Printer Friendly

Genetic architecture and candidate genes detected for chicken internal organ weight with a 600 K single nucleotide polymorphism array.


The internal organs are involved in many biological occurrences, including oxygen transportation, lipid metabolism, and digestion. Each organ has its own developmental process and gene expression profile [1]. Liver development affects abdominal fat, yolk deposit [2] which are all associated with economic performance, and provides a model of non-alcoholic steatohepatitis [3]. Heart is an organ for carrying oxygen [4], affect well-being of the chicken [5]. Physical and chemical digestion takes place in the proventriculus and gizzard, which regulate feed intake and energy balance [6]. Efficient utilization of feed is essential to reduce feed cost and improve economic benefit. Moreover, internal organ like gizzard and heart provide proximate composition and amino acid same as chicken meat [7]. In order to support an economical viable meat system, it is essential to efficient utilize these byproducts. Hence, improve these internal organs will not only improve well-being of the organism but also take full advantage of byproduct from chicken.

With the development of deep sequencing and statistical methods, it is possible to detect the quantitative trait loci (QTLs) of complex traits. In recent years, great advances have been gained in the selection of economic traits in chickens with molecular biology techniques. Nones et al [8] identified novel QTLs associated with gizzard, liver, heart, lung, and thigh located on Gallus gallus autosome (GGA) 1. Rosario et al [9] found novel QTLs on GGA1, GGA3, and GGA4 related with carcass traits. Boschiero et al [10] identified microsatellite markers associated with performance and carcass traits on GGA1 and GGA13. Most studies on chicken internal organ weight have focused on the early stages of development, which has hindered the progress of similar research in the layer hen, especially for spent layer hens in China.

In this study, we conducted a genome-wide association study (GWAS) on an F2 resource population generated from Dongxiang Blue-Shelled (DX) and White Leghorn (WL) chickens with a 600 K Affymetrix chip. Furthermore, a linkage disequilibrium (LD) analysis, variance estimated percentage, and gene annotation were performed by biotechnology. This research will provide the genetic architecture that underlies the weight of internal organs in layer chickens.


Ethical statement

Animal procedures and sample collections were carried out following the guidelines established by the Ministry of Agriculture of China. The Institution of Animal Care and Use Committee in the Poultry Institute, Chinese Academy of Agricultural Science, Yangzhou, China approved the procedures.

Phenotype collection and analyses

We constructed an [F.sub.2] resource population by reciprocal crosses between WL and DX chickens. The WL was from a commercial line and DX from a Chinese indigenous strain, therefore the two lines represent two very different genetic pools, including morphological, physiological, and production traits, such as feed efficiency [11], follicle number [12], and the yolk proportion [13]. The F2 resource population of 1512 birds from 49 half-sib and 590 full-sib families were produced from a WL/DX (25 [male]:407 [female]) and DX/WL (24 [male]:235 [female]) cross in the [F.sub.1] generation. The [F.sub.1] population was generated from WL (6 [male])xDX (133 [female]) and DX (6 [male])xWL (80 [female]) by initial reciprocal crosses, producing 1,029 and 552 chicks, respectively. The experimental animals were kept indoors under standardized conditions at the research base in the Jiangsu Institute of Poultry Science, Yangzhou, China. Birds had ad libitum access to feed and water that met all NRC requirements and were housed in single-hen cages under a 16 L:8 D lighting regime. The [F.sub.2] birds were slaughtered at 72 weeks; immediately after the slaughter we measured the weights of heart (HW), liver (LW), proventriculus (PW), and gizzard (GW). DNA was collected by standard venipuncture. Data including a descriptive phenotype and normality test statistics were summarized using R software.

Genotyping and quality control

The standard phenol/chloroform method was used to extract genomic DNA, which was genotyped against a 600 K Affymetrix Axiom Chicken Genotyping Array (Affymetrix, Inc. Santa Clara, CA, USA). Affymetrix Power Tools v1.16.0 (APT) ( software was used to analyze genotype calling and quality control (QC). A total of 1,512 hens and 532,299 single nucleotide polymorphisms (SNPs) remained valid after the application of APT for QC. An R script was used to calculate the SNP QC metrics and filter out individual SNPs falling below given thresholds. The effectiveness of the detecting quality was tested by PLINK v1.90 [14]. We imputed for some sporadic missing genotypes using the BEAGLE v4.0 package [15]. Briefly, samples with a dish QC of more than 0.82, call rate less than 97%, minor allele frequency (MAF) higher than 5%, Hardy-Weinberg equilibrium test P of more than 1xE[10.sup.-6] were kept for further analysis. After QC, 1,512 individuals and 435,867 SNPs were used in further GWA analysis.

Genome-wide association analysis

A principal component analysis was used to eliminate spurious associations resulting from the presence of cryptic relatedness or hidden population stratification. Then, the full SNP set to 41,130 independent SNPs were pruned via the -indep-pairwise 25 5 0.2 command (PLINK), and we included five principal components (PCs) as covariates in the mixed model. After simpleM was used to determine the threshold by correcting the number of multiple tests, the genome-wide significant and suggestive p-values were 8.43x[10.sup.-7] (0.05/59,308) and 1.69x[10.sup.-5] (1.00/59,308), respectively.

The GEMMA v0.94 package [16], with the exact mixed model approach, was implemented with the valid individuals and SNPs for univariate analysis. The independent SNPs were used to compute the centered relatedness matrix, and the significance p-value level between SNPs and phenotypes was calculated from the derived Wald test. The univariate linear mixed model was denoted as:

y = W[alpha]+x[beta]+u+[epsilon]

Where y represents an n-vector of phenotypic values for n individuals; W is a matrix of covariates (fixed effects with a column of 1s and top five PCs), [alpha] is a vector of the corresponding coefficients including the intercept; x is an n-vector of the genotypes of the SNP marker, [beta] is the effect size of the marker; [mu] is an n-vector of random effects; and e is an n-vector of errors.

The Manhattan plots and quantile-quantile (QQ) plots were drawn by "gap" packages in R software. The genomic inflation factor was calculated by using the GenABEL package in the R project.

Linkage disequilibrium analysis and gene annotation

The LD were implemented by Haploview v4.2 to test the association of significant SNPs with each other. When the distance between two consecutive genome-wide significant SNPs was greater than 10 Mb, they were considered as two separate QTL. The LD analyses would find a haplotype blocks because of these loci located in a strong LD region. But, a genuine causal locus could not be distinguished by present GWAS method. Therefore, a functional gene annotation on these significant SNPs were essential performed to characterize potential candidate genes. In our study, we searched candidate genes based on Galgal5 assembly, using BioMart system supported by Ensemble ( The genes nearest or harboring significant SNPs of internal organ weight were chosen as candidate locations.

SNP effects and chromosome heritability

The contributions to phenotypic variance (CPV) explained by the significant association SNPs was calculated with GCTA software [17] using univariate restricted maximum likelihood. The heritability explained by the eligible SNPs ([h.sup.2.sub.snp]) for GWAS or the variance contributed by each chromosome were also estimated with GCTA.


Phenotypic statistics and genetic parameters

Descriptive data for each organ are displayed in Table 1. After rank-based inverse normal transformation, all phenotypic values conformed to a normal distribution. The range of all internal organ is large, and with the coefficient of variance more than 19%. Table 2 displayed the genetic analysis, which showed that the heritability of internal organs was moderate to high, and the GW had the highest heritability (0.640). Genetic correlation between each organ weight was higher than 0.5, the highest genetic correlation was PW and GW (0.711).

Genome-wide association study results

The results of the GWAS were as follows. There are 380 significant SNPs on GGA1, 49 SNPs on GGA4 were related to internal organ weights, and six SNPs on GGA2 were related to HW (Table 2, Figure 1, Supplementary Figure S1, S2, S3). The regions ranged from 165.63 to 173.78 Mb on GGA1, and from 71.30 to 77.19 Mb on GGA4. After Venn diagram analysis, 17 and 6 significant SNPs located on GGA1 and GGA4, respectively, were found to affect all internal organ weights (Figure 2).

Estimation of contributions to phenotypic variance and plausible candidate genes

By using gene annotation analysis, it was obtained that the 17 significant SNPs on GGA1 were harbored in nine genes, general transcription factor IIF polypeptide 2 (GTF2F2), potassium channel tetramerization domain containing 4, spermatid associated, succinate-coA ligase ADP-forming [beta] subunit, 5-hydroxytryptamine receptor 2A (HTR2A), WD repeat and FYVE domain containing 2 (WDFY2), asparagine-linked glycosylation 11, NIMA related kinase 3, and cytoskeleton associated protein 2. The SNP on GGA1 possessing a MAF of more than 0.4 and the values of beta on GGA1 were all negative (Table 4). We then analyzed the allelic contribution to phenotypic variation (Table 4), the results show that the smaller the p value, the larger the CPV value. All significant SNPs on GGA1 explained over 5% of CPV for PW and GW, while the CPV for HW and LW was from 2.781% to 4.723%. The CPV of SNPs rs312726815, rs315120631, rs13972990 were highest for HW, LW, and PW/GW. The rs312726815 SNP located in the upstream (16,831 bp) of HTR2A, rs315120631 located on the intron of GTF2F2, rs13972990 was downstream (95,255 bp) of WDFY2.

Table 5 displayed that the SNPs on GGA4 possess a MAF more than 0.05 and less than 0.06, with values of beta being positive. These SNPs located in or near genes ligand dependent nuclear receptor corepressor like, non-SMC condensin I complex subunit G (NCAPG), leucine aminopeptidase 3, LIM domain binding 2, quinoid dihydropteridine reductase, prominin 1, and transmembrane anterior posterior transformation 1. We found one locus rs14491030 that located in the missense of the gene NCAPG and another SNP close to NCAPG.

Table 6 showed that the SNPs on GGA2 association with HW association with HW displayed a MAF of ~0.3, with the values of beta being positive. These SNPs located on the intron of sonic hedgehog (SHH) and close to genes 5-Hydroxytryptamine receptor 5A, insulin induced gene 1, and limb development membrane protein. The SNPs rs316413705 explained the highest CPV (3.109%) for HW.

Linkage disequilibrium analysis

The LD analysis was performed because of the potentially strong LD between neighboring variants. The LD analysis results revealed that significant SNPs on GGA1, GGA4, or GGA2 were all extremely strong in a LD status (Figure 3). Combined the aforementioned CPV calculation analysis, LD analysis and gene annotation, we considered rs312726815, rs315120631, and rs13972990 on GGA1 as primary candidate loci, with correspond to genes HTR2A, GTF2F2, and WDFY2 associated with HW LW and PW/GW respectively. The NCAPG gene harbored the missense locus rs14491030 on GGA4 were considered as a candidate gene for all internal organs, SHH gene with intron harbored the SNP rs14491030 on GGA2 had an association with HW.

Genome partitioning of genetic variation

The genetic architecture of internal organs was further illustrated by partitioning the genetic variation onto chromosome segments with an exploratory analysis. HW and PW in the joint model could not converge due to the relatively small sample size in the F2 population. Only the partitioning spectrum of LW and GW were estimated. The estimates of variance contributed by each chromosome exhibited a medium and strong linear relationship with the length of the chromosome for LW ([R.sup.2] = 0.493, 0.493 Figure 4A) and GW ([R.sup.2] = 0.753, 0.493 Figure 4B). For GGA1 this explained 18.29% and 10.03% of the phenotypic variance for GW and LW, respectively. While for GGA4 it explained 11.30% and 4.09% for GW and LW, respectively.


To accomplish the goal of "500 eggs in 100 weeks", it is essential to continuous select egg production trait as well as fitness trait (such as increasing heart capacity [18]) for which populations were not routinely selected. Because of the correlation between egg production and fitness showed positive or negative was not absolute, it is essential to detect molecular markers associated with economic traits as well as fitness traits to improve animals simultaneously. Moreover, some body composition trait like gizzard and proventriculus indirectly reflects feed efficiency [5]. Hence, improving these traits will not only make chicken well-being better but also take advantage of these organs. But, these traits are inconvenient or difficult to measure for animals should be killed. Application of DNA markers could simultaneously select economic trait and these traits.

Our results showed that phenotypic data displayed a large variation for all internal weights, which was probably because these traits were not chosen as selection index. The heritability showed a moderate to high, higher than that the previous reports. For LW, the heritability estimates were 0.31 [19] and 0.33 [20] at 6 weeks, and the heritability of HW and GW was 0.27 and 0.44, respectively, in a report by Venturini et al [20]. The difference of the heritability may be due to different breeds and ages [21]. In addition, we dissected genome partitioning of genetic variation for GW and LW. A positive linear correlation between the variance explained by each chromosome and its length were obtained in our work, which corresponds to previous reports [22]. In particular, GGA1 accounted for the largest genetic variance (18.29%) for GW, which is consistent with the highest heritability (0.640), whereas GGA1 accounted for 10.03% genetic variance for LW, the heritability of which is moderate (0.355). We then considered that with the higher heritability of the trait, the genetic variance of GGA will also be higher. This finding is in line with previous research especially in our resource population [23].

Using GWAS method, we totally found 17 SNPs on GGA1 and 6 SNPs on GGA4 associated with all internal organ weights, and 6 SNPs on GGA2 were associated with HW. Previous studies showed that QTL for internal organ weight, such as HW and LW, at 9 weeks was detected on GGA1 and GGA4, respectively [24]. HW and GW at 6 weeks were mapped to GGA13 [10]. The difference between the present study and previous reports might be due to differences in the age and population of the birds used in different studies. These studies used younger birds, the age of the birds in the present study was 72 weeks. In addition, one trait, such as internal organ weight, was controlled by more than one QTL, which will help us understand the genetic architecture underlying the quantitative traits that are controlled by polygenes.

The most significant loci in the GGA1 spanned from 166.99 to 169.66 Mb, which is also related to comb weight [25], egg weight [23], ovary weight [26], and feed intake [27] in our previous studies and growth traits in other reports [28-30]. The internal organs as fitness traits participate in metabolism that affects growth performance and economic traits, therefore the selection of body weight or egg weight will indirectly result in internal organ weight change. This has important ramifications for understanding the pleiotropic effects of the locus or gene [31].

After rigorous statistic and LD analysis, three SNPs rs3127 26815, rs315120631, and rs13972990, correspond to genes HTR2A, GTF2F2, and WDFY2, respectively, were selected as the candidate loci. HTR2A was considered as a candidate gene for HW. The mRNA of the gene, which encodes a receptor of 5-hydroxytryptamine (serotonin), is expressed within vascular smooth muscle, endothelial cells, and cardiomyocytes, and plays a role in vasoconstriction and cellular proliferation [32]. Also, the function of the HTR2A-encoded protein is associated with blood pressure and heart rate and is a contributing factor to cardiovascular disease in humans [33]. In birds, the size of heart was association with the rate of metabolism and indicates the capacity to move the blood [34]. Therefore, it is reasonable to propose that HTR2A was responsible for heart weight, the detailed information needs to be further investigated. Second gene GTF2F2 nearby SNP rs315120631 was considered as candidate gene for LW. GTF2F2 forms a heteromeric general transcription initiation factor [35]. GTF2F2 widely exists in many tissues and organs, especially in the liver, lung, and kidney. The liver is the main organ regulating body homeostasis [36], where a large amount of active transcription and translation takes place. Reports said that GTF2F2 is association with rat liver regeneration [37]. Liver is responsible for detoxification and energy balance, which involved in various of biological process including transcription and translation. Hence, further studies about the function of GTF2F2 in chicken internal organs are needed. The third SNP, rs13972990, is near gene WDFY2, which is associated with PW and GW WDFY2 (WD repeat and FYVE domain containing 2) is an endosomal protein, modulating the PI3K/AKT pathway, which is known to be involved in oncogenesis [38]. WDFY2 overexpression can increase adipogenesis, which may play a role in metabolic disorders such as diabetes in humans [39]. Proventriculus and gizzard modulated feed behavior [6] and feed overconsumption [5], indirectly regulating glycogen synthesis. It is tempting to speculate that WDFY2 may participate in proventriculus and gizzard growth though its involvement is not clear.

Moreover, the region on GGA4 was related to egg weight, eggshell weight, oviduct weight in our population [23,40]. NCAPG harbored the missense rs14491030 was a candidate gene for all internal weight. The gene was related to residual feed intake in cattle [41] and withers height in horses [42] in previous reports. Most breeders focus on egg production in layer chicken, resulting indirectly in selection for fitness trait or digestive organs. NCAPG was considered as candidate gene for all internal organ weight though no selection for these traits, pleiotropic effects of the locus on other traits that are under selection, or close linkage and LD with QTL that are under selection [43]. We provide an evidence that quantitative traits that may be controlled by polygenes and a single gene or mutation may lead to a host of alterations in multiple traits.

In addition, GWAS results showed that one QTL on GGA2 associated with HW. The SHH could be considered as a candidate gene. The SHH gene is one of the Hh proteins, which plays a crucial role in the development of all animals and regulates morphogenesis of a variety of tissues and organs in the embryo [44]. Previous studies showed that SHH is necessary for secondary heart field proliferation in humans [45]. Consequently, the function of SHH in the chicken heart may involve heart differentiation.


In present study, a GWAS strategy was performed to detect potential QTLs or genes association with internal organ traits. Our study provides evidence that the internal organ weight at 72-weeks old appears to moderate a high heritability and may share the similar genetic mechanisms. Five candidate genes were identified with significant effect on internal weights. These promising loci or genes could be helpful to simultaneous improve layer chicken economic trait and fitness trait to accomplish sustainable use of the chicken.


We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.


This work was supported in part by China Agriculture Research Systems (grant no. CARS-40-K01), Program for the Agricultural Major New Variety of the Jiangsu Province (grant no. PZCZ201729), National Natural Science Foundation of China (grant no. 31601938) and the Sanxin Project of Jiangsu Province, (grant no. SXGC[2017]252). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The authors also gratefully acknowledge the staff of the China Agricultural University for involvement in this program.


[1.] He Y, Li X, Zhang F, et al. Multi-breed genome-wide association study reveals novel loci associated with the weight of internal organs. Genet Sel Evol 2015;47:87.

[2.] Gruber M, Bos ES, Ab G. Hormonal control of vitellogenin synthesis in avian liver. Mol Cell Endocrinol 1976;5:41-50.

[3.] Ayala I, Castillo AM, Adanez G, et al. Hyperlipidemic chicken as a model of non-alcoholic steatohepatitis. Exp Biol Med (Maywood) 2009;234:10-6.

[4.] Zhang LF, Lian LS, Zhao CJ, et al. Expression pattern of HIF1alpha mRNA in brain, heart and liver tissues of Tibet chicken embryos in hypoxia revealed with quantitative real-time PCR. Animal 2007;1:1467-71.

[5.] Svihus B. The gizzard: function, influence of diet structure and effects on nutrient availability. Worlds Poult Sci J 2011;67:207-24.

[6.] Richards MP. Genetic regulation of feed intake and energy balance in poultry. Poult Sci 2003;82:907-16.

[7.] Maiti AK, Ahlawat SS. Effect of natural tenderizers on physicochemical properties of chicken gizzard and goat heart. Am J Food Technol 2011;6:80-6.

[8.] Nones K, Ledur M, Ruy D, et al. Mapping QTLs on chicken chromosome 1 for performance and carcass traits in a broiler x layer cross. Anim Genet 2006;37:95-100.

[9.] Rosario MFd, Ledur MC, Moura ASAMT, Coutinho LL, Garcia AAF. Genotypic characterization of microsatellite markers in broiler and layer selected chicken lines and their reciprocal F1s. Sci Agric 2009;66:150-8.

[10.] Boschiero C, Rosario M, Ledur M, et al. Associations between microsatellite markers and traits related to performance, carcass and organs in chickens. Int J Poult Sci 2009;8:615-20.

[11.] Yuan J, Dou T, Ma M, et al. Genetic parameters of feed efficiency traits in laying period of chickens. Poult Sci 2015;94:1470-5.

[12.] Shen M, Sun H, Qu L, et al. Genetic architecture and candidate genes identified for follicle number in chicken. Sci Rep 2017;7: 16412.

[13.] Jun G, Kehua W, Liang Q, et al. Random regression analyses to model the longitudinal measurements of yolk proportions in the laying chicken. Poult Sci 2017;96:3852-7.

[14.] Purcell S, Neale B, Todd-Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559-75.

[15.] Browning BL, Browning SR. A unified approach to genotype imputation and haplotype-phase inference for large data sets of trios and unrelated individuals. Am J Hum Genet 2009;84: 210-23.

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

[17.] Yang J, Lee S, Goddard M, Visscher P. GCTA: a tool for genomewide complex trait analysis. Am J Hum Genet 2011;88:76-82.

[18.] Deeb N, Lamont SJ. Genetic architecture of growth and body composition in unique chicken populations. J Hered 2002;93: 107-18.

[19.] Nunes Bdo N, Ramos SB, Savegnago RP, et al. Genetic parameters for body weight, carcass chemical composition and yield in a broiler-layer cross developed for QTL mapping. Genet Mol Biol 2011;34:429-34.

[20.] Venturini GC, Cruz VA, Rosa JO, et al. Genetic and phenotypic parameters of carcass and organ traits of broiler chickens. Genet Mol Res 2014;13:10294-300.

[21.] Nitsan Z, Dunnington EA, Siegel PB. Organ growth and digestive enzyme levels to fifteen days of age in lines of chickens differing in body weight. Poult Sci 1991;70:2040-8.

[22.] Gao G, Li Q, Zhao X, et al. Transcriptome profiling of the hypothalamus during prelaying and laying periods in Sichuan white geese (Anser cygnoides). Anim Sci J 2015;86:800-5.

[23.] Yi G, Shen M, Yuan J, et al. Genome-wide association study dissects genetic architecture underlying longitudinal egg weights in chickens. BMC Genomics 2015;16:746.

[24.] Navarro P, Visscher P, Knott S, et al. Mapping of quantitative trait loci affecting organ weights and blood variables in a broiler layer cross. Br Poult Sci 2005;46:430-42.

[25.] Shen M, Qu L, Ma M, et al. Genome-wide association studies for comb traits in chickens. PLoS One 2016;11:e0159081.

[26.] Sun C, Lu J, Yi G, et al. Promising loci and genes for yolk and ovary weight in chickens revealed by a genome-wide association study. PloS one 2015;10:e0137145.

[27.] Yuan J, Wang K, Yi G, et al. Genome-wide association studies for feed intake and efficiency in two laying periods of chickens. Genet Sel Evol 2015;47:82.

[28.] Zhang H, Zhang YD, Wang SZ, et al. Detection and fine mapping of quantitative trait loci for bone traits on chicken chromosome one. J Anim Breed Genet 2010;127:462-8.

[29.] Liu X, Zhang H, Li H, et al. Fine-mapping quantitative trait loci for body weight and abdominal fat traits: effects of marker density and sample size. Poult Sci 2008;87:1314-9.

[30.] Wang W, Wang J, Zhang T, et al. Genome-wide association study of growth traits in Jinghai Yellow chicken hens using SLAF-seq technology. Anim Genet 2015 Sept 13 [Epub]. https://

[31.] Hodgkin J. Seven types of pleiotropy. Int J Dev Biol 1998;42:5015.

[32.] Nagatomo T, Rashid M, Muntasir HA, Komiyama T. Functions of 5-HT 2A receptor and its antagonists in the cardiovascular system. Pharmacol Ther 2004;104:59-81.

[33.] Nichols CD. Serotonin 5-HT2A receptor function as a contributing factor to both neuropsychiatric and cardiovascular diseases. Cardiovasc Psychiatry Neurol 2009;2009:475108.

[34.] Hartman FA. Heart weight in birds. Condor 1955;57:221-38.

[35.] Sopta M, Burton ZF, Greenblatt J. Structure and associated DNA-helicase activity of a general transcription initiation factor that binds to RNA polymerase II. Nature 1989;341:410-4.

[36.] Milagro FI, Campion J, Martinez JA. Weight gain induced by high-fat feeding involves increased liver oxidative stress. Obesity (Silver Spring) 2006;14:1118-23.

[37.] Hu K, Xu C. Expression patterns analysis of the RNA biosynthesis, processing, transporting or catabolism-associated genes during rat liver regeneration. Life Sci J 2008;5:1-8.

[38.] Kannan K, Coarfa C, Rajapakshe K, et al. CDKN2D-WDFY2 is a cancer-specific fusion gene recurrent in high-grade serous ovarian carcinoma. PLoS Genet 2014;10:e1004216.

[39.] Fritzius T, Moelling K. Akt- and Foxo1-interacting WD-repeat-FYVE protein promotes adipogenesis. EMBO J 2008;27:1399410.

[40.] Shen M, Qu L, Ma M, et al. A genome-wide study to identify genes responsible for oviduct development in chickens. PLoS One 2017;12:e0189955.

[41.] Widmann P, Reverter A, Weikard R, et al. Systems biology analysis merging phenotype, metabolomic and genomic data identifies Non-SMC Condensin I Complex, Subunit G (NCAPG) and cellular maintenance processes as major contributors to genetic variability in bovine feed efficiency. PLoS One 2015;10: e0124574.

[42.] Tetens J, Widmann P, Kuhn C, Thaller G. A genome-wide association study indicates LCORL/NCAPG as a candidate locus for withers height in German Warmblood horses. Anim Genet 2013;44:467-71.

[43.] Abasht B, Dekkers J, Lamont S. Review of quantitative trait loci identified in the chicken. Poult Sci 2006;85:2079-96.

[44.] Balaskas N, Ribeiro A, Panovska J, et al. Gene regulatory logic for reading the Sonic Hedgehog signaling gradient in the vertebrate neural tube. Cell 2012;148:273-84.

[45.] Dyer LA, Kirby ML. Sonic hedgehog maintains proliferation in secondary heart field progenitors and is required for normal arterial pole formation. Dev Biol 2009;330:305-17.

Taocun Dou (1), (a), Manman Shen (1,2), (a), Meng Ma (1), Liang Qu (1,3), Yongfeng Li (1), Yuping Hu (1), Jian Lu (1), Jun Guo (1), Xingguo Wang (1), and Kehua Wang (1) *

* Corresponding Author: Kehua Wang

Tel: +86-38-0527-6606, Fax: +86-0514-8559-9012, E-mail:

(1) Jiangsu Institute of Poultry Science, Chinese Academy of Agricultural Sciences, Yangzhou, Jiangsu 225216, China

(2) College of Animal Science and Technology, Yangzhou University, Yangzhou, Jiangsu 225009, China

(3) College of Animal Science and Technology, Nanjing Agriculture University, Nanjing, Jiangsu 210095, China

(a) Taocun Dou and Manman Shen have contributed equally to this study.


Kehua Wang

Submitted Apr 9, 2018; Revised May 30, 2018; Accepted Jul 23, 2018

Caption: Figure 1. Manhattan plot and QQ-plot of genome-wide association analysis for heart weight. The left plot is the Manhattan plot, which shows the -log10 (observed p values) for association of single nucleotide polymorphisms (y-axis) plotted against their chromosomal positions on each chromosome (x-axis). The right plot is the QQ test for population structure, the x-axis indicates the expected -log10-transformed p values, and the y-axis shows the observed -log10-transformed p values.

Caption: Figure 2. Venn diagram of significant single nucleotide polymorphisms on GGA1 (A) and GGA4 (B) associated with four internal organ weights by univariate association test. Heart weight, liver weight, proventriculus weight, and gizzard weight are abbreviated as HW, LW, PW, and GW, respectively. GGA, Gallus gallus autosome.

Caption: Figure 3. LD plot for significant single nucleotide polymorphisms at GGA1, GGA4, and GGA2. Plot A represents LD analysis on GGA1, plot B represents LD analysis on GGA4, plot C represents LD analysis on GGA2. LD, linkage disequilibrium; GGA, Gallus gallus autosome.

Caption: Figure 4. Genome partitioning for LW and GW by joint analysis. The estimated proportion of variance captured by each chromosome against its size. The characters in the circles are the chromosome numbers. Plot A represents liver weight, plot B represents gizzard weight. LW, liver weight; GW, gizzard weight.
Table 1. Descriptive statistics for internal organ weight in the
[F.sub.2] population

Trait        Mean        SD        Min        Max       CV (%)

HW           5.42       1.03       3.00      12.20      19.00
LW          28.18       6.68      12.40      66.80      23.69
PW           1.64       0.79       0.10      14.30      48.36
GW           5.25       1.21       1.00      15.10      22.95

Mean, arithmetic mean; SD, standard deviation; Min, minimum; Max,
maximum; CV, coefficient of variance; HW, heart weight; LW, liver
weight; PW, proventriculus weight; GW, gizzard weight.

Table 2. Summary of genetic analysis for internal organs (1)

Trait         HW              LW              PW              GW

HW       0.283(0.043)    0.658(0.079)    0.666(0.043)    0.567(0.073)
LW           0.351       0.355(0.045)    0.645(0.075)    0.503(0.071)
PW           0.360           0.240       0.408(0.044)    0.711(0.049)
GW           0.389           0.358           0.475       0.640(0.037)

HW, heart weight; LW, liver weight; PW, proventriculus weight; GW,
gizzard weight.

(1) Diagonal, heritability estimates; Lower triangle, phenotypic
correlations; Upper triangle, genetic correlations. Standard errors of
the estimates are in parentheses.

Table 3. Genome-wide significant regions identified by genome-wide
association study for the weights of internal organ

Trait             GGA    No. snp     Region (Mb)     Trait     GGA

HW                 1        97      165.99-171.60      PW       1
                   2        6         7.79-8.34                 4
                   4        42       71.30-77.19
LW                 1        18      167.00-169.72      GW       1
                   4        7        75.48-76.07                4
Total              1       380      165.63-173.78               4
Venn diagram       1        17      166.99-169.66               4

Trait            No. snp     Region (Mb)

HW                 230      165.63-173.78
                    15       75.48-76.21

LW                 308      165.93-172.12
                    32       71.30-76.67
Total               49       71.30-77.19
Venn diagram        6        75.48-76.07

GGA, Gallus gallus autosome; SNP, single nucleotide polymorphisms; No.
snp, the number of significant SNPs; HW, heart weight; LW, liver
weight; PW, proventriculus weight; GW, gizzard weight.

Table 4. Information of 17 SNPs on GGA1 and genomic regions to
internal organ weight

SNP              Position (bp)     Location (bp)       Candidate/
                                                      nearest gene

rs314001986       166,994,822        intron/U-3          GTF2F2
rs15493356        166,998,353          intron            GTF2F2
rs31 5221322      167,002,252          intron            GTF2F2
rs31 51 20631     167,009,025          intron            GTF2F2
rs15493542        167,042,503          intron            GTF2F2
rs314998864       167,046,824          intron            GTF2F2
rs312642102       167,071,397         U_32764            SPERT
rs316853653       167,129,709          U_5263            SPERT
rs14914951        167,157,210          D_1359            SPERT
rs13969850        167,166,474         D_11512            SPERT
rs314631516       167,176,627         D_27951            SPERT
rs14914978        167,193,066         U_16438            HTR2A
rs31 2726815      167,689,410         U_16831            HTR2A
rs317496170       167,689,803         D_44012            WDFY2
rs13972990        169,562,091         D_95255            WDFY21
rs318027552       169,613,334      D_1533/D_37584     ALG11/CKAP2
rs13973177        169,663,934        intron/U-3       GTF2F2/KCTD4

SNP               EA/AA       MAF                Beta (SE)

                                             HW               LW

rs314001986        A/G        0.49     -0.277(0.051)    -0.274(0.054)
rs15493356         A/T       0.472     -0.256(0.051)    -0.287(0.053)
rs31 5221322       G/A       0.479     -0.291(0.051)    -0.284(0.054)
rs31 51 20631      A/G       0.479     -0.278(0.052)    -0.289(0.054)
rs15493542         G/A       0.436     -0.267(0.049)     -0.255(0.05)
rs314998864        G/A       0.471      -0.26(0.051)    -0.284(0.053)
rs312642102        T/C       0.478     -0.285(0.052)    -0.287(0.054)
rs316853653        G/C       0.365     -0.287(0.049)    -0.265(0.051)
rs14914951         T/C       0.491     -0.291(0.052)    -0.281(0.054)
rs13969850         C/T       0.474     -0.268(0.051)    -0.285(0.054)
rs314631516        G/A       0.442     -0.307(0.049)    -0.278(0.052)
rs14914978         A/T       0.427     -0.303(0.048)    -0.271(0.051)
rs31 2726815       A/C       0.402     -0.314(0.048)     -0.277(0.05)
rs317496170        T/A       0.401     -0.313(0.048)     -0.266(0.05)
rs13972990         T/C       0.456     -0.314(0.048)    -0.276(0.051)
rs318027552        G/A        0.34     -0.276(0.047)    -0.244(0.049)
rs13973177         G/C       0.451     -0.278(0.049)     -0.26(0.051)

SNP                         Beta (SE)               p-value

                       PW               GW             HW

rs314001986      -0.309(0.055)    -0.295(0.056)     8.23E-08
rs15493356       -0.315(0.054)    -0.278(0.054)     5.49E-07
rs31 5221322     -0.334(0.055)    -0.314(0.055)     1.9E-08
rs31 51 20631    -0.335(0.055)    -0.331(0.056)     8.52E-08
rs15493542        -0.31(0.051)    -0.292(0.049)     5.08E-08
rs314998864       -0.32(0.054)    -0.273(0.054)     4.06E-07
rs312642102      -0.323(0.055)    -0.339(0.055)     3.98E-08
rs316853653      -0.356(0.051)    -0.377(0.049)     6.01E-09
rs14914951        -0.33(0.055)    -0.345(0.056)     1.98E-08
rs13969850       -0.336(0.055)    -0.331(0.055)     2.18E-07
rs314631516      -0.363(0.052)    -0.366(0.052)     5.6E-10
rs14914978       -0.336(0.051)     -0.33(0.05)      4.74E-10
rs31 2726815      -0.383(0.05)    -0.399(0.049)     9.46E-11
rs317496170       -0.386(0.05)    -0.403(0.049)     1.13E-10
rs13972990       -0.386(0.051)     -0.44(0.05)      1.27E-10
rs318027552      -0.351(0.049)     -0.38(0.047)     4.82E-09
rs13973177       -0.358(0.051)    -0.399(0.051)     1.59E-08

SNP                            p-value                  CPV (%)

                     LW           PW           GW          HW

rs314001986       4.2E-07      2.41E-08     1.27E-07     3.691
rs15493356        6.2E-08      6.11E-09     2.84E-07     3.171
rs31 5221322      1.4E-07      1.29E-09     1.51E-08      4.04
rs31 51 20631     9.38E-08     1.54E-09     3.11E-09     3.725
rs15493542        3.98E-07     1.12E-09     3.32E-09     3.581
rs314998864       1.02E-07     3.97E-09     5.81E-07     3.269
rs312642102       1.14E-07     5.33E-09     1.24E-09     3.889
rs316853653       1.99E-07     5.06E-12     4.33E-14     3.883
rs14914951        2.38E-07     2.58E-09     7.65E-10     4.044
rs13969850        1.21E-07     1.1E-09      2.12E-09     3.481
rs314631516       7.97E-08     5.66E-12     1.76E-12     4.507
rs14914978        1.09E-07     8.28E-11     5.7E-11      4.403
rs31 2726815      4.56E-08     4.38E-14     5.29E-16     4.723
rs317496170       1.51E-07     2.64E-14     2.47E-16     4.677
rs13972990        7.3E-08      5.49E-14     6.9E-18       4.69
rs318027552       6.86E-07     8.44E-13     1.48E-15     3.341
rs13973177        4.75E-07     4.63E-12     6.07E-15     3.694

SNP                         CPV (%)

                    LW         PW         GW

rs314001986       3.819      5.042      4.964
rs15493356        4.151      5.264      4.417
rs31 5221322      4.066      5.786      5.496
rs31 51 20631     4.219      5.835      6.065
rs15493542        3.421      5.147      4.918
rs314998864       4.064      5.424       4.3
rs312642102       4.149       5.47      6.318
rs316853653       3.512      6.376      7.458
rs14914951        3.992      5.661      6.538
rs13969850        4.128      5.919      6.139
rs314631516       4.006      6.813      7.356
rs14914978        3.772      5.856      6.027
rs31 2726815       3.86      7.363      8.479
rs317496170        3.58      7.454      8.618
rs13972990         3.75      7.312      10.112
rs318027552       2.781      5.869      7.356
rs13973177        3.351      6.536      8.581

SNP, single nucleotide polymorphism; GGA, Gallus gallus autosome; EA,
effect allele (minor allele); AA, alternative allele (major allele);
MAF, minor allele frequency; Beta, estimated allelic substitution
effect per copy of the effect allele; SE, standard error of the beta;
CPV, contribution to phenotypic variance (%); HW, heart weight; LW,
liver weight; PW, proventriculus weight; GW, gizzard weight. GTF2F2,
general transcription factor IIF polypeptide 2; SPERT, spermatid
associated; HTR2A, 5-hydroxytryptamine receptor 2A; WDFY2, WD repeat
and FYVE domain containing 2; ALG11, asparagine-linked glycosylation
11; CKAP2, cytoskeleton associated protein 2; GTF2F2, general
transcription factor IIF polypeptide 2; KCTD4, potassium channel
tetramerization domain containing 4.

Table 5. Information of 6 SNPs on GGA4 and genomic regions to internal
organ weight

    SNP         Position       Location (bp)      Candidate/    EA/AA
                  (bp)                           nearest gene

 rs15619270    75,478,306     D_26163/D_1994      LCORUNCAPG     G/C
rs315201454    75,485,620     intron/D_33477     NCAPG/LCORL     A/G
 rs14491030    75,486,534    missense/D_34391    NCAPG/LCORL     A/G
rs316243629    75,495,451         intron            NCAPG        G/A
 rs14491074    75,653,956        U_102910            LDB2        T/C
 rs14491507    76,073,771         U_18141        PROM1/TAPT1     G/A

    SNP         MAF                     Beta (SE)

                             HW             LW             PW

 rs15619270    0.058    0.571(0.084)   0.437(0.085)   0.466(0.085)
rs315201454    0.058    0.581(0.084)   0.457(0.086)   0.465(0.086)
 rs14491030    0.059    0.549(0.083)   0.427(0.084)   0.441(0.085)
rs316243629    0.058    0.581(0.084)   0.457(0.086)   0.465(0.086)
 rs14491074    0.057    0.558(0.086)   0.479(0.087)   0.493(0.087)
 rs14491507    0.051    0.623(0.09)    0.468(0.092)   0.53(0.092)

    SNP         Beta (SE)                   p-value

                    GW            HW           LW           PW

 rs15619270    0.531(0.078)    1.58E-11     3.38E-07     5.5E-08
rs315201454    0.535(0.079)    7.75E-12     1.03E-07     6.7E-08
 rs14491030    0.527(0.077)    6.32E-11     4.86E-07     2.08E-07
rs316243629    0.535(0.079)    7.75E-12     1.03E-07     6.7E-08
 rs14491074    0.574(0.081)    1.16E-10     4.05E-08     1.79E-08
 rs14491507    0.61(0.085)     7.99E-12     3.94E-07     1.06E-08

    SNP         p-value                    CPV (%)

                   GW         HW       LW       PW       GW

 rs15619270     1.52E-11    3.503     1.98    2.369     3.29
rs315201454     1.52E-11    3.679    2.189    2.422    3.365
 rs14491030     1.4E-11     3.284    1.892    2.179     3.28
rs316243629     1.52E-11    3.679    2.189    2.422    3.365
 rs14491074     1.66E-12    3.333    2.414    2.636    3.821
 rs14491507     1.31E-12    3.787    2.075    2.794    3.981

SNP, single nucleotide polymorphism; GGA, Gallus gallus autosome; EA,
effect allele (minor allele); AA, alternative allele (major allele);
MAF, minor allele frequency; Beta, estimated allelic substitution
effect per copy of the effect allele; SE, standard error of the beta;
CPV, contribution to phenotypic variance (%); HW, heart weight; LW,
liver weight; PW, proventriculus weight; GW, gizzard weight. LCORL,
ligand dependent nuclear receptor corepressor like; NCAPG, non-SMC
condensin I complex subunit G; LDB2, LIM domain binding 2; PROMI,
prominin 1; TAPTI, transmembrane anterior posterior transformation 1.

Table 6. Information of 6 SNPs on GGA2 and genomic regions to heart

SNP               Position      Location       Candidate/      EA/AA
                     (bp)          (bp)       nearest gene

rs15060276        7,788,887      D_74,602        HTR5A          T/C
rs316413705       8,083,684       intron          SHH           T/C
rs313026677       8,096,208      U_6,355          SHH           A/T
rs14135507        8,159,096      U_69,243         SHH           G/A
rs315048247       8,175,293      U_85,440         SHH           C/A
rs316318736       8,393,250     U_303,397         SHH           A/G

SNP                MAF        Beta (SE)        p-value      CPV (%)

rs15060276        0.400      0.235(0.045)      2.50E-07      2.489
rs316413705       0.281      0.282(0.049)      1.35E-08      3.109
rs313026677       0.383      0.235(0.047)      8.08E-07      2.293
rs14135507        0.275      0.246(0.048)      4.89E-07      3.027
rs315048247       0.300      0.271(0.047)      1.49E-08      2.213
rs316318736       0.302      0.229(0.046)      7.57E-07      2.423

SNP, single nucleotide polymorphism; GGA, Gallus gallus autosome; EA,
effect allele (minor allele); AA, alternative allele (major allele);
MAF, minor allele frequency; Beta, estimated allelic substitution
effect per copy of the effect allele (EA); SE, standard error of the
beta, which means the effect size of minor alleles; CPV, contribution
to phenotypic variance (%); HTR5A, 5-Hydroxytryptamine receptor 5A;
SHH, sonic hedgehog.
COPYRIGHT 2019 Asian - Australasian Association of Animal Production Societies
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2019 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Dou, Taocun; Shen, Manman; Ma, Meng; Qu, Liang; Li, Yongfeng; Hu, Yuping; Lu, Jian; Guo, Jun; Wang,
Publication:Asian - Australasian Journal of Animal Sciences
Article Type:Report
Date:Mar 1, 2019
Previous Article:Estimation of the genetic milk yield parameters of Holstein cattle under heat stress in South Korea.
Next Article:Effects of exercise on myokine gene expression in horse skeletal muscles.

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |