Genome-wide association studies for feed intake and efficiency in two laying periods of chickens

Abstract

Background

Feed contributes to over 60 % of the total production costs in the poultry industry. Increasing feed costs prompt geneticists to include feed intake and efficiency as selection goals in breeding programs. In the present study, we used an F2 chicken population in a genome-wide association study (GWAS) to detect potential genetic variants and candidate genes associated with daily feed intake (FI) and feed efficiency, including residual feed intake (RFI) and feed conversion ratio (FCR).

Methods

A total of 1534 F2 hens from a White Leghorn and Dongxiang reciprocal cross were phenotyped for feed intake and efficiency between 37 and 40 weeks (FI1, RFI1, and FCR1) and between 57 and 60 weeks (FI2, RFI2, and FCR2), and genotyped using the chicken 600 K single nucleotide polymorphism (SNP) genotyping array. Univariate, bivariate, and conditional genome-wide association studies (GWAS) were performed with GEMMA, a genome-wide efficient mixed model association algorithm. The statistical significance threshold for association was inferred by the simpleM method.

Results

We identified eight genomic regions that each contained at least one genetic variant that showed a significant association with FI. Genomic regions on Gallus gallus (GGA) chromosome 4 coincided with known quantitative trait loci (QTL) that affect feed intake of layers. Of particular interest, eight SNPs on GGA1 in the region between 169.23 and 171.55 Mb were consistently associated with FI in both univariate and bivariate GWAS, which explained 3.72 and 2.57 % of the phenotypic variance of FI1 and FI2, respectively. The CAB39L gene can be considered as a promising candidate for FI1. For RFI, a haplotype block on GGA27 harbored a significant SNP associated with RFI2. The major allele of rs315135692 was favorable for a lower RFI, with a phenotypic difference of 3.35 g/day between opposite homozygous genotypes. Strong signals on GGA1 were detected in the bivariate GWAS for FCR.

Conclusions

The results demonstrated the polygenic nature of feed intake. GWAS identified novel variants and confirmed a QTL that was previously reported for feed intake in chickens. Genetic variants associated with feed efficiency may be used in genomic breeding programs to select more efficient layers.

Research Article

Open Access

References

Hume DA, Whitelaw CBA, Archibald AL. The future of animal production: improving productivity and sustainability. J Agr Sci. 2011;149:9-16.

Gunsett FC. Problems associated with selection for traits defined as a ratio of two component traits. In: Proceedings of the 3rd World Congress of Genetics Applied to Livestock Production: 16-22 July 1986, Lincoln. 1986;11:437-42.

Flock DK. Genetic-economic aspects of feed efficiency in laying hens. World Poultry Sci J. 1998;54:225-39.

Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157:1819-29.

Tabor HK, Risch NJ, Myers RM. Candidate-gene approaches for studying complex genetic traits: practical considerations. Nat Rev Genet. 2002;3:391-7.

Fan B, Du Z, Gorbach DM, Rothschild MF. Development and application of high-density SNP arrays in genomic studies of domestic animals. Asian Australas J Anim Sci. 2010;23:833-47.

Mackay TFC, Stone EA, Ayroles JF. The genetics of quantitative traits: challenges and prospects. Nat Rev Genet. 2009;10:565-77.

Zhang H, Wang Z, Wang S, Li H. Progress of genome wide association study in domestic animals. J Anim Sci Biotechnol. 2012;3:26.

Yuan J, Dou T, Ma M, Yi G, Chen S, Qu LJ, et al. Genetic parameters of feed efficiency traits in laying period of chickens. Poult Sci. 2015;94:1470-5.

Luiting P, Urff EM. Optimization of a model to estimate residual feed consumption in the laying hen. Livest Prod Sci. 1991;27:321-38.

Kranis A, Gheyas A, Boschiero C, Turner F, Yu L, Smith S, et al. Development of a high density 600 K SNP genotyping array for chicken. BMC Genomics. 2013;14:59.

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira Manuel AR, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559-75.

Marchini M, Howie B. Genotype imputation for genome-wide association studies. Nat Rev Genet. 2010;11:499-511.

Browning SR, Browning BL. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am J Hum Genet. 2007;81:1084-97.

Gu X, Feng C, Ma L, Song C, Wang Y, Yang D, et al. Genome-wide association study of body weight in chicken F2 resource population. PLoS One. 2011;6:e21872.

Price AL, Zaitlen NA, Reich D, Patterson N. New approaches to population stratification in genome-wide association studies. Nat Rev Genet. 2010;11:459-63.

Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012;44:821-4.

Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263-5.

Zhou X, Stephens M. Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat Methods. 2014;11:407-9.

Yang J, Ferreira T, Morris AP, Medland SE, Madden PAF, Heath AC, et al. Conditional and joint multiple-SNP analysis of GWAS summary statistics identifies additional variants influencing complex traits. Nat Genet. 2012;44:369-75.

Johnson RC, Nelson GW, Troyer JL, Lautenberger JA, Kessing BD, Winkler CA, et al. Accounting for multiple comparisons in a genome-wide association study (GWAS). BMC Genomics. 2010;11:724.

Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol. 2008;32:361-9.

Team RC. R: A language and environment for statistical computing. Vienna: R foundation for statistical computing; 2012. ISBN 3-900051-07-0.

Yang J, Benyamin B, McEvoy BP, Gordon S, Henders AK, Nyholt DR, et al. Estimation of pleiotropy between complex diseases using SNP-derived genomic relationships and restricted maximum likelihood. Bioinformatics. 2012;28:2540-2.

Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88:76-82.

Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31:3812-4.

Sherman EL, Nkrumah JD, Moore SS. Whole genome single nucleotide polymorphism associations with feed intake and feed efficiency in beef cattle. J Anim Sci. 2010;88:16-22.

Saatchi M, Garrick DJ, Taylor J. Genome-wide association study of feed efficiency in beef cattle. Animal Industry Report. 2014. http://lib.dr.iastate.edu/ans_air/vol660/iss1/20. Accessed 5 Jan 2015.

Do DN, Ostersen T, Strathe AB, Mark T, Jensen J, Kadarmideen HN. Genome-wide association and systems genetic analyses of residual feed intake, daily feed consumption, backfat and weight gain in pigs. BMC Genet. 2014;15:27.

Sahana G, Kadlecová V, Hornshøj H, Nielsen B, Christensen OF. A genome-wide association scan in pig identifies novel regions associated with feed efficiency trait. J Anim Sci. 2013;91:1041-50.

Onteru SK, Gorbach DM, Young JM, Garrick DJ, Dekkers JCM, Rothschild MF. Whole genome association studies of residual feed intake and related traits in the pig. PLoS One. 2013;8:e61756.

Bolormaa S, Pryce JE, Hayes BJ, Goddard ME. Multivariate analysis of a genome-wide association study in dairy cattle. J Dairy Sci. 2010;93:3818-33.

Kim S, Xing EP. Statistical estimation of correlated genome associations to a quantitative trait network. PLoS Genet. 2009;5:e1000587.

Zhu W, Zhang H. Why do we test multiple traits in genetic association studies? J Korean Stat Soc. 2009;38:1-10.

Yang J, Weedon MN, Purcell S, Lettre G, Estrada K, Willer CJ, Smith AV, et al. Genomic inflation factors under polygenic inheritance. Eur J Hum Genet. 2011;19:807-12.

Tuiskula-Haavisto M, Honkatukia M, Vilkki J, de Koning DJ, Schulman NF, Maki-Tanila A. Mapping of quantitative trait loci affecting quality and production traits in egg layers. Poult Sci. 2002;81:919-27.

Gudbjartsson DF, Walters GB, Thorleifsson G, Stefansson H, Halldorsson BV, Zusmanovich P, et al. Many sequence variants affecting diversity of adult human height. Nat Genet. 2008;40:609-15.

Pryce JE, Hayes BJ, Bolormaa S, Goddard ME. Polymorphic regions affecting human height also control stature in cattle. Genetics. 2011;187:981-4.

Lindholm-Perry AK, Sexten AK, Kuehn LA, Smith TPL, King DA, Shackelford SD, et al. Association, effects and validation of polymorphisms within the NCAPG-LCORL locus located on BTA6 with feed intake, gain, meat and carcass traits in beef cattle. BMC Genet. 2011;12:103.

Thorens B. Glucagon-like peptide-1 and control of insulin secretion. Diabete Metab. 1995;21:311-8.

Lavin JH, Wittert GA, Andrews J, Yeap B, Wishart JM, Morris HA, et al. Interaction of insulin, glucagon-like peptide 1, gastric inhibitory polypeptide, and appetite in response to intraduodenal carbohydrate. Am J Clin Nutr. 1998;68:591-8.

Tuiskula-Haavisto M, de Koning D-J, Honkatukia M, Schulman NF, Mäki-Tanila A, Vilkki J. Quantitative trait loci with parent-of-origin effects in chicken. Genet Res. 2004;84:57-66.

The mouse gene expression data in MGI. http://www.informatics.jax.org/. Accessed 20 Nov 2014.

Hawley SA, Boudeau J, Reid JL, Mustard KJ, Udd L, Makela TP, Alessi DR, et al. Complexes between the LKB1 tumor suppressor, STRADalpha/beta and MO25alpha/beta are upstream kinases in the AMP-activated protein kinase cascade. J Biol. 2003;2:28.

Minokoshi Y, Alquier T, Furukawa N, Kim YB, Lee A, Xue B, et al. AMP-kinase regulates food intake by responding to hormonal and nutrient signals in the hypothalamus. Nature. 2004;428:569-74.

Proszkowiec-Weglarz M, Richards MP, Ramachandran R, McMurtry JP. Characterization of the AMP-activated protein kinase pathway in chickens. Comp Biochem Physiol B Biochem Mol Biol. 2006;143:92-106.

Bolormaa S, Pryce JE, Reverter A, Zhang Y, Barendse W, Kemper K, et al. A Multi-trait, meta-analysis for detecting pleiotropic polymorphisms for stature, fatness and reproduction in beef cattle. PLoS Genet. 2014;10:e1004198.

Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in complex traits: challenges and strategies. Nat Rev Genet. 2013;14:483-95.

Visscher PM, Hill WG, Wray NR. Heritability in the genomics era-concepts and misconceptions. Nat Rev Genet. 2008;9:255-66.

Dascal N. Signalling via the G protein-activated K+ channels. Cell Signal. 1997;9:551-73.

Obici S, Feng Z, Arduini A, Conti R, Rossetti L. Inhibition of hypothalamic carnitine palmitoyltransferase-1 decreases food intake and glucose production. Nat Med. 2003;9:756-61.

Ronnett GV, Kleman AM, Kim EK, Landree LE, Tu Y. Fatty acid metabolism, the central nervous system, and feeding. Obesity (Silver Spring). 2006;14:201S-7S.

Wolfgang MJ, Kurama T, Dai Y, Suwa A, Asaumi M, Matsumoto S-I, Cha SH, Shimokawa T, Lane MD. The brain-specific carnitine palmitoyltransferase-1c regulates energy homeostasis. Proc Natl Acad Sci USA. 2006;103:7282-7.

Ka S, Markljung E, Ring H, Albert FW, Harun-Or-Rashid M, Wahlberg P, et al. Expression of carnitine palmitoyl-CoA transferase-1B is influenced by a cis-acting eQTL in two chicken lines selected for high and low body weight. Physiol Genomics. 2013;45:367-76.

Jin H, Oksenberg D, Ashkenazi A, Peroutka SJ, Duncan AM, Rozmahel R, et al. Characterization of the human 5-hydroxytryptamine1B receptor. J Biol Chem. 1992;267:5735-8.

Levitan RD, Kaplan AS, Masellis M, Basile VS, Walker ML, Lipson N, et al. Polymorphism of the serotonin 5-HT1B receptor gene (HTR1B) associated with minimum lifetime body mass index in women with bulimia nervosa. Biol Psychiatry. 2001;50:640-3.

Saudou F, Amara DA, Dierich A, LeMeur M, Ramboz S, Segu L, et al. Enhanced aggressive behavior in mice lacking 5-HT1B receptor. Science. 1994;265:1875-8.

Koch WJ, Hawes BE, Allen LF, Lefkowitz RJ. Direct evidence that Gi-coupled receptor stimulation of mitogen-activated protein kinase is mediated by G beta gamma activation of p21ras. Proc Natl Acad Sci USA. 1994;91:12706-10.

Chen Y, Gondro C, Quinn K, Herd RM, Parnell PF, Vanselow B. Global gene expression profiling reveals genes expressed differentially in cattle with high and low residual feed intake. Anim Genet. 2011;42:475-90.

Tall A. Plasma lipid transfer proteins. Annu Rev Biochem. 1995;64:235-57.

Sato K, Ohuchi A, Sato T, Schneider WJ, Akiba Y. Molecular characterization and expression of the cholesteryl ester transfer protein gene in chickens. Comp Biochem Physiol B Biochem Mol Biol. 2007;148:117-23.

Yamamoto-Hino M, Sugiyama T, Hikichi K, Mattei MG, Hasegawa K, Sekine S, et al. Cloning and characterization of human type 2 and type 3 inositol 1,4,5-trisphosphate receptors. Receptors Channels. 1994;2:9-22.

Jonchère V, Brionne A, Gautron J, Nys Y. Identification of uterine ion transporters for mineralisation precursors of the avian eggshell. BMC Physiol. 2012;12:10.