Genome-Wide Association Study for Feed Efficiency and Growth Traits in U.S
Total Page:16
File Type:pdf, Size:1020Kb
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Papers and Publications in Animal Science Animal Science Department 2017 Genome-wide association study for feed efficiency and owthgr traits in U.S. beef cattle Christopher M. Seabury Texas A&M University, [email protected] David L. Oldeschulte Texas A&M University, College Station, [email protected] Mahdi Saatchi Iowa State University, [email protected] Jonathan E. Beever University of Illinois at Urbana-Champaign, [email protected] Jared E. Decker University of Missouri, Columbia, [email protected] See next page for additional authors Follow this and additional works at: https://digitalcommons.unl.edu/animalscifacpub Part of the Genetics and Genomics Commons, and the Meat Science Commons Seabury, Christopher M.; Oldeschulte, David L.; Saatchi, Mahdi; Beever, Jonathan E.; Decker, Jared E.; Halley, Yvette A.; Bhattarai, Eric K.; Molael, Maral; Freetly, Harvey; Hansen, Stephanie L.; Yampara-Iquise, Helen; Johnson, Kristen A.; Kerley, Monty S.; Kim, JaeWoo; Loy, Daniel D.; Marques, Elisa; Neibergs, Holly L.; Schnabel, Robert D.; Shike, Daniel W.; Spangler, Matthew L.; Weaber, Bob; Garrick, D. J.; and Taylor, Jeremy F., "Genome-wide association study for feed efficiency and owthgr traits in U.S. beef cattle" (2017). Faculty Papers and Publications in Animal Science. 993. https://digitalcommons.unl.edu/animalscifacpub/993 This Article is brought to you for free and open access by the Animal Science Department at DigitalCommons@University of Nebraska - Lincoln. It has been accepted for inclusion in Faculty Papers and Publications in Animal Science by an authorized administrator of DigitalCommons@University of Nebraska - Lincoln. Authors Christopher M. Seabury, David L. Oldeschulte, Mahdi Saatchi, Jonathan E. Beever, Jared E. Decker, Yvette A. Halley, Eric K. Bhattarai, Maral Molael, Harvey Freetly, Stephanie L. Hansen, Helen Yampara-Iquise, Kristen A. Johnson, Monty S. Kerley, JaeWoo Kim, Daniel D. Loy, Elisa Marques, Holly L. Neibergs, Robert D. Schnabel, Daniel W. Shike, Matthew L. Spangler, Bob Weaber, D. J. Garrick, and Jeremy F. Taylor This article is available at DigitalCommons@University of Nebraska - Lincoln: https://digitalcommons.unl.edu/ animalscifacpub/993 Seabury et al. BMC Genomics (2017) 18:386 DOI 10.1186/s12864-017-3754-y RESEARCH ARTICLE Open Access Genome-wide association study for feed efficiency and growth traits in U.S. beef cattle Christopher M. Seabury1*, David L. Oldeschulte1, Mahdi Saatchi2, Jonathan E. Beever3, Jared E. Decker4,5, Yvette A. Halley1, Eric K. Bhattarai1, Maral Molaei1, Harvey C. Freetly6, Stephanie L. Hansen2, Helen Yampara-Iquise4, Kristen A. Johnson7, Monty S. Kerley4, JaeWoo Kim4, Daniel D. Loy2, Elisa Marques8, Holly L. Neibergs7, Robert D. Schnabel4,5, Daniel W. Shike3, Matthew L. Spangler9, Robert L. Weaber10, Dorian J. Garrick2,11 and Jeremy F. Taylor4* Abstract Background: Single nucleotide polymorphism (SNP) arrays for domestic cattle have catalyzed the identification of genetic markers associated with complex traits for inclusion in modern breeding and selection programs. Using actual and imputed Illumina 778K genotypes for 3887 U.S. beef cattle from 3 populations (Angus, Hereford, SimAngus), we performed genome-wide association analyses for feed efficiency and growth traits including average daily gain (ADG), dry matter intake (DMI), mid-test metabolic weight (MMWT), and residual feed intake (RFI), with marker-based heritability estimates produced for all traits and populations. Results: Moderate and/or large-effect QTL were detected for all traits in all populations, as jointly defined by the estimated proportion of variance explained (PVE) by marker effects (PVE ≥ 1.0%) and a nominal P-value threshold (P ≤ 5e-05). Lead SNPs with PVE ≥ 2.0% were considered putative evidence of large-effect QTL (n = 52), whereas those with PVE ≥ 1.0% but < 2.0% were considered putative evidence for moderate-effect QTL (n = 35). Identical or proximal lead SNPs associated with ADG, DMI, MMWT, and RFI collectively supported the potential for either pleiotropic QTL, or independent but proximal causal mutations for multiple traits within and between the analyzed populations. Marker-based heritability estimates for all investigated traits ranged from 0.18 to 0.60 using 778K genotypes, or from 0.17 to 0.57 using 50K genotypes (reduced from Illumina 778K HD to Illumina Bovine SNP50). An investigation to determine if QTL detected by 778K analysis could also be detected using 50K genotypes produced variable results, suggesting that 50K analyses were generally insufficient for QTL detection in these populations, and that relevant breeding or selection programs should be based on higher density analyses (imputed or directly ascertained). Conclusions: Fourteen moderate to large-effect QTL regions which ranged from being physically proximal (lead SNPs ≤ 3Mb) to fully overlapping for RFI, DMI, ADG, and MMWT were detected within and between populations, and included evidence for pleiotropy, proximal but independent causal mutations, and multi-breed QTL. Bovine positional candidate genes for these traits were functionally conserved across vertebrate species. Keywords: GWAS, QTL, Feed efficiency and growth, Beef Cattle * Correspondence: [email protected]; [email protected] 1Department of Veterinary Pathobiology, Texas A&M University, College Station 77843, USA 4Division of Animal Sciences, University of Missouri, Columbia 65211, USA Full list of author information is available at the end of the article © The Author(s). 2017 Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Seabury et al. BMC Genomics (2017) 18:386 Page 2 of 25 Background (DMI; lb/d), average daily gain on feed (ADG; lb/d), and Feed efficiency analyses for beef cattle and other import- mid-test metabolic body weight (MMWT; lb0.75)inU.S. ant food animals generally seek to relate measures of Angus, Hereford, and SimAngus (Simmental × Angus) feed intake with measures of animal productivity (i.e., beef cattle [13], with corresponding heritability estimates agricultural input versus output). Relevant to beef cattle produced for all investigated traits. Thereafter, we investi- production, and particularly among feedlot operations, gated whether QTL discovered using the Illumina the most costly component of the production process is BovineHD markers (hereafter 778K) were also detectable feed, with feed-based expenditures during certain stages using the Illumina BovineSNP50 Assay (hereafter 50K), of production representing more than 80% of the total which is important for modern genomic selection pro- costs [1]. The overall economic importance of feed effi- grams predicated on historic 50K analyses, and also ciency in beef cattle was initially recognized more than compared the marker-based heritability estimates pro- 40 years ago [2, 3], with modern improvements in the duced by the two arrays. Finally, we evaluated our 778K efficiency of feed utilization expected to yield positive results within the context of the established GWAS lit- economic returns across many facets of the beef industry erature, and found some positional candidate genes [1, 4]. Moreover, an estimated annual cost savings ex- underlying bovine feed efficiency and/or growth related ceeding $1 billion U.S. dollars could likely be achieved QTL to be conserved across divergent food-animal spe- by increasing the efficiency of feed utilization by 10% in cies, with a tangible proportion of our bovine GWAS U.S. beef cattle alone [5]. results also overlapping with specific loci implicated in Traditional measures of feed efficiency, such as the studies of obesity or metabolic syndrome, diabetes, and ratio of feed consumed to observed body weight gain insulin resistance among humans and/or mice. The re- (i.e., feed conversion ratio), are likely to be very success- sults of this study are immediately useful for enabling ful in selecting for increased growth rates in beef cattle genomic selection in beef cattle, but also suggest that (i.e., increased producer income), but selection for en- domestic cattle may be relevant models for some hu- hanced weaning or yearling weights may also lead to in- man biomedical studies. creased costs associated with maintaining larger mature cows (i.e., increased nutrient requirements and feed costs; Results and discussion increased calving difficulty due to larger birth weights, Heritability estimates for RFI, DMI, ADG, and MMWT etc) [4, 6, 7]. An alternative measure that has gained in Using a genomic relationship matrix [23] that was nor- popularity among livestock species is Residual Feed Intake malized using Gower’s centering approach [24, 25], to (RFI), with this trait often defined as the difference be- yield a sampling variance of 1.0, we estimated the pseudo- 2 2 2 2 tween an animal’s observed and expected feed intake in heritability [24, 25] [(i. e., ha = σa/(σa + σe); also repre- 2 relation to the animal’s