To date, research on poultry egg production performance has only been conducted within inter or intra-breed groups, while those combining both inter- and intra-breed groups are lacking. Egg production performance is known to differ markedly between Sichuan white goose (Anser cygnoides) and Landes goose (Anser anser). In order to understand the mechanism of egg production performance in geese, we undertook this study. Here, 18 ovarian stromal samples from both Sichuan white goose and Landes goose at the age of 145 days (3 individuals before egg production initiation for each breed) and 730 days (3 high- and low egg production individuals during non-laying periods for each breed) were collected to reveal the genome-wide expression profiles of ovarian mRNAs and lncRNAs between these two geese breeds at different physiological stages. Briefly, 58, 347, 797, 777, and 881 differentially expressed genes (DEGs) and 56, 24, 154, 105, and 224 differentially expressed long non-coding RNAs (DElncRNAs) were found in LLD vs. HLD (low egg production Landes goose vs. high egg production Landes goose), LSC vs. HSC (low egg production Sichuan White goose vs. high egg production Sichuan white goose), YLD vs. YSC (young Landes goose vs. young Sichuan white goose), HLD vs. HSC (high egg production Landes goose vs. high egg production Sichuan white goose), and LLD vs. LSC (low egg production Landes goose vs. low egg production Sichuan white goose) groups, respectively. Functional enrichment analysis of these DEGs and DElncRNAs suggest that the “neuroactive ligand–receptor interaction pathway” is crucial for egg production, and particularly, members of the 5-hydroxytryptamine receptor (HTR) family affect egg production by regulating ovarian metabolic function. Furthermore, the big differences in the secondary structures among HTR1F and HTR1B, HTR2B, and HTR7 may lead to their different expression patterns in goose ovaries of both inter- and intra-breed groups. These results provide novel insights into the mechanisms regulating poultry egg production performance.

MiRNAs are critical for steroidogenesis in granulosa cells (GCs) during ovarian follicular development. We have previously shown that miR-202-5p displays a stage-dependent expression pattern in GCs from goose follicles of different sizes, suggesting that this miRNA could be involved in the regulation of the functions of goose GCs; therefore, in this study, the effects of miR-202-5p on lipid metabolism and steroidogenesis in goose hierarchical follicular GCs (hGCs), as well as its mechanisms of action, were evaluated. Oil Red O staining and analyses of intracellular cholesterol and triglyceride contents showed that the overexpression of miR-202-5p significantly inhibited lipid deposition in hGCs; additionally, miR-202-5p significantly inhibited progesterone secretion in hGCs. A bioinformatics analysis and luciferase reporter assay indicated that Acyl-CoA synthetase long-chain family member 3 (ACSL3), which activates long-chain fatty acids for the synthesis of cellular lipids, is a potential target of miR-202-5p. ACSL3 silencing inhibited lipid deposition and estrogen secretion in hGCs. These data suggest that miR-202-5p exerts inhibitory effects on lipid deposition and steroidogenesis in goose hGCs by targeting the ACSL3 gene.

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8–10 mm) and those of hierarchical follicles (hGCs, F2–F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3′UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese’ egg production performance.

Currently, FRS and CRS are the two predominant dryland rearing systems in the goose industry. However, the effects of these two systems on goose growth performance and health, as well as the underlying mechanisms, have not been fully clarified. Thus, this study aimed to compare growth performance and immune status, as well as investigate the genome-wide transcriptomic profiles of spleen in geese, between CRS and FRS at 270 d of age. Phenotypically, the body weight and body size traits were higher in geese under FRS, while the weight and organ index of spleen were higher in geese under CRS (p < 0.05). Noticeably, the bursa of Fabricius of geese under FRS was degenerated, while that under CRS was retained. At the serum level, the immune globulin-G (IgG) and interleukin-6 (IL-6) levels were higher in geese under CRS (p < 0.05). At the transcriptomic level, we identified 251 differentially expressed genes (DEGs) in the spleen between CRS and FRS, which were mainly enriched in scavenger receptor activity, inflammatory response, immune response, neuroactive ligand–receptor interaction, phenylalanine metabolism, ECM receptor interaction, calcium signaling pathway, phenylalanine, tyrosine, and tryptophan biosynthesis, regulation of actin cytoskeleton, and MAPK signaling pathways. Furthermore, through protein–protein interaction (PPI) network analysis, ten candidate genes were identified, namely, VEGFA, FGF2, NGF, GPC1, NKX2-5, FGFR1, FGF1, MEIS1, CD36, and PAH. Further analysis demonstrated that geese in CRS could improve their immune ability through the “phenylalanine metabolism” pathway. Our results revealed that the FRS improved growth performance, whereas the CRS improved goose immune function by increasing levels of IL-6 and IgG in serum. Moreover, the phenylalanine metabolism pathway could exert positive effects on immune function of geese under CRS. These results can provide reliable references for understanding how floor and cage rearing systems affect goose growth performance and immune capacity.

The lower egg production of geese (20~60 eggs per year) compared with chicken and duck limits the development of the industry, while the yolk weight and fatty liver susceptibility of geese was higher than that of other poultry. Therefore, the relationship between lipid metabolism and the laying performance of geese remains to be explored. Phenotypically, we observed that the liver fat content of the high-, low-, and no-egg production groups decreased in turn, while the abdominal fat weight increased in turn. For transcriptional regulation, the KEGG pathways related to lipid metabolism were enriched in all pairwise comparisons of abdominal fat and liver through functional analysis. However, some KEGG pathways related to inflammation and the circadian rhythm pathway were enriched by DEGs only in abdominal fat and the liver, respectively. The DEGs in ovarian stroma among different groups enriched some KEGG pathways related to ovarian steroidogenesis and cell adhesion. Our research reveals that lipid metabolism regulated by the circadian rhythm of the liver may directly or indirectly affect ovarian function through the inflammation and hormone secretion of abdominal fat. These results offer new insights into the regulation mechanisms of goose reproductive traits.

(1) Background: The role of estrogen receptor gene 1 (ESR1) in female reproduction and lipid metabolism has been extensively investigated. However, its contribution to lipid metabolism during the development of poultry follicles remains unclear. (2) Methods: This study aimed to explore the function of ESR1 via overexpressing (ESR1^ov) and interfering (ESR1^si) with its expression in pre-hierarchical granulosa cells (phGCs) and hierarchical granulosa cells (poGCs). (3) Results: We successfully cloned and obtained an 1866 bp segment of the full-length CDS region of the Sichuan white goose ESR1 gene. In phGCs of the ESR1^ov and ESR1^si groups, there were no significant changes compared to the control group. However, in poGCs, the ESR1^ov group exhibited decreased lipid deposition, triglycerides, and cholesterol compared to the control group, while the ESR1^si group showed increased lipid deposition, triglycerides, and cholesterol. The expression of APOB and PPARα was significantly reduced in the ESR1^ov group compared to the ESR1^ov-NC group. Moreover, significant changes in the expression of ACCα, DGAT1, SCD, CPT1, and ATGL were observed between the ESR1^si and ESR1^si-NC group. (4) Conclusions: These findings shed light on the function and molecular mechanism of ESR1 in lipid metabolism in goose poGCs, providing a better understanding of the physiological process of goose follicular development.