Search

Coat colours and patterns are highly variable in cats and are determined mainly by several genes with Mendelian inheritance. A 2-bp deletion in agouti signalling protein (ASIP) is associated with melanism in domestic cats. Bengal cats are hybrids between domestic cats and Asian leopard cats (Prionailurus bengalensis), and the charcoal coat colouration/pattern in Bengals presents as a possible incomplete melanism. The complete coding region of ASIP was directly sequenced in Asian leopard, domestic and Bengal cats. Twenty-seven variants were identified between domestic and leopard cats and were investigated in Bengals and Savannahs, a hybrid with servals (Leptailurus serval). The leopard cat ASIP haplotype was distinguished from domestic cat by four synonymous and four non-synonymous exonic SNPs, as well as 19 intronic variants, including a 42-bp deletion in intron 4. Fifty-six of 64 reported charcoal cats were compound heterozygotes at ASIP, with leopard cat agouti (A(P) (be) ) and domestic cat non-agouti (a) haplotypes. Twenty-four Bengals had an additional unique haplotype (A2) for exon 2 that was not identified in leopard cats, servals or jungle cats (Felis chaus). The compound heterozygote state suggests the leopard cat allele, in combination with the recessive non-agouti allele, influences Bengal markings, producing a darker, yet not completely melanistic coat. This is the first validation of a leopard cat allele segregating in the Bengal breed and likely affecting their overall pelage phenotype. Genetic testing services need to be aware of the possible segregation of wild felid alleles in all assays performed on hybrid cats.

Chromosomal inversions are associated with local adaptation in many species. However, questions regarding how they are formed, maintained and impact various other evolutionary processes remain elusive. Here, using a large genomic dataset of long-read and short-read sequencing, we ask these questions in one of the most abundant vertebrates on Earth, the Atlantic herring. This species has four megabase-sized inversions associated with ecological adaptation that correlate with water temperature. The S and N inversion alleles at these four loci dominate in the southern and northern parts, respectively, of the species distribution in the North Atlantic Ocean. By determining breakpoint coordinates of the four inversions and the structural variations surrounding them, we hypothesize that these inversions are formed by ectopic recombination between duplicated sequences immediately outside of the inversions. We show that these are old inversions (>1 MY), albeit formed after the split between the Atlantic herring and its sister species, the Pacific herring. There is evidence for extensive gene flux between inversion alleles at all four loci. The large Ne of herring combined with the common occurrence of opposite homozygotes across the species distribution has allowed effective purifying selection to prevent the accumulation of genetic load and repeats within the inversions., (© 2024. The Author(s).)

The Greying with age phenotype in horses involves loss of hair pigmentation whereas skin pigmentation is not reduced, and a predisposition to melanoma. The causal mutation was initially reported as a duplication of a 4.6 kb intronic sequence in Syntaxin 17. The speed of greying varies considerably among Grey horses. Here we demonstrate the presence of two different Grey alleles, G2 carrying two tandem copies of the duplicated sequence and G3 carrying three. The latter is by far the most common allele, probably due to strong selection for the striking white phenotype. Our results reveal a remarkable dosage effect where the G3 allele is associated with fast greying and high incidence of melanoma whereas G2 is associated with slow greying and low incidence of melanoma. The copy number expansion transforms a weak enhancer to a strong melanocyte-specific enhancer that underlies hair greying (G2 and G3) and a drastically elevated risk of melanoma (G3 only). Our direct pedigree-based observation of the origin of a G2 allele from a G3 allele by copy number contraction demonstrates the dynamic evolution of this locus and provides the ultimate evidence for causality of the copy number variation of the 4.6 kb intronic sequence., (© 2024. The Author(s).)

Chromosomal abnormalities are a common cause of human miscarriage but rarely reported in any other species. As a result, there are currently inadequate animal models available to study this condition. Horses present one potential model since mares receive intense gynecological care. This allowed us to investigate the prevalence of chromosomal copy number aberrations in 256 products of conception (POC) in a naturally occurring model of pregnancy loss (PL). Triploidy (three haploid sets of chromosomes) was the most common aberration, found in 42% of POCs following PL over the embryonic period. Over the same period, trisomies and monosomies were identified in 11.6% of POCs and subchromosomal aberrations in 4.2%. Whole and subchromosomal aberrations involved 17 autosomes, with chromosomes 3, 4, and 20 having the highest number of aberrations. Triploid fetuses had clear gross developmental anomalies of the brain. Collectively, data demonstrate that alterations in chromosome number contribute to PL similarly in women and mares, with triploidy the dominant ploidy type over the key period of organogenesis. These findings, along with highly conserved synteny between human and horse chromosomes, similar gestation lengths, and the shared single greatest risk for PL being advancing maternal age, provide strong evidence for the first animal model to truly recapitulate many key features of human miscarriage arising due to chromosomal aberrations, with shared benefits for humans and equids., Competing Interests: Competing interests statement:The authors declare no competing interest.

Telomere-to-telomere (T2T) assemblies reveal new insights into the structure and function of the previously 'invisible' parts of the genome and allow comparative analyses of complete genomes across entire clades. We present here an open collaborative effort, termed the 'Ruminant T2T Consortium' (RT2T), that aims to generate complete diploid assemblies for numerous species of the Artiodactyla suborder Ruminantia to examine chromosomal evolution in the context of natural selection and domestication of species used as livestock., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Thoroughbred stallions that carry a double-homozygous genotype A/A-A/A for SNPs rs397316122 and rs69101140 in exon 5 of the FKBP6 gene (chr13; EquCab3.0) are uniquely subfertile due to impaired acrosomal exocytosis (IAE). In this study, the sperm proteome in frozen/thawed semen from subfertile Thoroughbred stallions was studied and compared to that of frozen/thawed sperm from fertile Thoroughbred stallions. A total of 2,220 proteins was identified, of which 140 proteins were found to be differentially abundant in sperm from the subfertile stallions compared to that of fertile stallions (83 less and 57 more abundant). Proteins of differential abundance in sperm from the subfertile stallions were mainly overrepresented in the "metabolism" and the "metabolism of lipids" pathways. One of these proteins, arylsulfatase F (ARSF), was studied by immunofluorescence. A lower proportion of sperm displaying ARSF signal at the acrosome region was observed in sperm from subfertile Thoroughbred stallions. In addition, heterologous zona pellucida binding assays revealed that sperm from subfertile Thoroughbred stallions bound at a lower proportion to zonae pellucidae than sperm from fertile Thoroughbred stallions. In conclusion, a group of differential abundance proteins, including some of acrosome origin, were identified in sperm from subfertile stallions with acrosome dysfunction., (© 2024. The Author(s).)

The employment of antibodies as a targeted drug delivery vehicle has proven successful which is exemplified by the emergence of antibody-drug conjugates (ADCs). However, ADCs are not without their shortcomings. Improvements may be made to the ADC platform by decoupling the cytotoxic drug from the delivery vehicle and conjugating an organometallic catalyst in its place. The resulting protein-metal catalyst conjugate was designed to uncage the masked cytotoxin administered as a separate entity. Macropinocytosis of albumin by cancerous cells suggests the potential of albumin acting as the tumor-targeting delivery vehicle. Herein reported are the first preparation and demonstration of ruthenium catalysts with cyclopentadienyl and quinoline-based ligands conjugated to albumin. The effective uncaging abilities were demonstrated on allyloxy carbamate (alloc)-protected rhodamine 110 and doxorubicin, providing a promising catalytic scaffold for the advancement of selective drug delivery methods in the future.

We generated single haplotype assemblies from a hinny hybrid which significantly improved the gapless contiguity for horse and donkey autosomal genomes and the X chromosomes. We added over 15 Mb of missing sequence to both X chromosomes, 60 Mb to donkey autosomes and corrected numerous errors in donkey and some in horse reference genomes. We resolved functionally important X-linked repeats: the DXZ4 macrosatellite and ampliconic Equine Testis Specific Transcript Y7 (ETSTY7). We pinpointed the location of the pseudoautosomal boundaries (PAB) and determined the size of the horse (1.8 Mb) and donkey (1.88 Mb) pseudoautosomal regions (PARs). We discovered distinct differences in horse and donkey PABs: a testis-expressed gene, XKR3Y, spans horse PAB with exons1-2 located in Y and exon3 in the X-Y PAR, whereas the donkey XKR3Y is Y-specific. DXZ4 had a similar ~ 8 kb monomer in both species with 10 copies in horse and 20 in donkey. We assigned hundreds of copies of ETSTY7, a sequence horizontally transferred from Parascaris and massively amplified in equids, to horse and donkey X chromosomes and three autosomes. The findings and products contribute to molecular studies of equid biology and advance research on X-linked conditions, sex chromosome regulation and evolution in equids., (© 2023. The Author(s).)

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer., (© 2023. The Author(s).)

The role of structurally dynamic genomic regions in speciation is poorly understood due to challenges inherent in diploid genome assembly. Here we reconstructed the evolutionary dynamics of structural variation in five cat species by phasing the genomes of three interspecies F1 hybrids to generate near-gapless single-haplotype assemblies. We discerned that cat genomes have a paucity of segmental duplications relative to great apes, explaining their remarkable karyotypic stability. X chromosomes were hotspots of structural variation, including enrichment with inversions in a large recombination desert with characteristics of a supergene. The X-linked macrosatellite DXZ4 evolves more rapidly than 99.5% of the genome clarifying its role in felid hybrid incompatibility. Resolved sensory gene repertoires revealed functional copy number changes associated with ecomorphological adaptations, sociality and domestication. This study highlights the value of gapless genomes to reveal structural mechanisms underpinning karyotypic evolution, reproductive isolation and ecological niche adaptation., (© 2023. The Author(s), under exclusive licence to Springer Nature America, Inc.)

We developed a highly contiguous chromosome-level reference genome for North American bison to provide a platform to evaluate the conservation, ecological, evolutionary, and population genomics of this species. Generated from a F1 hybrid between a North American bison dam and a domestic cattle bull, completeness and contiguity exceed that of other published bison genome assemblies. To demonstrate the utility for genome-wide variant frequency estimation, we compiled a genomic variant database consisting of 3 true albino bison and 44 wild-type pelage color bison. Through the examination of genomic variants fixed in the albino cohort and absent in the controls, we identified a nonsynonymous single nucleotide polymorphism (SNP) mutation on chromosome 29 in exon 3 of the tyrosinase gene (c.1114C>T). A TaqMan SNP Genotyping Assay was developed to genotype this SNP in a total of 283 animals across 29 herds. This assay confirmed the absence of homozygous variants in all animals except 7 true albino bison included in this study. In addition, the only heterozygous animals identified were 2 wild-type pelage color dams of albino offspring. Therefore, we propose that this new high-quality bison genome assembly and incipient variant database provides a highly robust and informative resource for genomics investigations for this iconic North American species., Competing Interests: Conflicts of interest statement The author(s) declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Background: The international Dog10K project aims to sequence and analyze several thousand canine genomes. Incorporating 20 × data from 1987 individuals, including 1611 dogs (321 breeds), 309 village dogs, 63 wolves, and four coyotes, we identify genomic variation across the canid family, setting the stage for detailed studies of domestication, behavior, morphology, disease susceptibility, and genome architecture and function., Results: We report the analysis of > 48 M single-nucleotide, indel, and structural variants spanning the autosomes, X chromosome, and mitochondria. We discover more than 75% of variation for 239 sampled breeds. Allele sharing analysis indicates that 94.9% of breeds form monophyletic clusters and 25 major clades. German Shepherd Dogs and related breeds show the highest allele sharing with independent breeds from multiple clades. On average, each breed dog differs from the UU_Cfam_GSD_1.0 reference at 26,960 deletions and 14,034 insertions greater than 50 bp, with wolves having 14% more variants. Discovered variants include retrogene insertions from 926 parent genes. To aid functional prioritization, single-nucleotide variants were annotated with SnpEff and Zoonomia phyloP constraint scores. Constrained positions were negatively correlated with allele frequency. Finally, the utility of the Dog10K data as an imputation reference panel is assessed, generating high-confidence calls across varied genotyping platform densities including for breeds not included in the Dog10K collection., Conclusions: We have developed a dense dataset of 1987 sequenced canids that reveals patterns of allele sharing, identifies likely functional variants, informs breed structure, and enables accurate imputation. Dog10K data are publicly available., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Introduction: Spontaneous rupture of tendons and ligaments is common in several species including humans. In horses, degenerative suspensory ligament desmitis (DSLD) is an important acquired idiopathic disease of a major energy-storing tendon-like structure. DSLD risk is increased in several breeds, including the Peruvian Horse. Affected horses have often been used for breeding before the disease is apparent. Breed predisposition suggests a substantial genetic contribution, but heritability and genetic architecture of DSLD have not been determined. Methods: To identify genomic regions associated with DSLD, we recruited a reference population of 183 Peruvian Horses, phenotyped as DSLD cases or controls, and undertook a genome-wide association study (GWAS), a regional window variance analysis using local genomic partitioning, a signatures of selection (SOS) analysis, and polygenic risk score (PRS) prediction of DSLD risk. We also estimated trait heritability from pedigrees. Results: Heritability was estimated in a population of 1,927 Peruvian horses at 0.22 ± 0.08. After establishing a permutation-based threshold for genome-wide significance, 151 DSLD risk single nucleotide polymorphisms (SNPs) were identified by GWAS. Multiple regions of enriched local heritability were identified across the genome, with strong enrichment signals on chromosomes 1, 2, 6, 10, 13, 16, 18, 22, and the X chromosome. With SOS analysis, there were 66 genes with a selection signature in DSLD cases that was not present in the control group that included the TGFB3 gene. Pathways enriched in DSLD cases included proteoglycan metabolism, extracellular matrix homeostasis, and signal transduction pathways that included the hedgehog signaling pathway. The best PRS predictive performance was obtained when we fitted 1% of top SNPs using a Bayesian Ridge Regression model which achieved the highest mean of R 2 on both the probit and logit liability scales, indicating a strong predictive performance. Discussion: We conclude that within-breed GWAS of DSLD in the Peruvian Horse has further confirmed that moderate heritability and a polygenic architecture underlies the trait and identified multiple DSLD SNP associations in novel tendinopathy candidate genes influencing disease risk. Pathways enriched with DSLD risk variants include ones that influence glycosaminoglycan metabolism, extracellular matrix homeostasis, signal transduction pathways., Competing Interests: PM, MM, SS, and SB are involved in setting up a genetic screening test for risk of DSLD in the Peruvian Horse at the University of Wisconsin-Madison. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Momen, Brauer, Patterson, Sample, Binversie, Davis, Cothran, Rosa, Brounts and Muir.)

The genomes of most vertebrates contain many V, D, and J gene segments within their Ig loci to construct highly variable CDR3 sequences through combinatorial diversity. This nucleotide variability translates into an antibody population containing extensive paratope diversity. Cattle have relatively few functional VDJ gene segments, requiring innovative approaches for generating diversity like the use of ultralong-encoding IGHV and IGHD gene segments that yield dramatically elongated CDR H3. Unique knob and stalk microdomains create protracted paratopes, where the antigen-binding knob sits atop a long stalk, allowing the antibody to bind both surface and recessed antigen epitopes. We examined genomes of twelve species of Bovidae to determine when ultralong-encoding IGHV and IGHD gene segments evolved. We located the 8-bp duplication encoding the unique TTVHQ motif in ultralong IGHV segments in six Bovid species (cattle, zebu, wild yak, domestic yak, American bison, and domestic gayal), but we did not find evidence of the duplication in species beyond the Bos and Bison genera. Additionally, we analyzed mRNA from bison spleen and identified a rich repertoire of expressed ultralong CDR H3 antibody mRNA, suggesting that bison use ultralong IGHV transcripts in their host defense. We found ultralong-encoding IGHD gene segments in all the same species except domestic yak, but again not beyond the Bos and Bison clade. Thus, the duplication event leading to this ultralong-encoding IGHV gene segment and the emergence of the ultralong-encoding IGHD gene segment appears to have evolved in a common ancestor of the Bos and Bison genera 5-10 million years ago., (© 2023. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Osteosarcoma prognosis has remained unchanged for the past three decades. In both humans and canines, treatment is limited to excision, radiation, and chemotherapy. Chemoresistance is the primary cause of treatment failure, and the trajectory of tumor evolution while under selective pressure from treatment is thought to be the major contributing factor in both species. We sought to understand the nature of platinum-based chemotherapy resistance by investigating cells that were subjected to repeated treatment and recovery cycles with increased carboplatin concentrations. Three HMPOS-derived cell lines, two resistant and one naïve, underwent single-cell RNA sequencing to examine transcriptomic perturbation and identify pathways leading to resistance and phenotypic changes. We identified the mechanisms of acquired chemoresistance and inferred the induced cellular trajectory that evolved with repeated exposure. The gene expression patterns indicated that acquired chemoresistance was strongly associated with a process similar to epithelial-mesenchymal transition (EMT), a phenomenon associated with the acquisition of migratory and invasive properties associated with metastatic disease. We conclude that the observed trajectory of tumor adaptability is directly correlated with chemoresistance and the phase of the EMT-like phenotype is directly affected by the level of chemoresistance. We infer that the EMT-like phenotype is a critical component of tumor evolution under treatment pressure and is vital to understanding the mechanisms of chemoresistance and to improving osteosarcoma prognosis.

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH , FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22 , IL26 , CSMD1 , LDHA , and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Binversie, Momen, Rosa, Davis and Muir.)

Structural rearrangements like copy number variations in the male-specific Y chromosome have been associated with male fertility phenotypes in human and mouse but have been sparsely studied in other mammalian species. Here, we designed digital droplet PCR assays for 7 horse male-specific Y chromosome multicopy genes and SRY and evaluated their absolute copy numbers in 209 normal male horses of 22 breeds, 73 XY horses with disorders of sex development and/or infertility, 5 Przewalski's horses and 2 kulans. This established baseline copy number for these genes in horses. The TSPY gene showed the highest copy number and was the most copy number variable between individuals and breeds. SRY was a single-copy gene in most horses but had 2-3 copies in some indigenous breeds. Since SRY is flanked by 2 copies of RBMY, their copy number variations were interrelated and may lead to SRY-negative XY disorders of sex development. The Przewalski's horse and kulan had 1 copy of SRY and RBMY. TSPY and ETSTY2 showed significant copy number variations between cryptorchid and normal males (P < 0.05). No significant copy number variations were observed in subfertile/infertile males. Notably, copy number of TSPY and ETSTY5 differed between successive male generations and between cloned horses, indicating germline and somatic mechanisms for copy number variations. We observed no correlation between male-specific Y chromosome gene copy number variations and male-specific Y chromosome haplotypes. We conclude that the ampliconic male-specific Y chromosome reference assembly has deficiencies and further studies with an improved male-specific Y chromosome assembly are needed to determine selective constraints over horse male-specific Y chromosome gene copy number and their relation to stallion reproduction and male biology., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Fluoroquinolones are a widely used class of chemotherapeutics within veterinary medicine, prized for their broad-spectrum bactericidal activity. These drugs present a known risk of retinal phototoxicity in domestic cats ( Felis catus ); therefore, using lower doses and alternative antibiotic classes is encouraged in this species. This adverse drug effect of fluoroquinolones, and enrofloxacin specifically, has been determined to be species-specific in domestic felids. Four feline-specific missense variants in ABCG2 result in four amino acid changes (E159M, S279L, H283Q, and T644I) that are unique to the domestic cat compared with multiple other nonfeline mammalian species. These changes alter the ABCG2 protein involved with the cellular transmembrane transport of drugs, including fluoroquinolones, making the protein functionally defective in domestic cats. The predisposition to fluoroquinolone-mediated phototoxicity in nondomestic felids was explored in this study. At least eight nondomestic felids share the four ABCG2 missense variants with domestic cats, and eleven other felids shared at least three of the four domestic cat variants. Taken together, these results suggest the genetic potential for nondomestic felids to also experience fluoroquinolone-induced retinal phototoxicity; therefore, cautions similar to those for domestic cats should be followed for these drugs in the entire feline taxon.

Degenerative suspensory ligament desmitis is a progressive idiopathic condition that leads to scarring and rupture of suspensory ligament fibers in multiple limbs in horses. The prevalence of degenerative suspensory ligament desmitis is breed related. Risk is high in the Peruvian Horse, whereas pony and draft breeds have low breed risk. Degenerative suspensory ligament desmitis occurs in families of Peruvian Horses, but its genetic architecture has not been definitively determined. We investigated contrasts between breeds with differing risk of degenerative suspensory ligament desmitis and identified associated risk variants and candidate genes. We analyzed 670k single nucleotide polymorphisms from 10 breeds, each of which was assigned one of the four breed degenerative suspensory ligament desmitis risk categories: control (Belgian, Icelandic Horse, Shetland Pony, and Welsh Pony), low risk (Lusitano, Arabian), medium risk (Standardbred, Thoroughbred, Quarter Horse), and high risk (Peruvian Horse). Single nucleotide polymorphisms were used for genome-wide association and selection signature analysis using breed-assigned risk levels. We found that the Peruvian Horse is a population with low effective population size and our breed contrasts suggest that degenerative suspensory ligament desmitis is a polygenic disease. Variant frequency exhibited signatures of positive selection across degenerative suspensory ligament desmitis breed risk groups on chromosomes 7, 18, and 23. Our results suggest degenerative suspensory ligament desmitis breed risk is associated with disturbances to suspensory ligament homeostasis where matrix responses to mechanical loading are perturbed through disturbances to aging in tendon (PIN1), mechanotransduction (KANK1, KANK2, JUNB, SEMA7A), collagen synthesis (COL4A1, COL5A2, COL5A3, COL6A5), matrix responses to hypoxia (PRDX2), lipid metabolism (LDLR, VLDLR), and BMP signaling (GREM2). Our results do not suggest that suspensory ligament proteoglycan turnover is a primary factor in disease pathogenesis., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

Small effective population sizes raise the probability of extinction by increasing the frequency of potentially deleterious alleles and reducing fitness. However, the extent to which cancers play a role in the fitness reduction of genetically depauperate wildlife populations is unknown. Santa Catalina island foxes ( Urocyon littoralis catalinae ) sampled in 2007-2008 have a high prevalence of ceruminous gland tumors, which was not detected in the population prior to a recent bottleneck caused by a canine distemper epidemic. The disease appears to be associated with inflammation from chronic ear mite ( Otodectes ) infections and secondary elevated levels of Staphyloccus pseudointermedius bacterial infections. However, no other environmental factors to date have been found to be associated with elevated cancer risk in this population. Here, we used whole genome sequencing of the case and control individuals from two islands to identify candidate loci associated with cancer based on genetic divergence, nucleotide diversity, allele frequency spectrum, and runs of homozygosity. We identified several candidate loci based on genomic signatures and putative gene functions, suggesting that cancer susceptibility in this population may be polygenic. Due to the efforts of a recovery program and weak fitness effects of late-onset disease, the population size has increased, which may allow selection to be more effective in removing these presumably slightly deleterious alleles. Long-term monitoring of the disease alleles, as well as overall genetic diversity, will provide crucial information for the long-term persistence of this threatened population.

Recent adaptive radiations are models for investigating mechanisms contributing to the evolution of biodiversity. An unresolved question is the relative importance of new mutations, ancestral variants, and introgressive hybridization for phenotypic evolution and speciation. Here, we address this issue using Darwin's finches and investigate the genomic architecture underlying their phenotypic diversity. Admixture mapping for beak and body size in the small, medium, and large ground finches revealed 28 loci showing strong genetic differentiation. These loci represent ancestral haplotype blocks with origins predating speciation events during the Darwin's finch radiation. Genes expressed in the developing beak are overrepresented in these genomic regions. Ancestral haplotypes constitute genetic modules for selection and act as key determinants of the unusual phenotypic diversity of Darwin's finches. Such ancestral haplotype blocks can be critical for how species adapt to environmental variability and change.

During the late nineteenth century North American bison underwent a significant population bottleneck resulting in a reduction in population size of over 99% and a species-level near-extinction event. Factors responsible for this destruction included indiscriminate killing, loss of access to suitable habitat, and diseases. At the nadir of this population crash, very few wild plains bison survived and were restricted to Yellowstone National Park, USA and a small number of wild wood bison remained in Wood Buffalo National Park, Canada. However, most surviving bison in the late 1800's were maintained by cattle ranchers in private herds where hybridization between bison with various breeds of domestic cattle was often encouraged. Over the last 20 years, the legacy of this introgression has been identified using mitochondrial DNA and limited nuclear microsatellite analyses. However, no genome-wide assessment has been performed, and some herds were believed to be free of introgression based on current genetic testing strategies. Herein, we report detailed analyses using whole genome sequencing from nineteen modern and six historical bison, chosen to represent the major lineages of bison, to identify and quantitate signatures of nuclear introgression in their recent (within 200 years) history. Both low and high coverage genomes provided evidence for recent introgression, including animals from Yellowstone, Wind Cave, and Elk Island National Parks which were previously thought to be free from hybridization with domestic cattle. We employed multiple approaches, including one developed for this work, to identify putative cattle haplotypes in each bison genome. These regions vary greatly in size and frequency by sample and herd, though we detected domestic cattle introgression in all bison genomes tested. Since our sampling strategy spanned across the diversity of modern bison populations, these finding are best explained by multiple historical hybridization events between these two species with significant genetic recombination over the last 200 years. Our results demonstrate that whole genome sequencing approaches are required to accurately quantitate cattle introgression in bison., (© 2022. The Author(s).)

Allergy to domestic cat affects up to 15% of the population, and sensitization to cat allergen is associated with asthma. Despite the pervasiveness of cat allergic disease, current treatments have limited impact. Here, we present a bioinformatics analysis of the major cat allergen, Fel d 1, and demonstrate proof of principle for CRISPR gene editing of the allergen. Sequence and structural analyses of Fel d 1 from 50 domestic cats identified conserved coding regions in genes CH1 and CH2 suitable for CRISPR editing. Comparative analyses of Fel d 1 and orthologous sequences from eight exotic felid species determined relatively low-sequence identities for CH1 and CH2 , and implied that the allergen may be nonessential for cats, given the apparent lack of evolutionary conservation. In vitro knockouts of domestic cat Fel d 1 using CRISPR-Cas9 yielded editing efficiencies of up to 55% and found no evidence of editing at predicted potential off-target sites. Taken together, our data indicate that Fel d 1 is both a rational and viable candidate for gene deletion, which may profoundly benefit cat allergy sufferers by removing the major allergen at the source.

The high mobility group AT-hook 2 (HMGA2) protein works as an architectural regulator by binding AT-rich DNA sequences to induce conformational changes affecting transcription. Genomic deletions disrupting HMGA2 coding sequences and flanking noncoding sequences cause dwarfism in mice and rabbits. Here, CRISPR/Cas9 was used in mice to generate an Hmga2 null allele that specifically disrupts only the coding sequence. The loss of one or both alleles of Hmga2 resulted in reduced body size of 20% and 60%, respectively, compared to wild-type littermates as well as an allometric reduction in skull length in Hmga2-/- mice. Both male and female Hmga2-/- mice are infertile, whereas Hmga2+/- mice are fertile. Examination of reproductive tissues of Hmga2-/- males revealed a significantly reduced size of testis, epididymis, and seminal vesicle compared to controls, and 70% of knock-out males showed externalized penis, but no cryptorchidism was observed. Sperm analyses revealed severe oligospermia in mutant males and slightly decreased sperm viability, increased DNA damage but normal sperm chromatin compaction. Testis histology surprisingly revealed a normal seminiferous epithelium, despite the significant reduction in testis size. In addition, Hmga2-/- mice showed a significantly reduced exploratory behavior. In summary, the phenotypic effects in mouse using targeted mutagenesis confirmed that Hmga2 is affecting prenatal and postnatal growth regulation, male reproductive tissue development, and presents the first indication that Hmga2 function is required for normal mouse behavior. No specific effect, despite an allometric reduction, on craniofacial development was noted in contrast to previous reports of an altered craniofacial development in mice and rabbits carrying deletions of both coding and noncoding sequences at the 5' part of Hmga2., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Melanotic ( Ml ) is a mutation in chickens that extends black (eumelanin) pigmentation in normally brown or red (pheomelanin) areas, thus affecting multiple within-feather patterns [J. W. Moore, J. R. Smyth Jr, J. Hered. 62, 215-219 (1971)]. In the present study, linkage mapping using a back-cross between Dark Cornish ( Ml / Ml ) and Partridge Plymouth Rock ( ml + /ml + ) chickens assigned Ml to an 820-kb region on chromosome 1. Identity-by-descent mapping, via whole-genome sequencing and diagnostic tests using a diverse set of chickens, refined the localization to the genomic region harboring GJA5 encoding gap-junction protein 5 (alias connexin 40) previously associated with pigmentation patterns in zebrafish. An insertion/deletion polymorphism located in the vicinity of the GJA5 promoter region was identified as the candidate causal mutation. Four different GJA5 transcripts were found to be expressed in feather follicles and at least two showed differential expression between genotypes. The results showed that Melanotic constitutes a cis -acting regulatory mutation affecting GJA5 expression. A recent study established the melanocortin-1 receptor ( MC1R) locus and the interaction between the MC1R receptor and its antagonist agouti-signaling protein as the primary mechanism underlying variation in within-feather pigmentation patterns in chickens. The present study advances understanding the mechanisms underlying variation in plumage color in birds because it demonstrates that the activity of connexin 40/GJA5 can modulate the periodic pigmentation patterns within individual feathers., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3 , which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds., (© 2021 Serres-Armero et al.; Published by Cold Spring Harbor Laboratory Press.)

In 2013, the U.S. Department of Agriculture Food Safety and Inspection Service (USDA-FSIS) began transitioning to whole genome sequencing (WGS) for foodborne disease outbreak- and recall-associated isolate identification of select bacterial species. While WGS offers greater precision, certain hurdles must be overcome before widespread application within the food industry is plausible. Challenges include diversity of sequencing platform outputs and lack of standardized bioinformatics workflows for data analyses. We sequenced DNA from USDA-FSIS approved, non-pathogenic E. coli surrogates and a derivative group of rifampicin-resistant mutants (rif R ) via both Oxford Nanopore MinION and Illumina MiSeq platforms to generate and annotate complete genomes. Genome sequences from each clone were assembled separately so long-read, short-read, and combined sequence assemblies could be directly compared. The combined sequence data approach provides more accurate completed genomes. The genomes from these isolates were verified to lack functional key E. coli elements commonly associated with pathogenesis. Genetic alterations known to confer rif R were also identified. As the food industry adopts WGS within its food safety programs, these data provide completed genomes for commonly used surrogate strains, with a direct comparison of sequence platforms and assembly strategies relevant to research/testing workflows applicable for both processors and regulators.

The Crest mutation in chicken shows incomplete dominance and causes a spectacular phenotype in which the small feathers normally present on the head are replaced by much larger feathers normally present only in dorsal skin. Using whole-genome sequencing, we show that the crest phenotype is caused by a 197 bp duplication of an evolutionarily conserved sequence located in the intron of HOXC10 on chromosome 33. A diagnostic test showed that the duplication was present in all 54 crested chickens representing eight breeds and absent from all 433 non-crested chickens representing 214 populations. The mutation causes ectopic expression of at least five closely linked HOXC genes, including HOXC10, in cranial skin of crested chickens. The result is consistent with the interpretation that the crest feathers are caused by an altered body region identity. The upregulated HOXC gene expression is expanded to skull tissue of Polish chickens showing a large crest often associated with cerebral hernia, but not in Silkie chickens characterized by a small crest, both homozygous for the duplication. Thus, the 197 bp duplication is required for the development of a large crest and susceptibility to cerebral hernia because only crested chicken show this malformation. However, this mutation is not sufficient to cause herniation because this malformation is not present in breeds with a small crest, like Silkie chickens., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

The unique evolutionary dynamics and complex structure make the Y chromosome the most diverse and least understood region in the mammalian genome, despite its undisputable role in sex determination, development, and male fertility. Here we present the first contig-level annotated draft assembly for the alpaca ( Vicugna pacos ) Y chromosome based on hybrid assembly of short- and long-read sequence data of flow-sorted Y. The latter was also used for cDNA selection providing Y-enriched testis transcriptome for annotation. The final assembly of 8.22 Mb comprised 4.5 Mb of male specific Y (MSY) and 3.7 Mb of the pseudoautosomal region. In MSY, we annotated 15 X-degenerate genes and two novel transcripts, but no transposed sequences. Two MSY genes, HSFY and RBMY, are multicopy. The pseudoautosomal boundary is located between SHROOM2 and HSFY. Comparative analysis shows that the small and cytogenetically distinct alpaca Y shares most of MSY sequences with the larger dromedary and Bactrian camel Y chromosomes. Most of alpaca X-degenerate genes are also shared with other mammalian MSYs, though WWC3Y is Y-specific only in alpaca/camels and the horse. The partial alpaca Y assembly is a starting point for further expansion and will have applications in the study of camelid populations and male biology.

Direct comparisons between historical and contemporary populations allow for detecting changes in genetic diversity through time and assessment of the impact of habitat fragmentation. Here, we determined the genetic architecture of both historical and modern lions to document changes in genetic diversity over the last century. We surveyed microsatellite and mitochondrial genome variation from 143 high-quality museum specimens of known provenance, allowing us to directly compare this information with data from several recently published nuclear and mitochondrial studies. Our results provide evidence for male-mediated gene flow and recent isolation of local subpopulations, likely due to habitat fragmentation. Nuclear markers showed a significant decrease in genetic diversity from the historical (HE = 0.833) to the modern (HE = 0.796) populations, whereas mitochondrial genetic diversity was maintained (Hd = 0.98 for both). Although the historical population appears to have been panmictic based on nDNA data, hierarchical structure analysis identified four tiers of genetic structure in modern populations and was able to detect most sampling locations. Mitogenome analyses identified four clusters: Southern, Mixed, Eastern, and Western and were consistent between modern and historically sampled haplotypes. Within the last century, habitat fragmentation caused lion subpopulations to become more geographically isolated as human expansion changed the African landscape. This resulted in an increase in fine-scale nuclear genetic structure and loss of genetic diversity as lion subpopulations became more differentiated, whereas mitochondrial structure and diversity were maintained over time., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Atlantic herring is widespread in North Atlantic and adjacent waters and is one of the most abundant vertebrates on earth. This species is well suited to explore genetic adaptation due to minute genetic differentiation at selectively neutral loci. Here, we report hundreds of loci underlying ecological adaptation to different geographic areas and spawning conditions. Four of these represent megabase inversions confirmed by long read sequencing. The genetic architecture underlying ecological adaptation in herring deviates from expectation under a classical infinitesimal model for complex traits because of large shifts in allele frequencies at hundreds of loci under selection., Competing Interests: FH, MJ, MP, LS, AF, BD, DB, FB, MC, GD, EF, AF, LA No competing interests declared, (© 2020, Han et al.)

The outcome of pathogen spillover from a reservoir to a novel host population can range from a "dead-end" when there is no onward transmission in the recipient population, to epidemic spread and even establishment in new hosts. Understanding the evolutionary epidemiology of spillover events leading to discrete outcomes in novel hosts is key to predicting risk and can lead to a better understanding of the mechanisms of emergence. Here we use a Bayesian phylodynamic approach to examine cross-species transmission and evolutionary dynamics during a canine distemper virus (CDV) spillover event causing clinical disease and population decline in an African lion population (Panthera leo) in the Serengeti Ecological Region between 1993 and 1994. Using 21 near-complete viral genomes from four species we found that this large-scale outbreak was likely  ignited by a single cross-species spillover event from a canid reservoir to noncanid hosts <1 year before disease detection and explosive spread of CDV in lions. Cross-species transmission from other noncanid species probably fuelled the high prevalence of CDV across spatially structured lion prides. Multiple lines of evidence suggest that spotted hyenas (Crocuta crocuta) could have acted as the proximate source of CDV exposure in lions. We report 13 nucleotide substitutions segregating CDV strains found in canids and noncanids. Our results are consistent with the hypothesis that virus evolution played a role in CDV emergence in noncanid hosts following spillover during the outbreak, suggest that host barriers to clinical infection can limit outcomes of CDV spillover in novel host species., (© 2020 John Wiley & Sons Ltd.)

Canine distemper virus (CDV) is a multi-host pathogen with variable clinical outcomes of infection across and within species. We used whole-genome sequencing (WGS) to search for viral markers correlated with clinical distemper in African lions. To identify candidate markers, we first documented single-nucleotide polymorphisms (SNPs) differentiating CDV strains associated with different clinical outcomes in lions in East Africa. We then conducted evolutionary analyses on WGS from all global CDV lineages to identify loci subject to selection. SNPs that both differentiated East African strains and were under selection were mapped to a phylogenetic tree representing global CDV diversity to assess if candidate markers correlated with documented outbreaks of clinical distemper in lions ( n = 3). Of 54 SNPs differentiating East African strains, ten were under positive or episodic diversifying selection and 20 occurred in the clinical strain despite strong purifying selection at those loci. Candidate markers were in functional domains of the RNP complex ( n = 19), the matrix protein ( n = 4), on CDV glycoproteins ( n = 5), and on the V protein ( n = 1). We found mutations at two loci in common between sequences from three CDV outbreaks of clinical distemper in African lions; one in the signaling lymphocytic activation molecule receptor (SLAM)-binding region of the hemagglutinin protein and another in the catalytic center of phosphodiester bond formation on the large polymerase protein. These results suggest convergent evolution at these sites may have a functional role in clinical distemper outbreaks in African lions and uncover potential novel barriers to pathogenicity in this species., Competing Interests: The authors declare no conflicts of interest.

The domestic cat (Felis catus) numbers over 94 million in the USA alone, occupies households as a companion animal, and, like humans, suffers from cancer and common and rare diseases. However, genome-wide sequence variant information is limited for this species. To empower trait analyses, a new cat genome reference assembly was developed from PacBio long sequence reads that significantly improve sequence representation and assembly contiguity. The whole genome sequences of 54 domestic cats were aligned to the reference to identify single nucleotide variants (SNVs) and structural variants (SVs). Across all cats, 16 SNVs predicted to have deleterious impacts and in a singleton state were identified as high priority candidates for causative mutations. One candidate was a stop gain in the tumor suppressor FBXW7. The SNV is found in cats segregating for feline mediastinal lymphoma and is a candidate for inherited cancer susceptibility. SV analysis revealed a complex deletion coupled with a nearby potential duplication event that was shared privately across three unrelated cats with dwarfism and is found within a known dwarfism associated region on cat chromosome B1. This SV interrupted UDP-glucose 6-dehydrogenase (UGDH), a gene involved in the biosynthesis of glycosaminoglycans. Importantly, UGDH has not yet been associated with human dwarfism and should be screened in undiagnosed patients. The new high-quality cat genome reference and the compilation of sequence variation demonstrate the importance of these resources when searching for disease causative alleles in the domestic cat and for identification of feline biomedical models., Competing Interests: The authors disclose there are no conflicts of interest. The funders had no role in study design, data collection, data analysis, interpretation of results, or decision to publish.