- Research
- Open access
- Published:
Identification and validation of senescence-related genes in polycystic ovary syndrome
Journal of Ovarian Research volume 17, Article number: 7 (2024)
Abstract
Background
Polycystic ovary syndrome (PCOS) is an exceedingly intractable issue affecting female endocrine and reproductive health. However, the etiology and intricate pathological mechanisms of PCOS remain unclear. Nowadays, aging was found to share multiple common pathological mechanisms with PCOS, which causes probing into the pathogenesis of PCOS from senescence. However, no bioinformatics analyses have specifically focused on connection between PCOS and ovarian aging.
Methods
Differentially expressed aging-related genes in PCOS were identified and then analyzed using function enrichment method. Hub genes were determined based on multiple algorithms, and expression validation of hub genes was performed in both datasets and experiments (human granulosa-like tumor cell line, KGN; human Granulosa Cell, hGCs). Finally, a transcription factor-miRNA-gene network of hub genes was constructed.
Results
Here, we identified 73 aging-related differential expression genes (ARDEGs) by intersecting DEGs in PCOS and senescence-related gene set. Furthermore, we performed biological functions and potential pathways of ARDEGs and potential hub genes were also screened by multiple algorithms. From the perspective of immune dysfunction, we analyzed the correlation between PCOS and immune cells. Finally, TF-miRNA-gene networks were constructed. Finally, TF-miRNA-gene networks were constructed.
Conclusions
Our work aimed to elucidate the relation between PCOS and cellular senescence based on bioinformatics strategy, deepening the understanding of mechanisms and to seek for novel therapy strategies for improving reproductive lifespan and female health. Exploring the potential molecular mechanism of cell aging in PCOS is expected to bring a new breakthrough for PCOS diagnosis and therapy strategies. And this, might deepen our understanding about intricate mechanisms of ovarian aging.
Background
Polycystic ovary syndrome (PCOS), a common female endocrine abnormality, is characterized by hyperandrogenism, ovulation dysfunction and ovarian polycystic morphology [1]. It’s not only the leading cause for female infertility but also contributes to a higher risk of obesity, hypertension, insulin resistance, type 2 diabetes, cardiovascular disease, and endometrial carcinoma [2]. However, the pathological mechanism of development and progression of polycystic ovary syndrome is still unclear.
Cellular senescence, also known as senescence-associated secretory phenotype (SASP), is a state of standing cell cycle arrest manifested by significant secretion of various bioactive molecules such as reactive oxygen species, pro-inflammatory cytokines, and chemokines [3]. Cellular senescence plays a central role in aging-related organ malfunction and it can be induced by aging in a variety of cells [3].
The primordial follicular pool, responsible for the women fertility and reproductive lifespan, has been established during the embryonic life and depleted with age in most of mammalians. Of note, granulosa cells (GCs) play an important part in maintaining the follicle pool [4]. Besides, GCs surround the oocyte, produce growth factors and sex steroids, and overall contribute to oocyte development [5]. With the process of GCs aging, meiosis progress is disturbed and led to oocyte dysfunction. And this, eventually causes the faliure in maturation and full developmental potential acquisition of oocyte [6]. In recent decades, with the progress of research on aging and PCOS, studies had revealed that PCOS and cellular senescence share some common pathogenesis, such as oxidative stress, inflammatory state, immune dysfunction, and so on [1,2,3, 7]. However, the pathophysiologies common to both between cellular senescence and PCOS is still unclear.
Gene expression analysis has become increasingly significant in biological research. Today, biomedicine is greatly facilitated by revolutionary developments in microarray technology and bioinformatics, which allow us to explore disease-related alterations in transcriptional and related biological pathways. Several published studies have used gene expression datasets obtained from the Gene Expression Omnibus (GEO) database to illuminate the biologic mechanisms underlying the pathogenesis of PCOS [2, 8, 9]. The results of those bioinformatic analyses provide ponderable hints to researchers for comprehending the pathogenesis of PCOS from different perspectives. However, no studies focus on the pathophysiology between aging and PCOS using bioinformatics analysis, and much potentially valuable information still waits to be revealed.
This study aims to analyze the PCOS-associated datasets from an angle of the aging-related genes (ARGs). The Aging-related differentially expression genes (DEGs, ARDEGs) were obtained by intersecting of DEGs and 1230 aging and senescence-associated genes. The potential biological functions and pathways of ARDEGs were analyzed by Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Set Enrichment Analysis (GSEA). Then, Weighted correlation network analysis (WGCNA) and Support vector machine-recursive feature elimination (SVM-RFE) algorithms were performed to identify potential hub genes involved in PCOS development. The expression level and diagnostic activity were validated using internal, external datasets, and receiver operating characteristic (ROC) curve. Additionally, the relevance between hub ARDEGs and immune infiltration within PCOS was explored. Finally, we conducted Transcription factor-microRNA (TF-miRNA) network of hub genes to explore the regulatory relationship between genes.
Results
Data processing and identification of DEGs
The flow diagram of the study is presented in Fig. 1. And the details of included datasets are listed in Table 1. We merged and normalized the expression matrix, then displayed it in principal component analysis (PCA) plots before and after the standardized processing (Fig. 2A,B). Then 446 up-regulated and 303 down-regulated genes were identified and displayed in volcano and heatmap plots (Fig. 2C,D).
Identification of ARDEGs and enrichment analysis
Totally, we obtained 1230 aging and senescence-associated genes after deleting duplicated genes; they are all listed in Supplementary Table S2. We intersected the DEGs and ARGs to screen out the ARDEGs. A total of 73 overlapped genes were revealed by venn plot (Fig. 3A), the genes were also visualized with a heatmap (Fig. 3B).
The BP results showed that the leukocyte cell–cell adhesion, leukocyte migration, leukocyte activation, cell activation involved in immune response, and cytokine-mediated signaling pathway were enriched (Fig. 3C). Besides, secretory granule membrane, protein complex involved in cell–matrix adhesion, external side of plasma membrane, specific granule, side, function including receptor ligand activity, signaling receptor activator activity, cytokine receptor binding, cytokine binding, and immune receptor activity were abundant in this category. The pathway analysis revealed that the cytokine-cytokine receptor interaction, Malaria, nuclear factor kappa-B (NF-κB) signal pathway, tumor necrosis factor (TNF) signal pathway, and Viral protein interaction with cytokine and cytokine receptor were the top 5 signal pathways enriched in KEGG terms (Fig. 3D). Moreover, we used GSEA to appraise the functional difference of ARDEGs between the PCOS and control samples. The GSEA result showed that ARDEGs were primarily involved in immune effector process, cell activation, leukocyte mediated immunity, cell activation involved in immune response, and myeloid leukocyte activation (Fig. 3E).
Co-expression network and hub aging-related DEGs
To identify aging-related hub genes of PCOS, two methods were applied. In the WGCNA algorithm, a soft-thresholding power of 8 was used to obtain approximate scale-free topology for the network (Fig. 4A). The genes with top 5000 variance were clustered into twelve modules (Fig. 4B), two modules (magent and purple) with a total of 183 genes were closely related to the disease. The clustering dendrogram showing the relation between eigengene and clinical characteristics can be seen in Fig. 4C. Correlation heatmap between genes was plotted based on topological overlap matrix (Fig. 4D).
In addition, the 73 ARDEGs were screened as crucial genes using the SVM-RFE algorithm. In total, 25 genes were identified as key candidate biomarkers including VEGFC, CXCR2, PECAM1, SMC5, DDIT3 and etc. (Fig. 4E,F). Besides, the genes overlapped in DEGs, genes from WGCNA key modules, and ARGs were obtained (Fig. 4G). The genes in WGCNA key modules and SVM-RFE algorithm were shown in Supplementary Table S3. Through the above two algorithms, we finally obtained the 4 key genes, namely ANXA3, CXCR2, IQGAP2 and LMNB1. The comprehensive information of the four molecules were shown in Table 2.
Predicted performance of hub biomarkers
To further identify the hub genes from the four candidate biomarkers, the ROC curve was performed for training dataset. The results showed that the all of the four biomarkers had an AUC of > 0.7 for predicting PCOS (Fig. 5A). Moreover, one extra dataset (GSE106724) was used to verify the above result as a testing group. The ROC analyses indicated that only CXCR2 and LMNB1 had an AUC of 0.78 and 0.84, respectively (Fig. 5B). We also validated the expression of four genes, differential expression of these genes, namely ANXA3,CXCR2, and LMNB1, was verified in the training group (Fig. 6A). Furthermore, only two genes (CXCR2 and LMNB1) showed the difference in gene expression in the testing dataset (Fig. 6B). As above-mentioned, ANXA3 and IQGAP2 showed unsatisfactory result for both diagnostic performance and expression difference. Therefore, we obtained two hub biomarkers, namely CXCR2 and LMNB1.
Validation of hub biomarkers
To verify the result of the bioinformatic analysis, the etoposide-induced aged granulosa cells (aGCs) model was established. The SA-β-gal assay validated the senescence phenotype (Fig. 7A). Using PCR analysis, we identified that CXCR2 was significantly upregulated and LMNB1 was significantly downregulated in aGCs when compared with young GCs (Fig. 7B). The baseline information for recruited women in this study was shown in Table 3.The above two genes were found up-regulated in PCOS patients when compared to the control group (Fig. 7C, D).
Immune cell infiltration analysis
The immune cells surrounding GCs were predicted. Compared to other immune cells, neutrophils and mast cells activated dominated in both PCOS and control groups (Fig. 8A). However, we detected no significant difference in immune cell infiltration among the groups (Fig. 8B). Moreover, we conducted correlation analysis to evaluate the correlation between hub genes and immune cells. The result revealed that the hub gene CXCR2 showed positive correlation with neutrophils (R = 0.630, P = 0.002), but negatively correlated with macrophage M1 cells (R = − 0.495, P = 0.019; Fig. 8C). Besides, we found that LMNB1 was negatively related with plasma cells (R = − 0.505, P = 0.016), macrophages M1 (R = − 0.498, P = 0.018), macrophages M2 (R = − 0.46, P = 0.031), and B cells naive (R = − 0.426, P = 0.048; Fig. 8D).
PPI network identification and construction of coregulatory network
To determine the regulatory mechanism, the proteins closely related to CXCR2 and LMNB1 were predicted using the STRING database. The proteins showed close relationship with CXCR2 including CXCL1, CXCL2, CXCL3, CXCL5, CXCL6, CXCL8, PPBP, ARB1,CD74, and VASP (Fig. 9A). For LMNB1 protein network, LEMD3, LMNA, TMPO, EMD, LBR, BANF1, SUN1, SUN2, PPKCA, and CDK1 were predicted (Fig. 9B). The interaction between TFs and hub genes were analyzed, the result showed four TF-genes (CEBP1, NFKB1, RELA, USF2) regulated the expression of CXCR2. Meanwhile, nine TF-genes play roles in TF-LMNB1 regulatory network. In order to comprehensively understand the connection that contributes to controll the hub genes expression, microRNA-gene network analyses were also performed. The TF-miRNA coregulatory network was showed by Cytoscape tool (Fig. 9C,D). 23 and 157 miRNAs were predicted to be involved in the expression regulation of hub genes. The information of edges was listed in supplementary Table S4.
Discussion
Polycystic ovary syndrome is a complex and intractable endocrinopathy, which exerts adverse influence on both reproduction and metabolism for women of childbearing age [10]. Cellular senescence is a condition of standing cell cycle arrest, which contribute great importance to aging and aging-associated diseases. Studies revealed that PCOS and ovarian aging share multiple common mechanisms in terms of pathogenesis such as oxidative stress, cellular metabolic process, mTOR pathway, and so on [11, 12]. Nonetheless, the underlying mechanism that promotes cell senescence and PCOS pathophysiology remain largely unknown. In this study, we aimed to reveal the intrinsic connection between PCOS and cell senescence, which may reveal the potential role of PCOS in female ovarian aging.
Recently, dozens of studies have been published to highlight the role of ARGs in various aging-related diseases. For instance, research by Xiao et al. probed into the potential role of aging-related genes in glioma [13]. Also, an aging-related signature genes for the prognosis of lung cancer was established [14]. Except for carcinomas, the role of ARGs in non-tumor conditions gained wide concern, such as asthma, pulmonary fibrosis, and so on [15, 16]. However, bioinformatic analysis focusing on ARGs in PCOS remains to be performed, and signature ARGs involved in PCOS pathogenesis are not determined.
In this study, we found the ARDEGs by intersecting DEGs and 1230 aging and senescence-associated genes. In addition, four hub ARDEGs (ANXA3, CXCR2, IQGAP2, and LMNB1) were identified by WGCNA and SVM-RFE algorithm, and the expression and diagnostic value was also validated using internal and external datasets.
The enrichment analysis of biological function for ARDEGs was performed. For GO analysis, it mainly focuses on immune-associated function, such as leukocyte cell migration, cytokine − mediated signaling pathway, and so. Similarly, GSEA method showed that those genes were mainly enriched in immune effector process, leukocyte mediated immunity, and cell adhesion. The function of the immune system declines during aging, actually, aging was also deemed as a chronic low-grade inflammatory state [17]. In the future, the molecular mechanism of immune in aging process needs further exploration.
The KEGG pathway investigation for ARDEGs indicated that those ARDEGs were primarily involved in NF-κB, PI3K-Akt, and JAK/STAT signal pathways. NF-κB pathway was revealed to be involved in the process of skin aging [18]. In osteoarthritis, a common aging-related degenerative joint disease, NF-κB signaling was also activated [19]. Previous research had linked the PI3K-Akt signal pathway to age-related diseases. Specially, it was found that the PI3K-Akt-mTOR signal pathway is active during hippocampus aging [20]. Furthermore, study showed that the drugs adjusting PI3K-Akt signal pathway played a role in anti-aging [21]. Above all, the result of functional enrichment showed that the ARDEGs in PCOS had close relationship within aging process, the role of those ARDEGs in ovarian aging deserves our attention and further studying.
To find the most important hub genes among the 73 ARDEGs, we performed WGCNA and SVM-RFE analysis to find 4 signature genes. Then the expression and diagnostic activity were assessed using both training and testing datasets. Besides, relative expression quantification was also performed in senescent GCs and GCs from PCOS patients by qRT-PCR. CXCR2 and LMNB1 were confirmed the excellent importance in PCOS and aging. Notably, CXCR2 showed significant up-regulation in both aged GCs and GCs of PCOS group, which indicated its potential role in promoting GCs senescence. However, expression trend of LMNB1 in PCOS GCs (up-regulated) is contrary to that in aGCs (down-regulated). This hinted us that PCOS might be a protective factor for ovarian aging and the underlying mechanism requires further exploration in the future.
We reviewed the role of CXCR2 and LMNB1 in PCOS and aging. Firstly, we found that CXCR2 showed great importance in senescence and tumor. Senescent cancer cells promoted CXCR2 expression and mediated the resistance against therapy, meanwhile, inhibition of CXCR2 was found to promote cancer cell senescence process, then improved the survival period of patients with Head-and-Neck and lung carcinoma [22, 23]. Also, the studies indicated that the gene was important in other cell aging-associated conditions, like depression, Alzheimer’s disease (AD), pathological healing, and atrophic Gastritis [24,25,26,27]. Recently, researchers found that inhibition of CXCL5-CXCR2 signal pathway in aged mouse embryos could improve the implantation rate, increasing the pregnancy rates and the delivery of normal pups [28]. In terms of mechanism of CXCR2, senescence-induced alterations in genes expression of CXCR2 may interacted with C-X-C motif chemokine ligand (CXCL) and promoted the proliferation and migration of pancreatic stellate cells [29]. Blockade of CXCR2 through toll-like receptor (TLR)2, and TLR5 signaling inhibited the secretion of interukin-8 (IL-8) or growth-regulated oncogene-alpha (GROa) by senescent mesenchymal stem cells [30]. Additionally, CXCR2 and p53 formed a positive feedback loop to constantly enhance cellular senescence in helicobacter pylori-induced atrophic gastritis [25]. CXCR2 knockout monocytes infused mice differentiated in TNF-α-releasing macrophages, leading to cell senescence [31]. IL-8/CXCR2 signaling pathway were found to accelerate senescence of placenta-derived mesenchymal stem cells [32].
Functional effect of nuclear LMNB1 was also studied, according to the cellular localization and function, the gene are crucially involved in nuclear processes like chromatin organization. This nuclear lamina degradation upon aberrant cellular stress induced cellular senescence, also, the loss of methylated LMNB1 played roles in aged human cells [33]. Besides, LMNB1 knockdown experiment revealed the implication in important cellular processes driving ER stress leading to cell senescence in malignant melanoma [34].
Besides, pathway analysis was performed for CXCR2 and LMNB1, and the result indicated an intrisic connection between CXCR Chemokine Receptors and Aging process (Table S5). Specifically, Aging is associated with changes in the immune system, a process known as immunosenescence. These changes can affect the function and expression of various immune receptors, including chemokine receptors like CXCR2. Also, these genes regulate the chemotaxis of neutrophils and leukocytes to some degree as identified by the PPI and pathway analysis above. The expression and function of chemokine receptors may be altered with age, and this can influence immune cell trafficking and responses to inflammation. Inflammation is directly associated with PCOS, which might help you tie things together.
Since the ARDEGs analysis of GO and GSEA were mostly enriched in immune pathway, which indicated the immune cells and immune dysfunction have relation to PCOS and cellular senescence. In this study, no different immune infiltration of immunocytes were found between groups, notably, correlation analysis revealed a significant connection between hub genes and immunocytes. Spearman analysis revealed that CXCR2 and LMNB1 were negatively correlated with macrophage M1, and CXCR2 was positively related to neutrophils. Furthermore, LMNB1 was negatively correlated with macrophage M2, B naïve, and plasma cells. PCOS is associated with increased neutrophils, which corelated positively with two hub biomarkers. Besides, PCOS is characterized by a shift from an anti-inflammatory M2 to a proinflammatory M1 state, yet, our finding of opposite connection of CXCR2 and LMNB1 with macrophage M1 indicated that the signature genes and the cell in PCOS need to be further clarified [35]. Besides, the chronic inflammation caused by imbalanced immunocytes may contribute to tissue fibrosis, one of the main mechanistic pillars of aging, and ultimately results in ovarian aging. As far as we know, CXCR2 and LMNB1 have not been studied in PCOS before; however, we speculated their latent roles in the immune process of PCOS. In the future, further research is required to reveal the intricate relationship between the hub genes and immune infiltration in both PCOS and aging.
To find the regulatory molecules among the 2 signature genes, we constructed a TF-miRNA-mRNA coregulatory network. miRNAs take part in cell senescence pathophysiology and have influence on the disease progression. According to the analysis result, the two hub genes were both regulated by hsa-mir-194-5p, hsa-mir-27a-3p, and hsa-mir-429. By examining the interaction between hub genes, TFs and miRNAs, we discovered that several miRNAs are implicated in PCOS (such as hsa-mir-29a-5p, hsa-mir-320, hsa-mir-429, hsa-mir-103a-3p, hsa-mir-192-5p, hsa-mir-145-5p, hsa-mir-218-5p, hsa-mir-1224-5p, hsa-mir-27a-3p) [36,37,38,39,40,41,42] and aging (such as hsa-mir-194-5p, hsa-mir-26b-5p, has-let-7e-5p, hsa-mir-181a-5p, hsa-mir-30a-3p, hsa-mir-200b-3p, hsa-mir-200c-3p, hsa-mir-30a-5p, hsa-mir-130a-3p, hsa-mir-506-3p, hsa-mir-93-5p, hsa-mir-148b-3p, hsa-mir-15b-5p, hsa-mir-17-3p) [43,44,45,46,47,48,49,50,51,52,53]. Some miRNAs have dual roles in cellular senescence and PCOS including hsa-mir-22-3p, hsa-mir-335-5p, 126-3p, 21-5p, 23a-3p, 124-3p, 122-5p, 34a-5p, 125a-5p, 125b-5p, 155-5p, 760, 199a-5p, and 424-5p [38, 45, 46, 54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72].
Several limitations of the study require to be noted. Firstly, the analysis is mainly based on the public datasets, which contained relatively small sample size. Thus, the results may be biased, and subsequent research with expand samples is required for obtaining more reliable conclusions. Secondly, screening of hub gene is mainly based on WGCNA and computational algorithm, and the results are hypothetical, in vitro and in vivo experiments are necessary to confirm the hub genes in PCOS and aging. Therefore, we can build on the existing research to provide a theoretical framework and foundation for further exploring the relation between PCOS and cellular senescence through gain or loss of function on biological models.
Conclusions
Altogether, our research determined two potential senescence-associated hub genes (CXCR2, and LMNB1) with good diagnostic properties in PCOS using bioinformatic analysis and machine learning methods. Furthermore, we constructed a TF-miRNA coregulatory pathway for the genes. The characteristic genes were remarkably linked with immunological functions, besides, CXCR2 and LMNB1 per se and their regulatory miRNAs showed close connection with multiple cellular senescence pathways. These genes may affect the occurrence and prognosis of PCOS by regulating senescence. The above findings deepen our understanding of the unique relationship between cellular senescence and PCOS, and indicate that in-depth study of PCOS may provide insight into the process of granulosa cell senescence. Based on the results from bioinformatic analysis, we can provide a theoretical framework and foundation for further exploring the relation between PCOS and cellular senescence through gain or loss of function on biological models. Further research is needed, which may lay the theoretical foundation for future studies in ovarian aging and improve reproductive lifespan for women.
Methods
Data collection and processing
Microarray data, including four transcription profiles [GSE34526, GSE102293, GSE114419, GSE106724] were downloaded from the NCBI GEO database. Inclusion criteria was as follows: (1) 2003 Rotterdam criteria was as the basis for PCOS diagnosis; (2) mural GCs were isolated from follicular fluid aspirates; (3) unaffected women were selected as a control group (notably, the control part included women of childbearing age with an infertile male partner, tubal obstruction, or unexplained infertility seeking for IVF); (4) quantitative real-time PCR (qRT-PCR) was performed for experimental validation. GSE34526, GSE102293, GSE114419 datasets were used as the training dataset and GSE106724 was used to validate the hub genes. The ComBat method was used to normalize the data. Raw data was processed by R statistical software (version 4.2.1) and “Bioconductor” tools. Genes without being simultaneously included in there datasets were excluded. The “sva” R package was applied to remove the batch effect.
Aging and senescence-associated genes
Aging and senescence-associated genes were searched from the Human Aging Genomic Resources2 (https://genomics.senescence.info/genes/index.html) and MSigDB gene sets3 [including datasets M16568, M24705, M14683, M27188, M13158, M14701, M11850, M40058, M39619] (https://www.gsea-msigdb.org/gsea/msigdb). Besides, a new senescence-related gene set [SenMayo] discovered by Mayo clinic was also included in this study [73].
Differential expression genes analysis
The “limma” package was used for analyzing differentially expressed genes in granulosa cells between control group and women with PCOS. And P-value < 0.05 with |Fold Change|> 1.5 were as the threshold for selecting DEGs. The significantly up-regulated and down-regulated genes were visualized through Heatmap and volcano plot using “Pheatmap” and “ggplot2” packages, respectively.
Identification and function enrichment analysis of ARDEGs
The Venn online tool (https://bioinfogp.cnb.csic.es/tools/venny/index.html) was applied for screening intersections between differential expression genes and cellular senescence-associated genes, namely, ARDEGs, which was visualized with heatmap plot using the “Pheatmap” package.
We used the “ClusterProfiler” R package to enrich the biological function of ARDEGs. GO analysis identified the biological phenomena of DEGs, which primarily consists of biological processes (BP), cellular components (CC), and molecular functions (MF). Moreover, pathway enrichment was done using KEGG. GSEA on cellular and senescence-associated DEGs was also performed using the “c5.all.v7.0.entrez.gmt” gene set, which was downloaded from the MSigDB database.
Construction of the co-expression network and hub ARDEGs identification
WGCNA is an algorithm to cluster genes into different modules and uncover the relationships between modules and disease traits [74]. To identify the key modules involved in the pathogenesis of PCOS, we used “WGCNA” package to perform the WGCNA. The analysis included following four steps: 1) Genes with top 5000 variance were selected to constructed the co-expressed network and gene module; the soft threshold power was set to the scale-free topology criterion (signed R2 = 0.85). 2) To assess network connectivity, acquired weighted adjacency matrix was transformed into a topological overlap matrix. 3) The matrix was then used to generate a clustering tree using the average-linkage hierarchal clustering method. 4) The genes were classified according to their expression patterns; in this way, the genes were separated into different modules.
SVM-RFE algorithm was also performed to screen out the aging-related crucial genes in PCOS patients. SVM is a machine learning method and has been widely used for classification and regression [75]. The SVM-RFE requires multiple classifiers to obtain the total ranking criterion score, which is used as the basis for feature rejection, so as to avoid over-fitting and search for the best features.
In this study, genes overlapped in DEGs, genes from WGCNA key modules, and ARGs were intersected with important genes identified by SVM-RFE algorithm to identify hub ARDEGs.
Expression validation and predicted performance of feature genes
To verify the above bioinformatics analyses, we draw the ROC curve to further evaluate the diagnostic capability of candidate biomarkers using “pROC” package in R software. P < 0.05 indicated the difference is statistically significant. Furthermore, internal and external datasets were also used for validating the expression of potential hub genes.
Review the role of hub genes in PCOS and aging
The basic and functional information for aging-related PCOS signature genes was achieved obtained from the BioGPS website (http:// biogps.org), the Human Protein Atlas (https://www.proteinatlas.org/), and GeneCards (https://www.genecards.org/). Besides, the Pubmed database was systematically searched with the keyword “CXCR2”, “LMNB1”, “PCOS”, “Polycystic Ovary Syndrome”, “aging”, and “senescence”. Articles that studied the relationship between hub genes and PCOS/aging were reviewed.
Immune infiltration analysis and its correlation with senescence-associated hub genes
Based on the microarray data, CIBERSORT research was performed to predict the differential abundance of 22 leukocyte subtypes between PCOS and control groups [76]. Bar plot of fractions and boxplot were draw by using the “ggplot2” package, respectively. Further, we used Spearman’s method for the correlation analysis between senescence-associated signature genes expression and the content of immune cells, p < 0.05 indicated significant difference.
Regulatory network analysis of CXCR2 and LMNB1
microRNA and TF usually take part in the gene transcription and post-transcription regulation. To explore the regulatory network and further identify the intrinsic associations of hub genes with PCOS and aging, the protein–protein interaction (PPI) network, TF, and microRNA that interacted with CXCR2 and LMNB1 were analyzed. The STRING online database (https://string-db.org/) was applicated to study hub gene and protein interactions. The Cistrome DB database (http://cistrome.org/db/#/) was searched for finding TFs of hub genes. MicroRNAs those targeting genes were also screened from miRNANet database (https://www.mirnet.ca/miRNet/home.xhtml). Then, we incorporated the TF-miRNA coregulated interactions. Cytoscape software (version 3.7.2) was used to visualize the PPI network and interactions between TFs, microRNAs, and signature genes.
Cell line culture and etoposide treatment
KGN, the human granulosa-like tumor cell line, was purchased from the Shanghai Fuheng Biotechnology. KGN cell line was cultured in DMEM/F12 medium containing 10% fetal bovine serum, 100 U/mL of penicillin and 100 μg/mL of streptomycin. Those cells were grown in 95% humidified air and 5% CO2 at 37°C. Cells were exposed to 5 μM etoposide for 48 h (in complete medium) to trigger cellular senescence. And cells were allowed a recovery period of 24 h prior to the following experiments.
Senescence-associated β-galactosidase staining (SA-β-gal assay)
SA-β-gal assay was applied to identify KGN senescence. KGN cells were seeded on a 24-well plate at a density of 2 ∗ 104 cells/well. After fixation, cells were stained with SA-β-gal staining solution overnight at 37°C and in CO2-free environment (C0602, Beyotime Biotechnology). Cells were imaged and photographed under a light microscope, and then we manually counted 3 fields of SA-β-gal-positive cells to determine the percentage.
Patients and sample acquisition
Between September 2021 and October 2022, we recruited 30 women who underwent in vitro fertilization and embryo transfer (IVF-ET) at the Reproductive Center of Qilu Hospital, Shandong University. 15 patients were diagnosed as PCOS based on the Rotterdam criteria, and 15 control patients were infertile due to either fallopian tube obstruction or husband’s subfertility.
After the ovarian stimulation, the follicular fluid was extracted under ultrasound echo guidance and immediately centrifuged for 10 min at 1000 g. We suspended the precipitates with lymphocyte separation medium (P8610, Solarbio), and extracted the intermediate layer cells (human Granulosa Cell, hGCs) after the centrifugation. hGCs were then stored at − 80 °C until RNA extraction.
Quantitative real‑time PCR (RT‑qPCR)
RNA extraction was performed from KGN cells and hGCs using RNA-easy Isolation Reagent (R701, Vazyme), and then total RNA was reverse-transcribed into cDNA with an RT Kit (11142ES60, Yesen). Each PCR system (10 μL) contained 5μL SYBR Green (11184ES03, Yesen), 30 ng cDNA and 300 nmol/L specific primers. We listed the primer sequence of each gene in Supplementary Table S1. And the PCR program consisted of 2 min at 95°C, followed by 40 cycles of 95°C for 10 s and 60°C for 30 s. We calculated the gene expression using the 2-ΔΔCT method and β-Actin was as the internal control.
Statistical analysis
Expression data for granulosa cells obtained by RT-qPCR was analysed with t-test using the GraphPad Prism version 8.0.2 software. Differences were considered statistically significant when the P-value was < 0.05.
Availability of data and materials
Microarray data in this work are available in the GEO online database (http://www.ncbi.nlm.nih.gov/geo).
Abbreviations
- PCOS:
-
Polycystic ovary syndrome
- DEGs:
-
Differentially expressed genes
- ARGs:
-
Aging-related genes
- ARDEGs:
-
Aging-related differential expression genes
- WGCNA:
-
Weighted gene co-expression network analysis
- SVM-RFE:
-
Support vector machine recursive feature elimination
- AUC:
-
Area under the curve
- GEO:
-
Gene Expression Omnibus
- GO:
-
Gene ontology
- BP:
-
Biological processes
- CC:
-
Cellular components
- MF:
-
Molecular functions
- KEGG:
-
Kyoto encyclopedia of genes and genomes
- GSEA:
-
Gene set enrichment analysis
- PPI:
-
Protein–protein interaction
- TF:
-
Transcriptional factor
- miRNA:
-
MicroRNA
- GCs:
-
Granulosa cells
- aGCs:
-
Aged granulosa cells
- yGCs:
-
Young granulosa cells
- SA-β-gal assay:
-
Senescence-associated β-galactosidase staining
- ETO:
-
Etoposide
References
Xia Q, Wang W, Liu Z, Xiao J, Qiao C, Zhao Y, et al. New insights into mechanisms of berberine in alleviating reproductive disorders of polycystic ovary syndrome: anti-inflammatory properties. Eur J Pharmacol. 2023;939:175433.
Yao X, Wang X. Bioinformatics searching of diagnostic markers and immune infiltration in polycystic ovary syndrome. Front Genet. 2022;13:937309.
Roy AL, Sierra F, Howcroft K, Singer DS, Sharpless N, Hodes RJ, Wilder EL, et al. A blueprint for characterizing senescence. Cell. 2020;183(5):1143–6.
Liu C, Peng J, Matzuk MM, Yao HHC. Lineage specification of ovarian theca cells requires multicellular interactions via oocyte and granulosa cells. Nat Commun. 2015;6:6934.
Jozkowiak M, Piotrowska-Kempisty H, Kobylarek D, Gorska N, Mozdziak P, Kempisty B, et al. Endocrine disrupting chemicals in polycystic ovary syndrome: the relevant role of the theca and granulosa cells in the pathogenesis of the ovarian dysfunction. Cells. 2022;12(1):174.
Rb G, Lj R, Dt A. Oocyte-Somatic Cell Interactions during Follicle Development in Mammals. Animal reproduction science. 2004; 82–83.
de Magalhães JP, Passos JF. Stress, cell senescence and organismal ageing. Mech Ageing Dev. 2018;170:2–9.
Wu T, Gao YY, Tang XN, Li Y, Dai J, Zhou S, et al. Construction of a competing endogenous RNA network to identify drug targets against polycystic ovary syndrome. Hum Reprod. 2022;37(12):2856–66.
Cozzolino M, Herraiz S, Titus S, Roberts L, Romeu M, Peinado I, et al. 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. 2022;37(6):1274–86.
Abolhassani N, Winterfeld U, Kaplan YC, Jaques C, Minder WB, Del GC, et al. Major malformations risk following early pregnancy exposure to metformin: a systematic review and meta-analysis. BMJ Open Diabetes Res Care. 2023;11(1):e002919.
Ho CH, Chang CM, Li HY, Shen HY, Lieu FK, Wang PSG. Dysregulated immunological and metabolic functions discovered by a polygenic integrative analysis for PCOS. Reprod Biomed Online. 2020;40(1):160–7.
Dong L, Teh DBL, Kennedy BK, Huang Z. Unraveling female reproductive senescence to enhance healthy longevity. Cell Res. 2023;33(1):11–29.
Xiao G, Zhang X, Zhang X, Chen Y, Xia Z, Cao H, et al. Aging-related genes are potential prognostic biomarkers for patients with gliomas. Aging (Albany NY). 2021;13(9):13239–63.
Xu Q, Chen Y. An aging-related gene signature-based model for risk stratification and prognosis prediction in lung adenocarcinoma. Front Cell Dev Biol. 2021;9:685379.
Yang Y, Yuan L, Yang M, Du X, Qin L, Wang L, et al. Aberrant methylation of aging-related genes in asthma. Front Mol Biosci. 2021;8:655285.
He J, Li X. Identification and validation of aging-related genes in idiopathic pulmonary fibrosis. Front Genet. 2022;13:780010.
Cisneros B, García-Aguirre I, Unzueta J, Arrieta-Cruz I, González-Morales O, Domínguez-Larrieta JM, et al. Immune system modulation in aging: molecular mechanisms and therapeutic targets. Front Immunol. 2022;13:1059173.
Trentini M, Zanolla I, Zanotti F, Tiengo E, Licastro D, Dal Monego S, et al. Apple derived exosomes improve collagen type i production and decrease MMPs during aging of the skin through downregulation of the NF-ΚB pathway as mode of action. Cells. 2022;11(24):3950.
Lu R, He Z, Zhang W, Wang Y, Cheng P, Lv Z, et al. Oroxin B alleviates osteoarthritis through anti-inflammation and inhibition of PI3K/AKT/MTOR signaling pathway and enhancement of autophagy. Front Endocrinol (Lausanne). 2022;13:1060721.
Meng S, Xia W, Pan M, Jia Y, He Z, Ge W. Proteomics profiling and pathway analysis of hippocampal aging in rhesus monkeys. BMC Neurosci. 2020;21(1):2.
Gong P, Wang D, Cui D, Yang Q, Wang P, Yang W, et al. Anti-aging function and molecular mechanism of radix astragali and radix astragali preparata via network pharmacology and PI3K/Akt signaling pathway. Phytomedicine. 2021;84:153509.
Schoetz U, Klein D, Hess J, Shnayien S, Spoerl S, Orth M, et al. Early senescence and production of senescence-associated cytokines are major determinants of radioresistance in head-and-neck squamous cell carcinoma. Cell Death Dis. 2021;12(12):1162.
Cheng Y, Mo F, Li Q, Han X, Shi H, Chen S, et al. Targeting CXCR2 inhibits the progression of lung cancer and promotes therapeutic effect of cisplatin. Mol Cancer. 2021;20(1):62.
Zhu YJ, Fan JJ, Wu FY, Zhang M, Song AQ, Li Y, et al. Aging promotes chronic stress-induced depressive-like behavior by activating NLRP1 inflammasome-driven inflammatory signaling in mice. Inflammation. 2022;45(6):2172–85.
Cai Q, Shi P, Yuan Y, Peng J, Ou X, Zhou W, et al. Inflammation-associated senescence promotes helicobacter pylori-induced atrophic gastritis. Cell Mol Gastroenterol Hepatol. 2021;11(3):857–80.
Patel D, Zhang X, Farrell JJ, Lunetta KL, Farrer LA. Set-based rare variant expression quantitative trait loci in blood and brain from alzheimer disease study participants. Genes (Basel). 2021;12(3):419.
Wilkinson HN, Hardman MJ. Wound senescence: a functional link between diabetes and ageing? Exp Dermatol. 2021;30(1):68–73.
Kawagoe Y, Kawashima I, Sato Y, Okamoto N, Matsubara K, Kawamura K. CXCL5-CXCR2 signaling is a senescence-associated secretory phenotype in preimplantation embryos. Aging Cell. 2020;19(10):e13240.
Takikawa T, Hamada S, Matsumoto R, Tanaka Y, Kataoka F, Sasaki A, Masamune A. Senescent human pancreatic stellate cells secrete CXCR2 agonist CXCLs to promote proliferation and migration of human pancreatic cancer AsPC-1 and MIAPaCa-2 cell lines. Int J Mol Sci. 2022;23(16):9275.
Kwon JH, Kim M, Um S, Lee HJ, Bae YK, Choi SJ, et al. Senescence-associated secretory phenotype suppression mediated by small-sized mesenchymal stem cells delays cellular senescence through TLR2 and TLR5 signaling. Cells. 2021;10(1):63.
Di Mitri D, Mirenda M, Vasilevska J, Calcinotto A, Delaleu N, Revandkar A, et al. Re-education of tumor-associated macrophages by CXCR2 blockade drives senescence and tumor inhibition in advanced prostate cancer. Cell Rep. 2019;28(8):2156-2168.e5.
Li JJ, Ma FX, Wang YW, Chen F, Lu SH, Chi Y, et al. Knockdown of IL-8 provoked premature senescence of placenta-derived mesenchymal stem cells. Stem Cells Dev. 2017;26(12):912–31.
Rao RA, Ketkar AA, Kedia N, Krishnamoorthy VK, Lakshmanan V, Kumar P, et al. KMT1 family methyltransferases regulate heterochromatin-nuclear periphery tethering via histone and non-histone protein methylation. EMBO Rep. 2019;20(5):e43260.
Lämmerhirt L, Kappelmann-Fenzl M, Fischer S, Pommer M, Zimmermann T, Kluge V, et al. K. Knockdown of Lamin B1 and the corresponding lamin B receptor leads to changes in heterochromatin state and senescence induction in malignant melanoma. Cells. 2022;11(14):2154.
Luan YY, Zhang L, Peng YQ, Li YY, Liu RX, Yin CH. Immune regulation in polycystic ovary syndrome. Clin Chim Acta. 2022;531:265–72.
Deswal R, Dang AS. Dissecting the role of Micro-RNAs as a diagnostic marker for polycystic ovary syndrome: a systematic review and meta-analysis. Fertil Steril. 2020;113(3):661-669.e2.
Mao Z, Li T, Zhao H, Qin Y, Wang X, Kang Y. Identification of epigenetic interactions between MicroRNA and DNA methylation associated with polycystic ovarian syndrome. J Hum Genet. 2021;66(2):123–37.
De Nardo MG, Baracat EC, Soares JM, Carvalho KC, Maciel GAR. Evaluation of circulating MicroRNA profiles in Brazilian women with polycystic ovary syndrome: a preliminary study. PLoS ONE. 2022;17(10):e0275031.
Zhang F, Li SP, Zhang T, Yu B, Zhang J, Ding HG, et al. High Throughput MicroRNAs sequencing profile of serum exosomes in women with and without polycystic ovarian syndrome. PeerJ. 2021;9:e10998.
Patil K, Joseph S, Shah J, Mukherjee S. An Integrated in silico analysis highlighted angiogenesis regulating MiRNA-MRNA network in PCOS pathophysiology. J Assist Reprod Genet. 2022;39(2):427–40.
Li Y, Yao N, Gao Y, Wang Y, Bai L, Xu J, et al. MiR-1224-5p attenuates polycystic ovary syndrome through inhibiting NOD-like receptor protein 3 inflammasome activation via targeting forkhead box O 1. Bioengineered. 2021;12(1):8555–69.
Wang M, Sun J, Xu B, Chrusciel M, Gao J, Bazert M, et al. Functional characterization of MicroRNA-27a-3p expression in human polycystic ovary syndrome. Endocrinology. 2018;159(1):297–309.
Liebig JK, Kuphal S, Bosserhoff AK. HuRdling Senescence: HuR Breaks BRAF-induced senescence in melanocytes and supports melanoma growth. Cancers (Basel). 2020;12(5):1299.
Liu L, Zhao C, Zhang H, Lu Y, Luo B, Yao Z, et al. Asporin regulated by MiR-26b-5p mediates chondrocyte senescence and exacerbates osteoarthritis progression via TGF-Β1/Smad2 pathway. Rheumatology (Oxford). 2022;61(6):2631–43.
Han J, Li W, Zhang J, Guan Y, Huang Y, Li X. Mechanism of CircHIPK3-MiRNA-124-3p/MiRNA-148b-3p-mediated inflammatory responses and cell senescence in Candida albicans-induced septic acute kidney injury. Gerontology. 2022;68(10):1145–65.
Markopoulos GS, Roupakia E, Tokamani M, Vartholomatos G, Tzavaras T, Hatziapostolou M, et al. Senescence-associated MicroRNAs target cell cycle regulatory genes in normal human lung fibroblasts. Exp Gerontol. 2017;96:110–22.
Terlecki-Zaniewicz L, Lämmermann I, Latreille J, Bobbili MR, Pils V, Schosserer M, Weinmüllner R, et al. Small extracellular vesicles and their MiRNA cargo are anti-apoptotic members of the senescence-associated secretory phenotype. Aging (Albany NY). 2018;10(5):1103–32.
Moimas S, Salton F, Kosmider B, Ring N, Volpe MC, Bahmed K, et al. MiR-200 family members reduce senescence and restore idiopathic pulmonary fibrosis type II alveolar epithelial cell transdifferentiation. ERJ Open Res. 2019;5(4):00138–2019.
Chevalier FP, Rorteau J, Ferraro S, Martin LS, Gonzalez-Torres A, Berthier A, et al. MiR-30a-5p alters epidermal terminal differentiation during aging by regulating BNIP3L/NIX-dependent mitophagy. Cells. 2022;11(5):836.
Jiang X, Ruan XL, Xue YX, Yang S, Shi M, Wang LN. Metformin reduces the senescence of renal tubular epithelial cells in diabetic nephropathy via the MBNL1/MiR-130a-3p/STAT3 pathway. Oxid Med Cell Longev. 2020;2020:8708236.
Borchardt H, Kogel A, Kalwa H, Weirauch U, Aigner A. Therapeutic MiR-506-3p replacement in pancreatic carcinoma leads to multiple effects including autophagy, apoptosis, senescence, and mitochondrial alterations in vitro and in vivo. Biomedicines. 2022;10(7):1692.
Choi JY, Shin HJ, Bae IH. MiR-93-5p suppresses cellular senescence by directly targeting Bcl-w and P21. Biochem Biophys Res Commun. 2018;505(4):1134–40.
Du WW, Li X, Li T, Li H, Khorshidi A, Liu F, et al. The MicroRNA MiR-17-3p inhibits mouse cardiac fibroblast senescence by targeting par4. J Cell Sci. 2015;128(2):293–304.
Yuan D, Luo J, Sun Y, Hao L, Zheng J, Yang Z. PCOS follicular fluid derived exosomal MiR-424-5p induces granulosa cells senescence by targeting CDCA4 expression. Cell Signal. 2021;85:110030.
Butler AE, Ramachandran V, Sathyapalan T, David R, Gooderham NJ, Benurwar M, et al. MicroRNA expression in women with and without polycystic ovarian syndrome matched for body mass index. Front Endocrinol (Lausanne). 2020;11:206.
Chou J, Kaller M, Jaeckel S, Rokavec M, Hermeking H. AP4 suppresses DNA damage, chromosomal instability and senescence via inducing MDC1/Mediator of DNA damage checkpoint 1 and repressing MIR22HG/MiR-22-3p. Mol Cancer. 2022;21(1):120.
Liu Y, Yang J, Yang X, Lai P, Mou Y, Deng J, et al. H2O2 down-regulates SIRT7’s protective role of endothelial premature dysfunction via MicroRNA-335-5p. Biosci Rep. 2022;42(5):BSR20211775.
Zhang S, Liu Y, Wang M, Ponikwicka-Tyszko D, Ma W, Krentowska A, et al. Role and mechanism of MiR-335-5p in the pathogenesis and treatment of polycystic ovary syndrome. Transl Res. 2022;S1931–5244(22):00172–4.
Jiang X, Li J, Zhang B, Hu J, Ma J, Cui L, et al. differential expression profile of plasma exosomal MicroRNAs in women with polycystic ovary syndrome. Fertil Steril. 2021;115(3):782–92.
Accardi G, Bono F, Cammarata G, Aiello A, Herrero MT, Alessandro R, et al. miR-126-3p and miR-21-5p as Hallmarks of bio-positive ageing; correlation analysis and machine learning prediction in young to ultra-centenarian sicilian population. Cells. 2022;11(9):1505.
Jo H-R, Jeong J-H. MicroRNA-mediated downregulation of HMGB2 contributes to cellular senescence in microvascular endothelial cells. Cells. 2022;11(3):584.
Chen L, Kong C. LINC00173 regulates polycystic ovarian syndrome progression by promoting apoptosis and repressing proliferation in ovarian granulosa cells via the MicroRNA-124-3p (MiR-124-3p)/Jagged Canonical Notch Ligand 1 (JAG1) pathway. Bioengineered. 2022;13(4):10373–85.
Huo Y, Ji S, Yang H, Wu W, Yu L, Ren Y, et al. differential expression of MicroRNA in the serum of patients with polycystic ovary syndrome with insulin resistance. Ann Transl Med. 2022;10(14):762.
Tian-Min Y, Suxia L, Shufang D, Dandan C, Long-Dan L, Shu BY, W. Combined transcriptomic and metabolomic analysis of women with polycystic ovary syndrome. Dis Markers. 2022;2022:4000424.
Romero-Ruiz A, Pineda B, Ovelleiro D, Perdices-Lopez C, Torres E, Vazquez MJ, et al. Molecular diagnosis of polycystic ovary syndrome in obese and non-obese women by targeted plasma MiRNA profiling. Eur J Endocrinol. 2021;185(5):637–52.
Zhang X, Xiao H, Zhang XEQ, Gong X, Li T, Han Y, et al. Decreased MicroRNA-125b-5p disrupts follicle steroidogenesis through targeting PAK3/ERK1/2 signalling in mouse preantral follicles. Metabolism. 2020;107:154241.
Cao J, Huo P, Cui K, Wei H, Cao J, Wang J, et al. 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. 2022;20(1):61.
Nguyen DDN, Zain SM, Kamarulzaman MH, Low TY, Chilian WM, Pan Y, et al. Intracellular and exosomal MicroRNAome profiling of human vascular smooth muscle cells during replicative senescence. Am J Physiol Heart Circ Physiol. 2021;321(4):H770–83.
Lee Y, Bae YS. Long non-coding RNA KCNQ1OT1 regulates protein kinase CK2 Via MiR-760 in senescence and calorie restriction. Int J Mol Sci. 2022;23(3):1888.
Deng L, Chen Q, Xie J, Wei W, Hui H. CircPUM1 promotes polycystic ovary syndrome progression by sponging to MiR-760. Gene. 2020;754:144903.
Tao W, Hong Y, He H, Han Q, Mao M, Hu B, et al. MicroRNA-199a-5p aggravates angiotensin II-induced vascular smooth muscle cell senescence by targeting sirtuin-1 in abdominal aortic aneurysm. J Cell Mol Med. 2021;25(13):6056–69.
Shao S, Wang H, Shao W, Liu N. MiR-199a-5p stimulates ovarian granulosa cell apoptosis in polycystic ovary syndrome. J Mol Endocrinol. 2020;65(4):187–201.
Zhang Q, Li J, Weng L. Identification and validation of aging-related genes in Alzheimer’s Disease. Front Neurosci. 2022;16:905722.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Sanz H, Valim C, Vegas E, Oller JM, Reverter F. SVM-RFE: selection and visualization of the most relevant features through non-linear kernels. BMC Bioinformatics. 2018;19(1):432.
Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, et al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12(5):453–7.
Acknowledgements
The authors acknowledge Zixiang Wang for his technical assistance.
Funding
This work was supported by the National Natural Science Foundation of China grant numbers 82071620.
Author information
Authors and Affiliations
Contributions
Study design: Dan‑ni Jiang, Yang Xu; Data retrieve: Dan‑ni Jiang, Yang Xu, Lin Yang, Xiao-juan Han; Statistical analysis: Lin Yang, Peng-fei Li, Qian‑ni Li; Result interpretation: Yang Yang; Writing‑original draft: Dan‑ni Jiang; Writing‑review & editing: Chao Lan; Funding acquisition: Chao Lan. The authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Medical Ethics Committee of Qilu Hospital, Shandong University (approval number: KYLL-202008–187).
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1:
Table S1. The primer sequences used in qRT-PCR.
Additional file 2:
Table S2. List of aging and senescence-associated genes.
Additional file 3:
Table S3. Gene lists of the WGCNA key modules and SVM-RFE algorithm.
Additional file 4:
Table S4. The edges information for TF-miRNA-mRNA regulatory network.
Additional file 5:
Table S5. The pathway information for CXCR2 and LMNB1.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jiang, D., Xu, Y., Yang, L. et al. Identification and validation of senescence-related genes in polycystic ovary syndrome. J Ovarian Res 17, 7 (2024). https://doi.org/10.1186/s13048-023-01338-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13048-023-01338-4