Search

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates. [ABSTRACT FROM AUTHOR]

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.)

Variant calling across diverse species remains challenging as most bioinformatics tools default to assumptions based on human genomes. DeepVariant (DV) excels without joint genotyping while offering fewer implementation barriers. However, the growing appeal of a "universal" algorithm has magnified the unknown impacts when used with non-human genomes. Here, we use bovine genomes to assess the limits of human-genome-trained models in other species. We introduce the first multi-species DV model that achieves a lower Mendelian Inheritance Error (MIE) rate during single-sample genotyping. Our novel approach, TrioTrain, automates extending DV for species without Genome In A Bottle (GIAB) resources and uses region shuffling to mitigate barriers for SLURM-based clusters. To offset imperfect truth labels for animal genomes, we remove Mendelian discordant variants before training, where models are tuned to genotype the offspring correctly. With TrioTrain, we use cattle, yak, and bison trios to build 30 model iterations across five phases. We observe remarkable performance across phases when testing the GIAB human trios with a mean SNP F1 score >0.990. In HG002, our phase 4 bovine model identifies more variants at a lower MIE rate than DeepTrio. In bovine F1-hybrid genomes, our model substantially reduces inheritance errors with a mean MIE rate of 0.03 percent. Although constrained by imperfect labels, we find that multi-species, trio-based training produces a robust variant calling model. Our research demonstrates that exclusively training with human genomes restricts the application of deep-learning approaches for comparative genomics., Competing Interests: COMPETING INTEREST STATEMENT The authors declare that they have no competing interests.

Directional selection alters the genome via hard sweeps, soft sweeps, and polygenic selection. However, mapping polygenic selection is difficult because it does not leave clear signatures on the genome like a selective sweep. In populations with temporally stratified genotypes, the Generation Proxy Selection Mapping (GPSM) method identifies variants associated with generation number (or appropriate proxy) and thus variants undergoing directional allele frequency changes. Here, we use GPSM on two large datasets of beef cattle to detect associations between an animal's generation and 11 million imputed SNPs. Using these datasets with high power and dense mapping resolution, GPSM detected a total of 294 unique loci actively under selection in two cattle breeds. We observed that GPSM has a high power to detect selection in the very recent past (<10 years), even when allele frequency changes are small. Variants identified by GPSM reside in genomic regions associated with known breed-specific selection objectives, such as fertility and maternal ability in Red Angus, and carcass merit and coat color in Simmental. Over 60% of the selected loci reside in or near (<50 kb) annotated genes. Using haplotype-based and composite selective sweep statistics, we identify hundreds of putative selective sweeps that likely occurred earlier in the evolution of these breeds; however, these sweeps have little overlap with recent polygenic selection. This makes GPSM a complementary approach to sweep detection methods when temporal genotype data are available. The selected loci that we identify across methods demonstrate the complex architecture of selection in domesticated cattle., Competing Interests: The authors declare no conflicts of interest related to this work., (© 2024 The Authors. Evolutionary Applications published by John Wiley & Sons Ltd.)

Creation of the bovine reference assembly paved the way to develop the high-throughput genotyping arrays of the single nucleotide polymorphisms (SNPs) based on the available map coordinates that facilitated major advances in gene mapping and selection programs. The assembly flaws, however, may cause false results in the downstream gene mapping studies. The most recent bovine reference genome (ARS-UCD1.2) is built on the long-read sequences that provides improved quality and continuity. By applying population genetic metrics in this study, we aimed to evaluate the map coordinates to which SNP markers were assigned. We employed a three-step approach by combining the recombination and linkage disequilibrium analyses to test if the markers fit into the assigned physical map coordinates. We applied the method to the bovine 50k array in a large pedigree of Holstein cattle and revealed a panel of 65 candidate markers, most of which were re-located either on a different chromosome or re-mapped as far as several million base pairs away on the same chromosome. This list of candidates accounts for 0.1% of the SNPs in the widely used 50k genotyping array and we foresee a reasonably larger set of markers being misplaced in the BovineHD 700K BeadChip. We suggest pre-removal of the candidate misplaced markers to reduce false signals in association mapping studies., (© 2022 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Background: Genotypic information produced from single nucleotide polymorphism (SNP) arrays has routinely been used to identify genomic regions associated with complex traits in beef and dairy cattle. Herein, we assembled a dataset consisting of 15,815 Red Angus beef cattle distributed across the continental U.S. and a union set of 836,118 imputed SNPs to conduct genome-wide association analyses (GWAA) for growth traits using univariate linear mixed models (LMM); including birth weight, weaning weight, and yearling weight. Genomic relationship matrix heritability estimates were produced for all growth traits, and genotype-by-environment (GxE) interactions were investigated., Results: Moderate to high heritabilities with small standard errors were estimated for birth weight (0.51 ± 0.01), weaning weight (0.25 ± 0.01), and yearling weight (0.42 ± 0.01). GWAA revealed 12 pleiotropic QTL (BTA6, BTA14, BTA20) influencing Red Angus birth weight, weaning weight, and yearling weight which met a nominal significance threshold (P ≤ 1e-05) for polygenic traits using 836K imputed SNPs. Moreover, positional candidate genes associated with Red Angus growth traits in this study (i.e., LCORL, LOC782905, NCAPG, HERC6, FAM184B, SLIT2, MMRN1, KCNIP4, CCSER1, GRID2, ARRDC3, PLAG1, IMPAD1, NSMAF, PENK, LOC112449660, MOS, SH3PXD2B, STC2, CPEB4) were also previously associated with feed efficiency, growth, and carcass traits in beef cattle. Collectively, 14 significant GxE interactions were also detected, but were less consistent among the investigated traits at a nominal significance threshold (P ≤ 1e-05); with one pleiotropic GxE interaction detected on BTA28 (24 Mb) for Red Angus weaning weight and yearling weight., Conclusions: Sixteen well-supported QTL regions detected from the GWAA and GxE GWAA for growth traits (birth weight, weaning weight, yearling weight) in U.S. Red Angus cattle were found to be pleiotropic. Twelve of these pleiotropic QTL were also identified in previous studies focusing on feed efficiency and growth traits in multiple beef breeds and/or their composites. In agreement with other beef cattle GxE studies our results implicate the role of vasodilation, metabolism, and the nervous system in the genetic sensitivity to environmental stress., (© 2022. The Author(s).)

This study is part of a concerted effort to identify and phenotype rare, deleterious mutations that adversely affect sperm quality, or convey high developmental and fertility potential to embryos and ensuing progeny. A rare, homozygous mutation in EML5 ( EML5 R1654W ), which encodes a microtubule-associated protein with high expression in testis and brain was identified in an Angus bull used extensively in artificial insemination (AI) for its outstanding progeny production traits. The bull's fertility was low in cross-breeding timed AI (TAI) (Pregnancy/TAI = 25.2%; n = 222) and, in general, AI breeding to Nellore cows (41%; n = 822). A search of the 1,000 Bull Genomes Run9 database revealed an additional 74 heterozygous animals and 8 homozygous animals harboring this exact mutation across several different breeds (0.7% frequency within the 6,191 sequenced animals). Phenotypically, spermatozoa from the homozygous Angus bull displayed prominent piriform and tapered heads, and outwardly protruding knobbed acrosomes. Additionally, an increased retention of EML5 was also observed in the sperm head of both homozygous and heterozygous Angus bulls compared to wild-type animals. This non-synonymous point mutation is located within a WD40 signaling domain repeat of EML5 and is predicted to be detrimental to overall protein function by genomic single nucleotide polymorphism (SNP) analysis and protein modeling. Future work will examine how this rare mutation affects field AI fertility and will characterize the role of EML5 in spermatogenesis., 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 Nogueira, Tirpák, Hamilton, Zigo, Kerns, Sutovsky, Kim, Volkmann, Jovine, Taylor, Schnabel and Sutovsky.)

The Neosho madtom (Noturus placidus) is a small catfish, generally less than 3 inches in length, unique to the Neosho-Spring River system within the Arkansas River Basin. It was federally listed as threatened in 1990, largely due to habitat loss. For conservation efforts, we generated whole-genome sequence data from 10 Neosho madtom individuals originating from 3 geographically separated populations to evaluate genetic diversity and population structure. A Neosho madtom genome was de novo assembled, and genome size and content were assessed. Single nucleotide polymorphisms were assessed from de Bruijn graphs, and via reference alignment with both the channel catfish (Ictalurus punctatus) reference genome and Neosho madtom reference genome. Principal component analysis and structure analysis indicated weak population structure, suggesting fish from the 3 locations represent a single population. Using a novel method, genome-wide conservation and divergence between the Neosho madtom, channel catfish, and zebrafish (Danio rerio) was assessed by pairwise contig alignment, which demonstrated that genes important to embryonic development frequently had conserved sequences. This research in a threatened species with no previously published genomic resources provides novel genetic information to guide current and future conservation efforts and demonstrates that using whole-genome sequencing provides detailed information of population structure and demography using only a limited number of rare and valuable samples., (Published by Oxford University Press on behalf of Genetics Society of America 2022. This work is written by US Government employees and is in the public domain in the US.)

Development of the American Breeds of beef cattle began in the 1920s as breeders and U. S. Experiment Station researchers began to create Bos taurus taurus × Bos taurus indicus hybrids using Brahman as the B. t. indicus source. By 1954, U.S. Breed Associations had been formed for Brangus (5/8 Angus × 3/8 Brahman), Beefmaster (½ Brahman × ¼ Shorthorn × ¼ Hereford), and Santa Gertrudis (5/8 Shorthorn × 3/8 Brahman). While these breeds were developed using mating designs expected to create base generation animals with the required genome contributions from progenitor breeds, each association has now registered advanced generation animals in which selection or drift may have caused the realized genome compositions to differ from initial expected proportions. The availability of high-density SNP genotypes for 9,161 Brangus, 3,762 Beefmaster, and 1,942 Santa Gertrudis animals allowed us to compare the realized genomic architectures of breed members to the base generation expectations. We used RFMix to estimate local ancestry and identify genomic regions in which the proportion of Brahman ancestry differed significantly from a priori expectations. For all three breeds, lower than expected levels of Brahman composition were found genome-wide, particularly in early-generation animals where we demonstrate that selection on beef production traits was likely responsible for the taurine enrichment. Using a proxy for generation number, we also contrasted the genomes of early- and advanced-generation animals and found that the indicine composition of the genome has increased with generation number likely due to selection on adaptive traits. Many of the most-highly differentiated genomic regions were breed specific, suggesting that differences in breeding objectives and selection intensities exist between the breeds. Global ancestry estimation is commonly performed in admixed animals to control for stratification in association studies. However, local ancestry estimation provides the opportunity to investigate the evolution of specific chromosomal segments and estimate haplotype effects on trait variation in admixed individuals. Investigating the genomic architecture of the American Breeds not only allows the estimation of indicine and taurine genome proportions genome-wide, but also the locations within the genome where either taurine or indicine alleles confer a selective advantage., 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 © 2021 Crum, Schnabel, Decker and Taylor.)

Selection on complex traits can rapidly drive evolution, especially in stressful environments. This polygenic selection does not leave intense sweep signatures on the genome, rather many loci experience small allele frequency shifts, resulting in large cumulative phenotypic changes. Directional selection and local adaptation are changing populations; but, identifying loci underlying polygenic or environmental selection has been difficult. We use genomic data on tens of thousands of cattle from three populations, distributed over time and landscapes, in linear mixed models with novel dependent variables to map signatures of selection on complex traits and local adaptation. We identify 207 genomic loci associated with an animal's birth date, representing ongoing selection for monogenic and polygenic traits. Additionally, hundreds of additional loci are associated with continuous and discrete environments, providing evidence for historical local adaptation. These candidate loci highlight the nervous system's possible role in local adaptation. While advanced technologies have increased the rate of directional selection in cattle, it has likely been at the expense of historically generated local adaptation, which is especially problematic in changing climates. When applied to large, diverse cattle datasets, these selection mapping methods provide an insight into how selection on complex traits continually shapes the genome. Further, understanding the genomic loci implicated in adaptation may help us breed more adapted and efficient cattle, and begin to understand the basis for mammalian adaptation, especially in changing climates. These selection mapping approaches help clarify selective forces and loci in evolutionary, model, and agricultural contexts., Competing Interests: The authors have declared that no competing interests exist.

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations., (Copyright © 2021. Published by Elsevier Inc.)

Background: Major advances in selection progress for cattle have been made following the introduction of genomic tools over the past 10-12 years. These tools depend upon the Bos taurus reference genome (UMD3.1.1), which was created using now-outdated technologies and is hindered by a variety of deficiencies and inaccuracies., Results: We present the new reference genome for cattle, ARS-UCD1.2, based on the same animal as the original to facilitate transfer and interpretation of results obtained from the earlier version, but applying a combination of modern technologies in a de novo assembly to increase continuity, accuracy, and completeness. The assembly includes 2.7 Gb and is >250× more continuous than the original assembly, with contig N50 >25 Mb and L50 of 32. We also greatly expanded supporting RNA-based data for annotation that identifies 30,396 total genes (21,039 protein coding). The new reference assembly is accessible in annotated form for public use., Conclusions: We demonstrate that improved continuity of assembled sequence warrants the adoption of ARS-UCD1.2 as the new cattle reference genome and that increased assembly accuracy will benefit future research on this species., (© The Author(s) 2020. Published by Oxford University Press.)

Background: During the last decade, the use of common-variant array-based single nucleotide polymorphism (SNP) genotyping in the beef and dairy industries has produced an astounding amount of medium-to-low density genomic data. Although low-density assays work well in the context of genomic prediction, they are less useful for detecting and mapping causal variants and the effects of rare variants are not captured. The objective of this project was to maximize the accuracies of genotype imputation from medium- and low-density assays to the marker set obtained by combining two high-density research assays (~ 850,000 SNPs), the Illumina BovineHD and the GGP-F250 assays, which contains a large proportion of rare and potentially functional variants and for which the assay design is described here. This 850 K SNP set is useful for both imputation to sequence-level genotypes and direct downstream analysis., Results: We found that a large multi-breed composite imputation reference panel that includes 36,131 samples with either BovineHD and/or GGP-F250 genotypes significantly increased imputation accuracy compared with a within-breed reference panel, particularly at variants with low minor allele frequencies. Individual animal imputation accuracies were maximized when more genetically similar animals were represented in the composite reference panel, particularly with complete 850 K genotypes. The addition of rare variants from the GGP-F250 assay to our composite reference panel significantly increased the imputation accuracy of rare variants that are exclusively present on the BovineHD assay. In addition, we show that an assay marker density of 50 K SNPs balances cost and accuracy for imputation to 850 K., Conclusions: Using high-density genotypes on all available individuals in a multi-breed reference panel maximized imputation accuracy for tested cattle populations. Admixed animals or those from breeds with a limited representation in the composite reference panel were still imputed at high accuracy, which is expected to further increase as the reference panel expands. We anticipate that the addition of rare variants from the GGP-F250 assay will increase the accuracy of imputation to sequence level.

Two littermate German Shorthaired Pointers, a male and a female, were adopted as puppies from an animal shelter. Both puppies developed normally until approximately 11 months of age when the male began to exhibit neurological signs including ataxia, vision loss, and behavioral changes indicative of cognitive decline. These signs increased in severity over time. The female remained neurologically normal and healthy. The affected dog was euthanized at approximately 21 months of age. Autofluorescent cytoplasmic storage bodies were detected in neurons in unstained tissue sections from the cerebellum, the cerebrum, and the retina. Electron micrographs of these storage bodies showed that they were membrane bound and that most contained tightly packed aggregates of membranous whorls along with a variety of other ultrastructural features. This ultrastructure, along with the autofluorescence and the clinical signs supported a diagnosis of neuronal ceroid lipofuscinosis (NCL). Unlike earlier investigated forms of canine NCL with causal alleles in ATP13A2 , TPP1 , MFSD8 and CLN5 that had autofluorescent cytoplasmic storage bodies in cardiac muscle, no autofluorescence was detected in cardiac muscle from the affected German Shorthaired Pointer. A 39-fold average coverage whole genome sequence indicated that the affected German Shorthaired Pointer was homozygous for the A allele of a G > A transversion at position 30,895,648 chromosome 37. This 37:30895648G > A mutation created a CLN8 termination codon that had been previously reported to cause NCL in a mixed breed dog with Australian Shepherd and Australian Cattle Dog ancestry. This nonsense allele was heterozygous in the clinically normal female sibling, while archived DNA samples from 512 other German Shorthaired Pointers were all homozygous for the reference allele. The affected German Shorthaired Pointer and the previously diagnosed mixed breed dog with the same nonsense mutation shaired an identical homozygous haplotype that extended for 4.41 Mb at the telomeric end of chromosome 37, indicating the both dogs inherited the nonsense mutation from a common ancestor., (© 2019 The Authors.)

Background: National genetic evaluations for disease resistance do not exist, precluding the genetic improvement of cattle for these traits. We imputed BovineHD genotypes to whole genome sequence for 2703 Holsteins that were cases or controls for Bovine Respiratory Disease and sampled from either California or New Mexico to construct and compare genomic prediction models. The sequence variation reference dataset comprised variants called for 1578 animals from Run 5 of the 1000 Bull Genomes Project, including 450 Holsteins and 29 animals sequenced from this study population. Genotypes for 9,282,726 variants with minor allele frequencies ≥5% were imputed and used to obtain genomic predictions in GEMMA using a Bayesian Sparse Linear Mixed Model., Results: Variation explained by markers increased from 13.6% using BovineHD data to 14.4% using imputed whole genome sequence data and the resolution of genomic regions detected as harbouring QTL substantially increased. Explained variation in the analysis of the combined California and New Mexico data was less than when data for each state were separately analysed and the estimated genetic correlation between risk of Bovine Respiratory Disease in California and New Mexico Holsteins was - 0.36. Consequently, genomic predictions trained using the data from one state did not accurately predict disease risk in the other state. To determine if a prediction model could be developed with utility in both states, we selected variants within genomic regions harbouring: 1) genes involved in the normal immune response to infection by pathogens responsible for Bovine Respiratory Disease detected by RNA-Seq analysis, and/or 2) QTL identified in the association analysis of the imputed sequence variants. The model based on QTL selected variants is biased but when trained in one state generated BRD risk predictions with positive accuracies in the other state., Conclusions: We demonstrate the utility of sequence-based and biology-driven model development for genomic selection. Disease phenotypes cannot be routinely recorded in most livestock species and the observed phenotypes may vary in their genomic architecture due to variation in the pathogen composition across environments. Elucidation of trait biology and genetic architecture may guide the development of prediction models with utility across breeds and environments.

Bovine respiratory disease (BRD) is the most common infectious disease of beef and dairy cattle and is characterized by a complex infectious etiology that includes a variety of viral and bacterial pathogens. We examined the global changes in mRNA abundance in healthy lung and lung lesions and in the lymphoid tissues bronchial lymph node, retropharyngeal lymph node, nasopharyngeal lymph node and pharyngeal tonsil collected at the peak of clinical disease from beef cattle experimentally challenged with either bovine respiratory syncytial virus, infectious bovine rhinotracheitis, bovine viral diarrhea virus, Mannheimia haemolytica or Mycoplasma bovis. We identified signatures of tissue-specific transcriptional responses indicative of tropism in the coordination of host's immune tissue responses to infection by viral or bacterial infections. Furthermore, our study shows that this tissue tropism in host transcriptional response to BRD pathogens results in the activation of different networks of response genes. The differential crosstalk among genes expressed in lymphoid tissues was predicted to be orchestrated by specific immune genes that act as 'key players' within expression networks. The results of this study serve as a basis for the development of innovative therapeutic strategies and for the selection of cattle with enhanced resistance to BRD.

Many genetic markers related to health or production traits are not evaluated in populations independent of the discovery population or related to phenotype. Here we evaluated 68 single nucleotide polymorphisms (SNP) in candidate genes previously associated with genetic merit for fertility and production traits for association with phenotypic measurements of fertility in a population of Holstein cows that was selected based on predicted transmitting ability (PTA) for daughter pregnancy rate (DPR; high, ≥1, n = 989; low, ≤ -1.0, n = 1,285). Cows with a high PTA for DPR had higher pregnancy rate at first service, fewer services per conception, and fewer days open than cows with a low PTA for DPR. Of the 68 SNP, 11 were associated with pregnancy rate at first service, 16 with services per conception, and 19 with days open. Single nucleotide polymorphisms in 12 genes (BDH2, BSP3, CAST, CD2, CD14, FUT1, FYB, GCNT3, HSD17B7, IBSP, OCLN, and PCCB) had significant associations with 2 fertility traits, and SNP in 4 genes (CSPP1, FCER1G, PMM2, and TBC1D24) had significant associations with each of the 3 traits. Results from this experiment were compared with results from 2 earlier studies in which the SNP were associated with genetic estimates of fertility. One study involved the same animals as used here, and the other study was of an independent population of bulls. A total of 13 SNP associated with 1 or more phenotypic estimates of fertility were directionally associated with genetic estimates of fertility in the same cow population. Moreover, 14 SNP associated with reproductive phenotype were directionally associated with genetic estimates of fertility in the bull population. Nine SNP (located in BCAS, BSP3, CAST, FUT1, HSD17B7, OCLN, PCCB, PMM2, and TBC1D24) had a directional association with fertility in all 3 studies. Examination of the function of the genes with SNP associated with reproduction in more than one study indicates the importance of steroid hormones and immune function as determinants of reproductive function. All but 1 of the 68 evaluated SNP were variable in 11 breeds besides Holstein, indicating the potential effects of these SNP on reproductive function across breeds of cattle., (Copyright © 2017 American Dairy Science Association. Published by Elsevier Inc. All rights reserved.)

Hereditary paroxysmal dyskinesias (PxD) are a heterogeneous group of movement disorders classified by frequency, duration, and triggers of the episodes. A young-adult onset canine PxD has segregated as an autosomal recessive trait in Soft-Coated Wheaten Terriers. The medical records and videos of episodes from 25 affected dogs were reviewed. The episodes of hyperkinesia and dystonia lasted from several minutes to several hours and could occur as often as >10/day. They were not associated with strenuous exercise or fasting but were sometimes triggered by excitement. The canine PxD phenotype most closely resembled paroxysmal non-kinesigenic dyskinesia (PNKD) of humans. Whole genome sequences were generated with DNA from 2 affected dogs and analyzed in comparison to 100 control canid whole genome sequences. The two whole genome sequences from dogs with PxD had a rare homozygous PIGN:c.398C > T transition, which predicted the substitution of an isoleucine for a highly conserved threonine in the encoded enzyme. All 25 PxD-affected dogs were PIGN:c.398T allele homozygotes, whereas there were no c.398T homozygotes among 1185 genotyped dogs without known histories of PxD. PIGN encodes an enzyme involved in the biosynthesis of glycosylphosphatidylinositol (GPI), which anchors a variety of proteins including CD59 to the cell surface. Flow cytometry of PIGN-knockout HEK239 cells expressing recombinant human PIGN with the c.398T variant showed reduced CD59 expression. Mutations in human PIGN have been associated with multiple congenital anomalies-hypotonia-seizures syndrome-1 (MCAHS1). Movement disorders can be a part of MCAHS1, but this is the first PxD associated with altered GPI anchor function., Competing Interests: The authors declare that they have no conflict of interest. Ethical approval All studies were approved by the University of Missouri Animal Care and Use Committee, and informed consent was obtained from the dogs’ owners.

The two living species of bison (European and American) are among the few terrestrial megafauna to have survived the late Pleistocene extinctions. Despite the extensive bovid fossil record in Eurasia, the evolutionary history of the European bison (or wisent, Bison bonasus) before the Holocene (<11.7 thousand years ago (kya)) remains a mystery. We use complete ancient mitochondrial genomes and genome-wide nuclear DNA surveys to reveal that the wisent is the product of hybridization between the extinct steppe bison (Bison priscus) and ancestors of modern cattle (aurochs, Bos primigenius) before 120 kya, and contains up to 10% aurochs genomic ancestry. Although undetected within the fossil record, ancestors of the wisent have alternated ecological dominance with steppe bison in association with major environmental shifts since at least 55 kya. Early cave artists recorded distinct morphological forms consistent with these replacement events, around the Last Glacial Maximum (LGM, ∼21-18 kya).

The diversity and population genetics of copy number variation (CNV) in domesticated animals are not well understood. In this study, we analysed 75 genomes of major taurine and indicine cattle breeds (including Angus, Brahman, Gir, Holstein, Jersey, Limousin, Nelore, and Romagnola), sequenced to 11-fold coverage to identify 1,853 non-redundant CNV regions. Supported by high validation rates in array comparative genomic hybridization (CGH) and qPCR experiments, these CNV regions accounted for 3.1% (87.5 Mb) of the cattle reference genome, representing a significant increase over previous estimates of the area of the genome that is copy number variable (∼2%). Further population genetics and evolutionary genomics analyses based on these CNVs revealed the population structures of the cattle taurine and indicine breeds and uncovered potential diversely selected CNVs near important functional genes, including AOX1, ASZ1, GAT, GLYAT, and KRTAP9-1 Additionally, 121 CNV gene regions were found to be either breed specific or differentially variable across breeds, such as RICTOR in dairy breeds and PNPLA3 in beef breeds. In contrast, clusters of the PRP and PAG genes were found to be duplicated in all sequenced animals, suggesting that subfunctionalization, neofunctionalization, or overdominance play roles in diversifying those fertility-related genes. These CNV results provide a new glimpse into the diverse selection histories of cattle breeds and a basis for correlating structural variation with complex traits in the future., (Published by Oxford University Press on behalf of Kazusa DNA Research Institute 2016. This work is written by (a) US Government employee(s) and is in the public domain in the US.)

An autosomal recessive disease of Black Russian Terriers was previously described as a juvenile-onset, laryngeal paralysis and polyneuropathy similar to Charcot Marie Tooth disease in humans. We found that in addition to an axonal neuropathy, affected dogs exhibit microphthalmia, cataracts, and miotic pupils. On histopathology, affected dogs exhibit a spongiform encephalopathy characterized by accumulations of abnormal, membrane-bound vacuoles of various sizes in neuronal cell bodies, axons and adrenal cells. DNA from an individual dog with this polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV) was used to generate a whole genome sequence which contained a homozygous RAB3GAP1:c.743delC mutation that was absent from 73 control canine whole genome sequences. An additional 12 Black Russian Terriers with POANV were RAB3GAP1:c.743delC homozygotes. DNA samples from 249 Black Russian Terriers with no known signs of POANV were either heterozygotes or homozygous for the reference allele. Mutations in human RAB3GAP1 cause Warburg micro syndrome (WARBM), a severe developmental disorder characterized by abnormalities of the eye, genitals and nervous system including a predominantly axonal peripheral neuropathy. RAB3GAP1 encodes the catalytic subunit of a GTPase activator protein and guanine exchange factor for Rab3 and Rab18 respectively. Rab proteins are involved in membrane trafficking in the endoplasmic reticulum, axonal transport, autophagy and synaptic transmission. The neuronal vacuolation and membranous inclusions and vacuoles in axons seen in this canine disorder likely reflect alterations of these processes. Thus, this canine disease could serve as a model for WARBM and provide insight into its pathogenesis and treatment., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

Population bottlenecks, inbreeding, and artificial selection can all, in principle, influence levels of deleterious genetic variation. However, the relative importance of each of these effects on genome-wide patterns of deleterious variation remains controversial. Domestic and wild canids offer a powerful system to address the role of these factors in influencing deleterious variation because their history is dominated by known bottlenecks and intense artificial selection. Here, we assess genome-wide patterns of deleterious variation in 90 whole-genome sequences from breed dogs, village dogs, and gray wolves. We find that the ratio of amino acid changing heterozygosity to silent heterozygosity is higher in dogs than in wolves and, on average, dogs have 2-3% higher genetic load than gray wolves. Multiple lines of evidence indicate this pattern is driven by less efficient natural selection due to bottlenecks associated with domestication and breed formation, rather than recent inbreeding. Further, we find regions of the genome implicated in selective sweeps are enriched for amino acid changing variants and Mendelian disease genes. To our knowledge, these results provide the first quantitative estimates of the increased burden of deleterious variants directly associated with domestication and have important implications for selective breeding programs and the conservation of rare and endangered species. Specifically, they highlight the costs associated with selective breeding and question the practice favoring the breeding of individuals that best fit breed standards. Our results also suggest that maintaining a large population size, rather than just avoiding inbreeding, is a critical factor for preventing the accumulation of deleterious variants.

Background: Next-generation sequencing projects commonly commence by aligning reads to a reference genome assembly. While improvements in alignment algorithms and computational hardware have greatly enhanced the efficiency and accuracy of alignments, a significant percentage of reads often remain unmapped., Results: We generated de novo assemblies of unmapped reads from the DNA and RNA sequencing of the Bos taurus reference individual and identified the closest matching sequence to each contig by alignment to the NCBI non-redundant nucleotide database using BLAST. As expected, many of these contigs represent vertebrate sequence that is absent, incomplete, or misassembled in the UMD3.1 reference assembly. However, numerous additional contigs represent invertebrate species. Most prominent were several species of Spirurid nematodes and a blood-borne parasite, Babesia bigemina. These species are either not present in the US or are not known to infect taurine cattle and the reference animal appears to have been host to unsequenced sister species., Conclusions: We demonstrate the importance of exploring unmapped reads to ascertain sequences that are either absent or misassembled in the reference assembly and for detecting sequences indicative of parasitic or commensal organisms.

Background: The neuronal ceroid lipofuscinoses are heritable lysosomal storage diseases characterized by progressive neurological impairment and the accumulation of autofluorescent storage granules in neurons and other cell types. Various forms of human neuronal ceroid lipofuscinosis have been attributed to mutations in at least 13 different genes. So far, mutations in the canine orthologs of 7 of these genes have been identified in DNA from dogs with neuronal ceroid lipofuscinosis. The identification of new causal mutations could lead to the establishment of canine models to investigate the pathogenesis of the corresponding human neuronal ceroid lipofuscinoses and to evaluate and optimize therapeutic interventions for these fatal human diseases., Case Presentation: We obtained blood and formalin-fixed paraffin-embedded brain sections from a rescue dog that was reported to be a young adult Chinese Crested. The dog was euthanized at approximately 19 months of age as a consequence of progressive neurological decline that included blindness, anxiety, and cognitive impairment. A diagnosis of neuronal ceroid lipofuscinosis was made based on neurological signs, magnetic resonance imaging of the brain, and fluorescence microscopic and electron microscopic examination of brain sections. We isolated DNA from the blood and used it to generate a whole genome sequence with 33-fold average coverage. Among the 7.2 million potential sequence variants revealed by aligning the sequence reads to the canine genome reference sequence was a homozygous single base pair deletion in the canine ortholog of one of 13 known human NCL genes: MFSD8:c.843delT. MFSD8:c.843delT is predicted to cause a frame shift and premature stop codon resulting in a truncated protein, MFSD8:p.F282Lfs13*, missing its 239 C-terminal amino acids. The MFSD8:c.843delT allele is absent from the whole genome sequences of 101 healthy canids or dogs with other diseases. The genotyping of archived DNA from 1478 Chinese Cresteds did not identify any additional MFSD8:c.843delT homozygotes and found only one heterozygote., Conclusion: We conclude that the neurodegenerative disease of the Chinese Crested rescue dog was neuronal ceroid lipofuscinosis and that homozygosity for the MFSD8:c.843delT sequence variant was very likely to be the molecular-genetic cause of the disease.

The neuronal ceroid lipofuscinoses (NCLs) are hereditary neurodegenerative diseases characterized by seizures and progressive cognitive decline, motor impairment, and vision loss accompanied by accumulation of autofluorescent lysosomal storage bodies in the central nervous system and elsewhere in the body. Mutations in at least 14 genes underlie the various forms of NCL. One of these genes, CLN8, encodes an intrinsic membrane protein of unknown function that appears to be localized primarily to the endoplasmic reticulum. Most CLN8 mutations in people result in a form of NCL with a late infantile onset and relatively rapid progression. A mixed breed dog with Australian Shepherd and Blue Heeler ancestry developed neurological signs characteristic of NCL starting at about 8months of age. The signs became progressively worse and the dog was euthanized at 21months of age due to seizures of increasing frequency and severity. Postmortem examination of the brain and retinas identified massive accumulations of intracellular autofluorescent inclusions characteristic of the NCLs. Whole genome sequencing of DNA from this dog identified a CLN8:c.585G>A transition that predicts a CLN8:p.Trp195* nonsense mutation. This mutation appears to be rare in both ancestral breeds. All of our 133 archived DNA samples from Blue Heelers, and 1481 of our 1488 archived Australian Shepherd DNA samples tested homozygous for the reference CLN8:c.585G allele. Four of the Australian Shepherd samples tested heterozygous and 3 tested homozygous for the mutant CLN8:c.585A allele. All 3 dogs homozygous for the A allele exhibited clinical signs of NCL and in 2 of them NCL was confirmed by postmortem evaluation of brain tissue. The occurrence of confirmed NCL in 3 of 4 CLN8:c.585A homozygous dogs, plus the occurrence of clinical signs consistent with NCL in the fourth homozygote strongly suggests that this rare truncating mutation causes NCL. Identification of this NCL-causing mutation provides the opportunity for identifying dogs that can be used to establish a canine model for the CLN8 disease (also known as late infantile variant or late infantile CLN8 disease)., (Copyright © 2014. Published by Elsevier Inc.)

Previous archeological and genetic research has shown that modern cattle breeds are descended from multiple independent domestication events of the wild aurochs (Bos primigenius) ∼10,000 y ago. Two primary areas of domestication in the Middle East/Europe and the Indian subcontinent resulted in taurine and indicine lines of cattle, respectively. American descendants of cattle brought by European explorers to the New World beginning in 1493 generally have been considered to belong to the taurine lineage. Our analyses of 47,506 single nucleotide polymorphisms show that these New World cattle breeds, as well as many related breeds of cattle in southern Europe, actually exhibit ancestry from both the taurine and indicine lineages. In this study, we show that, although European cattle are largely descended from the taurine lineage, gene flow from African cattle (partially of indicine origin) contributed substantial genomic components to both southern European cattle breeds and their New World descendants. New World cattle breeds, such as Texas Longhorns, provide an opportunity to study global population structure and domestication in cattle. Following their introduction into the Americas in the late 1400s, semiferal herds of cattle underwent between 80 and 200 generations of predominantly natural selection, as opposed to the human-mediated artificial selection of Old World breeding programs. Our analyses of global cattle breed population history show that the hybrid ancestry of New World breeds contributed genetic variation that likely facilitated the adaptation of these breeds to a novel environment.

We report a systematic study of gene expression during myogenesis and transdifferentiation in four bovine muscle tissues and of adipogenesis in three bovine fat tissues using DNA microarray analysis. One hundred hybridizations were performed and 7245 genes of known and unknown function were identified as being differentially expressed. Supervised hierarchical cluster analysis of gene expression patterns revealed the tissue specificity of genes. A close relationship in global gene expression observed for adipocyte-like cells derived from muscle and adipocytes derived from intramuscular fat suggests a common origin for these cells. The role of transthyretin in myogenesis is a novel finding. Different genes were highly induced during the transdifferentiation of myogenic satellite cells and in the adipogenesis of preadipocytes, indicating the involvement of different molecular mechanisms in these processes. Induction of CD36 and FABP4 expression in adipocyte-like cells and adipocytes may share a common pathway., (Copyright © 2012 Elsevier Inc. All rights reserved.)

A recessive, adult-onset neuronal ceroid-lipofuscinosis (NCL) occurs in Tibetan terriers. A genome-wide association study restricted this NCL locus to a 1.3Mb region of canine chromosome 2 which contains canine ATP13A2. NCL-affected dogs were homozygous for a single-base deletion in ATP13A2, predicted to produce a frameshift and premature termination codon. Homozygous truncating mutations in human ATP13A2 have been shown by others to cause Kufor-Rakeb syndrome (KRS), a rare neurodegenerative disease. These findings suggest that KRS is also an NCL, although analysis of KRS brain tissue will be needed to confirm this prediction. Generalized brain atrophy, behavioral changes, and cognitive decline occur in both people and dogs with ATP13A2 mutations; however, other clinical features differ between the species. For example, Tibetan terriers with NCL develop cerebellar ataxia not reported in KRS patients and KRS patients exhibit parkinsonism and pyramidal dysfunction not observed in affected Tibetan terriers. To see if ATP13A2 mutations could be responsible for some cases of human adult-onset NCL (Kufs disease), we resequenced ATP13A2 from 28 Kufs disease patients. None of these patients had ATP13A2 sequence variants likely to be causal for their disease, suggesting that mutations in this gene are not common causes of Kufs disease., (Copyright © 2011 Elsevier Inc. All rights reserved.)

Purpose: To identify the genetic cause of isolated canine ectopia lentis, a well-characterized veterinary disease commonly referred to as primary lens luxation (PLL) and to compare the canine disease with a newly described human Weill-Marchesani syndrome (WMS)-like disease of similar genetic etiology., Methods: Genomewide association analysis and fine mapping by homozygosity were used to identify the chromosomal segment harboring the PLL locus. The resequencing of a regional candidate gene was used to discover a mutation in a splice donor site predicted to cause exon skipping. Exon skipping was confirmed by reverse transcription-polymerase chain reaction amplification of RNA isolated from PLL-affected eyes and from skin fibroblast cultures from PLL-affected dogs. An allelic discrimination assay was used to genotype individual dogs at the splice donor site mutation., Results: The PLL locus was mapped to a 664-kb region of canine chromosome 3 containing regional candidate gene ADAMTS17. Resequencing ADAMTS17 revealed a GT-->AT splice-donor-site mutation at the 5' end of intron 10. The predicted exon 10 skipping and resultant frame shift were confirmed with RNA derived from PLL-affected dogs. The ADAMTS17 mutation was significantly associated with clinical PLL in three different dog breeds., Conclusions: A truncating mutation in canine ADAMTS17 causes PLL, a well-characterized veterinary disease, which can now be compared to a recently described rare WMS-like disease caused by truncating mutations of the human ADAMTS17 ortholog.

The Toll-like receptor (TLR) and peptidoglycan recognition protein 1 (PGLYRP1) genes play key roles in the innate immune systems of mammals. While the TLRs recognize a variety of invading pathogens and induce innate immune responses, PGLYRP1 is directly microbicidal. We used custom allele-specific assays to genotype and validate 220 diallelic variants, including 54 nonsynonymous SNPs in 11 bovine innate immune genes (TLR1-TLR10, PGLYRP1) for 37 cattle breeds. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and we were unable to differentiate between the specialized B. t. taurus beef and dairy breeds, despite an average polymorphism density of one locus per 219 bp. Ninety-nine tagSNPs and one tag insertion-deletion polymorphism were sufficient to predict 100% of the variation at all 11 innate immune loci in both subspecies and their hybrids, whereas 58 tagSNPs captured 100% of the variation at 172 loci in B. t. taurus. PolyPhen and SIFT analyses of nonsynonymous SNPs encoding amino acid replacements indicated that the majority of these substitutions were benign, but up to 31% were expected to potentially impact protein function. Several diversity-based tests provided support for strong purifying selection acting on TLR10 in B. t. taurus cattle. These results will broadly impact efforts related to bovine translational genomics.

The Pecorans (higher ruminants) are believed to have rapidly speciated in the Mid-Eocene, resulting in five distinct extant families: Antilocapridae, Giraffidae, Moschidae, Cervidae, and Bovidae. Due to the rapid radiation, the Pecoran phylogeny has proven difficult to resolve, and 11 of the 15 possible rooted phylogenies describing ancestral relationships among the Antilocapridae, Giraffidae, Cervidae, and Bovidae have each been argued as representations of the true phylogeny. Here we demonstrate that a genome-wide single nucleotide polymorphism (SNP) genotyping platform designed for one species can be used to genotype ancient DNA from an extinct species and DNA from species diverged up to 29 million years ago and that the produced genotypes can be used to resolve the phylogeny for this rapidly radiated infraorder. We used a high-throughput assay with 54,693 SNP loci developed for Bos taurus taurus to rapidly genotype 678 individuals representing 61 Pecoran species. We produced a highly resolved phylogeny for this diverse group based upon 40,843 genome-wide SNP, which is five times as many informative characters as have previously been analyzed. We also establish a method to amplify and screen genomic information from extinct species, and place Bison priscus within the Bovidae. The quality of genotype calls and the placement of samples within a well-supported phylogeny may provide an important test for validating the fidelity and integrity of ancient samples. Finally, we constructed a phylogenomic network to accurately describe the relationships between 48 cattle breeds and facilitate inferences concerning the history of domestication and breed formation.

The imprints of domestication and breed development on the genomes of livestock likely differ from those of companion animals. A deep draft sequence assembly of shotgun reads from a single Hereford female and comparative sequences sampled from six additional breeds were used to develop probes to interrogate 37,470 single-nucleotide polymorphisms (SNPs) in 497 cattle from 19 geographically and biologically diverse breeds. These data show that cattle have undergone a rapid recent decrease in effective population size from a very large ancestral population, possibly due to bottlenecks associated with domestication, selection, and breed formation. Domestication and artificial selection appear to have left detectable signatures of selection within the cattle genome, yet the current levels of diversity within breeds are at least as great as exists within humans.

Neonatal encephalopathy with seizures (NEWS) is a previously undescribed autosomal recessive disease of standard poodle puppies. Affected puppies are small and weak at birth. Many die in their first week of life. Those surviving past 1 week develop ataxia, a whole-body tremor, and, by 4 to 6 weeks of age, severe generalized clonic-tonic seizures. None have survived to 7 weeks of age. Cerebella from affected puppies were reduced in size and often contained dysplastic foci consisting of clusters of intermixed granule and Purkinje neurons. We used deoxyribonucleic acid samples from related standard poodles to map the NEWS locus to a 2.87-Mb segment of CFA36, which contains the canine ortholog of ATF2. This gene encodes activating transcription factor 2 (ATF-2), which participates in the cellular responses to a wide variety of stimuli. We amplified and sequenced all coding regions of canine ATF2 from a NEWS-affected puppy and identified a T > G transversion that predicts a methionine-to-arginine missense mutation at amino acid position 51. Methionine-51 lies within a hydrophobic docking site for mitogen-activated protein kinases that activate ATF-2 so the arginine substitution is likely to interfere with ATF-2 activation. All 20 NEWS-affected puppies in the standard poodle family were homozygous for the mutant G allele. The 58 clinically normal family members were either G/T heterozygotes or homozygous for the ancestral T allele. There are no previous reports of spontaneous ATF2 mutations in people or animals; however, atf2-knockout mice have cerebellar lesions that are similar to those in puppies with NEWS.

Bovine chromosome six (BTA6) harbors up to six quantitative trait loci (QTL) influencing the milk production of dairy cattle. In stark contrast to human, there is long-range linkage disequilibrium in dairy cattle, which has previously made it difficult to identify the mutations underlying these QTL. Using 38 microsatellite markers in a pedigree of 3,147 Holstein bulls, we fine mapped regions of BTA6 that had previously been shown to harbor QTL. Next, we sequenced a 12.3-kb region harboring Osteopontin, a positional candidate for the statistically most significant of the identified QTL. Nine mutations were identified, and only genotypes for the OPN3907 indel were concordant with the QTL genotypes of eight bulls that were established by segregation analysis. Four of these mutations were genotyped, and a joint linkage/linkage disequilibrium mapping analysis was used to demonstrate the existence of only two functionally distinct clusters of haplotypes within the QTL region, which were uniquely defined by OPN3907 alleles. We estimate a probability of 0.40 that no other mutation within this region is concordant with the QTL genotypes of these eight bulls. Finally, we demonstrate that the motif harboring OPN3907, which is upstream of the promoter and within a region known to harbor tissue-specific osteopontin regulatory elements, is moderately conserved among mammals. The motif was not retrieved from database queries and may be a novel regulatory element.