Dataset Information

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

ABSTRACT: BACKGROUND:Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value predictions and followed by gene identification and in silico functional analyses. The association results were further validated using five scenarios with different numbers of SNPs. RESULTS:Seven hundred and eighty-two SNPs were significantly associated with calving performance at a genome-wise false discovery rate (FDR) of 5%. Most of these significant SNPs were on chromosomes 18 (71.9%), 17 (7.4%), 5 (6.8%) and 7 (2.4%) and mapped to 675 genes, among which 142 included at least one significant SNP and 532 were nearby one (100 kbp). For body conformation traits, 607 SNPs were significant at a genome-wise FDR of 5% and most of them were located on chromosomes 5 (30%), 18 (27%), 20 (13%), 6 (6%), 7 (5%), 14 (5%) and 13 (3%). SNP enrichment functional analyses for calving traits at a FDR of 1% suggested potential biological processes including musculoskeletal movement, meiotic cell cycle, oocyte maturation and skeletal muscle contraction. Furthermore, pathway analyses suggested potential pathways associated with calving performance traits including tight junction, oxytocin signaling, and MAPK signaling (P < 0.10). The prediction ability of the 1206 significant SNPs was between 78 and 83% of the prediction ability of the BovineSNP50 SNPs for calving performance traits and between 35 and 79% for body conformation traits. CONCLUSIONS:Various SNPs that are significantly associated with calving performance are located within or nearby genes with potential roles in tight junction, oxytocin signaling, and MAPK signaling. Combining the significant SNPs or SNPs within or nearby gene(s) from the HD panel with the BovineSNP50 panel yielded a marginal increase in the accuracy of prediction of genomic estimated breeding values for all traits compared to the use of the BovineSNP50 panel alone.

SUBMITTER: Abo-Ismail MK

PROVIDER: S-EPMC6389134 | biostudies-literature | 2017 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

Abo-Ismail Mohammed K MK Brito Luiz F LF Miller Stephen P SP Sargolzaei Mehdi M Grossi Daniela A DA Moore Steve S SS Plastow Graham G Stothard Paul P Nayeri Shadi S Schenkel Flavio S FS

Genetics, selection, evolution : GSE 20171107 1

<h4>Background</h4>Our aim was to identify genomic regions via genome-wide association studies (GWAS) to improve the predictability of genetic merit in Holsteins for 10 calving and 28 body conformation traits. Animals were genotyped using the Illumina Bovine 50 K BeadChip and imputed to the Illumina BovineHD BeadChip (HD). GWAS were performed on 601,717 real and imputed single nucleotide polymorphism (SNP) genotypes using a single-SNP mixed linear model on 4841 Holstein bulls with breeding value ...[more]

PMID: 29115939

Similar Datasets

Project description:Simple Summary Udder conformation traits are economically important for dairy animals in the dairy industry. Milk production loss can be reduced via a better udder structure in dairy cattle. Udder traits are related to milk production and the somatic cell count, which is a sign of mammary infections (mastitis); as such, it is vital to understand the genetic architecture underlying udder traits in Holstein Friesian cattle for genetic development and long-term selection. Through a GWAS on udder structure traits in Chinese Holstein cattle, we identified several significant single-nucleotide polymorphisms (SNPs) and candidate genes associated with udder traits. The results could provide useful information regarding the genetics architecture of udder structure traits, thus improving the genetic information, health, longevity, and production of dairy cattle. Abstract Udder conformation traits are one of the most economic traits in dairy cows, greatly affecting animal health, milk production, and producer profitability in the dairy industry. Genetic analysis of udder structure and scores have been developed in Holstein cattle. In our research, we conducted a genome-wide association study for five udder traits, including anterior udder attachment (AUA), central suspensory ligament (CSL), posterior udder attachment height (PUAH), posterior udder attachment width (PUAW), and udder depth (UD), in which the fixed and random model circulating probability unification (FarmCPU) model was applied for the association analysis. The heritability and the standard errors of these five udder traits ranged from 0.04 ± 0.00 to 0.49 ± 0.03. Phenotype data were measured from 1000 Holstein cows, and the GeneSeek Genomic Profiler (GGP) Bovine 100 K SNP chip was used to analyze genotypic data in Holstein cattle. For GWAS analysis, 984 individual cows and 84,407 single-nucleotide polymorphisms (SNPs) remained after quality control; a total of 18 SNPs were found at the GW significant threshold (p < 5.90 × 10−7). Many candidate genes were identified within 200kb upstream or downstream of the significant SNPs, which include MGST1, MGST2, MTUS1, PRKN, STXBP6, GRID2, E2F8, CDH11, FOXP1, SLF1, TMEM117, SBF2, GC, ADGRB3, and GCLC. Pathway analysis revealed that 58 Gene Ontology (GO) terms and 18 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched with adjusted p values, and these GO terms and the KEGG pathway analysis were associated with biological information, metabolism, hormonal growth, and development processes. These results could give valuable biological information for the genetic architecture of udder conformation traits in dairy Holstein cattle.

Project description:BackgroundThe impact of additive-genetic relationships captured by single nucleotide polymorphisms (SNPs) on the accuracy of genomic breeding values (GEBVs) has been demonstrated, but recent studies on data obtained from Holstein populations have ignored this fact. However, this impact and the accuracy of GEBVs due to linkage disequilibrium (LD), which is fairly persistent over generations, must be known to implement future breeding programs.Materials and methodsThe data set used to investigate these questions consisted of 3,863 German Holstein bulls genotyped for 54,001 SNPs, their pedigree and daughter yield deviations for milk yield, fat yield, protein yield and somatic cell score. A cross-validation methodology was applied, where the maximum additive-genetic relationship (amax) between bulls in training and validation was controlled. GEBVs were estimated by a Bayesian model averaging approach (BayesB) and an animal model using the genomic relationship matrix (G-BLUP). The accuracy of GEBVs due to LD was estimated by a regression approach using accuracy of GEBVs and accuracy of pedigree-based BLUP-EBVs.ResultsAccuracy of GEBVs obtained by both BayesB and G-BLUP decreased with decreasing amax for all traits analyzed. The decay of accuracy tended to be larger for G-BLUP and with smaller training size. Differences between BayesB and G-BLUP became evident for the accuracy due to LD, where BayesB clearly outperformed G-BLUP with increasing training size.ConclusionsGEBV accuracy of current selection candidates varies due to different additive-genetic relationships relative to the training data. Accuracy of future candidates can be lower than reported in previous studies because information from close relatives will not be available when selection on GEBVs is applied. A Bayesian model averaging approach exploits LD information considerably better than G-BLUP and thus is the most promising method. Cross-validations should account for family structure in the data to allow for long-lasting genomic based breeding plans in animal and plant breeding.

Project description:BackgroundThe focus in dairy cattle breeding is gradually shifting from production to functional traits and genetic parameters of calving traits are estimated more frequently. However, across countries, various statistical models are used to estimate these parameters. This study evaluates different models for calving ease and stillbirth in United Kingdom Holstein-Friesian cattle.MethodsData from first and later parity records were used. Genetic parameters for calving ease, stillbirth and gestation length were estimated using the restricted maximum likelihood method, considering different models i.e. sire (-maternal grandsire), animal, univariate and bivariate models. Gestation length was fitted as a correlated indicator trait and, for all three traits, genetic correlations between first and later parities were estimated. Potential bias in estimates was avoided by acknowledging a possible environmental direct-maternal covariance. The total heritable variance was estimated for each trait to discuss its theoretical importance and practical value. Prediction error variances and accuracies were calculated to compare the models.Results and discussionOn average, direct and maternal heritabilities for calving traits were low, except for direct gestation length. Calving ease in first parity had a significant and negative direct-maternal genetic correlation. Gestation length was maternally correlated to stillbirth in first parity and directly correlated to calving ease in later parities. Multi-trait models had a slightly greater predictive ability than univariate models, especially for the lowly heritable traits. The computation time needed for sire (-maternal grandsire) models was much smaller than for animal models with only small differences in accuracy. The sire (-maternal grandsire) model was robust when additional genetic components were estimated, while the equivalent animal model had difficulties reaching convergence.ConclusionsFor the evaluation of calving traits, multi-trait models show a slight advantage over univariate models. Extended sire models (-maternal grandsire) are more practical and robust than animal models. Estimated genetic parameters for calving traits of UK Holstein cattle are consistent with literature. Calculating an aggregate estimated breeding value including direct and maternal values should encourage breeders to consider both direct and maternal effects in selection decisions.

Project description:Little is known about the correlations between the genetic susceptibility/resistance to Mycobacterium avium subsp. paratuberculosis (Map) infection and the estimated breeding values for type, production and functional traits. Previously, we identified 70 combinations of five single nucleotide polymorphisms (SNPs) in four bovine innate immune genes (SLC11A1, SP110, TLR2, CD209) that are associated with the genetic risk of paratuberculosis (PTB) or Johne's disease progression, which can be graded as low (LOWIN), latent (LATIN), or patent (PATIN) risk. Other possible combinations of these 5 SNPs were grouped in the average group (AVERIN). In the current study, differences in estimated breeding values (EBVs) for several traits were analyzed using linear models in a large cohort of Holstein cows (N = 15656) genotyped across Spain in 2016 or 2017. After the assignment of each genotyped cow to a risk group, cows within the PATIN risk group (N = 1448) had a superior combined genetic index (2797.57), type genetic index (524.62), milk yield (653.92 kg), protein yield (21.77 kg), fat yield (24.82 kg) and economic merit index (125 Euros) compared with the other three risk groups. Statistically significant differences in the longevity scores between the cows that were included in the PATIN risk group (108.85) and the LOWIN (107.82) and AVERIN (107.92) groups were also observed. The associations between the genetic risk groups and PTB diagnostic results were validated in a population of 99 cows from a Spanish farm with a high prevalence of PTB. Significant differences in ELISA readings between the PATIN (65.49 %) and the AVERIN (15.97 %), LATIN (2.11 %), and LOWIN (3.27 %) groups were observed. In addition, significant differences in Map DNA copies/gram of feces were observed between the PATIN and the other three risk groups. These results together with the substantial economic impact of PTB in dairy cattle support the selection of the animals with less susceptibility to PTB in the Spanish breeding program.

Project description:Reducing the incidence of both the degree of assistance required at calving, as well as the extent of perinatal mortality (PM) has both economic and societal benefits. The existence of heritable genetic variability in both traits signifies the presence of underlying genomic variability. The objective of the present study was to locate regions of the genome, and by extension putative genes and mutations, that are likely to be underpinning the genetic variability in direct calving difficulty (DCD), maternal calving difficulty (MCD), and PM. Imputed whole-genome single-nucleotide polymorphism (SNP) data on up to 8,304 Angus (AA), 17,175 Charolais (CH), 16,794 Limousin (LM), and 18,474 Holstein-Friesian (HF) sires representing 5,866,712 calving events from descendants were used. Several putative quantitative trait loci (QTL) regions associated with calving performance both within and across dairy and beef breeds were identified, although the majority were both breed- and trait-specific. QTL surrounding and encompassing the myostatin (MSTN) gene were associated (P < 5 × 10-8) with DCD and PM in both the CH and LM populations. The well-known Q204X mutation was the fifth strongest association with DCD in the CH population and accounted for 5.09% of the genetic variance in DCD. In contrast, none of the 259 segregating variants in MSTN were associated (P > × 10-6) with DCD in the LM population but a genomic region 617 kb downstream of MSTN was associated (P < 5 × 10-8). The genetic architecture for DCD differed in the HF population relative to the CH and LM, where two QTL encompassing ZNF613 on Bos taurus autosome (BTA)18 and PLAG1 on BTA14 were identified in the former. Pleiotropic SNP associated with all three calving performance traits were also identified in the three beef breeds; 5 SNP were pleiotropic in AA, 116 in LM, and 882 in CH but no SNP was associated with more than one trait within the HF population. The majority of these pleiotropic SNP were on BTA2 surrounding MSTN and were associated with both DCD and PM. Multiple previously reported, but also novel QTL, associated with calving performance were detected in this large study. These also included QTL regions harboring SNP with the same direction of allele substitution effect for both DCD and MCD thus contributing to a more effective simultaneous selection for both traits.

Project description:In crossbreeding programs, various genomic prediction models have been proposed for using phenotypic records of crossbred animals to increase the selection response for crossbred performance in purebred animals. A possible model is a model that assumes identical single nucleotide polymorphism (SNP) effects for the crossbred performance trait across breeds (ASGM). Another model is a genomic model that assumes breed-specific effects of SNP alleles (BSAM) for crossbred performance. The aim of this study was to derive and validate equations for predicting the reliability of estimated genomic breeding values for crossbred performance in both these models. Prediction equations were derived for situations when all (phenotyping and) genotyping data have already been collected, i.e. based on the genetic evaluation model, and for situations when all genotyping data are not yet available, i.e. when designing breeding programs.When all genotyping data are available, prediction equations are based on selection index theory. Without availability of all genotyping data, prediction equations are based on population parameters (e.g., heritability of the traits involved, genetic correlation between purebred and crossbred performance, effective number of chromosome segments). Validation of the equations for predicting the reliability of genomic breeding values without all genotyping data was performed based on simulated data of a two-way crossbreeding program, using either two closely-related breeds, or two unrelated breeds, to produce crossbred animals. The proposed equations can be used for an easy comparison of the reliability of genomic estimated breeding values across many scenarios, especially if all genotyping data are available. We show that BSAM outperforms ASGM for a specific breed, if the effective number of chromosome segments that originate from this breed and are shared by selection candidates of this breed and crossbred reference animals is less than half the effective number of all chromosome segments that are independently segregating in the same animals.The derived equations can be used to predict the reliability of genomic estimated breeding values for crossbred performance using ASGM or BSAM in many scenarios, and are thus useful to optimize the design of breeding programs. Scenarios can vary in terms of the genetic correlation between purebred and crossbred performances, heritabilities, number of reference animals, or distance between breeds.

Project description:BACKGROUND:Calving difficulty or dystocia has a great economic impact in the US dairy industry. Reported risk factors associated with calving difficulty are feto-pelvic disproportion, gestation length and conformation. Different dairy cattle breeds have different incidence of calving difficulty, with Holstein having the highest dystocia rates and Jersey the lowest. Genomic selection becomes important especially for complex traits with low heritability, where the accuracy of conventional selection is lower. However, for complex traits where a large number of genes influence the phenotype, genome-wide association studies showed limitations. Biological networks could overcome some of these limitations and better capture the genetic architecture of complex traits. In this paper, we characterize Holstein, Brown Swiss and Jersey breed-specific dystocia networks and employ them in genomic predictions. RESULTS:Marker association analysis identified single nucleotide polymorphisms explaining the largest average proportion of genetic variance on BTA18 in Holstein, BTA25 in Brown Swiss, and BTA15 in Jersey. Gene networks derived from the genome-wide association included 1272 genes in Holstein, 1454 genes in Brown Swiss, and 1455 genes in Jersey. Furthermore, 256 genes in Holstein network, 275 genes in the Brown Swiss network, and 253 genes in the Jersey network were within previously reported dystocia quantitative trait loci. The across-breed network included 80 genes, with 9 genes being within previously reported dystocia quantitative trait loci. The gene-gene interactions in this network differed in the different breeds. Gene ontology enrichment analysis of genes in the networks showed Regulation of ARF GTPase was very significant (FDR???0.0098) on Holstein. Neuron morphogenesis and differentiation was the term most enriched (FDR???0.0539) on the across-breed network. Genomic prediction models enriched with network-derived relationship matrices did not outperform regular GBLUP models. CONCLUSIONS:Regions identified in the genome were in the proximity of previously described quantitative trait loci that would most likely affect calving difficulty by altering the feto-pelvic proportion. Inclusion of identified networks did not increase prediction accuracy. The approach used in this paper could be extended to any instance with asymmetric distribution of phenotypes, for example, resistance to disease data.

Dataset Information

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

Publications

Genome-wide association studies and genomic prediction of breeding values for calving performance and body conformation traits in Holstein cattle.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets