Using bioinformatics and metabolomics to identify altered granulosa cells in patients with diminished ovarian reserve

Background During fertility treatment, diminished ovarian reserve (DOR) is a challenge that can seriously affect a patient's reproductive potential. However, the pathogenesis of DOR is still unclear and its treatment options are limited. This study aimed to explore DOR's molecular mechanisms. Methods We used R software to analyze the mRNA microarray dataset E-MTAB-391 downloaded from ArrayExpress, screen for differentially expressed genes (DEGs), and perform functional enrichment analyses. We also constructed the protein-protein interaction (PPI) and miRNA-mRNA networks. Ovarian granulosa cells (GCs) from women with DOR and the control group were collected to perform untargeted metabolomics analyses. Additionally, small molecule drugs were identified using the Connectivity Map database. Results We ultimately identified 138 DEGs. Our gene ontology (GO) analysis indicated that DEGs were mainly enriched in cytokine and steroid biosynthetic processes. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG), the DEGs were mainly enriched in the AGE (advanced glycation end-product)-RAGE (receptor for AGE) signaling pathway in diabetic complications and steroid biosynthesis. In the PPI network, we determined that JUN, EGR1, HMGCR, ATF3, and SQLE were hub genes that may be involved in steroid biosynthesis and inflammation. miRNAs also played a role in DOR development by regulating target genes. We validated the differences in steroid metabolism across GCs using liquid chromatography-tandem mass spectrometry (LC-MS/MS). We selected 31 small molecules with potentially positive or negative influences on DOR development. Conclusion We found that steroidogenesis and inflammation played critical roles in DOR development, and our results provide promising insights for predicting and treating DOR.