Next Article in Journal
Diverse Sphingolipid Species Harbor Different Effects on Ire1 Clustering
Next Article in Special Issue
Hepato-Protective Effects of Delta-Tocotrienol and Alpha-Tocopherol in Patients with Non-Alcoholic Fatty Liver Disease: Regulation of Circulating MicroRNA Expression
Previous Article in Journal
Discovery of Novel HIV Protease Inhibitors Using Modern Computational Techniques
Previous Article in Special Issue
Micro-Players of Great Significance—Host microRNA Signature in Viral Infections in Humans and Animals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LncGSAR Controls Ovarian Granulosa Cell Steroidogenesis via Sponging MiR-125b to Activate SCAP/SREBP Pathway

1
Laboratory of Small Ruminant Genetics, Breeding and Reproduction, College of Animal Science and Technology, Hebei Agricultural University, Baoding 071000, China
2
College of Life Science, Hebei Agricultural University, Baoding 071000, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(20), 12132; https://doi.org/10.3390/ijms232012132
Submission received: 14 September 2022 / Revised: 30 September 2022 / Accepted: 30 September 2022 / Published: 12 October 2022

Abstract

:
Long non-coding RNAs (lncRNAs) have been shown to play important roles in livestock fecundity, and many lncRNAs that affect follicular development and reproductive diseases have been identified in the ovary. However, only a few of them have been functionally annotated and mechanistically validated. In this study, we identified a new lncRNA (lncGSAR) and investigated its effects on the proliferation and steroidogenesis of ovine granulosa cells (GCs). High concentrations of glucose (add 33.6 mM glucose) caused high expression of lncGSAR in GCs by regulating its stability, and lncGSAR overexpression promoted GCs proliferation, estrogen secretion, and inhibited progesterone secretion, whereas interference with lncGASR had the opposite effect. Next, we found that the RNA molecules of lncGSAR act on MiR-125b as competitive endogenous RNA (ceRNA), and SREBP-cleavage-activating protein (SCAP) was verified as a target of MiR-125b. LncGASR overexpression increased the expression of SCAP, SREBP, and steroid hormone-related proteins, which can be attenuated by MiR-125b. Our results demonstrated that lncGSAR can act as a ceRNA to activate SCAP/SREBP signaling by sponging MiR-125b to regulate steroid hormone secretion in GCs. These findings provide new insights into the mechanisms of nutrient-regulated follicle development in ewes.

Graphical Abstract

1. Introduction

High reproductive efficiency (i.e., litter size) is the core of the sheep breeding industry, and the productivity of a breeding ewe determines its commercial value. Normal proliferation and steroidogenesis of follicle granulosa cells (GCs) are crucial for ovarian follicular development, oocyte maturation, and subsequent embryonic development [1,2]. Therefore, the discovery of genetic regulatory factors involved in GCs function is of great significance for improving ovulation rates and ewe fertility.
Normal glucose metabolism in GCs is essential for oocyte development and maturation as well as the protection of GCs development [3,4]. However, chronically high glucose levels have deleterious effects on the structure and function of the ovary, especially oocytes and GCs during folliculogenesis [5]. These complex cellular and development processes depend on the precise spatiotemporal expression of regulatory factors. Our previous studies found that glucose dose had a direct effect on GCs proliferation, apoptosis, and steroid hormone secretion, and glycolytic metabolites were positively correlated with steroid hormone secretions [6,7]. This confirms the connection between GCs steroidogenesis and glucose availability, but the molecular mechanisms underlying these changes require further study.
Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs that are longer than 200 nucleotides in length and no protein-coding potential, which are involved in numerous important biological processes [8]. In reproductive activities, lncRNAs played diverse roles in the regulation of follicular development [9], oocyte maturation [10], GCs differentiation [11], and reproductive diseases [12,13]. For example, NEAT1, a lncRNA, was demonstrated to facilitate ovarian GCs proliferation and inhibit apoptosis by regulating the IGF1 [14]. At the post-transcriptional level, lncRNA can serve as efficient microRNA (miRNA) sponges—termed competing endogenous RNAs (ceRNAs)—that interact with miRNA to regulate gene expression [15]. For example, FDNCR sponges miR-543-3p in GCs and prevents miR-543-3p from binding to the DCN 3′ UTR, resulting in DCN transactivation and TGF-b pathway inhibition and promotion of GCs apoptosis in Hu sheep [16]. In addition, lncRNA NORHA acts as a “sponge” that directly binds to the miR-183-96-182 cluster to induce apoptosis in porcine ovarian GCs [17]. Although functions of these lncRNAs have been partially characterized, most of their roles for steroid hormone secretion are still poorly understood.
In this study, to explore the role of lncRNAs in regulating glucose-stimulated ovarian GCs steroidogenesis, ovarian GCs cultured with different glucose doses (optimum glucose concentration (8.4 mM) and high glucose concentration (33.6 mM)) [6,7] were used for RNA sequencing (RNA-seq). Based on this result, we investigated the function of a novel lncRNA named granulosa cell steroidogenesis-associated RNA (GSAR). Moreover, the function and underlying regulatory mechanisms of lncGSAR were further explored, and our results elucidated that mechanism of lncGSAR served as a ceRNA regulating ovarian GCs proliferation and steroid hormone secretion. These results provide some clues to the lncRNAs mechanisms of regulation in nutrient-stimulated follicle development and ovulation.

2. Results

2.1. LncRNA-miRNA-mRNA Interaction Network

There are 461 (281 upregulated and 180 downregulated) differentially expressed lncRNAs (Figure 1A, Table S3), 94 (47 upregulated and 47 downregulated) differentially expressed miRNAs (Figure 1B, Table S4), and 796 (379 upregulated and 417 downregulated) differentially expressed genes (Figure 1C, Table S5) that were identified with a p-value < 0.05 as the cut-off. To determine the possible functional significance of observed changes in lncRNA levels in the optimum and high-glucose-induced GCs groups, a gene ontology (GO) term enrichment analysis was performed. We found the differentially expressed lncRNAs were found to be similar and significantly associated with steroid hormone stimulus in biological process term enrichment (Figure S1). Interestingly, the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of lncRNA target genes also found the ovarian steroidogenesis pathway and steroid hormone biosynthesis pathway were both enriched in comparison group (Figure S2). These findings suggest that steroidogenesis may play critical roles in glucose-induced function of GCs.
Protein-coding RNAs (mRNAs) and lncRNAs, which share the common miRNA response elements (MREs), can both compete for binding to miRNAs and regulate each other. In this study, we constructed a putative lncRNA–miRNA–mRNA crosstalk network involved in GCs steroidogenesis. A novel differentially expressed lncRNA (the intergenic lncRNA TCONS_00219351 is located on chromosome chr26: 25174136-25228108, and its RNA sequence is shown in Figure S3) was served as a candidate. The expression of lncRNA TCONS_00219351 was significantly up-regulated in the high-glucose group compared with the optimal-glucose concentration group, and the analysis of qRT-PCR also verified this result (Figure 1D). Thus, we suspected that lncRNA TCONS_00219351 plays an important role in steroid hormone synthesis of GCs and hereinafter refer to this lncRNA as lncGSAR for convenience. We further investigated the subcellular localization of lncGSAR, and the RT-PCR result confirmed that it is an RNA molecule present in the cytoplasm and nucleus (Figure 1E). The lncGSAR-binding miRNA were predicted as candidates using RNAhybrid and miRanda software. A ceRNAs network (lncGSAR–miRNA–mRNA) with 1 miRNA and 10 mRNAs was constructed in GCs (Figure 1F, Table S6). Notably, SREBP-cleavage-activating protein (SCAP) gene from lncRNA–miRNA–gene network was a critical sensor of glutamine, glucose, and sterol levels [18]. SCAP regulates sterol production by activating the sterol regulatory element-binding proteins (SREBPs) signaling pathway, which plays an important role in the subsequent steroid hormone synthesis [19].

2.2. High Glucose Regulates lncGSAR Expression by Affecting Its Stability

To validate the effect of different glucose on the expression of lncGSAR, we treated g GCs with high-concentration glucose with or without actinomycin D (ACTD), which is an RNA transcription inhibitor. As shown in Figure 2A, the lncGSAR level was increased after high-concentration glucose treatment; however, the enhanced lncGSAR fails to be rescued by ACTD, which suggests high-concentration glucose regulated lncGSAR expression dependent on enhancing its stability. To further explore the mechanism of glucose enhances lncGSAR stability, we cloned the full-length lncGSAR sequence after luciferase element of pMIR vector to obtain pMIR-luciferase lncGSAR and cloned the lncGSAR promoter before the luciferase element of Pgl3 vector to obtain Pgl3-lncGSAR promoter luciferase plasmids, respectively (Figure 2B). It is no surprise that the luciferase activity of Pgl3-lncGSAR promoter luciferase cannot be changed, yet the activity of pMIR-luciferase lncGSAR was increased under high-level glucose treatment (Figure 2C,D). Thus, these data indicate expression of lncGSAR was increased through reinforcing its stability under high-glucose stress.

2.3. LncGSAR Promotes Steroid Hormone Secretion and Proliferation of GCs

In order to assess the function of lncGSAR in GCs, lncGSAR overexpression vector (pcDNA3.1-lncGSAR) or lncGSAR small interfering RNA (si-lncGSAR) were constructed and transfected into GCs. The cell counting kit-8 (CCK-8) assay showed that overexpression of lncGSAR significantly increased GCs viability (Figure 3A). The 5′-bromo-2′-deoxyuridine (BrdU) detection also demonstrated that the proliferation rate of lncGSAR overexpression cells was significantly increased compared with that of the control (Figure 3C). Conversely, the opposite result was observed by lncGSAR interference (Figure 3B,C), indicating that lncGSAR can facilitate GCs proliferation.
To further investigate the role of lncGSAR in secretion of steroid hormones, we measured the concentrations of estradiol (E2) and progesterone (P4) in GCs by enzyme-linked immunosorbent assay (ELISA). The result revealed that overexpression of lncGSAR significantly increased concentrations of E2 (Figure 3D) but decreased concentrations of P4 (Figure 3F). In addition, we found that lncGSAR overexpression significantly upregulated expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs (Figure 3H). Western blot analysis also showed that overexpression of lncGSAR significantly improved the expression levels of steroidogenesis-related proteins and sterol regulatory element proteins (Figure 3J). On the contrary, lncGSAR interference inhibited E2 secretion (Figure 3E) but promoted P4 secretion (Figure 3G) as well as significantly downregulated mRNA and protein expression levels of steroidogenesis-related and sterol regulation-related genes (Figure 3I,J). These results suggest that lncGSAR might play a crucial role in the regulation of steroidogenesis in GCs.

2.4. LncGSAR Acts as a Molecular Sponge for MiR-125b

MiR-125b is the core component of lncRNA–miRNA–gene network and has been reported to play a critical role in reproductive hormone biosynthesis and follicle development [20]. In addition, MiR-125b regulated GCs apoptosis in the yak ovary by targeting BMPR1B [21]. In the present study, the lncRNA–miRNA–gene network suggests both the lncGSAR and MiR-125b exhibits the deeply potential to interact with each other (Figure 1F). Furthermore, the RNA-seq showed that the expression of MiR-125b was significantly downregulated in the high-glucose group, and the analysis of qRT-PCR also verified this result (Figure 4A). There was a significant negatively correlation between the MiR-125b expression and lncGSAR expression (Figure 4B).
To explore the regulatory relationship between lncGSAR and MiR-125b, granulosa cells were transfected with different doses of MiR-125b mimics (0, 300, 500, and 800 ng) for 48 h. We found that the expression levels of lncGSAR decreased with increasing doses of MiR-125b mimics (Figure 4D), while expression levels of MiR-125b gradually increased (Figure 4C). Next, we added MiR-125b inhibitor to the MiR-125b mimic group and MiR-NC group, respectively, and found that MiR-125b inhibitor could rescue MiR-125b-induced downregulation of lncGSAR expression level (Figure 4E). It is worth noting that addition of MiR-125b mimic reversed the upregulation of lncGSAR expression by high glucose so that the mRNA expression level of lncGSAR had no significant difference between the high-glucose group and the optimal-concentration group (p < 0.05) (Figure 4F). To further investigate the MiR-125b directly binding to lncGSAR, dual-luciferase reporter constructs containing the miRNA response element (MRE; wildtype (WT)) and mutant (MT) plasmid were co-transfected with MiR-125b mimics into GCs (Figure 4G). Luciferase activity of the lncGSAR-WT was dramatically decreased; nevertheless, MiR-125b failed to inhibit the lncGSAR-mut luciferase activity. (Figure 4H). Meanwhile, the RNA–RNA binding assay showed that MiR-125b directly contacts with lncGSAR and not the control RNA GAPDH (Figure 4I).

2.5. MiR-125b Suppresses GCs Proliferation and Steroidogenesis

To further explore the critical functions of MiR-125b in GCs, we transfected GCs with miRNA mimics or miRNA-negative control (MiR-NC) and then monitored the proliferation status of cells using CCK-8 assay, BrdU staining, and flow cytometry analysis. CCK-8 and BrdU staining demonstrated that the proliferation rate of MiR-125b-transfected GCs was significantly reduced compared with that of the MiR-NC-transfected cells (Figure 5A,B). Flow cytometry analysis of the cell cycle revealed that GCs transfected with the MiR-125b mimic could elevate the ratio of cells that progressed to the G1 phase and reduced the ratio of cells that progressed to the S phase (Figure 5C–H). After ELISA, we found MiR-125b overexpression significantly inhibited concentrations of E2 but significantly promoted production of P4 (Figure 5I,J). The expression levels of steroidogenesis- related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs were downregulated in MiR-125b overexpression GCs compared to MiR-NC transfected GCs (Figure 5K). In addition, overexpression of MiR-125b reduced the expression level of steroidogenesis-related proteins and sterol regulatory element proteins (Figure 5J). Together, these results suggest that MiR-125b inhibits GC proliferation and steroidogenesis.

2.6. MiR-125b Regulates GCs Proliferation and Steroidogenesis by Targeting SCAP-SREBP Axis

Based on the above ceRNA (lncGSAR–MiR-125b–mRNA) network, several genes (e.g., SCAP and GPRC5A) were predicted to be MiR-125b target genes. SCAP, a transmembrane structural protein located on the endoplasmic reticulum and functioning as a sterol sensor, attracts our attention [22]. Thus, a SCAP was selected as a candidate target of the MiR-125b for further study. In the present study, the RNA-seq and qRT-PCR both showed that the expression of SCAP was significantly upregulated in the high-glucose group (Figure 6A). There was a significant negatively correlation between the MiR-125b expression and SCAP expression (Figure 6B). Meanwhile, the proliferation rate and E2 concentrations of MiR-125b transfected GCs were significantly downregulated (Figure 6C,D), while the P4 concentration was upregulated (Figure 6E). However, overexpression of SCAP alleviated the reduced proliferation rate and concentrations of E2 or increased concentrations of P4 driven by MiR-125b (Figure 6C–E). Overexpression of SREBP also reversed MiR-125b-driven inhibition of proliferation rate and E2 concentration or promotion of P4 concentration (Figure 6F–H). In addition, we found MiR-125b overexpression significantly decreased the protein expression level of SCAP and simultaneously downregulated SREBP protein, and the SCAP and SREBP overexpression efficiency were shown in Figure 6I,J.
Next, luciferase reporter plasmids containing the SCAP-WT MRE motif or the mutated versions were constructed and co-transfected with MiR-125b mimics into GCs (Figure 6K). MiR-125b significantly reduced luciferase activity of the reporter containing the MRE-WT motif, while SCAP-MT luciferase reporter activity was unchanged under the MiR-125b overexpression (Figure 6L).

2.7. SCAP Is Involved in GCs Proliferation and Steroidogenesis via Activate SREBP

We investigated the potential functions of SCAP in GCs. Here, we transfected SCAP overexpression vector (pcDNA3.1-SCAP) or SCAP small interfering RNA (si-lncGSAR) into GCs. The CCK-8 assay and BrdU staining showed that overexpression of SCAP significantly promoted GC proliferation rate compared with that of the control (Figure 7A,C). Conversely, proliferation rate was significantly inhibited after SCAP KD (Figure 7B,C), indicating that SCAP promotes GCs proliferation.
To assess whether SCAP affects steroid hormone levels in GCs, we measured the concentrations of E2 and P4 in GCs after transfection pcDNA3.1-SCAP or vector control for 48 h by ELISA. The result revealed that overexpression of SCAP significantly increased concentrations of E2 (Figure 7D) but decreased concentrations of P4 (Figure 7F). In addition, we found that SCAP overexpression significantly upregulated expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs (Figure 7H). Western blot analysis also showed that overexpression of lncGSAR significantly improved the expression levels of steroidogenesis-related proteins and sterol regulatory element-binding proteins (Figure 7J). On the contrary, SCAP interference inhibited E2 secretion (Figure 7E) but promoted P4 secretion (Figure 7G) as well as significantly downregulated mRNA and protein expression levels of steroidogenesis-related and sterol regulation-related genes (Figure 7I,J). Together, these results suggest that SCAP is involved in GCs proliferation and steroidogenesis via activate SREBP.

2.8. LncRNA GSAR Promotes GCs Proliferation and Steroidogenesis via Activating the SCAP-SREBP Pathway by Acting as an MiR-125b Sponge

The established binding mechanism between MiR-125b/lncGSAR and MiR-125b/SCAP inspired us to further determine how lncGSAR regulates the direct targeting of MiR-125b to SCAP and whether lncGSAR regulated steroidogenesis-related proteins’ and sterol regulatory element-binding proteins’ expression is dependent on SCAP. Therefore, we firstly performed Western blot and dual-luciferase assays on GCs after co-transfection of lncGSAR, MiR-125b, and specific siRNA of SCAP. We found that overexpression of lncGSAR could restore the inhibition of steroidogenesis-related proteins’ and sterol regulatory element-binding proteins’ expression exerted by MiR-125b. Meanwhile, lncGSAR regulates these proteins’ expression, which is dependent on its regulation of SCAP (Figure 8A). Additionally, a dual-luciferase reporter analysis similarly demonstrated that lncGSAR blocks the direct binding of MiR-125b and its target SCAP (Figure 8B).
To further explore whether lncGSAR affects the binding of MiR-125b to its target SCAP, we labeled MiR-125b and lncGSAR with biotin and digoxigenin, respectively, as well as performed RNA-RNA binding assays. We found that biotin-labeled MiR-125b was associated with SCAP but not unlabeled MiR-125b. Furthermore, digoxigenin-labeled-sense-lncGSAR disrupted the binding between MiR-125b and SCAP but not for the digoxigenin-labeled antisense lncGSAR (Figure 8C). Taken together, all the results support the idea that lncGSAR acts as an MiR-125b sponge to protect the SCAP from the attack of MiR-125b, thereby promoting cell proliferation and E2 secretion and inhibiting P4 secretion.

3. Discussion

Ovarian follicular development is a highly coordinated and nutrition-sensitive process [23], and glucose is an important energy substrate for metabolism in the follicles of many species [24,25,26,27]. Previous studies confirm that short-term dietary supplementation of ewes during the luteal phase can increase fertility due to elevated glucose and steroid hormone concentrations in follicles [19,28,29], whereas glucose has deleterious effects on ovary structure and function at higher concentrations [5,30,31]. Our previous study also demonstrated that glucose dose had a significant effect on GCs function, with 8.4 mM representing the optimal glucose concentration for GCs to secrete steroid hormones in vitro and a higher concentration of glucose (33.6 mM) inhibiting GCs proliferation and steroidogenesis [6,7]. Therefore, there is a potential regulatory relationship between glucose metabolism and steroid hormone synthesis, and further investigation is required to validate the molecular mechanism of glucose-induced GCs differentiation and steroid hormone secretion. In the current study, we identified a candidate lncRNA (lncGSAR) involved in ovine follicular steroidogenesis and firstly reported the expression, function, and mechanism of lncRNAs in glucose-induced GCs steroidogenesis. Therefore, this highlighted the importance of lncRNAs in regulating ewe fertility, which may provide new perspectives on the molecular mechanisms by which nutrition regulates reproduction in ewes.
With the development of next-generation high-throughput sequencing, multiple lncRNAs have been identified and been shown to play different roles in regulating fecundity. However, studies on lncRNAs related to the regulation of follicle development mostly come from human reproductive diseases (such as polycystic ovary syndrome [13,32] and ovarian cancer [33]), and most of them focus on the effects of lncRNAs on GCs proliferation, apoptosis, and differentiation [34]. Our study demonstrated that lncGSAR was mainly expressed in high-glucose group GCs and controlled follicular development by regulating GCs steroidogenesis. This suggests that lncGSAR is involved in sheep nutrient (glucose) stimulated follicle development and regulates GCs hormone secretion. Studies found that lncRNAs containing MREs act as molecular sponges to effectively inhibit miRNA function [35,36]. For instance, lncRNA MLK7-AS1 regulates proliferation, metastasis, and EMT process in ovarian cancer cells by sponging miR-375 [37]. LncRNA OIP5-AS1 sponging miR-34a to promote ovarian carcinoma cell invasion and migration [38]. Similarly, using bioinformatics, dual-luciferase reporter assay, and RNA-RNA binding assay, we further confirmed that lncGSAR can adsorb MiR-125b in GCs to inhibit follicular function. These studies may provide insight into the mechanisms of lncRNA regulates fecundity in sheep.
As we all know, miRNAs directly or indirectly regulate follicular function [39] and hormone secretion [40]. MiR-125b is an important miRNA, and many studies have focused on its biological functions in cell proliferation and differentiation [41] as well as cancer cells [42]. Emerging evidence suggests that in both mouse and human GCs, androgens suppress follicular atresia by enhancing MiR-125b expression, which then targets proapoptotic proteins [43]. On this basis, our study revealed that MiR-125b is involved in the regulation of GCs hormone secretion function and further identified the sterol regulation-related gene SCAP as a novel target of MiR-125b. SCAP is a protein with a sterol-sensing domain (SSD) and seven WD domains. In the presence of cholesterol, this protein binds to sterol regulatory element-binding protein (SREBP) and mediates their transport from the endoplasmic reticulum to the Golgi apparatus. SREBP is then proteolytically cleaved and regulates sterol biosynthesis [44,45,46]. In present study, we also found that the SCAP/SREBP signaling pathway can regulate steroid hormone secretion in GCs by enhancing the expression of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) genes, which is consistent with the role of lncGSAR and inverse to that MiR-125b. Meanwhile, our results showed that overexpression of MiR-125b reduced the protein levels of SCAP and SREBP. Luciferase reporter gene analysis also confirmed that SCAP is a target gene of MiR-125b. These observations suggest that MiR-125b regulates GCs steroid hormone secretion by targeting SCAP. Furthermore, lncGSAR overexpression increased proteins’ expression of SCAP and SREBP, whereas this effect was eliminated by MiR-125b mimics. The dual-luciferase reporter assay and clips test also confirmed again that lncGSAR acts as an MiR-125b sponge to protect the SCAP from the attack of MiR-125b. This result further demonstrated that binding of lncGSAR to MiR-125b regulates GCs state and function by targeting SCAP. However, there are various models for lncRNAs’ functions, and the regulatory mechanism of lncGSAR in vivo requires further study.

4. Materials and Methods

4.1. Cell Collections and Culture

All procedures used in this study were approved by the Animal Care and Use Committee of Hebei Agricultural University (Hebei, P.R. China; permit number 2022100). Fresh ewe ovaries (from thin-tailed Han sheep, ages ranged from 1 to 1.5 years) were collected at the local abattoir (Baoding, Hebei, China) and transported to the laboratory within 3 h in a buffered saline solution supplemented with streptomycin/penicillin mixture (1%) maintained at 37 °C [47]. Small, immature follicles between 1 and 3 mm in diameter were punctured with a disposable syringe, with follicular fluid collected from numbers of ovaries (ovaries number > 50) to negate any individual animal effects. The follicle suspensions were pooled, and GCs were harvested immediately after centrifuging at 1000× g for 10 min. GCs were counted with a hemocytometer (Axio Vert. A1, Zeiss, Oberkochen, Germany). Then, the GCs were seeded in cell culture plates (Thermo Fisher Scientific, Waltham, MA, USA) at a density of 2 × 105/well and cultured in Dulbecco’s Modified Eagle Medium (DMEM/F12, Gibco, Carlsbad, CA, USA) supplemented with 10% fetal bovine serum (FBS) (Gibco, USA) and 1% streptomycin/penicillin mixture in a humidified atmosphere at 37 °C and 5% CO2 for 48 h with the medium changed every 24 h.

4.2. The Granulosa Cell (GC) Model of Glucose Treatment

When cell confluence reached 80% medium, the original medium is removed. Cells were washed with 1 × PBS, and then, all treatments were cultured in DMEM, which was prepared with no serum, glucose, pyruvate, and phenol red (Solarbio, Beijing, China) for 8 h. Subsequently, the cells were supplemented with 8.4 mM and 33.6 mM of glucose cultured for an additional 24 h. The 8.4 mM represented an optimum glucose concentration for the secretion of steroid hormones by ovine GCs in vitro [6,7], and the 33.6 mM group represents 30 times the physiological concentration of glucose in follicular fluid and was used to detect changes in steroid hormones at super-physiological concentrations [48,49,50]. This culture system was developed so that GC retains hormonally responsive aromatase activity and does not luteinize with time in culture [51,52,53,54]. The GCs were collected for subsequent RNA-seq after the 24 h treatment period.

4.3. RNA Sequencing (RNA-Seq) and Bioinformatics Analysis

Ovine GCs subjected to optimum groups (n = 3, add 8.4 mM glucose groups) and high groups (n = 3, add 33.6 mM glucose groups) were used for RNA-seq. The total RNA of each sample was isolated using TRIzol reagent (Invitrogen, Life Technologies, Carlsbad, CA, USA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 System (Agilent Technologies, Santa Rosa, CA, USA). Library preparation and Illumina sequencing analysis were performed as previously described, and putative lncRNAs were screened using unknown transcripts [7].
To investigate interactions between lncRNAs and mRNAs, we constructed a complementary pair network comprising mRNA and lncRNA using Cytoscape 3.6.1(By the National Institute of General Medical Sciences, USA) (https://cytoscape.org/, accessed on 11 September 2018) [55]. Heatmaps were generated by using the R package. The raw sequencing dataset supporting the results of this study was submitted to NCBI BioProject (PRJNA825818) (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE200668, accessed on 12 April 2022).

4.4. RNA Extraction, Complementary DNA (cDNA) Synthesis, and (Quantitative Real-Time PCR) qRT-PCR

Total RNA was extracted from cultured cells according to the manufacturer’s instructions and supplied with the TRIzol Reagent (Invitrogen, Life Technologies, Carlsbad, CA, USA); cDNA synthesis for RNA (mRNA and lncRNA) was carried out using the Prime Script RT Reagent Kit with gDNA Eraser (Perfect Real Time) (TaKaRa, Otsu, Japan). Primers were designed and synthesized by GenechemBio (Shanghai, China), and primers used for quantitative real-time PCR are listed in Table S1. Real-time quantitative PCR reactions were performed on a Bio-Rad CFX96 Real-Time Detection System using an iTaq Universal SYBR Green Supermix Kit (Bio-Rad Laboratories Inc., Hercules, CA, USA). Data analyses were performed using the 2ΔΔCt method. β-actin was used as an internal control for mRNA and lncRNA, and U6 was used as an internal control for miRNA.

4.5. Nuclear and Cytoplasmic RNA Fractionation Assay

For nuclear and cytoplasmic RNA separation, 1 × 106 cells were collected and extracted using a Paris kit (Life Technologies, Pleasanton, CA, USA), according to the manufacturer’s instructions. U6 and GAPDH were used as positive controls for the nucleus and cytoplasm, respectively.

4.6. Plasmid Construction, RNA Oligonucleotides, and Cell Transfection

The granulosa cells were incubated overnight (at 60–70%, confluence), and transfection or co-transfection was performed using Lipofectamine 2000 (Invitrogen, Shanghai, China) for 48 h. The SCAP overexpression construct was generated by amplifying the SCAP coding sequence, which was subsequently integrated into the HindIII/KpnI restriction sites of the pcDNA3.1 overexpression plasmid (named pcDNA3.1-SCAP). The SREBPs overexpression plasmid was generated by amplifying the SREBPs coding sequence, which was subsequently integrated into the HindIII/KpnI restriction sites of the pcDNA3.1 overexpression plasmid (named pcDNA3.1-SREBPs). The lncRNA GSAR sequence was amplified by PCR and then subcloned into the pCDNA3.1 vector, generating pCDNA3.1-lncGSAR. The empty pcDNA3.1 vector was used as control plasmid.
MiR-125b binding sites in SCAP 3′UTR or lncGSAR were amplified by PCR using a cDNA template synthesized from total RNA. Then, the PCR products were subcloned into XhoI/XbaI restriction sites in the pmirGLO dual-luciferase reporter vector to generate the pmirGLO-Luc-SCAP reporter and pmirGLO-Luc-lncGSAR reporter.
The MiR-125b mimics, mimic-negative control, siRNA target against the SCAP gene (si-SCAP), siRNA target against the lncGSAR gene (si-lncGSAR), and siRNA nonspecific control were designed and synthesized by RiboBio (Guangzhou, China).
The GCs were incubated overnight (at 60–70%, confluence). All transient transfections were performed with Lipofectamine 2000 reagent (Invitrogen, Carlsbad, CA, USA), according to manufacturer’s direction. The primers and oligonucleotide sequences used in this study are listed in Tables S1 and S2.

4.7. CCK-8 Assay

Primary GCs were seeded in 96-well plates and cultured in growth medium. After being transfected, cell proliferation was monitored using a TransDetect CCK (TransGen Biotech, Beijing, China), according to the manufacturer’s protocol. Absorbance was measured using a Model 680 Microplate Reader (Bio-Rad, Hercules, CA, USA) by optical density at a wavelength of 450 nm.

4.8. Flow Cytometric Analysis

For the flow cytometry analysis of the cell cycle, GCs were seeded in 12-well plates. When the cells grew to a density of 50% confluence, they were transfected with overexpression plasmids. After transfection for 24 h, the cells were collected and fixed overnight in 70% ethanol at 4 °C. Subsequently, the fixed cells were stained with a 50 µg/mL propidium iodide solution (Sigma Life Science, St. Louis, MO, USA) containing 10 µg/mL RNase A (Takara, Japan) and 0.2% (v/v) Triton X-100 (Sigma Life Science, St. Louis, MO, USA) and then incubated in the dark at 37 °C for 30 min. Flow cytometry analysis was performed on a BD Accuri C6 flow cytometer (BD Biosciences, San Jose, CA, USA), and data were processed using the FlowJo7.6 software (Stanford University, Stanford, CA, USA).

4.9. 5′-bromo-2′-deoxyuridine (BrdU) Detection

BrdU (Sigma, #B5002) incorporation assay was carried out by following the protocol. Briefly, BrdU was diluted to a final concentration of 0.03 mg/mL with fresh DMEM and then applied onto the cells grown on slices. Cells were incubated with 1.5 M HCl followed by 5-min fixation in 70% cold ethanol. Fluorescence staining with anti-BrdU antibody was then conducted. The slices images were captured by a Zeiss Axio Observer confocal microscope.

4.10. Quantification of E2 and P4 in Culture Medium

Granulosa cells were seeded in six-well plates and transfected with overexpression plasmid, siRNA, or miRNA mimics for 48 h, and the medium was collected for subsequent measurements. The E2 concentration of the culture medium was analyzed using a commercial radioimmunoassay (RIA) kit (H102, Nanjing Jiancheng Bioengineering Institute, Nanjing, China). The P4 concentration was also measured using RIA kit (H089, Nanjing Jiancheng Bioengineering Institute, Nanjing, China).

4.11. Western Blotting

Granulosa cells were seeded in six-well plates and transfected with overexpression plasmid, siRNA, or miRNA mimics for 48 h. Cells were harvested, washed with 1 × phosphate buffered saline (PBS), and lysed in RIPA lysis buffer. Equal amount of proteins were separated at sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). After incubation with the indicated primary and secondary antibodies, signals were visualized by ECL. Membrane was then ready for scanning by Image studio system. Protein quantification was conducted by ImageJ software. The primary antibodies used were anti-CYP11A1 (1:1000; catalogue no. ab67355, Abcam, Cambridge, UK), anti-CYP19A1 (1:1000; catalogue no. ab18995, Abcam, Cambridge, UK), anti-SCAP (1:1000; catalogue no. ab190103, Abcam, Cambridge, UK), anti-SREBP1 (1:1000; catalogue no. ab3259, Abcam, Cambridge, UK), anti-Argonaute-2 (1:1000; catalogue no. ab186733, Abcam, Cambridge, UK), and anti-β-actin (1:10,000; catalogue no. 66009-1-Ig, Proteintech, Chicago, IL, USA). The goat anti-rabbit IgG (H + L)-HRP (1:5000; catalogue no. ab6721, Abcam, Cambridge, UK) were used as a secondary antibody.

4.12. Dual-Luciferase Reporter Assay

Approximately 3 × 104 GCs were plated onto 24-well tissue culture plates 24 h before transfection. Cells were transfected with a mixture of Renilla luciferase and indicated luciferase reporters using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). Forty-eight hours after transfection, the cells were harvested and subjected to an assay by using the Dual Luciferase Reporter Assay system (Promega, Madison, WI, USA). The luciferase activity was detected by Fluorescence/MultiDetection Microplate Reader (BioTek, Winooski, VT, USA). The relative luciferase activities were normalized with the Renilla luciferase activities.

4.13. In Vitro RNA-RNA Interaction Pull-Down Assay

Biotin-labeled RNA was synthesized in vitro using Biotin RNA Labeling Mix (Roche, St Louis, MO, USA, 11685597910). After treatment with RNase-free DNase I, Biotin-labeled RNA was heated at 70 °C for 15 min followed by 2 min incubation on ice to recover the secondary structure of RNA. The RNA was then incubated with streptavidin agarose beads (Invitrogen, Carlsbad, CA, USA) overnight. The next day, the RNA-RNA complexes were pulled down and collected by streptavidin agarose beads. After that, the immunoprecipitated RNA was eluted, isolated, and reverse transcribed to cDNA for the subsequent qRT-PCR analysis.

4.14. Statistical Analysis

All experiments were performed at least three times. Data are presented as means ± standard error of the mean based on three independent experiments. All data were normally distributed, and variance was similar between the statistically compared groups. Statistical analyses were performed using SPSS version 22.0 (SPSS Inc., Chicago, IL, USA). Statistical differences were determined by one-way analysis of variance (ANOVA). Tukey’s test was used to estimate the significance of the results. A p-value < 0.05 were considered statistically significant.

5. Conclusions

In conclusion, we identified a novel lncGSAR that acts as a sponge for MiR-125b to activate the SCAP/SREBP pathway, resulting in promoting granulose cell proliferation and steroidogenesis. Therefore, this study identified a candidate lncRNA (lncGSAR) involved in ovine fecundity, providing insights into the regulatory mechanisms by which glucose regulates follicular development and a basis for new strategies regulating animal reproduction by nutrients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232012132/s1, Table S1: Information of Primers; Table S2: Oligonucleotide sequences in this study; Table S3: The information of differentially expressed 461 lncRNA transcripts (281 upregulated and 180 downregulated) in the optimum and high group; Table S4: The information of differentially expressed 94 miRNA transcripts in the optimum and high group; Table S5: The information of differentially expressed 796 mRNA transcripts in the optimum and high group; Table S6: The information of ceRNA (lncGSAR–miRNA–mRNA) Network; Figure S1: Gene ontology term enrichment analysis of differentially expressed lncRNA target genes in different glucose treatment groups; Figure S2: Kyoto Encyclopedia of Genes and Genomes enrichment analysis of differentially expressed lncRNA target genes in different glucose treatment groups; Figure S3: RNA sequence of lncGSAR.

Author Contributions

Conceptualization, Y.Z.; data curation, Y.W.; formal analysis, Y.W. and Y.L.; funding acquisition, Y.L. and Y.Z.; investigation, J.L., Y.G. and C.D.; methodology, Y.W. and Y.G.; supervision, Y.L. and Y.Z.; validation, C.D. and H.Y.; visualization, Y.W., Y.G. and S.J.; writing—original draft, Y.W.; writing—review and editing, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the China Agriculture Research System (CARS-38) and (CARS-39-23).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Administration of Animal Care and Use and was approved by the Animal Ethics Committee of Hebei Agricultural University, Baoding, China (2022100, 10 March 2022).

Data Availability Statement

The datasets used during the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jahromi, B.N.; Mosallanezhad, Z.; Matloob, N.; Davari, M.; Ghobadifar, M.A. The potential role of granulosa cells in the maturation rate of immature human oocytes and embryo development: A co-culture study. Clin. Exp. Reprod. Med. 2015, 42, 111–117. [Google Scholar] [CrossRef] [PubMed]
  2. McGee, E.A.; Hsueh, A.J. Initial and cyclic recruitment of ovarian follicles. Endocr. Rev. 2000, 21, 200–214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Dumesic, D.A.; Meldrum, D.R.; Katz-Jaffe, M.G.; Krisher, R.L.; Schoolcraft, W.B. Oocyte environment: Follicular fluid and cumulus cells are critical for oocyte health. Fertil. Steril. 2015, 103, 303–316. [Google Scholar] [CrossRef] [PubMed]
  4. Charlier, C.; Montfort, J.; Chabrol, O.; Brisard, D.; Nguyen, T.; Aurélie, L.C.; Laurent, R.P.; François, M.; Pierre, P.; Svetlana, U.; et al. Oocyte-somatic cells interactions, lessons from evolution. BMC Genom. 2012, 13, 560. [Google Scholar] [CrossRef] [Green Version]
  5. Fraunhoffer, N.A.; Abuelafia, M.A.; Barrientos, M.A.; Cimerman, K.V.; Olmos, M.F.; Chuluyan, E.; Barrios, M. Long-term apoptosis-related protein expression in the diabetic mouse ovary. PLoS ONE 2018, 13, e0203268. [Google Scholar] [CrossRef] [Green Version]
  6. Wang, Y.; Duan, C.H.; Guo, Y.X.; Li, J.J.; He, H.Y.; Li, R.T.; Liu, Y.Q.; Zhang, Y.J. Effects of glucose on glycolysis and steroidogenesis as well as related gene expression in ovine granulosa cells in vitro. Small Rumin. Res. 2022, 215, 106766. [Google Scholar] [CrossRef]
  7. Wang, Y.; Guo, Y.X.; Duan, C.H.; Yang, R.C.; Zhang, L.C.; Liu, Y.Q.; Zhang, Y.J. Long Non-Coding RNA GDAR Regulates Ovine Granulosa Cells Apoptosis by Affecting the Expression of Apoptosis-Related Genes. Int. J. Mol. Sci. 2022, 23, 5183. [Google Scholar] [CrossRef]
  8. Djebali, S.; Davis, C.A.; Merkel, A.; Dobin, A.; Lassmann, T.; Mortazavi, A.; Tanzer, A.; Lagarde, J.; Lin, W.; Schlesinger, F.; et al. Landscape of transcription in human cells. Nature 2012, 489, 101–108. [Google Scholar] [CrossRef] [Green Version]
  9. Zhao, Z.; Zou, X.; Lu, T.; Deng, M.; Li, Y.; Guo, Y.; Sun, B.L.; Liu, G.B.; Liu, D.W. Identification of mRNAs and lncRNAs Involved in the Regulation of Follicle Development in Goat. Front. Genet. 2020, 11, 589076. [Google Scholar] [CrossRef]
  10. Wang, J.; Niu, M.; Zhang, T.; Shen, W.; Cao, H. Genome-Wide Network of lncRNA–mRNA During Ovine Oocyte Development from Germinal Vesicle to Metaphase II in vitro. Front. Physiol. 2020, 11, 1019. [Google Scholar] [CrossRef]
  11. Luo, J.; Wang, K.; Yeh, S.; Sun, Y.; Liang, L.; Xiao, Y.; Xu, W.H.; Niu, Y.J.; Cheng, L.; Sankar, N.M.; et al. LncRNA-p21 alters the antiandrogen enzalutamide-induced prostate cancer neuroendocrine differentiation via modulating the EZH2/STAT3 signaling. Nat. Commun. 2019, 10, 2571. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Huang, J.Y.; Zhao, J.; Geng, X.Y.; Chu, W.W.; Li, S.; Chen, Z.J.; Du, Y.Z. Long non-coding RNA lnc-CCNL1-3:1 promotes granulosa cell apoptosis and suppresses glucose uptake in women with polycystic ovary syndrome. Mol. Ther. Nucleic Acid. 2020, 23, 614–628. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, G.L.; Liu, S.X.; Xing, G.L.; Wang, F. lncRNA PVT1/MicroRNA-17-5p/PTEN Axis Regulates Secretion of E2 and P4, Proliferation, and Apoptosis of Ovarian Granulosa Cells in PCOS. Mol. Ther. Nucleic Acid. 2020, 20, 205–216. [Google Scholar] [CrossRef] [PubMed]
  14. Zhen, J.R.; Li, J.L.; Li, X.; Wang, X.; Xiao, Y.L.; Sun, Z.Y.; Yu, Q. Downregulating lncRNA NEAT1 induces proliferation and represses apoptosis of ovarian granulosa cells in polycystic ovary syndrome via microRNA-381/IGF1 axis. J. Biomed. Sci. 2021, 28, 53. [Google Scholar] [CrossRef]
  15. Salmena, L.; Poliseno, L.; Tay, Y.; Kats, L.; Pandolfi, P.P. A ceRNA hypothesis: The Rosetta Stone of a hidden RNA language? Cell 2011, 146, 353–358. [Google Scholar] [CrossRef] [Green Version]
  16. Yao, X.L.; Gao, X.X.; Bao, Y.J.; El-Samahy, M.A.; Yang, J.Y.; Wang, Z.B.; Li, X.D.; Zhang, G.M.; Zhang, Y.L.; Liu, W.J.; et al. lncRNA FDNCR promotes apoptosis of granulosa cells by targeting the miR-543-3p/DCN/TGF-b signaling pathway in Hu sheep. Mol. Ther. Nucleic Acid. 2021, 24, 223–240. [Google Scholar] [CrossRef]
  17. Wang, Y.; Pan, Z.X.; Du, X.; Zhang, J.B.; Liu, H.L.; Li, Q.F. NORHA, a novel follicular atresia-related lncRNA, promotes porcine granulosa cell apoptosis via the miR-183-96-182 cluster and FoxO1 axis. J. Anim. Sci. Biotechnol. 2021, 12, 103. [Google Scholar] [CrossRef]
  18. Cheng, C.M.; Geng, F.; Li, Z.; Zhong, Y.G.; Wang, H.B.; Cheng, X.; Zhao, Y.; Mo, X.K.; Horbinski, C.; Duan, W.R.; et al. Ammonia stimulates SCAP/Insig dissociation and SREBP-1 activation to promote lipogenesis and tumour growth. Nat. Metab. 2022, 4, 575–588. [Google Scholar] [CrossRef]
  19. Shimizu-Albergine, M.; Yserloo, V.B.; Golkowski, M.G.; Ong, S.E.; Beavo, J.A.; Bornfeldt, K.E. SCAP/SREBP pathway is required for the full steroidogenic response to cyclic AMP. Proc. Natl. Acad. Sci. USA 2016, 113, E5685–E5693. [Google Scholar] [CrossRef] [Green Version]
  20. Song, P.Y.; Yue, Q.X.; Fu, Q.; Li, X.Y.; Li, X.J.; Zhou, R.Y.; Chen, X.R.; Tao, C.Y. Integrated Analysis of miRNA-mRNA interaction in ovaries of Turpan Black Sheep during follicular and luteal phases. Reprod. Domest. Anim. 2020, 56, 46–57. [Google Scholar] [CrossRef]
  21. Yao, Y.L.; Niu, J.Q.; Sizhu, S.L.; Li, B.J.; Chen, Y.; Li, R.Y.; Yangzong, Q.B.; Li, Q.F.; Xu, Y.F. microRNA-125b Regulates Apoptosis by Targeting Bone Morphogenetic Protein Receptor 1B in Yak Granulosa Cells. DNA Cell Biol. 2018, 37, 878–887. [Google Scholar] [CrossRef] [PubMed]
  22. Brown, A.J.; Sun, L.; Feramisco, J.D.; Brown, M.S.; Goldstein, J.L. Cholesterol addition to ER membranes alters conformation of SCAP, the SREBP escort protein that regulates cholesterol metabolism. Mol. Cell 2002, 10, 237–245. [Google Scholar] [CrossRef]
  23. Walsh, S.W.; Matthews, D.; Browne, J.A.; Forde, N.; Crowe, M.A.; Mihm, M.; Diskin, M.; Evans, A.C.O. Acute dietary restriction in heifers alters expression of genes regulating exposure and response to gonadotrophins and IGF in dominant follicles. Anim. Reprod. Sci. 2012, 133, 43–51. [Google Scholar] [CrossRef] [PubMed]
  24. Comin, A.; Peric, T.; Montillo, M.; Cappa, A.; Marchi, V.; Veronesi, M.C.; Prandi, A. Luteal activity and effect of dietary energy restriction on follicular development in lactating cows. Reprod. Domest. Anim. 2017, 52, 632–639. [Google Scholar] [CrossRef]
  25. Ying, S.J.; Wang, Z.Y.; Wang, C.L.; Nie, H.T.; He, D.Y.; Jia, R.X.; Wu, Y.C.; Wan, Y.J.; Zhou, Z.R.; Yan, Y.B.; et al. Effect of different levels of short-term feed intake on folliculogenesis and follicular fluid and plasma concentrations of lactate dehydrogenase, glucose, and hormones in Hu sheep during the luteal phase. Reproduction 2011, 142, 699–710. [Google Scholar] [CrossRef] [Green Version]
  26. Brownlee, M. Biochemicstry and molecular cell biology of diabetic complications. Nature 2001, 414, 813–820. [Google Scholar] [CrossRef]
  27. Mauvais-Jarvis, F.; Clegg, D.J.; Hevener, A.L. The role of estrogens in control of energy balanceand glucose homeostasis. Endocr. Rev. 2013, 34, 309–338. [Google Scholar] [CrossRef] [Green Version]
  28. El-Bahr, S.M.; Ghoneim, I.M.; Waheed, M.M. Biochemical and hormonal analysis of follicular fluid and serum of female dromedary camels (Camelus dromedarius) with different sized ovarian follicles. Anim. Reprod. Sci. 2015, 159, 98–103. [Google Scholar] [CrossRef]
  29. Guo, Y.X.; Duan, C.H.; Hao, Q.H.; Liu, Y.Q.; Li, T.; Zhang, Y.J. Effect of short-term nutritional supplementation on hormone concentrations in ovarian follicular fluid and steroid regulating gene mRNA abundances in granulosa cells of ewes. Anim. Reprod. Sci. 2019, 211, 106208. [Google Scholar] [CrossRef]
  30. Li, X.L.; Zhang, T.Y.; Li, S.X.; Deng, Y.Y.; Wang, L.H.; Tao, T.; Wang, S.J.; Gu, Y.Y.; Gu, W.Q.; Hong, J.; et al. Correlation between glucose metabolism and serum steroid hormones in patients with polycystic ovary syndrome. Clin. Endocrinol. 2020, 92, 350–357. [Google Scholar] [CrossRef]
  31. Chabrolle, C.; JeanPierre, E.; Tosca, L.; Ramé, C.; Dupont, J. Effects of high levels of glucose on the steroidogenesis and the expression of adiponectin receptors in rat ovarian cells. Reprod. Biol. Endocrinol. 2008, 6, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Jiao, J.; Shi, B.; Wang, T.R.; Fang, Y.Y.; Cao, T.F.; Zhou, Y.M.; Wang, X.X.; Li, D. Characterization of long non-coding RNA and messenger RNA profiles in follicular fluid from mature and immature ovarian follicles of healthy women and women with polycystic ovary syndrome. Hum. Reprod. 2018, 33, 1735–1748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Lou, W.Y.; Ding, B.S.; Zhong, G.S.; Du, C.Y.; Fan, W.M.; Fu, P.F. Dysregulation of pseudogene/lncRNA-hsa-miR-363-3p-SPOCK2 pathway fuels stage progression of ovarian cancer. Aging 2019, 11, 11416–11439. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, X.J.; Zhang, P.; Li, Y.M.; Zhao, N.; Han, H.Y. The AMPK-mTOR axis requires increased MALAT1 expression for promoting granulosa cell proliferation in endometriosis. Exp. Ther. Med. 2021, 21, 21. [Google Scholar] [CrossRef]
  35. Ebert, M.S.; Sharp, P.A. Emerging roles for natural microRNA sponges. Curr. Biol. 2010, 20, R858–R861. [Google Scholar] [CrossRef] [Green Version]
  36. Karreth, F.A.; Pandolfi, P.P. ceRNA cross-talk in cancer: When ce-bling rivalries go awry. Cancer Discov. 2013, 3, 1113–1121. [Google Scholar] [CrossRef] [Green Version]
  37. Yan, H.; Li, H.; Li, P.Y.; Li, X.; Lin, J.J.; Zhu, L.L.; Silva, M.A.; Wang, X.F.; Wang, P.; Zhang, Z. Long noncoding RNA MLK7-AS1 promotes ovarian cancer cells progression by modulating miR-375/YAP1 axis. J. Exp. Clin. Cancer Res. 2018, 37, 237. [Google Scholar] [CrossRef]
  38. Jiang, X.Z.; Ye, Z.X.; Jiang, Y.F.; Yu, W.; Fang, Q. LncRNA OIP5-AS1 upregulates snail expression by sponging miR-34a to promote ovarian carcinoma cell invasion and migration. Biol. Res. 2020, 53, 49. [Google Scholar] [CrossRef]
  39. Wang, X.M.; Meng, K.; Wang, H.Q.; Wang, Y.; Zhao, Y.Q.; Kang, J.; Zhang, Y.; Quan, F.S. Identification of small extracellular vesicle subtypes in follicular fluid: Insights into the function and miRNA profiles. J. Cell. Physiol. 2021, 236, 5633–5645. [Google Scholar] [CrossRef]
  40. Guo, Y.; Wu, Y.B.; Shi, J.H.; Zhuang, H.; Ci, L.; Huang, Q.; Wan, Z.P.; Yang, H.; Zhang, M.J.; Tan, Y.T.; et al. miR-29a/b 1 Regulates the Luteinizing Hormone Secretion and Affects Mouse Ovulation. Front. Endocrinol. 2021, 12, 636220. [Google Scholar] [CrossRef]
  41. Li, G.X.; Chen, Y.; Jin, W.J.; Zhai, B.; Li, Y.F.; Sun, G.R.; Li, H.; Kang, X.T.; Tian, Y.D. Effects of MiR-125b-5p on Preadipocyte Proliferation and Differentiation in Chicken. Mol. Biol. Rep. 2021, 48, 491–502. [Google Scholar] [CrossRef] [PubMed]
  42. Yu, T.Z.; Li, G.P.; Wang, C.G.; Gong, G.Q.; Wang, L.W.; Li, C.Y.; Chen, Y.; Wang, X.L. MIR210HG regulates glycolysis, cell proliferation, and metastasis of pancreatic cancer cells through MiR-125b-5p/HK2/PKM2 axis. RNA Biol. 2021, 18, 2513–2530. [Google Scholar] [CrossRef] [PubMed]
  43. Sen, A.; Prizant, H.; Light, A.; Biswas, A.; Hayes, E.; Lee, H.J.; Barad, D.; Gleicher, N.; Hammes, S.R. Androgens regulate ovarian follicular development by increasing follicle stimulating hormone receptor and microRNA-125b expression. Proc. Natl. Acad. Sci. USA 2014, 111, 8–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. DeBose-Boyd, R.A.; Brown, M.S.; Li, W.P.; Nohturfft, A.; Goldstein, J.L.; Espenshade, P.J. Transport-dependent proteolysis of SREBP: Relocation of site-1 protease from Golgi to ER obviates the need for SREBP transport to Golgi. Cell 1999, 99, 703–712. [Google Scholar] [CrossRef] [Green Version]
  45. Espenshade, P.J.; Cheng, D.; Goldstein, J.L.; Brown, M.S. Autocatalytic processing of site-1 protease removes propeptide and permits cleavage of sterol regulatory element-binding proteins. J. Biol. Chem. 1999, 274, 22795–22804. [Google Scholar] [CrossRef] [Green Version]
  46. Espenshade, P.J.; Li, W.P.; Yabe, D. Sterols block binding of COPII proteins to SCAP, thereby controlling SCAP sorting in ER. Proc. Natl. Acad. Sci. USA 2002, 99, 11694–11699. [Google Scholar] [CrossRef] [Green Version]
  47. Li, T.; Ma, A.T.; Zhang, Y.J.; Liu, Y.Q.; Wang, S.Y.; Shi, X.F.; Gao, K.; Duan, C.H. Effects of silencing INHα gene by RNAi on the mRNA expressions of StAR, FST, INHβB, and FSHR genes in cultured sheep granulosa cells. Small Rumin. Res. 2017, 157, 23–26. [Google Scholar] [CrossRef]
  48. Somchit, A.; Campbell, B.K.; Khalid, M.; Kendall, N.R.; Scaramuzzi, R.J. The effect of short-term nutritional supplementation of ewes with lupin grain (Lupinus luteus), during the luteal phase of the estrous cycle on the number of ovarian follicles and the concentrations of hormones and glucose in plasma and follicular fluid. Theriogenology 2007, 68, 1037–1046. [Google Scholar] [CrossRef]
  49. Scaramuzzi, R.J.; Campbell, B.K.; Souza, C.J.H.; Baird, D.T. Glucose uptake and lactate production by the autotransplanted ovary of the ewe during the luteal and follicular phases of the oestrous cycle. Theriogenology 2010, 73, 1061–1067. [Google Scholar] [CrossRef]
  50. Campbell, B.K.; Onions, V.; Kendall, N.R.; Guo, L.; Scaramuzzi, R.J. The effect of monosaccharide sugars and pyruvate on the differentiation and metabolism of sheep granulosa cells in vitro. Reproduction 2010, 140, 541–550. [Google Scholar] [CrossRef]
  51. Langhout, D.J.; Spicer, L.J.; Geisert, R.D. Development of a culture system for bovine granulosa cells: Effects of growth hormone, estradiol, and gonadotropins on cell proliferation, steroidogenesis, and protein synthesis. J. Anim. Sci. 1991, 69, 3321–3334. [Google Scholar] [CrossRef] [PubMed]
  52. Spicer, L.J.; Chamberlain, C.S. Influence of Cortisol on Insulin and Insulin-Like Growth Factor 1 (IGF-1)-Induced Steroid Production and on IGF-1 Receptors in Cultured Bovine Granulosa Cells and Thecal Cells. Endocrine 1998, 9, 153–162. [Google Scholar] [CrossRef]
  53. Spicer, L.J.; Chamberlain, C.S.; Maciel, S.M. Influence of gonadotropins on insulin and insulin-like growth factor-I (IGF-I)-induced steroid production by bovine granulosa cells. Domest. Anim. Endocrinol. 2002, 22, 237–254. [Google Scholar] [CrossRef]
  54. Spicer, L.J. Growth differentiation factor-9 has divergent effects on proliferation and steroidogenesis of bovine granulosa cells. J. Endocrinol. 2006, 189, 329–339. [Google Scholar] [CrossRef] [Green Version]
  55. Saito, R.; Smoot, M.E.; Ono, K.; Ruscheinski, J.; Wang, P.L.; Lotia, S.; Pico, A.R.; Bader, G.D.; Ideker, T. A travel guide to Cytoscape plugins. Nat. Method. 2012, 9, 1069–1076. [Google Scholar] [CrossRef]
Figure 1. LncRNA-miRNA-mRNA interaction network. (AC) Heatmap of long non-coding RNAs (lncRNAs) (A), microRNA (miRNA) (B), and mRNAs (C) showing hierarchical clustering of changed lncRNAs, miRNA, and mRNAs of granulosa cells (GCs) in different glucose treatment groups; up- and downregulated genes are colored in red and blue, respectively. (D) Quantitative real-time PCR (qRT-PCR) validations and RNA sequencing (RNA-seq) of lncGSAR in GCs from in different glucose treatment groups. (E) LncGSAR is localized in the cytoplasm and nucleus. GAPDH and U6 serve as cytoplasmic and nuclear localization controls, respectively. (F) LncRNA-miRNA-mRNA interaction network consists of one lncRNA (red circle), one miRNA (green arrow), and nine mRNAs (blue squares). Values represent means ± SEM for three individuals. * p < 0.05 and ** p < 0.01.
Figure 1. LncRNA-miRNA-mRNA interaction network. (AC) Heatmap of long non-coding RNAs (lncRNAs) (A), microRNA (miRNA) (B), and mRNAs (C) showing hierarchical clustering of changed lncRNAs, miRNA, and mRNAs of granulosa cells (GCs) in different glucose treatment groups; up- and downregulated genes are colored in red and blue, respectively. (D) Quantitative real-time PCR (qRT-PCR) validations and RNA sequencing (RNA-seq) of lncGSAR in GCs from in different glucose treatment groups. (E) LncGSAR is localized in the cytoplasm and nucleus. GAPDH and U6 serve as cytoplasmic and nuclear localization controls, respectively. (F) LncRNA-miRNA-mRNA interaction network consists of one lncRNA (red circle), one miRNA (green arrow), and nine mRNAs (blue squares). Values represent means ± SEM for three individuals. * p < 0.05 and ** p < 0.01.
Ijms 23 12132 g001
Figure 2. LncGSAR is transactivated upon glucose stimulation. (A) The GCs were treated with actinomycin D to explore the effect of glucose on the expression of lncGSAR. (B) Dual-luciferase reporter assay was conducted of lncGSAR promoter sequence and full-length sequence. (C) Dual-luciferase reporter for the lncGSAR promoter. (D) Detection of lncGSAR full-length luciferase activity. Values represent means ± SEM for three individuals. ** p < 0.01. N.S., not significant.
Figure 2. LncGSAR is transactivated upon glucose stimulation. (A) The GCs were treated with actinomycin D to explore the effect of glucose on the expression of lncGSAR. (B) Dual-luciferase reporter assay was conducted of lncGSAR promoter sequence and full-length sequence. (C) Dual-luciferase reporter for the lncGSAR promoter. (D) Detection of lncGSAR full-length luciferase activity. Values represent means ± SEM for three individuals. ** p < 0.01. N.S., not significant.
Ijms 23 12132 g002
Figure 3. LncGSAR promotes GCs steroidogenesis and proliferation. (A,B) CCK-8 assay was performed to assess the effect of lncGSAR overexpression and knockdown on GCs proliferation. (C) The 5′-bromo-2′-deoxyuridine (BrdU) analysis after transfection of pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids, i.e., si-lncGSAR and si-NC, in proliferating GCs; scale bars are 100 µm. (D,E) Determination of estradiol (E2) concentrations in GCs transfected with pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids, i.e., si-lncGSAR and si-NC. (F,G) The progesterone (P4) concentrations in GCs after transfection of pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids or si-lncGSAR and si-NC. (H,I) The qRT-PCR of relative expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs transfected with pcDNA3.1-lncGSAR and in GCs transfected with si-lncGSAR. (J) Overexpression and knockdown of lncGSAR affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. ** p < 0.01; *** p < 0.001.
Figure 3. LncGSAR promotes GCs steroidogenesis and proliferation. (A,B) CCK-8 assay was performed to assess the effect of lncGSAR overexpression and knockdown on GCs proliferation. (C) The 5′-bromo-2′-deoxyuridine (BrdU) analysis after transfection of pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids, i.e., si-lncGSAR and si-NC, in proliferating GCs; scale bars are 100 µm. (D,E) Determination of estradiol (E2) concentrations in GCs transfected with pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids, i.e., si-lncGSAR and si-NC. (F,G) The progesterone (P4) concentrations in GCs after transfection of pcDNA3.1-lncGSAR and pcDNA3.1 empty plasmids or si-lncGSAR and si-NC. (H,I) The qRT-PCR of relative expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs transfected with pcDNA3.1-lncGSAR and in GCs transfected with si-lncGSAR. (J) Overexpression and knockdown of lncGSAR affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. ** p < 0.01; *** p < 0.001.
Ijms 23 12132 g003
Figure 4. LncGSAR functions as a competing endogenous RNAs for MiR-125b. (A) The qRT-PCR validations and RNA-seq of MiR-125b in GCs from in different glucose treatment groups. (B) Pearson’s correlation was determined between lncGSAR and MiR-125b. (C,D) The expression levels of lncGSAR (C) and MiR-125b (D) in GCs after transfection of 0, 300, 500, or 800 ng of MiR-125b mimics. (E) Detection of lncGSAR expression levels after transfection of MiR-125b mimics or MiR-125b inhibitor. (F) Detection of lncGSAR expression levels after transfection of MiR-125b mimics in different glucose concentrations groups. (G) Schematic depicting the interactions of MiR-125b with wild-type lncGSAR (blue) and mutant lncGSAR (green). Red nucleotides indicate the seed sequence of MiR-125b. (H) The regulatory relationship between lncGSAR and MiR-125b was assessed using a dual-luciferase reporter gene assay. (I) The interaction of lncGSAR with MiR-125b was determined by RNA-RNA interaction pull-down assay. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01. N.S., not significant.
Figure 4. LncGSAR functions as a competing endogenous RNAs for MiR-125b. (A) The qRT-PCR validations and RNA-seq of MiR-125b in GCs from in different glucose treatment groups. (B) Pearson’s correlation was determined between lncGSAR and MiR-125b. (C,D) The expression levels of lncGSAR (C) and MiR-125b (D) in GCs after transfection of 0, 300, 500, or 800 ng of MiR-125b mimics. (E) Detection of lncGSAR expression levels after transfection of MiR-125b mimics or MiR-125b inhibitor. (F) Detection of lncGSAR expression levels after transfection of MiR-125b mimics in different glucose concentrations groups. (G) Schematic depicting the interactions of MiR-125b with wild-type lncGSAR (blue) and mutant lncGSAR (green). Red nucleotides indicate the seed sequence of MiR-125b. (H) The regulatory relationship between lncGSAR and MiR-125b was assessed using a dual-luciferase reporter gene assay. (I) The interaction of lncGSAR with MiR-125b was determined by RNA-RNA interaction pull-down assay. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01. N.S., not significant.
Ijms 23 12132 g004
Figure 5. MiR-125b inhibit GCs steroidogenesis and proliferation. (A) The cell counting kit-8 (CCK-8) assay of GCs at 48 h after transfection of MiR-125b mimics and MiR-NC. (B) BrdU analysis was performed to assess the effect of MiR-125b overexpression on GCs proliferation. (CH) Cell cycle analysis of GCs at 48 h after transfection of MiR-125b mimics plasmid. (I) Determination of E2 concentrations after 48 h of GCs transfected with MiR-125b mimics. (J) Determination of P4 concentrations after 48 h of GCs transfected with MiR-125b mimics. (K) The qRT-PCR of relative expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs transfected with MiR-125b mimics. (L) Overexpression of MiR-125b affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01. N.S., not significant.
Figure 5. MiR-125b inhibit GCs steroidogenesis and proliferation. (A) The cell counting kit-8 (CCK-8) assay of GCs at 48 h after transfection of MiR-125b mimics and MiR-NC. (B) BrdU analysis was performed to assess the effect of MiR-125b overexpression on GCs proliferation. (CH) Cell cycle analysis of GCs at 48 h after transfection of MiR-125b mimics plasmid. (I) Determination of E2 concentrations after 48 h of GCs transfected with MiR-125b mimics. (J) Determination of P4 concentrations after 48 h of GCs transfected with MiR-125b mimics. (K) The qRT-PCR of relative expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs transfected with MiR-125b mimics. (L) Overexpression of MiR-125b affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01. N.S., not significant.
Ijms 23 12132 g005
Figure 6. MiR-125b regulates GCs proliferation and steroidogenesis and by targeting SCAP-SREBP axis. (A) The qRT-PCR validations and RNA-seq of SCAP in GCs from in different glucose treatment groups. (B) Pearson’s correlation was determined between SCAP and MiR-125b. (C) MiR-125b or pcDNA3.1-lncGSAR and MiR-125 mimics were co-transfected into GCs for CCK-8 assay. (D,E) The concentrations of E2 and P4 determination in GCs transfected with MiR-125b mimics, i.e., pcDNA3.1-lncSCAP and MiR-125 mimics. (F) CCK-8 assay of GCs at 48 h after transfection of MiR-125b mimics plasmid or pcDNA3.1-lncSREBP and MiR-125 mimics. (G,H) The concentrations of E2 and P4 determination in GCs transfected with MiR-125b mimics, i.e., pcDNA3.1-lncSREBP and MiR-125 mimics. (I,J) Western blot detection of SCAP (I) and SREBP (J) in GCs after overexpression of MiR-125b. (K) Schematic depicting the interaction of MiR-125b with wild-type (blue) and mutant SCAP (green). Red nucleotides indicate the seed sequence of MiR-125b. (L) The regulatory relationship between MiR-125b and SCAP was assessed using a dual-luciferase reporter gene assay. * p < 0.05; ** p < 0.01. N.S., not significant.
Figure 6. MiR-125b regulates GCs proliferation and steroidogenesis and by targeting SCAP-SREBP axis. (A) The qRT-PCR validations and RNA-seq of SCAP in GCs from in different glucose treatment groups. (B) Pearson’s correlation was determined between SCAP and MiR-125b. (C) MiR-125b or pcDNA3.1-lncGSAR and MiR-125 mimics were co-transfected into GCs for CCK-8 assay. (D,E) The concentrations of E2 and P4 determination in GCs transfected with MiR-125b mimics, i.e., pcDNA3.1-lncSCAP and MiR-125 mimics. (F) CCK-8 assay of GCs at 48 h after transfection of MiR-125b mimics plasmid or pcDNA3.1-lncSREBP and MiR-125 mimics. (G,H) The concentrations of E2 and P4 determination in GCs transfected with MiR-125b mimics, i.e., pcDNA3.1-lncSREBP and MiR-125 mimics. (I,J) Western blot detection of SCAP (I) and SREBP (J) in GCs after overexpression of MiR-125b. (K) Schematic depicting the interaction of MiR-125b with wild-type (blue) and mutant SCAP (green). Red nucleotides indicate the seed sequence of MiR-125b. (L) The regulatory relationship between MiR-125b and SCAP was assessed using a dual-luciferase reporter gene assay. * p < 0.05; ** p < 0.01. N.S., not significant.
Ijms 23 12132 g006
Figure 7. SCAP promotes GCs steroidogenesis and proliferation. (A,B) CCK-8 assay was performed to assess the effect of SCAP overexpression and knockdown on GCs proliferation. (C) BrdU analysis was performed to assess the effect of SCAP overexpression and knockdown on GCs proliferation; scale bars are 100 µm. (D,E) Detection of E2 concentrations in GCs transfected with pcDNA3.1-SCAP and pcDNA3.1 empty plasmids, i.e., si-SCAP and si-NC. (F,G) Detection of P4 concentrations in GCs transfected with pcDNA3.1-SCAP and pcDNA3.1 empty plasmids, i.e., si-SCAP and si-NC. (H,I) RT-qPCR was used to examine the effect of overexpression and knockdown of SCAP on the expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs. (J) Overexpression and knockdown of SCAP affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01; *** p < 0.001. N.S., not significant.
Figure 7. SCAP promotes GCs steroidogenesis and proliferation. (A,B) CCK-8 assay was performed to assess the effect of SCAP overexpression and knockdown on GCs proliferation. (C) BrdU analysis was performed to assess the effect of SCAP overexpression and knockdown on GCs proliferation; scale bars are 100 µm. (D,E) Detection of E2 concentrations in GCs transfected with pcDNA3.1-SCAP and pcDNA3.1 empty plasmids, i.e., si-SCAP and si-NC. (F,G) Detection of P4 concentrations in GCs transfected with pcDNA3.1-SCAP and pcDNA3.1 empty plasmids, i.e., si-SCAP and si-NC. (H,I) RT-qPCR was used to examine the effect of overexpression and knockdown of SCAP on the expression levels of steroidogenesis-related (CYP11A1, CYP19A1, and 3β-HSD) and sterol regulation-related (SCAP and SREBP) mRNAs in GCs. (J) Overexpression and knockdown of SCAP affected the expression of steroidogenesis-related proteins and sterol regulatory element proteins. Values represent means ± SEM for three individuals. * p < 0.05; ** p < 0.01; *** p < 0.001. N.S., not significant.
Ijms 23 12132 g007
Figure 8. LncGSAR functions as a ceRNA to regulate SCAP-SREBP pathway by sponging MiR-125b, thus promoting GCs proliferation and steroidogenesis. (A) Western blot detection of steroidogenesis-related and sterol regulatory element proteins. (B) Dual-luciferase reporter demonstrated the binding relationship among lncGSAR, MiR-125b and SCAP. (C) The interactions of lncGSAR with MiR-125b and MiR-125b with SCAP were determined by RNA-RNA interaction pull-down assay. Values represent means ± SEM for three individuals. ** p < 0.01. N.S., not significant.
Figure 8. LncGSAR functions as a ceRNA to regulate SCAP-SREBP pathway by sponging MiR-125b, thus promoting GCs proliferation and steroidogenesis. (A) Western blot detection of steroidogenesis-related and sterol regulatory element proteins. (B) Dual-luciferase reporter demonstrated the binding relationship among lncGSAR, MiR-125b and SCAP. (C) The interactions of lncGSAR with MiR-125b and MiR-125b with SCAP were determined by RNA-RNA interaction pull-down assay. Values represent means ± SEM for three individuals. ** p < 0.01. N.S., not significant.
Ijms 23 12132 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Guo, Y.; Duan, C.; Li, J.; Ji, S.; Yan, H.; Liu, Y.; Zhang, Y. LncGSAR Controls Ovarian Granulosa Cell Steroidogenesis via Sponging MiR-125b to Activate SCAP/SREBP Pathway. Int. J. Mol. Sci. 2022, 23, 12132. https://doi.org/10.3390/ijms232012132

AMA Style

Wang Y, Guo Y, Duan C, Li J, Ji S, Yan H, Liu Y, Zhang Y. LncGSAR Controls Ovarian Granulosa Cell Steroidogenesis via Sponging MiR-125b to Activate SCAP/SREBP Pathway. International Journal of Molecular Sciences. 2022; 23(20):12132. https://doi.org/10.3390/ijms232012132

Chicago/Turabian Style

Wang, Yong, Yunxia Guo, Chunhui Duan, Junjie Li, Shoukun Ji, Huihui Yan, Yueqin Liu, and Yingjie Zhang. 2022. "LncGSAR Controls Ovarian Granulosa Cell Steroidogenesis via Sponging MiR-125b to Activate SCAP/SREBP Pathway" International Journal of Molecular Sciences 23, no. 20: 12132. https://doi.org/10.3390/ijms232012132

APA Style

Wang, Y., Guo, Y., Duan, C., Li, J., Ji, S., Yan, H., Liu, Y., & Zhang, Y. (2022). LncGSAR Controls Ovarian Granulosa Cell Steroidogenesis via Sponging MiR-125b to Activate SCAP/SREBP Pathway. International Journal of Molecular Sciences, 23(20), 12132. https://doi.org/10.3390/ijms232012132

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop