Integrated bioinformatics analysis elucidates granulosa cell whole-transcriptome landscape of PCOS in China
Journal of Ovarian Research volume 16, Article number: 154 (2023)
Polycystic ovary syndrome (PCOS) is a common reproductive, neuroendocrine, and metabolic disorder in women of reproductive age that affects up to 5–10% of women of reproductive age. The aetiology of follicle development arrest and critical issues regarding the abnormal follicular development in PCOS remain unclear. The present study aims to systematically evaluate granulosa cell whole-transcriptome sequencing data to gain more insights into the transcriptomic landscape and molecular mechanism of PCOS in China.
In the present study, the microarray datasets GSE138518, GSE168404, GSE193123, GSE138572, GSE95728, and GSE145296 were downloaded from the Gene Expression Omnibus (GEO) database. Subsequently, differential expression analysis was performed on the PCOS and control groups, followed by functional interaction prediction analysis to investigate gene-regulatory circuits in PCOS. Finally, hub genes and their associated ncRNAs were validated by qPCR in human-luteinized granulosa (hGL) cells and were correlated with the clinical characteristics of the patients.
A total of 200 differentially expressed mRNAs, 3 differentially expressed miRNAs, 52 differentially expressed lncRNAs, and 66 differentially expressed circRNAs were found in PCOS samples compared with controls. GO and KEGG enrichment analyses indicated that the DEGs were mostly enriched in phospholipid metabolic processes, steroid biosynthesis and inflammation related pathways. In addition, the upregulated miRNA hsa-miR-205-5p was significantly enriched in the ceRNA network, and two hub genes, MVD and PNPLA3, were regulated by hsa-miR-205-5p, which means that hsa-miR-205-5p may play a fundamental role in the pathogenesis of PCOS. We also found that MVD and PNPLA3 were related to metabolic processes and ovarian steroidogenesis, which may be the cause of the follicle development arrest in PCOS patients.
In summary, we systematically constructed a ceRNA network depicting the interactions between the ncRNAs and the hub genes in PCOS and control subjects and correlated the hub genes with the clinical characteristics of the patients, which provides valuable insights into the granulosa cell whole-transcriptome landscape of PCOS in China.
Polycystic ovary syndrome (PCOS) is a common reproductive, neuroendocrine, and metabolic disorder in women of reproductive age that affects up to 5–10% of reproductive aged women . Its main clinical manifestations are ovulatory dysfunction, hyperandrogenemia, and polycystic ovaries, which can lead to infertility [2, 3]. Additional clinical features include metabolic abnormalities such as insulin resistance, obesity and type 2 diabetes (T2D) . Although research on the aetiology of follicle development arrest in patients with PCOS has continuously emerged in recent years, critical issues regarding abnormal follicular development in PCOS and the precise network regulation mechanism remain unclear.
Granulosa cells (GCs) are an important cell type surrounding follicles that can interfere with follicle maturation and ovulation . Several studies have shown that GC dysfunction is associated with the disruption of follicle development, such as excessive follicular recruitment, obstruction of dominant follicular selection, follicular atresia, and anovulation and metabolic disorder in PCOS [6,7,8,9,10]. Folliculogenesis and steroidogenesis rely on oocyte-GC crosstalk, which can provide oocytes with nutrients and the removal of waste. The communication between oocytes and GCs is critical for normal follicular development and the secretion of steroid hormones by GCs . Moreover, studies have suggested that the metabolic process of GCs has an effect on follicular development [12, 13]. Until now, how GC dysfunction leads to the occurrence of PCOS has not been elucidated.
In recent years, an increasing number of studies have performed transcriptome sequencing of GCs to explore the pathogenesis of PCOS [14,15,16,17]. Competing endogenous RNAs (ceRNAs), natural decoys that compete for a common pool of microRNAs (miRNAs), represent a novel layer of gene regulation by systematically functionalizing miRNA response element (MRE)-harboring noncoding RNAs, such as long noncoding RNAs (lncRNAs), pseudogenes, and circular RNAs (circRNAs), and forming complex miRNA-mediated ceRNA networks [18,19,20,21]. Perturbation of the ceRNA crosstalk balance of cellular processes and functions leads to diseases such as PCOS [22,23,24]. Recently, a study reported that lncXIST inhibited human granulosa-like tumor cell viability and induced apoptosis by increasing the expression of Bcl2-like protein 11 (BCL2L11) via the sponging of miR-30c-5p by ceRNA . However, few studies have utilized GCs whole-transcriptome sequencing strategies, which allow accurate examination of global gene expression profiles, to elucidate the GCs transcriptome characteristics of PCOS in China.
In this study, we investigated the whole-transcriptome profiles in GCs of PCOS patients by using GEO datasets. Subsequently, differential mRNA, miRNA, lncRNA and circRNA expression analyses were performed between the control groups and PCOS groups, followed by functional interaction prediction analysis. The results showed that PNPLA3, MVD, MMP9, LCK, NCF1, OSM, C3, MLXIPL and TREM1 were closely related to the onset of PCOS. In addition, miR-205-5p, miR-210-5p, and miR-144-5p and their associated lncRNAs and circRNAs were also involved in the progression of PCOS. This study uncovered a reliable molecular basis of PCOS initiation and progression, and provided clues to investigate the onset and development of PCOS.
The National Center for Biotechnology Information (NCBI)’s Gene Expression Omnibus (GEO) repository (https://www.ncbi.nlm.nih.gov/geo/) was searched for gene expression profile datasets of GCs from women with or without PCOS. Gene expression datasets for 6 PCOS studies were downloaded (GSE138518, GSE168404, GSE193123, GSE138572, GSE95728, and GSE145296) from Gene Expression Omnibus (GEO, https://www.ncbi.nlm.nih.gov/geo/), including 34 control (control group), and 34 PCOS patients (case group) from China, as shown in Table 1. Background correction, and quantile normalization, were performed by RMA normalization.
Identification of differentially expressed mRNAs and noncoding RNAs
Differential gene analysis was performed by edgeR and the limma package. We defined RNAs with a P value < 0.05 and |log FC| > 0.5 as being differentially expressed. Overlapping differentially expressed mRNAs were identified in at least two mRNA datasets (GSE138518, GSE168404, and GSE193123). Overlapping differentially expressed miRNAs were identified in two miRNA datasets (GSE138572 and GSE168404). Furthermore, we identified differentially expressed lncRNAs and circRNAs in the lncRNA dataset (GSE95728) and circRNA dataset (GSE145296).
Protein-Protein Interaction (PPI) network, module extraction and hub genes identification
To further investigate the function of dif-mRNAs at the protein level, we constructed a PPI network using the STRING database (https://string-db.org/) and visualized it by Cytoscape. Interactions with a combined score ≥ 0.4 were considered statistically significant. Then, to identify highly interacting hub mRNA clustering, we established “Molecular Complex Detection” (MCODE), a clustering algorithm that identifies locally densely connected regions in a large PPI network based on node-weighting arithmetic with degree cut-off score = 2, k-core = 2 and max depth = 100. The cytoHubba app in Cytoscape was used to disclose the hub genes in the PPI network. In the whole PPI network, the top 50 hub genes ranked by maximal clique centrality (MCC) were obtained.
Functional annotation and enrichment analysis
Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were conducted to investigate the roles of the differentially expressed mRNAs. GO enrichment analysis and KEGG enrichment analysis were carried out using clusterProfiler. A GO tree was constructed to summarize the affected functions. The mutual regulatory relationships between enriched KEGG pathways were illustrated by Pathway-Act networks.
Competing endogenous RNA (ceRNA) network construction
The predicted lncRNA‒miRNA pairs and miRNA‒mRNA pairs were collected from the miRcode and TargetScan databases. The lncRNA‒miRNA‒mRNA networks were visualized by Cytoscape. The predicted circRNA‒miRNA pairs were collected from the circbank database. The circRNA‒miRNA‒mRNA networks were visualized by Cytoscape. Furthermore, according to the lncRNA‒miRNA‒mRNA and circRNA‒miRNA‒mRNA networks, differentially expressed circRNAs, lncRNAs, and mRNAs that were regulated by the same miRNA were further screened to construct a ceRNA Network.
Real-time quantitative PCR analysis
Two DEGs, four lncRNAs, four circRNAs and one miRNA from the ceRNA network were selected, and their differential expression was validated by real-time PCR. Total RNA was derived from human-luteinized granulosa (hGL) cells of PCOS and control. PCOS was diagnosed according to Rotterdam’s diagnostic criteria, which required the presence of any two of the following three conditions while excluding other causes of excessive androgen: 1) Oligo-ovulation or anovulation; 2) clinical manifestations of high androgen levels or hyperandrogenism; 3) polycystic ovaries identified by ultrasound, with one ovary or bilateral ovary showing a diameter of 2-9 mm follicles ≥ 12, or ovarian volume ≥ 10ml. The control group consisted of patients with tubal infertility. The reverse transcription of RNA was performed with the PrimeScript RT Reagent Kit (Perfect Real Time, Takara), following the manufacturer’s instructions. Then, qPCR was conducted to amplify cDNA samples using SYBR Green PCR Master Mix (#MR101; Vazyme, China; Q711; Vazyme, China). The expression levels of mRNAs, circRNAs and lncRNAs were normalized to GAPDH, and miRNA expression levels were normalized to U6 and analysed using the 2-ΔΔCt method. Sequences of the specific real-time PCR primers for selected genes and ncRNAs were presented in Supplemental Table 4.
A GC tumor-derived cell line (KGN) was utilized to explore the miRNA‒mRNA network. KGN cells were seeded into six-well plates at a density of 2 × 105 cells/well in DMEM supplemented with 5% foetal bovine serum. All the cell models used in this study were cultured in a humidified atmosphere of 5% CO2 and 95% air at 37 °C, and the cell culture medium was changed every 2 days in all experiments.
Cells were cultured to approximately 70% density and then transfected with 100 nM hsa-miR-205-5p mimics and negative control (RIBBIO, Guangzhou, China) using Lipofectamine RNA iMAX according to the manufacturer’s instructions (Life Technologies). The efficiency of the hsa-miR-205-5p mimics was detected by quantitative real-time RT‒PCR.
Western blot analysis
After the treatment, all the cells were lysed in cell lysis buffer (Cell Signaling Technology), and the protein concentration was determined using a Pierce Rapid Gold BCA kit following the manufacturer’s instructions (Thermo Fisher, USA). Equal amounts of protein were loaded and separated using sodium dodecyl sulfate‒polyacrylamide gel electrophoresis (SDS‒PAGE) analysis. Afterwards, the proteins were transferred onto polyvinylidene difluoride (PVDF) membranes (Bio-Rad, USA), followed by blocking with Tris-buffered saline (TBS) containing 5% nonfat dry milk for 1 h at room temperature and incubated overnight at 4 °C with corresponding primary antibodies. The next day, the membranes were washed with TBS for 1 h and then incubated in the appropriate HRP-conjugated secondary antibody for 30 min. Similarly, the membranes were washed with TBS for 1 h after secondary antibody incubation. Finally, the immunoreactive bands were detected using an enhanced chemiluminescent substrate (Bio-Rad) and X-ray film. The intensities of the bands were quantified with Image-Pro Plus software (v.4.5; Media Cybernetics, USA).
All experiments were repeated three times. Data analysis was performed using GraphPad Prism version 9 (GraphPad Software Inc, San Diego, CA, USA). All data are presented as the mean ± standard error of the mean (SEM). The Mann‒Whitney U test or unpaired Student’s t test was performed to compare mRNA, lncRNA, circRNA and miRNA expression between groups. Statistical significance was defined as a two-tailed P value < 0.05.
Differential expression analysis and functional enrichment analysis
According to the screening criteria, a total of 286 overlapping differentially expressed mRNAs (dif-mRNAs) in at least two mRNA datasets were obtained (Fig. 1A). Six overlapping differentially expressed miRNAs (dif-miRNAs) were identified (Fig. 1B). A total of 4256 differentially expressed lncRNAs (dif-lncRNAs) were obtained, of which 2634 were upregulated and 1622 were downregulated (Fig. 1C). A total of 3614 differentially expressed circRNAs (dif-circRNAs) were obtained, of which 1823 were upregulated and 1791 were downregulated (Fig. 1D). Then, functional enrichment analysis of dif-mRNAs was performed and the top 30 GO terms or KEGG pathways were displayed. We found that the GO terms that differentially expressed genes were mainly involved in were the regulation of phospholipid metabolic process, regulation of phosphatidylcholine metabolic process, and hydrogen peroxide catabolic process (Fig. 1E). Additionally, the KEGG pathways that the differentially expressed genes were significantly associated with included steroid biosynthesis, fluid shear stress and atherosclerosis, and complement and coagulation cascades (Fig. 1F).
Protein-Protein Interaction (PPI) network, module extraction and hub genes identification
The PPI network based on dif-mRNA consisted of 202 nodes and 630 interaction pairs (Fig. 2A). Nodes with high topological scores can be regarded as key nodes of the network. Using the Cytoscape plug-in MCODE (score ≥ 4), four subnetwork modules were aggregated and extracted from the PPI network (Fig. 2B). Module A (score = 6) contained 9 nodes and 24 interaction pairs, in which inflammatory factors were most included, such as interleukin 10 (IL 10, degree = 70), LCK proto-oncogene, Src family tyrosine kinase (LCK, degree = 36), and CD14 molecule (CD14, degree = 20). Module B (score = 4.8) contained 21 nodes and 48 interaction pairs, in which the genes were associated with lipid metabolism, such as patatin-like phospholipase domain containing 3 (PNPLA3, degree = 12), mevalonate diphosphate decarboxylase (MVD, degree = 24), and fatty acid desaturase 2 (FADS2, degree = 20). Module C (score = 4) contained 4 nodes and 6 interaction pairs, including acyl-CoA synthetase short chain family member 2 (ACSS2, degree = 26) and stearoyl-CoA desaturase (SCD, degree = 32). Module D (score = 4) contained 5 nodes and 8 interaction pairs, including matrix metallopeptidase 9 (MMP9, degree = 68) and interleukin 1 beta (IL 1B, degree = 106). Moreover, genes in the modules were subjected to GO enrichment analysis. According to the significance order, the top 10 terms for each module were selected for display (Fig. 2C). The genes in Module A were significantly involved in the chronic inflammatory response, Toll-like receptor binding, and RAGE receptor binding. Genes in Module B were significantly associated with cholesterol catabolic process and sterol catabolic process. Genes in Module C were concerned with oxidoreductase activity and cellular lipid biosynthetic process. Genes in Module D were enriched in superoxide-generating NAD(P)H oxidase activity and immature T-cell proliferation in thymus.
To further identify the hub genes, we identified hub mRNAs by cytoHubba, a plug in Cytoscape, and the top 50 hub genes were identified and selected from the PPI network (Fig. 3A). GO enrichment analysis showed that hub mRNAs were mainly enriched in cholesterol metabolic process, sterol metabolic process, and regulation of inflammatory response (Fig. 3B). GO-Tree analysis of GO biological process (GO-BP) terms was performed to further understand the core BPs associated with PCOS based on their subordinate and interaction relationships (Fig. 3C). A Pathway-Act network was constructed to further investigate the mutual interactions of pathways and to obtain the hub pathways that may play a vital role in PCOS (Fig. 3D). The top pathways that showed interactions with other surrounding pathways were the Toll-like receptor signaling pathway, cholesterol metabolism, ovarian steroidogenesis, and the AMPK signaling pathway. These results indicated that lipid metabolism, especially the cholesterol metabolic pathway, and regulation of the inflammatory response, especially the Toll-like receptor signaling pathway, might play key roles in PCOS initiation and progression.
ceRNA network construction
Based on the regulatory relationship of dif-miRNA–dif-mRNA and dif-miRNA–dif-lncRNA, the lncRNAs and mRNAs that were significantly differentially expressed and regulated by the same miRNA were screened. In total, 457 lncRNA–miRNA–mRNA interactions were finally obtained (Fig. 4 and Supplement Table 1), including 50 upregulated and 5 downregulated lncRNAs, 91 upregulated and 110 downregulated mRNAs, and 2 upregulated and 2 downregulated miRNAs.
Based on the regulatory relationship of dif-miRNA–dif-mRNA and dif-miRNA–dif-circRNA, dif-circRNA and mRNA regulated by the same miRNAs were screened, resulting in 707 interaction relationships of circRNA–miRNA–mRNA. There were 39 upregulated and 49 downregulated circRNAs, 117 upregulated mRNAs and 101 downregulated mRNAs, and 2 upregulated miRNAs and 4 downregulated miRNAs. The circRNA–miRNA–mRNA network is shown in Fig. 5 and Supplemental Table 2.
Furthermore, according to the lncRNA–miRNA–mRNA and circRNA–miRNA–mRNA networks, differentially expressed circRNAs, lncRNAs, and mRNAs that were regulated by the same miRNA were further screened. Finally, 504 interaction pairs were obtained (Fig. 6 and Supplement Table 3), of which 31 upregulated and 35 downregulated circRNAs, 48 upregulated and 4 downregulated lncRNAs, 91 upregulated and 109 downregulated mRNAs, and 1 upregulated and 2 downregulated miRNAs (hsa-miR-144-5p, and hsa-miR-210-5p, downregulated; hsa-miR-205-5p, upregulated) were included.
To further narrow down the hub dif-mRNAs to identify those that are the core dif-mRNAs, we merged the hub dif-mRNAs from the PPI network, MCODE cluster mRNAs, 11 overlapping mRNAs, and dif-mRNAs in the ceRNA network and found two core dif-mRNAs (Fig. 7A). These core mRNAs were all downregulated in PCOS and regulated by hsa-miR-205-5p, including the triacylglycerol lipase PNPLA3 and cholesterol biosynthesis enzyme MVD. To fully validate these most differentiated mRNAs, lncRNAs, miRNAs and circRNAs, we performed RT–qPCR using hGL cells from Women’s Hospital, School of Medicine, Zhejiang University. The primers for real-time quantitative PCR are shown in Supplemental Table 4. The baseline characteristics of the patients are shown in Table 2. AFC, BMI, duration of infertility, LH level, ratio of LH to FSH, T level, AMH, triglyceride, Total cholesterol, LDL-C, HOMA-IR, number of oocytes retrieved, and number of MII oocytes were higher in the PCOS group than in the control group. The number of 2PN fertilized oocytes was lower in the PCOS group. We observed downregulation of PNPLA3 and MVD (Fig. 7B and C), and upregulation of hsa-miR-205-5p in hGL (Fig. 7D). Moreover, transfection of KGN cells with hsa-miR-205-5p mimics led to downregulation of PNPLA3 and MVD at both mRNA and protein levels, supporting the regulatory relationship between hsa-miR-205-5p and PNPLA3 and MVD (Fig. 7E-I). Additionally, we found that two lncRNAs (SAA2 and SLC30A10) and four circRNAs (hsa-circ-0020555, hsa-circ-0027651, hsa-circ-0086809, and hsa-circ-0118448) which regulate hsa-miR-205-5p were all downregulated in hGL (Fig. 7J-M, N-Q).
To further explore the relationship between PNPLA3 and MVD and the occurrence and progression of PCOS, we correlated the mRNA expression of these two genes with the clinical characteristics of patients (Table 3). We found that PNPLA3 was negatively correlated with AFC, number of oocytes retrieved, number of MII oocytes, LH levels, and AMH levels, and positively correlated with the rate of day 3 good-quality embryos, P levels and HDL-C levels (Fig. 8A-H, Figure S1). MVD was negatively correlated with total cholesterol levels, LDL-C levels and HDL-C levels, and positively correlated with the rate of day 3 good-quality embryos, LH levels, estradiol levels and triglyceride levels (Fig. 9A-H, Figure S2). MVD is a key enzyme in cholesterol synthesis, and cholesterol is an important precursor for steroid hormone synthesis, hence, we speculate that reduced MVD levels in PCOS affect cholesterol synthesis and further affect the synthesis of steroid hormones such as estrogen, which in turn affects oocyte development, but the mechanism still needs to be further explored.
The current study aimed to identify common DEGs of PCOS from multiple datasets and to identify the potential genes responsible for the occurrence and progression of PCOS. In the present study, we analysed the whole-transcriptome characteristics of PCOS GCs. The GO and KEGG pathway analyses showed that cholesterol catabolic process, sterol catabolic process, and chronic inflammatory response were significantly enriched. To further explore the interactions between mRNA and ncRNA during the pathogenesis of PCOS, we constructed PPI and ceRNA networks. It turned out that mRNA (PNPLA3 and MVD), miRNA (hsa-miR-205-5p, hsa-miR-210-5p, has-and miR-144-5p), and miRNA-associated lncRNAs and circRNAs play key roles in the initiation and progression of PCOS. Overall, we delineated the GCs whole-transcriptomic landscape and identified possible changes in ovarian steroidogenesis, metabolic disorders, and immune dysfunction in PCOS.
The cross-talk between GCs and follicles plays a key role in follicular development and maturation . PCOS patients show blocked follicle development in the small sinus follicular phase and polycystic ovaries . In the early stage of follicle development, follicles are mainly composed of immature oocytes and GCs surrounding them . Studies have shown that in the process of primordial follicle activation and follicle development, the metabolism of oocytes and surrounding GCs is significantly enhanced [28, 29]. Metabolism disorders such as lipid metabolism disorders or hyperandrogenemia can cause follicular development arrest . In the current study, we found that metabolic processes and ovarian steroidogenesis were significantly impaired. Cholesterol catabolic processes and sterol catabolic processes were significantly enriched in PCOS GCs. Among the genes we identified through coexpression of three datasets, PNPLA3, MVD, MMP9, oncostatin M (OSM), LCK, triggering receptor expressed on myeloid cells 1 (TREM1), FADS2, proprotein convertase subtilisin/kexin type 9 (PCSK9), and C3 are involved in lipid metabolism and ovarian dysfunction. FADS2 is decreased in the transcriptome sequencing of PCOS patients with androgen excess and regulates lipid metabolism of PCOS patients [14, 31]. In addition, PNPLA3 may regulate oocyte development capacity in PCOS . Abnormal high expression of proprotein PCSK9 may be involved in pathogenesis of PCOS by affecting lipid metabolism and ovarian function, and the inhibition of PCSK9 partly reverses the pathological changes in PCOS [33,34,35]. MMP9 and TREM1 were identified as hub genes in PCOS patients with non-alcoholic fatty liver (NAFLD). Upregulation of MMP9 could contribute to excess abdominal adiposity and hyperandrogenism, which might be related to increased cardiovascular risk in PCOS [36,37,38]. Metformin could decrease MMP9 levels to improve the lipid metabolism of PCOS [39, 40]. OSM, a recently identified adipokine, could negatively affect oocyte maturation and fertilization rates . The tyrosine kinase LCK may have an effect on PCOS patients with insulin resistance (IR) and may play an important role in PCOS pathogenesis . On all accounts, our results showed that abnormal lipid metabolism was closely related to the occurrence and development of PCOS; however, the specific mechanism underlying how these metabolic disorders contribute to ovulatory disorders and follicular atresia in PCOS remains elusive. These genes involved in metabolic disorders in PCOS GCs modulated cholesterol synthesis and ovarian steroidogenesis, which indicated that fatty acid and cholesterol metabolism disorders in PCOS GCs may contribute to the impairment in oocyte maturation.
In recent years, several studies have shown that women with PCOS present with chronic low-grade inflammation, indicating abnormal expression of the proinflammatory cytokines interleukin-1 (IL-1), CRP, and interleukin-18 (IL-18), as well as endothelial dysfunction and increased oxidative stress [43,44,45,46]. Anti-inflammatory therapy can improve the symptoms of PCOS [47, 48]. We found that the chronic inflammatory response was enriched and that IL-18 production was increased, which indicated an immune imbalance in PCOS patients. In addition, C3 was involved in the regulation of inflammation. Higher C3 levels also had a stronger association with IR and might be an inflammatory marker of IR in women with PCOS [49, 50]. Women with PCOS exhibited higher MMP9 levels, which could contribute to chronic low-grade inflammation and an atherothrombotic state . Moreover, in PCOS offspring, MMP9 was also increased, suggesting that these children may exhibit increased chronic low-grade inflammation . In the physiological state, GCs excrete prostaglandin E2 (PGE2) and some inflammatory cytokines and chemokines to promote ovulation . These genes are involved in immune regulation to maintain the immune balance of GCs in PCOS patients, which could support follicle development and ovulation. The dysregulation of genes associated with immune regulation could be a molecular clue regarding the abnormal folliculogenesis and anovulation in PCOS patients.
The circRNA dataset (GSE145296) was used cumulus cells (CCs), while the others were used mural granulosa cells (MGCs). This dataset is currently the only one in the GEO database performing circRNA sequencing on granulosa cells from PCOS patients. CircRNAs play an essential role in whole transcriptome analysis and have significant physiological functions, such as regulating mRNA expression by acting as competitive endogenous RNAs for miRNAs, which can impact cellular function and contribute to PCOS development [21, 22, 53, 54]. In addition, undifferentiated GCs differentiate into MGCs and CCs during follicle antrum formation based on their location distribution . The oocyte within the follicle is surrounded by CCs, which connect with MGCs that are responsible for secreting steroid hormones. Crosstalk between oocytes and surrounding somatic cells is crucial for mammalian oogenesis [25, 56]. Both cell types originate from undifferentiated GCs, the transcriptomes of these cell types differ before ovulation [57, 58], but they converge after ovulation such that their gene expression profiles become completely overlapping. Notably, CCs transition to be more like MGCs after human chorionic gonadotropin (HCG) . The CCs used in the circRNA dataset (GSE145296) was retrieved 36 h after HCG administration. Hence, the CCs and MGCs could have similar gene expression profiles. Thus, we kept this dataset for integrated analysis. We also validated the differential circRNAs identified in the dataset using hGL. The results showed that all four circRNAs were downregulated. Thus, using this dataset of CCs for complete transcriptome analysis provides valuable insights into changes in the whole transcriptome of GCs in PCOS.
In summary, this study found several DEGs related to metabolic processes, ovarian steroidogenesis and the immune response, including MVD, PNPLA3, MMP9, OSM, and C3. Additionally, this study identified some ncRNAs, especially hsa-miR-205-5p and its associated circRNAs and lncRNAs. The relationship among these molecules may contribute to the onset and development of PCOS. We delineated the GC whole-transcriptomic landscape and provided a valuable direction aiming to improve the fertility of PCOS in China.
Availability of data and materials
The microarray datasets GSE138518, GSE168404, GSE193123, GSE138572, GSE95728, GSE145296 were downloaded from GEO database (https://www.ncbi.nlm.nih.gov/geo/).
Polycystic ovary syndrome
In vitro fertilization
Intracytoplasmic sperm injection
Fertilized oocytes with two primary pronucleus
Follicle stimulating hormone
Antral follicle count
Body mass index
Homeostatic model assessment of insulin resistance
Human chorionic gonadotropin
Mural granulosa cells
Norman RJ, Dewailly D, Legro RS, Hickey TE. Polycystic ovary syndrome. Lancet. 2007;370(9588):685–97. https://doi.org/10.1016/s0140-6736(07)61345-2.
Fauser BC, Tarlatzis BC, Rebar RW, Legro RS, Balen AH, Lobo R, et al. Consensus on women’s health aspects of polycystic ovary syndrome (PCOS): the Amsterdam ESHRE/ASRM-Sponsored 3rd PCOS Consensus Workshop Group. Fertil Steril. 2012;97(1):28-38.e25. https://doi.org/10.1016/j.fertnstert.2011.09.024.
Tata B, Mimouni NEH, Barbotin AL, Malone SA, Loyens A, Pigny P, et al. Elevated prenatal anti-Müllerian hormone reprograms the fetus and induces polycystic ovary syndrome in adulthood. Nat Med. 2018;24(6):834–46. https://doi.org/10.1038/s41591-018-0035-5.
Azziz R, Carmina E, Chen Z, Dunaif A, Laven JS, Legro RS, et al. Polycystic ovary syndrome. Nat Rev Dis Primers. 2016;2:16057. https://doi.org/10.1038/nrdp.2016.57.
Li X, Zhu Q, Wang W, Qi J, He Y, Wang Y, et al. Elevated chemerin induces insulin resistance in human granulosa-lutein cells from polycystic ovary syndrome patients. FASEB J. 2019;33(10):11303–13. https://doi.org/10.1096/fj.201802829R.
Jonard S, Dewailly D. The follicular excess in polycystic ovaries, due to intra-ovarian hyperandrogenism, may be the main culprit for the follicular arrest. Hum Reprod Update. 2004;10(2):107–17. https://doi.org/10.1093/humupd/dmh010.
Franks S, Stark J, Hardy K. Follicle dynamics and anovulation in polycystic ovary syndrome. Hum Reprod Update. 2008;14(4):367–78. https://doi.org/10.1093/humupd/dmn015.
El-Hayek S, Yang Q, Abbassi L, FitzHarris G, Clarke HJ. Mammalian oocytes locally remodel follicular architecture to provide the foundation for germline-soma communication. Curr Biol. 2018;28(7):1124-1131e3. https://doi.org/10.1016/j.cub.2018.02.039.
Li R, Albertini DF. The road to maturation: somatic cell interaction and self-organization of the mammalian oocyte. Nat Rev Mol Cell Biol. 2013;14(3):141–52. https://doi.org/10.1038/nrm3531.
Ding L, Gao F, Zhang M, Yan W, Tang R, Zhang C, et al. Higher PDCD4 expression is associated with obesity, insulin resistance, lipid metabolism disorders, and granulosa cell apoptosis in polycystic ovary syndrome. Fertil Steril. 2016;105(5):1330-1337e3. https://doi.org/10.1016/j.fertnstert.2016.01.020.
Chang HM, Qiao J, Leung PC. Oocyte-somatic cell interactions in the human ovary-novel role of bone morphogenetic proteins and growth differentiation factors. Hum Reprod Update. 2016;23(1):1–18. https://doi.org/10.1093/humupd/dmw039.
Gioacchini G, Notarstefano V, Sereni E, Zacà C, Coticchio G, Giorgini E, et al. Does the molecular and metabolic profile of human granulosa cells correlate with oocyte fate? New insights by Fourier transform infrared microspectroscopy analysis. Mol Hum Reprod. 2018;24(11):521–32. https://doi.org/10.1093/molehr/gay035.
Yilmaz B, Vellanki P, Ata B, Yildiz BO. Metabolic syndrome, hypertension, and hyperlipidemia in mothers, fathers, sisters, and brothers of women with polycystic ovary syndrome: a systematic review and meta-analysis. Fertil Steril. 2018;109(2):356-364e32. https://doi.org/10.1016/j.fertnstert.2017.10.018.
Liao B, Qi X, Yun C, Qiao J, Pang Y. Effects of Androgen excess-related metabolic disturbances on Granulosa cell function and Follicular Development. Front Endocrinol (Lausanne). 2022;13:815968. https://doi.org/10.3389/fendo.2022.815968.
Zhao R, Jiang Y, Zhao S, Zhao H. Multiomics analysis reveals molecular abnormalities in granulosa cells of women with polycystic ovary syndrome. Front Genet. 2021;12: 648701. https://doi.org/10.3389/fgene.2021.648701.
Li J, Chen H, Gou M, Tian C, Wang H, Song X, et al. Molecular features of polycystic ovary syndrome revealed by transcriptome analysis of oocytes and cumulus cells. Front Cell Dev Biol. 2021;9: 735684. https://doi.org/10.3389/fcell.2021.735684.
Zhou J, Huang X, Xue B, Wei Y, Hua F. Bioinformatics analysis of the molecular mechanism of obesity in polycystic ovary syndrome. Aging. 2021;13(9):12631–40. https://doi.org/10.18632/aging.202938.
Zhao Y, Tao M, Wei M, Du S, Wang H, Wang X. Mesenchymal stem cells derived exosomal mir-323-3p promotes proliferation and inhibits apoptosis of cumulus cells in polycystic ovary syndrome (PCOS). Artif Cells Nanomed Biotechnol. 2019;47(1):3804–13. https://doi.org/10.1080/21691401.2019.1669619.
Geng X, Zhao J, Huang J, Li S, Chu W, Wang WS, et al. lnc-MAP3K13-7:1 inhibits ovarian GC proliferation in PCOS via DNMT1 downregulation-mediated CDKN1A promoter hypomethylation. Mol Ther. 2021;29(3):1279–93. https://doi.org/10.1016/j.ymthe.2020.11.018.
Huang J, Zhao J, Geng X, Chu W, Li S, Chen ZJ, et al. Long non-coding RNA lnc-CCNL1-3:1 promotes granulosa cell apoptosis and suppresses glucose uptake in women with polycystic ovary syndrome. Mol Ther Nucleic Acids. 2021;23:614–28. https://doi.org/10.1016/j.omtn.2020.12.008.
Lu X, Gao H, Zhu B, Lin G. Circular RNA circ_RANBP9 exacerbates polycystic ovary syndrome via microRNA-136-5p/XIAP axis. Bioengineered. 2021;12(1):6748–58. https://doi.org/10.1080/21655979.2021.1964157.
Xu L, Xiong F, Bai Y, Xiao J, Zhang Y, Chen J, et al. Circ_0043532 regulates miR-182/SGK3 axis to promote granulosa cell progression in polycystic ovary syndrome. Reprod Biol Endocrinol. 2021;19(1):167. https://doi.org/10.1186/s12958-021-00839-5.
Guo Y, Peng X, Cheng R, Chen H, Luo X. Long non-coding RNA-X-inactive specific transcript inhibits cell viability, and induces apoptosis through the microRNA-30c-5p/Bcl2-like protein 11 signaling axis in human granulosa-like tumor cells. Bioengineered. 2022;13(6):14107–17. https://doi.org/10.1080/21655979.2022.2080366.
Xia H, Zhao Y. miR-155 is high-expressed in polycystic ovarian syndrome and promotes cell proliferation and migration through targeting PDCD4 in KGN cells. Artif Cells Nanomed Biotechnol. 2020;48(1):197–205. https://doi.org/10.1080/21691401.2019.1699826.
Zhang Y, Wang Y, Feng X, Zhang S, Xu X, Li L, et al. Oocyte-derived microvilli control female fertility by optimizing ovarian follicle selection in mice. Nat Commun. 2021;12(1):2523. https://doi.org/10.1038/s41467-021-22829-2.
Webber LJ, Stubbs S, Stark J, Trew GH, Margara R, Hardy K, et al. Formation and early development of follicles in the polycystic ovary. Lancet. 2003;362(9389):1017–21. https://doi.org/10.1016/s0140-6736(03)14410-8.
Adhikari D, Liu K. Molecular mechanisms underlying the activation of mammalian primordial follicles. Endocr Rev. 2009;30(5):438–64. https://doi.org/10.1210/er.2008-0048.
Cinco R, Digman MA, Gratton E, Luderer U. Spatial characterization of bioenergetics and metabolism of primordial to preovulatory follicles in whole ex vivo murine ovary. Biol Reprod. 2016;95(6):129. https://doi.org/10.1095/biolreprod.116.142141.
Harris SE, Leese HJ, Gosden RG, Picton HM. Pyruvate and oxygen consumption throughout the growth and development of murine oocytes. Mol Reprod Dev. 2009;76(3):231–8. https://doi.org/10.1002/mrd.20945.
Han Q, Wang J, Li W, Chen ZJ, Du Y. Androgen-induced gut dysbiosis disrupts glucolipid metabolism and endocrinal functions in polycystic ovary syndrome. Microbiome. 2021;9(1):101. https://doi.org/10.1186/s40168-021-01046-5.
Tian Y, Zhang W, Zhao S, Sun Y, Bian Y, Chen T, et al. FADS1-FADS2 gene cluster confers risk to polycystic ovary syndrome. Sci Rep. 2016;6: 21195. https://doi.org/10.1038/srep21195.
Celik N, Aydin S, Ugur K, Yardim M, Acet M, Yavuzkir S, et al. Patatin-like phospholipase domain containing 3-gene (adiponutrin), preptin, kisspeptin and amylin regulates oocyte developmental capacity in PCOS. Cell Mol Biol (Noisy-le-grand). 2018;64(15):7–12.
Wang M, Zhao D, Xu L, Guo W, Nie L, Lei Y, et al. Role of PCSK9 in lipid metabolic disorders and ovarian dysfunction in polycystic ovary syndrome. Metabolism. 2019;94:47–58. https://doi.org/10.1016/j.metabol.2019.02.002.
Bizoń A, Franik G, Madej P. The role of proprotein convertase subtilisin/kexin type-9 concentration and paraoxonase 1 activities in the blood of women with polycystic ovary syndrome. Environ Toxicol Pharmacol. 2021;84:103612. https://doi.org/10.1016/j.etap.2021.103612.
Guo WJ, Wang YC, Ma YD, Cui ZH, Zhang LX, Nie L, et al. Contribution of high-fat diet-induced PCSK9 upregulation to a mouse model of PCOS is mediated partly by SREBP2. Reproduction. 2021;162(6):397–410. https://doi.org/10.1530/rep-21-0164.
Lewandowski KC, Komorowski J, O’Callaghan CJ, Tan BK, Chen J, Prelevic GM, et al. Increased circulating levels of matrix metalloproteinase-2 and – 9 in women with the polycystic ovary syndrome. J Clin Endocrinol Metab. 2006;91(3):1173–7. https://doi.org/10.1210/jc.2005-0648.
González F, Kirwan JP, Rote NS, Minium J. Glucose ingestion stimulates atherothrombotic inflammation in polycystic ovary syndrome. Am J Physiol Endocrinol Metab. 2013;304(4):E375-83. https://doi.org/10.1152/ajpendo.00491.2012.
Shalev E, Goldman S, Ben-Shlomo I. The balance between MMP-9 and MMP-2 and their tissue inhibitor (TIMP)-1 in luteinized granulosa cells: comparison between women with PCOS and normal ovulatory women. Mol Hum Reprod. 2001;7(4):325–31. https://doi.org/10.1093/molehr/7.4.325.
Chen Z, Wei H, Zhao X, Xin X, Peng L, Ning Y, et al. Metformin treatment alleviates polycystic ovary syndrome by decreasing the expression of MMP-2 and MMP-9 via H19/miR-29b-3p and AKT/mTOR/autophagy signaling pathways. J Cell Physiol. 2019;234(11):19964–76. https://doi.org/10.1002/jcp.28594.
Xiong F, Xiao J, Bai Y, Zhang Y, Li Q, Lishuang X. Metformin inhibits estradiol and progesterone-induced decidualization of endometrial stromal cells by regulating expression of progesterone receptor, cytokines and matrix metalloproteinases. Biomed Pharmacother. 2019;109:1578–85. https://doi.org/10.1016/j.biopha.2018.10.128.
Nikanfar S, Hamdi K, Haiaty S, Samadi N, Shahnazi V, Fattahi A, et al. Oncostatin M and its receptor in women with polycystic ovary syndrome and association with assisted reproductive technology outcomes. Reprod Biol. 2022;22(2): 100633. https://doi.org/10.1016/j.repbio.2022.100633.
Shen H, Xu X, Fu Z, Xu C, Wang Y. The interactions of CAP and LYN with the insulin signaling transducer CBL play an important role in polycystic ovary syndrome. Metabolism. 2022;131: 155164. https://doi.org/10.1016/j.metabol.2022.155164.
Popovic M, Sartorius G, Christ-Crain M. Chronic low-grade inflammation in polycystic ovary syndrome: is there a (patho)-physiological role for interleukin-1? Semin Immunopathol. 2019;41(4):447–59. https://doi.org/10.1007/s00281-019-00737-4.
Duleba AJ, Dokras A. Is PCOS an inflammatory process? Fertil Steril. 2012;97(1):7–12. https://doi.org/10.1016/j.fertnstert.2011.11.023.
Liu Y, Liu H, Li Z, Fan H, Yan X, Liu X, et al. The release of peripheral immune inflammatory cytokines promote an inflammatory cascade in PCOS patients via altering the follicular microenvironment. Front Immunol. 2021;12: 685724. https://doi.org/10.3389/fimmu.2021.685724.
Rudnicka E, Suchta K, Grymowicz M, Calik-Ksepka A, Smolarczyk K, Duszewska AM, et al. Chronic low grade inflammation in pathogenesis of PCOS. Int J Mol Sci. 2021;22(7): 3789. https://doi.org/10.3390/ijms22073789.
Jensterle M, Janez A, Fliers E, DeVries JH, Vrtacnik-Bokal E, Siegelaar SE. The role of glucagon-like peptide-1 in reproduction: from physiology to therapeutic perspective. Hum Reprod Update. 2019;25(4):504–17. https://doi.org/10.1093/humupd/dmz019.
Chugh RM, Park HS, El Andaloussi A, Elsharoud A, Esfandyari S, Ulin M, et al. Mesenchymal stem cell therapy ameliorates metabolic dysfunction and restores fertility in a PCOS mouse model through interleukin-10. Stem Cell Res Ther. 2021;12(1):388. https://doi.org/10.1186/s13287-021-02472-w.
Yang S, Li Q, Song Y, Tian B, Cheng Q, Qing H, et al. Serum complement C3 has a stronger association with insulin resistance than high-sensitivity C-reactive protein in women with polycystic ovary syndrome. Fertil Steril. 2011;95(5):1749–53. https://doi.org/10.1016/j.fertnstert.2011.01.136.
Wu Y, Zhang J, Wen Y, Wang H, Zhang M, Cianflone K. Increased acylation-stimulating protein, C-reactive protein, and lipid levels in young women with polycystic ovary syndrome. Fertil Steril. 2009;91(1):213–9. https://doi.org/10.1016/j.fertnstert.2007.11.031.
Daan NM, Koster MP, de Wilde MA, Dalmeijer GW, Evelein AM, Fauser BC, et al. Biomarker profiles in women with PCOS and PCOS offspring; a pilot study. PLoS One. 2016;11(11):e0165033. https://doi.org/10.1371/journal.pone.0165033.
Duffy DM. Novel contraceptive targets to inhibit ovulation: the prostaglandin E2 pathway. Hum Reprod Update. 2015;21(5):652–70. https://doi.org/10.1093/humupd/dmv026.
Ma Z, Zhao H, Zhang Y, Liu X, Hao C. Novel circular RNA expression in the cumulus cells of patients with polycystic ovary syndrome. Arch Gynecol Obstet. 2019;299(6):1715–25. https://doi.org/10.1007/s00404-019-05122-y.
Jia C, Wang S, Yin C, Liu L, Zhou L, Ma Y. Loss of hsa_circ_0118530 inhibits human granulosa-like tumor cell line KGN cell injury by sponging miR-136. Gene. 2020;744: 144591. https://doi.org/10.1016/j.gene.2020.144591.
Turathum B, Gao EM, Chian RC. The function of cumulus cells in oocyte growth and maturation and in subsequent ovulation and fertilization. Cells. 2021;10(9): 2292. https://doi.org/10.3390/cells10092292.
Aydos A, Gurel A, Oztemur Islakoglu Y, Noyan S, Gokce B, Ecemis T, et al. Identification of polycystic ovary syndrome (PCOS) specific genes in cumulus and mural granulosa cells. PLoS One. 2016;11(12): e0168875. https://doi.org/10.1371/journal.pone.0168875.
Wigglesworth K, Lee KB, Emori C, Sugiura K, Eppig JJ. Transcriptomic diversification of developing cumulus and mural granulosa cells in mouse ovarian follicles. Biol Reprod. 2015;92(1):23. https://doi.org/10.1095/biolreprod.114.121756.
Grøndahl ML, Andersen CY, Bogstad J, Borgbo T, Boujida VH, Borup R. Specific genes are selectively expressed between cumulus and granulosa cells from individual human pre-ovulatory follicles. Mol Hum Reprod. 2012;18(12):572–84. https://doi.org/10.1093/molehr/gas035.
Chaffin CL, Lee YS, VandeVoort CA, Patel BG, Latham KE. Rhesus monkey cumulus cells revert to a mural granulosa cell state after an ovulatory stimulus. Endocrinology. 2012;153(11):5535–45. https://doi.org/10.1210/en.2012-1733.
Natural Science Foundation of China (No. 81873819); Zhejiang Provincial Science and Technology project of China (No. 2021C03100).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The lncRNA-miRNA-mRNA Network.
The circRNA-miRNA-mRNA Network.
The Competing Endogenous RNA (ceRNA) Network.
The primers for real-time quantitative PCR.
The PNPLA3 mRNA levels with patients’ clinical characteristics. The correlation analysis of PNPLA3 mRNA levels with BMI (A), No. of 2PN fertilized (B), No. of day 3 good-quality embryos (C), 2PN fertilization rate (D), FSH levels (E), Estradiol (F) TT levels (G), Triglyceride (H), Total cholesterol (I), LDL-C (J), HOMA-IR (K).
The MVD mRNA levels with patients’ clinical characteristics. The correlation analysis of MVD mRNA levels with BMI (A), AFC (B), No. of oocytes retrieved (C), No. of MII oocytes (D), No. of 2PN fertilized (E), No. of day 3 good-quality embryos (F), 2PN fertilization rate (G), FSH levels (H), P levels (I), TT levels (J), AMH levels (K), HOMA-IR (L).
About this article
Cite this article
Li, Q., Sang, Y., Chen, Q. et al. Integrated bioinformatics analysis elucidates granulosa cell whole-transcriptome landscape of PCOS in China. J Ovarian Res 16, 154 (2023). https://doi.org/10.1186/s13048-023-01223-0