- Institute of Reproductive Health, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
Polycystic ovarian syndrome (PCOS) is one of the most common gynecological endocrine disorders. MicroRNAs (miRNAs) play extensive roles in the pathogenesis of PCOS and can serve as potential diagnostic markers. However, most studies focused on the regulatory mechanisms of individual miRNAs, and the combined regulatory effects of multiple miRNAs remain unclear. The aim of this study was to identify the common targets of miR-223-3p, miR-122-5p, and miR-93-5p; and assess the transcript levels of some of these targets in PCOS rat ovaries. Transcriptome profiles of granulosa cells from PCOS patients were obtained from the Gene Expression Omnibus (GEO) database to identify differentially expressed genes (DEGs). A total of 1,144 DEGs were screened, 204 of which were upregulated and 940 were downregulated. According to the miRWalk algorithm, 4,284 genes were targeted by all three miRNAs at the same time, and intersection with DEGs was used to obtain candidate target genes. A total of 265 candidate target genes were screened, and the detected target genes were subjected to Gene ontology (GO) and KEGG pathway enrichment, followed by PPI network analysis. Then, qRT-PCR was used to determine the levels of 12 genes in PCOS rat ovaries. The expressions of 10 of these genes were found to be consistent with our bioinformatics results. In conclusion, JMJD1C, PLCG2, SMAD3, FOSL2, TGFB1, TRIB1, GAS7, TRIM25, NFYA, and CALCRL may participate in the development of PCOS. Our findings contribute to the identification of biomarkers that may promote the effective prevention and treatment of PCOS in the future.
Introduction
PCOS is a common endocrine and metabolic disease in women of reproductive age. According to the Rotterdam diagnostic criteria, it encompasses all combinations of ovulation dysfunction, hyperandrogenemia, and polycystic ovary morphology (Ehrmann, 2005; Rosenfield and Ehrmann, 2016). Patients with PCOS often have a range of other health problems, including infertility, insulin resistance (IR), and obesity (Patel, 2018). Several studies suggest that genetic, environmental, and epigenetic factors may play important roles in the pathogenesis of PCOS (Escobar-Morreale, 2018), but the exact mechanism remains largely unclear.
MiRNAs are small endogenous single-stranded non-coding RNA molecules consisting of 18–25 nucleotides that regulate gene expression at the post-transcriptional level (Fabian et al., 2010). It has been reported that miRNAs are involved in various signaling pathways in PCOS, including amino acid metabolism, hormone regulation, cell differentiation, etc. (Imbar and Eisenberg, 2014; Mu et al., 2021). Differentially expressed miRNAs play important roles in PCOS pathogenesis and serve as potential diagnostic markers (Deswal and Dang, 2020). In addition, miRNAs regulate many biological processes associated with obesity, including adipogenesis, insulin secretion, and glucose uptake (Butler et al., 2020). Inflammation of adipose tissue in obese patients contributes to obesity-related metabolic dysfunction, such as insulin resistance and type 2 diabetes (Wu and Ballantyne, 2020).
Upregulation of miR-223 and miR-93 was detected in adipose tissue of PCOS patients (Chen et al., 2013; Udesen et al., 2020), while miR-93 was also significantly upregulated in the ovarian cortex and follicular fluid of PCOS patients (Jiang et al., 2015; Butler et al., 2019). Interestingly, miR-223 and miR-93 were not only overexpressed in adipose tissue of PCOS patients, but also in control patients with IR, with which it was positively correlated in vivo (Chen et al., 2013; Chuang et al., 2015). MiR-223 and miR-93 downregulate GLUT4 expression and inhibit insulin-stimulated glucose uptake in adipocytes, suggesting that they may play an important role in other IR-related diseases such as T2DM and obesity. Increased serum miR-122 levels were found in PCOS patients with impaired glucose metabolism (Jiang et al., 2016). Elevated miR-122 in circulation was positively associated with obesity and IR in young adults (Wang et al., 2015). Furthermore, miR-122 and miR-223 were found to be increased in obesity or hyperglycemia, and their intracellular roles are related to the development of IR (Murri et al., 2018; Udesen et al., 2020).
Experimental studies confirmed that the expression of miR-223-3p, miR-122-5p and miR-93-5p was significantly upregulated in PCOS (Chen et al., 2013; Jiang et al., 2016). This may indicate that these differential miRNAs are key molecules involved in the pathological process of PCOS. However, most studies focused the regulation mechanism of individual miRNAs, while the combined regulatory effect of multiple miRNAs remained unclear.
Each miRNA may affect hundreds of targets, while itself being regulated by several distinct miRNAs. Since the efficacy of single markers is limited, and multi-marker-based models can provide more reliable information for the diagnosis and therapeutic management of PCOS, we aimed to identify the common targets of miR-223-3p, miR-122-5p, and miR-93-5p; and assess the transcript levels of some of these targets in PCOS rat ovaries.
Materials and methods
Microarray data and DEGs screening
The gene expression profile GSE34526 was acquired from the GEO website. GSE34526 contained 3 granulosa cells samples (GSM850527-GSM850529) from female controls and 7 granulosa cells samples (GSM850530-GSM850536) from PCOS patients, based on the GPL570 [HG-U133_Plus_2] Affymetrix Human Genome U133 Plus 2.0 Array. DEGs between PCOS and normal samples were obtained from the GEO database by GEO2R analysis (https://www.ncbi.nlm.nih.gov/geo/geo2r). A p-value <0.05 and |logFC| > 2 as used as the DEGs cut-off criteria. A volcano plot was generated to visualize DEG expression changes using GraphPad Prism 8.0 software (GraphPad, United States).
Prediction of target genes
The miRWalk Version 3.0 (https://mirwalk.umm.uni-heidelberg.de), which is linked to three online databases (TargetScan, miRDB, and miRTarBase), stores predicted data including experimentally verified miRNA-target interactions. The targets genes of miR-223-3p, miR-122-5p, and miR-93-5p were downloaded from the miRWalk 3.0 database, and the intersection genes were selected for further analysis. Then, the overlapping genes among DEGs and targets genes of the miRNAs were obtained using the Venn diagrams tool (https://bioinformatics.psb.ugent.be).
Functional enrichment analysis of target genes
GO functional enrichment and KEGG pathway analyses were performed for the gene overlaps using Metascape (https://metascape.org). Metascape always uses the latest data, which integrates data sources from GO, KEGG, UniProt, and DrugBank to achieve pathway enrichment and biological process annotation. p < 0.01 was considered statistically significant.
PPI network analysis of target genes
The Search Tool for the Retrieval of Interacting Genes (STRING) online database (https://string-db.org) was used to predict the protein-protein interaction (PPI) network specific to target genes. Then, PPI networks were visualized using Cytoscape version 3.8.2 (https://cytoscape.org) and the clusters (highly interconnected regions) in the PPI network were identified using the Cytoscape plugin MCODE.
Genes related to obesity with PCOS
The DEGs related to obesity in PCOS were identified in the dataset of GSE80432, including four granulosa cell samples (GSM2127203-GSM2127204, GSM2127215-GSM2127216) from normal weight PCOS patients and four granulosa cell samples (GSM2127209, GSM2127211-GSM2127212, GSM2127214) from obese PCOS patients. The datasets were analyzed on the GPL6244 [HuGene-1_0-st] Affymetrix Human Gene 1.0 ST Array [transcript (gene) version]. The p-value <0.05 and Fold Change >1.5 were set as the cut-off criteria for DEGs. The Venn diagram tool was used to identify the genes shared between DEGs and targets genes of miRNAs.
TFs related to miR-223-3p, miR-122-5p and miR-93-5p
The prediction of TFs related to miR-223-3p, miR-122-5p, and miR-93-5p was conducted using TransmiR v2.0 (https://www.cuilab.cn/transmir), which is a database for TF-miRNA regulatory association. Subsequently, the TFs involved in PCOS shared between DEGs and predicted TFs were screened using the Venn diagram tool.
Animal model
Female Sprague–Dawley rats (n = 20; 21 days old) were obtained from the Hubei Provincial Center for Disease Control and Prevention (Wuhan, China). Animal experiments were approved by the Animal Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (Wuhan, China). PCOS was induced by injection of DHEA (60 mg/kg body weight) dissolved in 0.2 mL sesame oil. The control group was injected with the vehicle (0.2 mL sesame oil). Daily treatment was continued for up to 21 days. In this study, we used a rat model that exhibits reproductive and metabolic abnormalities similar to human PCOS to uncover molecular mechanisms.
Quantitative real-time polymerase chain reaction (qRT-PCR)
Total RNA was extracted from the ovaries of PCOS rats using TRIzol reagent (Ambion). The cDNA as generated using the PrimeScript™ RT reagent kit or miRNA First-Strand Synthesis kit (Takara). The qRT-PCR was performed using qPCR SYBR Green Master Mix (Vazyme). GAPDH and U6 were used as internal controls for mRNA and miRNA expression, respectively. The relative expression of the mRNAs and miRNAs was calculated using the 2−(ΔΔCt) method.
Statistical analysis
Data are presented as the means ± SD. Data were analyzed using Student’s t-test. Differences with p-values <0.05 were considered statistically significant. The results were analyzed and visualized using GraphPad Prism 8.0 software (GraphPad Inc, United States).
Results
Identification of DEGs and target genes in PCOS
After normalization of the microarray data from GSE34526 (Figure 1A), a total of 1,144 DEGs were identified in PCOS patients compared with female controls, consisting of 204 upregulated genes and 940 downregulated genes (Figure 1B). The downstream target genes associated with miRNAs were identified, including 6,653 for miR-223-3p, 12,132 for miR-122-5p, and 15,170 for miR-93-5p. Among them, 4,284 target genes were shared among all three miRNAs. Finally, a total of 265 overlapping genes shared between DEGs and target genes were obtained, including 33 upregulated genes and 232 downregulated genes (Figure 1C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1097706/fgene-14-1097706-HTML/image_m/fgene-14-1097706-g001.jpg)
FIGURE 1. (A) Box plot of gene expression profiles after standardization. (B) Volcano plot of DEGs in PCOS. Red dots represent upregulation, blue dots represent downregulation, and gray dots represent no differential expression. (C) Venn diagram of DEGs and miRNA target genes.
GO and KEGG analysis of candidate target genes in PCOS
To further explore the function of the identified target genes, GO term and KEGG pathway enrichment analyses were performed. GO analysis was divided into the MF, CC, and BP categories (Figure 2). The results showed that 265 genes were mainly enriched in the regulation of cytokine production, actin filament-based processes, positive regulation of the immune response, cell junction organization, and myeloid leukocyte differentiation. Moreover, 12 KEGG pathways were overrepresented, including salmonella infection, morphine addiction, osteoclast differentiation, cell adhesion molecules, adherens junction, and endocytosis (Figure 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1097706/fgene-14-1097706-HTML/image_m/fgene-14-1097706-g002.jpg)
FIGURE 2. GO enrichment analysis of overlapping genes obtained from the intersection between DEGs and miRNA target genes. (A) Bar chart of the most highly enriched terms. (B) Network of the most highly enriched terms.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1097706/fgene-14-1097706-HTML/image_m/fgene-14-1097706-g003.jpg)
FIGURE 3. KEGG enrichment analysis of overlapping genes obtained from the intersection between DEGs and miRNA target genes. (A) Genes involved in enriched signaling pathways. (B) Network of the most highly enriched signaling pathways.
PPI network construction and hub gene selection
To explore the interactions among proteins encoded by the identified target genes, a PPI network was constructed, including 265 nodes and 389 edges (Figure 4). The three downregulated genes JMJD1C, PLCG2, and SERPINA1 were selected as hub genes based on enrichment degree ≥2 as the cutoff criterion.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1097706/fgene-14-1097706-HTML/image_m/fgene-14-1097706-g004.jpg)
FIGURE 4. Protein-protein interaction network analysis of the overlapping genes. The PPI network consists of 265 nodes and 389 edges. Red indicates upregulation, and blue indicates downregulation.
TFs associated with miR-223-3p, miR-122-5p and miR-93-5p
A total of 229 TFs associated with miR-223-3p, miR-122-5p, and miR-93-5p were downloaded from the TransmiR database. Compared with normal granulosa cells, 7 TFs were differentially downregulated in PCOS, including TGFB1, SMAD3, FOSL2, JMJD1C, CREB1, TRIM25, and NFYA.
DEGs related to obesity in PCOS
To further clarify the potential mechanisms of miRNAs related to adipogenesis in obese patients with PCOS, the genes that were differentially expressed in granulosa cells and associated with obesity were selected for further analysis. Based on the intersection between 265 genes and DEGs from the GSE80432 dataset, one upregulated and four downregulated genes were obtained, including CALCRL, TGFB1, TRIB1, GAS7, and FOSL2. TGFB1 and FOSL2 are both major transcription factors and closely related to the pathological process of obesity in PCOS.
Validation of the identified DEGs in a rat model of PCOS
To validate the results of bioinformatics analysis, we detected the expression of three miRNAs and 12 predicted target genes using qRT-PCR in PCOS model rats. As shown in Figure 5, miR-223-3p, miR-122-5p, and miR-93-5p were significantly upregulated in the PCOS group compared with control group. JMJD1C, PLCG2, SMAD3, FOSL2, TGFB1, TRIB1, GAS7, TRIM25, and NFYA were downregulated, while CALCRL was highly expressed in PCOS rats. Furthermore, SERPINA1 and CREB1 were not validated in PCOS rats.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1097706/fgene-14-1097706-HTML/image_m/fgene-14-1097706-g005.jpg)
FIGURE 5. The qPCR verification results of three miRNAs (A) and twelve candidate target genes (B) in the rat model of PCOS. An asterisk indicates a significant difference between the two groups (∗p < 0.05) according to Student’s t-test.
Discussion
PCOS is a common endocrine and metabolic disorder that can easily lead to ovulation dysfunction in women of childbearing age, resulting in infertility. The development of oocytes is inseparable from the role of follicular granulosa cells, which provide nutrients for oocytes through gap junctions and regulate oocyte development through paracrine signals (Clarke, 2018). In this study, a total of 265 candidate target genes were screened by integrating the information of DEGs and miRNA target genes in granulosa cells of PCOS patients with insulin resistance, among which 7 TF genes and five target genes closely related to obesity were identified. The function of candidate target genes was analyzed and the enriched GO terms included regulation of cytokine production, actin filament-based process, positive regulation of immune response, cell junction organization, and myeloid leukocyte differentiation.
We selected three DEGs as hub genes, the top three of which were JMJD1C, SERPINA1, and PLCG2. It stands to reason that miRNAs may play a major role in PCOS through these key target genes. We verified their transcriptional levels in a rat model of PCOS and found that the expression of JMJD1C and PLCG2 was downregulated, while the expression of SERPINA1 was unchanged. The mechanism may be that miRNA mediates the attenuation of JMJD1C and PLCG2 mRNA, hindering the translation of SERPINA1 mRNA without affecting its stability. Among them, JMJD1C and SERPINA1 play important roles in adipogenesis. JMJD1C regulates the induction of adipogenic transcription factors via H3K9me2 (Buerger et al., 2017) and promotes adipogenesis in vivo to increase liver and plasma triglyceride levels (Viscarra et al., 2020). SERPINA1, which may serve as an important molecular marker of obesity, affects energy expenditure by regulating the AMPK pathway and promotes the development of obesity-related metabolic complications (Mansuy-Aubert et al., 2013). PLCG2 regulates oxidative stress in hypoxic-ischemic encephalopathy, and activated PLCG2 signaling to attenuate oxidative stress-induced neuronal degeneration and apoptosis (Hu et al., 2020). Ovarian oxidative stress imbalance is a key feature of PCOS, suggesting that downregulation of PLCG2 is closely related to PCOS (Murri et al., 2013).
The selected transcription factors included TGFB1, SMAD3, FOSL2, JMJD1C, CREB1, TRIM25, and NFYA. The mRNA levels of TGFB1, SMAD3, FOSL2, TRIM25, and NFYA were decreased in our animal model. TGFB1 partially inhibits adipogenesis through SMAD3 (Tsurutani et al., 2011), and downregulation of the TGFB1 signaling pathway may promote adipogenesis by altering the expression of adipogenic genes to change the PCOS phenotype (Li and Wu, 2020). Furthermore, androgens were found to induce ovarian fibrosis through the TGFB1 signaling pathway in a rat model of PCOS (Wang et al., 2018). It is speculated that TGFB1 mainly regulates adipogenesis in PCOS and has less effect on ovarian fibrosis. FOSL2 is known to be a key regulator of adipokine LEP expression in obese mice and humans (Wrann et al., 2012). TRIM25 enhances the antioxidant defense by activating Nrf2 (Liu et al., 2020), and downregulation of TRIM25 may contribute to the oxidative stress imbalance phenotype of PCOS. NFYA regulates the gene expression of adiponectin in adipose tissue, one of the adipokines secreted by adipocytes that regulate energy homeostasis related to insulin sensitivity (Park et al., 2004).
Although CALCRL gene expression was reported to be lower in the adipose tissue of obese patients (Aguilera et al., 2015; Kim et al., 2020), CALCRL transcription was increased in the rat PCOS model. This difference may be caused by differences between species. TRIB1 knockout mice exhibit obesity and increased lipid accumulation in their livers (Bauer et al., 2015). TRIB1 can promote adipose tissue thermogenesis by regulating mitochondrial function (Zhang et al., 2021). PCOS patients have lower brown adipose tissue (BAT) activity and reduced BAT thermogenesis, which is associated with increased insulin resistance (Robinson et al., 1992; Shorakae et al., 2019). The expression of the GAS7 gene was found to be decreased in the ovaries of obese women, and its transcript was also downregulated in the ovaries of DHT-treated rats (Ruebel et al., 2017), which was consistent with our prediction. These results demonstrate the accuracy of predicting miRNA target genes and demonstrate that miRNAs regulate the pathological process of PCOS.
Although this study revealed the common targets of three miRNAs in PCOS using bioinformatics analyses and validated them using a rat model, there are still some limitations in the present study. Firstly, the GSE34526 and GSE80432 datasets were collected from granulosa cells samples, which may be the reason why the expression of SERPINA1 and CREB1 in PCOS rats was inconsistent with the bioinformatics analyses. In addition, the miRNA–mRNA relationships in PCOS were based on target prediction. However, the fact that the predicted targets are dysregulated in the PCOS rat model does not mean that they were targeted by the three miRNAs, which needs further experimental validation. Thus, a knockdown or overexpression of the miRNAs in granulosa cells and assessment of the mRNA and protein levels of the predicted targets would have yielded more reliable information.
In conclusion, JMJD1C, PLCG2, SMAD3, FOSL2, TGFB1, TRIB1, GAS7, TRIM25, NFYA, and CALCRL were found to be potential common targets of the three miRNAs. These possible targets were downregulated in PCOS rat ovaries, and so may serve as biomarkers for PCOS diagnosis and treatment. Instead of relying solely on the target gene prediction of a single miRNA, we used multiple miRNAs identified in clinical samples to predict common target genes, which provides strong evidence for the involvement of miRNAs in PCOS and broadens our understanding of gene expression changes in PCOS.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was reviewed and approved by Animal Ethics Committee of Tongji Medical College, Huazhong University of Science and Technology (Wuhan, China).
Author contributions
LZ conceived the study, analyzed the data, and wrote the manuscript. QF performed experiments. CZ and WX reviewed the paper.
Funding
This work was supported by the grants from the National Key Research and Development Program of China (2022YFC2702705).
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.
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.2023.1097706/full#supplementary-material
References
Aguilera, C. M., Gomez-Llorente, C., Tofe, I., Gil-Campos, M., Cañete, R., and Gil, Á. (2015). Genome-wide expression in visceral adipose tissue from obese prepubertal children. Int. J. Mol. Sci. 16 (4), 7723–7737. doi:10.3390/ijms16047723
Bauer, R. C., Sasaki, M., Cohen, D. M., Cui, J., Smith, M. A., Yenilmez, B. O., et al. (2015). Tribbles-1 regulates hepatic lipogenesis through posttranscriptional regulation of C/EBPα. J. Clin. Invest. 125 (10), 3809–3818. doi:10.1172/jci77095
Buerger, F., Müller, S., Ney, N., Weiner, J., Heiker, J. T., Kallendrusch, S., et al. (2017). Depletion of Jmjd1c impairs adipogenesis in murine 3T3-L1 cells. Biochim. Biophys. Acta Mol. Basis Dis. 1863 (7), 1709–1717. doi:10.1016/j.bbadis.2017.05.011
Butler, A. E., Ramachandran, V., Hayat, S., Dargham, S. R., Cunningham, T. K., Benurwar, M., et al. (2019). Expression of microRNA in follicular fluid in women with and without PCOS. Sci. Rep. 9 (1), 16306. doi:10.1038/s41598-019-52856-5
Butler, A. E., Ramachandran, V., Sathyapalan, T., David, R., Gooderham, N. J., Benurwar, M., et al. (2020). microRNA expression in women with and without polycystic ovarian syndrome matched for body mass index. Front. Endocrinol. (Lausanne) 11, 206. doi:10.3389/fendo.2020.00206
Chen, Y. H., Heneidi, S., Lee, J. M., Layman, L. C., Stepp, D. W., Gamboa, G. M., et al. (2013). miRNA-93 inhibits GLUT4 and is overexpressed in adipose tissue of polycystic ovary syndrome patients and women with insulin resistance. Diabetes 62 (7), 2278–2286. doi:10.2337/db12-0963
Chuang, T. Y., Wu, H. L., Chen, C. C., Gamboa, G. M., Layman, L. C., Diamond, M. P., et al. (2015). MicroRNA-223 expression is upregulated in insulin resistant human adipose tissue. J. Diabetes Res. 2015, 943659. doi:10.1155/2015/943659
Clarke, H. J. (2018). Regulation of germ cell development by intercellular signaling in the mammalian ovarian follicle. Wiley Interdiscip. Rev. Dev. Biol. 7 (1). doi:10.1002/wdev.294
Deswal, R., and Dang, A. S. (2020). Dissecting the role of micro-RNAs as a diagnostic marker for polycystic ovary syndrome: A systematic review and meta-analysis. Fertil. Steril. 113 (3), 661–669.e2. doi:10.1016/j.fertnstert.2019.11.001
Ehrmann, D. A. (2005). Polycystic ovary syndrome. N. Engl. J. Med. 352 (12), 1223–1236. doi:10.1056/NEJMra041536
Escobar-Morreale, H. F. (2018). Polycystic ovary syndrome: Definition, aetiology, diagnosis and treatment. Nat. Rev. Endocrinol. 14 (5), 270–284. doi:10.1038/nrendo.2018.24
Fabian, M. R., Sonenberg, N., and Filipowicz, W. (2010). Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 79, 351–379. doi:10.1146/annurev-biochem-060308-103103
Hu, X., Li, S., Doycheva, D. M., Huang, L., Lenahan, C., Liu, R., et al. (2020). Rh-CSF1 attenuates oxidative stress and neuronal apoptosis via the CSF1R/PLCG2/PKA/UCP2 signaling pathway in a rat model of neonatal HIE. Oxid. Med. Cell Longev. 2020, 6801587. doi:10.1155/2020/6801587
Imbar, T., and Eisenberg, I. (2014). Regulatory role of microRNAs in ovarian function. Fertil. Steril. 101 (6), 1524–1530. doi:10.1016/j.fertnstert.2014.04.024
Jiang, L., Huang, J., Chen, Y., Yang, Y., Li, R., Li, Y., et al. (2016). Identification of several circulating microRNAs from a genome-wide circulating microRNA expression profile as potential biomarkers for impaired glucose metabolism in polycystic ovarian syndrome. Endocrine 53 (1), 280–290. doi:10.1007/s12020-016-0878-9
Jiang, L., Huang, J., Li, L., Chen, Y., Chen, X., Zhao, X., et al. (2015). MicroRNA-93 promotes ovarian granulosa cells proliferation through targeting CDKN1A in polycystic ovarian syndrome. J. Clin. Endocrinol. Metab. 100 (5), E729–E738. doi:10.1210/jc.2014-3827
Kim, J., Lee, S. K., Kim, D., Choe, H., Jang, Y. J., Park, H. S., et al. (2020). Altered expression of adrenomedullin 2 and its receptor in the adipose tissue of obese patients. J. Clin. Endocrinol. Metab. 105 (1), dgz066. doi:10.1210/clinem/dgz066
Li, S. N., and Wu, J. F. (2020). TGF-β/SMAD signaling regulation of mesenchymal stem cells in adipocyte commitment. Stem Cell Res. Ther. 11 (1), 41. doi:10.1186/s13287-020-1552-y
Liu, Y., Tao, S., Liao, L., Li, Y., Li, H., Li, Z., et al. (2020). TRIM25 promotes the cell survival and growth of hepatocellular carcinoma through targeting Keap1-Nrf2 pathway. Nat. Commun. 11 (1), 348. doi:10.1038/s41467-019-14190-2
Mansuy-Aubert, V., Zhou, Q. L., Xie, X., Gong, Z., Huang, J. Y., Khan, A. R., et al. (2013). Imbalance between neutrophil elastase and its inhibitor α1-antitrypsin in obesity alters insulin sensitivity, inflammation, and energy expenditure. Cell Metab. 17 (4), 534–548. doi:10.1016/j.cmet.2013.03.005
Mu, L., Sun, X., Tu, M., and Zhang, D. (2021). Non-coding RNAs in polycystic ovary syndrome: A systematic review and meta-analysis. Reprod. Biol. Endocrinol. 19 (1), 10. doi:10.1186/s12958-020-00687-9
Murri, M., Insenser, M., Fernández-Durán, E., San-Millán, J. L., Luque-Ramírez, M., and Escobar-Morreale, H. F. (2018). Non-targeted profiling of circulating microRNAs in women with polycystic ovary syndrome (PCOS): Effects of obesity and sex hormones. Metabolism 86, 49–60. doi:10.1016/j.metabol.2018.01.011
Murri, M., Luque-Ramírez, M., Insenser, M., Ojeda-Ojeda, M., and Escobar-Morreale, H. F. (2013). Circulating markers of oxidative stress and polycystic ovary syndrome (PCOS): A systematic review and meta-analysis. Hum. Reprod. Update 19 (3), 268–288. doi:10.1093/humupd/dms059
Park, S. K., Oh, S. Y., Lee, M. Y., Yoon, S., Kim, K. S., and Kim, J. W. (2004). CCAAT/enhancer binding protein and nuclear factor-Y regulate adiponectin gene expression in adipose tissue. Diabetes 53 (11), 2757–2766. doi:10.2337/diabetes.53.11.2757
Patel, S. (2018). Polycystic ovary syndrome (PCOS), an inflammatory, systemic, lifestyle endocrinopathy. J. Steroid Biochem. Mol. Biol. 182, 27–36. doi:10.1016/j.jsbmb.2018.04.008
Robinson, S., Chan, S. P., Spacey, S., Anyaoku, V., Johnston, D. G., and Franks, S. (1992). Postprandial thermogenesis is reduced in polycystic ovary syndrome and is associated with increased insulin resistance. Clin. Endocrinol. (Oxf) 36 (6), 537–543. doi:10.1111/j.1365-2265.1992.tb02262.x
Rosenfield, R. L., and Ehrmann, D. A. (2016). The pathogenesis of polycystic ovary syndrome (PCOS): The hypothesis of PCOS as functional ovarian hyperandrogenism revisited. Endocr. Rev. 37 (5), 467–520. doi:10.1210/er.2015-1104
Ruebel, M. L., Cotter, M., Sims, C. R., Moutos, D. M., Badger, T. M., Cleves, M. A., et al. (2017). Obesity modulates inflammation and lipid metabolism oocyte gene expression: A single-cell transcriptome perspective. J. Clin. Endocrinol. Metab. 102 (6), 2029–2038. doi:10.1210/jc.2016-3524
Shorakae, S., Jona, E., de Courten, B., Lambert, G. W., Lambert, E. A., Phillips, S. E., et al. (2019). Brown adipose tissue thermogenesis in polycystic ovary syndrome. Clin. Endocrinol. (Oxf) 90 (3), 425–432. doi:10.1111/cen.13913
Tsurutani, Y., Fujimoto, M., Takemoto, M., Irisuna, H., Koshizaka, M., Onishi, S., et al. (2011). The roles of transforming growth factor-β and Smad3 signaling in adipocyte differentiation and obesity. Biochem. Biophys. Res. Commun. 407 (1), 68–73. doi:10.1016/j.bbrc.2011.02.106
Udesen, P. B., Glintborg, D., Sørensen, A. E., Svendsen, R., Nielsen, N. L. S., Wissing, M. L. M., et al. (2020). Metformin decreases miR-122, miR-223 and miR-29a in women with polycystic ovary syndrome. Endocr. Connect. 9 (11), 1075–1084. doi:10.1530/ec-20-0195
Viscarra, J. A., Wang, Y., Nguyen, H. P., Choi, Y. G., and Sul, H. S. (2020). Histone demethylase JMJD1C is phosphorylated by mTOR to activate de novo lipogenesis. Nat. Commun. 11 (1), 796. doi:10.1038/s41467-020-14617-1
Wang, D., Wang, W., Liang, Q., He, X., Xia, Y., Shen, S., et al. (2018). DHEA-induced ovarian hyperfibrosis is mediated by TGF-β signaling pathway. J. Ovarian Res. 11 (1), 6. doi:10.1186/s13048-017-0375-7
Wang, R., Hong, J., Cao, Y., Shi, J., Gu, W., Ning, G., et al. (2015). Elevated circulating microRNA-122 is associated with obesity and insulin resistance in young adults. Eur. J. Endocrinol. 172 (3), 291–300. doi:10.1530/eje-14-0867
Wrann, C. D., Eguchi, J., Bozec, A., Xu, Z., Mikkelsen, T., Gimble, J., et al. (2012). FOSL2 promotes leptin gene expression in human and mouse adipocytes. J. Clin. Invest. 122 (3), 1010–1021. doi:10.1172/jci58431
Wu, H., and Ballantyne, C. M. (2020). Metabolic inflammation and insulin resistance in obesity. Circ. Res. 126 (11), 1549–1564. doi:10.1161/circresaha.119.315896
Keywords: PCOS, miRNA, bioinformatics analysis, biomarkers, granulosa cells
Citation: Zou L, Feng Q, Xia W and Zhu C (2023) Bioinformatics analysis of the common targets of miR-223-3p, miR-122-5p, and miR-93-5p in polycystic ovarian syndrome. Front. Genet. 14:1097706. doi: 10.3389/fgene.2023.1097706
Received: 14 November 2022; Accepted: 06 February 2023;
Published: 16 February 2023.
Edited by:
Xiao Wang, Konge Larsen ApS, DenmarkReviewed by:
Enoch Appiah Adu Gyamfi, The State University of New York (SUNY), United StatesZhiming Zhao, Second Hospital of Hebei Medical University, China
Copyright © 2023 Zou, Feng, Xia and Zhu. 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: Changhong Zhu, reprodcentre@163.com