Polycystic Ovary Syndrome and Epithelial-mesenchymal Transition: Mendelian Randomization and Single-cell Analysis Insights

Overview

Journal J Ovarian Res

Publisher Biomed Central

Specialties Gynecology & Obstetrics
Reproductive Medicine

Date 2025 Feb 19

PMID 39972362

Authors

Dong Liu

Dan Liu

Kunyan Zhou

Affiliations

Soon will be listed here.

Abstract

Background: The process of epithelial-mesenchymal transition (EMT) may promote fibrosis in ovarian tissue related to polycystic ovary syndrome (PCOS), thus affecting ovarian function and hormonal balance.

Objective: This study aimed to explore key genes associated with EMT in PCOS and their potential molecular regulatory mechanisms, exclusively from the perspective of transcriptomics and single-cell RNA sequencing (scRNA-seq), combined with Mendelian Randomization (MR) analysis.

Methods: The dataset for PCOS and EMT-related genes (EMT-RGs) were sourced from public databases. The key genes in this study were identified via differential expression analysis, MR, and evaluation of expression levels. Enrichment analysis and a series of functional analyses were conducted on these genes to further elucidate their potential mechanisms. Subsequently, using scRNA-seq data and validation of the expression of key genes, key cell group in PCOS were identified, followed by pseudo-time and cell communication analyses to provide deeper insights.

Results: Three key genes, NUCB2 [odds ratio (OR) = 0.8634, 95% confidence interval (CI): 0.8145-0.9152, P < 0.0001], PGF (OR = 0.8393, 95% CI: 0.7185-0.9805, P < 0.05), and CRIM1 (OR = 0.7539, 95% CI: 0.6556-0.670, P < 0.0001), were identified as having a unidirectional causal association with PCOS and were associated with a reduced risk of PCOS. In public datasets, NUCB2 exhibited significantly increased expression in PCOS samples, while PGF and CRIM1 showed the opposite trends. These three genes were enriched in pathways related to cellular functions, metabolic processes, and the operation of the nervous system, and they were co-expressed in smooth muscle. Additionally, five cell clusters were annotated, among which fibroblasts were identified as key cells due to their highest expression of all three key genes. Further analysis revealed a bifurcation event occurring during the mid-development stage of fibroblasts, with PCOS samples displaying a higher abundance of fibroblasts. In PCOS samples, fibroblasts exhibited more extensive communication with secretory epithelial cells, indicating a more complex intercellular interaction within this condition.

Conclusion: This study identified three EMT-RGs: NUCB2, PGF, and CRIM1, which were associated with a reduced risk of PCOS, with fibroblast identified as a key cell group in the disease's pathology. This provides new insights for PCOS research.

References

Li M, Youngren J, Dunaif A, Goldfine I, Maddux B, Zhang B . Decreased insulin receptor (IR) autophosphorylation in fibroblasts from patients with PCOS: effects of serine kinase inhibitors and IR activators. J Clin Endocrinol Metab. 2002; 87(9):4088-93. DOI: 10.1210/jc.2002-020363. View

Escobar-Morreale H . Polycystic ovary syndrome: definition, aetiology, diagnosis and treatment. Nat Rev Endocrinol. 2018; 14(5):270-284. DOI: 10.1038/nrendo.2018.24. View

Zhang H, Meltzer P, Davis S . RCircos: an R package for Circos 2D track plots. BMC Bioinformatics. 2013; 14:244. PMC: 3765848. DOI: 10.1186/1471-2105-14-244. View

Yu G, Wang L, Han Y, He Q . clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16(5):284-7. PMC: 3339379. DOI: 10.1089/omi.2011.0118. View

Armanini D, Boscaro M, Bordin L, Sabbadin C . Controversies in the Pathogenesis, Diagnosis and Treatment of PCOS: Focus on Insulin Resistance, Inflammation, and Hyperandrogenism. Int J Mol Sci. 2022; 23(8). PMC: 9030414. DOI: 10.3390/ijms23084110. View

Zheng Y, Gao W, Zhang Q, Cheng X, Liu Y, Qi Z . Ferroptosis and Autophagy-Related Genes in the Pathogenesis of Ischemic Cardiomyopathy. Front Cardiovasc Med. 2022; 9:906753. PMC: 9279674. DOI: 10.3389/fcvm.2022.906753. View

Gao S, Long F, Jiang Z, Shi J, Ma D, Yang Y . The complex metabolic interactions of liver tissue and hepatic exosome in PCOS mice at young and middle age. Front Physiol. 2022; 13:990987. PMC: 9531160. DOI: 10.3389/fphys.2022.990987. View

Satija R, Farrell J, Gennert D, Schier A, Regev A . Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015; 33(5):495-502. PMC: 4430369. DOI: 10.1038/nbt.3192. View

Sadeghi H, Adeli I, Calina D, Docea A, Mousavi T, Daniali M . Polycystic Ovary Syndrome: A Comprehensive Review of Pathogenesis, Management, and Drug Repurposing. Int J Mol Sci. 2022; 23(2). PMC: 8775814. DOI: 10.3390/ijms23020583. View

10.

Hu J, Lin F, Yin Y, Shang Y, Xiao Z, Xu W . Adipocyte-derived exosomal miR-30c-5p promotes ovarian angiogenesis in polycystic ovary syndrome via the SOCS3/STAT3/VEGFA pathway. J Steroid Biochem Mol Biol. 2023; 230:106278. DOI: 10.1016/j.jsbmb.2023.106278. View

11.

Bowden J, Davey Smith G, Haycock P, Burgess S . Consistent Estimation in Mendelian Randomization with Some Invalid Instruments Using a Weighted Median Estimator. Genet Epidemiol. 2016; 40(4):304-14. PMC: 4849733. DOI: 10.1002/gepi.21965. View

12.

Du Y, Li F, Li S, Ding L, Liu M . Causal relationship between polycystic ovary syndrome and chronic kidney disease: A Mendelian randomization study. Front Endocrinol (Lausanne). 2023; 14:1120119. PMC: 10050750. DOI: 10.3389/fendo.2023.1120119. View

13.

Gustavsson E, Zhang D, Reynolds R, Garcia-Ruiz S, Ryten M . ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2. Bioinformatics. 2022; 38(15):3844-3846. PMC: 9344834. DOI: 10.1093/bioinformatics/btac409. View

14.

Kalluri R, Weinberg R . The basics of epithelial-mesenchymal transition. J Clin Invest. 2009; 119(6):1420-8. PMC: 2689101. DOI: 10.1172/JCI39104. View

15.

Miao C, Chen Y, Fang X, Zhao Y, Wang R, Zhang Q . Identification of the shared gene signatures and pathways between polycystic ovary syndrome and endometrial cancer: An omics data based combined approach. PLoS One. 2022; 17(7):e0271380. PMC: 9278750. DOI: 10.1371/journal.pone.0271380. View

16.

Love M, Huber W, Anders S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15(12):550. PMC: 4302049. DOI: 10.1186/s13059-014-0550-8. View

17.

Hartwig F, Davey Smith G, Bowden J . Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017; 46(6):1985-1998. PMC: 5837715. DOI: 10.1093/ije/dyx102. View

18.

Harris R, McAllister J, Strauss 3rd J . Single-Cell RNA-Seq Identifies Pathways and Genes Contributing to the Hyperandrogenemia Associated with Polycystic Ovary Syndrome. Int J Mol Sci. 2023; 24(13). PMC: 10341507. DOI: 10.3390/ijms241310611. View

19.

Szczuko M, Kikut J, Szczuko U, Szydlowska I, Nawrocka-Rutkowska J, Zietek M . Nutrition Strategy and Life Style in Polycystic Ovary Syndrome-Narrative Review. Nutrients. 2021; 13(7). PMC: 8308732. DOI: 10.3390/nu13072452. View

20.

Hemani G, Zheng J, Elsworth B, Wade K, Haberland V, Baird D . The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018; 7. PMC: 5976434. DOI: 10.7554/eLife.34408. View