Skip to main content

Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens

Abstract

Background

Wenchang chickens are one of the most popular local chicken breeds in the Chinese chicken industry. However, the low feed efficiency is the main shortcoming of this breed. Therefore, there is a need to find a more precise breeding method to improve the feed efficiency of Wenchang chickens. In this study, we explored important candidate genes and variants for feed efficiency and growth traits through genome-wide association study (GWAS) analysis.

Results

Estimates of genomic heritability for growth and feed efficiency traits, including residual feed intake (RFI) of 0.05, average daily food intake (ADFI) of 0.21, average daily weight gain (ADG) of 0.24, body weight (BW) at 87, 95, 104, 113 days of age (BW87, BW95, BW104 and BW113) ranged from 0.30 to 0.44. Important candidate genes related to feed efficiency and growth traits were identified, such as PLCE1, LAP3, MED28, QDPR, LDB2 and SEL1L3 genes.

Conclusion

The results identified important candidate genes for feed efficiency and growth traits in Wenchang chickens and provide a theoretical basis for the development of new molecular breeding technology.

Peer Review reports

Introduction

Poultry production is an important entreprise worldwide. Chicken is considered a healthy white meat source that is lower in fat, calories and cholesterol than other red meat sources [1]. In recent years, the promotion of white meat consumption has gradually become a trend [2]. As an important part of the white meat market, chickens represent an efficient and inexpensive source of animal protein [3]. Since the 1980s, with the continuous development of modern broiler production, China has become the second largest country in the world in terms of chicken production and consumption [4]. The main breeds of chickens are native breeds and broilers [5]. Wenchang chicken (Fig. 1) is one of the most popular chicken breeds in the local chicken industry, originating from Hainan Island in the South China Sea and has been raised for 400 years. It is famous for its excellent meat quality and is one of the four famous dishes from Hainan [6]. It’s an economic mainstay of animal husbandry in Hainan Province, with an annual output of nearly 100 million chickens and a total output value of 1.78 billion dollars in 2020 [7].

Fig. 1
figure 1

Picture for Wenchang chicken. (a), roosters (b), hens

However, the primary shortcoming of Wenchang chickens is their low feed efficiency [8]. Therefore, there is an urgent need to find a more reliable breeding method to improve the feed efficiency of Wenchang chickens. Feed represents more than 70% of the total cost of poultry production, and improving feed efficiency has consistently been the goal of any chicken breeding strategy [9]. Feed efficiency is contingent upon the relation between the feed intake (FI) and the growth (or bodyweight gain) of an animal and is quantified by several indexes, such as RFI and feed conversion rate (FCR) [10,11,12]. Feed efficiency is influenced by several factors, including the breed and its sex, age, diet, and management [13, 14].The RFI is an important index measuring feed efficiency, and it is defined as the difference between the actual and expected feed intake [15].In 1963, the concept of RFI was first proposed in beef cattle research [16], and it was first applied to chickens by Luiting in 1991 [17], covering the calculation methods of RFI, heritability calculation, and the phenotypic and genetic correlations with relevant traits. Since RFI reflects the variation in feed efficiency, it appears to be independent of growth traits [15]. Studies have shown that RFI is moderately heritable in poultry [12, 18]. The study by Bai et al. demonstrates that selecting for low RFI can improve poultry feed efficiency without compromising growth performance [19]. Growth traits are key in poultry breeding, and the properties of growth traits need to be considered in breeding for feed efficiency. In recent years, GWAS has been applied in poultry to explore the association between host genetics and economic traits [20]. Earlier studies identified candidate genes including NSUN3, EPHA6, and AGK for broilers RFI, while LAP3 was identified as a candidate gene for broiler body weight [21,22,23,24,25,26]. By deeply studying the functions and regulatory mechanisms of these genes, our aim is to provide more efficient and sustainable solutions for the chicken breeding technologies [4].

In the breeding process of this study, the focus will be on maintaining the excellent meat quality while moderately improving the feed efficiency and growth rate. The objectives of this study were (1) to determine the inheritance pattern of feed efficiency and growth traits to provide a basis for formulating a better breeding method for Wenchang chickens, and (2) to find key variants affecting feed efficiency and growth traits of Wenchang chickens and also to elucidate the molecular mechanism of feed efficiency and growth traits.

Materials and methods

Experimental birds

The study uses 1,547 chickens hatched and raised by Hainan Tanniu Wenchang Chicken Co., Ltd. These birds are from a commercial line of Wenchang chicken selected for 18 generations. All birds were raised in three-tier battery cages (one bird per cage) under the same management and nutritional conditions. The diet was formulated based on the National Research Council (NRC) requirements [27] and the Feeding Standards of Chickens established by the Ministry of Agriculture, Beijing, China [28]. The chickens were raised from 87 to 113 days to collect data on individual phenotypes.

Phenotypes

The phenotypes measured in this study included body weight and feed intake of chickens.

For body weight measurements, an electronic scale with a precision of 0.1 g was utilized. The trough was emptied of any remaining feed 12 h prior to weighing, and body weights were recorded at 87, 95, 104, and 113 days of age.

Feed intake in this experiment was recorded for 26 days. Fresh feed was added at a fixed time every day and recorded. In addition to the above measured phenotypes, the phenotypes calculated from the measured phenotypes included FI, ADFI, ADG, Metabolic weight in the middle of the test (MWT), and RFI.

FI, g.

$$\eqalign{& FI\hspace{0.17em}=\hspace{0.17em}feed\,amount\,during the\,test\hspace{0.17em}- \cr & \hspace{0.17em}Weight\,of\,feed\,remaining\,at\,the\,end\,of \,the\,test\cr}$$
(1)

ADFI, g/d.

$$\text{A}\text{D}\text{F}\text{I}=\text{F}\text{I} / 26$$
(2)

ADG, g/d.

$$\text{A}\text{D}\text{G}= (\text{B}\text{W}113-\text{B}\text{W}87) / 26$$
(3)

MWT, g.

$$\text{M}\text{W}\text{T}= \left[\right(\text{B}\text{W}87 + \text{B}\text{W}113) / {2]}^{0.75}$$
(4)

RFI, g/d: The RFI was obtained from the multiple regression equation of average daily feed intake and metabolic weight in the middle of the test and average daily weight gain [18]. The equation was as follows:

$$\text{R}\text{F}\text{I}=\text{A}\text{D}\text{F}\text{I}-({\upmu } + \text{s}\text{e}\text{x} + {\upbeta }1 \text{M}\text{W}\text{T} + {\upbeta }2 \text{A}\text{D}\text{G})$$
(5)

where ADFI represents the mean average daily feed intake, µ is the intercept, sex is a fixed effect, MWT and ADG are as defined above, β1 and β2 represent the partial regression coefficient.

Genotyping, imputation, and quality control (QC)

At 112 days of age, 2 mL of blood was collected from the wing vein, placed in an anticoagulant tube containing EDTA-K2, mixed, and stored at -20℃ until later analysis. Genomic DNA was extracted from blood samples with the phenol-chloroform method. Genotyping was conducted with a customized chicken 55 K SNP array (Beijing Compass Biotechnology Co., Ltd., Beijing, China) [29]. QC of the generated genotype data was achieved using PLINK (V2.0) (https://www.cog-genomics.org/plink/2.0/) software. The specific process setting the individual genotype detection rate at ≥ 90%, single SNP detection rate at ≥ 90%, minimum allele frequency (MAF) of 95%, and retention of SNPs on autosomes 1–28. A total of 45,278 SNPs in 1,479 chickens (762 males and 717 females) passed the QC. Whole genome resequencing of 247 individuals from the 1,479 Wenchang chickens, including 25 males and 222 females was carried out. The sequencing generated 150 bp paired end reads on an Illumina NovaSeq 6000 platform with the average depth of approximately 10×, at the Shenzhen BGI Co., Ltd. The QC standards were as follows: setting the detection rate of individual genotype at ≥ 90%, the detection rate of single SNP site at ≥ 90%, a MAF of 95%, and Hardy-Weinberg equilibrium (HWE) at P < 0.000001. A total of 12,590,784 autosomal SNPs in 247 individuals were passed the QC.

Beagle 5.2 software was used to impute the 55 K chip data to the whole-genome sequence (WGS) level [30]. Before imputation, inconsistencies between the target panel and the reference panel were checked using conform-gt software (http://faculty.washington.edu/browning/conform-gt.html). Then, the 55 K SNP chip data were populated to the resequencing level using Beagle 5.2 software. The QC condition: HWE at P < 1.00e-6, setting the detection rate of individual genotype at ≥ 90%, the detection rate of single SNP site at ≥ 90%, and a MAF of ≥ 0.05. A total of 12,184,765 autosomal SNPs from 1,479 samples passed QC.

Estimation of heritability and genetic correlations

Phenotypic correlation coefficients were calculated using ggpairs within the GGally R package, and then genetic parameter estimation was performed using restricted maximum likelihood (REML) of GCTA v1.93.2 beta software [31, 32]. The statistical model used was:

$$y = Xb+Z\alpha +e$$
$$\eqalign{y &= Xb + Za + e \cr Var\left[ {\matrix{\alpha \cr e \cr } } \right] &= \left[ {\matrix{{G\sigma _a^2} & 0 \cr 0 & {1\sigma _e^2} \cr } } \right] \cr}$$
(6)

where, \(\text{y}\) is a vector of observations, \(\text{b}\) is a vector of fixed effects (i.e., sex), \({\upalpha }\) is the random vector representing the genomic effects,\(e\)is the vector of random residual effects, \(\text{X}\)and \(\text{Z}\) are incidence matrices. The distribution of the random animal effect \(\alpha\) is \(\alpha \ \sim N(0,\,{\rm{G}}\sigma _a^2)\) with \(\text{G}\) is being the genomic relationship matrix, and \({\sigma }_{a}^{2}\) being the additive genetic variance.

Genome-wide association study

A GWAS analysis was carried out between all the genotyped SNPs and feed efficiency and growth traits using a mixed linear model (MLM). The MLM for feed efficiency and growth traits was performed using sex (female or male) as a fixed effect and the top three principal components (PCs) as covariates. All association tests were performed using the MLM option in GCTA based on the following model [33]:

$$y = Xb + Z\mu + e$$
(7)

where y is a vector of observations, X and Z are incidence matrices for the vectors for parameters b and µ, b is a vector of fixed effects including the sex and three eigenvectors from principal component analysis (PCA), µ is the vector of the additive genetic effect of the candidate SNP to be tested for association, and e is the vector of the residual effect.

The Bonferroni correction method was used in this study to determine the significance thresholds, and the formula for performing the Bonferroni corrected multiple tests was as follows:

$$P=\alpha /N$$
(8)

Where \(P\) is the corrected significance threshold,\(\alpha\)represents the significance threshold for a single test, and N represents the number of multiple hypothesis tests, i.e., the number of SNPs analysed by GWAS. We calculated the number of genome-wide independent markers using the PLINK (V1.9) command -indep-pairwise, with a window size of 25 SNPs, a step of five SNPs, and an r2 threshold of 0.2. Manhattan and quantile‒quantile (Q-Q) plots were derived from the GWAS results using the qqman (https://cran.r-project.org/web/packages/qqman/) and Cairo (http://www.rforge.net/Cairo/) packages within R software (http://www.r-project.org/). LD blocks of target regions were performed using Haploview v4.2 software [34]. For additive and dominance effects of important SNPs on traits, the calculation process in this study was done in ASReml v4.1 software (https://asreml.kb.vsni.co.uk/knowledge-base/asreml_documentation). The SNP positions were updated according to the newest release from Ensembl (https://asia.ensembl.org/index.html). Identification of the closest genes to genome-wide significant and suggestive variants was obtained using the ChIPpeakAnno package (https://www.bioconductor.org/packages/devel/bioc/vignettes/ChIPpeakAnno/inst/doc/pipeline.html). Gene function enrichment analysis was performed using bioinformatics (https://www.bioinformatics.com.cn).

Differential expression analysis

In this study, the RFI or BW of 127 broiler chickens was ranked from low to high. Using individuals with the highest (n = 15) and lowest RFI (n = 15) phenotypes, and individuals with the highest (n = 15) and lowest (n = 15) BW. The Wenchang chickens used in this study were a mix of males and females, so we did not differentiate between genders in the 30 chickens selected. Gene expression of important candidate genes in the different groups can demonstrate the reliability of the GWAS results. Figures were generated using GraphPad Prism 8 [4].

Results

Description of phenotypic traits

Descriptive statistics of the feed efficiency traits of Wenchang chickens are shown in Table 1; Fig. 2. According to the body weights at 87, 95, 104 and 113 days, Wenchang chickens exhibited slow growth, with an ADFI of 82.52 g/d and an ADG of 8.66 g/d. The RFI ranged from − 39.14 to 36.25 g/d.

Table 1 Descriptive statistics for feed efficiency and growth traits of Wenchang chickens (N1 = 1532)
Fig. 2
figure 2

Phenotypic data and correlation analysis of Wenchang chicken

Genetic parameters’ estimates

Comparison of density markers before and after the imputation in Fig. 3 revealed a significant increase in loci after imputation. The estimated genetic parameters of residual feed intake and body weight are shown in Table 2. The heritability of the RFI, BW, ADFI and ADG ranged from moderate and low level (0.05–0.44). RFI had the highest genetic association with ADFI and ADG, at 0.92 and 0.86, respectively (P < 0.001). Regarding phenotypic correlation, the genetic correlation between RFI and ADFI was 0.63 and decreased with body weight to a minimum of 0.20. Based on the result, the heritability of RFI in Wenchang chicken is low, while the heritability of body weight is not significantly different from that of other breeds.

Fig. 3
figure 3

Comparison of density markers before and after the imputation

Table 2 Estimates of genetic parameters for feed efficiency and growth traits of Wenchang chickens

Genome-wide association study of feed efficiency and growth traits

The GWAS results showed significant SNPs on GGA 2, 6 and 26 were associated with RFI (expansion coefficient λ of 0.987) (Fig. 4a-b). One of significant SNPs 6_21123592 on GGA6 was located on the intron of the PLCE1 gene, with this SNP explaining 2.46% of the genetic variation. Additionally, significant SNPs 2_45795056 and 26_2851843 were located in the lncRNA introns of the ENSGALG00000052614 and ENSGALG00000001264 genes, respectively (Table 3). The GWAS results showed significant SNPs on GGA 2, 4, 6, 19, and 28 were associated with ADFI (expansion coefficient λ of 0.994) (Fig. 4c-d). Specifically, the significant SNPs 4_75971941, 4_85488222, 19_2763444, and 28_388715 were found on the genes LAP3, IMMT, GATSL2, and FBN3, respectively (Table 3). The SNP 6_21123592 for RFI in the CT genotype exhibited significantly higher values than those in the CC genotype (Table 4). Another SNP 2_45795056 showed a significant additive effect on ADFI but a nonsignificant effect on the two feed efficiency traits assessed in this study. The estimated additive effect of candidate SNP 6_21123592 on RFI and ADFI were − 0.81 and 2.67 respectively, while the dominance effect was − 2.99 and 0.93 respectively. Significant dominance effects were observed on both RFI and ADFI, while the additive effect was not significant (Table 5).

Fig. 4
figure 4

Manhattan and Q-Q plots of GWAS for feed efficiency traits. The horizontal red lines indicate the thresholds for genome-wide significance (P = 3.47E-08), and the horizontal blue lines indicate the thresholds for suggestive significance (P = 6.93E-07)

Table 3 The significant SNPs information of the feed efficiency traits in Wenchang chicken
Table 4 Effect of candidate SNPs on feed efficiency traits
Table 5 Additive and dominance effects of candidate SNPs on feed efficiency traits1

The GWAS analysis results for growth traits at different ages are presented in Fig. 5; Table 6. Significant QTLs affecting BW, ADG, and ADFI were observed in the 73.15–76.55 Mb interval of GGA4. Specifically, BW87 (Fig. 5a, Table S1) had 661 SNPs in the significant interval, BW95 had 933 SNPs (Fig. 5c, Table S2), BW104 had 1018 SNPs (Fig. 5e, Table S3), BW113 had 1150 SNPs (Fig. 5g, Table S4), and ADFI and ADG had 3 and 7 SNPs (Table 3; Fig. 4c and e). The most prominent SNP among these traits was 4_75971941, located in the intron of the LAP3 gene. There were 13 important candidate genes in the colocalization interval of different day-ages, including LAP3, KCNIP4, NCAPG, LDB2, FAM184B, SEL1L3, ZCCHC4, PPARGC1A, PACRGL, SLIT2, LCORL, MED28 and QDPR (Table 6). Through gene function enrichment analysis (Figure S1), we identified significant enrichment primarily in functions such as RNA biosynthetic process, regulation of RNA metabolic process, and regulation of cellular macromolecule biosynthetic process.

Fig. 5
figure 5

Manhattan and Q-Q plots of GWAS for growth traits. The horizontal red lines indicate the thresholds for genome-wide significance (P = 3.47E-08), and the horizontal blue lines indicate the thresholds for suggestive significance (P = 6.93E-07)

Table 6 The significant SNP information of growth traits in Wenchang chicken

The four LD blocks were detected within the common location of growth traits on GGA4:75971.28-75974.42 kb (Fig. 6), each containing 3–10 SNPs. The most significant SNP 4_75971941 was located on block1, which is located in the intron of the LAP3 gene. Genotype analysis showed that the body weight of the CC genotype at different ages was significantly higher than that of the CT genotype. Specifically, the weight of the CC genotype for BW87 was 54.89 g heavier than that of the CT genotype, and the weight of the CC genotype for BW113 was 82.80 g heavier than that of the CT genotype (Fig. 6b-g). The estimated additive effects of candidate SNP 4_75971941 on body weight at different ages were 37.97, 46.67, 53.38 and 60.53, while the dominance effects were − 16.92, -16.19, -17.98 and − 22.26. SNP 4_75971941 had significant additive effects on BW and ADG ranging from 37.97 to 60.53. The additive effects increased with age, while the locus had significant dominance effects on BW and nonsignificant dominance effects on ADG and ADFI (Table 7).

Fig. 6
figure 6

Association results of the candidate SNPs on GGA4. (a), LD analysis of the 25 significant SNPs on GGA4. (b, c, d, e, f, g), The phenotypic differences of individuals with different genotypes at rs80610898 on GGA4. *, P < 0.05; **, P < 0.01; ***, P < 0.001; ns, no significance

Table 7 Additive and dominance effects of candidate SNPs on ADFI and growth traits

The candidate sites for body weight were determined to be located in distinct LD blocks. In this study, haplotype association analysis was conducted using the glm model with sex as a fixed effect (Table 8). The analysis indicated that block1 had a highly significant effect (P < 0.01) on BW at various ages, while block2 exhibited a significant effect (P < 0.05) on BW across different age groups. The results revealed that block2, block3, and block4 all had a significant impact on ADFI (P < 0.05), while none of the four LD blocks showed a significant impact on RFI (P > 0.05).

Table 8 The effect of different LD blocks on feed efficiency and growth traits in Wenchang chickens

Differential expression analysis

We conducted gene differential expression analysis on all genes mapped to the feed efficiency and growth traits, respectively. Significant differences in phenotypes were observed between the high and low RFI groups (Fig. 7a) and between the high and low BW groups (Fig. 7c). Between the high and low RFI groups, the PLCE1 and IMMT genes showed significantly (P < 0.01) higher expression levels in the high RFI group than in the low RFI group (Fig. 7b). Similarly, expression level of LAP3, MED28, and QDPR genes was significantly higher in the high BW group than that in the low BW group (P < 0.01), while expression of LDB2 and SEL1L3 genes was significantly (P < 0.05) higher in the low BW group than the high BW group (Fig. 7d).

Fig. 7
figure 7

The expressions of candidate genes for feed efficiency and growth traits. (a), Analysis of RFI Differences between high and low RFI groups. (b), Differential expression of related candidate genes between high and low RFI groups. (c), Analysis of BW Differences between high and low BW groups. (d), Differential expression of related candidate genes between high and low weight groups. Data were expressed as the mean ± SEM (n = 15), *** P < 0.001, ** P < 0.01, * P < 0.05

Discussion

Phenotypic and estimation of the genetic parameters of Wenchang chickens

Body weight traits (BW87-BW113) of Wenchang chickens were selected for the GWAS analysis because Wenchang chickens are listing age within that age range. Through different trait measurements at this stage, our study revealed that the BW87 was 1,409.48 g, that of BW113 was 1,894.36 g, the ADG was 18.66 g/d, and the relative growth rate was 25.60%. These findings are consistent with previous studies on yellow-feathered chickens in China [35, 36]. The low heritability of RFI may be caused by the lack of systematic breeding efforts and the chickens are in the latter stage of growth, similar to the results of previous studies [37, 38]. Heritability of BW, ADFI and ADG revealed in the current study ranged from medium to low (0.21–0.44), which was similar to that of other breeds [39,40,41]. The phenotypic and genetic correlations between RFI and ADFI were 0.92 and 0.63, similar to the research of Shirali et al [42]. To optimize breeding outcomes, it is recommended to integrate additional metrics, such as ADFI, and implement a multi-trait selection approach throughout the breeding process, ultimately enhancing breeding results in Wenchang chickens.

Genome-wide association study of feed efficiency and growth traits

To investigate the genetic architecture of feed efficiency and growth traits in Wenchang chickens, a mixed linear model was used for GWAS analysis of related traits. Through GWAS analysis, we identified significant SNP 6_21123592 on GGA6 was located within the intron of the PLCE1 gene associated with RFI. Previous studies have shown that PLCE1 gene is highly expressed in the nervous system and belongs to the phosphoinositide-specific phospholipase C family. The production of second messenger molecules such as diacylglycerol is regulated by activated phosphatidylinositol-specific phospholipase C enzymes, which mediate small GTPases of the Ras superfamily through the activity of its Ras guanine exchange factor. As the effector of heterotrimer and small G protein, PLCE1 is involved in regulating cell growth, T-cell activation, actin organization and cell survival. Mapping the PLCE1 gene function were mostly related to nervous system activity, which regulates the function of the brain in different ways, and the brain was key to regulating diet behavior and body energy homeostasis [43,44,45,46].

The significant QTL affecting BW, ADG, and ADFI was located in the 73.15–76.55 Mb interval of GGA4.There were 13 important candidate genes in the colocalization interval related to BW, including LAP3, KCNIP4, NCAPG, LDB2, FAM184B, SEL1L3, ZCCHC4, PPARGC1A, PACRGL, SLIT2, LCORL, MED28 and QDPR. Among these candidates LAP3, MED28, QDPR, LDB2 and SEL1L3 demonstrated differential expression between high and low groups. Comparison with the Animal QTL database (Chicken QTL Database at (animalgenome.org)) reveals a total of 305 QTLs related to BW and 221 QTLs associated with ADG within 73.15–76.55 Mb interval of GGA4. LAP3 has been shown to catalyze the hydrolysis of the amino-terminal leucine residues of protein or peptide substrates, with diverse functions in mammals, invertebrates, microbes, and plants [47]. The primary function of LAP3 lies in protein maturation and degradation, processes crucial for metabolism, development, adaptation, and repair [48]. LAP3 gene variation may underlie variations in growth rates among species and significant genetic polymorphism of traits of interest in breeding, potentially leading to applications in animal breeding. Furthermore, studies have been conducted on its SNP and its association with growth traits in mammals, such as bovine [49]. Another study found that the LAP3 gene may have a potential function affecting muscle development in sheep [50]. Prenatal development stages are directly related to the growth and development of individual skeletal muscle, which determines the number of muscle fibers and postnatal muscle mass and further exerts long-term effects on the postnatal growth of animals [51,52,53]. Related research that tracked LAP3 mutations in Hu sheep populations reportedly linked, the mutations with body weight at different growth stages [54]. In poultry LAP3 gene was found to be associated with chicken growth traits [22].

Related study linked the MED28 gene with live weight in sheep [55].We found that the MED28 gene was related to muscle development in pig [56]. QDPR for an enzyme that regulates tetrahydrobiopterin (BH4), a cofactor for enzymes involved in neurotransmitter synthesis and blood pressure regulation. Therefore, QDPR gene are also important genes in the regulation of growth [57]. Previous studies have shown that the LDB2 gene located at GGA4 is important for chicken growth traits, and a 31-bp indel was significantly correlated with multiple growth and carcass traits in the F2 population and affected the expression of the LDB2 gene in muscle tissue [26, 58, 59]. It was also identified as an important candidate gene for rapid growth in chickens and had the strongest association with late body weight in Jinghai yellow chicken hens [22, 60, 61]. According to relevant studies, the KCNIP4 gene is located on GGA4, and the candidate gene belongs to the potassium channel interaction protein family and has a wide range of physiological functions, including heart rate regulation, insulin secretion, neurotransmitter release, and smooth muscle contraction. It was considered to be an important candidate gene for growth traits of chickens. In addition, it was reported in different breeds and different growth stages, which also verified that the KCNIP4 and FAM184B genes can affect the growth and development of chickens [58, 62, 63]. Studies have shown that the NCAPG gene is an important candidate gene for mammalian body size growth traits in growth trait association analysis of horses, sheep, and domestic donkeys and is involved in chromosome condensation and methylation [64,65,66,67]. The FAM184B gene has been found in previous studies to be associated with cattle carcass weight [68].

Differential expression analysis

Transcriptomic data enable the quantification of DNA or RNA abundance and expression levels [69]. Differential analysis between different groups is conducted by measuring the expression levels of gene RNA, thereby identifying distinct patterns and variations in gene expression among different groups. The results of related studies also confirm gene expression data of important candidate genes in different groups can demonstrate the reliability of the GWAS results [70, 71]. We showed that expression of PLCE1 in the high and low RFI groups of broilers, validatied it a key candidate gene for RFI. We found that the candidate gene LAP3 was related to the BW, ADFI and ADG traits, as demonstrated by the significant difference in the expression of LAP3 in high and low body recombination in broilers, and relevant research reports also indicated that this might be an important candidate gene affecting growth traits. This study also found that the candidate genes LDB2, KCNIP4, FAM184B, and NCAPG were located on GGA4 and were related to growth traits, which provides an important reference value for subsequent research on the growth traits of Wenchang chickens.

The validity of the SNPs and candidate genes obtained in this study is worth further extensive verification. The significant SNPs and candidate genes identified in this study can be incorporated into the chip markers in future research. By utilizing genomic selection breeding techniques, these markers can be used to breed and improve the target traits of WenChang chickens, ultimately enhancing breeding efficiency.

Conclusions

In this study, we identified that the significant SNP 6_21123592 was located in candidate gene PLCE1 for feed efficiency traits of Wenchang chickens, and the significant SNP 4_75971941 was located in candidate gene LAP3. Other candidate genes including MED28, QDPR, LDB2, and SEL1L3 were identified for growth traits of Wenchang chickens. This provides a good theoretical basis for developing methods of Wenchang chickens breeding, and by further studying the functions and regulatory mechanisms of these genes, we could provide more efficient solutions for breeding of these chickens.

Data availability

The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2021) in National Genomics Data Center (Nucleic Acids Res 2022), China National Center for Bioinformation / Beijing Institute of Genomics, Chinese Academy of Sciences that are publicly accessible at https://bigd.big.ac.cn/gsa/browse/CRA016976.

References

  1. Yang XT, Sun JH, Zhao GP, Li W, Tan XD, Zheng MQ, Feng FR, Liu DW, Wen J, Liu RR. Identification of major loci and candidate genes for meat production-related traits in broilers. Front Genet 2021, 12.

  2. Kong F, Zhao G, He Z, Sun J, Wang X, Liu D, Zhu D, Liu R, Wen J. Serum creatine kinase as a Biomarker to predict wooden breast. Front Physiol. 2021;12:711711.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Liu Y, Li HJ, Wang M, Zhang XH, Yang L, Zhao CJ, Wu CX. Genetic architectures and selection signatures of body height in Chinese indigenous donkeys revealed by next-generation sequencing. Anim Genet. 2022;53(4):487–97.

    Article  CAS  PubMed  Google Scholar 

  4. Tan X, Liu R, Zhao D, He Z, Li W, Zheng M, Li Q, Wang Q, Liu D, Feng F et al. Large-scale genomic and transcriptomic analyses elucidate the genetic basis of high meat yield in chickens. J Adv Res 2023.

  5. Havenstein GB. Poultry breeding and Genetics. Poult Sci. 1991;70(3):662–3.

    Article  Google Scholar 

  6. Tan Z, Luo LL, Wang XZ, Wen Q, Zhou L, Wu KB. Characterization of the cecal microbiome composition of Wenchang chickens before and after fattening. PLoS ONE 2019, 14(12).

  7. Gu LH, Jiang QC, Chen YY, Zheng XL, Zhou HL, Xu TS. Transcriptome-wide study revealed m6A and miRNA regulation of embryonic breast muscle development in Wenchang chickens. Front Vet Sci 2022, 9.

  8. Tian S, Tang W, Zhong Z, Wang Z, Xie X, Liu H, Chen F, Liu J, Han Y, Qin Y et al. Identification of Runs of Homozygosity Islands and Functional Variants in Wenchang Chicken. Animals (Basel) 2023, 13(10).

  9. Aggrey SE, Karnuah AB, Sebastian B, Anthony NB. Genetic properties of feed efficiency parameters in meat-type chickens. Genet Sel Evol 2010, 42.

  10. Zampiga M, Calini F, Sirri F. Importance of feed efficiency for sustainable intensification of chicken meat production: implications and role for amino acids, feed enzymes and organic trace minerals. World Poult Sci J. 2021;77(3):639–59.

    Article  Google Scholar 

  11. Titus HW, Mehring AL, Brumbaugh JH. Variation of feed Conversion. Poult Sci. 1953;32(6):1074–7.

    Article  Google Scholar 

  12. Li W, Zheng MQ, Zhao GP, Wang J, Liu J, Wang SL, Feng FR, Liu DW, Zhu D, Li QH et al. Identification of QTL regions and candidate genes for growth and feed efficiency in broilers. Genet Sel Evol 2021, 53(1).

  13. Sell-Kubiak E, Wimmers K, Reyer H, Szwaczkowski T. Genetic aspects of feed efficiency and reduction of environmental footprint in broilers: a review. J Appl Genet. 2017;58(4):487–98.

    Article  PubMed  PubMed Central  Google Scholar 

  14. He Z, Liu R, Wang M, Wang Q, Zheng J, Ding J, Wen J, Fahey AG, Zhao G. Combined effect of microbially derived cecal SCFA and host genetics on feed efficiency in broiler chickens. Microbiome. 2023;11(1):198.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Wen CL, Yan W, Zheng JX, Ji CL, Zhang DX, Sun CJ, Yang N. Feed efficiency measures and their relationships with production and meat quality traits in slower growing broilers. Poult Sci. 2018;97(7):2356–64.

    Article  CAS  PubMed  Google Scholar 

  16. Koch RM, Swiger LA, Chambers D, Gregory KE. Efficiency of feed use in beef cattle. J Anim Sci. 1963;22(2):486–94.

    Article  Google Scholar 

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

    Article  Google Scholar 

  18. Case LA, Wood BJ, Miller SP. The genetic parameters of feed efficiency and its component traits in the Turkey (Meleagris gallopavo). Genet Sel Evol. 2012;44(1):2.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bai H, Guo Q, Yang B, Dong Z, Li X, Song Q, Jiang Y, Wang Z, Chang G, Chen G. Effects of residual feed intake divergence on growth performance, carcass traits, meat quality, and blood biochemical parameters in small-sized meat ducks. Poult Sci. 2022;101(9):101990.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tan X, He Z, Fahey AG, Zhao G, Liu R, Wen J. Research progress and applications of genome-wide association study in farm animals. Anim Res One Health. 2023;1(1):56–77.

    Article  Google Scholar 

  21. Hu ZL, Park CA, Wu XL, Reecy JM. Animal QTLdb: an improved database tool for livestock animal QTL/association data dissemination in the post-genome era. Nucleic Acids Res. 2013;41(D1):D871–9.

    Article  CAS  PubMed  Google Scholar 

  22. Liu R, Sun Y, Zhao G, Wang F, Wu D, Zheng M, Chen J, Zhang L, Hu Y, Wen J. Genome-wide association study identifies loci and candidate genes for body composition and meat quality traits in Beijing-You chickens. PLoS ONE. 2013;8(4):e61172.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Moreira GCM, Salvian M, Boschiero C, Cesar ASM, Reecy JM, Godoy TF, Ledur MC, Garrick D, Mourão GB, Coutinho LL. Genome-wide association scan for QTL and their positional candidate genes associated with internal organ traits in chickens. BMC Genomics. 2019;20(1):669.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Li J, Akanno EC, Valente TS, Abo-Ismail M, Karisa BK, Wang Z, Plastow GS. Genomic heritability and Genome-Wide Association Studies of Plasma Metabolites in crossbred beef cattle. Front Genet. 2020;11:538600.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Dou D, Shen L, Zhou J, Cao Z, Luan P, Li Y, Xiao F, Guo H, Li H, Zhang H. Genome-wide association studies for growth traits in broilers. BMC Genomic Data. 2022;23(1):1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Li Y-d, Bai X, Liu X, Wang W-j, Li Z-w, Wang N, Xiao F, Gao H-h, Guo H-s, Li H, et al. Integration of genome-wide association study and selection signatures reveals genetic determinants for skeletal muscle production traits in an F2 chicken population. J Integr Agric. 2022;21(7):2065–75.

    Article  CAS  Google Scholar 

  27. National Research Council (NRC). In. The grants Register 2023: the complete guide to Postgraduate Funding Worldwide. London: Palgrave Macmillan UK; 2022. pp. 783–4.

    Google Scholar 

  28. Liu L, Liu XJ, Cui HX, Liu RR, Zhao GP, Wen J. Transcriptional insights into key genes and pathways controlling muscle lipid metabolism in broiler chickens. In: Bmc Genomics vol. 20; 2019.

  29. Liu R, Xing S, Wang J, Zheng M, Cui H, Crooijmans R, Li Q, Zhao G, Wen J. A new chicken 55K SNP genotyping array. BMC Genomics. 2019;20(1):410.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Browning BL, Zhou Y, Browning SR. A one-penny Imputed Genome from Next-Generation reference panels. Am J Hum Genet. 2018;103(3):338–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Kumar SK, Feldman MW, Rehkopf DH, Tuljapurkar S. Limitations of GCTA as a solution to the missing heritability problem (113, pg E61, 2015). P Natl Acad Sci USA. 2016;113(6):E813–813.

    Google Scholar 

  32. Kumar SK, Feldman MW, Rehkopf DH, Tuljapurkar S. Limitations of GCTA as a solution to the missing heritability problem. P Natl Acad Sci USA. 2016;113(1):E61–70.

    CAS  Google Scholar 

  33. Yang JA, Lee SH, Goddard ME, Visscher PM. GCTA: A Tool for Genome-wide Complex Trait Analysis. Am J Hum Genet. 2011;88(1):76–82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

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

    Article  CAS  PubMed  Google Scholar 

  35. Tang H, Gong YZ, Wu CX, Jiang J, Wang Y, Li K. Variation of meat quality traits among five genotypes of chicken. Poult Sci. 2009;88(10):2212–8.

    Article  CAS  PubMed  Google Scholar 

  36. Huang Q, Wen C, Yan W, Sun C, Gu S, Zheng J, Yang N. Comparative analysis of the characteristics of digestive organs in broiler chickens with different feed efficiencies. Poult Sci. 2022;101(12):102184.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Yuan J, Dou T, Ma M, Yi G, Chen S, Qu L, Shen M, Qu L, Wang K, Yang N. Genetic parameters of feed efficiency traits in laying period of chickens. Poult Sci. 2015;94(7):1470–5.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Marchesi JAP, Ono RK, Cantao ME, Ibelli AMG, Peixoto JD, Moreira GCM, Godoy TF, Coutinho LL, Munari DP, Ledur MC. Exploring the genetic architecture of feed efficiency traits in chickens. Sci Rep-Uk 2021, 11(1).

  39. Venturini GC, Stafuzza NB, Cardoso DF, Baldi F, Ledur MC, Peixoto JO, El Faro L, Munari DP. Association between ACTA1 candidate gene and performance, organs and carcass traits in broilers. Poult Sci. 2015;94(12):2863–9.

    Article  CAS  PubMed  Google Scholar 

  40. Mebratie W, Shirali M, Madsen P, Sapp RL, Hawken R, Jensen J. The effect of selection and sex on genetic parameters of body weight at different ages in a commercial broiler chicken population. Livest Sci. 2017;204:78–87.

    Article  Google Scholar 

  41. Mebratie W, Madsen P, Hawken R, Romé H, Marois D, Henshall J, Bovenhuis H, Jensen J. Genetic parameters for body weight and different definitions of residual feed intake in broiler chickens. Genet Sel Evol. 2019;51(1):53.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Shirali M, Varley PF, Jensen J. Bayesian estimation of direct and correlated responses to selection on linear or ratio expressions of feed efficiency in pigs. Genet Sel Evol. 2018;50(1):33.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Rao J, Ashraf S, Tan W, van der Ven AT, Gee HY, Braun DA, Feher K, George SP, Esmaeilniakooshkghazi A, Choi WI, et al. Advillin acts upstream of phospholipase C ϵ1 in steroid-resistant nephrotic syndrome. J Clin Invest. 2017;127(12):4257–69.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Strazza M, Adam K, Smrcka AV, Lerrer S, Mor A. PLCepsilon1 suppresses tumor growth by regulating murine T cell mobilization. Clin Exp Immunol. 2020;200(1):53–60.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Yu S, Choi WI, Choi YJ, Kim HY, Hildebrandt F, Gee HY. PLCE1 regulates the migration, proliferation, and differentiation of podocytes. Exp Mol Med. 2020;52(4):594–603.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Huang LH, Liao CD, Wu HH, Huang PW. PLCE1 is a poor prognostic marker and may promote immune escape from osteosarcoma by the CD70-CD27 signaling pathway. Bosnian J Basic Med. 2022;22(6):992–1004.

    Article  CAS  Google Scholar 

  47. Yao HH, Ren FZ, Bao YB, Dong YH, Lin ZH. Molecular characterization and expression of the LAP3 gene and its Association with Growth traits in the blood clam Tegillarca Granosa. Fishes-Basel 2021, 6(4).

  48. de Seny D, Baiwir D, Bianchi E, Cobraiville G, Deroyer C, Poulet C, Malaise O, Paulissen G, Kaiser MJ, Hauzeur JP et al. New proteins contributing to Immune Cell infiltration and pannus formation of synovial membrane from Arthritis diseases. Int J Mol Sci 2021, 23(1).

  49. Zheng X, Ju ZH, Wang J, Li QL, Huang JM, Zhang AW, Zhong JF, Wang CF. Single nucleotide polymorphisms, haplotypes and combined genotypes of LAP3 gene in bovine and their association with milk production traits. Mol Biol Rep. 2011;38(6):4053–61.

    Article  CAS  PubMed  Google Scholar 

  50. Ge L, Su PW, Wang S, Gu YF, Cao XK, Lv XY, Wang SH, Getachew T, Mwacharo JM, Haile A et al. New Insight into the Role of the Leucine Aminopeptidase 3 (LAP3) in Cell Proliferation and Myogenic Differentiation in Sheep Embryonic Myoblasts. Genes-Basel 2022, 13(8).

  51. Rehfeldt C, Kuhn G. Consequences of birth weight for postnatal performance and carcass quality in pigs as related to myogenesis. J Anim Sci. 2006;84(Supplsuppl):E113–123.

    Article  PubMed  Google Scholar 

  52. Du M, Tong J, Zhao J, Underwood KR, Zhu M, Ford SP, Nathanielsz PW. Fetal programming of skeletal muscle development in ruminant animals. J Anim Sci. 2010;88(13 Suppl):E51–60.

    Article  CAS  PubMed  Google Scholar 

  53. Du M, Zhao JX, Yan X, Huang Y, Nicodemus LV, Yue W, McCormick RJ, Zhu MJ. Fetal muscle development, mesenchymal multipotent cell differentiation, and associated signaling pathways. J Anim Sci. 2011;89(2):583–90.

    Article  CAS  PubMed  Google Scholar 

  54. La YF, Zhang XX, Li FD, Zhang DY, Li C, Mo FT, Wang WM. Molecular characterization and expression of SPP1, LAP3 and LCORL and their Association with Growth traits in Sheep. Genes-Basel 2019, 10(8).

  55. Ramos Z, Garrick DJ, Blair HT, Vera B, Ciappesoni G, Kenyon PR. Genomic regions Associated with wool, Growth and Reproduction traits in Uruguayan Merino Sheep. Genes-Basel. 2023;14(1):167.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Zhuang X, Xie F, Lin Z, Luo J, Chen T, Xi Q, Zhang Y, Sun J. Effect of mir-493-5p on proliferation and differentiation of myoblast by targeting ANKRD17. Cell Tissue Res. 2023;393(1):119–32.

    Article  CAS  PubMed  Google Scholar 

  57. Girish A, Sutar S, Murthy TPK, Premanand SA, Garg V, Patil L, Shreyas S, Shukla R, Yadav AK, Singh TR. Comprehensive bioinformatics analysis of structural and functional consequences of deleterious missense mutations in the human QDPR gene. J Biomol Struct Dynamics 2023:1–17.

  58. Wang S, Wang Y, Li Y, Xiao F, Guo H, Gao H, Wang N, Zhang H, Li H. Genome-Wide Association Study and Selective Sweep Analysis Reveal the Genetic Architecture of body weights in a Chicken F2 Resource Population. Front Vet Sci. 2022;9:875454.

    Article  PubMed  PubMed Central  Google Scholar 

  59. Zhang Y, Wang Y, Li Y, Wu J, Wang X, Bian C, Tian Y, Sun G, Han R, Liu X, et al. Genome-wide association study reveals the genetic determinism of growth traits in a Gushi-Anka F2 chicken population. Heredity. 2021;126(2):293–307.

    Article  CAS  PubMed  Google Scholar 

  60. Wang WH, Wang JY, Zhang T, Wang Y, Zhang Y, Han K. Genome-wide association study of growth traits in Jinghai Yellow chicken hens using SLAF-seq technology. Anim Genet. 2019;50(2):175–6.

    Article  CAS  PubMed  Google Scholar 

  61. Wei C, Hou D, Feng Y, Li T, Jing Z, Li W, Han R, Li G, Sun G, Tian Y, et al. Molecular characterization and a duplicated 31-bp indel within the LDB2 gene and its associations with production performance in chickens. Gene. 2020;761:145046.

    Article  CAS  PubMed  Google Scholar 

  62. Jin CF, Chen YJ, Yang ZQ, Shi K, Chen CK. A genome-wide association study of growth trait-related single nucleotide polymorphisms in Chinese Yancheng chickens. Genet Mol Res. 2015;14(4):15783–92.

    Article  CAS  PubMed  Google Scholar 

  63. Cha J, Choo H, Srikanth K, Lee SH, Son JW, Park MR, Kim N, Jang GW, Park JE. Genome-Wide Association Study Identifies 12 Loci Associated with Body Weight at Age 8 weeks in Korean native chickens. Genes (Basel) 2021, 12(8).

  64. Makvandi-Nejad S, Hoffman GE, Allen JJ, Chu E, Gu E, Chandler AM, Loredo AI, Bellone RR, Mezey JG, Brooks SA, et al. Four loci explain 83% of size variation in the horse. PLoS ONE. 2012;7(7):e39929.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Bouwman AC, Daetwyler HD, Chamberlain AJ, Ponce CH, Sargolzaei M, Schenkel FS, Sahana G, Govignon-Gion A, Boitard S, Dolezal M, et al. Meta-analysis of genome-wide association studies for cattle stature identifies common genes that regulate body size in mammals. Nat Genet. 2018;50(3):362–7.

    Article  CAS  PubMed  Google Scholar 

  66. Posbergh CJ, Huson HJ. All sheeps and sizes: a genetic investigation of mature body size across sheep breeds reveals a polygenic nature. Anim Genet. 2021;52(1):99–107.

    Article  CAS  PubMed  Google Scholar 

  67. Liu R, Kong F, Xing S, He Z, Bai L, Sun J, Tan X, Zhao D, Zhao G, Wen J. Dominant changes in the breast muscle lipid profiles of broiler chickens with wooden breast syndrome revealed by lipidomics analyses. J Anim Sci Biotechno. 2022;13(1):93.

    Article  CAS  Google Scholar 

  68. Naserkheil M, Mehrban H, Lee D, Park MN. Genome-wide Association Study for Carcass Primal cut yields using single-step bayesian Approach in Hanwoo Cattle. Front Genet. 2021;12:752424.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  69. Wu Z, Wang Z, Xie Y, Liu G, Shang X, Zhan N. Transcriptome and metabolome profiling provide insights into flavonoid synthesis in Acanthus Ilicifolius Linn. Genes (Basel) 2023, 14(3).

  70. Shao M, Shi K, Zhao Q, Duan Y, Shen Y, Tian J, He K, Li D, Yu M, Lu Y et al. Transcriptome analysis reveals the differentially expressed genes Associated with Growth in Guangxi Partridge Chickens. Genes (Basel) 2022, 13(5).

  71. Li X, Yang J, Shen M, Xie XL, Liu GJ, Xu YX, Lv FH, Yang H, Yang YL, Liu CB, et al. Whole-genome resequencing of wild and domestic sheep identifies genes associated with morphological and agronomic traits. Nat Commun. 2020;11(1):2815.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thanks the Institute of Animal Science of CAAS, Chinese Academy of Agricultural Sciences for the support of this experiment.

Funding

This research was supported by the STI 2030-Major Project (2023ZD04072), the Wenchang Chicken superiority characteristic industrial cluster project (WCSCICP20211106), the Project of Sanya Yazhou Bay Science and Technology City (SKJC-2022-PTDX-002), the Special Project for Southern Propagation of the Chinese Academy of Agricultural Sciences (YYLH04), the National chicken industry technology system project (CARS-41), the Study of the Key Genetic Resources [JBGS (2021) 107].

Author information

Authors and Affiliations

Authors

Contributions

KQC, Conceptualization, Data curation and Writing-original draft. RRL, Conceptualization and Validation. LMW, Investigation. XPW, Investigation. HXC, Conceptualization and Data curation. NL, Software. JW, Conceptualization, Validation and Supervision. YXC, Writing-review & editing. GPZ, Writing-review & editing.

Corresponding authors

Correspondence to Yuxiao Chang or Guiping Zhao.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in this study followed the guidelines for the care and use of experimental animals (HNSYY20230203). The research protocol was approved by the Ethics Committee of Sanya Research Institute, Hainan Academy of Agricultural Sciences. All methods were carried out in accordance with relevant guidelines and regulations. This study was carried out in compliance with the ARRIVE guidelines (Animal Research: Reporting of In Vivo Experiments).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Conflict of interest

All authors disclosed no relevant relationships.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cai, K., Liu, R., Wei, L. et al. Genome-wide association analysis identify candidate genes for feed efficiency and growth traits in Wenchang chickens. BMC Genomics 25, 645 (2024). https://doi.org/10.1186/s12864-024-10559-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10559-w

Keywords