Explore the vast range of titles available.

GWAS conducted directly on imputed whole genome sequence have led to the identification of numerous genetic variants associated with agronomic traits in cattle. However, such variants are often simply markers in linkage disequilibrium with the actual causal variants, which is a limiting factor for the development of accurate genomic predictions. It is possible to identify causal variants by integrating information on how variants impact gene expression into GWAS output. RNA splicing plays a major role in regulating gene expression. Thus, assessing the effect of variants on RNA splicing may explain their function. Here, we use a high-throughput strategy to functionally analyse putative splice-disrupting variants in the bovine genome. Using GWAS, massively parallel reporter assay and deep learning algorithms designed to predict splice-disrupting variants, we identify 38 splice-disrupting variants associated with complex traits in cattle, three of which could be classified as causal. Our results indicate that splice-disrupting variants are widely found in the quantitative trait loci related to these phenotypes. Using our combined approach, we also assess the validity of splicing predictors originally developed to analyse human variants in the context of the bovine genome. The effects of GWAS-identified variants on splicing may explain their impact on gene expression. Here, the authors use GWAS, bioinformatic and functional analyses to highlight the role splicing variants play in complex traits in 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.

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.

Background Genome-wide association studies (GWAS) were performed at the sequence level to identify candidate mutations that affect the expression of six major milk proteins in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) dairy cattle. Whey protein (α-lactalbumin and β-lactoglobulin) and casein (αs1, αs2, β, and κ) contents were estimated by mid-infrared (MIR) spectrometry, with medium to high accuracy (0.59 ≤ R 2  ≤ 0.92), for 848,068 test-day milk samples from 156,660 cows in the first three lactations. Milk composition was evaluated as average test-day measurements adjusted for environmental effects. Next, we genotyped a subset of 8080 cows (2967 MON, 2737 NOR, and 2306 HOL) with the BovineSNP50 Beadchip. For each breed, genotypes were first imputed to high-density (HD) using HD single nucleotide polymorphisms (SNPs) genotypes of 522 MON, 546 NOR, and 776 HOL bulls. The resulting HD SNP genotypes were subsequently imputed to the sequence level using 27 million high-quality sequence variants selected from Run4 of the 1000 Bull Genomes consortium (1147 bulls). Within-breed, multi-breed, and conditional GWAS were performed. Results Thirty-four distinct genomic regions were identified. Three regions on chromosomes 6, 11, and 20 had very significant effects on milk composition and were shared across the three breeds. Other significant effects, which partially overlapped across breeds, were found on almost all the autosomes. Multi-breed analyses provided a larger number of significant genomic regions with smaller confidence intervals than within-breed analyses. Combinations of within-breed, multi-breed, and conditional analyses led to the identification of putative causative variants in several candidate genes that presented significant protein–protein interactions enrichment, including those with previously described effects on milk composition ( SLC37A1 , MGST1 , ABCG2 , CSN1S1 , CSN2 , CSN1S2 , CSN3 , PAEP , DGAT1 , AGPAT6 ) and those with effects reported for the first time here ( ALPL , ANKH , PICALM ). Conclusions GWAS applied to fine-scale phenotypes, multiple breeds, and whole-genome sequences seems to be effective to identify candidate gene variants. However, although we identified functional links between some candidate genes and milk phenotypes, the causality between candidate variants and milk protein composition remains to be demonstrated. Nevertheless, the identification of potential causative mutations that underlie milk protein composition may have immediate applications for improvements in cheese-making.

Background Due to their contribution to global warming, methane emissions from ruminants have been the subject of considerable scientific interest. It has been proposed that such emissions might be reduced using genetic selection; proposed phenotypes differ in the measurement methods used (direct or predicted methane emissions) and in the unit under consideration (g/d, g/kg of milk, g/kg of intake, residual methane emissions). Identifying the quantitative trait loci (QTLs) and candidate genes responsible for genetic variation in methane emissions allows a better understanding of the underlying genetic architecture of these phenotypes. Therefore, the aim of this study was to identify the genomic regions associated with six methane traits predicted from milk mid-infrared (MIR) spectra (0.33 ≤ R 2  ≤ 0.88) in French Holstein dairy cows using genome-wide association studies at the whole-genome-sequence level. Results Six methane emission traits—in g/d, in g/kg of fat- and protein-corrected milk, and in g/kg of dry matter intake—were predicted from milk MIR spectra routinely collected by French milk recording companies. A genome-wide association study of the predicted methane emissions of 40,609 primiparous Holstein cows was conducted using imputed whole-genome-sequence data. This analysis revealed 57 genomic regions of interest; between 1 and 8 QTLs were identified on each of the autosomes except 4, 12, 21, 24 and 26. We identified multiple genomic regions that were shared by two or more predicted methane traits, illustrating their common genetic basis. Functional annotation revealed potential candidate genes, in particular FASN , DGAT1 , ACSS2 , and KCNIP4 , which could be involved in biological pathways possibly related to methane production. Conclusions The methane traits studied here, which were predicted from milk MIR spectra, appear to be highly polygenic. Several genomic regions associated with these traits contain candidate genes previously associated with milk traits. Functional annotation and comparisons with studies using direct methane measurements support some potential candidate genes involved in biological pathways related to methane production. However, the overlap with genes influencing milk traits highlights the challenge of distinguishing whether these regions genuinely influence methane emissions or reflect the use of milk MIR spectra to predict the phenotypes.

Background Sequence-based genome-wide association studies (GWAS) provide high statistical power to identify candidate causal mutations when a large number of individuals with both sequence variant genotypes and phenotypes is available. A meta-analysis combines summary statistics from multiple GWAS and increases the power to detect trait-associated variants without requiring access to data at the individual level of the GWAS mapping cohorts. Because linkage disequilibrium between adjacent markers is conserved only over short distances across breeds, a multi-breed meta-analysis can improve mapping precision. Results To maximise the power to identify quantitative trait loci (QTL), we combined the results of nine within-population GWAS that used imputed sequence variant genotypes of 94,321 cattle from eight breeds, to perform a large-scale meta-analysis for fat and protein percentage in cattle. The meta-analysis detected (p ≤ 10 −8 ) 138 QTL for fat percentage and 176 QTL for protein percentage. This was more than the number of QTL detected in all within-population GWAS together (124 QTL for fat percentage and 104 QTL for protein percentage). Among all the lead variants, 100 QTL for fat percentage and 114 QTL for protein percentage had the same direction of effect in all within-population GWAS. This indicates either persistence of the linkage phase between the causal variant and the lead variant across breeds or that some of the lead variants might indeed be causal or tightly linked with causal variants. The percentage of intergenic variants was substantially lower for significant variants than for non-significant variants, and significant variants had mostly moderate to high minor allele frequencies. Significant variants were also clustered in genes that are known to be relevant for fat and protein percentages in milk. Conclusions Our study identified a large number of QTL associated with fat and protein percentage in dairy cattle. We demonstrated that large-scale multi-breed meta-analysis reveals more QTL at the nucleotide resolution than within-population GWAS. Significant variants were more often located in genic regions than non-significant variants and a large part of them was located in potentially regulatory regions.

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.

Background Linear models that are commonly used to predict breeding values in livestock species consider paternal influence solely as a genetic effect. However, emerging evidence in several species suggests the potential effect of non-genetic semen-mediated paternal effects on offspring phenotype. This study contributes to such research by analyzing the extent of non-genetic paternal effects on the performance of Holstein, Montbéliarde, and Normande dairy cows. Insemination data, including semen Batch Identifier (BI, a combination of bull identification and collection date), was associated with various traits measured in cows born from the insemination. These traits encompassed stature, milk production (milk, fat, and protein yields), udder health (somatic cell score and clinical mastitis), and female fertility (conception rates of heifers and cows). We estimated (1) the effects of age at collection and heat stress during spermatogenesis, and (2) the variance components associated with BI or Weekly aggregated BI (WBI). Results Overall, the non-genetic paternal effect estimates were small and of limited biological importance. However, while heat stress during spermatogenesis did not show significant associations with any of the traits studied in daughters, we observed significant effects of bull age at semen collection on the udder health of daughters. Indeed, cows born from bulls collected after 1500 days of age had higher somatic cell scores compared to those born from bulls collected at a younger age (less than 400 days old) in both Holstein and Normande breeds (+ 3% and + 5% of the phenotypic mean, respectively). In addition, across all breeds and traits analyzed, the estimates of non-genetic paternal variance were consistently low, representing on average 0.13% and 0.09% of the phenotypic variance for BI and WBI, respectively (ranging from 0 to 0.7%). These estimates did not significantly differ from zero, except for milk production traits (milk, fat, and protein yields) in the Holstein breed and protein yield in the Montbéliarde breed when WBI was considered. Conclusions Our findings indicate that non-genetic paternal information transmitted through semen does not substantially influence the offspring phenotype in dairy cattle breeds for routinely measured traits. This lack of substantial impact may be attributed to limited transmission or minimal exposure of elite bulls to adverse conditions.

Background Mastitis is a disease that incurs significant costs in the dairy industry. A promising approach to mitigate its negative effects is to genetically improve the resistance of dairy cattle to mastitis. A meta-analysis of genome-wide association studies (GWAS) across multiple breeds for clinical mastitis (CM) and its indicator trait, somatic cell score (SCS), is a powerful method to identify functional genetic variants that impact mastitis resistance. Results We conducted meta-analyses of eight and fourteen GWAS on CM and SCS, respectively, using 30,689 and 119,438 animals from six dairy cattle breeds. Methods for the meta-analyses were selected to properly account for the multi-breed structure of the GWAS data. Our study revealed 58 lead markers that were associated with mastitis incidence, including 16 loci that did not overlap with previously identified quantitative trait loci (QTL), as curated at the Animal QTLdb. Post-GWAS analysis techniques such as gene-based analysis and genomic feature enrichment analysis enabled prioritization of 31 candidate genes and 14 credible candidate causal variants that affect mastitis. Conclusions Our list of candidate genes can help to elucidate the genetic architecture underlying mastitis resistance and provide better tools for the prevention or treatment of mastitis, ultimately contributing to more sustainable animal production.

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.