LncPrep + 96kb 2.2 kb Inhibits Estradiol Secretion From Granulosa Cells by Inducing EDF1 Translocation

Introduction

Ovarian follicles are the basic functional units of ovaries. In mammals, ovarian follicular development is a continuous process throughout the reproductive life span. Follicles progress through the primordial, primary, and secondary stages before developing an antral cavity. With further growth and differentiation, preovulatory follicles form, and oocytes are released after being stimulated with luteinizing hormone (LH) (Zhang et al., 2011).

The growth and maturation of follicles are always accompanied by granulosa cell proliferation and differentiation (Johnson, 2015). Granulosa cells can acquire specific functional characteristics after differentiation, including the expression of specific receptors and the secretion of progesterone and estrogen to coordinate oocyte maturation and subsequent events such as ovulation, fertilization, and early embryo development (Sasson et al., 2004; Wang Q. et al., 2013).

Estrogens are the main steroids produced by granulosa cells in mature follicles. Ovarian estrogens play key roles in the female reproductive system and secondary sexual characteristics. In the reproductive system, estrogens regulate follicular development (Baba et al., 2017). It is well known that ovarian follicles can develop until the preantral and early antral stages independently of follicle-stimulating hormone (FSH). Estrogens stimulate the growth of primary and secondary follicles until the preantral stage (Laven and Fauser, 2006). Studies have shown that estrogens increase the number and size of ovarian follicles in vivo and the size of rat, mouse, and bovine follicles in vitro (Rosenfeld et al., 2001). In a steroid-depleted milieu, estradiol also improved secondary follicle survival, growth, antrum formation, and oocyte health in adult rhesus macaques (Ting et al., 2015). However, other studies showed that estrogens had inhibitory effect on follicle development in some species. In vivo diethylstilbestrol (DES) exposure induces atresia of the dominant follicles in rats and monkeys (Hutz et al., 1986). In vitro high-dose estradiol exposure delayed or inhibited oocyte meiotic maturation during murine secondary follicle culture (Tarumi et al., 2014). Therefore, the effects of estrogens can be stimulatory or inhibitory, depending on the species, dose, and duration of exposure. Studies have also shown that estrogens can lead to synchronization of follicular development in prepubertal mice (Wu et al., 2015). Studies of the underlying mechanism have shown that estrogens stimulate the proliferation of granulosa cells, have antiapoptotic effects on granulosa cells, and control follicular neovascularization (Li et al., 2013; Chou and Chen, 2018). Besides the effects on follicle development, estrogens exert negative feedback on LH production in the early stage of the menstrual cycle. When estrogen levels reach a critical level, they exert positive feedback on LH production and ovulation (Drummond, 2006; Holesh et al., 2020).

Estradiol, one main estrogen, is synthesized by granulosa cells. Aromatase in granulosa cells converts androgens (from theca cells) into estradiol. The loss of this enzyme is responsible for mouse infertility (Fisher et al., 1998). The expression of aromatase is controlled in a cell-specific, temporal, and spatial manner, which limits estradiol production. Studies have shown that aromatase is expressed at low levels in the granulosa cells of small growing follicles during the fetal and neonatal periods, and estradiol production is limited at this stage. In 21-day-old rats, aromatase is no longer expressed in growing follicles, but is expressed in only healthy large antral follicles and preovulatory follicles. Steroidogenesis can occur in only the postpubertal period. Therefore, the granulosa cells of large antral follicles and preovulatory follicles are the main cells that produce estradiol (Guigon et al., 2003; Sakurada et al., 2006; Stocco, 2008; Dewailly et al., 2016).

It is well accepted that FSH is the major inducer of aromatase activity in granulosa cells. In addition, other endocrine factors, paracrine factors, and autocrine factors all contribute to estrogen production (Stocco, 2008; Hobeika et al., 2019). For example, insulinlike growth factor 1 (IGF1) alone increased estradiol production to a level that was comparable to the level induced by FSH. Follicle-stimulating hormone and IGF1 have a synergistic effect on estradiol production (Sahmi et al., 2014). The proximal ends of the aromatase promoter contain a number of fairly well-conserved regulatory elements, including steroidogenic factor 1 (SF1, NR5A1), liver receptor homolog 1 (LRH-1, NR5A2), adenosine 3,5-cyclic monophosphoric acid (cAMP), and forkhead transcription factor 2 (Foxl2), which contribute to the regulation of aromatase and estradiol production (Stocco, 2008).

ENCODE (Encyclopedia of DNA Elements) shows that the vast majority (80.4%) of the human genome can be transcribed into RNA. Non-coding RNAs account for most of the genome (Consortium, 2012; Dhahbi, 2014) and include small non-coding RNAs and long non-coding RNAs (lncRNAs). MiRNAs have been demonstrated to play pivotal roles in regulating the sequential recruitment, selection, and development of follicles, atresia, and ovulation (Hossain et al., 2012; Imbar and Eisenberg, 2014; Li et al., 2015). Over the last few years, lncRNAs have been demonstrated to have important biological functions, which are involved in gene expression, and interact directly with biological macromolecules (DNAs, RNAs, and proteins) (Mercer and Mattick, 2013). Long non-coding RNAs have notably decreased sequence conservation compared to protein coding genes, but their secondary structure, rather than their sequence, usually dictates their function (Pang et al., 2006; Struhl, 2007; Mercer and Mattick, 2013). Long non-coding RNAs are implicated in a variety of biological processes, including dosage compensation, imprinting, cell cycle control, cell development, and gametogenesis (Yang et al., 2014). Moreover, the spatial and temporal-specific expression of lncRNAs shows a strong association with the cell differentiation process (Cesana et al., 2011; Wang Y. et al., 2013). Studies have shown that some lncRNAs are highly expressed in the ovaries. For instance, Neat1 plays an important role in corpus luteum formation, as Neat1 knockout mice fail to become pregnant (Nakagawa et al., 2014). Long non-coding RNA SRA was also demonstrated to promote growth, inhibit apoptosis, and induce the production of estradiol and progesterone in granulosa cells. Abnormal expression of lncRNA SRA may be a risk factor for polycystic ovary syndrome (Li et al., 2018). Long non-coding RNA-Amhr2 regulates Amhr2 gene activation in mouse ovarian granulosa cells (Kimura et al., 2017).

A novel lncRNA, termed lncPrep + 96kb, has been identified as the product of conserved non-coding sequences in the intergenic region adjacent to prolyl oligopeptidase (Matsubara et al., 2014). LncPrep + 96kb consists of two transcripts of different sizes, 2.2 and 2.8 kb. Matsubara et al. (2014) found that this lncRNA was highly expressed in murine granulosa cells. However, the physiological function of lncPrep + 96kb in granulosa cells is still unclear. In this study, we found that lncPrep + 96kb was mainly expressed in granulosa cells of small growing follicles and identified the inhibitory effect of lncPrep + 96kb 2.2 kb on aromatase expression and estradiol production in granulosa cells. We also further explored the possible underlying mechanism.

Materials and Methods

Animals

Kunming mice were purchased from the Animal Facility of Nanchang University and housed in a temperature- and light-controlled facility with free access to water and food. The experimental protocols were approved by the ethical committee of Nanchang University (SYXK2015-0001, December 29, 2015).

Chemicals

An antibody against endothelial differentiation-related factor 1 (EDF1) (BS70093) was obtained from Bio World, Nanjing, China. An antiaromatase antibody (BA3704) was obtained from BOSTER, Wuhan, China. Real-time polymerase chain reaction (PCR) mix (208054) was obtained from QIAGEN, Germany. ANTI-FLAG^®M2 magnetic beads (M8823) were obtained from Sigma, St. Louis, MO, United States. Phusion^® high-fidelity DNA polymerase and endonuclease were obtained from New England Biolab, Ipswich, MA, United States. FuGENE 6 transfection reagent (E2691) was purchased from Promega, Madison, WI, United States. An anti–β-actin antibody (HC201), anti-Flag antibody (HT201), and a CCK8 kit (FC101) were obtained from TransGen Biotech, Beijing, China. The digoxigenin (DIG) RNA labeling mix (11277073910) and anti-DIG-Fab antibody (110932749100) were obtained from Roche, Mannheim, Germany. pLV[shRNA]-EGFP:T2A:Puro-U6 > mEdf1 plasmid was obtained from Vector Builder, Guangzhou, China.

In situ Hybridization

The probe sequence of lncPrep + 96kb was amplified using Phusion^® high-fidelity DNA polymerase and ligated into the pGEMT vector. The DIG-labeled RNA probe was synthesized as previously reported (Zhang et al., 2011). Ovaries were fixed in 4% paraformaldehyde for 12 h and embedded in paraffin. Then, 5 μm sections were cut. After deparaffinization and rehydration, the sections were permeabilized with proteinase K (10 mg/L) in proteinase K buffer for 25 min at 37°C, incubated in prehybridization buffer for 4 h at 42°C, and incubated with the hybridization solutions, which were composed of DIG-labeled antisense RNA probes and hybridization buffer, at 56°C for 18 h. Detection of the hybridized probe was performed using an alkaline phosphatase–conjugated anti-DIG antibody (Fab fragment). After washing, the color was developed with color-generating solution (4-nitroblue tetrazolium chloride and 5-bromo-4-chloro-3-indolyl-phosphate,4-toluidinesalt) for 2 h.

CatRAPID Prediction of RNA–Protein Interactions

We used the CatRAPID online server to evaluate the interaction between lncPrep + 96kb 2.2 kb and EDF1 based on the secondary structure, hydrogen bonding, and van der Waals interactions, which rapidly predict the interactions between RNA and proteins (Agostini et al., 2013; Livi et al., 2016; Shang et al., 2019). The binding capacity of the EDF1 protein was estimated via algorithms using the CatRAPID signature module. The binding sites between lncPrep + 96kb 2.2 kb and EDF1 were evaluated by the CatRAPID fragments module.

Plasmid Constructs

To construct the lncPrep + 96kb 2.2 kb overexpression plasmid, the whole sequence of lncPrep + 96kb 2.2 kb was amplified by PCR, digested with the KpnI and EcoRI restriction enzymes, and cloned into the PCDNA3.0 plasmid vector. Bridge PCR was employed to construct the lncPrep + 96kb 2.2 kb mutation plasmids. The primers are shown in Table 1. The guide RNA (gRNA) targeted against lncPrep + 96kb was cloned into the pSpCas9(BB)-2A-GFP vector (PX458). We designed the gRNAs online using a web server¹, and the oligonucleotide pairs for each gRNA were subsequently synthesized (Table 2). The PX458 plasmid was digested with BbsI and annealed to the gRNAs. DNA sequencing was used to confirm successful insertion.

TABLE 1

Table 1. Oligonucleotides used for constructing plasmids.

TABLE 2

Table 2. Sequences of gRNAs designed online.

Granulosa Cell Culture and Transfection

The isolation of granulosa cells was performed as follows. Female mice (21 days) were injected with 5 IU pregnant mare serum gonadotropin (PMSG), and the ovaries were collected 48 h later. The fat and capsule tissues adhering to the ovaries were removed. Following two washes, the ovaries were mechanically dissected with a syringe needle in sterile phosphate-buffered saline (PBS). The suspension was collected and filtered through 75 μm strainers to remove debris and was centrifuged at 500 g for 5 min. The cells were resuspended in Dulbecco medium eagle medium/F12 culture medium supplemented with 5% fetal bovine serum, 100 IU/mL penicillin, and 100 μg/mL streptomycin sulfate. The cells were seeded and cultured overnight until they adhered to the cell culture plate. The plasmids were transfected into the granulosa cells using FuGENE 6 when the cell confluence reached 70%. After 8 h, the medium was changed. Then 100 nM androstenedione was added into the medium as a substrate of estradiol. The medium was harvested for estradiol examination, and the cells were lysed for mRNA and protein extraction at the indicated times.

Cell Viability Assay

Granulosa cells were seeded in 96-well plates at 5 × 10³ cells per well. The cells were transfected with the indicated plasmids for 48 h. Ten microliters of CCK-8 reagent was added to each well and incubated for 2 h. The absorbance of reduced WST-8(2-(2-methoxy-4-nitrophenyl)-3-(4-nitrophenyl)-5-(2,4-disulfophenyl)-2H-tetrazolium) was measured at 450 nm using an enzyme-linked immunosorbent assay plate reader.

Radioimmunoassay

The cultured granulosa cells were transfected with different plasmids. The supernatant was harvested at different times and stored at −80°C for subsequent hormone measurement. Estradiol was assessed by a commercial laboratory (Beijing Sino-uk Institute of Biological Technology, Beijing, China) using a commercial radioimmunoassay kit. For the kit, the cross-reactivities of steroid hormones and other peptides are less than 4%. The intra-assay and inter-assay variation coefficients did not exceed 10%. The sensitivity of estradiol detection was 2 pg/mL.

Immunofluorescence

Granulosa cells were fixed in 4% paraformaldehyde for 30 min. After washing with PBS, the slides were incubated with 5% bovine serum albumin in PBS for 30 min at room temperature to reduce non-specific binding and then incubated with a primary antibody against EDF1 (2.5 μg/mL) at 4°C overnight. After three rinses in PBS, the slides were incubated for 1 h at room temperature with goat anti–rabbit immunoglobulin G conjugated to fluorescein isothiocyanate. DAPI (4′,6-diamidino-2-phenylindole) was used to stain the nuclei. For the negative control, the slides were incubated with PBS instead of the primary antibody. Other procedures were exactly the same.

RNA Extraction and Real-Time PCR

Total RNA was isolated with TRIzol reagent according to the manufacturer’s instructions. Two micrograms of total RNA was used to synthesize cDNA using reverse transcriptase. Real-time PCR was carried out in a 20 μL reaction volume that consisted of 10 μL 2X Brilliant SYBR Green qPCR Master Mix, 1 μL cDNA, 0.5 μL primers, and 300 nM ROX reference dye using an ABI thermal cycler 7500 under the following conditions: 94°C for 30 s and 40 cycles of 94°C for 5 s and 60°C for 60 s. Melting curve analysis was conducted to analyze the purity of the products after real-time PCR. Experiments were repeated three times, and each gene was analyzed in parallel. The sequences for the specific primers are listed in Table 3. The threshold cycle (Ct) values were normalized to the mRNA expression level of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and calculated with the ΔΔCt method, which was then expressed as 2^–ΔΔCt for the sake of analyzing the fold change.

TABLE 3

Table 3. Oligonucleotides used for real-time PCR.

Protein Extraction and Western Blotting

Granulosa cells were washed with PBS and lysed in RIPA (radioimmune precipitation assay) lysis buffer containing complete Mini Protease Inhibitor Cocktail tablets (Roche, Mannheim, Germany). A Bradford assay (Bio-Rad Laboratories, Hercules, CA, United States) was employed to estimate the protein concentration. The proteins were separated on 12% sodium dodecyl sulfate–polyacrylamide gel electrophoresis gels and transferred to nitrocellulose membranes using a 200 mA constant current for 2 h. The blots were blocked with 5% skim milk for 1 h and incubated overnight at 4°C with a primary antibody (500 ng/mL for aromatase and EDF1). The membranes were washed in Tris-buffered saline Tween-20 three times and incubated with a horseradish peroxidase–conjugated secondary antibody at room temperature for 1 h. ECL substrate was used to visualize the target bands, and the band intensity was quantified by densitometry (Bio-Rad Image Lab). β-Actin served as an internal control.

RNA Immunoprecipitation and Mass Spectrometry

The interaction between lncPrep + 96kb 2.2 kb and protein was detected by RNA immunoprecipitation (RIP) based on the specific interaction between the MS2 coat protein and MS2 binding sites (MS2bs) as reported (Li et al., 2017). The pcDNA3.0-12 × MS2bs and pcDNA3.0-Flag-2 × MS2 plasmids were gifts from Professor Zheng Xiaofei (Beijing Institute of Radiation Medicine, Beijing, China). We amplified the full length of lncPrep + 96kb 2.2 kb and cloned it into pcDNA3.0-12 × MS2bs, which was named pcDNA3.0-lncPrep + 96kb 2.2 kb-12 × MS2bs. The pcDNA3.0-Flag-2 × MS2 and pcDNA3.0-lncPrep + 96kb 2.2 kb-12 × MS2bs plasmids were cotransfected into granulosa cells. After 48 h, the cells were washed with PBS and then cross-linked with UV light for 15 min. Then the cells were lysed in Pierce IP Lysis buffer containing protease, RNA inhibitors, and DNase I for 30 min on ice, and centrifuged at 12,000 g for 10 min. ANTI-FLAG^® M2 magnetic beads were washed with TBS solution, and the cell lysates were added and incubated overnight at 4°C. The FLAG-tagged proteins were immunoprecipitated from the supernatants with ANTI-FLAG^® M2 magnetic beads. After three washes with a low-salt wash buffer (50 mM Tris-Cl, pH 7.4, and 150 mM NaCl), the magnetic beads were boiled in loading buffer, and the protein samples were further analyzed through mass spectrometry (Changya Biotech Co., Ltd., Shanghai, China). Western blot analysis was further used to confirm the results of mass spectrometry.

Statistical Analysis

All statistical analyses were performed using GraphPad Prism v6.01 (GraphPad Software Inc., San Diego, CA, United States). Data were presented as the mean and standard error of the mean. An independent-samples t-test was used for statistical comparisons between two groups. One-way analysis of variance followed by the Student–Newman–Keuls test was used for statistical comparisons among multiple groups. A p-value less than 0.05 was deemed statistically significant.

Results

Lncprep + 96kb Is Spatiotemporally Expressed in Murine Ovaries

The expression of lncPrep + 96kb in ovaries was detected through in situ hybridization. We found that the lncPrep + 96kb signal was mainly located in the granulosa cells of primary and secondary follicles. A weak signal was also detected in the granulosa cells of antral follicles (Figure 1).

FIGURE 1

Figure 1. The expression of lncPrep + 96kb in murine ovaries. The negative control is shown in (A). (B,C) Show the positive signal of lncPrep + 96kb. Scale bar = 100 μm.

LncPrep + 96kb 2.2 kb Decreases the Production of Estradiol and Inhibits the Expression of Aromatase

To observe the function of lncPrep + 96kb 2.2 kb in granulosa cells, we overexpressed and knocked down the expression of lncPrep + 96kb 2.2 kb. The CCK-8 assay results showed that the proliferation of granulosa cells was not affected (Figure 2A). Overexpression of lncPrep + 96kb 2.2 kb inhibited the production of estradiol, and knockdown of lncPrep + 96kb 2.2 kb promoted the secretion of estradiol (Figures 2B,C). The expression of aromatase was negatively regulated by lncPrep + 96kb 2.2 kb at the mRNA and protein levels (Figures 2D–F).

FIGURE 2

Figure 2. Estradiol production and aromatase expression are regulated by lncPrep + 96kb 2.2 kb. (A) Shows the effect of lncPrep + 96kb 2.2 kb on cell proliferation. (B,C) Show estradiol production after the overexpression and knockdown of lncPrep + 96kb 2.2 kb. (D,E) Show the mRNA expression of aromatase after overexpression and knockdown of lncPrep + 96kb 2.2 kb. (F) Shows the protein expression of aromatase after overexpression and knockdown of lncPrep + 96kb 2.2 kb. The experiment was independently repeated for three times. *p < 0.05, which was a specific two-group comparison.

Lncprep + 96kb 2.2 kb Binds to EDF1 and Induces the Translocation of EDF1

To identify lncPrep + 96kb 2.2 kb interacting proteins, we conducted RIP based on the specific interaction between the MS2 coat protein and the MS2 binding sites (MS2bs) and mass spectrometry analysis. Endothelial differentiation-related factor 1 was identified as a binding partner of lncprep + 96kb 2.2 kb. Western blotting with an anti-EDF1 antibody confirmed the existence of EDF1 within lncPrep + 96kb 2.2 kb pull-down samples (Figure 3A). We also observed the effect of lncPrep + 96kb 2.2 kb on EDF1 and found that overexpression of lncPrep + 96kb 2.2 kb induced the translocation of EDF1 from the nucleus to the cytoplasm, however, knockdown of lncPrep + 96kb 2.2 kb had no influence on the location of EDF1 (Figure 3B).

FIGURE 3

Figure 3. The effect of lncPrep + 96kb 2.2 kb on EDF1. (A) Western blot was used to determine the specific interaction of lncPrep + 96kb 2.2 kb with EDF1. (B) Immunofluorescence was used to examine the location of EDF1 in granulosa cells. Scale bar = 100 μm.

EDF1-Specific shRNA Downregulates the Expression of Aromatase and Inhibits the Secretion of Estradiol

To confirm the function of EDF1 in murine ovaries, we knocked down the expression of EDF1 by EDF1-specific shRNA. After knockdown of EDF1, estradiol production and aromatase expression were decreased, suggesting that EDF1 promoted the secretion of estradiol (Figure 4).

FIGURE 4

Figure 4. Estradiol production and aromatase expression are regulated by EDF1. (A) Shows the effect of EDF1-specific shRNA on estradiol production. (B–D) Shows the effect of EDF1-specific shRNA on aromatase expression at the mRNA and protein levels. The experiment was independently repeated for three times. *p < 0.05, which was a specific two-group comparison.

The 1,979–2,077 Site of lncPrep + 96kb 2.2 kb Contributes to the Binding of EDF1

To prove the effect of lncPrep + 96kb 2.2 kb on aromatase and estradiol through EDF1, we first evaluated the possible interaction between lncPrep + 96kb 2.2 kb and EDF1 through bioinformatics approaches. The results from the determination of the CatRAPID signature revealed that the overall EDF1/RNA interaction score was 0.75 (Figure 5A). Then, the CatRAPID fragment algorithm, which is based on the interactions between individual polypeptides and nucleotide fragments, was used for analysis. As shown in Figure 5B, the 1,979–2,077 and 603–690 nucleotide positions of lncPrep + 96kb 2.2 kb have high capabilities to bind to EDF1. Then, we deleted these binding sites and found that mutation of 1,979–2,077 in lncPrep + 96kb 2.2 kb rescued the effects of lncPrep + 96kb 2.2 kb on estradiol production and aromatase expression. However, the 603–690 site did not contribute to the effect of lncPrep + 96kb 2.2 kb (Figures 6A–D).

FIGURE 5

Figure 5. Prediction of the interaction between lncPrep + 96kb 2.2 kb and EDF1. (A) The RNA-binding capability of the EDF1 protein was predicted by the CatRAPID signature module. Overall interaction scores greater than 50% indicate binding capability. (B) CatRAPID fragment module prediction of the interaction profile and matrix for lncPrep + 96kb 2.2 kb and EDF1.

FIGURE 6

Figure 6. The effects of the mutating the lncPrep + 96kb 2.2 kb binding sites on estradiol production and aromatase expression. (A) Shows estradiol production after mutating the lncPrep + 96kb 2.2 kb binding sites. (B) Shows the expression of aromatase after mutating the lncPrep + 96kb 2.2 kb binding sites. (C,D) Show the protein expression of aromatase after mutating the lncPrep + 96kb 2.2 kb binding sites. *p < 0.05. The experiment was independently repeated for three times.

Discussion

Spatially and temporally specific expression of lncRNAs has a strong association with the cell differentiation process (Cesana et al., 2011; Wang Y. et al., 2013). LncPrep + 96kb is a novel lncRNA that is highly expressed in murine granulosa cells. We examined the spatially and temporally specific expression of lncPrep + 96kb through in situ hybridization. The results showed that lncPrep + 96kb was highly and specifically expressed in the granulosa cells of primary follicles and secondary follicles, and a low expression level was identified in antral follicles.

Early growth in follicular development is accompanied by an increase in oocyte diameter and the proliferation of granulosa cells. Early follicular development requires the appropriate expression of numerous genes at different developmental stages (Bonnet et al., 2013, 2015; Ernst et al., 2018). LncPrep + 96kb was specifically expressed in granulosa cells of early stage follicles, so we inferred that it may influence the proliferation of granulosa cells. However, the results did not support this hypothesis. Compared with the granulosa cells of early growing follicles, the granulosa cells of antral follicles express different differentiation markers, such as FSH receptor and aromatase, and have the capacity to secrete estradiol. We cultured granulosa cells and observed the effect of lncPrep + 96kb 2.2 kb on estradiol production. We found that lncPrep + 96kb 2.2 kb inhibited estradiol production and aromatase expression. These results suggest that high expression of lncPrep + 96kb 2.2 kb may maintain the granulosa cells of early stage follicles in an undifferentiated state.

As the key enzyme of estradiol synthesis, aromatase controls the production of estradiol. In adult rats, aromatase is no longer expressed by small growing follicles, but is confined to healthy large antral follicles and preovulatory follicles (Guigon et al., 2003; Sakurada et al., 2006; Stocco, 2008; Dewailly et al., 2016). The spatially and temporally specific expression of aromatase was strictly controlled. The proximal ends of aromatase promoters contain a number of fairly well-conserved regulatory elements. The ovary-specific PII promoter binds to cAMP-response element-binding protein (CREB), LRH1/NR5A2, and SF1/NR5A1 (Stocco, 2008; Lai et al., 2014). LRH1 is primarily expressed in granulosa cells and luteal cells, whereas SF1 is expressed in theca/interstitial cells and granulosa cells (Falender et al., 2003; Hinshelwood et al., 2003; Liu et al., 2003; Lai et al., 2014). Therefore, in granulosa cells, SF1 and LRH1 contribute to the regulation of aromatase.

As unstable biomacromolecules, lncRNAs bind to specific proteins to form RNA/protein complexes to ensure their stability in cells. In addition, the biological function of lncRNAs cannot be achieved without the reversible or irreversible binding of related protein molecules. RNA–protein interactions are the basis of many cellular biological processes, such as transcription, posttranscriptional regulation, protein translation, RNA stability, and RNA transport and localization. Long non-coding RNAs can interact with proteins in a variety of forms to serve as protein decoys and protein scaffolds to perform functions related to protein targeting and intracellular signaling (Wang and Chang, 2011). We employed RIP based on the specific interaction between the MS2 coat protein and MS2 binding sites (MS2bs) to identify the proteins that bind to lncPrep + 96kb 2.2 kb. Mass spectrometry and Western blotting identified EDF1, which can bind to lncPrep + 96kb 2.2 kb. Further study showed that overexpression of lncPrep + 96kb 2.2 kb induced the translocation of EDF1 from the nucleus to the cytoplasm.

Endothelial differentiation-related factor 1, a highly conserved intracellular protein, was initially identified as an evolutionarily conserved transcriptional coactivator that links regulatory factors and the TATA element-binding protein (Takemaru et al., 1997). Studies have shown that EDF1 can act as a bridging factor, enabling the interaction of SF1, LRH1, liver-X-receptor-α (LXR-α), and peroxisome proliferator activated receptor-γ with the transcription machinery (Brendel et al., 2002; Leidi et al., 2009). As a transcriptional coactivator, EDF1 functions in the nucleus of cells and interacts with TATA box-binding proteins. Some studies also showed that EDF1 was localized in the cytosol, where it binds to calmodulin and regulates CaM availability (Ballabio et al., 2004; Cazzaniga et al., 2018). Therefore, EDF1 serves different functions depending on the location of EDF1 in cells. Recent studies have shown that the lncRNA Blnc1 can bind to EDF1 and contribute to LXR transcriptional activation. According to the RIP results, we also showed that EDF1 can bind to lncPrep + 96kb 2.2 kb, and lncPrep + 96kb 2.2 kb promoted the translocation of EDF1 from the nucleus to the cytoplasm. We tested the function of EDF1 in granulosa cells and found that EDF1-specific shRNA indeed inhibited estradiol production and aromatase expression. Endothelial differentiation-related factor 1 is a coactivator of SF1 and LRH1, which are the main regulators of the aromatase promoter. Therefore, we infer that EDF1 contributes to estradiol production in granulosa cells as a coactivator. Before prepuberty, lncPrep + 96kb 2.2 kb binds to EDF1 and remains in the cytoplasm to prevent the activity of EDF1 as a coactivator.

Subsequently, the binding sites of lncPrep + 96kb 2.2 kb and EDF1 were predicted by CatRAPID. Two locations in lncPrep + 96kb 2.2 kb (1,979–2,077 and 603–690) had high scores. After mutation of these two locations, we found that the loss of 1,979–2,077 but not 603–690 reduced the effect of lncPrep + 96kb 2.2 kb on estradiol production and aromatase expression. These results suggested that the binding of lncPrep + 96kb 2.2 kb with EDF1 contributed to the localization of EDF1 and inhibited aromatase expression. We will explore the function of EDF1 in the cytoplasm in future studies.

Conclusion

Taken together, the results of this study demonstrate the function of lncPrep + 96kb 2.2 kb in granulosa cells. This lncRNA was highly expressed in granulosa cells of early growing follicles. It inhibited aromatase expression and estradiol production in granulosa cells by inducing the translocation of EDF1 from the nucleus to the cytoplasm.

Data Availability Statement

All datasets generated for this study are included in the article/Supplementary Material.

Ethics Statement

The animal study was reviewed and approved by the Ethical Committee of Nanchang University.

Author Contributions

FF, JW, and RB performed the experiments and wrote the manuscript. LL, XT, SH, HZ, WW, and LX performed the experiments and analyzed the data. CZ conceived and designed the experiments. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (Grant nos. 81601242 and 81960272) and the Natural Science Foundation of Jiangxi Province (Grant no. 20161BAB215199).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2020.00481/full#supplementary-material

FIGURE S1 | The purity of granulosa cells. The expression of FSHR was detected with immunofluorescence to confirm the purity of granulosa cells. We examined the FSHR signal in granulosa cells isolated from PMSG treated mice through mechanical isolation. Scale bar = 100 μm.

FIGURE S2 | The effect of lncPrep + 96kb 2.2 kb on StAR, P450scc, SF and GATA4. Western blotting was used to examine the protein level of StAR, P450scc, SF and GATA4 after overexpression and knockdown of lncPrep + 96kb 2.2 kb.

FIGURE S3 | The effect of lncPrep + 96kb 2.2 kb on estradiol production and aromatase expression in granulosa cells from untreated mice. The granulosa cells were isolated from 21 days female mice that were not treated with PMSG. After the overexpression and knockdown of lncPrep + 96kb 2.2 kb, we examined the estradiol production at 72 h and aromatase expression at 48 h. (A) shows estradiol production after the overexpression and knockdown of lncPrep + 96kb 2.2 kb. (B) shows the mRNA expression of aromatase after overexpression and knockdown of lncPrep + 96kb 2.2 kb. (C) shows the protein expression of aromatase after overexpression and knockdown of lncPrep + 96kb 2.2 kb. ^∗p < 0.05.

Footnotes

1. ^ http://crispr.mit.edu

References

Agostini, F., Zanzoni, A., Klus, P., Marchese, D., Cirillo, D., and Tartaglia, G. G. (2013). catRAPID omics: a web server for large-scale prediction of protein-RNA interactions. Bioinformatics 29, 2928–2930. doi: 10.1093/bioinformatics/btt495

PubMed Abstract | CrossRef Full Text | Google Scholar

Baba, T., Ting, A. Y., Tkachenko, O., Xu, J., and Stouffer, R. L. (2017). Direct actions of androgen, estrogen and anti-Mullerian hormone on primate secondary follicle development in the absence of FSH in vitro. Hum. Reprod. 32, 2456–2464. doi: 10.1093/humrep/dex322

PubMed Abstract | CrossRef Full Text | Google Scholar

Ballabio, E., Mariotti, M., De Benedictis, L., and Maier, J. A. (2004). The dual role of endothelial differentiation-related factor-1 in the cytosol and nucleus: modulation by protein kinase A. Cell. Mol. Life Sci. 61, 1069–1074. doi: 10.1007/s00018-004-4016-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonnet, A., Cabau, C., Bouchez, O., Sarry, J., Marsaud, N., Foissac, S., et al. (2013). An overview of gene expression dynamics during early ovarian folliculogenesis: specificity of follicular compartments and bi-directional dialog. BMC Genomics 14:904. doi: 10.1186/1471-2164-14-904

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonnet, A., Servin, B., Mulsant, P., and Mandon-Pepin, B. (2015). Spatio-temporal gene expression profiling during in vivo early ovarian folliculogenesis: integrated transcriptomic study and molecular signature of early follicular growth. PLoS One 10:e0141482. doi: 10.1371/journal.pone.0141482

PubMed Abstract | CrossRef Full Text | Google Scholar

Brendel, C., Gelman, L., and Auwerx, J. (2002). Multiprotein bridging factor-1 (MBF-1) is a cofactor for nuclear receptors that regulate lipid metabolism. Mol. Endocrinol. 16, 1367–1377. doi: 10.1210/mend.16.6.0843

PubMed Abstract | CrossRef Full Text | Google Scholar

Cazzaniga, A., Locatelli, L., Castiglioni, S., and Maier, J. (2018). The contribution of EDF1 to PPARgamma transcriptional activation in VEGF-treated human endothelial cells. Int. J. Mol. Sci. 19:1830. doi: 10.3390/ijms19071830

PubMed Abstract | CrossRef Full Text | Google Scholar

Cesana, M., Cacchiarelli, D., Legnini, I., Santini, T., Sthandier, O., Chinappi, M., et al. (2011). A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell 147, 358–369. doi: 10.1016/j.cell.2011.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Chou, C. H., and Chen, M. J. (2018). The effect of steroid hormones on ovarian follicle development. Vitam. Horm. 107, 155–175. doi: 10.1016/bs.vh.2018.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Consortium, E. P. (2012). An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74. doi: 10.1038/nature11247

PubMed Abstract | CrossRef Full Text | Google Scholar

Dewailly, D., Robin, G., Peigne, M., Decanter, C., Pigny, P., and Catteau-Jonard, S. (2016). Interactions between androgens, FSH, anti-Mullerian hormone and estradiol during folliculogenesis in the human normal and polycystic ovary. Hum. Reprod. Update 22, 709–724. doi: 10.1093/humupd/dmw027

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhahbi, J. M. (2014). Circulating small noncoding RNAs as biomarkers of aging. Ageing Res. Rev. 17, 86–98. doi: 10.1016/j.arr.2014.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Drummond, A. E. (2006). The role of steroids in follicular growth. Reprod. Biol. Endocrinol. 4:16. doi: 10.1186/1477-7827-4-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Ernst, E. H., Franks, S., Hardy, K., Villesen, P., and Lykke-Hartmann, K. (2018). Granulosa cells from human primordial and primary follicles show differential global gene expression profiles. Hum. Reprod. 33, 666–679. doi: 10.1093/humrep/dey011

PubMed Abstract | CrossRef Full Text | Google Scholar

Falender, A. E., Lanz, R., Malenfant, D., Belanger, L., and Richards, J. S. (2003). Differential expression of steroidogenic factor-1 and FTF/LRH-1 in the rodent ovary. Endocrinology 144, 3598–3610. doi: 10.1210/en.2002-0137

PubMed Abstract | CrossRef Full Text | Google Scholar

Fisher, C. R., Graves, K. H., Parlow, A. F., and Simpson, E. R. (1998). Characterization of mice deficient in aromatase (ArKO) because of targeted disruption of the cyp19 gene. Proc. Natl. Acad. Sci. U.S.A. 95, 6965–6970. doi: 10.1073/pnas.95.12.6965

PubMed Abstract | CrossRef Full Text | Google Scholar

Guigon, C. J., Mazaud, S., Forest, M. G., Brailly-Tabard, S., Coudouel, N., and Magre, S. (2003). Unaltered development of the initial follicular waves and normal pubertal onset in female rats after neonatal deletion of the follicular reserve. Endocrinology 144, 3651–3662. doi: 10.1210/en.2003-0072

PubMed Abstract | CrossRef Full Text | Google Scholar

Hinshelwood, M. M., Repa, J. J., Shelton, J. M., Richardson, J. A., Mangelsdorf, D. J., and Mendelson, C. R. (2003). Expression of LRH-1 and SF-1 in the mouse ovary: localization in different cell types correlates with differing function. Mol. Cell. Endocrinol. 207, 39–45. doi: 10.1016/s0303-7207(03)00257-0

CrossRef Full Text | Google Scholar

Hobeika, E., Armouti, M., Kala, H., Fierro, M. A., Winston, N. J., Scoccia, B., et al. (2019). Oocyte-secreted factors synergize with FSH to promote aromatase expression in primary human cumulus cells. J. Clin. Endocrinol. Metab. 104, 1667–1676. doi: 10.1210/jc.2018-01705

PubMed Abstract | CrossRef Full Text | Google Scholar

Holesh, J. E., Hazhirkarzar, B., and Lord, M. (2020). Physiology, Ovulation. Treasure Island, FL: StatPearls.

Google Scholar

Hossain, M. M., Salilew-Wondim, D., Schellander, K., and Tesfaye, D. (2012). The role of microRNAs in mammalian oocytes and embryos. Anim. Reprod. Sci. 134, 36–44. doi: 10.1016/j.anireprosci.2012.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Hutz, R. J., Dierschke, D. J., and Wolf, R. C. (1986). Markers of atresia in ovarian follicular components from rhesus monkeys treated with estradiol-17 beta. Biol. Reprod. 34, 65–70. doi: 10.1095/biolreprod34.1.65

PubMed Abstract | CrossRef Full Text | Google Scholar

Imbar, T., and Eisenberg, I. (2014). Regulatory role of microRNAs in ovarian function. Fertil. Steril. 101, 1524–1530. doi: 10.1016/j.fertnstert.2014.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, A. L. (2015). Ovarian follicle selection and granulosa cell differentiation. Poult. Sci. 94, 781–785. doi: 10.3382/ps/peu008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimura, A. P., Yoneda, R., Kurihara, M., Mayama, S., and Matsubara, S. (2017). A long noncoding RNA, lncRNA-Amhr2, plays a role in Amhr2 gene activation in mouse ovarian granulosa cells. Endocrinology 158, 4105–4121. doi: 10.1210/en.2017-00619

PubMed Abstract | CrossRef Full Text | Google Scholar

Lai, W. A., Yeh, Y. T., Fang, W. L., Wu, L. S., Harada, N., Wang, P. H., et al. (2014). Calcineurin and CRTC2 mediate FSH and TGFbeta1 upregulation of Cyp19a1 and Nr5a in ovary granulosa cells. J. Mol. Endocrinol. 53, 259–270. doi: 10.1530/JME-14-0048

PubMed Abstract | CrossRef Full Text | Google Scholar

Laven, J. S., and Fauser, B. C. (2006). What role of estrogens in ovarian stimulation. Maturitas 54, 356–362. doi: 10.1016/j.maturitas.2006.04.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Leidi, M., Mariotti, M., and Maier, J. A. (2009). Transcriptional coactivator EDF-1 is required for PPARgamma-stimulated adipogenesis. Cell. Mol. Life Sci. 66, 2733–2742. doi: 10.1007/s00018-009-0069-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G., Jiang, H., Zheng, C., Zhu, G., Xu, Y., Sheng, X., et al. (2017). Long noncoding RNA MRAK009713 is a novel regulator of neuropathic pain in rats. Pain 158, 2042–2052. doi: 10.1097/j.pain.0000000000001013

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Q., He, H., Zhang, Y. L., Li, X. M., Guo, X., Huo, R., et al. (2013). Phosphoinositide 3-kinase p110delta mediates estrogen- and FSH-stimulated ovarian follicle growth. Mol. Endocrinol. 27, 1468–1482. doi: 10.1210/me.2013-1082

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Fang, Y., Liu, Y., and Yang, X. (2015). MicroRNAs in ovarian function and disorders. J. Ovarian Res. 8:51. doi: 10.1186/s13048-015-0162-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Wang, H., Zhou, D., Shuang, T., Zhao, H., and Chen, B. (2018). Up-regulation of long noncoding RNA SRA promotes cell growth, inhibits cell apoptosis, and induces secretion of estradiol and progesterone in ovarian granular cells of mice. Med. Sci. Monit. 24, 2384–2390. doi: 10.12659/msm.907138

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, D. L., Liu, W. Z., Li, Q. L., Wang, H. M., Qian, D., Treuter, E., et al. (2003). Expression and functional analysis of liver receptor homologue 1 as a potential steroidogenic factor in rat ovary. Biol. Reprod. 69, 508–517. doi: 10.1095/biolreprod.102.011767

PubMed Abstract | CrossRef Full Text | Google Scholar

Livi, C. M., Klus, P., Delli Ponti, R., and Tartaglia, G. G. (2016). catRAPID signature: identification of ribonucleoproteins and RNA-binding regions. Bioinformatics 32, 773–775. doi: 10.1093/bioinformatics/btv629

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsubara, S., Kurihara, M., and Kimura, A. P. (2014). A long non-coding RNA transcribed from conserved non-coding sequences contributes to the mouse prolyl oligopeptidase gene activation. J. Biochem. 155, 243–256. doi: 10.1093/jb/mvt113

PubMed Abstract | CrossRef Full Text | Google Scholar

Mercer, T. R., and Mattick, J. S. (2013). Structure and function of long noncoding RNAs in epigenetic regulation. Nat. Struct. Mol. Biol. 20, 300–307. doi: 10.1038/nsmb.2480

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakagawa, S., Shimada, M., Yanaka, K., Mito, M., Arai, T., Takahashi, E., et al. (2014). The lncRNA Neat1 is required for corpus luteum formation and the establishment of pregnancy in a subpopulation of mice. Development 141, 4618–4627. doi: 10.1242/dev.110544

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, K. C., Frith, M. C., and Mattick, J. S. (2006). Rapid evolution of noncoding RNAs: lack of conservation does not mean lack of function. Trends Genet. 22, 1–5. doi: 10.1016/j.tig.2005.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosenfeld, C. S., Wagner, J. S., Roberts, R. M., and Lubahn, D. B. (2001). Intraovarian actions of oestrogen. Reproduction 122, 215–226. doi: 10.1530/rep.0.1220215

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahmi, F., Nicola, E. S., Zamberlam, G. O., Goncalves, P. D., Vanselow, J., and Price, C. A. (2014). Factors regulating the bovine, caprine, rat and human ovarian aromatase promoters in a bovine granulosa cell model. Gen. Comp. Endocrinol. 200, 10–17. doi: 10.1016/j.ygcen.2014.02.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakurada, Y., Shirota, M., Inoue, K., Uchida, N., and Shirota, K. (2006). New approach to in situ quantification of ovarian gene expression in rat using a laser microdissection technique: relationship between follicle types and regulation of inhibin-alpha and cytochrome P450aromatase genes in the rat ovary. Histochem. Cell Biol. 126, 735–741. doi: 10.1007/s00418-006-0205-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Sasson, R., Rimon, E., Dantes, A., Cohen, T., Shinder, V., Land-Bracha, A., et al. (2004). Gonadotrophin-induced gene regulation in human granulosa cells obtained from IVF patients. Modulation of steroidogenic genes, cytoskeletal genes and genes coding for apoptotic signalling and protein kinases. Mol. Hum. Reprod. 10, 299–311. doi: 10.1093/molehr/gah041

PubMed Abstract | CrossRef Full Text | Google Scholar

Shang, Z., Yu, J., Sun, L., Tian, J., Zhu, S., Zhang, B., et al. (2019). LncRNA PCAT1 activates AKT and NF-kappaB signaling in castration-resistant prostate cancer by regulating the PHLPP/FKBP51/IKKalpha complex. Nucleic Acids Res. 47, 4211–4225. doi: 10.1093/nar/gkz108

PubMed Abstract | CrossRef Full Text | Google Scholar

Stocco, C. (2008). Aromatase expression in the ovary: hormonal and molecular regulation. Steroids 73, 473–487. doi: 10.1016/j.steroids.2008.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Struhl, K. (2007). Transcriptional noise and the fidelity of initiation by RNA polymerase II. Nat. Struct. Mol. Biol. 14, 103–105. doi: 10.1038/nsmb0207-103

PubMed Abstract | CrossRef Full Text | Google Scholar

Takemaru, K., Li, F. Q., Ueda, H., and Hirose, S. (1997). Multiprotein bridging factor 1 (MBF1) is an evolutionarily conserved transcriptional coactivator that connects a regulatory factor and TATA element-binding protein. Proc. Natl. Acad. Sci. U.S.A. 94, 7251–7256. doi: 10.1073/pnas.94.14.7251

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarumi, W., Itoh, M. T., and Suzuki, N. (2014). Effects of 5alpha-dihydrotestosterone and 17beta-estradiol on the mouse ovarian follicle development and oocyte maturation. PLoS One 9:e99423. doi: 10.1371/journal.pone.0099423

PubMed Abstract | CrossRef Full Text | Google Scholar

Ting, A. Y., Xu, J., and Stouffer, R. L. (2015). Differential effects of estrogen and progesterone on development of primate secondary follicles in a steroid-depleted milieu in vitro. Hum. Reprod. 30, 1907–1917. doi: 10.1093/humrep/dev119

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K. C., and Chang, H. Y. (2011). Molecular mechanisms of long noncoding RNAs. Mol. Cell 43, 904–914. doi: 10.1016/j.molcel.2011.08.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Q., Leader, A., and Tsang, B. K. (2013). Inhibitory roles of prohibitin and chemerin in FSH-induced rat granulosa cell steroidogenesis. Endocrinology 154, 956–967. doi: 10.1210/en.2012-1836

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Xu, Z., Jiang, J., Xu, C., Kang, J., Xiao, L., et al. (2013). Endogenous miRNA sponge lincRNA-RoR regulates Oct4, Nanog, and Sox2 in human embryonic stem cell self-renewal. Dev. Cell 25, 69–80. doi: 10.1016/j.devcel.2013.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, B., Shi, Y., Gong, X., Yu, L., Chen, Q., Wang, J., et al. (2015). Evaluation of follicular synchronization caused by estrogen administration and its reproductive outcome. PLoS One 10:e0127595. doi: 10.1371/journal.pone.0127595

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, L., Froberg, J. E., and Lee, J. T. (2014). Long noncoding RNAs: fresh perspectives into the RNA world. Trends Biochem. Sci. 39, 35–43. doi: 10.1016/j.tibs.2013.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C. P., Yang, J. L., Zhang, J., Li, L., Huang, L., Ji, S. Y., et al. (2011). Notch signaling is involved in ovarian follicle development by regulating granulosa cell proliferation. Endocrinology 152, 2437–2447. doi: 10.1210/en.2010-1182

PubMed Abstract | CrossRef Full Text | Google Scholar