Transcriptome comparative analysis of ovarian follicles reveals the key genes and signaling pathways implicated in hen egg production

DUSP8-attenuated ERK1/2 signaling mediates lipogenesis and steroidogenesis in chicken granulosa cells

The lipogenesis and steroidogenesis of granulosa cells are crucial during follicular development, yet it remains unclear whether dual-specificity phosphatase 8 (DUSP8) is involved. In this study, the specific role of DUSP8 in lipogenesis and steroidogenesis was investigated through culturing chicken granulosa cells in vitro. The results revealed that the expression levels of adipogenic genes were elevated after DUSP8 overexpression and reduced after knockdown. The same was observed for lipid deposition in granulosa cells. Meanwhile, the steroidogenic gene expression and progesterone synthesis were promoted after DUSP8 overexpression and inhibited after knockdown. In addition, we also found that DUSP8 blocked the phosphorylation of extracellular regulatory kinase 1/2 (ERK1/2). Based on the previous results that activated ERK1/2 signaling inhibited lipid deposition and progesterone synthesis in chicken granulosa cells, we demonstrated that DUSP8 promoted lipid deposition and progesterone synthesis through mediating the ERK1/2 signaling pathway. The results will improve our understanding of the molecular regulatory mechanisms regarding lipid metabolism and progesterone synthesis in chicken granulosa cells.

CircRNAs involved in the red light of effect on follicle selection in pigeons

Red light (RL) can enhance egg production in poultry. CircRNAs play a crucial role by serving as transcriptional regulators. However, their role in influencing follicle development in White King pigeons remains unexplored. In this study, 54 paired White King pigeons were chosen and divided into RL and white light (WL) groups, each with 3 subgroups. The egg production of paired pigeons in each replicate was recorded for 45 d, and the characteristics of follicle development were monitored during the laying interval (LI). The granulosa cell layer from follicles of the second-largest follicle (F2) was collected, and high-throughput sequencing was performed to elucidate the molecular mechanism of follicle development in pigeons. The study confirmed that RL enhances egg production in pigeons. Additionally, under RL, the F2 follicle was selected, while under WL, small follicles were kept on the third day (LI3). A total of 5,510 circRNAs were identified across all samples, revealing differentially expressed circRNAs (DECs) in various comparisons: 627 in RF1 vs. WF1, 900 in RF2 vs. WF2, 606 in RF1 vs. RF2, and 937 in WF1 vs. WF2. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that host genes of DECs were enriched in pathways like steroid hormone biosynthesis, oocyte meiosis, GnRH signaling pathway, and apoptosis pathway. Moreover, circRNA_5497, circRNA_2016, and circRNA_3328 were common DECs across 4 groups, sharing miRNA binding sites with follicle selection-associated genes. In conclusion, our findings suggest that RL promotes egg production by stimulating follicle selection during LI, offering insights into the regulatory mechanisms of circRNAs in follicle selection under RL. This knowledge can help enhance the reproductive performance of pigeons.

Unraveling the mysteries of chicken proteomics: Insights into follicle development and reproduction

Chicken proteomics is a valuable method for comprehending the many mechanisms involved in follicle growth and reproduction in birds. This study offers a thorough summary of the latest progress in chicken proteomics research, specifically highlighting the knowledge obtained regarding follicle development and reproductive physiology. Proteomic studies have revealed essential proteins and pathways that play a role in follicle development, including those that control oocyte size, maturation, and ovulation. Proteomic investigations have provided insight into the molecular pathways that govern reproductive processes. By utilizing advanced proteomic technologies, including mass spectrometry and protein microarray analysis, we have been able to identify and measure many proteins in chicken follicles at their different developmental stages. The utilization of proteomic methods has enabled the identification of previously unknown biomarkers for reproductive efficiency that expedited the creation of innovative diagnostic instruments for monitoring reproductive health in chicken. Chicken proteomics not only offers insights into follicle growth and reproduction but also uncovers the effects of environmental influences on reproductive function. This provides new opportunities for exploring the molecular pathways that cause these effects. The integration of current data with upcoming proteomic technologies offers the potential for innovative strategies to enhance chicken reproduction.

The potential influence of genome-wide adaptive divergence on conservation translocation outcome in an isolated greater sage-grouse population

Conservation translocations are an important conservation tool commonly employed to augment declining or reestablish extirpated populations. One goal of augmentation is to increase genetic diversity and reduce the risk of inbreeding depression (i.e., genetic rescue). However, introducing individuals from significantly diverged populations risks disrupting coadapted traits and reducing local fitness (i.e., outbreeding depression). Genetic data are increasingly more accessible for wildlife species and can provide unique insight regarding the presence and retention of introduced genetic variation from augmentation as an indicator of effectiveness and adaptive similarity as an indicator of source and recipient population suitability. We used 2 genetic data sets to evaluate augmentation of isolated populations of greater sage-grouse (Centrocercus urophasianus) in the northwestern region of the species range (Washington, USA) and to retrospectively evaluate adaptive divergence among source and recipient populations. We developed 2 statistical models for microsatellite data to evaluate augmentation outcomes. We used one model to predict genetic diversity after augmentation and compared these predictions with observations of genetic change. We used the second model to quantify the amount of observed reproduction attributed to transplants (proof of population integration). We also characterized genome-wide adaptive divergence among source and recipient populations. Observed genetic diversity (HO = 0.65) was higher in the recipient population than predicted had no augmentation occurred (HO = 0.58) but less than what was predicted by our model (HO = 0.75). The amount of shared genetic variation between the 2 geographically isolated resident populations increased, which is evidence of periodic gene flow previously assumed to be rare. Among candidate adaptive genes associated with elevated fixation index (FST) (143 genes) or local environmental variables (97 and 157 genes for each genotype-environment association method, respectively), we found clusters of genes with related functions that may influence the ability of transplants to use local resources and navigate unfamiliar environments and their reproductive potential, all possible reasons for low genetic retention from augmentation.

Integrating genomics and transcriptomics to identify candidate genes for high-altitude adaptation and egg production in Nixi chicken

1. This study combined genome-wide selection signal analysis with RNA-sequencing to identify candidate genes associated with high altitude adaptation and egg production performance in Nixi chickens (NXC).2. Based on the whole-genome data from 20 NXC (♂:10; ♀:10), the population selection signal was analysed by sliding window analysis. The selected genes were screened by combination with the population differentiation statistic (FST). The sequence diversity statistic (θπ). RNA-seq was performed on the ovarian tissues of NXC (n = 6) and Lohmann laying hens (n = 6) to analyse the differentially expressed genes (DEGs) between the two groups. The functional enrichment analysis of the selected genes and differentially expressed genes was performed.3. There were 742 genes under strong positive selection and 509 differentially expressed genes screened in NXC. Integrated analysis of the genome and transcriptome revealing 26 overlapping genes. The candidate genes for adaptation to a high-altitude environment, as well as for egg production, disease resistance, vision and pigmentation in NXC were preliminarily screened.4. The results provided theoretical guidance for further research on the genetic resource protection and utilisation of NXC.

NR5A1 and NR5A2 regulate follicle development in chicken (Gallus gallus) by altering proliferation, apoptosis, and steroid hormone synthesis of granulosa cells

Chicken ovarian follicle development is regulated by complex and dynamic gene expression. Nuclear receptor 5A1 and 5A2 (NR5A1 and NR5A2, respectively) are key genes that regulate steroid hormone production and gonadal development in mammals; however, studies on follicular development in the chicken ovary are scarce. In this study, we investigated the functions of NR5A1 and NR5A2 on follicle development in chickens. The results showed that the expression of NR5A1 and NR5A2 was significantly higher in small yellow follicles and F5. Furthermore, the expression of NR5A1 and NR5A2 was significantly higher in follicular tissues of peak-laying hens (30 wk) than in follicular tissues of late-laying hens (60 wk), with high expression abundance in granulosa cells (GC). The overexpression of NR5A1 and NR5A2 significantly promoted proliferation and inhibited apoptosis of cultured GC; upregulated STAR, CYP11A1, and CYP19A1 expression and estradiol (E2) and progesterone (P4) synthesis in GC from preovulatory follicles (po-GC); and increased STAR, CYP11A1, and CYP19A1 promoter activities. In addition, follicle-stimulating hormone treatment significantly upregulated NR5A1 and NR5A2 expression in po-GC and significantly promoted FSHR, CYP11A1, and HSD3B1 expression in GC from pre-hierarchical follicles and po-GC. The core promoter region of NR5A1 was identified at the -1,095- to -483-bp and -2,054- to -1,536-bp regions from the translation start site (+1), and the core promoter region of NR5A2 was at -998 to -489 bp. Two single nucleotide polymorphisms (SNP) were identified in the core promoter region of the NR5A1 gene, which differed between high- and low-yielding chicken groups. Our study suggested that NR5A1 and NR5A2 promoted chicken follicle development by promoting GC proliferation and E2 and P4 hormone synthesis and inhibiting apoptosis. Moreover, we identified the promoter core region or functional site that regulates NR5A1 and NR5A2 expression.

Cellular atlas of the human ovary using morphologically guided spatial transcriptomics and single-cell sequencing

The reproductive and endocrine functions of the ovary involve spatially defined interactions among specialized cell populations. Despite the ovary's importance in fertility and endocrine health, functional attributes of ovarian cells are largely uncharacterized. Here, we profiled >18,000 genes in 257 regions from the ovaries of two premenopausal donors to examine the functional units in the ovary. We also generated single-cell RNA sequencing data for 21,198 cells from three additional donors and identified four major cell types and four immune cell subtypes. Custom selection of sampling areas revealed distinct gene activities for oocytes, theca, and granulosa cells. These data contributed panels of oocyte-, theca-, and granulosa-specific genes, thus expanding the knowledge of molecular programs driving follicle development. Serial samples around oocytes and across the cortex and medulla uncovered previously unappreciated variation of hormone and extracellular matrix remodeling activities. This combined spatial and single-cell atlas serves as a resource for future studies of rare cells and pathological states in the ovary.

Identification and functional analysis of ovarian lncRNAs during different egg laying periods in Taihe Black-Bone Chickens

Introduction: Long non-coding RNA (lncRNA) refers to a category of non-coding RNA molecules exceeding 200 nucleotides in length, which exerts a regulatory role in the context of ovarian development. There is a paucity of research examining the involvement of lncRNA in the regulation of ovary development in Taihe Black-Bone Chickens. In order to further investigate the egg laying regulation mechanisms of Taihe Black-Bone Chickens at different periods, transcriptome analysis was conducted on the ovarian tissues at different laying periods. Methods: This study randomly selected ovarian tissues from 12 chickens for RNA-seq. Four chickens were selected for each period, including the early laying period (102 days, Pre), the peak laying period (203 days, Peak), and the late laying period (394 days, Late). Based on our previous study of mRNA expression profiles in the same ovarian tissue, we identified three differentially expressed lncRNAs (DE lncRNAs) at different periods and searched for their cis- and trans-target genes to draw an lncRNA-mRNA network. Results and discussion: In three groups of ovarian tissues, we identified 136 DE lncRNAs, with 8 showing specific expression during the early laying period, 10 showing specific expression during the peak laying period, and 4 showing specific expression during the late laying period. The lncRNA-mRNA network revealed 16 pairs of lncRNA-target genes associated with 7 DE lncRNAs, and these 14 target genes were involved in the regulation of reproductive traits. Furthermore, these reproductive-related target genes were primarily associated with signaling pathways related to follicle and ovary development in Taihe Black-Bone Chickens, including cytokine-cytokine receptor interaction, TGF-beta signaling pathway, tyrosine metabolism, ECM-receptor interaction, focal adhesion, neuroactive ligand-receptor interaction, and cell adhesion molecules (CAMs). This study offers valuable insights for a comprehensive understanding of the influence of lncRNAs on poultry reproductive traits.

Transcriptome analysis of ovarian tissues highlights genes controlling energy homeostasis and oxidative stress as potential drivers of heterosis for egg number and clutch size in crossbred laying hens

Heterosis is the major benefit of crossbreeding and has been exploited in laying hens breeding for a long time. This genetic phenomenon has been linked to various modes of nonadditive gene action. However, the molecular mechanism of heterosis for egg production in laying hens has not been fully elucidated. To fill this research gap, we sequenced mRNAs and lncRNAs of the ovary stroma containing prehierarchical follicles in White Leghorn, Rhode Island Red chickens as well as their reciprocal crossbreds that demonstrated heterosis for egg number and clutch size. We further delineated the modes of mRNAs and lncRNAs expression to identify their potential functions in the observed heterosis. Results showed that dominance was the principal mode of nonadditive expression exhibited by mRNAs and lncRNAs in the prehierarchical follicles of crossbred hens. Specifically, low-parent dominance was the main mode of mRNA expression, while high-parent dominance was the predominant mode of lncRNA expression. Important pathways enriched by genes that showed higher expression in crossbreds compared to either one or both parental lines were cell adhesion molecules, tyrosine and purine metabolism. In contrast, ECM-receptor interaction, focal adhesion, PPAR signaling, and ferroptosis were enriched in genes with lower expression in the crossbred. Protein network interaction identified nonadditively expressed genes including apolipoprotein B (APOB), transferrin, acyl-CoA synthetase medium-chain family member (APOBEC) 3, APOBEC1 complementation factor, and cathepsin S as hub genes. Among these potential hub genes, APOB was the only gene with underdominance expression common to the 2 reciprocal crossbred lines, and has been linked to oxidative stress. LncRNAs with nonadditive expression in the crossbred hens targeted natriuretic peptide receptor 1, epidermal differentiation protein beta, spermatogenesis-associated gene 22, sperm-associated antigen 16, melanocortin 2 receptor, dolichol kinase, glycine amiinotransferase, and prolactin releasing hormone receptor. In conclusion, genes with nonadditive expression in the crossbred may play crucial roles in follicle growth and atresia by improving follicle competence and increasing oxidative stress, respectively. These 2 phenomena could underpin heterosis for egg production in crossbred laying hens.

Whole-Transcriptome Sequencing of Ovary Reveals the ceRNA Regulation Network in Egg Production of Gaoyou Duck

To investigate the regulatory mechanism of the competing endogenous RNAs (ceRNAs) on the egg performance of Gaoyou ducks, full transcriptome sequencing was performed to analyze the ovarian tissues in Gaoyou ducks. The ducks were categorized into high- and low-yield groups based on the individual in-cage egg production records and the hematoxylin-eosin (HE) staining results. The differentially expressed genes (DEGs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) were further processed by GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses. In total, 72 DEmRNAs; 23 DElncRNAs; 4 DEcircRNAs; and 5 signaling pathways, including the ovarian steroidogenesis, PI3K-Akt, hedgehog, tryptophan metabolism, and oocyte meiosis signaling pathways, were significantly enriched. These results suggest that they could be associated with the Gaoyou duck's ovarian function and affect the total egg production or double-yolked egg production. Furthermore, a coregulation network based on the related candidate ceRNAs across the high- and low-yield egg production groups was constructed. Our findings provide new insights into the mechanisms underlying the molecular regulation of related circRNA/lncRNA-miRNA-mRNA in the egg production and double-yolked egg traits of Gaoyou ducks.

The DNA methylation status of the serotonin metabolic pathway associated with reproductive inactivation induced by long-light exposure in Magang geese

BACKGROUND

Domestic geese are seasonal breeders and have the lowest reproductive capacity among all poultry species. Magang geese is a topical short-day breeder, short photoperiod exposure stimulates its reproductive activity while long photoperiod inhibits. To explore epigenetic change that could influence reproductive activity, we performed whole genome bisulfite sequencing and transcriptome sequencing in the hypothalamus at three reproductive stages during long-light exposure in male Magang geese.

RESULTS

A total number of 10,602 differentially methylated regions (DMRs) were identified among three comparison groups. We observed that the vast majority of DMRs were enriched in intron regions. By integrating the BS-sequencing and RNA-seq data, the correlation between methylation changes of CG DMRs and expression changes of their associated genes was significant only for genes containing CG DMRs in their intron. A total of 278 DMR-associated DEGs were obtained among the three stages. KEGG analysis revealed that the DMR-associated DEGs were mainly involved in 11 pathways. Among them, the neuroactive ligand-receptor interaction pathway was significantly enriched in both two comparisons (RA vs.RD and RD vs.RI); the Wnt signaling pathway, apelin signaling pathway, melanogenesis, calcium signaling pathway, focal adhesion, and adherens junction were significantly enriched in the RA vs. RI comparison. In addition, the expression level of two serotonin-metabolic genes was significantly altered during reproductive axis inactivation by the methylation status of their promoter region (TPH2) and intron region (SLC18A2), respectively. These results were confirmed by Bisulfite sequencing PCR (BSP), pyrosequencing, and real-time qPCR, indicating that serotonin metabolic signaling may play a key role in decreasing the reproductive activity of Magang geese induced by long-light exposure. Furthermore, we performed a metabolomics approach to investigate the concentration of neurotransmitters among the three stages, and found that 5-HIAA, the last product of the serotonin metabolic pathway, was significantly decreased in the hypothalamus during RI.

CONCLUSIONS

Our study reveals that the methylation status of the serotonin metabolic pathway in the hypothalamus is associated with reproductive inactivation, and provided new insight into the effect of DNA methylation on the reproductive regulation of the hypothalamus in Magang geese.

Mating modifies the expression of crucial oxidative-reductive transcripts in the pig oviductal sperm reservoir: is the female ensuring sperm survival?

Background

Mating induces large changes in the female genital tract, warranting female homeostasis and immune preparation for pregnancy, including the preservation of crucial oxidative status among its pathways. Being highly susceptible to oxidative stress, sperm survival and preserved function depend on the seminal plasma, a protection that is removed during sperm handling but also after mating when spermatozoa enter the oviduct. Therefore, it is pertinent to consider that the female sperm reservoir takes up this protection, providing a suitable environment for sperm viability. These aspects have not been explored despite the increasing strategies in modulating the female status through diet control and nutritional supplementation.

Aims

To test the hypothesis that mating modifies the expression of crucial oxidative-reductive transcripts across the entire pig female genital tract (cervix to infundibulum) and, particularly in the sperm reservoir at the utero-tubal junction, before ovulation, a period dominated by estrogen stimulation of ovarian as well as of seminal origin.

Methods

The differential expression of estrogen (ER) and progesterone (PR) receptors and of 59 oxidative-reductive transcripts were studied using a species-specific microarray platform, in specific segments of the peri-ovulatory sow reproductive tract in response to mating.

Results

Mating induced changes along the entire tract, with a conspicuous downregulation of both ER and PR and an upregulation of superoxide dismutase 1 (SOD1), glutaredoxin (GLRX3), and peroxiredoxin 1 and 3 (PRDX1, PRDX3), among other NADH Dehydrogenase Ubiquinone Flavoproteins, in the distal uterus segment. These changes perhaps helped prevent oxidative stress in the area adjacent to the sperm reservoir at the utero-tubal junction. Concomitantly, there were a downregulation of catalase (CAT) and NADH dehydrogenase (ubiquinone) oxidoreductases 1 beta subcomplex, subunit 1 (NDUFB1) in the utero-tubal junction alongside an overall downregulation of CAT, SOD1, and PRDX3 in the ampullar and infundibulum segments.

Conclusions

Natural mating is an inducer of changes in the expression of female genes commanding antioxidant enzymes relevant for sperm survival during sperm transport, under predominant estrogen influence through the bloodstream and semen. The findings could contribute to the design of new therapeutics for the female to improve oxidative-reductive balance.

Characterization of ovarian follicles, serum steroid hormone concentration and steroidogenic gene expression profiles in the developing ovarian follicles in White King pigeons

Paired pigeons only lay 2 eggs in a laying period, which is closely related to ovarian follicle development, but this process is not well understood. In this study, 60 pairs of 12-mo-old White King pigeons were selected and serum and follicles were collected at 4 stages of laying interval (LI), including the first (LI1), the third (LI3), the fifth (LI5), and the seventh day (LI7). Morphological results showed that paired pigeons normally had 2 preovulatory follicles and the second-largest follicle (F2) developed from LI3 and had been selected in LI5. Prehierarchical follicles were coupled and hierarchical, which was in accordance with its clutch size. The P₄ concentration increased gradually from LI1 to LI5, reaching a maximum of 30.67 ng/mL in LI5 and decreasing to 27.83 ng/mL in LI7 (P < 0.05). The levels of T in LI1 and LI5 were higher than LI3 and LI7 (P < 0.05), although there was no significant difference in E₂ in LI (P > 0.05), but it stayed at high levels. In the TCs of the largest follicle (F1), HSD3B1 mRNA and HSD17B1 mRNA levels peaked in LI7. The expression pattern of CYP17A1 and CYP19A1 was similar, increasing from LI3 to LI5 and then decreasing. In the TCs of F2, the expressions of HSD3B1 and CYP17A1 had no significant difference between LI5 and LI7 (P > 0.05), while the expression pattern of HSD17B1 and CYP19A1 was the opposite. In TCs of SF1, HSD3B1 mRNA level peaked in LI3 while CYP19A1 mRNA levels peaked in LI7. The expression of CYP17A1 had a minor change (P > 0.05) and the expression pattern of HSD17B1 was similar to F1. It was concluded that the morphological characteristics of follicles during the LI for the first time, including the number and diameter of small follicles (SFs) and hierarchical follicles in pigeon and the concentrations of steroid hormones and expressions of steroidogenic genes in TCs of different follicles could explain the growth and selection of 2 preovulatory follicles. This study facilitates further research into the regulation of ovulation and egg production in pigeons.

DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens

BACKGROUND

The quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation.

RESULTS

The results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group.

CONCLUSIONS

Overall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.

Transcriptomic analysis of ovarian follicles uncovers the crucial genes relevant to follicle selection and preovulatory hierarchy in hens

Follicle selection and preovulatory hierarchy of hen ovaries were important stages of follicle development and crucially determining egg-laying performance. The selected follicles with a higher expression level of follicle-stimulating hormone receptor (FSHR) mRNA that facilitates response to FSH, and rapidly develops into preovulatory follicles with distinctive characteristics of granulosa cells (GCs) proliferation and differentiation. Identification of the key genes involved in these developmental events is helpful for elucidation of the molecular mechanism underlying egg-laying traits in chicken and other domestic fowl. Herein, the comparative transcriptomic analysis of ovarian prehierarchical follicles before selection (BSF), follicles at selection stage (ASF), and hierarchical follicles (HF) were implemented in the Jilin Black chicken (JB) and Lohmann Brown layer (LB) with the divergences in their egg-laying performance by RNA-sequencing. The results showed that nine deferentially expressed genes (DEGs), including STMN4, FABP3, ROBO2, RSPO4, and DMRT1 were revealed between follicles BSF and ASF; and seventeen DEGs, such as SLC6A15, SLITRK3, PRKG2 and TMC3 were mined between ASF and HF. These two group DEGs being co-expressed between BSF and ASF, and between ASF and HF were compared and substantiated in the JB and LB layers, respectively. Furthermore, 10 signaling pathways, such as cAMP signaling, PPAR signaling pathway, AMPK(Adenosine 5'-monophosphate (AMP)-activated protein kinase) pathway, and estrogen signaling pathway were also identified. Moreover, the roles of two representative candidates ROBO2 and PRKG2 genes presented as downregulated mRNA expression pattern in the transcriptomic profiles were further verified in vitro. The results demonstrated that downregulation of ROBO2 or PRKG2 significantly increased the expression levels of FSHR mRNA and protein with the boosted expression of CCND1, STAR, and BCL-2, whereas remarkably inhibited the expression of Caspase-3, consequently, brought about the decrease of GC apoptosis in the ovarian follicles, but increase of GC proliferation and differentiation serving as the hallmarks for follicle selection. It indicated that ROBO2 and PRKG2 may play indispensable roles in follicle selection and preovulatory hierarchy of hen ovaries separately. Our findings provided a comparative transcriptomic evidence for clarifying the molecular mechanism of the follicle development underlying egg-laying traits in chicken.

Transcriptomic analysis of the thyroid and ovarian stroma reveals key pathways and potential candidate genes associated with egg production in ducks

The importance of thyroid-related genes has been repeatedly mentioned in the transcriptome studies of poultry with different laying performance, yet there are few systematic studies to unravel the regulatory mechanisms of the thyroid-ovary axis in the poultry egg production process. In this study, we compared the transcriptome profiles in the thyroid and ovarian stroma between high egg production (GP) and low egg production (DP) ducks, and then revealed the pathways and candidate genes involved in the process. We identified 1,114 and 733 differentially expressed genes (DEGs) in the thyroid and ovarian stroma, separately. The Gene Ontology (GO) analysis showed that a total of 504 and 189 GO terms were identified in the thyroid and ovarian stroma (P < 0.05). Three common GO terms were identified from the top 5 GO terms with the highest significant level in two tissues, including extracellular space, calcium ion binding, and integral component of plasma membrane. The enrichment analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that 15 and 14 KEGG pathways were significantly (P < 0.05) enriched in the thyroid and ovarian stroma, respectively. And, there were 8 common pathways, including neuroactive ligand-receptor interaction, calcium signaling pathway, ECM-receptor interaction, PPAR signaling pathway, melanogenesis, wnt signaling pathway, vascular smooth muscle contraction, and cytokine-cytokine receptor interaction. Notably, the neuroactive ligand-receptor interaction pathway was the most significantly enriched by the DEGs both in the thyroid and ovarian stroma. The interaction among DEGs enriched in the neuroactive ligand-receptor interaction and ECM-receptor interaction suggested that the thyroid may regulate ovarian development by these genes. Through integrated analysis of the protein-protein interaction (PPI) network and KEGG pathway maps, 9 key DEGs (PTH, THBS2, THBS4, CD36, ADIPOQ, ACSL6, PRKAA2, CRH, and PCK1) were identified, which could play crucial roles in the thyroid to regulate ovarian function and then affect egg-laying performance between GP and DP. This study serves as a basis to explore the molecular mechanism of the thyroid affecting ovarian function and egg production in female ducks and may help to identify molecular markers that can be used for duck genetic selection.

Comparative Transcriptome Analysis Provided a New Insight into the Molecular Mechanisms of Epididymis Regulating Semen Volume in Drakes

Semen volume is an important factor in artificial insemination (AI) of ducks. In drakes, seminal plasma that is produced by the epididymis determines the semen volume. However, the mechanism of epididymis regulating semen volume of drakes remains unclear. Therefore, the aim of the present study was to preliminarily reveal the mechanism regulating the semen volume through comparing the epididymal histomorphology and mRNA expression profiles between drakes with high-volume semen (HVS) and low-volume semen (LVS). Phenotypically, drakes in the HVS group produced more sperm than drakes in the LVS group. In addition, compared with the HVS group, the ductal square of ductuli conjugentes (DC) and dutus epididymidis (DE) in epididymis was significantly smaller in the LVS group, and the lumenal diameter and epithelial thickness of DC/DE were significantly shorter in the LVS group. In transcriptional regulation, 72 different expression genes (DEGs) were identified from the epididymis between HVS and LVS groups. Gene Ontology (GO) analysis indicated that the DEGs were mainly related to hormone secretion, neurotransmitter synthesis/transport, transmembrane signal transduction, transmembrane transporter activity, and nervous system development (p < 0.05). Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis showed that the DEGs were significantly enriched in pathways associated with hormone and neurotransmitter transmission (p < 0.05). In addition, further analysis of the top five pathways enriched by KEGG, nine key candidate genes (including SLC18A2, SNAP25, CACNA1B, GABRG2, DRD3, CAMK2A, NR5A1, and STAR) were identified, which could play a crucial role in the formation of semen. These data provide new insights into the molecular mechanism regulating semen volume of drakes and make feasible the breeding of drakes by semen volume.

Transcriptomics and Metabolomics Analysis of the Ovaries of High and Low Egg Production Chickens

Simple Summary

The ovarian tissues of different breeds of hens during egg production were investigated through transcriptomics and metabolomics to provide a more comprehensive understanding of the molecular mechanisms of the ovary during egg production. Four genes involved in egg production were predicted by the transcriptome, including P2RX1, INHBB, VIPR2, and FABP3, and several close metabolites associated with reproduction were identified in the metabolome, including 17α-hydroxyprogesterone, iloprost, spermidine and adenosine. Correlation analysis of specific differential genes and differential metabolites identified important gene-metabolite pairs VIPR2–Spermidine and P2RX1–Spermidine in the reproductive process.

Abstract

Egg production is a pivotal indicator for evaluating the fertility of poultry, and the ovary is an essential organ for egg production and plays an indispensable role in poultry production and reproduction. In order to investigate different aspects of egg production mechanisms in different poultry, in this study we performed a metabolomic analysis of the transcriptomic combination of the ovaries of two chicken breeds, the high-production Ninghai indigenous chickens and the low-production Wuliangshan black-boned chickens, to analyze the biosynthesis and potential key genes and metabolic pathways in the ovaries during egg production. We predicted four genes in the transcriptomic that are associated with egg production, namely P2RX1, INHBB, VIPR2, and FABP3, and identified three important pathways during egg production, “Calcium signaling pathway”, “Neuroactive ligand–receptor interaction” and “Cytokine–cytokine receptor interaction”, respectively. In the metabolomic 149 significantly differential metabolites were identified, 99 in the negative model and 50 in the positive model, of which 17α-hydroxyprogesterone, iloprost, spermidine, and adenosine are important metabolites involved in reproduction. By integrating transcriptomics and metabolomics, the correlation between specific differential genes and differential metabolites identified important gene-metabolite pairs “VIPR2-Spermidine” and “P2RX1-Spermidine” in egg production. In conclusion, these data provide a better understanding of the molecular differences between the ovaries of low- and high-production hens and provide a theoretical basis for further studies on the mechanics of poultry egg production.

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.