Landscape of MicroRNA Regulatory Network Architecture and Functional Rerouting in Cancer

Abstract

Somatic mutations are a major source of cancer development, and many driver mutations have been identified in protein coding regions. However, the function of mutations located in miRNA and their target binding sites throughout the human genome remains largely unknown. Here, we built detailed cancer-specific miRNA regulatory networks across 30 cancer types to systematically analyze the effect of mutations in miRNAs and their target sites in 30 untranslated region (30 UTR), coding sequence (CDS), and 50 UTR regions. A total of 3,518,261 mutations from 9,819 samples were mapped to miRNA-gene interactions (mGI). Mutations in miRNAs showed a mutually exclusive pattern with mutations in their target genes in almost all cancer types. A linear regression method identified 148 candidate driver mutations that can significantly perturb miRNA regulatory networks. Driver mutations in 30UTRs played their roles by altering RNA binding energy and the expression of target genes. Finally, mutated driver gene targets in 30 UTRs were significantly downregulated in cancer and functioned as tumor suppressors during cancer progression, suggesting potential miRNA candidates with significant clinical implications. A user-friendly, open-access web portal (mGI-map) was developed to facilitate further use of this data resource. Together, these results will facilitate novel noncoding biomarker identification and therapeutic drug design targeting the miRNA regulatory networks. Significance: A detailed miRNA-gene interaction map reveals extensive miRNA-mediated gene regulatory networks with mutation-induced perturbations across multiple cancers, serving as a resource for noncoding biomarker discovery and drug development.