Abstract
This study aimed to identify microRNAs (miRNAs) whose expression levels are altered by high-risk human papillomavirus (HR-HPV) infection in women with epithelial ovarian neoplasms. MiRNA expression was quantified by real-time polymerase chain reaction, while HR-HPV DNA was quantified using digital-droplet PCR. Analysis of 11 miRNAs demonstrated significantly lower hsa-miR-25-5p expression in HPV-infected compared to uninfected ovarian tissues (p = 0.0405), while differences in miRNA expression in corresponding serum were statistically insignificant. The expression of hsa-miR-218-5p in ovarian tumors was significantly higher in high-grade serous ovarian carcinoma (HGSOC) cases than in other neoplasms (p = 0.0166). In addition, hsa-miR-218-5p was significantly upregulated, whereas hsa-miR-191-5p was significantly downregulated in tissues with stage III/IV FIGO (p = 0.0009 and p = 0.0305, respectively). Using unsupervised clustering, we identified three unique patient groups with significantly varied frequencies of HPV16/18-positive samples and varied miRNA expression profiles. In multivariate analysis, high expression of hsa-miR-16-5p was an independent prognostic factor for poor overall survival (p = 0.0068). This preliminary analysis showed the changes in miRNA expression in ovarian neoplasms during HPV infection and those collected from HGSOCs or patients with advanced disease. This prospective study can provide new insights into the pathogenesis of ovarian neoplasms and host-virus interactions.