LncRNA-m18as1竞争性结合miR-18a-5p通过Smad2/3通路调控大鼠原代腺垂体细胞FSH分泌

LncRNA-m18as1 competitively binds with miR-18a-5p to regulate follicle-stimulating hormone secretion through the Smad2/3 pathway in rat primary pituitary cells LncRNA-m18as1竞争性结合miR-18a-5p通过Smad2/3通路调控大鼠原代腺垂体细胞FSH分泌 1 equal first-author https://orcid.org/0000-0002-6047-4916 2 equal https://orcid.org/0000-0001-7052-4967 1 1 1 1 1 1 1 corresp Wei GAO, gaowei81@jlu.edu.cn Wei GAO, gaowei81@jlu.edu.cn gaowei81@jlu.edu.cn https://orcid.org/0000-0001-8093-0555 2 corresp corresp Wei GAO, gaowei81@jlu.edu.cn Wei GAO, gaowei81@jlu.edu.cn Bao YUAN, yuan_bao@jlu.edu.cn Bao YUAN, yuan_bao@jlu.edu.cn yuan_bao@jlu.edu.cn https://orcid.org/0000-0003-3490-0755 Department of Laboratory Animals, College of Animal Sciences, Jilin University, Changchun 130062, China Department of Laboratory Animals, College of Animal Sciences, Jilin University, Changchun 130062, China 吉林大学动物科学学院，实验动物系，中国长春市，130062 吉林大学动物科学学院，实验动物系，中国长春市，130062 Jilin Provincial Model Animal Engineering Research Center, Jilin University, Changchun 130062, China Jilin Provincial Model Animal Engineering Research Center, Jilin University, Changchun 130062, China 吉林大学吉林省模型动物工程研究中心，中国长春市，130062 吉林大学吉林省模型动物工程研究中心，中国长春市，130062 Long noncoding RNAs (lncRNAs) are expressed in different species and different tissues, and perform different functions, but little is known about their involvement in the synthesis or secretion of follicle-stimulating hormone (FSH). In general, we have revealed lncRNA‍‒‍microRNA (miRNA)‍‒‍‍messenger RNA (mRNA) interactions that may play important roles in rat primary pituitary cells. In this study, a new lncRNA was identified for the first time. First, we analyzed the gene expression of lncRNA-m18as1 in different tissues and different stages by reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and observed the localization of lncRNA-m18as1 with fluorescence in situ hybridization, which indicated that this lncRNA was distributed mainly in the cytoplasm. Next, we used RT-qPCR and enzyme-linked immunosorbent assay (ELISA) to analyze the regulation of FSH synthesis and secretion after overexpression or knockdown of lncRNA-m18as1 and found that lncRNA-m18as1 was positively correlated with FSH synthesis and secretion. In addition, mothers against decapentaplegic homolog 2 (Smad2) was highly expressed in our sequencing results. We also screened miR-18a-5p from our sequencing results as a miRNA that may bind to lncRNA-m18as1 and Smad2. We used RNA immunoprecipitation-qPCR (RIP-qPCR) and/or dual luciferase assays to confirm that lncRNA-m18as1 interacted with miR-18a-5p and miR-18a-5p interacted with Smad2. Fluorescence in situ hybridization (FISH) showed that lncRNA-m18as1 and miR-18a-5p were localized mainly in the cytoplasm. Finally, we determined the relationship among lncRNA-m18as1, miR-18a-5p, and the Smad2/3 pathway. Overall, we found that lncRNA-m18as1 acts as a molecular sponge of miR-18a-5p to regulate the synthesis and secretion of FSH through the Smad2/3 pathway. 目的： 2 探索长链非编码RNA（lncRNA）-微小RNA（miRNA）-信使RNA（mRNA）网络对促卵泡素（FSH）合成与分泌的作用，为调控FSH的分子机制提供理论基础。 创新点： 2 筛选出了一个新的lncRNA作为研究对象；首次鉴定并分析了该lncRNA对FSH合成与分泌的调控作用；确定了该lncRNA调控FSH合成与分泌的机制。 方法： 2 我们通过逆转录定量聚合酶链反应（RT-qPCR）筛选出了一个新的lncRNA，并根据功能将其最终命名为lncRNA-m18as1。经过敲降或过表达lncRNA-m18as1后，我们采用RT-qPCR与酶联免疫吸附剂测定（ELISA）分析了lncRNA-m18as1对 Fshβ mRNA以及FSH分泌的调控作用。我们预测并确定了lncRNA-m18as1发挥作用的lncRNAm18as1/miR-18a-5p/Smad2轴。我们使用RNA结合蛋白免疫沉淀测定-逆转录定量聚合酶链反应（RIP-qPCR）和/或双荧光素酶报告分析方法分析了miR-18a-5p与lncRNA-m18as1、Smad2的靶向关系。此外，我们使用RT-qPCR和蛋白质印迹法（western blot）分析了lncRNA-m18as1和miR-18a-5p在轴中对向下游因子的调控作用，同时通过荧光原位杂交技术（FISH）观察了lncRNA-m18as1以及miR-18a-5p在细胞核与细胞质中的分布。 结论： 2 研究结果表明在大鼠腺垂体中高表达且阶段性表达的lncRNA-m18as1促进 Fshβ mRNA的表达与FSH的分泌。进一步的分子机制研究表明lncRNA-m18as1竞争性结合miR-18a-5p调控Smad2蛋白的表达，进而调控FSH的合成与分泌。 Long noncoding RNA (lncRNA) MicroRNA (miRNA) Competitive endogenous RNA (ceRNA) Follicle-stimulating hormone (FSH) Mothers against decapentaplegic homolog 2/3 (Smad2/3) 长链非编码RNA（lncRNA） 微小RNA（miRNA） 竞争性内源RNA（ceRNA） 促卵泡素（FSH） Smad家族成员2/3（Smad2/3） 1 Introduction 1 Follicle-stimulating hormone (FSH) is a key hormone in the regulation of mammalian reproduction, and is synthesized and secreted by the pituitary gland ( Acevedo-Rodriguez et al., 2018 Acevedo-Rodriguez et al., 2018 ). FSH is a heterodimeric glycoprotein hormone. In addition to FSH, the glycoprotein family also includes luteinizing hormone (LH), thyroid-stimulating hormone (TSH), and chorionic gonadotropin (CG) ( Pierce and Parsons, 1981 Pierce and Parsons, 1981 ). As an intermediate factor of the hypothalamic‍‍–‍‍pituitary‍‍–‍‍gonadal (HPG) axis, FSH regulates downstream spermatogenesis and oogenesis by binding to FSH receptors located on testicular supporting cells and ovarian granulosa cells ( Heckert and Griswold, 1993 Heckert and Griswold, 1993 ; George et al., 2011 George et al., 2011 ; Kaprara and Huhtaniemi, 2018 Kaprara and Huhtaniemi, 2018 ; Kumar, 2018 Kumar, 2018 ). Secretion of FSH is regulated by gonadotropin-releasing hormone (GnRH) synthesized by the hypothalamus as well as other factors, such as single-nucleotide polymorphisms (SNPs), follistatin, progesterone, activin, and microRNAs (miRNAs) ( Popovics et al., 2011 Popovics et al., 2011 ; Han et al., 2017b Han et al., 2017b ; Schubert et al., 2019 Schubert et al., 2019 ; Stamatiades et al., 2019 Stamatiades et al., 2019 ; Motta et al., 2020 Motta et al., 2020 ). MiRNAs represent a class of small biological molecules with a length of 19–23 nucleotides ( Lei et al., 2017 Lei et al., 2017 ). There are reports in the literature that miR-7a2 acts as a positive regulator of FSH and LH synthesis ( Ahmed et al., 2017 Ahmed et al., 2017 ) and that miR-‍132/212 mediates GnRH to regulate FSH expression ( Lannes et al., 2016 Lannes et al., 2016 ). In addition, several miRNAs, such as miR-21-3p and miR-433, downregulate expression of the Fshβ gene in rat anterior pituitary cells ( Han et al., 2017b Han et al., 2017b ). However, to date, our understanding of the mechanisms regulating FSH secretion is still limited. Long noncoding RNAs (lncRNAs), which represent a recently discovered type of regulatory RNA, are noncoding RNAs (ncRNAs) with a length of more than 200 nucleotides that can still function through a variety of mechanisms even though they cannot be translated into proteins ( Ma et al., 2013 Ma et al., 2013 ; Cipriano and Ballarino, 2018 Cipriano and Ballarino, 2018 ). For instance, lncRNAs regulate the stability and/or translation of target messenger RNAs (mRNAs) to play roles at the posttranscriptional level ( Nohetal., 2018 Noh et al., 2018 ). Some examples include the lncRNA β‍-‍site amyloid precursor protein-cleaving enzyme 1 (BACE1) antisense (BACE1-AS), which can enhance the stability of BACE1 mRNA, and the lncRNA growth arrest-specific 5 (GAS5), which binds to c-Myc mRNA to regulate c-Myc translation through lncRNA–mRNA interaction ( Faghihietal., 2008 Faghihi et al., 2008 ; Huet al., 2014 Hu et al., 2014 ). LncRNAs can also act as miRNA "sponges" or competitive endogenous RNAs (ceRNAs) to regulate target miRNA genes ( Salmenaet al., 2011 Salmena et al., 2011 ; Xu et al., 2020 Xu et al., 2020 ). For example, lncRNA-RNA component of mitochondrial RNA-processing endoribonuclease (RMRP) regulates miR-206 as a molecular sponge to promote proliferation, migration, and invasion in bladder cancer (BC) cell lines ( Caoet al., 2019 Cao et al., 2019 ); the lncRNA plasmacytoma variant translocation 1 (PVT1) serves as a ceRNA to negatively regulate the expression of miR-424 in cervical cancer cells ( Gaoet al., 2017 Gao et al., 2017 ); and lncRNA-maternally expressed gene 3 (MEG3) can upregulate B-cell lymphoma 2 (Bcl-2) expression by competitively binding with miR-181 family members to regulate the occurrence of gastric cancer ( Peng et al., 2015 Peng et al., 2015 ). However, there have been relatively few studies on interactions between lncRNAs and on ceRNA mechanisms in the context of animal reproduction regulation. Mothers against decapentaplegic homologs (Smads), which are eukaryotic transcription factors, play significant roles in the transforming growth factor-β (TGF‍‍-‍‍β) signaling cascade ( Massagueet al., 2005 Massague et al., 2005 ). In addition to regulating the transcription factor TGF‍-β, Smad proteins are affected by some ncRNAs. For instance, the lncRNA definitive endoderm-associated lncRNA 1 (DEANR1) induces activation of the endoderm factor forkhead box A2 (FOXA2) by interacting with Smad2/3 and positively regulating FOXA2 ( Jianget al., 2015 Jiang et al., 2015 ). The lncRNA LINC00458 directly interacts with Smad2/3 to regulate the endothelial formation process ( Chenet al., 2020 Chen et al., 2020 ). Overexpression of miR-155-5p can inhibit development of vascular calcification by inhibiting activation of the Smad2/3 signaling pathway ( Zhaoet al., 2020 Zhao et al., 2020 ). The lncRNA oxidative stress responsive serine-rich 1-antisense RNA 1 (OSER1-AS1), a molecular sponge of miR-433-3p, promotes the occurrence and development of non-small cell lung cancer by increasing Smad2 expression ( Liuet al., 2020 Liu et al., 2020 ). A role for Smad2 in regulating reproduction and development has also been reported ( Bernard, 2004 Bernard, 2004 ). However, there have been few studies on the mechanisms by which ceRNAs regulate FSH synthesis and secretion through the Smad2/3 pathway. In this study, we identified the new lncRNA MSTRG.‍37997, which we have named lncRNA-miR-18a-5p sponge one (lncRNA-m18as1). We explored whether lncRNA-m18sa1, as a molecular sponge of miR-18a-5p, regulates FSH expression through the Smad2/3 pathway. 2 Materials and methods 1 2.1 Animals and cell culture 2 All Sprague-Dawley rats used in this experiment were of specific pathogen-free (SPF) grade and were obtained from Liaoning Changsheng Biotechnology Co., Ltd. (Benxi, Liaoning, China). Samples of different tissues (heart, liver, spleen, lung, kidney, brain, pituitary, pineal, muscle, and testis) were taken from three 8-week-old male rats to detect gene expression of lncRNA-m18as1. Pituitary samples were taken from three male rats in different age groups (2, 4, 6, and 8 weeks old) to analyze gene expression of lncRNA-m18as1. Rat primary pituitary cells were obtained as described in our previous study ( Han et al., 2017b Han et al., 2017b ). The rat primary pituitary cells were cultured in Dulbecco's modified Eagle's medium/F12 (DMEM/F12; HyClone, Logan, UT, USA) containing 100 U/mL penicillin/streptomycin (HyClone) and 10% fetal bovine serum (FBS; Gibco, Grand Island, NY, USA) in a 37 ℃ humidified incubator (Lishen, Shanghai, China) with 5% CO 2 . The rat pituitary adenoma cell lines GH3 and MMQ were obtained from the National Infrastructure of Cell Line Resource (resource numbers: 3111C0001CCC000008 and 3111C0001CCC000081) and cultured under the same conditions. 2.2 siRNA synthesis and plasmid construction 2 LncRNA-m18as1 small interfering RNAs (siRNAs) and the negative control siRNA were designed and synthesized by GeneCreate Biological Engineering Co., Ltd. (Wuhan, China). LncRNA-m18as1 complementary DNA (cDNA) was inserted into the pReceiver plasmid (pcDNA3.1) to construct a stable expression plasmid carrying lncRNA-m18as1 (this plasmid was also synthesized by GeneCreate Biological Engineering Co., Ltd.). Plasmid extraction was conducted according to the protocol of an EndoFree Maxi Plasmid kit, which was obtained from Tiangen Biotech Co., Ltd. (Beijing, China). 2.3 Synthesis of miRNA mimic 2 The miR-18a-5p mimic and negative control were synthesized by RiboBio Biotech Co., Ltd. (Guangzhou, China). 2.4 Cell transfection 2 Rat primary pituitary cells were seeded on a 24-well plate at a density of 3×10 5 cells/well. A Lipofectamine 2000 Transfection kit was purchased from Invitrogen (Carlsbad, California, USA), and the recommended protocol was used to perform all transfection experiments in this study. The treatment concentrations of miR-18a-5p mimic, miR-18a-5p negative control, lncRNA-m18as1 siRNA, and siRNA negative control were 100 nmol/L. The final concentrations of the lncRNA-m18as1 overexpression plasmid were 500, 1000, and 1500 ng/mL. Serum-free basal medium was used throughout the transfection process. After transfection, the cells were incubated for 48 h for gene expression analysis. 2.5 Fluorescence in situ hybridization 2 Rat primary pituitary cells were seeded on a 24-well plate with cell slides at a density of 3×10 5 cells/well. Specific probes for lncRNA-m18as1 and miR-18a-5p were synthesized by Suzhou Genepharma Co., Ltd. (Suzhou, China). Fluorescence in situ hybridization (FISH) was performed in accordance with the manufacturer's recommended instructions to determine the localization of lncRNA-m18as1 and miR-18a-5p. Fluorescence microscopy (Nicon, Tokyo, Japan) was used to capture images. 2.6 Immunofluorescence analysis 2 Rat primary pituitary cells were seeded on a 24-well plate with cell slides at a density of 3×10 5 cells/well to determine the nuclear and cytoplasmic distribution of phosphorylated Smad4 after treatment. Notably, the whole process had to be carried out in a dark and humid box at room temperature to prevent the cells from drying out and to prevent fluorescence quenching. Fluorescence microscopy (Nicon) was used to capture images. 2.7 RNA immunoprecipitation-qPCR (RIP-qPCR) and luciferase reporter analysis 2 First, rat primary pituitary cells were transfected with the miR-‍‍18a-5p mimic and its negative control, and a lysate was prepared with RIP Lysis Buffer. Next, the cell lysate was centrifuged, and the supernatant was divided into two groups. An anti-argonaute 2 (anti-AGO2) antibody or a nonspecific immunoglobulin G (IgG) antibody (Abcam, Cambridge, MA, UK) was added, and the samples were incubated overnight at 4 ℃. Then, protein A/G magnetic beads were added, and the samples were incubated for 1 h at 4 ℃. After that, the pellets were centrifuged and washed using RIP Wash Buffer. Enrichment of lncRNA-m18as1 and miR-18a was analyzed by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). More detailed experimental steps were carried out in accordance with the provided protocol of a Magna RIP RNA Binding Protein Immunoprecipitation kit (Millipore, Billerica, MA, USA). We used RNAhybrid 2.2 ( https://bibiserv.cebitec.uni-bielefeld. https://bibiserv.cebitec.uni-bielefeld. ‍‍de/rnahybrid) and TargetScan ( http://www. http://www. ‍‍targetscan.‍‍org/vert_71) to predict the potential binding site between lncRNA-m18as1 and miR-18a-5p. A sequence from approximately 200 base pairs upstream to 200 base pairs downstream of the lncRNA-m18as1 binding site for miR-18a-5p was cloned into the pmirGLO plasmid to obtain lncRNA-m18as1-WT; and a sequence with the lncRNA-m18a-5p binding site for miR‍-‍‍18a-5p was mutated to obtain lncRNA-m18as1-MuT. After cotransfection with miR-18a-5p mimic or its negative control, four experimental treatment groups were formed: the lncRNA-m18as1-WT+miR-control group, lncRNA-m18as1-WT+miR-18a-5p group, lncRNA-m18as1-MuT+miR-control group, and lncRNA-m18as1-MuT+miR-18a-5p group. The same method was used to predict the binding site between miR-18a-5p and Smad2 mRNA and obtain Smad2-WT and Smad2-MuT. After cotransfection with miR-18a-5p mimic or its negative control, four experimental treatment groups were formed: Smad2-WT+miR-control group, Smad2-WT+miR-18a-5p group, Smad2-MuT+miR-control group, and Smad2-MuT+miR-18a-5p group. The intensities of firefly luciferin and renilla luciferin were detected with a Dual Luciferase Reporter Gene Assay kit (Beyotime, Shanghai, China). 2.8 RNA isolation and RT-qPCR 2 Total RNA was extracted using TRIzol reagent (Invitrogen) according to the corresponding protocol. Reverse transcription of RNA was performed with a FastKing gDNA Dispelling RT SuperMix kit (Tiangen Biotech Co., Ltd.). RT-qPCR was performed using SuperReal PreMix Plus (SYBR Green, Tiangen Biotech Co., Ltd.), and the primers were designed on the National Center for Biotechnology Information (NCBI) website ( https://blast. https://blast. ‍‍ncbi.‍nlm.‍nih.‍gov/Blast.‍cgi) and verified with Primer 3 website ( https://primer3. https://primer3. ‍‍ut.‍ee). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was employed as an endogenous reference for mRNAs, lncRNA and U6 were employed as endogenous references for miRNA. All primers used in this study are listed in Table S1. 2.9 ELISA 2 After 48 h of transfection, the culture supernatant was collected and the FSH secretion level was analyzed with a rat FSH enzyme-linked immunosorbent assay (ELISA) kit (Enzyme-linked Biotechnology Co., Ltd., Shanghai, China) following the manufacturer's protocol. The concentration gradient of the standard curve was 0, 0.75, 1.5, 3, 6, and 12 IU/L. 2.10 Western blot analysis 2 After 48 h of transfection, rat primary pituitary cells were lysed in radioimmunoprecipitation assay (RIPA) buffer containing 1% (volume fraction) phenylmethylsulfonyl fluoride (PMSF) and 1% (volume fraction) protein phosphatase inhibitor Mix (Beijing Solarbio Science & Technology Co., Ltd., China). The concentration of each protein sample was determined with a bicinchoninic acid (BCA) assay kit (Millipore). After separation with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE), the proteins were transferred to polyvinylidene difluoride (PVDF) membranes (Millipore). Then, the PVDF membranes were incubated in blocking solution prepared with 5% (volume fraction) bull serum albumin (BSA) dissolved in Tris-buffered saline Tween (TBST) for 1 h at room temperature and then incubated with specific primary antibodies at 4 ℃ overnight. On the second day, the membranes were washed three times in TBST for 10 min each before being incubated with secondary antibodies for 1 h at room temperature. The immunoreactive bands were detected using Millipore Sigma TM Immobilon TM Western Chemiluminescent horseradish peroxidase (HRP) Substrate (ECL) (WBKL0500, Millipore). The antibodies used in this experiment were anti-GAPDH (Cell Signaling Technology, Boston, MA, USA), anti-Smad4 (Cell Signaling Technology), anti-Smad2/3 (Cell Signaling Technology), anti-phospho-Smad2/3 (Abcam), and anti-phospho-Smad4 (Affinity, Golden, CO, USA). The primary antibody and blocking solution were diluted 1:‍‍1000 (volume ratio) for subsequent use, while the secondary antibody and blocking solution were diluted 1:5000 (volume ratio). 2.11 Statistical analysis 2 Western blot bands were grayscale-analyzed with ImageJ (Version 1.52a, National Institutes of Health, Bethesda, MD, USA) and the data from luciferase reporter analysis and western blot were normalized. Then all data were presented as mean±standard deviation (SD) of three independent experiments. Data were subjected to multiple comparisons for significant differences by SPSS 19.0 (SPSS, Chicago, IL, USA). The whole graphs were generated using GraphPad Prism 8 software (Prism for Mac OS X, GraphPad Software Inc., San Diego, California, USA). P < 0.05 was considered to indicate statistical significance. 3 Results 1 3.1 Identification of lncRNA-m18as1 2 In our previous study, we identified lncRNAs in immature and mature rat pituitary glands ( Hanet al., 2017a Han et al., 2017a ). Based on these data, we selected the lncRNA MSTRG.37997 as our candidate research object and renamed it lncRNA-m18as1. LncRNA-m18as1 is an intergenic lncRNA on chromosome 14 composed of two introns ( Fig. 1 Fig. 1 a). We used Coding Potential Calculator (CPC) 2.0 ( http://cpc2.gao-lab.org http://cpc2.gao-lab.org ) to prove that lncRNA-m18as1 did not have coding ability ( Fig. 1 Fig. 1 b). We further proved this lack of coding potential through related tools on the NCBI and University of California Santa Cruz (UCSC) websites (data not shown). We further investigated the expression of lncRNA-m18as1 in different tissues of male rats. The results suggested that lncRNA-m18as1 was indeed expressed in these tissues. The pituitary had the highest expression of lncRNA-m18as1, followed by the pineal gland, spleen, lung, testis, liver, and brain, while the heart, kidney, and muscle had similar expression that was significantly lower than the expression in other tissues ( Fig. 1 Fig. 1 c). In addition, we detected the gene expression of lncRNA-m18as1 in the pituitary tissues of rats of different ages (2, 4, 6, and 8 weeks old). When the expression of lncRNA-m18as1 in 2-week-old rats was used as a control, the expression of lncRNA-m18as1 in 4-week-old rats was decreased in comparison, and then expression gradually increased with increasing age ( Fig. 1 Fig. 1 d). Quantitative analysis of lncRNA-m18as1 in different pituitary cells showed that lncRNA-m18as1 expression was lowest in the MMQ cell line and was almost the same in primary pituitary cells and the GH3 cell line ( Fig. 1 Fig. 1 e). The FISH results indicated that lncRNA-m18as1 was located mainly in the cytoplasm ( Fig. 1 Fig. 1 f). 3.2 Effect of lncRNA-m18as1 on secretion of FSH by upregulating expression of Fshβ 2 To determine transfection efficiency, we transfected siRNAs and overexpression plasmids coupled with a fluorescent marker into pituitary cells and found that the cells were marked with red fluorescence, indicating that they were successfully transfected (Fig. S1). To determine which siRNA had the best knockdown effect, we transfected three siRNAs into rat primary pituitary cells and tested the expression of the lncRNA-m18as1 gene. The three siRNAs downregulated lncRNA-m18as1 to a similar degree, so we ultimately chose siRNA-3 for follow-up experimental treatment according to the RT-qPCR result (hereafter referred to as the siRNA) ( Fig. 2 Fig. 2 a). We transfected the lncRNA-m18as1 overexpression plasmid into rat primary pituitary cells at three different concentrations (500, 1000, and 1500 ng/mL) and screened 500 ng/mL as the optimal concentration ( Fig. 2 Fig. 2 b). To determine whether lncRNA-m18as1 regulated the Fshβ gene, we transfected cells with the lncRNA-m18as1 siRNA and overexpression plasmid and analyzed Fshβ regulation. We also detected the effects of the above treatments on Lhβ and Cga . The results suggested that knockdown of lncRNA-m18as1 suppressed the mRNA expression of Fshβ ( P < 0.01), Lhβ ( P < 0.05), and Cga ( P < 0.05) and reduced secretion of FSH ( P < 0.01) (Figs. 2c 2c and 2d 2d ). However, overexpression of lncRNA-m18as1 had the opposite effects: it increased the mRNA expression of Fshβ , Lhβ , and Cga , as well as secretion of FSH (Figs. 2e 2e and 2f 2f ). 3.3 Prediction and confirmation of the lncRNA-m18as1 target gene 2 LncRNA-m18as1 was highly expressed in our sequencing data, and Smad2 was also highly expressed. Therefore, to explore the mechanism mediated by lncRNA-m18as1 and Smad2, we used RNAhybrid 2.2 and TargetScan to predict miRNAs that might be able to bind with lncRNA-m18as1 or Smad2 in our miRNA sequencing results. Interestingly, miR-18a-5p had a base pair that matched lncRNA-m18as1 and a base pair that matched Smad2 (Figs. 3a 3a and 3b 3b ). Therefore, we hypothesized that lncRNA-miRNA-mRNA network existed, which could regulate the synthesis and secretion of FSH, namely, an lncRNA-m18as1/miR-18a-5p/Smad2 network. To confirm the interaction between lncRNA-m18as1 and miR-18a-5p, we detected the enrichment of miR-18a-5p and lncRNA-m18as1 in cells transfected with the miR-18a-5p mimic or negative control by AGO2 RIP-qPCR. In the treatment group, lncRNA-m18as1 was enriched along with miR-18a-5p ( Fig. 3 Fig. 3 e). The dual luciferase reporter assay indicated that the binding site was no longer capable of binding miR-18a-5p after mutation ( Fig. 3 Fig. 3 d). In general, these results suggested that lncRNA-m18as1 and miR-18a-5p had the ability to bind. Additionally, we observed the colocalization of lncRNA-m18as1 and miR-18a-5p in cells by FISH, which showed that these molecules were located mainly in the cytoplasm ( Fig. 3 Fig. 3 f). To confirm the interaction between miR-18a-5p and Smad2, we performed a dual luciferase reporter gene experiment. This experiment showed that fluorescence increased significantly after binding site mutation ( P < 0.05), indicating that this site was the target of miR-18a-5p ( Fig. 3 Fig. 3 c). These results suggested that lncRNA-m18as1 may have competitively bound with miR-18a-5p to regulate expression of Smad2. 3.4 Competitive binding of lncRNA-m18as1 to miR-18a-5p to regulate the protein level of Smad2 2 To analyze the relationship between this ceRNA mechanism and Smad2, we treated rat primary pituitary cells with an overexpression plasmid or siRNA of lncRNA-m18as1. The results showed that overexpression of lncRNA-m18as1 decreased expression of miR-18a-5p ( P < 0.05) and increased the protein level of Smad2 (60 kDa) ( P < 0.05), as shown in Figs. 4 Figs. 4 a‍‍‒‍‍4d. Knockdown of lncRNA-m18as1 exhibited the opposite result ( Figs. 4 Figs. 4 e‍‒‍4h). Then, we cotransfected lncRNA-m18as1 and miR-18a-5p in a rescue experiment. Overexpression of miR-18a-5p inhibited the lncRNA-m18as1 overexpression-induced increases in Smad2 protein levels ( P < 0.05; Figs. 4i 4i and 4j 4j ).These findings suggested that competitive binding of lncRNA-m18as1 to miR-18a-5p may have regulated secretion of FSH through the Smad2/3 pathway. 3.5 Regulation of FSH by lncRNA-m18as1 through the Smad2/3 pathway 2 To determine whether lncRNA-m18as1 competitively bound miR-18a-5p through the Smad2/3 pathway to regulate FSH, we transfected cells with the lncRNA-m18as1 overexpression plasmid and added the Smad2/3 pathway inhibitor SB431542 ( Giraldez et al., 2006 Giraldez et al., 2006 ; Krol et al., 2010 Krol et al., 2010 ). We found that treatment with the inhibitor SB431542 inhibited upregulation of phosphorylated Smad2 and phosphorylated Smad4 (Figs. 5a 5a and 5b 5b ), and the ELISA results demonstrated that inhibition of the Smad2/3 pathway reduced secretion of FSH ( Fig. 5 Fig. 5 c). We also observed nucleocytoplasmic distribution of phosphorylated Smad4. We found that phosphorylated Smad4 was distributed mainly in the nuclei of rat primary pituitary cells ( Fig. 5 Fig. 5 d). These findings indicated that lncRNA-m18as1 sponged miR-18a-5p to regulate FSH through the Smad2/3 pathway. 4 Discussion 1 With the deepening of research, interest in lncRNAs is growing. Although lncRNAs are usually expressed at lower levels than other RNAs, they often exhibit more significant tissue specificity ( Ransohoff et al., 2018 Ransohoff et al., 2018 ). In this study, lncRNA-m18as1, an lncRNA located on chromosome 14, had the highest expression in the pituitary gland, while its expression in the heart, kidneys, and muscles was lower. Moreover, expression levels of lncRNA-m18as1 differed in two pituitary tumor cell lines, which suggests that this lncRNA exhibits cell type-specific expression. The detected levels of lncRNA-m18as1 in rats of different ages were consistent with the concept that lncRNAs have stage specificity ( Wanget al., 2021 Wang et al., 2021 ). LncRNAs are also known to exhibit spatial specificity. It is this spatial specificity that enables lncRNAs to have different functions in the cytoplasm and the nucleus. In the nucleus, lncRNAs function as chromatin regulators to regulate transcriptional procedures and as scaffolds to participate in the establishment and spatial organization of nuclear compartments. In the cytoplasm, lncRNAs regulate the transduction and transformation of pathway signals and subsequent modification of gene expression ( Yu and Li, 2020 Yu and Li, 2020 ; Bridgeset al., 2021 Bridges et al., 2021 ). We observed through FISH that lncRNA‍-m18as1 was located mainly in the cytoplasm, suggesting that it may be involved in a mechanism regulating pathway activation or gene expression. FSH secreted by the pituitary gland plays a vital role in reproduction and development in mammals, and many related articles have recently been published ( Dewailly et al., 2016 Dewailly et al., 2016 ; Stilley and Segaloff, 2018 Stilley and Segaloff, 2018 ; Santi et al., 2020 Santi et al., 2020 ). The network mechanism that regulates the synthesis and secretion of FSH is very complicated. For example, hormones including GnRH, steroid hormones, and pituitary adenylate cyclase activating polypeptide (PACAP) play a very important role in regulating the synthesis and secretion of FSH ( Marshall and Kelch, 1986 Marshall and Kelch, 1986 ; Fujii et al., 2002 Fujii et al., 2002 ; Kanasaki et al., 2012 Kanasaki et al., 2012 ). Studies have also shown that the activator protein-1 ( Ap-1 ) gene, forkhead box L2 ( Foxl2 ) gene, and others, play a certain regulatory role in the transcriptional stage of FSH synthesis and secretion ( Coss et al., 2004 Coss et al., 2004 ; Lamba et al., 2009 Lamba et al., 2009 ). Multiple pathways including the cyclic adenosine monophosphate (cAMP) signaling pathway and mitogen-activated protein kinase (MAPK) signaling pathway are constantly enriched with more extensive research ( Garrel et al., 2010 Garrel et al., 2010 ; Tsutsumi et al., 2010 Tsutsumi et al., 2010 ; Sun et al., 2015 Sun et al., 2015 ). However, there have been few studies on the lncRNAs that regulate the synthesis and secretion of FSH. In this study, we found that lncRNA-m18as1 could promote FSH synthesis and secretion. The regulatory mechanism was not specific, because the gene expression of Lhβ and Cga was similarly regulated. A complicated network may have caused this phenomenon. MiRNAs, which are highly conserved ncRNAs composed of 19‒23 nucleotides, control many developmental and cellular processes in animals, plants, and protozoa as posttranscriptional regulators ( Krol et al., 2010 Krol et al., 2010 ). In recent years, explorations of the roles of miRNAs in embryonic development, tumor progression, metabolic diseases, and other processes have opened up new dimensions of miRNAs ( Bhaskaran and Mohan, 2014 Bhaskaran and Mohan, 2014 ). It has been reported that miRNAs, like lncRNAs, are spatially specific ( Long et al., 2018 Long et al., 2018 ). Therefore, while observing lncRNA-m18as1 by FISH, we also observed miR-18a-5p and found that miR-18a-5p was also located mainly in the cytoplasm. This finding hints that lncRNA-m18as1 and miR-18a-5p may interact. One of the most well-known mechanisms of lncRNAs is molecular sponging, in which lncRNAs bind to miRNAs to regulate downstream processes through the ceRNA mechanism ( Chan and Tay, 2018 Chan and Tay, 2018 ). For example, lncRNAs can competitively bind to miRNAs to participate in the pathological development of diseases ( Li et al., 2019 Li et al., 2019 ; Wang et al., 2019 Wang et al., 2019 ). In our study, the results showed that lncRNA-m18as1 completely bound with miR-18a-5p, which is consistent with the view that lncRNAs regulate gene expression by binding to miRNAs ( Cesana et al., 2011 Cesana et al., 2011 ). Our results demonstrated that Smad2 was the target gene of miR-18a-5p in rat primary pituitary cells, which means that lncRNA-m18as1 may completely bind with miR-18a-5p to regulate expression of Smad2. Indeed, miR-18a-5p weakened the upregulation of the total Smad2/3 protein level induced by lncRNA-m18as1, suggesting that lncRNA-m18as1 sponged miR-18a-5p to regulate the activity of the Smad2/3 signaling pathway. Then, we added SB431542, an inhibitor of the Smad2/3 pathway, together with lncRNA-m18as1 overexpression plasmid to rat primary pituitary cells. After inhibition of this pathway, the concentration of FSH in the culture supernatant was reduced, indicating that lncRNA-m18as1 and miR-18a-5p form a ceRNA network to regulate the synthesis and secretion of FSH by stabilizing the content of Smad2. Phosphorylated Smad4 was located mainly in the nucleus, and there was export from or import into the nucleus. Previous studies have also demonstrated the distribution of phosphorylated Smad4 ( Massague et al., 2005 Massague et al., 2005 ; Arjunan et al., 2016 Arjunan et al., 2016 ). In conclusion, we found for the first time that lncRNA-m18as1, as a molecular sponge of miR-18a-5p, regulates expression of FSH through the Smad2/3 pathway ( Fig. 6 Fig. 6 ). Our results add to the existing knowledge about the relationship between ceRNA mechanisms and mammalian reproduction. References Acevedo-Rodriguez A , Kauffman AS , Cherrington BD , et al. , 2018 . Emerging insights into hypothalamic-pituitary-gonadal axis regulation and interaction with stress signalling . J Neuroendocrinol , 30 ( 10 ): e12590 . https://doi.org/10.1111/jne.12590 https://doi.org/10.1111/jne.12590 Emerging insights into hypothalamic-pituitary-gonadal axis regulation and interaction with stress signalling Ahmed K , LaPierre MP , Gasser E , et al. , 2017 . Loss of microRNA-7a2 induces hypogonadotropic hypogonadism and infertility . J Clin Invest , 127 ( 3 ): 1061 - 1074 . https://doi.org/10.1172/JCI90031 https://doi.org/10.1172/JCI90031 Loss of microRNA-7a2 induces hypogonadotropic hypogonadism and infertility Arjunan P , Gnanaprakasam JP , Ananth S , et al. , 2016 . Increased retinal expression of the pro-angiogenic receptor GPR91 via BMP6 in a mouse model of juvenile hemochromatosis . Invest Ophthalmol Vis Sci , 57 ( 4 ): 1612 - 1619 . https://doi.org/10.1167/iovs.15-17437 https://doi.org/10.1167/iovs.15-17437 Increased retinal expression of the pro-angiogenic receptor GPR91 via BMP6 in a mouse model of juvenile hemochromatosis Bernard DJ , 2004 . Both SMAD2 and SMAD3 mediate activin-stimulated expression of the follicle-stimulating hormone β subunit in mouse gonadotrope cells . Mol Endocrinol , 18 ( 3 ): 606 - 623 . https://doi.org/10.1210/me.2003-0264 https://doi.org/10.1210/me.2003-0264 Both SMAD2 and SMAD3 mediate activin-stimulated expression of the follicle-stimulating hormone β subunit in mouse gonadotrope cells Bhaskaran M , Mohan M , 2014 . MicroRNAs: history, biogenesis, and their evolving role in animal development and disease . Vet Pathol , 51 ( 4 ): 759 - 774 . https://doi.org/10.1177/0300985813502820 https://doi.org/10.1177/0300985813502820 MicroRNAs: history, biogenesis, and their evolving role in animal development and disease Bridges MC , Daulagala AC , Kourtidis A , 2021 . LNCcation: lncRNA localization and function . J Cell Biol , 220 ( 2 ): e202009045 . https://doi.org/10.1083/jcb.202009045 https://doi.org/10.1083/jcb.202009045 LNCcation: lncRNA localization and function Cao HL , Liu ZJ , Huang PL , et al. , 2019 . lncRNA-RMRP promotes proliferation, migration and invasion of bladder cancer via miR-206 . Eur Rev Med Pharmacol Sci , 23 ( 3 ): 1012 - 1021 . https://doi.org/10.26355/eurrev_201902_16988 https://doi.org/10.26355/eurrev_201902_16988 lncRNA-RMRP promotes proliferation, migration and invasion of bladder cancer via miR-206 Cesana M , Cacchiarelli D , Legnini I , et al. , 2011 . A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA . Cell , 147 ( 2 ): 358 - 369 . https://doi.org/10.1016/j.cell.2011.09.028 https://doi.org/10.1016/j.cell.2011.09.028 A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA Chan JJ , Tay Y , 2018 . Noncoding RNA: RNA regulatory networks in cancer . Int J Mol Sci , 19 ( 5 ): 1310 . https://doi.org/10.3390/ijms19051310 https://doi.org/10.3390/ijms19051310 Noncoding RNA: RNA regulatory networks in cancer Chen YF , Li YJ , Chou CH , et al. , 2020 . Control of matrix stiffness promotes endodermal lineage specification by regulating SMAD2/3 via lncRNA LINC00458 . Sci Adv , 6 ( 6 ): eaay0264 . https://doi.org/10.1126/sciadv.aay0264 https://doi.org/10.1126/sciadv.aay0264 Control of matrix stiffness promotes endodermal lineage specification by regulating SMAD2/3 via lncRNA LINC00458 Cipriano A , Ballarino M , 2018 . The ever-evolving concept of the gene: the use of RNA/protein experimental techniques to understand genome functions . Front Mol Biosci , 5 : 20 . https://doi.org/10.3389/fmolb.2018.00020 https://doi.org/10.3389/fmolb.2018.00020 The ever-evolving concept of the gene: the use of RNA/protein experimental techniques to understand genome functions Coss D , Jacobs SBR , Bender CE , et al. , 2004 . A novel AP-1 site is critical for maximal induction of the follicle-stimulating hormone β gene by gonadotropin-releasing hormone . J Biol Chem , 279 ( 1 ): 152 - 162 . https://doi.org/10.1074/jbc.M304697200 https://doi.org/10.1074/jbc.M304697200 A novel AP-1 site is critical for maximal induction of the follicle-stimulating hormone β gene by gonadotropin-releasing hormone Dewailly D , Robin G , Peigne M , et al. , 2016 . Interactions between androgens, FSH, anti-Müllerian hormone and estradiol during folliculogenesis in the human normal and polycystic ovary . Hum Reprod Update , 22 ( 6 ): 709 - 724 . https://doi.org/10.1093/humupd/dmw027 https://doi.org/10.1093/humupd/dmw027 Interactions between androgens, FSH, anti-Müllerian hormone and estradiol during folliculogenesis in the human normal and polycystic ovary Faghihi MA , Modarresi F , Khalil AM , et al. , 2008 . Expression of a noncoding RNA is elevated in Alzheimer’s disease and drives rapid feed-forward regulation of β‍-secretase . Nat Med , 14 ( 7 ): 723 - 730 . https://doi.org/10.1038/nm1784 https://doi.org/10.1038/nm1784 Expression of a noncoding RNA is elevated in Alzheimer’s disease and drives rapid feed-forward regulation of β‍-secretase Fujii Y , Okada Y , Moore JP , et al. , 2002 . Evidence that PACAP and GnRH down-regulate follicle-stimulating hormone-‍β mRNA levels by stimulating follistatin gene expression: effects on folliculostellate cells, gonadotrophs and LβT2 gonadotroph cells . Mol Cell Endocrinol , 192 ( 1-2 ): 55 - 64 . https://doi.org/10.1016/s0303-7207(02)00109-0 https://doi.org/10.1016/s0303-7207(02)00109-0 Evidence that PACAP and GnRH down-regulate follicle-stimulating hormone-‍β mRNA levels by stimulating follistatin gene expression: effects on folliculostellate cells, gonadotrophs and LβT2 gonadotroph cells Gao YL , Zhao ZS , Zhang MY , et al. , 2017 . Long noncoding RNA PVT1 facilitates cervical cancer progression via negative regulating of miR-424 . Oncol Res Featur Preclin Clin Cancer Ther , 25 ( 8 ): 1391 - 1398 . https://doi.org/10.3727/096504017X14881559833562 https://doi.org/10.3727/096504017X14881559833562 Long noncoding RNA PVT1 facilitates cervical cancer progression via negative regulating of miR-424 Garrel G , Simon V , Thieulant ML , et al. , 2010 . Sustained gonadotropin-releasing hormone stimulation mobilizes the cAMP/PKA pathway to induce nitric oxide synthase type 1 expression in rat pituitary cells in vitro and in vivo at proestrus . Biol Reprod , 82 ( 6 ): 1170 - 1179 . https://doi.org/10.1095/biolreprod.109.082925 https://doi.org/10.1095/biolreprod.109.082925 Sustained gonadotropin-releasing hormone stimulation mobilizes the cAMP/PKA pathway to induce nitric oxide synthase type 1 expression in rat pituitary cells in vitro and in vivo at proestrus George JW , Dille EA , Heckert LL , 2011 . Current concepts of follicle-stimulating hormone receptor gene regulation . Biol Reprod , 84 ( 1 ): 7 - 17 . https://doi.org/10.1095/biolreprod.110.085043 https://doi.org/10.1095/biolreprod.110.085043 Current concepts of follicle-stimulating hormone receptor gene regulation Giraldez AJ , Mishima Y , Rihel J , et al. , 2006 . Zebrafish miR-430 promotes deadenylation and clearance of maternal mRNAs . Science , 312 ( 5770 ): 75 - 79 . https://doi.org/10.1126/science.1122689 https://doi.org/10.1126/science.1122689 Zebrafish miR-430 promotes deadenylation and clearance of maternal mRNAs Han DX , Sun XL , Fu Y , et al. , 2017a . Identification of long non-coding RNAs in the immature and mature rat anter ior pituitary. Sci Rep , 7 : 17780 . https://doi.org/10.1038/s41598-017-17996-6 https://doi.org/10.1038/s41598-017-17996-6 Identification of long non-coding RNAs in the immature and mature rat anter Han DX , Sun XL , Xu MQ , et al. , 2017b . Roles of differential expression of microRNA-21-3p and microRNA-433 in FSH regulation in rat anterior pituitary cells . Oncotarget , 8 ( 22 ): 36553 - 36565 . https://doi.org/10.18632/oncotarget.16615 https://doi.org/10.18632/oncotarget.16615 Roles of differential expression of microRNA-21-3p and microRNA-433 in FSH regulation in rat anterior pituitary cells Heckert L , Griswold MD , 1993 . Expression of the FSH receptor in the testis . In: Bardin CW (Ed.), Recent Progress in Hormone Research . Academic Press, Manhattan, p. 61 - 77 . https://doi.org/10.1016/b978-0-12-571148-7.50006-3 https://doi.org/10.1016/b978-0-12-571148-7.50006-3 Expression of the FSH receptor in the testis Hu GZ , Lou ZK , Gupta M , 2014 . The long non-coding RNA GAS5 cooperates with the eukaryotic translation initiation factor 4E to regulate c-Myc translation . PLoS ONE , 9 ( 9 ): e107016 . https://doi.org/10.1371/journal.pone.0107016 https://doi.org/10.1371/journal.pone.0107016 The long non-coding RNA GAS5 cooperates with the eukaryotic translation initiation factor 4E to regulate c-Myc translation Jiang W , Liu YT , Liu R , et al. , 2015 . The lncRNA DEANR1 facilitates human endoderm differentiation by activating FOXA2 expression . Cell Rep , 11 ( 1 ): 137 - 148 . https://doi.org/10.1016/j.celrep.2015.03.008 https://doi.org/10.1016/j.celrep.2015.03.008 The lncRNA DEANR1 facilitates human endoderm differentiation by activating FOXA2 expression Kanasaki H , Purwana IN , Mijiddorj T , et al. , 2012 . Effects of estradiol and progesterone on gonadotropin LHβ‍- and FSHβ‍-subunit promoter activities in gonadotroph LβT2 cells . Neuro Endocrinol Lett , 33 ( 6 ): 608 - 613 . Effects of estradiol and progesterone on gonadotropin LHβ‍- and FSHβ‍-subunit promoter activities in gonadotroph LβT2 cells Kaprara A , Huhtaniemi IT , 2018 . The hypothalamus-pituitary-gonad axis: tales of mice and men . Metabolism , 86 : 3 - 17 . https://doi.org/10.1016/j.metabol.2017.11.018 https://doi.org/10.1016/j.metabol.2017.11.018 The hypothalamus-pituitary-gonad axis: tales of mice and men Krol J , Loedige I , Filipowicz W , 2010 . The widespread regulation of microRNA biogenesis, function and decay . Nat Rev Genet , 11 ( 9 ): 597 - 610 . https://doi.org/10.1038/nrg2843 https://doi.org/10.1038/nrg2843 The widespread regulation of microRNA biogenesis, function and decay Kumar TR , 2018 . Extragonadal actions of FSH: a critical need for novel genetic models . Endocrinology , 159 ( 1 ): 2 - 8 . https://doi.org/10.1210/en.2017-03118 https://doi.org/10.1210/en.2017-03118 Extragonadal actions of FSH: a critical need for novel genetic models Lamba P , Fortin J , Tran S , et al. , 2009 . A novel role for the forkhead transcription factor FOXL2 in activin A-regulated follicle-stimulating hormone β subunit transcription . Mol Endocrinol , 23 ( 7 ): 1001 - 1013 . https://doi.org/10.1210/me.2008-0324 https://doi.org/10.1210/me.2008-0324 A novel role for the forkhead transcription factor FOXL2 in activin A-regulated follicle-stimulating hormone β subunit transcription Lannes J , L'hôte D , Fernandez-Vega A , et al. , 2016 . A regulatory loop between miR-132 and miR-125b involved in gonadotrope cells desensitization to GnRH . Sci Rep , 6 : 34676 . https://doi.org/10.1038/srep34676 https://doi.org/10.1038/srep34676 A regulatory loop between miR-132 and miR-125b involved in gonadotrope cells desensitization to GnRH Lei KC , Liang X , Gao YW , et al. , 2017 . Lnc-ATB contributes to gastric cancer growth through a miR-141-3p/TGFβ2 feedback loop . Biochem Biophys Res Commun , 484 ( 3 ): 514 - 521 . https://doi.org/10.1016/j.bbrc.2017.01.094 https://doi.org/10.1016/j.bbrc.2017.01.094 Lnc-ATB contributes to gastric cancer growth through a miR-141-3p/TGFβ2 feedback loop Li M , Duan LW , Li YX , et al. , 2019 . Long noncoding RNA/circular noncoding RNA-miRNA-mRNA axes in cardiovascular diseases . Life Sci , 233 : 116440 . https://doi.org/10.1016/j.lfs.2019.04.066 https://doi.org/10.1016/j.lfs.2019.04.066 Long noncoding RNA/circular noncoding RNA-miRNA-mRNA axes in cardiovascular diseases Liu XM , Huang SS , Guan Y , et al. , 2020 . Long noncoding RNA OSER1-AS1 promotes the malignant properties of non-small cell lung cancer by sponging microRNA-433-3p and thereby increasing Smad2 expression . Oncol Rep , 44 ( 2 ): 599 - 610 . https://doi.org/10.3892/or.2020.7645 https://doi.org/10.3892/or.2020.7645 Long noncoding RNA OSER1-AS1 promotes the malignant properties of non-small cell lung cancer by sponging microRNA-433-3p and thereby increasing Smad2 expression Long H , Wang X , Chen YJ , et al. , 2018 . Dysregulation of microRNAs in autoimmune diseases: pathogenesis, biomarkers and potential therapeutic targets . Cancer Lett , 428 : 90 - 103 . https://doi.org/10.1016/j.canlet.2018.04.016 https://doi.org/10.1016/j.canlet.2018.04.016 Dysregulation of microRNAs in autoimmune diseases: pathogenesis, biomarkers and potential therapeutic targets Ma LN , Bajic VB , Zhang Z , 2013 . On the classification of long non-coding RNAs . RNA Biol , 10 ( 6 ): 925 - 933 . https://doi.org/10.4161/rna.24604 https://doi.org/10.4161/rna.24604 On the classification of long non-coding RNAs Marshall JC , Kelch RP , 1986 . Gonadotropin-releasing hormone: role of pulsatile secretion in the regulation of reproduction . N Engl J Med , 315 ( 23 ): 1459 - 1468 . https://doi.org/10.1056/NEJM198612043152306 https://doi.org/10.1056/NEJM198612043152306 Gonadotropin-releasing hormone: role of pulsatile secretion in the regulation of reproduction Massague J , Seoane J , Wotton D , 2005 . Smad transcription factors . Genes Dev , 19 ( 23 ): 2783 - 2810 . https://doi.org/10.1101/gad.1350705 https://doi.org/10.1101/gad.1350705 Smad transcription factors Motta JCL , Madureira G , Silva LO , et al. , 2020 . Interactions of circulating estradiol and progesterone on changes in endometrial area and pituitary responsiveness to GnRH . Biol Reprod , 103 ( 3 ): 643 - 653 . https://doi.org/10.1093/biolre/ioaa065 https://doi.org/10.1093/biolre/ioaa065 Interactions of circulating estradiol and progesterone on changes in endometrial area and pituitary responsiveness to GnRH Noh JH , Kim KM , McClusky WG , et al. , 2018 . Cytoplasmic functions of long noncoding RNAs . WIREs RNA , 9 ( 3 ): e1471 . https://doi.org/10.1002/wrna.1471 https://doi.org/10.1002/wrna.1471 Cytoplasmic functions of long noncoding RNAs Peng WZ , Si S , Zhang QX , et al. , 2015 . Long non-coding RNA MEG3 functions as a competing endogenous RNA to regulate gastric cancer progression . J Exp Clin Cancer Res , 34 : 79 . https://doi.org/10.1186/s13046-015-0197-7 https://doi.org/10.1186/s13046-015-0197-7 Long non-coding RNA MEG3 functions as a competing endogenous RNA to regulate gastric cancer progression Pierce JG , Parsons TF , 1981 . Glycoprotein hormones: structure and function . Annu Rev Biochem , 50 : 465 - 495 . https://doi.org/10.1146/annurev.bi.50.070181.002341 https://doi.org/10.1146/annurev.bi.50.070181.002341 Glycoprotein hormones: structure and function Popovics P , Rekasi Z , Stewart AJ , et al. , 2011 . Regulation of pituitary inhibin/activin subunits and follistatin gene expression by GnRH in female rats . J Endocrinol , 210 ( 1 ): 71 - 79 . https://doi.org/10.1530/JOE-10-0485 https://doi.org/10.1530/JOE-10-0485 Regulation of pituitary inhibin/activin subunits and follistatin gene expression by GnRH in female rats Ransohoff JD , Wei YN , Khavari PA , 2018 . The functions and unique features of long intergenic non-coding RNA . Nat Rev Mol Cell Biol , 19 ( 3 ): 143 - 157 . https://doi.org/10.1038/nrm.2017.104 https://doi.org/10.1038/nrm.2017.104 The functions and unique features of long intergenic non-coding RNA Salmena L , Poliseno L , Tay Y , et al. , 2011 . A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Cell , 146 ( 3 ): 353 - 358 . https://doi.org/10.1016/j.cell.2011.07.014 https://doi.org/10.1016/j.cell.2011.07.014 A ceRNA hypothesis: the rosetta stone of a hidden RNA language? Santi D , Crepieux P , Reiter E , et al. , 2020 . Follicle-stimulating hormone (FSH) action on spermatogenesis: a focus on physiological and therapeutic roles . J Clin Med , 9 ( 4 ): 1014 . https://doi.org/10.3390/jcm9041014 https://doi.org/10.3390/jcm9041014 Follicle-stimulating hormone (FSH) action on spermatogenesis: a focus on physiological and therapeutic roles Schubert M , Pérez Lanuza L , Gromoll J , 2019 . Pharmacogenetics of FSH action in the male . Front Endocrinol , 10 : 47 . https://doi.org/10.3389/fendo.2019.00047 https://doi.org/10.3389/fendo.2019.00047 Pharmacogenetics of FSH action in the male Stamatiades GA , Carroll RS , Kaiser UB , 2019 . GnRH—a key regulator of FSH . Endocrinology , 160 ( 1 ): 57 - 67 . https://doi.org/10.1210/en.2018-00889 https://doi.org/10.1210/en.2018-00889 GnRH—a key regulator of FSH Stilley JAW , Segaloff DL , 2018 . FSH actions and pregnancy: looking beyond ovarian FSH receptors . Endocrinology , 159 ( 12 ): 4033 - 4042 . https://doi.org/10.1210/en.2018-00497 https://doi.org/10.1210/en.2018-00497 FSH actions and pregnancy: looking beyond ovarian FSH receptors Sun Y , Liu WZ , Liu T , et al. , 2015 . Signaling pathway of MAPK/ERK in cell proliferation, differentiation, migration, senescence and apoptosis . J Recept Sig Transd , 35 ( 6 ): 600 - 604 . https://doi.org/10.3109/10799893.2015.1030412 https://doi.org/10.3109/10799893.2015.1030412 Signaling pathway of MAPK/ERK in cell proliferation, differentiation, migration, senescence and apoptosis Tsutsumi R , Mistry D , Webster NJG , 2010 . Signaling responses to pulsatile gonadotropin-releasing hormone in LβT2 gonadotrope cells . J Biol Chem , 285 ( 26 ):‍ 20262 - 20272 . https://doi.org/10.1074/jbc.M110.132662 https://doi.org/10.1074/jbc.M110.132662 Signaling responses to pulsatile gonadotropin-releasing hormone in LβT2 gonadotrope cells Wang LM , Song FB , Zhu WB , et al. , 2021 . The stage-specific long non-coding RNAs and mRNAs identification and analysis during early development of common carp, Cyprinus carpio . Genomics , 113 ( 1 ): 20 - 28 . https://doi.org/10.1016/j.ygeno.2020.11.025 https://doi.org/10.1016/j.ygeno.2020.11.025 The stage-specific long non-coding RNAs and mRNAs identification and analysis during early development of common carp, Cyprinus carpio Wang LY , Cho KB , Li Y , et al. , 2019 . Long noncoding RNA (lncRNA)‍‍-mediated competing endogenous RNA networks provide novel potential biomarkers and therapeutic targets for colorectal cancer. Int J Mol Sci , 20 ( 2 2) :575 8. https://doi.org/10.3390/ijms202257 https://doi.org/10.3390/ijms202257 58 :575 Xu J , Wu KJ , Jia QJ , et al. , 2020 . Roles of miRNA and lncRNA in triple-negative breast cancer . J Zhejiang Univ-Sci B (Biomed & Biotechnol) , 21 ( 9 ): 673 - 689 . https://doi.org/10.1631/jzus.B1900709 https://doi.org/10.1631/jzus.B1900709 Roles of miRNA and lncRNA in triple-negative breast cancer Yu H , Li SB , 2020 . Role of LINC00152 in non-small cell lung cancer . J Zhejiang Univ-Sci B (Biomed & Biotechnol) , 21 ( 3 ): 179 - 191 . https://doi.org/10.1631/jzus.B1900312 https://doi.org/10.1631/jzus.B1900312 Role of LINC00152 in non-small cell lung cancer Zhao FL , Wu Y , Yang W , et al. , 2020 . Inhibition of vascular calcification by microRNA-155-5p is accompanied by the inactivation of TGF-β1/Smad2/3 signaling pathway . Acta Histochem , 122 : 151551 . https://doi.org/10.1016/j.acthis.2020.151551 https://doi.org/10.1016/j.acthis.2020.151551 Inhibition of vascular calcification by microRNA-155-5p is accompanied by the inactivation of TGF-β1/Smad2/3 signaling pathway Wenzhi REN, Bao YUAN, and Wei GAO were responsible for the main experimental concept and design. Weidi ZHANG, Guokun ZHAO, Haoqi WANG, and Haixiang GUO performed rat pituitary gland dissection. Weidi ZHANG, Dongxu HAN, and Guokun ZHAO performed other experiments. Weidi ZHANG, Haoqi WANG, Yi ZHENG, and Zhonghao JI performed the data analyses. Weidi ZHANG and Bao YUAN wrote the manuscript. All the authors have read and approved the final version, and therefore, have full access to all the data in the study and take responsibility for the integrity and security of the data. This work was supported by the National Natural Science Foundation of China (No. 31872349). This work was supported by the National Natural Science Foundation of China (No. 31872349). 10.1631/jzus.B2101052 1673-1581（2022）06-0502-13 This work was supported by the National Natural Sci‐ence Foundation of China (No. 31872349). 2021-12-13 2022-02-04 2022-06-15 2022-12-27 2022-05-13 Copyright © Zhejiang University and Springer-Verlag GmbH Germany, part of Springer Nature 2022 2022 浙江大学出版社 Weidi ZHANG, Wenzhi REN, Dongxu HAN, Guokun ZHAO, Haoqi WANG, Haixiang GUO, Yi ZHENG, Zhonghao JI, Wei GAO, and Bao YUAN declared that they have no conflict of interest. All animals were fed at the Jilin Provincial Key Laboratory of Animal Models, and the animal procedures were carried out in accordance with the relevant guidelines of the Jilin University Guide for the Care and Use of Laboratory Animals. All animal procedures were performed in accordance with the protocol approved by the Institutional Animal Care and Use Committee (IACUC) of Jilin University (permit number: SY201912021), China. Compliance with ethics guidelines Supplementary information 1 Fig. S1; Table S1 浙江大学学报（英文版）（B辑：生物医学和生物技术） 6 23