- 1International Peace Maternity and Child Health Hospital, School of Medicine, Shanghai Jiao Tong University, Shanghai, China
- 2Shanghai Key Laboratory of Embryo Original Disease, Shanghai, China
- 3Obstetrics and Gynecology Hospital, Institute of Reproduction and Development, Fudan University, Shanghai, China
- 4Research Units of Embryo Original Diseases, Chinese Academy of Medical Sciences, Shanghai, China
Polycystic ovary syndrome (PCOS), a common and frustrating syndrome in women of reproductive age, is characterized by symptoms including hyperandrogenemia, ovulation dysfunction, and polycystic ovaries. The role of competitive endogenous RNA (ceRNA) networks is receiving increasing attention and has been reported in multiple complicated diseases, such as various carcinomas, endometriosis, and tubal factor infertility. However, the association of ceRNA networks with the pathogenesis of PCOS remains unclear. This study aimed to construct a ceRNA network orchestrated by exosomal lnRNA and circRNA in PCOS. We screened RNA data of 34 samples from the Gene Expression Omnibus (GEO) database for differentially expressed lncRNAs (DELs), miRNAs (DEMs), mRNAs (DEGs), and circRNA associated with the progression of PCOS (PCOS, n = 17 vs. normal, n = 17). A protein–protein interaction (PPI) network, gene set enrichment analysis (GSEA), and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were conducted. Importantly, the function of the ceRNA network was explored using GO and KEGG enrichment analyses. We identified 46 DELs (25 upregulated and 21 downregulated), 31 DEMs (20 upregulated and 11 downregulated), 165 DEGs (52 upregulated and 113 downregulated), and 1 differentially expressed circRNA. The PPI network had 79 nodes and 112 edges. The GSEA results showed that these genes were mainly related to oxidative phosphorylation; TNF signaling pathways; and valine, leucine, and isoleucine degradation. GO and KEGG analyses revealed that the DEGs were significantly enriched in lipid metabolism, peroxisome proliferator-activated receptor (PPAR) signaling pathways, and fatty acid metabolism. Additionally, we constructed a novel PCOS-associated lncRNA–miRNA–mRNA ceRNA triple network and a circRNA-related network. Thereafter, we described the potential roles played by follicular fluid exosomes in PCOS. Our present study describes the molecular pathogenesis of PCOS in human ovarian granulosa cells at the post-transcriptional level, which provides new insights for the clinical diagnosis and treatment of PCOS and further scientific research.
1 Introduction
Polycystic ovary syndrome (PCOS) is a common disease in women and is associated with endocrine and metabolic disorders. Approximately 10.1% Chinese women of reproductive age are affected by PCOS (Wu et al., 2021). The disease has two typical characteristics: androgen excess (hyperandrogenism) and ovarian dysfunction [oligoovulation and/or polycystic ovary pattern (PCOM)]. In addition, most patients exhibit significant metabolic abnormalities, mainly insulin resistance and compensatory hyperinsulinemia (Ezeh et al., 2021). The most widely used classification of PCOS is the Rotterdam definition, which is currently supported by most scientific and health sectors (Legro et al., 2013). The definition suggests that any woman with at least two of the following three diseases can be diagnosed with PCOS: clinical or biochemical hyperandrogenism, ovulation disorders, and PCOM (Rotterdam ESHRE/ASRM-Sponsored PCOS Consensus Workshop Group, 2004). PCOS is caused by various factors, including androgen imbalance, increased luteinizing hormone (LH) levels, inflammatory response, and oxidative stress, and its pathophysiological mechanism is complex (Li et al., 2017).
Granulosa cells, a layer of cells occurring on the surface of follicles, are an important component of the ovary. Mounting evidence indicates that the disruption of follicle development and ovulation in PCOS is associated with granulosa cell dysfunction (Liu et al., 2022; Zhang et al., 2022). In addition, mural granulosa cells participate in oocyte meiosis arrest by secreting oocyte maturation inhibitory factor (OMI) (Kawamura et al., 2011). Studies have also found that granulosa cell proliferation is increased in PCOS patients (Bao et al., 2021). Additionally, the abnormal proliferation of granulosa cells may increase the number of preantral follicles and ovarian volume, which may lead to the formation of large follicles in the ovaries of PCOS patients. Therefore, the abnormal state of granulosa cells is associated with the occurrence of PCOS, suggesting that improving the abnormal state of granulosa cells has potential therapeutic value.
Exosomes are membrane-wrapped vesicles, 30–150 nm in diameter, and are secreted by multiple cells into the extracellular microenvironment. They contain various types of nucleic acids, including mRNA, microRNA (miRNA), and long non-coding RNA (lncRNA), and are indispensable for extracellular communication (Kalluri and LeBleu, 2020; Donoso-Quezada et al., 2021). Recent literature shows that exosomes can be synthesized and released in different parts of the female reproductive tract, such as in the tubal epithelial cells, endometrium, uterus, and follicular fluid, which is the key microenvironment for follicular development (Machtinger et al., 2016; Gurung et al., 2020; Li et al., 2020a; Giacomini et al., 2021). Exosomes can transmit information to granulosa cells (Pegtel and Gould, 2019). Abnormalities in either exosomes or granulosa cells play an important role in the development of PCOS. However, the mechanism of communication between exosomes and granulosa cells, as well as their therapeutic targets, have not been well understood. Therefore, an in-depth exploration of exosomes and granulosa cells and their role in the occurrence and development of PCOS is required.
The vast majority of mammalian genomes are transcribed as non-coding RNA, for instance, miRNA, lncRNA, and circular RNA (circRNA) (Yao et al., 2019). Recently, lncRNA has captured the attention of many scientists. A study on the lncRNA of follicular fluid exosomes using high-throughput sequencing screened 1,866 differentially expressed lncRNAs. Functional analysis showed that they were closely related to pathways such as endocytosis, Hippo, and MAPK. This suggests that lncRNAs are important in the pathogenesis of PCOS (Wang et al., 2020). Studies have identified differential expression profiles of lncRNAs and mRNAs in exosomes and have revealed some biological processes and regulatory networks related to PCOS (Zhou et al., 2022). CircRNAs are covalently closed structures (Chen, 2016). At present, a few studies have found that circRNAs are involved in the occurrence and development of PCOS. For example, Jia et al. reported that circ_0118530 is increased in PCOS patients (Jia et al., 2020).
miRNAs consist of 9–25 nucleotides (Wang et al., 2015) and are also crucial for the development of PCOS (Abdalla et al., 2020; Bao et al., 2021). For example, one study found that miR-200b expression was significantly increased in PCOS patients compared to that in controls. Moreover, the proliferation of a granulosa-like tumor cell line was inhibited by the overexpression of miR-200b and miR-200c (He et al., 2019). The expression of miR-93 and miR-21 was increased in the granulosa cells of hyperandrogenic PCOS patients, and the expression of miR-93 and miR-21 was positively correlated with free testosterone and free androgen index (Tu et al., 2019). However, different sequencing platforms, different statistical methods, and small sample sizes can all lead to erratic study results. Further comprehensive analyses are necessary to reliably identify the molecular mechanisms of disease development.
The theory about competitive endogenous RNA (ceRNA) was first put forward by Salmena et al. (2011). Subsequently, experiments have confirmed that some RNAs act as ceRNA, such as cIRS-7 (also known as CDR1-AS, a circRNA) regulating the activity of miR-7 in the central nervous system (Memczak et al., 2013). The ability of pseudogenes, such as PTENP1, to act as ceRNAs has also been experimentally demonstrated (Yu et al., 2014). Furthermore, lncRNAs can interact with miRNAs, messenger RNAs (mRNAs), and proteins for performing their biological functions. Additionally, lncRNAs act as ceRNAs that communicate with mRNAs through competitively shared miRNAs (Ergun and Oztuzcu, 2015; Piao et al., 2020); the specific mechanism is as follows. When ceRNA is transcriptionally silenced, the parental mRNA is transcribed and exported to the cytoplasm, where it is targeted by the microRNA-guided RNA-induced silencing complex (miRNA-RISC), resulting in accelerated degradation, translation blockade, and reduced expression. When the ceRNA with competing target sites becomes transcriptional active, it will compete for miRNA targeting and binding to the RISC complex. This isolates the miRNA-RISC complex from the parent gene and results in the increased expression of the parent gene (Thomson and Dinger, 2016).
The role of ceRNA networks is receiving increasing attention and has been reported in several complicated diseases, including carcinomas (Chen et al., 2022; Xiong et al., 2022), endometriosis (Li et al., 2022), and tubal factor infertility (Li J. Z. et al., 2021a). Some previous studies have constructed ceRNA networks for the development of PCOS. For example, a study constructed an lncRNA–mRNA network associated with PCOS on the basis of granulosa cells, and key lncRNAs were found to be significantly associated with processes and pathways related to inflammation, oxidative stress, and apoptosis (Ma et al., 2021). Other studies have also constructed ceRNA networks and identified ceRNA axes such as lncRNAs miR-483-5p-GOT2; most of the axes are closely associated with steroid biosynthesis and metabolic pathways (Zhao et al., 2021). However, few studies have investigated PCOS-associated lncRNA–miRNA–mRNA ceRNA triple networks between follicular fluid and granulosa cells. In the current study, to explore the possible mechanisms of PCOS pathogenesis in human ovarian granulosa cells at the post-transcriptional level, we conducted investigations (described in the flowchart in Figure 1 and Supplementary Figure S1), identified DELs in exosomes and DEMs and DEGs in granulosa cells, and constructed an integrated analysis of the lncRNA–miRNA–mRNA ceRNA network.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/979924/fgene-13-979924-HTML/image_m/fgene-13-979924-g001.jpg)
FIGURE 1. Workflow of data processing and analysis. DEL, Differentially expressed lncRNA; DEM, differentially expressed miRNA; DEG, differentially expressed mRNA; lncRNA, long non-coding RNA; miRNA, microRNA; mRNA, messenger RNA; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; GSEA, gene set enrichment analysis; ceRNA, competing endogenous RNA.
2 Materials and methods
2.1 Data collection
We searched the Gene Expression Omnibus (GEO) database for datasets related to PCOS and exosomes using the keywords “polycystic ovary syndrome” and “exosomes” and downloaded four datasets, namely, GSE145461, GSE138572, GSE193123, and GSE155489. To obtain data on lncRNA and circRNA expressions in the follicular fluid exosomes of patients with PCOS, the GSE145461 dataset (platform: GPL20301 Illumina HiSeq 4,000) was used. This dataset had been obtained by deep sequencing of the transcriptome of follicular fluid from five PCOS patients and five control patients using Illumina GAIIx. To obtain data on miRNA expression in granulosa cells of patients with PCOS, the GSE138572 dataset was used (platform: GPL11154 Illumina HiSeq, 2000). This dataset had been developed using ovarian granulosa cells from five PCOS patients and five healthy volunteers. The dataset GSE193123 (platform: GPL24676 Illumina NovaSeq 6,000) contains gene/mRNA expression data obtained from ovarian granulosa cells from three healthy women and three women with PCOS. The dataset GSE155489 (platform: GPL20795 HiSeq X Ten) contains gene/mRNA expression data obtained from granulosa cells from four healthy women and four women with PCOS.
2.2 Identification of differentially expressed genes
Differentially expressed lncRNAs (DELs) in exosome samples and differentially expressed miRNAs (DEMs) and differentially expressed mRNAs (DEGs) in granulosa cells between patients with PCOS and healthy individuals were screened using DEseq2 R packages (Love et al., 2014). Differentially expressed circRNAs in exosome samples between patients with PCOS and healthy individuals were screened using limma and edgeR packages (Ritchie et al., 2015), respectively. For all screening analyses, p-value < 0.05 was deemed statistically significant.
2.3 Construction of a protein–protein interaction network
We used STRING 11.0 (https://string-db.org/), an online database, to study important genes associated with PCOS in granulosa cells. A PPI network on DEGs was constructed with a confidence score ≥ 0.4. We generated the results using Cytoscape software (V2.8.3) (Shannon et al., 2003).
2.4 Functional enrichment analysis
To explore the possible biological role of granulosa cells in PCOS, GSEA was performed by the clusterProfiler R package. Moreover, to forecast the possible functions of these DEGs, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (Kanehisa et al., 2012) and Gene Ontology (GO) analyses of DEGs were conducted using the clusterProfiler R package. GO analysis was conducted from three distinct aspects: biological process (BP), cellular component (CC), and molecular function (MF). p-value < 0.05 was considered statistically significant. To predict the function of these DEGs with high accuracy, we also performed GO analysis and KEGG analysis on DEGs using ClueGO v 2.5.8. ClueGO is a plugin for Cytoscape that allows for the classification of non-redundant GO terms into categories and the visualization of functionally related genes in a clustering network (Bindea et al., 2009). The threshold was set at p-value < 0.05.
2.5 Construction of lncRNA–miRNA–mRNA ceRNA network
lncRNA–miRNA pairs and miRNA–mRNA pairs can form lncRNA–miRNA–mRNA triplets. miRNA can bind to targeted mRNA to promote the degradation of mRNA, whereas lncRNA can bind to targeted miRNA to inhibit the degradation of mRNA. Here, we used the ggalluvial R package (Rosvall and Bergstrom, 2010) to construct lncRNA–miRNA–mRNA triplets on the basis of DEGs and DELs through miRcode v 11(Jeggari et al., 2012; http://www.mircode.org/mircode/), miRDB v 7.0 (Wong et al., 2015; http://mirdb.org/), miRTarBase (Huang et al., 2020; https://mirtarbase.cuhk.edu.cn/∼miRTarBase/miRTarBase_2022/php/index.php), and TargetScan v 7.2 (Enright et al., 2003; http://targetscan.org/vert_72/).
First, we predicted the lncRNA–miRNA pairs through the miRcode database on the base of DELs. miRcode is a human microRNA target prediction database based on GENCODE gene annotation. Then, we used three databases, miRDB, miRTarBase, and TargetScan, to predict the target genes of miRNAs. miRDB is an online database for target prediction and functional annotation of miRNAs through the MirTarget bioinformatics tool. miRTarBase is a microRNA target database with experimentally confirmed data, and TargetScan is a software for predicting miRNA binding sites in mammals. We set the parameters such that genes appearing in any two of these databases were considered target genes for these miRNAs. Subsequently, we compared the predicted target genes with DEGs and used overlapping genes and the pairs of roles corresponding to them to build lncRNA–miRNA–mRNA triplets (ceRNA network).
2.6 GO and KEGG enrichment analysis of the ceRNA network
To further understand the function of genes in the ceRNA network, we accomplished GO and KEGG path analyses with the clusterProfiler R package on our ceRNA network to predict the function of this network. p-value < 0.05 was considered statistically significant.
2.7 CircRNA-related network
First, considering the differentially expressed circRNAs, we predicted the circRNA–miRNA pairs using the circinteractome (Dudekula et al., 2016) and the cancer-specific circRNA database (CSCD) (Xia et al., 2018). We then set the parameters such that genes appearing in two of these databases were considered target miRNAs for these circRNAs. Subsequently, we used TargetScan, miRDB, and miRTarBase databases to predict mRNAs as previously described to build the circRNA–miRNA–mRNA network.
3 Results
3.1 Identification of differentially expressed genes
We identified 46 DELs (25 upregulated and 21 downregulated) in GSE145461 (Figure 2A). Additionally, 31 DEMs (20 upregulated and 11 downregulated) were identified in GSE138572. Further, 1,328 differentially expressed mRNAs (601 upregulated and 727 downregulated) were identified in GSE193123, and 5,268 differentially expressed mRNAs (2,530 upregulated and 2,738 downregulated) were identified in GSE155489. Next, 165 DEGs (52 upregulated and 113 downregulated) were identified (Figure 2B). However, only one differentially expressed circRNA was identified (Supplementary Tables S1, S2, Supplementary Figure S2). The identified DELs, DEMs, and DEGs were used in the next analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/979924/fgene-13-979924-HTML/image_m/fgene-13-979924-g002.jpg)
FIGURE 2. Differential analysis of lncRNA, miRNA, mRNA, mRNA crossover genes, and PPI network analysis. (A) Difference analysis. Volplot of DELs, DEMs, and DEGs. (B) Overlapped DEGs screening. (C) PPI network of DEGs. Red nodes represent upregulated genes in PCOS, and green nodes represent downregulated genes in PCOS.
3.2 PPI network construction
We utilized the STRING database to build the PPI network of DEGs. One hundred and sixty-five DEGs were delivered to the STRING database, and 79 genes were mapped to the PPI network. The PPI network of DEGs contained 79 nodes and 112 edges, and their PPI enrichment p-values were below 0.05 (Figure 2C). In this network, several important genes are downregulated (e.g., SCD, ACADM, and APOE), whereas only a few genes are upregulated (e.g., ACSL6, RELA, and PISD).
3.3 Functional enrichment analysis
As seen in Figure 3A, GSEA showed that these genes were mainly associated with the oxidative phosphorylation pathway; TNF signaling pathway; and valine, leucine, and isoleucine degradation. As illustrated in Figure 3B, Figure 3C, and Supplementary Tables S1, the results of the GO analysis demonstrated that the BP of DEGs was significantly enriched in carboxylic acid decomposition, organic acid catabolism, small-molecule catabolism, monocarboxylic acid biosynthesis, and fatty acid metabolism. Changes in the CC of DEGs were mainly observed in the mitochondrial matrix, peroxisomal matrix, microbody lumen, peroxisome, and microbody. Changes in the MF were mainly observed in monocarboxylic acid, carboxylic acid, amide, and peptide hormone binding and in oxidoreductase activity, acting on the CH-CH group of donors. GO analysis from ClueGO showed that DEGs were mainly enriched in the monocarboxylic acid biosynthetic process, CoA carboxylase activity, monocarboxylic acid binding, and carboxylic acid catabolic process. Similarly, KEGG pathway analysis revealed that the DEGs were primarily enriched in valine, leucine, and isoleucine degradation; the PPAR signaling pathway; biosynthesis of unsaturated fatty acids; ovarian steroidogenesis; and fatty acid metabolism. KEGG analysis from ClueGO indicated that the DEGs were mainly enriched in valine, leucine, and isoleucine degradation; the PPAR signaling pathway; ovarian steroidogenesis; adherens junction; GnRH signaling pathway; and glycine, serine, and threonine metabolism (Figure 4 and Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/979924/fgene-13-979924-HTML/image_m/fgene-13-979924-g003.jpg)
FIGURE 3. GSEA enrichment analysis and GO analysis from ClueGO. (A) GSEA analysis. (B) GO analysis of DEGs from clusterprofiler. (C) GO analysis of DEGs from ClueGO (version: v.2.5.8).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/979924/fgene-13-979924-HTML/image_m/fgene-13-979924-g004.jpg)
FIGURE 4. KEGG enrichment analysis and KEGG analysis from ClueGO. (A) KEGG analysis of DEGs from clusterprofiler. (B) KEGG analysis of DEGs from ClueGO (version: v.2.5.8).
3.4 ceRNA network of lncRNA–miRNA–mRNA construction
As illustrated in Figure 5A, a ceRNA network was constructed including seven lncRNAs (e.g., SAMD12 and GNAQ), two miRNAs (hsa-miR-508-3p and hsa-miR-20b-5p), and nine mRNAs (e.g., CNKSR3, ITGB8, and HPS5). In the ceRNA network that we mapped, we found seven lncRNAs, of which six were upregulated and one was downregulated. The two identified miRNAs identified were upregulated. Of the nine mRNAs identified, two were upregulated and seven were downregulated. These results suggest that the lncRNA carried by exosomes may regulate the expression of some mRNAs in PCOS by hsa-miR-508-3p and hsa-miR-20b-5p.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/979924/fgene-13-979924-HTML/image_m/fgene-13-979924-g005.jpg)
FIGURE 5. Construction of ceRNA network and its functional enrichment analysis. (A) ceRNA network. Red nodes represent upregulated lncRNA/miRNA/mRNA in PCOS, and green nodes represent downregulated lncRNA/miRNA/mRNA genes in PCOS. (B) GO analysis of ceRNA network. (C) KEGG analysis of ceRNA network.
3.5 Functional enrichment analysis of ceRNA network
We performed enrichment analysis of GO function and the KEGG pathway. The top 10 enriched GO items are listed in Figure 5B and Table 2. GO analysis revealed that the ceRNA network was significantly enriched in transforming growth factor beta activation, urate metabolic process, response to the redox state, myeloid dendritic cell differentiation, regulation of Rac protein signal transduction, secondary palate development, ganglioside metabolic process, positive regulation of behavior, myeloid dendritic cell activation, and negative regulation of peptidyl-serine phosphorylation. KEGG pathway enrichment analysis was further conducted for this ceRNA network. As listed in Figure 5C and Table 3 (all six signal pathways), this ceRNA network suggested significant enrichment in the circadian rhythm, adherens junction, ECM–receptor interaction, etc.
3.6 Construction of circRNA-related network
We used circinteractome and CSCD database to predict miRNAs of differentially expressed circRNA, and only two miRNAs were found after intersection (Supplementary Tables S3, S4, Supplementary Figure S3). However, both these miRNAs had no intersection with DEM (Supplementary Figure S4). We then used miRDB, miRTarBase, and TargetScan database to predict the mRNA of the two miRNAs identified (Supplementary Table S5). The target mRNA and DEG had an intersection gene, HMGA1 (Supplementary Figure S5). Using the results obtained, we constructed the circRNA-related network (Supplementary Figure S6).
4 Discussion
Exosomes are a newly discovered component of the follicular fluid in ovaries. They deliver miRNAs, mRNAs, and proteins related to ovarian follicle development and oocyte maturation (Uzbekova et al., 2020; Cao et al., 2022). Exosomes in the follicular fluid can transmit genetic information to target cells, thereby changing the function of target cells, which is a novel mechanism for explaining the etiology of PCOS (Pegtel and Gould, 2019). As one of the “goods” delivered by exosomes, lncRNAs potentially play various important regulatory roles in the pathogenesis of PCOS (Tu et al., 2021). Importantly, in the ceRNA hypothesis, lncRNAs are called miRNA sponges because they can bind miRNAs and prevent miRNA function, and lncRNAs take advantage of this by playing a regulatory role (Tan et al., 2015; Duan et al., 2020). Therefore, we need to further explore the relationship between the lncRNAs of exosomes and the miRNA and mRNA in granulosa cells to understand the disease progression of PCOS. Some studies have targeted exosomes or granulosa cells to construct ceRNA networks. The novelty of our study is that it is based on the construction of ceRNA networks employing lncRNA in exosomes and miRNA and mRNA in granulosa cells. We found that the pathway was mainly enriched in the TNF signaling pathway, PPAR signaling pathway, fatty acid metabolism, circadian rhythm, etc.
The PPI networks proposed that some key genes, such as SCD, APOE, and ACADM, were downregulated in patients with PCOS. From the KEGG analysis graph plotted by ClueGO, we can see that SCD is in the PPAR pathway, and ACADM is in the PPAR pathway and valine, leucine, and isoleucine degradation pathway. SCD in humans has two subtypes and is closely associated with lipid metabolism and insulin resistance. One study used non-targeted metabolomics to study the serum metabolic profile and found increased stearoyl coenzyme A desaturase (SCD) activity in PCOS patients (Zhao et al., 2014). However, we analyzed genes from granulosa cells showing reduced SCD expression. Gene expression in different tissues cannot be compared directly. The APOE ε2 allele appears to be associated with abdominal obesity, insulin resistance, and metabolic syndrome in a group of women from southwest China who have PCOS (Liu et al., 2013). The ACADM gene encodes a medium-chain acyl coenzyme A dehydrogenase, which is involved in mitochondrial medium-chain fatty acid β-oxidation. At present, little research has focused on SCD, APOE, and ACADM in PCOS. Therefore, these genes are relatively recently discovered molecules in PCOS that should be researched further.
Through enrichment analysis, we conclude that PCOS is mainly enriched in the following pathways: oxidative phosphorylation; glycogen biosynthesis; TNF, PPAR, and MAPK signaling; fatty acid metabolic process; lipid catabolic process, etc. One report consistent with our study showed that DEGs in granulosa cells of PCOS are mainly involved in oxidative phosphorylation (Cozzolino et al., 2022). In addition, our observations are similar to the results of a previous study that found that DEGs in granulosa cells are mainly enriched in steroid biosynthesis and metabolism-related signaling, such as glycolysis or gluconeogenesis (Zhao et al., 2021). Moreover, studies from various countries have shown that patients with PCOS tend to have abnormal glucose metabolism, and some of them meet the WHO diagnostic criteria for type 2 diabetes mellitus (T2DM) (Ganie et al., 2016).
PCOS is also considered to be involved in chronic low-grade inflammation, and long-term inflammation may lead to hyperandrogenism unrelated to obesity (Rudnicka et al., 2021; Zhang et al., 2021). Women with PCOS have elevated levels of inflammatory cytokines such as IL-18, IL-1, TNF-α, and IL-6 (Liu et al., 2021; Rudnicka et al., 2021). Moreover, it was recently documented that TNF-α, a member of the adipokine family, may be associated with hyperandrogenism, increased insulin resistance, and obesity, which are common factors in the pathogenesis of PCOS (Escobar-Morreale et al., 2001; Alzamil, 2020; Uzdogan et al., 2021). In addition, in patients with PCOS, several dietary factors can induce oxidative stress and stimulate an inflammatory response, which may be related to insulin resistance in the disease. This oxidative stress may promote hyperandrogenism by reducing levels of the sex hormone-binding globulin (Sun et al., 2021). Besides, peroxisome proliferator-activated receptor-γ (PPAR-γ), a ligand-activated nuclear transcription factor, has been reported to be engaged in the inflammatory response in PCOS (San-Millán and Escobar-Morreale, 2010). One study suggested that the link between obesity and insulin resistance appears to be inflammatory, as pro-inflammatory cytokines are mainly produced by adipose tissue, and alterations in inflammation serum markers appear to be a feature of PCOS itself, unrelated to the PPAR variant (Knebel et al., 2008). Studies by other researchers support these findings, demonstrating that differentially expressed genes are engaged in the MAPK signaling pathway and other signaling pathways related to inflammation (Zhou et al., 2022). MAPK is the key to ovarian follicle development and ovulation. It is well known that activity of MAPK pathway is crucial for meiotic resumption and oocyte expansion in mouse oocytes (Su et al., 2002). Other studies have also found that MAPK is involved in regulating ovulation dysfunction in PCOS patients (Zhou et al., 2020). Another report shows that sterol regulatory element binding protein (SREBP-1), as the promoter of toll-like receptor (TLR)-4 signal transduction, can reprogram fatty acid metabolism, thereby inducing a pro-inflammatory response in PCOS (Wang et al., 2018). Our present results are consistent with these studies.
We constructed a ceRNA network based on the above results. Liu et al. (2019) observed that lncRNA SAMD12-AS1 was highly expressed in liver cancer tissues, and SAMD12-AS1 promoted cell growth and prevented apoptosis. However, if SAMD12-AS1 was knocked down, cell proliferation was inhibited, and apoptosis was enhanced. High lncRNA SAMD12-AS1 expression in glioma cells may increase the invasive and migratory capacity of cells via p53, and its high expression in patients may put them at increased risk of developing lymph node or distant metastases (Yu et al., 2019). In addition, lnc-GNAQ-6:1 expression is low in serum exosomes of gastric cancer patients and could be used as a novel diagnostic marker for gastric cancer by increasing sample size (Li et al., 2020b). Although lncRNAs have been studied in different cancerous cells and tissues, similar studies pertaining to PCOS are lacking. Therefore, in-depth studies should be conducted to address the research gap in the future. In addition, it has been shown that miRNAs mediate competition between ceRNAs and thus constitute an additional level of post-transcriptional regulation, which is essential in plenty of biological settings (Figliuzzi et al., 2013). For hsa-miR-508-3p, several studies have discovered that HSPA8 is highly expressed in acute myeloid leukemia, and HSPA8 is associated with hsa-miR-508-3p (Li and Ge, 2021). Another study showed that MiR-508-3p inhibits apoptosis and promotes the proliferation of middle ear cholesteatoma cells via the PTEN/PI3K/AKT pathway (Liu and Ma, 2021). For hsa-miR-20b-5p, studies have found that the MALAT1/hsa-mir-20b-5p/TXNIP axis plays an important role in liver inflammation in patients with chronic hepatitis B virus infection and non-alcoholic fatty liver disease (Li J. et al., 2021b). In addition, many studies have investigated the relationship between hsa-miR-20b-5p and cancer (Drobna et al., 2020; Filetti et al., 2021). However, we have not seen hsa-miR-508-3p and hsa-miR-20b-5p reported in studies on PCOS. Therefore, we may be able to further study these two miRNAs in vivo and in vitro in the future.
Through functional pathway analysis, we found that the ceRNA network potentially regulates multiple signaling pathways, such as the circadian rhythm. There is indeed some evidence to indicate that the development of PCOS is associated with a person’s circadian rhythm. Levels of melatonin in serum (Terzieva et al., 2013; Alizadeh et al., 2021), follicular fluid (Xie et al., 2021), and urine (Shreeve et al., 2013) and its metabolites are elevated in PCOS patients. Physiologically, the female reproductive system is governed by the hypothalamic–pituitary–gonadal axis in a time-dependent manner; thus, an abnormal daily rhythm would contribute to PCOS (Shao et al., 2021). In addition, one study showed night-shift work to be significantly associated with PCOS (Wang et al., 2021).
However, due to the limitation of data sets, we only found one differential circRNA, and the miRNA predicted by circRNA was not one of the DEMs from GSE138572; this needs to be further studied by increasing the datasets. Nevertheless, we identified a gene, HMAG1. HMGA1 has a wide range of roles in insulin resistance, glucose homeostasis, lipid metabolism, and atherosclerosis, and has recently been proposed as a convincing molecule involved in the pathological features of overlapping two diseases, PCOS and cardiovascular disease (De Rosa et al., 2018).
Although bioinformatics analysis is a powerful tool for understanding the pathogenesis of PCOS, further validation of these results at the molecular, cellular, and tissue levels is necessary. In the future, specific mechanisms involved in the regulation of disease development by these lncRNAs, miRNA, and mRNA, and intercellular communication would require further study. Other researchers may develop a more comprehensive ceRNA network based on a larger sample size, better database, and better algorithm. Importantly, more in-depth studies of these RNAs in vivo and in vitro are required.
Taken together, this study demonstrates the expression profiles of lncRNAs in exosomes, miRNAs, and mRNAs in granulosa cells, and presents ceRNA networks of exosomes and granulosa cells. The aberrant expression of these lncRNAs, miRNAs, and genes may be associated with multiple pathogeneses in patients with PCOS. This novel ceRNA network may prove useful for the study of the pathogenesis and treatment of PCOS.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
LB and XL conceived the idea and designed the research. LB and YL analyzed the data. LB, YG, and JG drafted the first version of the manuscript. XL and HH reviewed the manuscript. XL was responsible for the funding application and the supervision and management. All authors reviewed the intellectual content of the manuscript and approved the submission of the final version of the manuscript.
Funding
This work was sponsored by the National Natural Science Foundation of China (82171687), National Natural Science Funds (82192873), the CAMS Innovation Fund for Medical Sciences (2019-I2M-5-064), the National Natural Science Foundation of China (82088102), the Shanghai Frontiers Science Center of Reproduction and Development, and the Collaborative Innovation Program of Shanghai Municipal Health Commission (2020CXJQ01).
Acknowledgments
We would like to thank the editors (www.editage.com) for English language editing.
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.
The reviewers XL and JY declared a shared affiliation with the authors LB, YG, and TL to the handling editor at time of review
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.979924/full#supplementary-material
References
Abdalla, M., Deshmukh, H., Atkin, S. L., and Sathyapalan, T. (2020). miRNAs as a novel clinical biomarker and therapeutic targets in polycystic ovary syndrome (PCOS): a review. Life Sci. 259, 118174. doi:10.1016/j.lfs.2020.118174
Alizadeh, M., Karandish, M., Asghari Jafarabadi, M., Heidari, L., Nikbakht, R., Babaahmadi Rezaei, H., et al. (2021). Metabolic and hormonal effects of melatonin and/or magnesium supplementation in women with polycystic ovary syndrome: A randomized, double-blind, placebo-controlled trial. Nutr. Metab. 18, 57. doi:10.1186/s12986-021-00586-9
Alzamil, H. (2020). Elevated serum TNF-α is related to obesity in type 2 diabetes mellitus and is associated with glycemic control and insulin resistance. J. Obes. 2020, 5076858. doi:10.1155/2020/5076858
Bao, D., Li, M., Zhou, D., Zhuang, C., Ge, Z., Wei, Q., et al. (2021). miR-130b-3p is high-expressed in polycystic ovarian syndrome and promotes granulosa cell proliferation by targeting SMAD4. J. Steroid Biochem. Mol. Biol. 209, 105844. doi:10.1016/j.jsbmb.2021.105844
Bindea, G., Mlecnik, B., Hackl, H., Charoentong, P., Tosolini, M., Kirilovsky, A., et al. (2009). ClueGO: A cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25, 1091–1093. doi:10.1093/bioinformatics/btp101
Cao, J., Huo, P., Cui, K., Wei, H., Cao, J., Wang, J., et al. (2022). Correction: Follicular fluid-derived exosomal miR-143-3p/miR-155-5p regulate follicular dysplasia by modulating glycolysis in granulosa cells in polycystic ovary syndrome. Cell. Commun. Signal. 20, 116. doi:10.1186/s12964-022-00938-9
Chen, L. L. (2016). The biogenesis and emerging roles of circular RNAs. Nat. Rev. Mol. Cell. Biol. 17, 205–211. doi:10.1038/nrm.2015.32
Chen, X., Guo, J., Zhou, F., Ren, W., Huang, X., Pu, J., et al. (2022). Long non-coding RNA AL139385.1 as a novel prognostic biomarker in lung adenocarcinoma. Front. Oncol. 12, 905871. doi:10.3389/fonc.2022.905871
Cozzolino, M., Herraiz, S., Titus, S., Roberts, L., Romeu, M., Peinado, I., et al. (2022). Transcriptomic landscape of granulosa cells and peripheral blood mononuclear cells in women with PCOS compared to young poor responders and women with normal response. Hum. Reprod. 37 (6), 1274–1286. doi:10.1093/humrep/deac069
De Rosa, S., Arcidiacono, B., Chiefari, E., Brunetti, A., Indolfi, C., and Foti, D. P. (2018). Type 2 diabetes mellitus and cardiovascular disease: Genetic and epigenetic links. Front. Endocrinol. 9, 2. doi:10.3389/fendo.2018.00002
Donoso-Quezada, J., Ayala-Mar, S., and González-Valdez, J. (2021). The role of lipids in exosome biology and intercellular communication: Function, analytics and applications. Traffic 22, 204–220. doi:10.1111/tra.12803
Drobna, M., Szarzyńska, B., Jaksik, R., Sędek, Ł., Kuchmiy, A., Taghon, T., et al. (2020). hsa-miR-20b-5p and hsa-miR-363-3p affect expression of PTEN and BIM tumor suppressor genes and modulate survival of T-ALL cells in vitro. Cells 9, 1137. doi:10.3390/cells9051137
Duan, X., Wu, Y., Zhang, Z., and Lu, Z. (2020). Identification and analysis of dysregulated lncRNA and associated ceRNA in the pathogenesis of keloid. Ann. Transl. Med. 8, 222. doi:10.21037/atm.2020.01.07
Dudekula, D. B., Panda, A. C., Grammatikakis, I., De, S., Abdelmohsen, K., and Gorospe, M. (2016). CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol. 13, 34–42. doi:10.1080/15476286.2015.1128065
Enright, A. J., John, B., Gaul, U., Tuschl, T., Sander, C., and Marks, D. S. (2003). MicroRNA targets in Drosophila. Genome Biol. 5 (1), R1. doi:10.1186/gb-2003-5-1-r1
Ergun, S., and Oztuzcu, S. (2015). Oncocers: ceRNA-mediated cross-talk by sponging miRNAs in oncogenic pathways. Tumour Biol. 36, 3129–3136. doi:10.1007/s13277-015-3346-x
Escobar-Morreale, H. F., Calvo, R. M., Sancho, J., and San Millán, J. L. (2001). TNF-Α and hyperandrogenism: A clinical, biochemical, and molecular genetic study. J. Clin. Endocrinol. Metab. 86, 3761–3767. doi:10.1210/jcem.86.8.7770
Ezeh, U., Ezeh, C., Pisarska, M. D., and Azziz, R. (2021). Menstrual dysfunction in polycystic ovary syndrome: Association with dynamic state insulin resistance rather than hyperandrogenism. Fertil. Steril. 115, 1557–1568. doi:10.1016/j.fertnstert.2020.12.015
Figliuzzi, M., Marinari, E., and De Martino, A. (2013). MicroRNAs as a selective channel of communication between competing RNAs: A steady-state theory. Biophys. J. 104, 1203–1213. doi:10.1016/j.bpj.2013.01.012
Filetti, V., Loreto, C., Falzone, L., Lombardo, C., Cannizzaro, E., Castorina, S., et al. (2021). Diagnostic and prognostic value of three microRNAs in environmental asbestiform fibers-associated malignant mesothelioma. J. Pers. Med. 11, 1205. doi:10.3390/jpm11111205
Ganie, M. A., Dhingra, A., Nisar, S., Sreenivas, V., Shah, Z. A., Rashid, A., et al. (2016). Oral glucose tolerance test significantly impacts the prevalence of abnormal glucose tolerance among Indian women with polycystic ovary syndrome: Lessons from a large database of two tertiary care centers on the Indian subcontinent. Fertil. Steril. 105, 194–201. doi:10.1016/j.fertnstert.2015.09.005
Giacomini, E., Scotti, G. M., Vanni, V. S., Lazarevic, D., Makieva, S., Privitera, L., et al. (2021). Global transcriptomic changes occur in uterine fluid-derived extracellular vesicles during the endometrial window for embryo implantation. Hum. Reprod. 36, 2249–2274. doi:10.1093/humrep/deab123
Gurung, S., Greening, D. W., Catt, S., Salamonsen, L., and Evans, J. (2020). Exosomes and soluble secretome from hormone-treated endometrial epithelial cells direct embryo implantation. Mol. Hum. Reprod. 26, 510–520. doi:10.1093/molehr/gaaa034
He, T., Sun, Y., Zhang, Y., Zhao, S., Zheng, Y., Hao, G., et al. (2019). MicroRNA-200b and microRNA-200c are up-regulated in PCOS granulosa cell and inhibit KGN cell proliferation via targeting PTEN. Reprod. Biol. Endocrinol. 17, 68. doi:10.1186/s12958-019-0505-8
Huang, H. Y., Lin, Y. C., Li, J., Huang, K. Y., Shrestha, S., Hong, H. C., et al. (2020). miRTarBase 2020: updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 48 (D1), D148–d154. doi:10.1093/nar/gkz896
Jeggari, A., Marks, D. S., and Larsson, E. (2012). miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics 28 (15), 2062–2063. doi:10.1093/bioinformatics/bts344
Jia, C., Wang, S., Yin, C., Liu, L., Zhou, L., and Ma, Y. (2020). Loss of hsa_circ_0118530 inhibits human granulosa-like tumor cell line KGN cell injury by sponging miR-136. Gene 744, 144591. doi:10.1016/j.gene.2020.144591
Kalluri, R., and LeBleu, V. S. (2020). The biology, function, and biomedical applications of exosomes. Science 367, eaau6977. doi:10.1126/science.aau6977
Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 40, D109–D114. doi:10.1093/nar/gkr988
Kawamura, K., Cheng, Y., Kawamura, N., Takae, S., Okada, A., Kawagoe, Y., et al. (2011). Pre-ovulatory LH/hCG surge decreases C-type natriuretic peptide secretion by ovarian granulosa cells to promote meiotic resumption of pre-ovulatory oocytes. Hum. Reprod. 26, 3094–3101. doi:10.1093/humrep/der282
Knebel, B., Janssen, O. E., Hahn, S., Jacob, S., Gleich, J., Kotzka, J., et al. (2008). Increased low grade inflammatory serum markers in patients with Polycystic ovary syndrome (PCOS) and their relationship to PPARgamma gene variants. Exp. Clin. Endocrinol. Diabetes 116, 481–486. doi:10.1055/s-2008-1058085
Legro, R. S., Arslanian, S. A., Ehrmann, D. A., Hoeger, K. M., Murad, M. H., Pasquali, R., et al. (2013). Diagnosis and treatment of polycystic ovary syndrome: An endocrine society clinical practice guideline. J. Clin. Endocrinol. Metab. 98, 4565–4592. doi:10.1210/jc.2013-2350
Li, H., Huang, X., Chang, X., Yao, J., He, Q., Shen, Z., et al. (2020a). S100-A9 protein in exosomes derived from follicular fluid promotes inflammation via activation of NF-κB pathway in polycystic ovary syndrome. J. Cell. Mol. Med. 24, 114–125. doi:10.1111/jcmm.14642
Li, J., and Ge, Z. (2021). High HSPA8 expression predicts adverse outcomes of acute myeloid leukemia. BMC Cancer 21, 475. doi:10.1186/s12885-021-08193-w
Li, J., Ren, L., Li, M., Yang, C., Chen, J., and Chen, Q. (2021b). Screening of potential key genes related to tubal factor infertility based on competitive endogenous RNA network. Genet. Test. Mol. Biomarkers 25 (5), 325–333. doi:10.1089/gtmb.2020.0083
Li, J. Z., Ye, L. H., Wang, D. H., Zhang, H. C., Li, T. Y., Liu, Z. Q., et al. (2021a). The identify role and molecular mechanism of the MALAT1/hsa-mir-20b-5p/TXNIP axis in liver inflammation caused by CHB in patients with chronic HBV infection complicated with NAFLD. Virus Res. 298, 198405. doi:10.1016/j.virusres.2021.198405
Li, S., Zhang, M., Zhang, H., Hu, K., Cai, C., Wang, J., et al. (2020b). Exosomal long noncoding RNA lnc-GNAQ-6:1 may serve as a diagnostic marker for gastric cancer. Clin. Chim. Acta. 501, 252–257. doi:10.1016/j.cca.2019.10.047
Li, S., Zhao, L., and Wan, X. H. (2017). A missense variant rs4645843 in TNF-α gene is a risk factor of polycystic ovary syndrome in the Uygur population. Tohoku J. Exp. Med. 243, 95–100. doi:10.1620/tjem.243.95
Li, X., Guo, L., Zhang, W., He, J., Ai, L., Yu, C., et al. (2022). Identification of potential molecular mechanism related to infertile endometriosis. Front. Vet. Sci. 9, 845709. doi:10.3389/fvets.2022.845709
Liu, D., and Ma, X. (2021). MiR-508-3p promotes proliferation and inhibits apoptosis of middle ear cholesteatoma cells by targeting PTEN/PI3K/AKT pathway. Int. J. Med. Sci. 18, 3224–3235. doi:10.7150/ijms.60907
Liu, H. W., Zhang, F., Fan, P., Bai, H., Zhang, J. X., and Wang, Y. (2013). Effects of apolipoprotein E genotypes on metabolic profile and oxidative stress in southwest Chinese women with polycystic ovary syndrome. Eur. J. Obstet. Gynecol. Reprod. Biol. 170, 146–151. doi:10.1016/j.ejogrb.2013.04.016
Liu, Q., Liu, N., Shangguan, Q., Zhang, F., Chai, W., Tong, X., et al. (2019). LncRNA SAMD12-AS1 promotes cell proliferation and inhibits apoptosis by interacting with NPM1. Sci. Rep. 9, 11593. doi:10.1038/s41598-019-48116-1
Liu, Y., Li, Z., Wang, Y., Cai, Q., Liu, H., Xu, C., et al. (2022). IL-15 participates in the pathogenesis of polycystic ovary syndrome by affecting the activity of granulosa cells. Front. Endocrinol. (Lausanne) 13, 787876. doi:10.3389/fendo.2022.787876
Liu, Y., Liu, H., Li, Z., Fan, H., Yan, X., Liu, X., et al. (2021). The release of peripheral immune inflammatory cytokines promote an inflammatory cascade in PCOS patients via altering the follicular microenvironment. Front. Immunol. 12, 685724. doi:10.3389/fimmu.2021.685724
Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550. doi:10.1186/s13059-014-0550-8
Ma, Y., Ma, L., Cao, Y., and Zhai, J. (2021). Construction of a ceRNA-based lncRNA-mRNA network to identify functional lncRNAs in polycystic ovarian syndrome. Aging 13, 8481–8496. doi:10.18632/aging.202659
Machtinger, R., Laurent, L. C., and Baccarelli, A. A. (2016). Extracellular vesicles: Roles in gamete maturation, fertilization and embryo implantation. Hum. Reprod. Update 22, 182–193. doi:10.1093/humupd/dmv055
Memczak, S., Jens, M., Elefsinioti, A., Torti, F., Krueger, J., Rybak, A., et al. (2013). Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 495, 333–338. doi:10.1038/nature11928
Pegtel, D. M., and Gould, S. J. (2019). Exosomes. Annu. Rev. Biochem. 88, 487–514. doi:10.1146/annurev-biochem-013118-111902
Piao, Z., Zou, R., Lin, Y., Li, Z., Bai, Z., Zhou, L., et al. (2020). Overexpression of lncRNA H19 leads to reduced proliferation in TSCC cells through miR-675-5p/GPR55. Transl. Cancer Res. 9, 891–900. doi:10.21037/tcr.2019.12.45
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, e47. doi:10.1093/nar/gkv007
Rosvall, M., and Bergstrom, C. T. (2010). Mapping change in large networks. PLOS ONE 5, e8694. doi:10.1371/journal.pone.0008694
Rotterdam ESHRE/ASRM-Sponsored PCOS Consensus Workshop Group (2004). Revised 2003 consensus on diagnostic criteria and long-term health risks related to polycystic ovary syndrome (PCOS). Hum. Reprod. 19, 41–47. doi:10.1093/humrep/deh098
Rudnicka, E., Suchta, K., Grymowicz, M., Calik-Ksepka, A., Smolarczyk, K., Duszewska, A. M., et al. (2021). Chronic low grade inflammation in pathogenesis of PCOS. Int. J. Mol. Sci. 22, 3789. doi:10.3390/ijms22073789
Salmena, L., Poliseno, L., Tay, Y., Kats, L., and Pandolfi, P. P. (2011). A ceRNA hypothesis: The rosetta stone of a hidden RNA language? Cell. 146, 353–358. doi:10.1016/j.cell.2011.07.014
San‐Millán, J. L., and Escobar‐Morreale, H. F. (2010). The role of genetic variation in peroxisome proliferator-activated receptors in the polycystic ovary syndrome (PCOS): An original case–control study followed by systematic review and meta-analysis of existing evidence. Clin. Endocrinol. 72, 383–392. doi:10.1111/j.1365-2265.2009.03679.x
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Shao, S., Zhao, H., Lu, Z., Lei, X., and Zhang, Y. (2021). Circadian rhythms within the female HPG axis: From physiology to etiology. Endocrinology 162, bqab117. doi:10.1210/endocr/bqab117
Shreeve, N., Cagampang, F., Sadek, K., Tolhurst, M., Houldey, A., Hill, C. M., et al. (2013). Poor sleep in PCOS; is melatonin the culprit? Hum. Reprod. 28, 1348–1353. doi:10.1093/humrep/det013
Su, Y. Q., Wigglesworth, K., Pendola, F. L., O’Brien, M. J., and Eppig, J. J. (2002). Mitogen-activated protein kinase activity in cumulus cells is essential for gonadotropin-induced oocyte meiotic resumption and cumulus expansion in the mouse. Endocrinology 143, 2221–2232. doi:10.1210/endo.143.6.8845
Sun, Y., Li, S., Liu, H., Bai, H., Hu, K., Zhang, R., et al. (2021). Oxidative stress promotes hyperandrogenism by reducing sex hormone-binding globulin in polycystic ovary syndrome. Fertil. Steril. 116, 1641–1650. doi:10.1016/j.fertnstert.2021.07.1203
Tan, J. Y., Sirey, T., Honti, F., Graham, B., Piovesan, A., Merkenschlager, M., et al. (2015). Extensive microRNA-mediated crosstalk between lncRNAs and mRNAs in mouse embryonic stem cells. Genome Res. 25, 655–666. doi:10.1101/gr.181974.114
Terzieva, D. D., Orbetzova, M. M., Mitkov, M. D., and Mateva, N. G. (2013). Serum melatonin in women with polycystic ovary syndrome. Folia Med. 55, 10–15. doi:10.2478/folmed-2013-0012
Thomson, D. W., and Dinger, M. E. (2016). Endogenous microRNA sponges: Evidence and controversy. Nat. Rev. Genet. 17, 272–283. doi:10.1038/nrg.2016.20
Tu, J., Cheung, A. H. H., Chan, C. L. K., and Chan, W. Y. (2019). The role of microRNAs in ovarian granulosa cells in health and disease. Front. Endocrinol. 10, 174. doi:10.3389/fendo.2019.00174
Tu, M., Wu, Y., Mu, L., and Zhang, D. (2021). Long non-coding RNAs: Novel players in the pathogenesis of polycystic ovary syndrome. Ann. Transl. Med. 9, 173. doi:10.21037/atm-20-5044
Uzbekova, S., Almiñana, C., Labas, V., Teixeira-Gomes, A. P., Combes-Soia, L., Tsikis, G., et al. (2020). Protein cargo of extracellular vesicles from bovine follicular fluid and analysis of their origin from different ovarian cells. Front. Vet. Sci. 7, 584948. doi:10.3389/fvets.2020.584948
Uzdogan, A., Kuru Pekcan, M., Cil, A. P., Kisa, U., and Akbiyik, F. (2021). Progranulin and tumor necrosis factor-alpha in lean polycystic ovary syndrome patients. Gynecol. Endocrinol. 37, 925–929. doi:10.1080/09513590.2021.1958311
Wang, F., Xie, N., Wu, Y., Zhang, Q., Zhu, Y., Dai, M., et al. (2021). Association between circadian rhythm disruption and polycystic ovary syndrome. Fertil. Steril. 115, 771–781. doi:10.1016/j.fertnstert.2020.08.1425
Wang, L., Fan, H., Zou, Y., Yuan, Q., Hu, X., Chen, X., et al. (2020). Aberrant expression of long non-coding RNAs in exosomes in follicle fluid from PCOS patients. Front. Genet. 11, 608178. doi:10.3389/fgene.2020.608178
Wang, X., Chen, X., Meng, Q., Jing, H., Lu, H., Yang, Y., et al. (2015). MiR-181b regulates cisplatin chemosensitivity and metastasis by targeting TGFβR1/Smad signaling pathway in NSCLC. Sci. Rep. 5, 17618. doi:10.1038/srep17618
Wang, Y., He, J., and Yang, J. (2018). Eicosapentaenoic acid improves polycystic ovary syndrome in rats via sterol regulatory element-binding protein 1 (SREBP-1)/toll-like receptor 4 (TLR4) pathway. Med. Sci. Monit. 24, 2091–2097. doi:10.12659/MSM.909098
Wong, N., and Wang, X. (2015). miRDB: an online resource for microRNA target prediction and functional annotations. Nucleic Acids Res. 43, D146–D152. doi:10.1093/nar/gku1104
Wu, Q., Gao, J., Bai, D., Yang, Z., and Liao, Q. (2021). The prevalence of polycystic ovarian syndrome in Chinese women: A meta-analysis. Ann. Palliat. Med. 10, 74–87. doi:10.21037/apm-20-1893
Xia, S., Feng, J., Chen, K., Ma, Y., Gong, J., Cai, F., et al. (2018). Cscd: A database for cancer-specific circular RNAs. Nucleic Acids Res. 46, D925–D929. doi:10.1093/nar/gkx863
Xie, F., Zhang, J., Zhai, M., Liu, Y., Hu, H., Yu, Z., et al. (2021). Melatonin ameliorates ovarian dysfunction by regulating autophagy in PCOS via the PI3K-Akt pathway. Reproduction 162, 73–82. doi:10.1530/REP-20-0643
Xiong, Y., Pang, M., Du, Y., Yu, X., Yuan, J., Liu, W., et al. (2022). The LINC01929/miR-6875-5p/ADAMTS12 axis in the ceRNA network regulates the development of advanced bladder cancer. Front. Oncol. 12, 856560. doi:10.3389/fonc.2022.856560
Yao, R. W., Wang, Y., and Chen, L. L. (2019). Cellular functions of long noncoding RNAs. Nat. Cell. Biol. 21, 542–551. doi:10.1038/s41556-019-0311-8
Yu, B. X., Zou, L., Li, S., and Du, Y. L. (2019). LncRNA SAMD12-AS1 down-regulates P53 to promote malignant progression of glioma. Eur. Rev. Med. Pharmacol. Sci. 23, 8456–8467. doi:10.26355/eurrev_201910_19158
Yu, G., Yao, W., Gumireddy, K., Li, A., Wang, J., Xiao, W., et al. (2014). Pseudogene PTENP1 functions as a competing endogenous RNA to suppress clear-cell renal cell carcinoma progression. Mol. Cancer Ther. 13, 3086–3097. doi:10.1158/1535-7163.MCT-14-0245
Zhang, Q., Ren, J., Wang, F., Pan, M., Cui, L., Li, M., et al. (2022). Mitochondrial and glucose metabolic dysfunctions in granulosa cells induce impaired oocytes of polycystic ovary syndrome through Sirtuin 3. Free Radic. Biol. Med. 187, 1–16. doi:10.1016/j.freeradbiomed.2022.05.010
Zhang, Y., Ran, Y., Kong, L., Geng, L., Huang, H., Zhang, H., et al. (2021). Decreased SFRP5 correlated with excessive metabolic inflammation in polycystic ovary syndrome could be reversed by metformin: Implication of its role in dysregulated metabolism. J. Ovarian Res. 14, 97. doi:10.1186/s13048-021-00847-4
Zhao, R., Jiang, Y., Zhao, S., and Zhao, H. (2021). Multiomics analysis reveals molecular abnormalities in granulosa cells of women with polycystic ovary syndrome. Front. Genet. 12, 648701. doi:10.3389/fgene.2021.648701
Zhao, X., Xu, F., Qi, B., Hao, S., Li, Y., Li, Y., et al. (2014). Serum metabolomics study of polycystic ovary syndrome based on liquid chromatography–mass spectrometry. J. Proteome Res. 13, 1101–1111. doi:10.1021/pr401130w
Zhou, R., Li, S., Liu, J., Wu, H., Yao, G., Sun, Y., et al. (2020). Up-regulated FHL2 inhibits ovulation through interacting with androgen receptor and ERK1/2 in polycystic ovary syndrome. EBiomedicine 52, 102635. doi:10.1016/j.ebiom.2020.102635
Keywords: ceRNA network, exosomes, granulosa cells, polycystic ovary syndrome, post-transcriptional regulation
Citation: Bai L, Gong J, Guo Y, Li Y, Huang H and Liu X (2022) Construction of a ceRNA network in polycystic ovary syndrome (PCOS) driven by exosomal lncRNA. Front. Genet. 13:979924. doi: 10.3389/fgene.2022.979924
Received: 28 June 2022; Accepted: 17 October 2022;
Published: 04 November 2022.
Edited by:
Jing-Woei Li, Prince of Wales Hospital, ChinaReviewed by:
Yufang Liu, Hebei University of Engineering, ChinaXiaoyong Li, Shanghai Jiao Tong University, China
Jin Haixia, First Affiliated Hospital of Zhengzhou University, China
Jin Yu, Shanghai Jiao Tong University, China
Copyright © 2022 Bai, Gong, Guo, Li, Huang and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinmei Liu, liuxinmei@fudan.edu.cn, nanlilac@hotmail.com
†These authors have contributed equally to this work and share first authorship