Residual feed intake (RFI) is a powerful indicator for energy utilization efficiency and responds to selection. Low RFI selection enables a reduction in feed intake without affecting growth performance. However, the effective variants or major genes dedicated to phenotypic differences in RFI in quality chickens are unclear. Therefore, a genome-wide association study (GWAS) and RNA sequencing were performed on RFI to identify genetic variants and potential candidate genes associated with energy improvement. A lower average daily feed intake was found in low-RFI birds compared to high-RFI birds. The heritability of RFI measured from 44 to 83 d of age was 0.35. GWAS showed that 32 of the significant single nucleotide polymorphisms (SNPs) associated with the RFI (P < 10−4) accounted for 53.01 % of the additive genetic variance. More than half of the effective SNPs were located in a 1 Mb region (16.3–17.3 Mb) of chicken (Gallus gallus) chromosome (GGA) 12. Thus, focusing on this region should enable a deeper understanding of energy utilization. RNA sequencing was performed to profile the liver transcriptomes of four male chickens selected from the high and low tails of the RFI. One hundred and sixteen unique genes were identified as differentially expressed genes (DEGs). Some of these genes were relevant to appetite, cell activities, and fat metabolism, such as CCKAR, HSP90B1, and PCK1. Some potential genes within the 500 Kb flanking region of the significant RFI-related SNPs detected in GWAS (i.e., MGP, HIST1H110, HIST1H2A4L3, OC3, NR0B2, PER2, ST6GALNAC2, and G0S2) were also identified as DEGs in chickens with divergent RFIs. The GWAS findings showed that the 1 Mb narrow region of GGA12 should be important because it contained genes involved in energy-consuming processes, such as lipogenesis, social behavior, and immunity. Similar results were obtained in the transcriptome sequencing experiments. In general, low-RFI birds seemed to optimize energy employment by reducing energy expenditure in cell activities, immune responses, and physical activity compared to eating.

The age of onset of sexual maturity is an important reproductive trait in chickens. In this study, we explored candidate genes associated with sexual maturity and ovary development in chickens. We performed DGE RNA-sequencing analyses of ovaries of pre-laying (P-F-O1, L-F-O1) and laying (P-F-O2, L-F-O2) hens of two sub-breeds of Ningdu Yellow chicken. A total of 3197 genes were identified in the two comparisons, and 966 and 1860 genes were detected exclusively in comparisons of P-F-O1 vs. P-F-O2 and L-F-O1 vs. L-F-O2, respectively. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that genes involved in transmembrane signaling receptor activity, cell adhesion, developmental processes, the neuroactive ligand–receptor interaction pathway, and the calcium signaling pathway were enriched in both comparisons. Genes on these pathways, including growth hormone (GH), integrin subunit beta 3 (ITGB3), thyroid stimulating hormone subunit beta (TSHB), prolactin (PRL), and transforming growth factor beta 3 (TGFB3), play indispensable roles in sexual maturity. As a gene unique to poultry, hen egg protein 21 kDa (HEP21) was chosen as the candidate gene. Differential expression and association analyses were performed. RNA-seq data and qPCR showed that HEP21 was significantly differentially expressed in pre-pubertal and pubertal ovaries. A total of 23 variations were detected in HEP21. Association analyses of single nucleotide polymorphisms (SNPs) in HEP21 and reproductive traits showed that rs315156783 was significantly related to comb height at 84 and 91 days. These results indicate that HEP21 is a candidate gene for sexual maturity in chickens. Our results contribute to a more comprehensive understanding of sexual maturity and reproduction in chickens.

Methyltransferase 3 (METTL3), which has been demonstrated to play a crucial role in a variety of biological processes, is the key enzyme for catalyzing m⁶A modification in RNA. However, the complete protein sequence of METTL3 in quail has not been annotated, and its function in skeletal muscle of quails remains unknown. In the current study, the full-length coding sequence of the quail METTL3 was obtained through the 3′ rapid amplification of cDNA ends (3’ RACE) and its homology with that of other species was predicted based on a generated phylogenetic tree. A Cell Counting Kit-8 assay and flow cytometry in a quail myoblast cell line (QM7) demonstrated that METTL3 promotes myoblast proliferation. The overexpression of METTL3 in QM7 cells significantly increased the expression levels of the myoblast differentiation markers myogenin (MYOG), myogenic differentiation 1 (MYOD1), and myocyte enhancer factor 2C (MEF2C), further demonstrating that METTL3 promotes myoblast differentiation. Additionally, transcriptome sequencing following METTL3 overexpression revealed that METTL3 controls the expression of various genes involved in RNA splicing and the regulation of gene expression, as well as pathways such as the MAPK signaling pathway. Taken together, our findings demonstrated that METTL3 plays a vital function in quail myoblast proliferation and differentiation and that the METTL3-mediated RNA m⁶A modification represents an important epigenetic regulatory mechanism in poultry skeletal muscle development.