Lipid metabolism participates in regulating the functions of granulosa cells (GCs), which is important for follicular development. In this experiment, goose GCs from pre-hierarchical follicles and hierarchical follicles were selected to be the model for studying the putative regulatory role of lipid metabolism in apoptosis and steroidogenesis, through overexpression and interference with fatty acid synthase (FASN). When FASN was overexpressed, the lipid accumulation was increased in hierarchical GCs (hGCs) and it was increased in the two categorized GCs when FASN was interfered. In addition, the apoptosis of the two categorized GCs was increased when FASN was overexpressed, and their progesterone production was decreased when FASN was interfered. The results of qRT-PCR showed that, when FASN was overexpressed, the expression level of CYP11A1 was decreased in pre-hierarchical GCs (phGCs), while the expression levels of SCD1, DGAT2, APOB, and StAR were increased in hGCs. When FASN was interfered, the expression levels of CPT-1, DGAT2, and StAR were decreased whereas the expression level of CYP11A1 was increased in phGCs, and the expression levels of CPT-1, SCD1, and StAR were decreased in hGCs. These results not only identify the different effects of manipulated FASN expression on lipid metabolism of goose phGCs and hGCs but also demonstrate that FASN-mediated lipid metabolism plays an important role in regulating apoptosis and steroidogenesis of in vitro cultured goose GCs.

Long non-coding RNAs (lncRNAs) and mRNAs are temporally expressed during chicken follicle development. However, follicle transcriptome studies in chickens with timepoints relating to changes in luteinizing hormone (LH) levels are rare. In this study, gene expression in Rohman layers was investigated at three distinct stages of the ovulatory cycle: zeitgeber time 0 (ZT0, 9:00 a.m.), zeitgeber time 12 (ZT12, 9:00 p.m.), and zeitgeber time 20 (ZT20, 5:00 a.m.) representing the early, middle, and LH surge stages, respectively, of the ovulatory cycle. Gene expression profiles were explored during follicle development at ZT0, ZT12, and ZT20 using Ribo-Zero RNA sequencing. The three stages were separated into two major stages, including the pre-LH surge and the LH surge stages. A total of 12,479 mRNAs and 7528 lncRNAs were identified among the three stages, and 4531, 523 differentially expressed genes (DEGs) and 2367, 211 differentially expressed lncRNAs (DELs) were identified in the ZT20 vs. ZT12, and ZT12 vs. ZT0, comparisons. Functional enrichment analysis revealed that genes involved in cell proliferation and metabolism processes (lipid-related) were mainly enriched in the ZT0 and ZT12 stages, respectively, and genes related to oxidative stress, steroids regulation, and inflammatory process were enriched in the ZT20 stage. These findings provide the basis for further investigation of the specific genetic and molecular functions of follicle development in chickens.

Ornithine decarboxylase (ODC) plays an indispensable role in the process of polyamine biosynthesis. Polyamines are a pivotal part of living cells and have diverse roles in the regulation of cell proliferation and apoptosis, aging and reproduction. However, to date, there have been no reports about ODC regulating follicular development in goose ovaries. Here, we constructed ODC siRNA and overexpression plasmids and transfected them into goose primary granulosa cells (GCs) to elucidate the effects of ODC interference and overexpression on the polyamine metabolism, hormone levels, cell apoptosis and proliferation of granulosa cells. After interfering with ODC in GCs, the mRNA and protein levels of ODC and the content of putrescine were greatly decreased (P < 0.05). When ODC was overexpressed, ODC mRNA and protein levels and putrescine content were greatly increased (P < 0.05). The polyamine-metabolizing enzyme genes ornithine decarboxylase antizyme 1 (OAZ1) and spermidine / spermine-N1-acetyltransferase (SSAT) were significantly increased, and spermidine synthase (SPDS) was significantly decreased when ODC was downregulated (P < 0.05). OAZ1, SPDS and SSAT were significantly increased when ODC was upregulated (P < 0.05). In addition, after interference with ODC, progesterone (P4) levels in the culture medium of GCs increased greatly (P < 0.05), while the overexpression of ODC caused the P4 level to decrease significantly (P < 0.05). After ODC downregulation, granulosa cell activity was significantly reduced, the apoptosis rate was significantly increased, and the BCL-2 / BAX ratio was downregulated (P < 0.05). Under ODC overexpression, the activity of GCs was notably increased, the apoptosis rate was significantly reduced, and the BCL-2 / BAX protein ratio was upregulated (P < 0.05). Our study successfully induced ODC interference and overexpression in goose ovarian GCs, and ODC regulated mainly putrescine content in GCs with a slight influence on spermidine and spermine. Moreover, ODC participated in the adjustment of P4 levels in the culture medium of GCs, promoted granulosa cell proliferation and inhibited granulosa cell apoptosis.

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.

Although miR-27b-3p has been evidenced to regulate the proliferation, apoptosis, and differentiation of a variety of mammalian cell types, its actions and mechanisms on ovarian cell steroidogenesis remains largely unknown in both mammalian and avian species. In this study, we aimed to determine the expression profiles of miR-27b-3p in granulosa cell layers during goose ovarian follicle development and to reveal its actions on estrogen (E2) secretion of goose granulosa cells as well as the underlying regulatory mechanisms. It was observed that miR-27b-3p was ubiquitously expressed throughout follicle development but exhibited much higher levels in hierarchical- than in prehierarchical follicles. In cultured granulosa cells from the fourth through second largest preovulatory (F4-F2) follicles of goose, up- and downregulation of miR-27b-3p by using its mimic and inhibitor significantly decreased and increased E2 secretion, respectively. Meanwhile, the mRNA levels of STAR and CYP19A1 were significantly reduced while those of CYP11A1 and 3βHSD were elevated in the mimic-transfected granulosa cells. By comparison, downregulation of miR-27b-3p enhanced the mRNA levels of STAR but had no significant effects on those of CYP19A1, CYP11A1, and 3βHSD. Results from bioinformatic prediction and luciferase reporter assay demonstrated that CYP1B1 was a downstream target of miR-27b-3p. Although the siRNA-mediated downregulation of CYP1B1 did not significantly change E2 secretion by goose granulosa cells, it reduced the mRNA levels of STAR and CYP19A1 as well as those of LKB1 and AMPKα, which are involved in the AMPK signaling pathway. Taken together, these data suggest that miR-27b-3p plays an inhibitory role in E2 secretion by goose F4-F2 granulosa cells, at least in part, by targeting CYP1B1 through the AMPK signaling pathway.

Leptin was known as a pivotal regulator for the control of food intake and energy expenditure. However, leptin has also been found to be involved in the regulation of female reproductive system through interactions with pathways in the hypothalamic-hypophyseal axis and direct action at the ovarian level. In the present study, granulosa cells from goose ovarian preovulatory (F1-F3) follicles were cultured with leptin (0, 1, 10 or 100ng/ml). The proliferative and anti-apoptotic actions of leptin in granulosa cells were revealed by CCK-8, BrdU and TUNEL assays. Quantitative real-time PCR and Western blot analyses further indicated that leptin treatment led to increased expression of cyclin D1, cyclin D2, cyclin D3 and bcl-2, and decreased expression of p21 and caspase-3. The effects were involved in the activation of the PI3K/Akt/mTOR signaling pathway, as leptin treatment enhanced the expression of PI3K, Akt1, Akt2, Raptor, mTOR, S6K and p-S6K. Moreover, blockade of the PI3K/Akt/mTOR pathway attenuated the influences of leptin on proliferation and apoptosis of granulosa cells, considering that activated factors by leptin were inhibited in the presence of either 20μM LY294002 (a PI3K inhibitor) or 10μM rapamycin (an mTOR inhibitor). In addition, leptin had a modulatory effect on the expression of its receptor at the transcriptional and translational levels, and blockade of PI3K/Akt/mTOR inhibited both basal and leptin-induced Lepr gene and protein expression. These findings suggest that leptin exerts its proliferative and anti-apoptotic effects on goose granulosa cells through the PI3K/Akt/mTOR signaling pathway via interaction with its receptor.

Oxidative stress is the major culprits responsible for ovarian dysfunction by damaging granulosa cells (GCs). Ferritin heavy chain (FHC) may participate in the regulation of ovarian function by mediating GCs apoptosis. However, the specific regulatory function of FHC in follicular GCs remains unclear. Here, 3-nitropropionic acid (3-NPA) was utilized to establish an oxidative stress model of follicular GCs of Sichuan white geese. To explore the regulatory effects of FHC on oxidative stress and apoptosis of primary GCs in geese by interfering or overexpressing FHC gene. After transfection of siRNA-FHC to GCs for 60 h, the expressions of FHC gene and protein decreased significantly (P < 0.05). After FHC overexpression for 72 h, the expressions of FHC mRNA and protein upregulated considerably (P < 0.05). The activity of GCs was impaired after interfering with FHC and 3-NPA coincubated (P < 0.05). When overexpression of FHC combined with 3-NPA treatment, the activity of GCs was remarkably enhanced (P < 0.05). After interference FHC and 3-NPA treatment, NF-κB and NRF2 gene expression decreased (P < 0.05), the intracellular reactive oxygen species (ROS) level increased greatly (P < 0.05), BCL-2 expression reduced, BAX/BCL-2 ratio intensified (P < 0.05), the mitochondrial membrane potential decreased notably (P < 0.05), and the apoptosis rate of GCs aggravated (P < 0.05). While overexpression of FHC combined with 3-NPA treatment could promote BCL-2 protein expression and reduce BAX/BCL-2 ratio, indicating that FHC regulated the mitochondrial membrane potential and apoptosis of GCs by mediating the expression of BCL-2. Taken together, our research manifested that FHC alleviated the inhibitory effect of 3-NPA on the activity of GCs. FHC knockdown could suppress the expression of NRF2 and NF-κB genes, reduce BCL-2 expression and augment BAX/BCL-2 ratio, contributing to the accumulation of ROS and jeopardizing mitochondrial membrane potential, as well as exacerbating GCs apoptosis.

Bone morphogenetic protein 4 (BMP4) has an important role in regulating cellular proliferation, differentiation and apoptosis. It, however, is still unclear as to the mechanisms by which BMP4 regulates the apoptosis of granulosa cells (GCs) in geese. In the present study, there was cloning of the full-length coding sequence of goose BMP4 gene, which consisted of 1212 nucleotides encoding 403 amino acids. Its deduced amino acid sequence comprised one signal peptide, one TGFβ pro-peptide and one mature peptide domain. Results from conducting the quantitative real-time PCR (qPCR) indicated the relative abundances of BMP4 mRNA in geese GCs increased gradually from the relative abundances in pre-hierarchical follicles that were 4 to 6 mm in diameter to that in the fifth largest (F5) follicle and then relative abundances of BMP4 mRNA decreased with further development as the largest (F1) follicle. Results from use of the TUNEL assay indicated that overexpression of the goose BMP4 gene suppressed GC apoptosis and this was confirmed when relative abundances of the CAD, Caspase-9 and Caspase-3 proteins were determined using western blotting. In addition, overexpression of the BMP4 gene induced phosphorylation of AKT, which was inhibited with use of the PI3K inhibitor, LY294002. Co-transfection of BMP4 and LY294002 resulted in increased relative abundances of Caspase-9 and CAD proteins but had no effect on that of Caspase-3. Taken together, these results suggested that expression of the BMP4 gene resulted in a reduction in Caspase-9 protein leading to inhibition of GC apoptosis via the PI3K/AKT signaling pathway in geese.

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.

Despite their important functions and nearly ubiquitous presence in cells, an understanding of the biology of intracellular lipid droplets (LDs) in goose follicle development remains limited. An integrated study of lipidomic and transcriptomic analyses was performed in a cellular model of stearoyl-CoA desaturase (SCD) function, to determine the effects of intracellular LDs on follicle development in geese.

The circadian clock is reported to play a role in the ovaries in a variety of vertebrate species, including the domestic hen. However, the ovary is an organ that changes daily, and the laying hen maintains a strict follicular hierarchy. The aim of this study was to examine the spatial-temporal expression of several known canonical clock genes in the granulosa and theca layers of six hierarchy follicles. We demonstrated that the granulosa cells (GCs) of the F1-F3 follicles harbored intrinsic oscillatory mechanisms in vivo. In addition, cultured granulosa cells (GCs) from F1 follicles exposed to luteinizing hormone (LH) synchronization displayed Per2 mRNA oscillations, whereas, the less mature GCs (F5 plus F6) displayed no circadian change in Per2 mRNA levels. Cultures containing follicle-stimulating hormone (FSH) combined with LH expressed levels of Per2 mRNA that were 2.5-fold higher than those in cultures with LH or FSH alone. These results show that there is spatial specificity in the localization of clock cells in hen preovulatory follicles. In addition, our results support the hypothesis that gonadotropins provide a cue for the development of the functional cellular clock in immature GCs.

Very low density lipoprotein receptor (VLDLR)-mediated endocytosis of plasma lipoproteins into the ovary is essential for ovarian follicle development. Two splice variants of VLDLR have been identified in several species, yet little is known about their distinctive roles in ovarian developing follicles. In the present study, the full-length cDNAs of two splice isoforms of VLDLR were obtained from geese (Anser cygnoide) ovaries using the RACE method. The longer isoform (TypeI VLDLR) is 3141bp and contains five conserved structural domains, while the other (TypeII VLDLR) lacks 90bp encoding for the O-linked sugar domain. TypeII VLDLR was predominantly expressed in the ovary, with greater amounts of mRNA in theca and granulosa cells from early stages of follicle development but decreased during vitellogenesis. However, there was minimal expression of the TypeI VLDLR gene in theca cells and expression was almost undetectable in granulosa cells throughout follicle development. Yolk VLDL concentrations decreased as stage of development advanced while yolk triglyceride and cholesterol concentrations increased in a follicular size-dependent manner. The significant correlations between transcripts of TypeII VLDLR and yolk lipids supported its important role on yolk lipid deposition. In addition, in vitro experiments suggested that exogenous cholesterol, 25-hydroxycholesterol and mevinolin (a highly potent competitive inhibitor of HMG-CoA) treatment could significantly alter TypeII VLDLR gene expression in granulosa cells from both pre-hierarchical and pre-ovulatory follicles. Collectively, data from the present study indicate that TypeII VLDLR is more important for the transport of plasma lipoproteins into developing follicles than TypeI VLDLR, and provide new evidence about the influence of steroids in modulating VLDLR gene expression in ovarian cells.

The transcription factors GATA-4 and GATA-6, members of the GATA family, play an important role in ovarian cell proliferation, differentiation and apoptosis. In this study, the full-length coding sequences of goose GATA-4 and GATA-6 were cloned and characterized. GATA-4 and GATA-6 consist of 1236 and 1104 nucleotides encoding proteins with 411 and 367 amino acids, respectively. The deduced amino acid sequences of both proteins include two adjacent zinc finger domains with the distinctive form (CVNC-X17-CNAC)-X29-(CANC-X17-CNAC) and share 84.76% identity within this domain. In silico prediction together with matching of the high affinity RRXS(T)Y motif revealed that the GATA-4 protein might be phosphorylated predominantly at S(233), but no phosphorylation site was found in the GATA-6 protein. Real-time quantitative PCR analysis showed that GATA-4 and GATA-6 mRNAs were co-expressed in goose follicles, moderately expressed in granulosa cells and weakly expressed in theca cells. The expression level of GATA-4 mRNA in healthy follicles was significantly higher than in atretic follicles or postovulatory follicles (P<0.01), and the expression level of GATA-6 mRNA in healthy follicles was significantly lower than in atretic follicles or postovulatory follicles (P<0.01). The expression level of GATA-4 mRNA in granulosa cells was downregulated during follicle development; the peak of expression occurred in the 8-10 mm follicles, and the lowest expression was in the F1 follicles. GATA-6 was upregulated and reached its peak expression in the F1 follicles. These results indicate that the molecular structural differences in goose GATA-4 and GATA-6 may be related to their different roles during follicle development.

Anti-Müllerian hormone (AMH) is recognized as a reliable marker of ovarian reserve. However, the regulatory mechanism of goose AMH gene remains poorly understood. In the present study, both the full-length coding sequence (CDS) and promoter sequence of goose AMH have been cloned. Its CDS consisted of 2013 nucleotides encoding 670 amino acids and the amino acid sequence contained two structural domain: AMH-N and transforming growth factor beta (TGF-β) domain. The obtained promoter sequence spanned from the -2386 bp to its transcription start site (ATG). Core promoter regions and regulatory elements were identified as well as transcription factors were predicted in its promoter sequence. The luciferase activity was the highest spanning from the -331 to -1 bp by constructing deletion promoter reporter vectors. In CHO cells, the luciferase activity significantly increased by co-expression of AMH and GATA binding protein 4 (GATA-4), while that significantly decreased by mutating the binding sites of GATA-4 located in the -778 and -1477 bp. Results from quantitative real-time polymerase chain reaction (qPCR) indicated that levels of AMH mRNA in geese granulosa layers decreased gradually with the increasing follicular diameter. Taken together, it could be concluded that the transcriptional activity of AMH was activated by GATA-4 to inhibit the development of small follicles in goose.

Stearoyl-CoA desaturase (SCD) is known to be an important rate-limiting enzyme in the production of monounsaturated fatty acids (MUFAs). However, the role of this enzyme in goose follicular development is poorly understood. To investigate the metabolic mechanism of SCD during goose follicular development, we observed its expression patterns in vivo and in vitro using quantitative reverse-transcription (qRT)-PCR. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) was used to determine a cellular model of SCD function in granulosa cells (GCs) via SCD overexpression and knockdown. qRT-PCR analysis showed that SCD was abundantly expressed in the GC layer, and was upregulated in preovulatory follicles. Peak expression was found in F1 and prehierarchal follicles with diameters of 4-6 mm and 8-10 mm, respectively. We further found that mRNA expression and corresponding enzyme activity occur in a time-dependent oscillation pattern in vitro, beginning on the first day of GC culture. By LC-MS/MS, we identified numerous changes in metabolite activation and developed an overview of multiple metabolic pathways, 10 of which were associated with lipid metabolism and enriched in both the overexpressed and knockdown groups. Finally, we confirmed cholesterol and pantothenol or pantothenate as potential metabolite biomarkers to study SCD-related lipid metabolism in goose GCs.