- Research Article
- Open Access
Genome-wide association studies for feed intake and efficiency in two laying periods of chickens
- Jingwei Yuan†1,
- Kehua Wang†2,
- Guoqiang Yi1,
- Meng Ma2,
- Taocun Dou2,
- Congjiao Sun1,
- Lu-Jiang Qu1,
- Manman Shen2,
- Liang Qu2 and
- Ning Yang1Email author
© Yuan et al. 2015
- Received: 28 January 2015
- Accepted: 7 October 2015
- Published: 16 October 2015
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).
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.
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.
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.
- Quantitative Trait Locus
- Feed Intake
- Haplotype Block
- Feed Conversion Ratio
- Significant SNPs
The long-term challenge for animal breeders is to improve the productivity of major livestock species to meet the growing demands for livestock products and minimize their impact on the environment and global natural resources . An effective way to face this challenge is to improve feed intake and efficiency that play a decisive role in the economic benefit of livestock husbandry. Therefore, genetic improvement through the selection of animals with a greater ability to use feed could be one of the primary breeding goals. Generally, feed conversion ratio (FCR) and residual feed intake (RFI) were used to measure feed efficiency. FCR is widely used but not a suitable selection trait because of its complex correlations with growth and production traits . RFI is defined as the difference between the observed and expected feed intake given a certain production, which is a sensitive and accurate measure for feed efficiency in chickens. Although traditional selection for RFI has made substantial genetic progress (see review ), it did not maintain the potential of current egg production of layers. However, selection based on genomic information could be a promising alternative .
In the previous decades, quantitative trait loci (QTL) for many traits in chicken have been studied. For feed intake (FI), RFI, and FCR, 37 QTL have been detected on various chicken chromosomes (http://www.animalgenome.org/cgi-bin/QTLdb/GG/index). However, these QTL mapping studies and candidate gene approaches are insufficient because of the low power of linkage analyses and bias in the detection of biologically plausible candidates for complex traits [5, 6].
Compared to previous studies, genome-wide association studies (GWAS) can involve larger genomic regions and detect smaller associated chromosomal regions, and provide more precise estimates of the size and direction of the effects of alleles at the associated loci . Recently, GWAS have revealed many important findings associated with production traits, disease resistance, and morphological characteristics in chickens . However, to date, no GWAS has been performed to evaluate the genetic architecture of feed intake and efficiency in chickens. In addition, the availability of the chicken 600 K high-density single nucleotide polymorphism (SNP) chip permits interrogation of the entire chicken genome at resolutions that were previously unattainable, making the genetic analysis more powerful. The purpose of our study was to identify genetic variants and candidate genes for feed intake and efficiency in chicken using the genome-wide association approach in an F2 resource population genotyped with the chicken 600 K SNP array.
This study was conducted in accordance with the Guidelines for Experimental Animals established by the Animal Welfare Committee of the Chinese Agricultural University.
Birds and phenotypes
An F2 resource population was produced by reciprocal crossing since 2011 between White Leghorn (WL) chicken that originated from the Shanghai Poultry Breeding Co., Ltd, with selection on egg production and quality, and Dongxiang Blue-shelled (DBS) chicken, an indigenous breed introduced from the Jiangxi province, China, with selection on egg production, eggshell color, and egg quality since 1998 at the research farm in the Jiangsu Institute of Poultry Science. WL is a global commercial layer breed that is recognized for its high egg production efficiency, whereas DBS is recognized for its blue eggshell color and meat quality. The two parental lines differed in various characters including morphological, physiological, and production traits. Reciprocal mating between unrelated individuals, i.e., 6 WL (♂) × 133 DBS (♀) and 6 DBS (♂) × 80 WL (♀), was used to produce the F1 generation. Then, 25 F1 males and 406 F1 females from WL (♂) × DBS (♀) and 24 F1 males and 233 F1 females from DBS (♂) × WL (♀) were randomly selected to produce the F2 generation, yielding 1856 males and 1893 females from 590 full-sib families in a single hatch as described previously .
Genotyping, imputation, and quality control
Genomic DNA was extracted from blood samples using the phenol–chloroform method. The eligible DNA was quantified, and then genotyped using the 600 K Affymetrix Axiom high-density array. The SNP chip contained 580,961 SNPs across 28 autosomes, two linkage groups (LGE64 and LGE22C19W28_E50C23), and two sex chromosomes . The Axiom™ GT1 algorithm was used for genotype calling of the 580,961 SNPs from the 1534 samples. The quality control of genotype data was performed with Affymetrix Power Tools (APT) developed by Affymetrix (http://affymetrix.com/): poor quality markers (for which genotyping failed in more than 5 % of the samples) and poor quality DNA samples (for which genotyping failed for more than 3 % of the markers) were removed. The SNP data quality was then controlled by PLINK software  to exclude from further analysis SNPs with a minor allele frequency less than 0.01 and those that deviated from the Hardy–Weinberg equilibrium (HWE) (P value < 1e−6). Finally, 435,867 autosomal SNPs and 1512 samples passed the quality control. To perform association analysis with GEMMA and infer the independent tests with SimpleM, and to increase the power of the GWAS analysis, genotype imputation  was performed. Missing genotypes were imputed based on information from pedigree and remaining SNP genotypes for the F2 population, as implemented in the Beagle Version 4 program with R-squared >0.3 . Quality control of the imputed data was performed using PLINK software, which excluded SNPs that deviated from HWE. Finally, 435,243 SNPs and 1512 samples remained for the GWAS.
All autosomal SNPs were pruned using the indep-pairwise option in PLINK, with a window size of 25 SNPs, a step of five SNPs, and a r2 threshold of 0.2; this resulted in 22,343 independent SNPs. Pairwise identity-by-state (IBS) distances were calculated between all individuals using the independent SNPs, and then multidimensional scaling (MDS) components were obtained using the MDS-plot option based on the IBS matrix. All options and parameters are similar to those previously used  and as implemented in PLINK software.
To further characterize candidate regions that affect traits, we performed linkage disequilibrium (LD) analysis for the chromosomal regions with multiple clustered significant SNPs using the solid spin algorithm implemented in Haploview version 4.2 software . Given the high genetic correlation for each trait between the two laying periods, we performed an association study using a bivariate linear mixed model (bivariate GWAS) , implemented in GEMMA software to detect the genes that were involved in both traits. Finally, a conditional analysis, including the most significant SNPs of a genome-wide scan as covariates into the model in a stepwise manner, was performed to identify secondary association signals .
Since Bonferroni correction is overly conservative due to the high LD in genetic data, it may produce false negative results . Therefore, we calculated the number of effectively independent tests based on the SimpleM method  implemented in R software . It includes three steps: (1) using the cor () function in R to derive the composite LD (CLD) correlation matrix from the SNP dataset; (2) using the R function eigen () to calculate the eigenvalues; and (3) inferring the effective number of independent tests through principal component analysis (PCA), which amounted to 59,286 independent tests. Then, the threshold P-value of the 5 % genome-wide significance was adjusted to 8.43 × 10−7 (0.05/59,286), and 1.69 × 10−5 (1.00/59,286) for suggestive significance.
In addition, SNP-based genetic correlations between traits recorded during the two laying periods, and phenotypic variances explained by the significantly associated SNPs were calculated using GCTA software [24, 25]. The annotated genes that were closest to significant SNPs were identified using Ensembl and NCBI annotation of the Gallus gallus 4.0 genome version.
Descriptive statistics and genetic correlations (rg) between two laying periods for feed intake and efficiency
Univariate and bivariate genome-wide association studies
SNPs significantly associated with daily feed intake detected in both univariate and bivariate GWAS analyses
We performed conditional analysis by fitting the most significant SNP in the model to test for additional associations. The resulting Manhattan and Q–Q plots are in Additional file 2: Figure S1. After conditional analysis, the previously significant SNPs for FI1 on GGA1 disappeared, whereas five SNPs on GGA4 between 75.5 and 76.2 Mb, which were in close proximity to the significant SNP that was detected in the univariate analysis, surpassed the genome-wide significance threshold. These five SNPs were clustered in two blocks by linkage analysis (Fig. 2), and they all had a low MAF (<0.1), thus, we chose SNP rs314523951 previously detected in the univariate analysis as the leading SNP in the following analyses. Interestingly, a region on GGA1 that corresponded to a 30 kb haplotype block between 159.51 and 159.54 Mb included six SNPs (Fig. 2) for FI2 after conditional analysis, which suggested that this region harbors another QTL that affects FI2. Nine SNPs on GGA27 between 3.08 and 3.50 Mb were also significantly associated with FI2, of which eight were in a 418 kb long LD block (Fig. 2). After conditional analysis, genomic control inflation factors decreased to 1.02 for FI1 and 1.03 for FI2, respectively. Moreover, no significant association was found on GGA1 after conditional analysis of the results of the bivariate GWAS, and the genomic control inflation factor was reduced to 1.06.
Leading SNPs associated with feed intake and efficiency in two laying periods
CPV % (se)f
Residual feed intake
SNPs significantly associated with residual feed intake and feed conversion ratio detected in the bivariate GWAS analysis
Feed conversion ratio
Univariate GWAS analysis identified 17 SNPs that were associated with FCR (P value < 1.69e−5), including six SNPs for FCR1 and 11 SNPs for FCR2. Together, the six SNPs for FCR1 explained 3.00 % (SE = 0.02) of its phenotypic variance. Five of these six SNPs were on GGA7 in a 1.94 Mb region between 15.86 and 17.80 Mb, and three of them, i.e. rs317101341, rs317398212, and rs313786079 were in a haplotype block that spanned 324 kb (see Additional file 4: Figure S2). The most significant association was found between FCR1 and SNP rs317101341 (P value = 2.59e−6), which was located at 42.52 kb upstream of the CDCA7 (cell division cycle associated 7) gene (Table 3). Notably, all SNPs on GGA7 had a consistently favorable effect on FCR with a low MAF (<0.1), which indicates that most of the birds in this population carried the dominant genotype. The sixth SNP associated with FCR1 was located on GGA1 at 56.00 Mb in the TBXAS1 (thromboxane A synthase 1 (platelet)) gene. Together the 11 SNPs associated with FCR2 explained 5.14 % (SE = 0.03) of the phenotypic variance. The detailed annotations of these SNPs are in Additional file 3: Table S5, among which SNP rs312620976 with the strongest association signal for FCR2 (P value = 1.86e−6) was situated 16.04 kb downstream of the CETP (cholesteryl ester transfer protein) gene on GGA11 (Table 3).
GWAS is a powerful tool for the genetic analysis of important production traits in farm animals. Recently, GWAS have identified a number of QTL, candidate regions, and SNPs associated with daily FI, RFI, and FCR in beef cattle [27, 28] and swine [29–31]. To our knowledge, the work presented here is the first GWAS of feed intake and efficiency during the laying period in chickens.
Similar to previous findings in dairy cattle , the bivariate analysis detected more associations than the univariate analysis, including genetic variants that affected the same trait expressed in each of the laying periods and genetic variants that affected only expression of a trait in one of the laying periods. These results confirm that multivariate analysis increases the power of detection when genetic correlations occur between traits [32–34].
In practice, it is usually considered that a genomic control inflation factor λ less than 1.05 indicates no population stratification . In the present study, slight population stratification was found in the univariate and bivariate GWAS analyses for FI. Generally, genomic inflation is caused by family structure, cryptic relatedness, and genotyping errors, which lead to spurious associations. However, the powerful mixed model approach that we used here can account for the genomic inflation effectively . Moreover, no strong evidence of population stratification was found for RFI and FCR, for which a limited number of significant genetic variants were detected, which indirectly indicates the polygenic inheritance of FI. Thus, we attributed the observed inflation in our study to the presence of polygenic inheritance , since higher λ are obtained when the number of significantly associated SNPs is large. Furthermore, after conditional analysis, the corresponding associated signals disappeared and λ decreased, which demonstrated that genomic inflation was caused by large number of associated SNPs.
Associations with feed intake
The region associated with FI1 on GGA4 was concordant with a previously reported QTL for a layer population . Interestingly, the five SNPs that were detected by conditional analysis were located within the NCAPG-LCORL (non-SMC condensin I complex, subunit G, and ligand dependent nuclear receptor corepressor-like) locus. This locus is associated with stature in humans, cattle, and horses [37, 38] and is significantly associated with FI in beef cattle . Based on this evidence, it is likely that this locus also affects FI in chickens but further validation on a larger population is necessary. For FI2, the most significant SNP rs316723909 was situated upstream of the GIP (gastric inhibitory polypeptide) gene, which encodes an incretin hormone inducing insulin secretion  and has a role in mediating appetite and energy intake .
The region that we detected on GGA1 for FI was located downstream of a previously reported QTL . The most significant SNP associated with FI1 was within the CAB39L (calcium binding protein 39-like) gene. The protein CAB39L, also known as MO25beta encoding calcium-binding protein 39-like, is expressed in the nervous system and sensory organs in mice . It activates and interacts with the serine/threonine kinase 11/liver kinase (B1STK11/LKB1), which catalyzes the process of phosphorylation to activate AMP-activated protein kinase (AMPK), thus plays a role in the regulation of food intake [44, 45]. Moreover, studies in chickens have suggested that an AMPK pathway similar to that in mammals exists in chickens . Therefore, CAB39L may be a promising candidate gene for FI1, and further functional approaches are necessary to validate this hypothesis. Moreover, the eight significant SNPs that were identified for both FI1 and FI2 on GGA1 had a consistently positive effect, which indicates that this candidate region has a similar effect on daily feed intake in each period and is responsible for the high genetic correlation for feed intake between the two laying periods [7, 47, 48].
Associations with residual feed intake
RFI is a complex trait that is affected by many factors including feed-, growth- and egg-related traits. The changes observed in these traits between the two laying periods lead to different associated genetic variants between RFI1 and RFI2. In addition, environmental conditions can differ between the two laying periods and impact the genetic effect of these variants . This might explain why the identified variants differed between RFI1 and RFI2, and significant SNPs were detected only for RFI2. The region on GGA27 that is associated with RFI2 can be considered as a novel candidate region, since no QTL or SNP was previously reported in this region. The highly significant SNP is located in a novel gene LOC429785 that has no known function in chickens. In humans, LOC429785 is a member of the GIRK (G protein-coupled inwardly-rectifying potassium channel) family and mediates several functions in the central nervous system . Onteru et al.  reported that the GNG4 (guanine nucleotide binding protein 4) gene was associated with RFI in pigs, which acted by activating inward rectifier potassium channels, a mechanism that needs to be further studied in chickens. In the present study, the most significant SNP for RFI1 was in a region for which no gene information was available. The closest annotated gene to this SNP was carnitine palmitoyltransferase 1A (CPT1A) located 21.57 kb downstream. CPT1 is inhibited by malonyl-coenzyme A, which is involved in the hypothalamic regulation of energy expenditure and food intake [51–53]. Ka et al.  demonstrated differential expression of CPT1B in two chicken lines that were selected for high (HWS) and low (LWS) body weight; furthermore, higher levels in the hypothalamus but lower levels in muscle in the HWS compared with that in the LWS line were correlated with increased appetite and food intake.
Unlike the results of the bivariate analysis for FI, the SNPs that were significantly associated with RFI2 on GGA27 did not reach genome-wide significance in the bivariate analysis for RFI. This may be due to a discrepancy between the effects of SNPs on FI and RFI, i.e., the phenotypic variance explained by the SNPs associated with FI on GGA1 was much larger than that explained by the SNPs associated with RFI on GGA27. Four highly significant SNPs were found to be associated with RFI in the bivariate analysis, of which two, rs312873273 and rs312274422, were located at 79.25 Mb on GGA3 and at 1.37 Mb on GGA17, which corresponds to two interesting candidate genes, HTR1B and RGS3, respectively. The HTR1B protein is a G-protein-linked receptor for serotonin  and is associated with eating behavior and appetite control [56, 57]. RGS3 interacts directly with G-proteins and some evidence suggests that RGS family members may indirectly affect proteins in the MAPK (mitogen-activated protein kinase) signal transduction pathways . Furthermore, it has been suggested that MAPK is associated with residual feed intake in pigs and beef cattle [31, 59]. Further research is needed on these two candidate genes in chickens.
Associations with feed conversion ratio
In the univariate GWAS analysis, no SNP reached the genome-wide significance level for FCR1 and FCR2, which indicated that FCR was a complex trait controlled by several or numerous genes, each with a small effect. Thus, a larger sample size is needed to detect associations with FCR. In addition, applying a looser statistical significance threshold by using a false discovery rate method or a Bayesian model may identify more interesting candidate genes. Nevertheless, based on genome information for the significant SNPs, we were able to suggest several promising candidate genes. For FCR2, the most significant SNP rs312620976 is within a promising candidate gene CETP, which is involved in the transfer of cholesteryl ester from high-density lipoproteins to other lipoproteins . Previously, Sato et al. reported that impaired CETP activity played an important role in egg production of laying hens , and egg production is one of important components affecting FCR . Moreover, ITPR2, a promising candidate gene that was identified in the bivariate analysis, encodes a receptor for inositol triphosphate and a calcium channel . Since the inositol trisphosphate receptor type 2 is involved in ion transfer in the hen’s uterus , it is possible that this gene will impact FCR by affecting egg formation.
This study identified eight genomic regions that are significantly associated with FI. Seven of these associations are reported for the first time, which demonstrates the power of high-density SNP arrays and distinct parental lines. One of the significant allele substitution effects was observed for a SNP associated with RFI in the laying period between 57 and 60 weeks and caused a 3.35 g/day phenotypic difference. It is located in LOC429785, and is a candidate gene that possibly exerts a potential function in the control of RFI. Furthermore, we identified the CETP and ITPR2 genes that influenced FCR by affecting egg production. Our results provide valuable knowledge on SNP and candidate genes that are involved in the genetic architecture of FI, RFI, and FCR in two laying periods of chickens. The SNP or regions that are associated with feed intake and efficiency can be used as fundamental information in marker-assisted or genomic selection.
JY contributed to genomic DNA extraction, performed GWAS analyses and drafted the manuscript. KW designed the study and contributed to population construction, phenotype collection and data analysis. GY contributed to statistical analysis and helped to draft the manuscript. MM, TD, CS, LJQ, MS and LQ contributed to population construction, collection of phenotypic data, and statistical analysis. NY supervised and designed the study and contributed to drafting the manuscript. All authors read and approved the final manuscript.
The authors gratefully acknowledge staff at Jiangsu Institute of Poultry Science for their technical assistance. The work was supported in parts by the National High Technology Development Plan of China (2013AA102501), National Scientific Supporting Projects of China (2011BAD28B03), China Agriculture Research Systems (CARS-41), and Program for Changjiang Scholars and Innovative Research in University (IRT1191).
The authors declare that they have no competing interests.
Open AccessThis 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.
- Hume DA, Whitelaw CBA, Archibald AL. The future of animal production: improving productivity and sustainability. J Agr Sci. 2011;149:9–16.View ArticleGoogle Scholar
- 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.Google Scholar
- Flock DK. Genetic-economic aspects of feed efficiency in laying hens. World Poultry Sci J. 1998;54:225–39.View ArticleGoogle Scholar
- Meuwissen THE, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157:1819–29.PubMed CentralPubMedGoogle Scholar
- Tabor HK, Risch NJ, Myers RM. Candidate-gene approaches for studying complex genetic traits: practical considerations. Nat Rev Genet. 2002;3:391–7.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Mackay TFC, Stone EA, Ayroles JF. The genetics of quantitative traits: challenges and prospects. Nat Rev Genet. 2009;10:565–77.View ArticlePubMedGoogle Scholar
- Zhang H, Wang Z, Wang S, Li H. Progress of genome wide association study in domestic animals. J Anim Sci Biotechnol. 2012;3:26.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Luiting P, Urff EM. Optimization of a model to estimate residual feed consumption in the laying hen. Livest Prod Sci. 1991;27:321–38.View ArticleGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Marchini M, Howie B. Genotype imputation for genome-wide association studies. Nat Rev Genet. 2010;11:499–511.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012;44:821–4.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.View ArticlePubMedGoogle Scholar
- Zhou X, Stephens M. Efficient multivariate linear mixed model algorithms for genome-wide association studies. Nat Methods. 2014;11:407–9.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Team RC. R: A language and environment for statistical computing. Vienna: R foundation for statistical computing; 2012. ISBN 3-900051-07-0.Google Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Ng PC, Henikoff S. SIFT: predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31:3812–4.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Kim S, Xing EP. Statistical estimation of correlated genome associations to a quantitative trait network. PLoS Genet. 2009;5:e1000587.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu W, Zhang H. Why do we test multiple traits in genetic association studies? J Korean Stat Soc. 2009;38:1–10.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Pryce JE, Hayes BJ, Bolormaa S, Goddard ME. Polymorphic regions affecting human height also control stature in cattle. Genetics. 2011;187:981–4.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Thorens B. Glucagon-like peptide-1 and control of insulin secretion. Diabete Metab. 1995;21:311–8.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Solovieff N, Cotsapas C, Lee PH, Purcell SM, Smoller JW. Pleiotropy in complex traits: challenges and strategies. Nat Rev Genet. 2013;14:483–95.PubMed CentralView ArticlePubMedGoogle Scholar
- Visscher PM, Hill WG, Wray NR. Heritability in the genomics era–concepts and misconceptions. Nat Rev Genet. 2008;9:255–66.View ArticlePubMedGoogle Scholar
- Dascal N. Signalling via the G protein-activated K+ channels. Cell Signal. 1997;9:551–73.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Tall A. Plasma lipid transfer proteins. Annu Rev Biochem. 1995;64:235–57.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar