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

The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds

peer-reviewedDomestication and the subsequent selection of animals for either economic or morphological features can leave a variety of imprints on the genome of a population. Genomic regions subjected to high selective pressures often show reduced genetic diversity and frequent runs of homozygosity (ROH). Therefore, the objective of the present study was to use 42,182 autosomal SNPs to identify genomic regions in 3,191 sheep from six commercial breeds subjected to selection pressure and to quantify the genetic diversity within each breed using ROH. In addition, the historical effective population size of each breed was also estimated and, in conjunction with ROH, was used to elucidate the demographic history of the six breeds. ROH were common in the autosomes of animals in the present study, but the observed breed differences in patterns of ROH length and burden suggested differences in breed effective population size and recent management. ROH provided a sufficient predictor of the pedigree inbreeding coefficient, with an estimated correlation between both measures of 0.62. Genomic regions under putative selection were identified using two complementary algorithms; the fixation index and hapFLK. The identified regions under putative selection included candidate genes associated with skin pigmentation, body size and muscle formation; such characteristics are often sought after in modern-day breeding programs. These regions of selection frequently overlapped with high ROH regions both within and across breeds. Multiple yet uncharacterised genes also resided within putative regions of selection. This further substantiates the need for a more comprehensive annotation of the sheep genome as these uncharacterised genes may contribute to traits of interest in the animal sciences. Despite this, the regions identified as under putative selection in the current study provide an insight into the mechanisms leading to breed differentiation and genetic variation in meat production.This work was supported by MultiGS Research Stimulus Fund (11/S/112) and OviGen project (14/S/849) which are funded by the Department of Agriculture, Food and Marine, Ireland

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

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 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

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

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

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

Genome-wide association analyses of carcass traits using copy number variants and raw intensity values of single nucleotide polymorphisms in cattle

peer-reviewedBackground The carcass value of cattle is a function of carcass weight and quality. Given the economic importance of carcass merit to producers, it is routinely included in beef breeding objectives. A detailed understanding of the genetic variants that contribute to carcass merit is useful to maximize the efficiency of breeding for improved carcass merit. The objectives of the present study were two-fold: firstly, to perform genome-wide association analyses of carcass weight, carcass conformation, and carcass fat using copy number variant (CNV) data in a population of 923 Holstein-Friesian, 945 Charolais, and 974 Limousin bulls; and secondly to perform separate association analyses of carcass traits on the same population of cattle using the Log R ratio (LRR) values of 712,555 single nucleotide polymorphisms (SNPs). The LRR value of a SNP is a measure of the signal intensity of the SNP generated during the genotyping process. Results A total of 13,969, 3,954, and 2,805 detected CNVs were tested for association with the three carcass traits for the Holstein-Friesian, Charolais, and Limousin, respectively. The copy number of 16 CNVs and the LRR of 34 SNPs were associated with at least one of the three carcass traits in at least one of the three cattle breeds. With the exception of three SNPs, none of the quantitative trait loci detected in the CNV association analyses or the SNP LRR association analyses were also detected using traditional association analyses based on SNP allele counts. Many of the CNVs and SNPs associated with the carcass traits were located near genes related to the structure and function of the spliceosome and the ribosome; in particular, U6 which encodes a spliceosomal subunit and 5S rRNA which encodes a ribosomal subunit. Conclusions The present study demonstrates that CNV data and SNP LRR data can be used to detect genomic regions associated with carcass traits in cattle providing information on quantitative trait loci over and above those detected using just SNP allele counts, as is the approach typically employed in genome-wide association analyses

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