A recurrent de novo missense mutation in COL1A1 causes osteogenesis imperfecta type II and preterm delivery in Normande cattle.

BACKGROUND: Nine male and eight female calves born to a Normande artificial insemination bull named "Ly" were referred to the French National Observatory of Bovine Abnormalities for multiple fractures, shortened gestation, and stillbirth or perinatal mortality. RESULTS: Using Illumina BovineSNP50 array genotypes from affected calves and 84 half-sib controls, the associated locus was mapped to a 6.5-Mb interval on chromosome 19, assuming autosomal inheritance with germline mosaicism. Subsequent comparison of the whole-genome sequences of one case and 5116 control genomes, followed by genotyping in the affected pedigree, identified a de novo missense substitution within the NC1 domain of the COL1A1 gene (Chr19 g.36,473,965G > A; p.D1412N) as unique candidate variant. Interestingly, the affected residue was completely conserved among 243 vertebrate orthologs, and the same substitution in humans has been reported to cause type II osteogenesis imperfecta (OI), a connective tissue disorder that is characterized primarily by bone deformity and fragility. Moreover, three COL1A1 mutations have been described to cause the same syndrome in cattle. Necropsy, computed tomography, radiology, and histology confirmed the diagnosis of type II OI, further supporting the causality of this variant. In addition, a detailed analysis of gestation length and perinatal mortality in 1387 offspring of Ly and more than 160,000 progeny of 63 control bulls allowed us to statistically confirm in a large pedigree the association between type II OI and preterm delivery, which is probably due to premature rupture of fetal membranes and has been reported in several isolated cases of type II OI in humans and cattle. Finally, analysis of perinatal mortality rates and segregation distortion supported a low level of germ cell mosaicism in Ly, with an estimate of 4.5% to 7.7% of mutant sperm and thus 63 to 107 affected calves born. These numbers contrast with the 17 cases reported and raise concerns about the underreporting of congenital defects to heredo-surveillance platforms, even for textbook genetic syndromes. CONCLUSIONS: In conclusion, we describe a large animal model for a recurrent substitution in COL1A1 that is responsible for type II OI in humans. More generally, this study highlights the utility of such datasets and large half-sib families available in livestock species to characterize sporadic genetic defects.

Effect of temperature-humidity index on the evolution of trade-offs between fertility and production in dairy cattle.

BACKGROUND: In the current context of climate change, livestock production faces many challenges to improve the sustainability of systems. Dairy farming, in particular, must find ways to select animals that will be able to achieve sufficient overall production while maintaining their reproductive ability in environments with increasing temperatures. With future forecasted climate conditions in mind, this study used data from Holstein and Montbeliarde dairy cattle to: (1) estimate the genetic-by-temperature-humidity index (THI) interactions for female fertility, and (2) evaluate the production-fertility trade-off with increasing values of THI. RESULTS: Two-trait random regression models were fitted for conception rate (fertility) and test-day protein yield (production). For fertility, genetic correlations between different THI values were generally above 0.75, suggesting weak genotype-by-THI interactions for conception rate in both breeds. However, the genetic correlations between the conception rate breeding values at the current average THI (THI = 50, corresponding to a 24-h average temperature of 8 °C at 50% relative humidity) and their slopes (i.e., potential reranking) for heat stress scenarios (THI > 70), were different for each breed. For Montbeliarde, this correlation tended to be positive (i.e., overall the best reproducers are less affected by heat stress), whereas for Holstein it was approximately zero. Finally, our results indicated a weak antagonism between production and fertility, although for Montbeliarde this antagonism intensified with increasing THI. CONCLUSIONS: Within the range of weather conditions studied, increasing temperatures are not expected to exacerbate the fertility-production trade-off. However, our results indicated that the animals with the best breeding values for production today will be the most affected by temperature increases, both in terms of fertility and production. Nonetheless, these animals should remain among the most productive ones during heat waves. For Montbeliarde, the current selection program for fertility seems to be adequate for ensuring the adaptation of fertility traits to temperature increases, without adverse effects on production. Such a conclusion cannot be drawn for Holstein. In the future, the incorporation of a heat tolerance index into dairy cattle breeding programs would be valuable to promote the selection of animals adapted to future climate conditions.

Expected values for the accuracy of predicted breeding values accounting for genetic differences between reference and target populations.

BACKGROUND: Genetic merit, or breeding values as referred to in livestock and crop breeding programs, is one of the keys to the successful selection of animals in commercial farming systems. The developments in statistical methods during the twentieth century and single nucleotide polymorphism (SNP) chip technologies in the twenty-first century have revolutionized agricultural production, by allowing highly accurate predictions of breeding values for selection candidates at a very early age. Nonetheless, for many breeding populations, realized accuracies of predicted breeding values (PBV) remain below the theoretical maximum, even when the reference population is sufficiently large, and SNPs included in the model are in sufficient linkage disequilibrium (LD) with the quantitative trait locus (QTL). This is particularly noticeable over generations, as we observe the so-called erosion of the effects of SNPs due to recombinations, accompanied by the erosion of the accuracy of prediction. While accurately quantifying the erosion at the individual SNP level is a difficult and unresolved task, quantifying the erosion of the accuracy of prediction is a more tractable problem. In this paper, we describe a method that uses the relationship between reference and target populations to calculate expected values for the accuracies of predicted breeding values for non-phenotyped individuals accounting for erosion. The accuracy of the expected values was evaluated through simulations, and a further evaluation was performed on real data. RESULTS: Using simulations, we empirically confirmed that our expected values for the accuracy of PBV accounting for erosion were able to correctly determine the prediction accuracy of breeding values for non-phenotyped individuals. When comparing the expected to the realized accuracies of PBV with real data, only one out of the four traits evaluated presented accuracies that were significantly higher than the expected, approaching h 2 . CONCLUSIONS: We defined an index of genetic correlation between reference and target populations, which summarizes the expected overall erosion due to differences in allele frequencies and LD patterns between populations. We used this correlation along with a trait's heritability to derive expected values for the accuracy ( R ) of PBV accounting for the erosion, and demonstrated that our derived E R | erosion is a reliable metric.

Sequence-based GWAS meta-analyses for beef production traits.

BACKGROUND: Combining the results of within-population genome-wide association studies (GWAS) based on whole-genome sequences into a single meta-analysis (MA) is an accurate and powerful method for identifying variants associated with complex traits. As part of the H2020 BovReg project, we performed sequence-level MA for beef production traits. Five partners from France, Switzerland, Germany, and Canada contributed summary statistics from sequence-based GWAS conducted with 54,782 animals from 15 purebred or crossbred populations. We combined the summary statistics for four growth, nine morphology, and 15 carcass traits into 16 MA, using both fixed effects and z-score methods. RESULTS: The fixed-effects method was generally more informative to provide indication on potentially causal variants, although we combined substantially different traits in each MA. In comparison with within-population GWAS, this approach highlighted (i) a larger number of quantitative trait loci (QTL), (ii) QTL more frequently located in genomic regions known for their effects on growth and meat/carcass traits, (iii) a smaller number of genomic variants within the QTL, and (iv) candidate variants that were more frequently located in genes. MA pinpointed variants in genes, including MSTN, LCORL, and PLAG1 that have been previously associated with morphology and carcass traits. We also identified dozens of other variants located in genes associated with growth and carcass traits, or with a function that may be related to meat production (e.g., HS6ST1, HERC2, WDR75, COL3A1, SLIT2, MED28, and ANKAR). Some of these variants overlapped with expression or splicing QTL reported in the cattle Genotype-Tissue Expression atlas (CattleGTEx) and could therefore regulate gene expression. CONCLUSIONS: By identifying candidate genes and potential causal variants associated with beef production traits in cattle, MA demonstrates great potential for investigating the biological mechanisms underlying these traits. As a complement to within-population GWAS, this approach can provide deeper insights into the genetic architecture of complex traits in beef cattle.

Large-scale detection and characterization of interchromosomal rearrangements in normozoospermic bulls using massive genotype and phenotype data sets.

In this paper, we developed a highly sensitive approach to detect interchromosomal rearrangements in cattle by searching for abnormal linkage disequilibrium patterns between markers located on different chromosomes in large paternal half-sib families genotyped as part of routine genomic evaluations. We screened 5571 families of artificial insemination sires from 15 breeds and revealed 13 putative interchromosomal rearrangements, 12 of which were validated by cytogenetic analysis and long-read sequencing. These consisted of one Robertsonian fusion, 10 reciprocal translocations, and the first case of insertional translocation reported in cattle. Taking advantage of the wealth of data available in cattle, we performed a series of complementary analyses to define the exact nature of these rearrangements, investigate their origins, and search for factors that may have favored their occurrence. We also evaluated the risks to the livestock industry and showed significant negative effects on several traits in the sires and in their balanced or aneuploid progeny compared with wild-type controls. Thus, we present the most comprehensive and thorough screen for interchromosomal rearrangements compatible with normal spermatogenesis in livestock species. This approach is readily applicable to any population that benefits from large genotype data sets, and will have direct applications in animal breeding. Finally, it also offers interesting prospects for basic research by allowing the detection of smaller and rarer types of chromosomal rearrangements than GTG banding, which are interesting models for studying gene regulation and the organization of genome structure.

X-linked genes influence various complex traits in dairy cattle.

BACKGROUND: The search for quantitative trait loci (QTL) affecting traits of interest in mammals is frequently limited to autosomes, with the X chromosome excluded because of its hemizygosity in males. This study aimed to assess the importance of the X chromosome in the genetic determinism of 11 complex traits related to milk production, milk composition, mastitis resistance, fertility, and stature in 236,496 cows from three major French dairy breeds (Holstein, Montbéliarde, and Normande) and three breeds of regional importance (Abondance, Tarentaise, and Vosgienne). RESULTS: Estimates of the proportions of heritability due to autosomes and X chromosome (h²X) were consistent among breeds. On average over the 11 traits, h²X=0.008 and the X chromosome explained ~ 3.5% of total genetic variance. GWAS was performed within-breed at the sequence level (~ 200,000 genetic variants) and then combined in a meta-analysis. QTL were identified for most breeds and traits analyzed, with the exception of Tarentaise and Vosgienne and two fertility traits. Overall, 3, 74, 59, and 71 QTL were identified in Abondance, Montbéliarde, Normande, and Holstein, respectively, and most were associated with the most-heritable traits (milk traits and stature). The meta-analyses, which assessed a total of 157 QTL for the different traits, highlighted new QTL and refined the positions of some QTL found in the within-breed analyses. Altogether, our analyses identified a number of functional candidate genes, with the most notable being GPC3, MBNL3, HS6ST2, and DMD for dairy traits; TMEM164, ACSL4, ENOX2, HTR2C, AMOT, and IRAK1 for udder health; MAMLD1 and COL4A6 for fertility; and NRK, ESX1, GPR50, GPC3, and GPC4 for stature. CONCLUSIONS: This study demonstrates the importance of the X chromosome in the genetic determinism of complex traits in dairy cattle and highlights new functional candidate genes and variants for these traits. These results could potentially be extended to other species as many X-linked genes are shared among mammals.

The Bovine Pangenome Consortium: democratizing production and accessibility of genome assemblies for global cattle breeds and other bovine species.

The Bovine Pangenome Consortium (BPC) is an international collaboration dedicated to the assembly of cattle genomes to develop a more complete representation of cattle genomic diversity. The goal of the BPC is to provide genome assemblies and a community-agreed pangenome representation to replace breed-specific reference assemblies for cattle genomics. The BPC invites partners sharing our vision to participate in the production of these assemblies and the development of a common, community-approved, pangenome reference as a public resource for the research community ( https://bovinepangenome.github.io/ ). This community-driven resource will provide the context for comparison between studies and the future foundation for cattle genomic selection.

Integrin alpha 6 homozygous splice-site mutation causes a new form of junctional epidermolysis bullosa in Charolais cattle.

BACKGROUND: Inherited epidermolysis bullosa (EB) is a group of painful and life-threatening genetic disorders that are characterized by mechanically induced blistering of the skin and mucous membranes. Congenital skin fragility resembling EB was recently reported in three Charolais calves born in two distinct herds from unaffected parents. Phenotypic and genetic analyses were carried out to describe this condition and its molecular etiology. RESULTS: Genealogical, pathological and histological investigations confirmed the diagnosis of recessive EB. However, the affected calves showed milder clinical signs compared to another form of EB, which was previously reported in the same breed and is caused by a homozygous deletion of the ITGB4 gene. Homozygosity mapping followed by analysis of the whole-genome sequences of two cases and 5031 control individuals enabled us to prioritize a splice donor site of ITGA6 (c.2160 + 1G > T; Chr2 g.24112740C > A) as the most compelling candidate variant. This substitution showed a perfect genotype-phenotype correlation in the two affected pedigrees and was found to segregate only in Charolais, and at a very low frequency (f = 1.6 × 10-4) after genotyping 186,154 animals from 15 breeds. Finally, RT-PCR analyses revealed increased retention of introns 14 and 15 of the ITGA6 gene in a heterozygous mutant cow compared with a matched control. The mutant mRNA is predicted to cause a frameshift (ITGA6 p.I657Mfs1) that affects the assembly of the integrin α6ß4 dimer and its correct anchoring to the cell membrane. This dimer is a key component of the hemidesmosome anchoring complex, which ensures the attachment of basal epithelial cells to the basal membrane. Based on these elements, we arrived at a diagnosis of junctional EB. CONCLUSIONS: We report a rare example of partial phenocopies observed in the same breed and due to mutations that affect two members of the same protein dimer, and provide the first evidence of an ITGA6 mutation that causes EB in livestock species.

A Bos taurus sequencing methods benchmark for assembly, haplotyping, and variant calling.

Inspired by the production of reference data sets in the Genome in a Bottle project, we sequenced one Charolais heifer with different technologies: Illumina paired-end, Oxford Nanopore, Pacific Biosciences (HiFi and CLR), 10X Genomics linked-reads, and Hi-C. In order to generate haplotypic assemblies, we also sequenced both parents with short reads. From these data, we built two haplotyped trio high quality reference genomes and a consensus assembly, using up-to-date software packages. The assemblies obtained using PacBio HiFi reaches a size of 3.2 Gb, which is significantly larger than the 2.7 Gb ARS-UCD1.2 reference. The BUSCO score of the consensus assembly reaches a completeness of 95.8%, among highly conserved mammal genes. We also identified 35,866 structural variants larger than 50 base pairs. This assembly is a contribution to the bovine pangenome for the "Charolais" breed. These datasets will prove to be useful resources enabling the community to gain additional insight on sequencing technologies for applications such as SNP, indel or structural variant calling, and de novo assembly.

Estimation of genotype by temperature-humidity index interactions on milk production and udder health traits in Montbeliarde cows.

BACKGROUND: Heat stress negatively influences cattle welfare, health and productivity. To cope with the forecasted increases in temperature and heat waves frequency, identifying high-producing animals that are tolerant to heat is of capital importance to maintain milk production. This study, based on the joint analysis of on-farm performance and weather data, had two objectives: (1) to determine the response in production performances (milk, fat and protein yields, fat and protein contents) and udder health (somatic cell score) to temperature-humidity index (THI) variations in Montbeliarde cows, and (2) to estimate the interactions between genotype and THI, to enable the identification of the most adapted animals for facing the expected increases in temperature. RESULTS: Test-day records from first and second lactations from 2016 to 2020 were associated with the average THI during the three days before the test-day record. In total, 446,717 test-day records from 55,650 cows in first lactation and 457,516 test-day records from 58,229 cows in second lactation were analysed. The optimal THI was below 55 (i.e. ~ 12-13 °C) for all traits. Individual responses to THI were estimated by random regression models, which also included individual responses to days in milk. Regardless of the stage of lactation, genetic correlations along the THI gradient were above 0.80, which suggests that genotype-by-THI interactions were weak for production and udder health traits. Nevertheless, a variability in the individual slope of decay could be highlighted at high THI. The genetic correlation between production level at moderate THI and the slope at high THI was negative, while for somatic cell score, it was positive, indicating that heat stress amplifies the susceptibility to mastitis. CONCLUSIONS: The optimal THI for French Montbeliarde cows is below 55 for production and udder health traits. Genetic-by-THI interactions are weak in French Montbeliarde cows for production and udder health traits, but not all animals react in the same way to high temperatures. Even if there is little room for improvement, using a heat tolerance index in cattle selection would be relevant to anticipate the expected increases in temperature. Further investigations are needed to interpret this variability on production traits. However, the current selection for mastitis resistance seems appropriate to adapt cattle to rising temperatures.

New insights into the genetic resistance to paratuberculosis in Holstein cattle via single-step genomic evaluation.

BACKGROUND: Bovine paratuberculosis, or Johne's disease (JD), is a contagious and incurable disease caused by Mycobacterium avium subsp. paratuberculosis (MAP). It has adverse effects on animal welfare and is very difficult to control, leading to serious economic consequences. An important line of defense to this disease is host genetic resistance to MAP, which, when it will be more fully understood, could be improved through selective breeding. Using a large dataset of Holstein cows (161,253 animals including 56,766 cows with ELISA serological phenotypes and 12,431 animals with genotypes), we applied a single-step single nucleotide polymorphism (SNP) best linear unbiased prediction approach to investigate the genetic determinism underlying resistance to this disease (heritability estimate and identification of relevant genomic regions) and estimated genetic trends, reliability, and relative risk factors associated with genomic predictions. RESULTS: Resistance to JD was moderately heritable (0.14) and 16 genomic regions were detected that accounted for at least 0.05% of the breeding values variance (GV) in resistance to JD, and were located on chromosomes 1, 3, 5, 6, 7, 19, 20, 21, 23, 25, and 27, with the highest percentage of variance explained by regions on chromosomes 23 (0.36% GV), 5 (0.22% GV), 1 (0.14% GV), and 3 (0.13% GV). When estimated for the whole chromosomes, the autosomes with the largest overall contributions were chromosomes 3 (5.3% GV), 10 (4.8%), 23 (4.7%), 1 (3.6%), 7 (3.4%), 5 (2.9%), 12 (2.5%), 11 (2.2%), and 13 (2%). We estimated a slightly favorable genetic trend in resistance to JD over the last two decades, which can be explained by a low positive genetic correlation between resistance to JD and total merit index (+ 0.06). Finally, in a validation population of 907 cows, relatively reliable genomic predictions (reliability = 0.55) were obtained, which allowed the identification of cows at high risk of infection. CONCLUSIONS: This study provides new insights into the genetic determinism of resistance to JD and shows that this trait can be predicted from SNP genotypes. It has led to the implementation of a single-step genomic evaluation that should rapidly become an effective tool for controlling paratuberculosis on French Holstein farms.

Mutation of the MYH3 gene causes recessive cleft palate in Limousine cattle.

BACKGROUND: The palate is a structure separating the oral and nasal cavities and its integrity is essential for feeding and breathing. The total or partial opening of the palate is called a cleft palate and is a common malformation in mammals with environmental or hereditary aetiologies. Generally, it compromises life expectancy in the absence of surgical repair. A new form of non-syndromic cleft palate arose recently in Limousine cattle, with animals referred to the French National Observatory of Bovine Abnormalities since 2012. Since the number of affected animals has increased steadily, this study was undertaken to identify the cause of this disease. RESULTS: Based on pedigree analysis, occurrence of cleft palate in Limousine cattle was concordant with an autosomal recessive mode of inheritance. Genotyping of 16 affected animals and homozygosity mapping led to the identification of a single disease-associated haplotype on Bos taurus chromosome (BTA)19. The genome of two affected animals was sequenced, and their sequences were compared to the ARS-UCD1.2 reference genome to identify variants. The likely causal variants were compared to the variant database of the 1000 bull genome project and two fully linked mutations in exon 24 of the MYH3 (myosin heavy chain) gene were detected: a 1-bp non-synonymous substitution (BTA19:g.29609623A>G) and a 11-bp frameshift deletion (BTA19:g.29609605-29609615del). These two mutations were specific to the Limousine breed, with an estimated allele frequency of 2.4% and are predicted to be deleterious. The frameshift leads to a premature termination codon. Accordingly, mRNA and protein analyses in muscles from wild-type and affected animals revealed a decrease in MYH3 expression in affected animals, probably due to mRNA decay, as well as an absence of the MYH3 protein in these animals. MYH3 is mostly expressed in muscles, including craniofacial muscles, during embryogenesis, and its absence may impair palate formation. CONCLUSIONS: We describe a new form of hereditary cleft palate in Limousine cattle. We identified two fully linked and deleterious mutations, ultimately leading to the loss-of-function of the MYH3 protein. The mutations were included on the Illumina EuroG10k v8 and EuroGMD v1 SNP chips and are used to set up a reliable eradication strategy in the French Limousine breed.

Standardized Whole Blood Assay and Bead-Based Cytokine Profiling Reveal Commonalities and Diversity of the Response to Bacteria and TLR Ligands in Cattle.

Recent developments in multiplex technologies enable the determination of a large nu\mber of soluble proteins such as cytokines in various biological samples. More than a one-by-one determination of the concentration of immune mediators, they permit the establishment of secretion profiles for a more accurate description of conditions related to infectious diseases or vaccination. Cytokine profiling has recently been made available for bovine species with the development of a Luminex® technology-based 15-plex assay. Independently from the manufacturer, we evaluated the bovine cytokine/chemokine multiplex assay for limits of detection, recovery rate, and reproducibility. Furthermore, we assessed cytokine secretion in blood samples from 107 cows upon stimulation with heat-killed bacteria and TLR2/4 ligands compared to a null condition. Secretion patterns were analyzed either using the absolute concentration of cytokines or using their relative concentration with respect to the overall secretion level induced by each stimulus. Using Partial Least Square-Discriminant Analysis, we show that the 15-cytokine profile is different under Escherichia coli, Staphylococcus aureus, and Streptococcus uberis conditions, and that IFN-γ, IL-1ß, and TNF-α contribute the most to differentiate these conditions. LPS and E. coli induced largely overlapping biological responses, but S. aureus and S. uberis were associated with distinct cytokine profiles than their respective TLR ligands. Finally, results based on adjusted or absolute cytokine levels yielded similar discriminative power, but led to different stimuli-related signatures.

Predicting male fertility from the sperm methylome: application to 120 bulls with hundreds of artificial insemination records.

BACKGROUND: Conflicting results regarding alterations to sperm DNA methylation in cases of spermatogenesis defects, male infertility and poor developmental outcomes have been reported in humans. Bulls used for artificial insemination represent a relevant model in this field, as the broad dissemination of bull semen considerably alleviates confounding factors and enables the precise assessment of male fertility. This study was therefore designed to assess the potential for sperm DNA methylation to predict bull fertility. RESULTS: A unique collection of 100 sperm samples was constituted by pooling 2-5 ejaculates per bull from 100 Montbéliarde bulls of comparable ages, assessed as fertile (n = 57) or subfertile (n = 43) based on non-return rates 56 days after insemination. The DNA methylation profiles of these samples were obtained using reduced representation bisulfite sequencing. After excluding putative sequence polymorphisms, 490 fertility-related differentially methylated cytosines (DMCs) were identified, most of which were hypermethylated in subfertile bulls. Interestingly, 46 genes targeted by DMCs are involved in embryonic and fetal development, sperm function and maturation, or have been related to fertility in genome-wide association studies; five of these were further analyzed by pyrosequencing. In order to evaluate the prognostic value of fertility-related DMCs, the sperm samples were split between training (n = 67) and testing (n = 33) sets. Using a Random Forest approach, a predictive model was built from the methylation values obtained on the training set. The predictive accuracy of this model was 72% on the testing set and 72% on individual ejaculates collected from an independent cohort of 20 bulls. CONCLUSION: This study, conducted on the largest set of bull sperm samples so far examined in epigenetic analyses, demonstrated that the sperm methylome is a valuable source of male fertility biomarkers. The next challenge is to combine these results with other data on the same sperm samples in order to improve the quality of the model and better understand the interplay between DNA methylation and other molecular features in the regulation of fertility. This research may have potential applications in human medicine, where infertility affects the interaction between a male and a female, thus making it difficult to isolate the male factor.

Genome-wide association study with imputed whole-genome sequence variants including large deletions for female fertility in 3 Nordic dairy cattle breeds.

Fertility is an economically important trait in livestock. Poor fertility in dairy cattle can be due to loss-of-function variants affecting any essential gene that causes early embryonic mortality in homozygotes. To identify fertility-associated quantitative trait loci, we performed single-marker association analyses for 8 fertility traits in Holstein, Jersey, and Nordic Red Dairy cattle using imputed whole-genome sequence variants including SNPs, indels, and large deletion. We then performed stepwise selection of independent markers from GWAS loci using conditional and joint association analyses. From single-marker analyses for fertility traits, we reported genome-wide significant associations of 30,384 SNPs, 178 indels, and 3 deletions in Holstein; 23,481 SNPs, 189 indels, and 13 deletions in Nordic Red; and 17 SNPs in Jersey cattle. Conditional and joint association analyses identified 37 and 23 independent associations in Holstein and Nordic Red Dairy cattle, respectively. Fertility-associated GWAS loci were enriched for developmental and cellular processes (Gene Ontology enrichment, false discovery rate < 0.05). For these quantitative trait loci regions (top marker and 500 kb of surrounding regions), we proposed several candidate genes with functional annotations corresponding to embryonic lethality and various fertility-related phenotypes in mouse and cattle. The inclusion of these top markers in future releases of the custom SNP chip used for genomic evaluations will enable their validation in independent populations and improve the accuracy of genomic predictions.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets.

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

Sequence-based GWAS and post-GWAS analyses reveal a key role of SLC37A1, ANKH, and regulatory regions on bovine milk mineral content.

The mineral composition of bovine milk plays an important role in determining its nutritional and cheese-making value. Concentrations of the main minerals predicted from mid-infrared spectra produced during milk recording, combined with cow genotypes, provide a unique opportunity to decipher the genetic determinism of these traits. The present study included 1 million test-day predictions of Ca, Mg, P, K, Na, and citrate content from 126,876 Montbéliarde cows, of which 19,586 had genotype data available. All investigated traits were highly heritable (0.50-0.58), with the exception of Na (0.32). A sequence-based genome-wide association study (GWAS) detected 50 QTL (18 affecting two to five traits) and positional candidate genes and variants, mostly located in non-coding sequences. In silico post-GWAS analyses highlighted 877 variants that could be regulatory SNPs altering transcription factor (TF) binding sites or located in non-coding RNA (mainly lncRNA). Furthermore, we found 47 positional candidate genes and 45 TFs highly expressed in mammary gland compared to 90 other bovine tissues. Among the mammary-specific genes, SLC37A1 and ANKH, encoding proteins involved in ion transport were located in the most significant QTL. This study therefore highlights a comprehensive set of functional candidate genes and variants that affect milk mineral content.

Confirmed effects of candidate variants for milk production, udder health, and udder morphology in dairy cattle.

BACKGROUND: Over the last years, genome-wide association studies (GWAS) based on imputed whole-genome sequences (WGS) have been used to detect quantitative trait loci (QTL) and highlight candidate genes for important traits. However, in general this approach does not allow to validate the effects of candidate mutations or determine if they are truly causative for the trait(s) in question. To address these questions, we applied a two-step, within-breed GWAS approach on 15 traits (5 linked with milk production, 2 with udder health, and 8 with udder morphology) in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) cattle. We detected the most-promising candidate variants (CV) using imputed WGS of 2515 MON, 2203 NOR, and 6321 HOL bulls, and validated their effects in three younger populations of 23,926 MON, 9400 NOR, and 51,977 HOL cows. RESULTS: Bull sequence-based GWAS detected 84 QTL: 13, 10, and 30 for milk production traits; 3, 0, and 2 for somatic cell score (SCS); and 8, 2 and 16 for udder morphology traits, in MON, NOR, and HOL respectively. Five genomic regions with effects on milk production traits were shared among the three breeds whereas six (2 for production and 4 for udder morphology and health traits) had effects in two breeds. In 80 of these QTL, 855 CV were highlighted based on the significance of their effects and functional annotation. The subsequent GWAS on MON, NOR, and HOL cows validated 8, 9, and 23 QTL for production traits; 0, 0, and 1 for SCS; and 4, 1, and 8 for udder morphology traits, respectively. In 47 of the 54 confirmed QTL, the CV identified in bulls had more significant effects than single nucleotide polymorphisms (SNPs) from the standard 50K chip. The best CV for each validated QTL was located in a gene that was functionally related to production (36 QTL) or udder (9 QTL) traits. CONCLUSIONS: Using this two-step GWAS approach, we identified and validated 54 QTL that included CV mostly located within functional candidate genes and explained up to 6.3% (udder traits) and 37% (production traits) of the genetic variance of economically important dairy traits. These CV are now included in the chip used to evaluate French dairy cattle and can be integrated into routine genomic evaluation.

Weighted single-step genomic best linear unbiased prediction integrating variants selected from sequencing data by association and bioinformatics analyses.

BACKGROUND: Sequencing data enable the detection of causal loci or single nucleotide polymorphisms (SNPs) highly linked to causal loci to improve genomic prediction. However, until now, studies on integrating such SNPs using a single-step genomic best linear unbiased prediction (ssGBLUP) model are scarce. We investigated the integration of sequencing SNPs selected by association (1262 SNPs) and bioinformatics (2359 SNPs) analyses into the currently used 54K-SNP chip, using three ssGBLUP models which make different assumptions on the distribution of SNP effects: a basic ssGBLUP model, a so-called featured ssGBLUP (ssFGBLUP) model that considered selected sequencing SNPs as a feature genetic component, and a weighted ssGBLUP (ssWGBLUP) model in which the genomic relationship matrix was weighted by the SNP variances estimated from a Bayesian whole-genome regression model, with every 1, 30, or 100 adjacent SNPs within a chromosome region sharing the same variance. We used data on milk production and female fertility in Danish Jersey. In total, 15,823 genotyped and 528,981‬ non-genotyped females born between 1990 and 2013 were used as reference population and 7415 genotyped females and 33,040 non-genotyped females born between 2014 and 2016 were used as validation population. RESULTS: With basic ssGBLUP, integrating SNPs selected from sequencing data improved prediction reliabilities for milk and protein yields, but resulted in limited or no improvement for fat yield and female fertility. Model performances depended on the SNP set used. When using ssWGBLUP with the 54K SNPs, reliabilities for milk and protein yields improved by 0.028 for genotyped animals and by 0.006 for non-genotyped animals compared with ssGBLUP. However, with the SNP set that included SNPs selected from sequencing data, no statistically significant difference in prediction reliability was observed between the three ssGBLUP models. CONCLUSIONS: In summary, when using 54K SNPs, a ssWGBLUP model with a common weight on the SNPs in a given region is a feasible approach for single-trait genetic evaluation. Integrating relevant SNPs selected from sequencing data into the standard SNP chip can improve the reliability of genomic prediction. Based on such SNP data, a basic ssGBLUP model was suggested since no significant improvement was observed from using alternative models such as ssWGBLUP and ssFGBLUP.