Dataset Information

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

ABSTRACT:

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.

SUBMITTER: Tribout T

PROVIDER: S-EPMC7529513 | biostudies-literature | 2020 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

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

Tribout Thierry T Croiseau Pascal P Lefebvre Rachel R Barbat Anne A Boussaha Mekki M Fritz Sébastien S Boichard Didier D Hoze Chris C Sanchez Marie-Pierre MP

Genetics, selection, evolution : GSE 20201001 1

<h4>Background</h4>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 ...[more]

PMID: 32998688

Similar Datasets

Project description:BackgroundGenome wide association studies (GWAS) in most cattle breeds result in large genomic intervals of significant associations making it difficult to identify causal mutations. This is due to the extensive, low-level linkage disequilibrium within a cattle breed. As there is less linkage disequilibrium across breeds, multibreed GWAS may improve precision of causal variant mapping. Here we test this hypothesis in a Holstein and Jersey cattle data set with 17,925 individuals with records for production and functional traits and 632,003 SNP markers.ResultsBy using a cross validation strategy within the Holstein and Jersey data sets, we were able to identify and confirm a large number of QTL. As expected, the precision of mapping these QTL within the breeds was limited. In the multibreed analysis, we found that many loci were not segregating in both breeds. This was partly an artefact of power of the experiments, with the number of QTL shared between the breeds generally increasing with trait heritability. False discovery rates suggest that the multibreed analysis was less powerful than between breed analyses, in terms of how much genetic variance was explained by the detected QTL. However, the multibreed analysis could more accurately pinpoint the location of the well-described mutations affecting milk production such as DGAT1. Further, the significant SNP in the multibreed analysis were significantly enriched in genes regions, to a considerably greater extent than was observed in the single breed analyses. In addition, we have refined QTL on BTA5 and BTA19 to very small intervals and identified a small number of potential candidate genes in these, as well as in a number of other regions.ConclusionWhere QTL are segregating across breed, multibreed GWAS can refine these to reasonably small genomic intervals. However, such QTL appear to represent only a fraction of the genetic variation. Our results suggest a significant proportion of QTL affecting milk production segregate within rather than across breeds, at least for Holstein and Jersey cattle.

Project description:The premises for the potential success of molecular breeding is the ability to identify major genes associated with important dairy related traits. The present study was taken up with the objectives to identify single nucleotide polymorphism (SNP) of bovine MASP2 and SIRT1 genes and its effect on estimated breeding values (EBVs) and to estimate genetic parameters for lactation milk yield (LMY), 305-day milk yield (305dMY), 305-day fat yield (305dFY), 305-day solid not fat yield (305dSNFY) and lactation length (LL) in Sahiwal dairy cattle to devise a promising improvement strategy. Genetic parameters and breeding values of milk production traits were estimated from 935 Sahiwal cattle population (1979-2019) reared at National Dairy Research Institute at Karnal, India. A total of 7 SNPs, where one SNP (g.499C>T) in exon 2 and four SNPs (g.576G>A, g.609T>C, g.684G>T and g.845A>G) in exon 3 region of MASP2 gene and 2 SNPs (g.-306T>C and g.-274G>C) in the promoter region of SIRT1 gene were identified in Sahiwal cattle population. Five of these identified SNPs were chosen for further genotyping by PCR-RFLP and association analysis. Association analysis was performed using estimated breeding values (n = 150) to test the effect of SNPs on LMY, 305dMY, 305dFY, 305dSNFY and LL. Association analysis revealed that, three SNP markers (g.499C>T, g.609T>C and g.-306T>C) were significantly associated with all milk yield traits. The estimates for heritability using repeatability model for LMY, 305dMY, 305dFY, 305dSNFY and LL were low, however the corresponding estimates from first parity were 0.20±0.08, 0.17±0.08, 0.13±0.09, 0.13±0.09 and 0.24, respectively. The repeatability estimates were moderate to high indicating consistency of performance over the parities and hence reliability of first lactation traits. Genetic correlations among the traits of first parity were high (0.55 to 0.99). From the results we could conclude that optimum strategy to improve the Sahiwal cattle further would be selecting the animals based on their first lactation 305dMY. Option top include the significant SNP in selection criteria can be explored. Taken together, a 2-stage selection approach, select Sahiwal animals early for the SNP and then on the basis of first lactation 305dMY will help to save resources.

Project description:An annual pattern of milk composition has been well recognized in dairy cattle, with the highest milk fat and protein concentration observed during the winter and lowest occurring in the summer; however, rhythms of milk yield and composition have not been well quantified. Cosinor rhythmometry is commonly used to model repeating daily and annual rhythms and allows determination of the amplitude (peak to mean), acrophase (time at peak), and period (time between peaks) of the rhythm. The objective of this study was to use cosinor rhythmometry to characterize the annual rhythms of milk yield and milk fat and protein concentration and yield using both national milk market and cow-level data. First, 10 yr of monthly average milk butterfat and protein concentration for each Federal Milk Marketing Order were obtained from the US Department of Agriculture Agricultural Marketing Service database. Fat and protein concentration fit a cosine function with a 12-mo period in all milk markets. We noted an interaction between milk marketing order and milk fat and protein concentration. The acrophase (time at peak) of the fat concentration rhythm ranged from December 4 to January 19 in all regions, whereas the rhythm of protein concentration peaked between December 27 and January 6. The amplitude (peak to mean) of the annual rhythm ranged from 0.07 to 0.14 percentage points for milk fat and from 0.08 to 0.12 percentage points for milk protein. The amplitude of the milk fat rhythm generally was lower in southern markets and higher in northern markets. Second, the annual rhythm of milk yield and milk fat and protein yield and concentration were analyzed in monthly test day data from 1,684 cows from 11 tiestall herds in Pennsylvania. Fat and protein concentration fit an annual rhythm in all herds, whereas milk and milk fat and protein yield only fit rhythms in 8 of the 11 herds. On average, milk yield peaked in April, fat and protein yield peaked in February, fat concentration peaked in January, and protein concentration peaked in December. Amplitudes of milk, fat, and protein yield averaged 0.82 kg, 55.3 g, and 30.4 g, respectively. Milk fat and protein concentration had average amplitudes of 0.12 and 0.07, respectively, similar to the results of the milk market data. Generally, milk yield and milk components fit annual rhythm regardless of parity or diacylglycerol O-acyltransferase 1 (DGAT1) K232A polymorphism, with only cows of the low-frequency AA genotype (5.2% of total cows) failing to fit rhythm of milk yield. In conclusion, the yearly rhythms of milk yield and fat and protein concentration and yield consistently occur regardless of region, herd, parity, or DGAT1 genotype and supports generation by a conserved endogenous annual rhythm.

Project description:BackgroundGenome-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 ≤ R2 ≤ 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.ResultsThirty-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).ConclusionsGWAS 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.

Dataset Information

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

Background

Results

Conclusions

Publications

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

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets