Abstract

The remarkable growth of multi-platform genomic profiles has led to the multiomics data integration challenge. The effective integration of such data provides a comprehensive view of the molecular complexity of cancer tumors and can significantly improve clinical out-come predictions. In this study, we present a novel network-based integration method of multiomics data as well as a clustering technique involving the Wasserstein (Earth Movers) distance from the theory of optimal mass transport. We applied our proposed method of integrative Wasserstein-based clustering (iWCluster) to invasive breast carcinoma from The Cancer Genome Atlas (TCGA) project. The subtypes were characterized by the concordant effect of mRNA expression, DNA copy number alteration, and DNA methylation as well as the interaction network connectivity of the gene products. iW-Cluster is substantially more effective in distinguishing clusters with different survival rates as compared to isolated one-dimensional conventional omics analysis. Applying iWCluster to breast cancer TCGA data successfully recovered the known PAM50 molecular subtypes. In addition, iWCluster preserves the gene-specific data, which enables us to interpret the results and perform further analysis of significant genes for a specific cluster. The gene ontology enrichment analysis of significant genes in our substantially low survival sub-group leads to the well-known phenomenon of tumor hypoxia and the transcription factor ETS1 whose expression is induced by hypoxia. Increased expression of ETS1 is associated with an increased risk of recurrence and worse prognosis in breast cancer. Consequently, we believe iWCluster has the potential to discover novel subtypes by accentuating the genes that have concordant multiomics measurements in their interaction network, which are challenging to find without the network inference or with single omics analysis.