We aimed to develop a diagnostic model identifying ovarian cancer (OC) from benign ovarian tumors using metagenomic data from serum microbe-derived extracellular vesicles (EVs). We obtained serum samples from 166 patients with pathologically confirmed OC and 76 patients with benign ovarian tumors. For model construction and validation, samples were randomly divided into training and test sets in the ratio 2:1. Isolation of microbial EVs from serum samples of the patients and 16S rDNA amplicon sequencing were carried out. Metagenomic and clinicopathologic data-based OC diagnostic models were constructed in the training set and then validated in the test set. There were significant differences in the metagenomic profiles between the OC and benign ovarian tumor groups; specifically, genus Acinetobacter was significantly more abundant in the OC group. More importantly, Acinetobacter was the only common genus identified by seven different statistical analysis methods. Among the various metagenomic and clinicopathologic data-based OC diagnostic models, the model consisting of age, serum CA-125 levels, and relative abundance of Acinetobacter showed the best diagnostic performance with the area under the receiver operating characteristic curve of 0.898 and 0.846 in the training and test sets, respectively. Thus, our findings establish a metagenomic analysis of serum microbe-derived EVs as a potential tool for the diagnosis of OC.

Recent advances in high-throughput technology and the emergence of large-scale genomic datasets have enabled detection of genomic features that affect clinical outcomes. Although many previous computational studies have analysed the effect of each single gene or the additive effects of multiple genes on the clinical outcome, less attention has been devoted to the identification of gene-gene interactions of general type that are associated with the clinical outcome. Moreover, the integration of information from multiple molecular profiles adds another challenge to this problem. Recently, network-based approaches have gained huge popularity. However, previous network construction methods have been more concerned with the relationship between features only, rather than the effect of feature interactions on clinical outcome. We propose a mutual information-based integrative network analysis framework (MINA) that identifies gene pairs associated with clinical outcome and systematically analyses the resulting networks over multiple genomic profiles. We implement an efficient non-parametric testing scheme that ensures the significance of detected gene interactions. We develop a tool named MINA that automates the proposed analysis scheme of identifying outcome-associated gene interactions and generating various networks from those interacting pairs for downstream analysis. We demonstrate the proposed framework using real data from ovarian cancer patients in The Cancer Genome Atlas (TCGA). Statistically significant gene pairs associated with survival were identified from multiple genomic profiles, which include many individual genes that have weak or no effect on survival. Moreover, we also show that integrated networks, constructed by merging networks from multiple genomic profiles, demonstrate better topological properties and biological significance than individual networks. We have developed a simple but powerful analysis tool that is able to detect gene-gene interactions associated with clinical outcome on multiple genomic profiles. By being network-based, our approach provides a better insight into the underlying gene-gene interaction mechanisms that affect the clinical outcome of cancer patients.