Specific expression profile of follicular fluid-derived exosomal microRNAs in patients with diminished ovarian reserve

Abstract Background Diminished ovarian reserve (DOR) is defined as a reduction in ovarian reserve and oocyte quality. The pathophysiology of DOR has not been completely explained as of yet. Scholars have uncovered a large number of exosomes that have been detected in follicular fluid, and exosomal miRNAs have been proven to play a critical role in controlling ovarian disorders and follicle formation. We focused on the expression profile of follicular fluid-derived exosomal microRNAs (miRNAs) and attempted to understand if their role is connected to the pathomechanism of DOR. Methods The follicular fluid-derived differentially expressed exosomal miRNAs (DEmiRs) between patients with DOR and those with normal ovarian function were investigated using the next-generation sequencing (NGS) method. The main metabolic and signaling pathways of DEmiRs were identified using the KEGG pathway database, disease ontology (DO) analysis, and gene ontology (GO) analysis. In the end, a Protein-Protein Interaction (PPI) network was built to search for exosomal miRNAs and their target genes that were potentially strongly connected with DOR. Results In comparison to normal controls, 52 DEmiRs were discovered in follicular fluid-derived exosomes of DOR patients, of which 19 were up-regulated and 33 were down-regulated (|log2(fold change) |>2, P < 0.05). GO, DO analysis, and the KEGG pathway database revealed that many of these DEmiRs have broad biological roles that are connected to ovarian function and disorders. The top ten DEmiRs in terms of expression were then chosen for miRNA-mRNA interaction analysis. Totally, 8 experimentally supported miRNAs (hsa-miR-1246, hsa-miR-483-3p, hsa-miR-122-5p, hsa-miR-130b-3p, hsa-miR-342-3p, hsa-miR-625-3p, hsa-miR-675-3p, and hsa-miR-134-5p) and 126 target genes were filtrated by utilizing Cytoscape software. The module analysis findings of the PPI network showed that the main module cluster with a score > 6.0 (MCODE score = 15) had six hub genes, including IGFR, VEGFA, KRAS, ERBB2, RHOA, and PTEN (MCODE score = 11.472). Conclusion Our data suggested a special expression profile of follicular fluid-derived exosomal miRNAs in patients with DOR, which was probably correlated to ovarian dysfunction and follicle formation. These results may give a unique insight into a better understanding of the molecular process in the pathogenesis of DOR or other ovarian diseases.