Transcriptomic Analyses of the Hypothalamic-Pituitary-Gonadal Axis Identify Candidate Genes Related to Egg Production in Xinjiang Yili Geese

Transcriptome analyses of shell color and egg production traits between the uteruses of blue-green eggshell chickens and Hy-Line brown layers

Blue-green eggs exhibit unique shell color; however, compared to commercial layers, blue-green eggshell chickens have lower egg production and lack uniform shell colors. Aiming to confirm the molecular mechanisms that affect shell color and egg production, this study collected the uteruses of 12 blue-green eggshell chickens (BG group) and six Hy-Line layers (Brown group), which had significantly different shell color indexes (SCI) and egg numbers at 300 days of age (EN300). Transcriptome sequencing and comparative analyses were subsequently performed. BG hens were divided into two groups for comparative analysis (BGblue vs. BGgreen and BGlow vs. BGhigh, respectively), based on the differences in SCI and EN300, respectively. The result of weighted gene co-expression network (WGCNA) analysis showed that the sequenced and mapped 16,785 genes were clustered into 18 modules, among which six modules with a total of 4270 genes were highly correlated with SCI and EN300 traits. Five hundred and eleven differentially expressed genes (DEGs) belonged to the six key modules. Through KEGG mapping, GO enrichment, and Cytoscape network analysis, nine Hub genes were tightly associated with SCI and EN300. The up-regulated genes were CCR2, CCR8, CD40LG, IL18RAP, INHBA, and P2RY13, while the down-regulated genes were ABCA13, ADCY2, and GRM1. Co-analyses with the results of comparisons between the BG subgroups revealed that the expression of solute carrier (SLC) proteins and ABC transporters were highly related to eggshell color, while cytokine-cytokine receptor interactions and neuroactive ligand-receptor interactions were key pathways affecting egg production. The expression of extracellular cytokines and membrane receptors were significantly up-regulated in low-yield chickens. The candidate DEGs and pathways found in the study will assist in clarifying the molecular mechanisms affecting shell color and egg production, and improve the breeding of blue-green eggshell chickens.

The hypothalamic transcriptome reveals the importance of visual perception on the egg production of Wanxi white geese

Egg performance significantly impacts the development of the local goose industry. The hypothalamus plays an essential role in the egg production of birds. However, few potential candidate genes and biological functions related to egg production in geese have been identified in hypothalamus tissue. In this study, 115 geese were raised and observed for 5 months during the laying period. To understand the regulation mechanism of egg production, the hypothalamus transcriptome profiles of these geese were sequenced using RNA-seq. The hypothalamus samples of four high egg production (HEP) and four low egg production (LEP) geese were selected and collected, respectively. A total of 14,679 genes were identified in the samples. After multiple bioinformatics analyses, Gene Ontology (GO) annotations indicated that genes related to egg production were mainly enriched in biological processes of "response to light stimulus," "sensory system development," and "visual perception." Six potential candidate genes (PDE6C, RHO, MFRP, F2, APOB, and IL6) based on their corresponding GO terms and interaction networks were identified. These identified candidate genes can be used as selection markers to improve the egg production of Wanxi white geese. Our study highlights how visual perception may affect the regulation of geese egg production.

Genome-wide association study identified candidate genes for egg production traits in the Longyan Shan-ma duck

Egg production is an important economic trait in layer ducks and understanding the genetics basis is important for their breeding. In this study, a genome-wide association study (GWAS) for egg production traits in 303 female Longyan Shan-ma ducks was performed based on a genotyping-by-sequencing strategy. Sixty-two single nucleotide polymorphisms (SNPs) associated with egg weight traits were identified (P < 9.48 × 10-5), including 8 SNPs at 5% linkage disequilibrium (LD)-based Bonferroni-corrected genome-wide significance level (P < 4.74 × 10-6). One hundred and nineteen SNPs were associated with egg number traits (P < 9.48 × 10-5), including 13 SNPs with 5% LD-based Bonferroni-corrected genome-wide significance (P < 4.74 × 10-6). These SNPs annotated 146 target genes which contained known candidate genes for egg production traits, such as prolactin and prolactin releasing hormone receptor. This study identified that these associated genes were significantly enriched in egg production-related pathways (P < 0.05), such as the oxytocin signaling, MAPK signaling, and calcium signaling pathways. It was notable that 18 genes were differentially expressed in ovarian tissues between higher and lower egg production in Shan-ma ducks. The identified potential candidate genes and pathways provide insight into the genetic basis underlying the egg production trait of layer ducks.

The Omics Revolution in Understanding Chicken Reproduction: A Comprehensive Review

Omics approaches have significantly contributed to our understanding of several aspects of chicken reproduction. This review paper gives an overview of the use of omics technologies such as genomics, transcriptomics, proteomics, and metabolomics to elucidate the mechanisms of chicken reproduction. Genomics has transformed the study of chicken reproduction by allowing the examination of the full genetic makeup of chickens, resulting in the discovery of genes associated with reproductive features and disorders. Transcriptomics has provided insights into the gene expression patterns and regulatory mechanisms involved in reproductive processes, allowing for a better knowledge of developmental stages and hormone regulation. Furthermore, proteomics has made it easier to identify and quantify the proteins involved in reproductive physiology to better understand the molecular mechanisms driving fertility, embryonic development, and egg quality. Metabolomics has emerged as a useful technique for understanding the metabolic pathways and biomarkers linked to reproductive performance, providing vital insights for enhancing breeding tactics and reproductive health. The integration of omics data has resulted in the identification of critical molecular pathways and biomarkers linked with chicken reproductive features, providing the opportunity for targeted genetic selection and improved reproductive management approaches. Furthermore, omics technologies have helped to create biomarkers for fertility and embryonic viability, providing the poultry sector with tools for effective breeding and reproductive health management. Finally, omics technologies have greatly improved our understanding of chicken reproduction by revealing the molecular complexities that underpin reproductive processes.

Transcriptomic analysis of the HPG axis-related tissues reveals potential candidate genes and regulatory pathways associated with testicular size in Hu sheep

Testicular size is an excellent proxy for selecting high-fertility rams. The hypothalamus-pituitary-gonadal (HPG) axis plays an important role in regulating reproductive capacity in vertebrates, while key genes and regulatory pathways within the HPG axis associated with testicular size remain largely unknown in sheep. This study comprehensively compared the transcriptomic profiles in the hypothalamus, pituitary and testis of rams after sexual maturity between the large-testis group (LTG, testicular weight = 454.29 ± 54.24 g) and the small-testis group (STG, testicular weight = 77.29 ± 10.76 g). In total, 914, 795 and 10518 differentially expressed genes (DEGs) were identified in the hypothalamus, pituitary and testis between LTG and STG, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DEGs were mainly involved in the biological processes of reproduction, biological regulation, and development process. Notably, the neuroactive ligand-receptor interaction and cAMP signaling pathways, commonly enriched by the DEGs in the hypothalamus and pituitary between two groups, were considered as two key signal pathways regulating testicular development through the HPGs axis. Weighted gene co-expression network analysis (WGCNA) identified two modules that were significantly associated with testicular size, and 97 key genes were selected with high module membership (MM) and gene significance (GS) in these two modules. Finally, a protein-protein interaction (PPI) network was constructed, and ten genes with the highest degree were represented as hub genes, including FOS, NPY, SST, F2, AGT, NTS, OXT, EDN1, VIP and TAC1. Taken together, these results provide new insights into the molecular mechanism underlying the HPG axis regulating testicular size of Hu sheep.

Analysis of Transcriptomic Differences in the Ovaries of High- and Low-Laying Ducks

The egg-laying performance of Shan Ma ducks (Anas Platyrhynchos) is a crucial economic trait. Nevertheless, limited research has been conducted on the egg-laying performance of this species. We examined routine blood indicators and observed higher levels of metabolic and immune-related factors in the high-egg-production group compared with the low-egg-production group. Furthermore, we explored the ovarian transcriptome of both high- and low-egg-production groups of Shan Ma ducks using Illumina NovaSeq 6000 sequencing. A total of 1357 differentially expressed genes (DEGs) were identified, with 686 down-regulated and 671 up-regulated in the high-egg-production (HEP) ducks and low-egg-production (LEP) ducks. Several genes involved in the regulation of ovarian development, including neuropeptide Y (NPY), cell cycle protein-dependent kinase 1 (CDK1), and transcription factor 1 (E2F1), exhibited significant differential expressions at varying stages of egg production. Pathway functional analysis revealed that the DEGs were primarily associated with the steroid biosynthesis pathway, and the neuroactive ligand-receptor interaction pathway exhibited higher activity in the HEP group compared to the LEP group. This study offers valuable information about and novel insights into high egg production.

Liver Transcriptome Profiling Identifies Key Genes Related to Lipid Metabolism in Yili Geese

The Yili goose is the only indigenous goose breed that originates from Anser anser in China, known for its adaptability, strong flying ability, and tender meat with a low body lipid content. The liver plays a crucial role in lipid and glucose metabolism, including the intake, secretion, transportation, and storage of fatty acids (FAs). In this study, RNA-sequencing (RNA-seq) technology was performed to analyze the liver differentially expressed genes of Yili geese and their hybrid geese to investigate differences in liver lipid and glucose metabolism. A total of 452 differentially expressed genes (Q-value < 0.05) were identified. Notably, in KEGG enrichment analysis, four pathways (Q-value < 0.05) were enriched to be associated with lipid and glucose metabolism, including the metabolic pathway, PI3K-Akt signaling pathway, glycolysis/gluconeogenesis, and steroid biosynthesis. This study provides insights into potential candidate genes and metabolic pathways that affect the liver lipid metabolism of Yili goose. These findings provide a better understanding of animal liver lipid deposition and metabolism.

Identification of potential candidate genes and regulatory pathways related to reproductive capacity in hypothalamus and pituitarium of male ducks (Anas platyrhynchos) by differential transcriptome analysis

The improvement of reproductive capacity of poultry is important for the poultry industry. The existing studies on reproductive capacity mainly focus on the testis tissue, but few reports on regulationary effect of brain neuroendocrime on reproductive capacity have been available. The hypothalamus-pituitarium-gonad (HPG) axis is an important pathway regulating spermatogenesis and sexual behavior. This study analyzed the gene expression in the hypothalamus and pituitary tissues of male ducks in high-semen-quality group (DH), low-semen-quality group (DL), and non-response group (DN) by RNA-sequencing. A total of 1980 differentially expressed genes (DEGs) were identified, and significantly less DEGs were found in pituitary gland than in hypothalamus. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that these DEGs were mainly enriched in nerve-related and synapse-related biological processes, mitochondrial inner membrane formation pathway, and ribosome structure pathway. Notably, the neuroactive ligand-receptor interaction pathway significantly enriched in all three comparisons (DH vs. DL, DH vs. DN, and DL vs. DN) was related to different reproductive performance such as semen quality and sexual response. Furthermore, six genes, including POMC, CPLX2, HAPLN2, EGR4, TOX3, and MSH4, were identified as candidate genes regulating reproductive capacity. Our findings provide new insights into the regulation mechanisms underlying the reproductive performance of male poultry, and offer a valuable reference for duck breeding programs aimed at promoting reproductive capacity.

Identification of differentially expressed genes and signaling pathways in Gaoyou duck ovary at different physiological stages

Introduction

Gaoyou duck is famous in China and abroad for its good production of double-yolk eggs. However, there has been no systematic research on the egg-laying characteristics of the Gaoyou duck, which limits the development and utilization of breed resource.

Methods

To identify the essential genes related to ovarian development, the transcriptome profiles of the ovaries of Gaoyou ducks at different physiological stages were analyzed. The transcriptome profiles of the ovaries of Gaoyou ducks at 150 d (before laying), 240 d (egg laying) and 500 d (nesting) were constructed, and the differentially expressed genes (DEGs) underwent GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses.

Results

The 6 randomly selected DEGs were verified by real-time fluorescent quantitative PCR that their relative expression was consistent with the transcriptional expression profile. Furthermore, KEGG analysis found that 8 candidate signaling pathways were essential for ovarian development, including the MAPK signaling pathway, Progesterone-mediated oocyte maturation, Cell adhesion molecules (CAMs), NOD-like receptor signaling pathway, ECM-receptor interaction, Focal adhesion, TGF-beta signaling path-way and Phagosome. Finally, 5 key DEGs were identified to participate in ovarian development, including TGIF1, TGFBR2, RAF1, PTK2 and FGF10.

Discussion

Our findings reveal the mechanisms under-lying the molecular regulation of related genes in Gaoyou duck ovarian development.

Gonadal Transcriptome Analysis and Sequence Characterization of Sex-Related Genes in Cranoglanis bouderius

In China, the Cranoglanis bouderius is classified as a national class II-protected animal. The development of C. bouderius populations has been affected by a variety of factors over the past few decades, with severe declines occurring. Considering the likelihood of continued population declines of the C. bouderius in the future, it is critical to investigate the currently unknown characteristics of gonadal differentiation and sex-related genes for C. bouderius conservation. In this study, the Illumina sequencing platform was used to sequence the gonadal transcriptome of the C. bouderius to identify the pathways and genes related to gonadal development and analyze the expression differences in the gonads. A total of 12,002 DEGs were identified, with 7220 being significantly expressed in the ovary and 4782 being significantly expressed in the testis. According to the functional enrichment results, the cell cycle, RNA transport, apoptosis, Wnt signaling pathway, p53 signaling pathway, and prolactin signaling pathway play important roles in sex development in the C. bouderius. Furthermore, the sequence characterization and evolutionary analysis revealed that AMH, DAX1, NANOS1, and AR of the C. bouderius are highly conserved. Specifically, the qRT-PCR results from various tissues showed significant differences in AMH, DAX1, NANOS1, and AR expression levels in the gonads of both sexes of C. bouderius. These analyses indicated that AMH, DAX1, NANOS1, and AR may play important roles in the differentiation and development of C. bouderius gonads. To our best knowledge, this study is the first to analyze the C. bouderius gonadal transcriptome and identify the structures of sex-related genes, laying the foundation for future research.

RNA-seq identifies differentially expressed genes involved in csal1 overexpression in granulosa cells of prehierarchical follicles in Chinese Dagu hens

The transcription factor csal1 is an important molecule that plays a critical regulatory function in ovarian follicle development, as confirmed by our previous data. However, the candidate genes of csal1 and its regulatory mechanism remain poorly understood in the granulosa cells (GCs) of chicken prehierarchical follicles (PFs). Six transcriptomes of csal1 and empty vector were analyzed in Chinese Dagu hens by RNA sequencing. Six cDNA libraries were constructed, with more than 42 million clean reads and 16,779 unigenes. Of these 16,779 unigenes, 2,762 differentially expressed genes (DEGs) were found in GCs, including 1,605 upregulated and 1,157 downregulated unigenes. Fourteen genes, including BMP5, TACR2, AMH, PLAG1, MYOD1, BOP1, SIPA1, NOTCH1, BCL2L1, SOX9, ADGRA2, WNT5A, SLC7A11, and GATAD2B, were related to GC proliferation and differentiation, hormone production, ovarian follicular development, regulation of reproductive processes, and signaling pathways in the PFs. Further analysis demonstrated the DEGs in GCs of ovarian follicles were enriched in neuroactive ligand-receptor interaction, cell adhesion molecules, and pathways related to cytochrome P450, indicating a critical function for csal1 in the generation of egg-laying features by controlling ovarian follicle development. For the first time, the current study represents the transcriptome analysis with ectopic csal1 expression. These findings provide significant evidence for investigating the molecular mechanism by which csal1 controls PF development in the hen ovary.

MicroRNAs with non-additive expression in the ovary of hybrid hens target genes enriched in key reproductive pathways that may influence heterosis for egg laying traits

Heterosis has been extensively exploited in chicken breeding to improve laying traits in commercial hybrid stock. However, the molecular mechanisms underlying it remains elusive. This study characterizes the miRNAome in the pre-hierarchical follicles of purebred and hybrid laying hens, and investigate the functions of miRNAs with non-additive expression in the pre-hierarchical follicles as they modulate heterosis for egg number and clutch size. To achieve that aim, White Leghorn and Rhode Island Red chicken lines were reciprocally crossed to generate hybrids. The crossbreds demonstrated heterosis for egg number and clutch size, and pre-hierarchical follicles from 4 birds of each genotype were collected at 53 weeks of age. Mode of miRNA expression was characterized after miRNA sequencing. A total of 50 miRNAs including 30 novel ones, were found to exhibit non-additive expression. Dominance was the predominant mode of expression exhibited by majority of the miRNAs. Functional analysis of target genes of the known miRNAs with non-additive expression revealed Gene Ontology terms related to regulation of transcription, metabolic processes and gene expression. KEGG and REACTOME pathways including hedgehog, cellular senescence, wnt, TGF-β, progesterone-mediated oocyte maturation, oocyte meiosis, GnRH signaling, signal transduction and generic transcription, which can be linked to primordial follicle activation, growth and ovulation, were significantly enriched by target genes of miRNAs with non-additive expression. Majority of the genes enriched in these biological pathways were targeted by gga-miR-19a, gga-miR-19b, gga-miR-375, gga-miR-135a, and gga-miR-7 and 7b, thus, revealing their synergistic roles in enhancing processes that could influence heterosis for egg number and clutch size in hybrid hens.

Comprehensive analysis of differentially expressed profiles of mRNA, lncRNA, and miRNA of Yili geese ovary at different egg-laying stages

Background

The development of the ovaries is an important factor that affects egg production performance in geese. Ovarian development is regulated by genes that are expressed dynamically and stage-specifically. The transcriptome profile analysis on ovarian tissues of goose at different egg laying stages could provide an important basis for screening and identifying key genes regulating ovarian development.

Results

In this study, 4 ovary tissues at each breeding period of pre-laying (PP), laying (LP), and ceased-laying period (CP), respectively, with significant morphology difference, were used for RNA extraction and mRNAs, lncRNAs, and miRNAs comparison in Yili geese. CeRNA regulatory network was constructed for key genes screening. A total of 337, 1136, and 525 differentially expressed DE mRNAs, 466, 925, and 742 DE lncRNAs and 258, 1131 and 909 DE miRNAs were identified between PP and LP, between CP and LP, and between CP and PP groups, respectively. Functional enrichment analysis showed that the differentially expressed mRNAs and non-coding RNA target genes were mainly involved in the cell process, cytokine-cytokine receptor interaction, phagosome, calcium signaling pathway, steroid biosynthesis and ECM-receptor interaction. Differential genes and non-coding RNAs, PDGFRB, ERBB4, LHCGR, MSTRG.129094.34, MSTRG.3524.1 and gga-miR-145–5p, related to reproduction and ovarian development were highly enriched. Furthermore, lncRNA-miRNA-mRNA regulatory networks related to ovary development were constructed.

Conclusions

Our study found dramatic transcriptomic differences in ovaries of Yili geese at different egg-laying stages, and a differential lncRNA-miRNA-mRNA regulatory network related to cell proliferation, differentiation and apoptosis and involved in stromal follicle development were established and preliminarily validated, which could be regarded as a key regulatory pathway of ovarian development in Yili geese.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08774-4.

Ovarian Transcriptomic Analysis of Ninghai Indigenous Chickens at Different Egg-Laying Periods

Egg production is an essential indicator of poultry fertility. The ovary is a crucial organ involved in egg production; however, little is known about the key genes and signaling pathways involved in the whole egg-laying cycle of hens. In order to explore the mechanism of egg production at different stages of the egg-laying process, ovarian tissues from four chickens were randomly selected for transcriptome analysis at each of the three ages (145 d, 204 d, and 300 d in the early, peak, and late stages of egg laying). A total of 12 gene libraries were constructed, and a total of 8433 differential genes were identified from NH145d vs. NH204d, NH145d vs. NH300d and NH300d vs. NH204d (Ninghai 145-day-old, Ninghai 204-day-old, and Ninghai 300-day-old), with 1176, 1653 and 1868 up-regulated genes, and 621, 1955 and 1160 down-regulated genes, respectively. In each of the two comparison groups, 73, 1004, and 1030 differentially expressed genes were found to be co-expressed. We analyzed the differentially expressed genes and predicted nine genes involved in egg production regulation, including LRP8, BMP6, ZP4, COL4A1, VCAN, INHBA, LOX, PTX3, and IHH, as well as several essential egg production pathways, such as regulation adhesion molecules (CAMs), calcium signaling pathways, neuroactive ligand–receptor interaction, and cytokine–cytokine receptor interaction. Transcriptional analysis of the chicken ovary during different phases of egg-lay will provide a useful molecular basis for study of the development of the egg-laying ovary.

Transcriptome analysis reveals key genes and pathways associated with egg production in Nandan-Yao domestic chicken

BACKGROUND

Egg production is a very important economic trait in chicken breeding, but its molecular mechanism is unclear until now. Nandan-Yao chicken (Gallus gallus domesticus) is a native breed in Guangxi province, China, which is famous for good meet quality, but with low egg production.

METHODS

To explore the molecular regulation related to egg production, high egg production (HEP) and low egg production (LEP) were divided according to the total egg number at 55 weeks, and the concentration of serum sex hormones was tested to evaluate the physiological function of ovary and uterus. RNA sequencing (RNA-Seq) was used to explore the transcriptome from the ovary and uterus of Nandan-Yao chicken.

RESULTS

The levels of serum sex hormone showed that concentrations of estradiol (E2), follicle-stimulating hormone (FSH), and luteotropic hormone (LH) were significantly higher in HEP than those in LEP (P < 0.01), while the concentration of testosterone (T) was significantly lower in HEP (P < 0.01). RNA-Seq analysis identified 901 and 2763 differentially expressed genes (DEGs) in ovary and uterus, respectively. Enrichment analysis showed that DEGs were significantly involved in the regulation of tight junction in the ovary (P < 0.05), while in uterus, DEGs were mainly enriched in the phagosome, ECM-receptor interaction, cell adhesion molecules (CAMs), focal adhesion, cardiac muscle contraction, cytokine-cytokine receptor interaction, and the regulation of MAPK signaling pathway (P < 0.05). Protein network interaction and function analyses revealed that FN1, FGF7, SOX2 identified from the ovary, and UQCRH, COX5A, FN1 from the uterus might be key candidate genes for egg production in Nandan-Yao chicken.

CONCLUSIONS

Our study provided key candidate genes and pathways involved in the egg-laying process of Nandan-Yao chicken and could help to further understand the molecular mechanisms of chicken reproduction.

Comparative analysis of hypothalamus transcriptome between laying hens with different egg-laying rates

Egg-laying performance is one of the most important economic traits in the poultry industry. Commercial layers can lay one egg almost every day during their peak-laying period. However, many Chinese indigenous chicken breeds show a relatively low egg-laying rate, even during their peak-laying period. To understand what makes the difference in egg production, we compared the hypothalamus transcriptome profiles of Lushi blue-shelled-egg chickens (LBS), a Chinese indigenous breed with low egg-laying rate and Rhode Island Red chickens (RIR), a commercial layer with relatively high egg-laying rate using RNA-seq. A total of 753 differentially expressed genes (DEGs) were obtained. Of these DEGs, 38 genes were enriched in 2 Gene Ontology (GO) terms, namely reproduction term and the reproductive process term, and 6 KEGG pathways, namely Wnt signaling pathway, Oocyte meiosis, GnRH signaling pathway, Thyroid hormone signaling pathway, Thyroid hormone synthesis and MAPK signaling pathway, which have been long known to be involved in egg production regulation. To further determine the core genes from the 38 DEGs, protein-protein interaction (PPI) network, co-expression network and transcriptional regulatory network analyses were carried out. After integrated analysis and experimental validation, 4 core genes including RAC1, MRE11A, MAP7 and SOX5 were identified as the potential core genes that are responsible for the laying-rate difference between the 2 breeds. These findings paved the way for future investigating the mechanism of egg-laying regulation and enriched the chicken reproductive regulation theory.

Transcriptome Analysis of Ovarian Follicles Reveals Potential Pivotal Genes Associated With Increased and Decreased Rates of Chicken Egg Production

Egg production is an important economic trait in the commercial poultry industry. Ovarian follicle development plays a pivotal role in regulation of laying hen performance and reproductive physiology. However, the key genes and signaling pathways involved in the various-stages of laying hen follicular development remain poorly understood. In this study, transcriptomes of ovarian follicles at three developmental stages, the large white follicle (LWF), small yellow follicle (SYF), and large yellow follicle (LYF), were comparatively analyzed in hens with high (HR) and low (LR) egg-laying rates by RNA-sequencing. Eighteen cDNA libraries were constructed and a total of 236, 544, and 386 unigenes were significantly differentially expressed in the LWF, SYF, and LYF follicles of HR and LR hens, respectively. Among them, 47 co-transcribed differentially expressed genes (DEGs) in LWF and SYF, 68 co-expressed DEGs in SYF and LYF, and 54 co-expressed DEGs in LWF and LYF were mined. Thirteen co-expressed DEGs were found in LWF, SYF, and LYF follicles. Eighteen candidate genes, including P2RX1, CAB39L, BLK, CSMD3, GPR65, ADRB2, CSMD1, PLPP4, ATF3, PRLL, STMN3, RORB, PIK3R1, PERP1, ACSBG1, MRTO4, CDKN1A, and EDA2R were identified to be potentially related to egg production. Furthermore, Kyoto Encyclopedia of Genes and Genomes analysis indicated neuroactive ligand-receptor interaction, cell adhesion molecules, peroxisome proliferator-activated receptor pathway, and cAMP signaling pathway might elicit an important role in formation of egg-laying traits by influencing ovarian follicle development. This study represents the first transcriptome analysis of various-sized follicles between HR and LR hens. These results provide useful molecular evidence for elucidating the genetic mechanism underlying ovarian follicle development associated with egg production in chicken.