Genome Sequencing and RNA-Motif Analysis Reveal Novel Damaging Noncoding Mutations in Human Tumors.

A major challenge in cancer research is to determine the biological and clinical significance of somatic mutations in noncoding regions. This has been studied in terms of recurrence, functional impact, and association to individual regulatory sites, but the combinatorial contribution of mutations to common RNA regulatory motifs has not been explored. Therefore, we developed a new method, MIRA (mutation identification for RNA alterations), to perform an unbiased and comprehensive study of significantly mutated regions (SMR) affecting binding sites for RNA-binding proteins (RBP) in cancer. Extracting signals related to RNA-related selection processes and using RNA sequencing (RNA-seq) data from the same specimens, we identified alterations in RNA expression and splicing linked to mutations on RBP binding sites. We found SRSF10 and MBNL1 motifs in introns, HNRPLL motifs at 5' UTRs, as well as 5' and 3' splice-site motifs, among others, with specific mutational patterns that disrupt the motif and impact RNA processing. MIRA facilitates the integrative analysis of multiple genome sites that operate collectively through common RBPs and aids in the interpretation of noncoding variants in cancer. MIRA is available at https://github.com/comprna/miraImplications: The study of recurrent cancer mutations on potential RBP binding sites reveals new alterations in introns, untranslated regions, and long noncoding RNAs that impact RNA processing and provide a new layer of insight that can aid in the interpretation of noncoding variants in cancer genomes. Mol Cancer Res; 16(7); 1112-24. ©2018 AACR.

A cloud-based workflow to quantify transcript-expression levels in public cancer compendia.

Public compendia of sequencing data are now measured in petabytes. Accordingly, it is infeasible for researchers to transfer these data to local computers. Recently, the National Cancer Institute began exploring opportunities to work with molecular data in cloud-computing environments. With this approach, it becomes possible for scientists to take their tools to the data and thereby avoid large data transfers. It also becomes feasible to scale computing resources to the needs of a given analysis. We quantified transcript-expression levels for 12,307 RNA-Sequencing samples from the Cancer Cell Line Encyclopedia and The Cancer Genome Atlas. We used two cloud-based configurations and examined the performance and cost profiles of each configuration. Using preemptible virtual machines, we processed the samples for as little as $0.09 (USD) per sample. As the samples were processed, we collected performance metrics, which helped us track the duration of each processing step and quantified computational resources used at different stages of sample processing. Although the computational demands of reference alignment and expression quantification have decreased considerably, there remains a critical need for researchers to optimize preprocessing steps. We have stored the software, scripts, and processed data in a publicly accessible repository (https://osf.io/gqrz9).