Feed contributes to over 60 % of the total production costs in the poultry industry. Increasing feed costs prompt geneticists to include feed intake and efficiency as selection goals in breeding programs. In the present study, we used an F2 chicken population in a genome-wide association study (GWAS) to detect potential genetic variants and candidate genes associated with daily feed intake (FI) and feed efficiency, including residual feed intake (RFI) and feed conversion ratio (FCR). A total of 1534 F2 hens from a White Leghorn and Dongxiang reciprocal cross were phenotyped for feed intake and efficiency between 37 and 40 weeks (FI1, RFI1, and FCR1) and between 57 and 60 weeks (FI2, RFI2, and FCR2), and genotyped using the chicken 600 K single nucleotide polymorphism (SNP) genotyping array. Univariate, bivariate, and conditional genome-wide association studies (GWAS) were performed with GEMMA, a genome-wide efficient mixed model association algorithm. The statistical significance threshold for association was inferred by the simpleM method. We identified eight genomic regions that each contained at least one genetic variant that showed a significant association with FI. Genomic regions on Gallus gallus (GGA) chromosome 4 coincided with known quantitative trait loci (QTL) that affect feed intake of layers. Of particular interest, eight SNPs on GGA1 in the region between 169.23 and 171.55 Mb were consistently associated with FI in both univariate and bivariate GWAS, which explained 3.72 and 2.57 % of the phenotypic variance of FI1 and FI2, respectively. The CAB39L gene can be considered as a promising candidate for FI1. For RFI, a haplotype block on GGA27 harbored a significant SNP associated with RFI2. The major allele of rs315135692 was favorable for a lower RFI, with a phenotypic difference of 3.35 g/day between opposite homozygous genotypes. Strong signals on GGA1 were detected in the bivariate GWAS for FCR. The results demonstrated the polygenic nature of feed intake. GWAS identified novel variants and confirmed a QTL that was previously reported for feed intake in chickens. Genetic variants associated with feed efficiency may be used in genomic breeding programs to select more efficient layers.

Circular RNA (circRNA) is a type of noncoding RNA involved in a variety of biological processes, especially in post-transcriptional regulation. The granulosa cells of follicles play a determining role in ovarian development. However, the function of circRNA in chicken follicles is unclear. To better understand the molecular mechanism underlying follicular development and granulosa cell function, we performed a strategy of second-generation sequencing and linear RNA depletion for granulosa cells from small yellow follicles (SYF, 5–8 mm), the smallest hierarchal follicles (F6, 9–12 mm), and the largest hierarchal follicles (F1, ~ 40 mm). We predicted a total of 11,642 circRNAs that distributed on almost all chromosomes. The majority of the splice lengths of circRNAs were 200–500 nt and mainly produced from intron and CDS regions. During follicle growth, differentially expressed (DE) circRNAs showed dynamic changes which were tissue- and stage-specific. The host genes of DE circRNAs were functionally enriched in GTPase activity and several pathways involved in reproduction. Moreover, bioinformatic prediction analysis for circRalGPS2 demonstrated that circRNAs from the same genes may share common miRNA to act as a sponge. The predicted target genes were enriched in various biological processes including cognition, cell communication, and regulation of signaling, and several pathways related to reproduction such as tight junction, oocyte meiosis, progesterone-mediated oocyte maturation, and GnRH signaling. This study provides a starting point for further experimental investigations into chicken circRNAs and casts a light on the understanding of follicle development.

As a major economic trait in chickens, egg weight (EW) receives widespread interests in breeding, production and consumption. However, limited information is available for underlying genetic architecture of longitudinal trend in EW. Herein, we measured EWs at nine time points from onset of laying to 60 week of age, and conducted comprehensive genome-wide association studies (GWAS) in 1,534 F2 hens derived from reciprocal crosses between White Leghorn and Dongxiang chickens. Egg weights at all ages except the first egg weight (FEW) exhibited high SNP-based heritability estimates (0.47 ~ 0.60). Strong pair-wise genetic correlations (0.77 ~ 1.00) were found among all EWs. Nine separate univariate genome-wide screens suggested 73 signals showing significant associations with longitudinal EWs. After multivariate and conditional analyses, four variants on three chromosomes remained independent contributions. The minor alleles at two loci exerted consistent and positive substitution effects on EWs, and other two were negative. The four loci together accounted for 3.84 % of the phenotypic variance for FEW and 7.29 ~ 11.06 % for EWs from 32 to 60 week of age. We obtained five candidate genes, of which NCAPG harbors a non-synonymous SNP (rs14491030) causing a valine-to-alanine amino-acid substitution. Genome partitioning analysis indicated a strong linear correlation between the variance explained by each chromosome and its length, which provided evidence that EW follows a highly polygenic nature of inheritance. Identification of significant genetic causes that together implicate EWs at different ages will greatly advance our understanding of the genetic basis behind longitudinal EWs, and would be helpful to illuminate the future breeding direction on how to select desired egg size.

Eggshell is subject to quality loss with aging process of laying hens, and damaged eggshells result in economic losses of eggs. However, the genetic architecture underlying the dynamic eggshell quality remains elusive. Here, we measured eggshell quality traits, including eggshell weight (ESW), eggshell thickness (EST) and eggshell strength (ESS) at 11 time points from onset of laying to 72 weeks of age and conducted comprehensive genome-wide association studies (GWAS) in 1534 F2 hens derived from reciprocal crosses between White Leghorn (WL) and Dongxiang chickens (DX). ESWs at all ages exhibited moderate SNP-based heritability estimates (0.30 ~ 0.46), while the estimates for EST (0.21 ~ 0.31) and ESS (0.20 ~ 0.27) were relatively low. Eleven independent univariate genome-wide screens for each trait totally identified 1059, 1026 and 1356 significant associations with ESW, EST and ESS, respectively. Most significant loci were in a region spanning from 57.3 to 71.4 Mb of chromosome 1 (GGA1), which together account for 8.4 ~ 16.5 % of the phenotypic variance for ESW from 32 to 72 weeks of age, 4.1 ~ 6.9 % and 2.95 ~ 16.1 % for EST and ESS from 40 to 72 weeks of age. According to linkage disequilibrium (LD) and conditional analysis, the significant SNPs in this region were in extremely strong linkage disequilibrium status. Ultimately, two missense SNPs in GGA1 and one in GGA4 were considered as promising loci on three independent genes including ITPR2, PIK3C2G, and NCAPG. The homozygotes of advantageously effective alleles on PIK3C2G and ITPR2 possessed the best eggshell quality and could partly counteract the negative effect of aging process. NCAPG had certain effect on eggshell quality for young hens. Identification of the promising region as well as potential candidate genes will greatly advance our understanding of the genetic basis underlying dynamic eggshell quality and has the practical significance in breeding program for the improvement of eggshell quality, especially at the later part of laying cycle.

The growth traits are important traits in chickens. Compared to white feather broiler breeds, Chinese local broiler breeds have a slow growth rate. The main genes affecting the growth traits of local chickens in China are still unclear and need to be further explored. This experiment used fast-growth and slow-growth groups of the Jinghai Yellow chicken as the research objects. Three males and three females with similar body weights were selected from the two groups at four weeks old and eight weeks old, respectively, with a total of 24 individuals selected. After slaughter, their chest muscles were taken for transcriptome sequencing. In the differentially expressed genes screening, all of the genes obtained were screened by fold change ≥ 2 and false discovery rate (FDR) < 0.05. For four-week-old chickens, a total of 172 differentially expressed genes were screened in males, where there were 68 upregulated genes and 104 downregulated genes in the fast-growth group when compared with the slow-growth group. A total of 31 differentially expressed genes were screened in females, where there were 11 upregulated genes and 20 downregulated genes in the fast-growth group when compared with the slow-growth group. For eight-week-old chickens, a total of 37 differentially expressed genes were screened in males. The fast-growth group had 28 upregulated genes and 9 downregulated genes when compared with the slow-growth group. A total of 44 differentially expressed genes were screened in females. The fast-growth group had 13 upregulated genes and 31 downregulated genes when compared with the slow-growth group. Through gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, many genes were found to be related to cell proliferation and differentiation, muscle growth, and cell division such as SNCG, MCL1, ARNTL, PLPPR4, VAMP1, etc. Real-time PCR results were consistent with the RNA-Seq data and validated the findings. The results of this study will help to understand the regulation mechanism of the growth and development of Jinghai Yellow chicken and provide a theoretical basis for improving the growth rate of Chinese local chicken breeds.