Genome-wide association study for calving performance using high-density genotypes in dairy and beef cattle

peer-reviewedBackground Calving difficulty and perinatal mortality are prevalent in modern-day cattle production systems. It is well-established that there is a genetic component to both traits, yet little is known about their underlying genomic architecture, particularly in beef breeds. Therefore, we performed a genome-wide association study using high-density genotypes to elucidate the genomic architecture of these traits and to identify regions of the bovine genome associated with them. Results Genomic regions associated with calving difficulty (direct and maternal) and perinatal mortality were detected using two statistical approaches: (1) single-SNP (single nucleotide polymorphism) regression and (2) a Bayesian approach. Data included high-density genotypes on 770 Holstein-Friesian, 927 Charolais and 963 Limousin bulls. Several novel or previously identified genomic regions were detected but associations differed by breed. For example, two genomic associations, one each on chromosomes 18 and 2 explained 2.49 % and 3.13 % of the genetic variance in direct calving difficulty in the Holstein-Friesian and Charolais populations, respectively. Imputed Holstein-Friesian sequence data was used to refine the genomic regions responsible for significant associations. Several candidate genes on chromosome 18 were identified and four highly significant missense variants were detected within three of these genes (SIGLEC12, CTU1, and ZNF615). Nevertheless, only CTU1 contained a missense variant with a putative impact on direct calving difficulty based on SIFT (0.06) and Polyphen (0.95) scores. Using imputed sequence data, we refined a genomic region on chromosome 4 associated with maternal calving difficulty in the Holstein-Friesian population and found the strongest association with an intronic variant in the PCLO gene. A meta-analysis was performed across the three breeds for each calving performance trait to identify common variants associated with these traits in the three breeds. Our results suggest that a portion of the genetic variation in calving performance is common to all three breeds. Conclusion The genomic architecture of calving performance is complex and mainly influenced by many polymorphisms of small effect. We identified several associations of moderate effect size but the majority were breed-specific, indicating that breed-specific alleles exist for calving performance or that the linkage phase between genotyped allele and causal mutation varies between breeds

Genomic Regions Associated With Gestation Length Detected Using Whole-Genome Sequence Data Differ Between Dairy and Beef Cattle

peer-reviewedWhile many association studies exist that have attempted to relate genomic markers to phenotypic performance in cattle, very few have considered gestation length as a phenotype, and of those that did, none used whole genome sequence data from multiple breeds. The objective of the present study was therefore to relate imputed whole genome sequence data to estimated breeding values for gestation length using 22,566 sires (representing 2,262,706 progeny) of multiple breeds [Angus (AA), Charolais (CH), Holstein-Friesian (HF), and Limousin (LM)]. The associations were undertaken within breed using linear mixed models that accounted for genomic relatedness among sires; a separate association analysis was undertaken with all breeds analysed together but with breed included as a fixed effect in the model. Furthermore, the genome was divided into 500 kb segments and whether or not segments harboured a single nucleotide polymorphism (SNP) with a P ≤ 1 × 10-4 common to different combinations of breeds was determined. Putative quantitative trait loci (QTL) regions associated with gestation length were detected in all breeds; significant associations with gestation length were only detected in the HF population and in the across-breed analysis of all 22,566 sires. Twenty-five SNPs were significantly associated (P ≤ 5 × 10-8) with gestation length in the HF population. Of the 25 significant SNPs, 18 were located within three QTLs on Bos taurus autosome number (BTA) 18, six were in two QTL on BTA19, and one was located within a QTL on BTA7. The strongest association was rs381577268, a downstream variant of ZNF613 located within a QTL spanning from 58.06 to 58.19 Mb on BTA18; it accounted for 1.37% of the genetic variance in gestation length. Overall there were 11 HF animals within the edited dataset that were homozygous for the T allele at rs381577268 and these had a 3.3 day longer (P < 0.0001) estimated breeding value (EBV) for gestation length than the heterozygous animals and a 4.7 day longer (P < 0.0001) EBV for gestation length than the homozygous CC animals. The majority of the 500 kb windows harboring a SNP with a P ≤ 1 × 10-4 were unique to a single breed and no window was shared among all four breeds for gestation length, suggesting any QTLs identified are breed-specific associations

Reaffirmation of known major genes and the identification of novel candidate genes associated with carcass-related metrics based on whole genome sequence within a large multi-breed cattle population

peer-reviewedBackground: The high narrow sense heritability of carcass traits suggests that the underlying additive genetic potential of an individual should be strongly correlated with both animal carcass quality and quantity, and therefore, by extension, carcass value. Therefore, the objective of the present study was to detect genomic regions associated with three carcass traits, namely carcass weight, conformation and fat cover, using imputed whole genome sequence in 28,470 dairy and beef sires from six breeds with a total of 2,199,926 phenotyped progeny. Results: Major genes previously associated with carcass performance were identified, as well as several putative novel candidate genes that likely operate both within and across breeds. The role of MSTN in carcass performance was re-affirmed with the segregating Q204X mutation explaining 1.21, 1.11 and 5.95% of the genetic variance in carcass weight, fat and conformation, respectively in the Charolais population. In addition, a genomic region on BTA6 encompassing the NCAPG/LCORL locus, which is a known candidate locus associated with body size, was associated with carcass weight in Angus, Charolais and Limousin. Novel candidate genes identified included ZFAT in Angus, and SLC40A1 and the olfactory gene cluster on BTA15 in Charolais. Although the majority of associations were breed specific, associations that operated across breeds included SORCS1 on BTA26, MCTP2 on BTA21 and ARL15 on BTA20; these are of particular interest due to their potential informativeness in across-breed genomic evaluations. Genomic regions affecting all three carcass traits were identified in each of the breeds, although these were mainly concentrated on BTA2 and BTA6, surrounding MSTN and NCAPG/LCORL, respectively. This suggests that although major genes may be associated with all three carcass traits, the majority of genes containing significant variants (unadjusted p-value < 10− 4 ) may be trait specific associations of small effect. Conclusions: Although plausible novel candidate genes were identified, the proportion of variance explained by these candidates was minimal thus reaffirming that while carcass performance may be affected by major genes in the form of MSTN and NCAPG/LCORL, the majority of variance is attributed to the additive (and possibly multiplicative) effect of many polymorphisms of small effect

The relationship between serum insulin-like growth factor-1 (IGF-1) concentration and reproductive performance, and genome-wide associations for serum IGF-1 in Holstein cows

peer-reviewedThe objectives of this study were to determine (1) factors associated with serum concentration of insulin-like growth factor-1 (IGF-1); (2) the relationship between serum IGF-1 concentration during the first week postpartum and ovarian cyclicity status by 35 d postpartum (DPP); (3) an optimum serum IGF-1 concentration threshold predictive of pregnancy to first artificial insemination (P/AI), including its diagnostic values; (4) the associations among categories of serum IGF-1 concentration and reproductive outcomes (P/AI and pregnancy risk up to 150 and 250 DPP); and (5) single nucleotide polymorphisms (SNP) associated with phenotypic variation in serum IGF-1 concentration in dairy cows. Serum IGF-1 concentration was determined at 7 (±2.4; ±standard error of the mean) DPP in 647 lactating Holstein cows (213 primiparous, 434 multiparous) from 7 herds in Alberta, Canada. The overall mean, median, minimum, and maximum serum IGF-1 concentrations during the first week postpartum were 37.8 (±1.23), 31.0, 20.0, and 225.0 ng/mL, respectively. Herd, age, parity, precalving body condition score, and season of blood sampling were all identified as factors associated with serum IGF-1 concentrations. Although serum IGF-1 concentration during the first week postpartum had no association with ovarian cyclicity status by 35 DPP in primiparous cows, it was greater in cyclic than in acyclic multiparous cows (32.2 vs. 27.4 ng/mL, respectively). The optimum serum IGF-1 thresholds predictive of P/AI were 85.0 ng/mL (sensitivity = 31.9%; specificity = 89.1%) and 31.0 ng/mL (sensitivity = 45.5%; specificity = 66.9%) for primiparous and multiparous cows, respectively. When cows were grouped into either high or low IGF-1 categories (greater or less than or equal to 85.0 ng/mL for primiparous cows and greater or less than or equal to 31.0 ng/mL for multiparous cows, respectively), primiparous cows with high IGF-1 had 4.43 times greater odds of P/AI and a tendency for higher pregnancy risk up to 150 DPP than those with low IGF-1, but not up to 250 DPP. Likewise, multiparous cows with high IGF-1 had 1.61 times greater odds of P/AI than those with low IGF-1. Pregnancy risk up to 150 and 250 DPP, however, did not differ between IGF-1 categories in multiparous cows. Moreover, 37 SNP across 10 Bos taurus autosomes were associated with variation in serum IGF-1 concentration, and 4 previously identified candidate genes related to fertility that were in linkage disequilibrium with some of these SNP were also identified

Genome-wide association study of endo-parasite phenotypes using imputed whole-genome sequence data in dairy and beef cattle

peer-reviewedBackground: Quantitative genetic studies suggest the existence of variation at the genome level that affects the ability of cattle to resist to parasitic diseases. The objective of the current study was to identify regions of the bovine genome that are associated with resistance to endo-parasites. Methods: Individual cattle records were available for Fasciola hepatica-damaged liver from 18 abattoirs. Deregressed estimated breeding values (EBV) for F. hepatica-damaged liver were generated for genotyped animals with a record for F. hepatica-damaged liver and for genotyped sires with a least one progeny record for F. hepatica-damaged liver; 3702 animals were available. In addition, individual cow records for antibody response to F. hepatica on 6388 genotyped dairy cows, antibody response to Ostertagia ostertagi on 8334 genotyped dairy cows and antibody response to Neospora caninum on 4597 genotyped dairy cows were adjusted for non-genetic effects. Genotypes were imputed to whole-sequence; after edits, 14,190,141 single nucleotide polymorphisms (SNPs) and 16,603,644 SNPs were available for cattle with deregressed EBV for F. hepatica-damaged liver and cows with an antibody response to a parasitic disease, respectively. Association analyses were undertaken using linear regression on one SNP at a time, in which a genomic relationship matrix accounted for the relationships between animals. Results: Genomic regions for F. hepatica-damaged liver were located on Bos taurus autosomes (BTA) 1, 8, 11, 16, 17 and 18; each region included at least one SNP with a p value lower than 10−6. Five SNPs were identified as significant (q value < 0.05) for antibody response to N. caninum and were located on BTA21 or 25. For antibody response to F. hepatica and O. ostertagi, six and nine quantitative trait loci (QTL) regions that included at least one SNP with a p value lower than 10−6 were identified, respectively. Gene set enrichment analysis revealed a significant association between functional annotations related to the olfactory system and QTL that were suggestively associated with endo-parasite phenotypes. Conclusions: A number of novel genomic regions were suggestively associated with endo-parasite phenotypes across the bovine genome and two genomic regions on BTA21 and 25 were associated with antibody response to N. caninum

Breed- and trait-specific associations define the genetic architecture of calving performance traits in cattle

peer-reviewedReducing 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.Science Foundation Irelan

The relationship between serum anti-Müllerian hormone concentrations and fertility, and genome-wide associations for anti-Müllerian hormone in Holstein cows

peer-reviewedThe objectives of this study were to (1) evaluate factors associated with variation in circulating anti-Müllerian hormone (AMH) concentrations, (2) establish an optimum AMH threshold predictive of pregnancy to first artificial insemination (P/AI), (3) examine the relationship between AMH and fertility (P/AI, pregnancy loss between 30 and 60 d after artificial insemination, and pregnancy risk up to 250 d postpartum), and (4) identify quantitative trait loci associated with phenotypic variation of AMH concentrations in dairy cows. Serum AMH concentrations (pg/mL) were determined at 7 ± 2.4 d postpartum in 647 lactating Holstein cows (213 primiparous, 434 multiparous) from 1 research and 6 commercial dairy herds in Alberta, Canada. Of these, 589 cows were genotyped on the 26K Bovine BeadChip (Neogen Inc., Lincoln, NE) and subsequently imputed to the Illumina Bovine High Density BeadChip (Illumina, San Diego, CA) for genome-wide association analysis for variation in serum AMH concentrations. Factors associated with variation in serum AMH concentrations and the relationship between categories of AMH and aforementioned fertility outcomes were evaluated only in a subset of 460 cows that had a complete data set available. The overall mean (±standard error of the mean), median, minimum, and maximum AMH concentrations were 191.1 ± 6.3, 151.7, 13.9, and 1,879.0 pg/mL, respectively. The AMH concentrations were not associated with herd, precalving body condition score, postpartum week, and season of sampling; the lactation number, however, had a quadratic relationship with serum AMH concentrations (116.2, 204.9 204.5, and 157.9 pg/mL for first, second, third, and ≥fourth lactation, respectively). The optimum AMH threshold predictive of P/AI could not be established because the receiver operating characteristic curve analysis model was nonsignificant. Categories of AMH [low (285.0 pg/mL; n = 92) based on lowest 20%, intermediate 60%, and highest 20% serum AMH) had no associations with P/AI (34, 43, and 40%), pregnancy loss between 30 and 60 d after artificial insemination (20, 12, and 8%), or pregnancy risk up to 250 d postpartum. One candidate gene associated with AMH production [AMH gene on Bos taurus autosome (BTA) 7] and 4 candidate genes related to embryo development (SCAI and PPP6C genes on BTA11 and FGF18 and EEF2K genes on BTA20 and BTA25, respectively) were in linkage disequilibrium with single nucleotide polymorphisms associated with phenotypic variation in serum AMH in dairy cows

Genomic regions associated with muscularity in beef cattle differ in five contrasting cattle breeds

peer-reviewedBackground Linear type traits, which reflect the muscular characteristics of an animal, could provide insight into how, in some cases, morphologically very different animals can yield the same carcass weight. Such variability may contribute to differences in the overall value of the carcass since primal cuts vary greatly in price; such variability may also hinder successful genome-based association studies. Therefore, the objective of our study was to identify genomic regions that are associated with five muscularity linear type traits and to determine if these significant regions are common across five different breeds. Analyses were carried out using linear mixed models on imputed whole-genome sequence data in each of the five breeds, separately. Then, the results of the within-breed analyses were used to conduct an across-breed meta-analysis per trait. Results We identified many quantitative trait loci (QTL) that are located across the whole genome and associated with each trait in each breed. The only commonality among the breeds and traits was a large-effect pleiotropic QTL on BTA2 that contained the MSTN gene, which was associated with all traits in the Charolais and Limousin breeds. Other plausible candidate genes were identified for muscularity traits including PDE1A, PPP1R1C and multiple collagen and HOXD genes. In addition, associated (gene ontology) GO terms and KEGG pathways tended to differ between breeds and between traits especially in the numerically smaller populations of Angus, Hereford, and Simmental breeds. Most of the SNPs that were associated with any of the traits were intergenic or intronic SNPs located within regulatory regions of the genome. Conclusions The commonality between the Charolais and Limousin breeds indicates that the genetic architecture of the muscularity traits may be similar in these breeds due to their similar origins. Conversely, there were vast differences in the QTL associated with muscularity in Angus, Hereford, and Simmental. Knowledge of these differences in genetic architecture between breeds is useful to develop accurate genomic prediction equations that can operate effectively across breeds. Overall, the associated QTL differed according to trait, which suggests that breeding for a morphologically different (e.g. longer and wider versus shorter and smaller) more efficient animal may become possible in the future

Genomic Regions Associated With Skeletal Type Traits in Beef and Dairy Cattle Are Common to Regions Associated With Carcass Traits, Feed Intake and Calving Difficulty

Linear type traits describing the skeletal characteristics of an animal are moderately to strongly genetically correlated with a range of other performance traits in cattle including feed intake, reproduction traits and carcass merit; thus, type traits could also provide useful insights into the morphological differences among animals underpinning phenotypic differences in these complex traits. The objective of the present study was to identify genomic regions associated with five subjectively scored skeletal linear traits, to determine if these associated regions are common in multiple beef and dairy breeds, and also to determine if these regions overlap with those proposed elsewhere to be associated with correlated performance traits. Analyses were carried out using linear mixed models on imputed whole genome sequence data separately in 1,444 Angus, 1,129 Hereford, 6,433 Charolais, 8,745 Limousin, 1,698 Simmental, and 4,494 Holstein-Friesian cattle, all scored for the linear type traits. There was, on average, 18 months difference in age at assessment of the beef versus the dairy animals. While the majority of the identified quantitative trait loci (QTL), and thus genes, were both trait-specific and breed-specific, a large-effect pleiotropic QTL on BTA6 containing the NCAPG and LCORL genes was associated with all skeletal traits in the Limousin population and with wither height in the Angus. Other than that, little overlap existed in detected QTLs for the skeletal type traits in the other breeds. Only two QTLs overlapped the beef and dairy breeds; both QTLs were located on BTA5 and were associated with height in both the Angus and the Holstein-Friesian, despite the difference in age at assessment. Several detected QTLs in the present study overlapped with QTLs documented elsewhere that are associated with carcass traits, feed intake, and calving difficulty. While most breeding programs select for the macro-traits like carcass weight, carcass conformation, and feed intake, the higher degree of granularity with selection on the individual linear type traits in a multi-trait index underpinning the macro-level goal traits, presents an opportunity to help resolve genetic antagonisms among morphological traits in the pursuit of the animal with optimum performance metrics.</p

The relationship between anogenital distance and fertility, and genome-wide associations for anogenital distance in Irish Holstein-Friesian cows

peer-reviewedThe evaluation of anogenital distance (AGD), the distance from the center of the anus to base of the clitoris, as a potential fertility trait for genetic selection in dairy cows has generated recent interest. The objectives of this cross-sectional observational study were to (1) characterize the distribution and variability of AGD, (2) determine factors associated with AGD, (3) estimate heritability for AGD, (4) identify single nucleotide polymorphisms (SNP) associated with phenotypic variation of AGD, and (5) validate the relationship between categories of AGD and fertility in Irish Holstein-Friesian cows. Anogenital distance was measured using digital calipers in 1,180 Holstein cows (mean ± standard deviation: 225 ± 79 d in milk) from 10 dairy herds located in Munster, Ireland. In addition, age (yr), weight (kg), height at hip (cm), and body condition score (BCS) at the time of AGD measurement were determined in a subset of 281 cows. Genotype information available from 908 cows was subsequently imputed to the Illumina Bovine High Density BeadChip (Illumina Inc., San Diego, CA) for genome-wide association analysis of phenotypic variation in AGD. Overall, AGD had a normal distribution and high variability (mean ± standard deviation; 119.2 ± 11.6 mm). Anogenital distance was weakly but positively associated with cow age, hip height, and body weight, and negatively associated with BCS; the phenotypic variation in AGD that was explainable by these variables was small (coefficient of determination; R2 = 0.09, 0.06, 0.10, and 0.02, respectively). The estimated heritability for AGD was 0.37 (standard error of mean ± 0.08). Six SNP of suggestive significance were identified on Bos taurus autosomes 6, 15, 20, and 26; however, none of these SNP was related to previously identified candidate genes for fertility. Cows were categorized into quartiles (Q1; 86 to 111 mm; n = 311, Q2; 112 to 120 mm; n = 330; Q3; 121 to 127 mm; n = 265, and Q4; 128 to 160 mm; n = 274) based on AGD and the association with reproductive outcomes examined (21-d submission rate, pregnancy to first AI, pregnancy rate within 21, 42 and 84-d after the farm mating start date, and number of times bred). None of the reproductive variables differed significantly between AGD categories. In summary, despite identification of high variability and moderate heritability for AGD in Irish Holstein-Friesian cows, reproductive outcomes did not differ between categories of AGD. This latter result differs from our previous finding of an inverse relationship between AGD and pregnancy outcomes in first- and second-parity Canadian Holstein cows, emphasizing the need to test and validate this new phenotype in diverse cow populations