Screening of potential biomarkers for polycystic ovary syndrome and identification of expression and immune characteristics

Background: Polycystic ovary syndrome (PCOS) seriously affects the fertility and health of women of childbearing age. We look forward to finding potential biomarkers for PCOS that can aid clinical diagnosis. Methods: We acquired PCOS and normal granulosa cell (GC) expression profiles from the Gene Expression Omnibus (GEO) database. After data preprocessing, differentially expressed genes (DEGs) were screened by limma package, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis and Gene Set Enrichment Analysis (GSEA) were performed. Recursive feature elimination (RFE) algorithm and the least absolute shrinkage and selection operator (LASSO) Cox regression analysis were used to acquire feature genes as potential biomarkers. Time-dependent receiver operator characteristic curve (ROC curve) and Confusion matrix were used to verify the classification performance of biomarkers. Then, the expression characteristics of biomarkers in PCOS and normal cells were analyzed, and the insulin resistance (IR) score of samples was computed by ssGSEA. Immune characterization of biomarkers was evaluated using MCP counter and single sample gene set enrichment analysis (ssGSEA). Finally, the correlation between biomarkers and the scores of each pathway was assessed. Results: We acquired 93 DEGs, and the enrichment results indicated that most of DEGs in PCOS group were significantly enriched in immune-related biological pathways. Further screening results indicated that JDP2 and HMOX1 were potential biomarkers. The area under ROC curve (AUC) value and Confusion matrix of the two biomarkers were ideal when separated and combined. In the combination, the training set AUC = 0.929 and the test set AUC = 0.917 indicated good diagnostic performance of the two biomarkers. Both biomarkers were highly expressed in the PCOS group, and both biomarkers, which should be suppressed in the preovulation phase, were elevated in PCOS tissues. The IR score of PCOS group was higher, and the expression of JDP2 and HMOX1 showed a significant positive correlation with IR score. Most immune cell scores and immune infiltration results were significantly higher in PCOS. Comprehensive analysis indicated that the two biomarkers had strong correlation with immune-related pathways. Conclusion: We acquired two potential biomarkers, JDP2 and HMOX1. We found that they were highly expressed in the PCOS and had a strong positive correlation with immune-related pathways.