A missense SNP in the tumor suppressor SETD2 reduces H3K36me3 and mitotic spindle integrity in Drosophila.

Mutations in SETD2 are among the most prevalent drivers of renal cell carcinoma (RCC). We identified a novel single nucleotide polymorphism (SNP) in SETD2, E902Q, within a subset of RCC patients, which manifests as both an inherited or tumor-associated somatic mutation. To determine if the SNP is biologically functional, we used CRISPR-based genome editing to generate the orthologous mutation within the Drosophila melanogaster Set2 gene. In Drosophila, the homologous amino acid substitution, E741Q, reduces H3K36me3 levels comparable to Set2 knockdown, and this loss is rescued by reintroduction of a wild-type Set2 transgene. We similarly uncovered significant defects in spindle morphogenesis, consistent with the established role of SETD2 in methylating α-Tubulin during mitosis to regulate microtubule dynamics and maintain genome stability. These data indicate the Set2 E741Q SNP affects both histone methylation and spindle integrity. Moreover, this work further suggests the SETD2 E902Q SNP may hold clinical relevance.

Systematic discovery of mutation-directed neo-protein-protein interactions in cancer.

Comprehensive sequencing of patient tumors reveals genomic mutations across tumor types that enable tumorigenesis and progression. A subset of oncogenic driver mutations results in neomorphic activity where the mutant protein mediates functions not engaged by the parental molecule. Here, we identify prevalent variant-enabled neomorph-protein-protein interactions (neoPPI) with a quantitative high-throughput differential screening (qHT-dS) platform. The coupling of highly sensitive BRET biosensors with miniaturized coexpression in an ultra-HTS format allows large-scale monitoring of the interactions of wild-type and mutant variant counterparts with a library of cancer-associated proteins in live cells. The screening of 17,792 interactions with 2,172,864 data points revealed a landscape of gain of interactions encompassing both oncogenic and tumor suppressor mutations. For example, the recurrent BRAF V600E lesion mediates KEAP1 neoPPI, rewiring a BRAFV600E/KEAP1 signaling axis and creating collateral vulnerability to NQO1 substrates, offering a combination therapeutic strategy. Thus, cancer genomic alterations can create neo-interactions, informing variant-directed therapeutic approaches for precision medicine.

NetMix: A Network-Structured Mixture Model for Reduced-Bias Estimation of Altered Subnetworks.

A classic problem in computational biology is the identification of altered subnetworks: subnetworks of an interaction network that contain genes/proteins that are differentially expressed, highly mutated, or otherwise aberrant compared with other genes/proteins. Numerous methods have been developed to solve this problem under various assumptions, but the statistical properties of these methods are often unknown. For example, some widely used methods are reported to output very large subnetworks that are difficult to interpret biologically. In this work, we formulate the identification of altered subnetworks as the problem of estimating the parameters of a class of probability distributions that we call the Altered Subset Distribution (ASD). We derive a connection between a popular method, jActiveModules, and the maximum likelihood estimator (MLE) of the ASD. We show that the MLE is statistically biased, explaining the large subnetworks output by jActiveModules. Based on these insights, we introduce NetMix, an algorithm that uses Gaussian mixture models to obtain less biased estimates of the parameters of the ASD. We demonstrate that NetMix outperforms existing methods in identifying altered subnetworks on both simulated and real data, including the identification of differentially expressed genes from both microarray and RNA-seq experiments and the identification of cancer driver genes in somatic mutation data.

Germline Features Associated with Immune Infiltration in Solid Tumors.

The immune composition of the tumor microenvironment influences response and resistance to immunotherapies. While numerous studies have identified somatic correlates of immune infiltration, germline features that associate with immune infiltrates in cancers remain incompletely characterized. We analyze seven million autosomal germline variants in the TCGA cohort and test for association with established immune-related phenotypes that describe the tumor immune microenvironment. We identify one SNP associated with the amount of infiltrating follicular helper T cells; 23 candidate genes, some of which are involved in cytokine-mediated signaling and others containing cancer-risk SNPs; and networks with genes that are part of the DNA repair and transcription elongation pathways. In addition, we find a positive association between polygenic risk for rheumatoid arthritis and amount of infiltrating CD8+ T cells. Overall, we identify multiple germline genetic features associated with tumor-immune phenotypes and develop a framework for probing inherited features that contribute to differences in immune infiltration.

Pathway and network analysis of more than 2500 whole cancer genomes.

The catalog of cancer driver mutations in protein-coding genes has greatly expanded in the past decade. However, non-coding cancer driver mutations are less well-characterized and only a handful of recurrent non-coding mutations, most notably TERT promoter mutations, have been reported. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancer across 38 tumor types, we perform multi-faceted pathway and network analyses of non-coding mutations across 2583 whole cancer genomes from 27 tumor types compiled by the ICGC/TCGA PCAWG project that was motivated by the success of pathway and network analyses in prioritizing rare mutations in protein-coding genes. While few non-coding genomic elements are recurrently mutated in this cohort, we identify 93 genes harboring non-coding mutations that cluster into several modules of interacting proteins. Among these are promoter mutations associated with reduced mRNA expression in TP53, TLE4, and TCF4. We find that biological processes had variable proportions of coding and non-coding mutations, with chromatin remodeling and proliferation pathways altered primarily by coding mutations, while developmental pathways, including Wnt and Notch, altered by both coding and non-coding mutations. RNA splicing is primarily altered by non-coding mutations in this cohort, and samples containing non-coding mutations in well-known RNA splicing factors exhibit similar gene expression signatures as samples with coding mutations in these genes. These analyses contribute a new repertoire of possible cancer genes and mechanisms that are altered by non-coding mutations and offer insights into additional cancer vulnerabilities that can be investigated for potential therapeutic treatments.

Hierarchical HotNet: identifying hierarchies of altered subnetworks.

Motivation: The analysis of high-dimensional 'omics data is often informed by the use of biological interaction networks. For example, protein-protein interaction networks have been used to analyze gene expression data, to prioritize germline variants, and to identify somatic driver mutations in cancer. In these and other applications, the underlying computational problem is to identify altered subnetworks containing genes that are both highly altered in an 'omics dataset and are topologically close (e.g. connected) on an interaction network. Results: We introduce Hierarchical HotNet, an algorithm that finds a hierarchy of altered subnetworks. Hierarchical HotNet assesses the statistical significance of the resulting subnetworks over a range of biological scales and explicitly controls for ascertainment bias in the network. We evaluate the performance of Hierarchical HotNet and several other algorithms that identify altered subnetworks on the problem of predicting cancer genes and significantly mutated subnetworks. On somatic mutation data from The Cancer Genome Atlas, Hierarchical HotNet outperforms other methods and identifies significantly mutated subnetworks containing both well-known cancer genes and candidate cancer genes that are rarely mutated in the cohort. Hierarchical HotNet is a robust algorithm for identifying altered subnetworks across different 'omics datasets. Availability and implementation: http://github.com/raphael-group/hierarchical-hotnet. Supplementary information: Supplementary material are available at Bioinformatics online.

A weighted exact test for mutually exclusive mutations in cancer.

MOTIVATION: The somatic mutations in the pathways that drive cancer development tend to be mutually exclusive across tumors, providing a signal for distinguishing driver mutations from a larger number of random passenger mutations. This mutual exclusivity signal can be confounded by high and highly variable mutation rates across a cohort of samples. Current statistical tests for exclusivity that incorporate both per-gene and per-sample mutational frequencies are computationally expensive and have limited precision. RESULTS: We formulate a weighted exact test for assessing the significance of mutual exclusivity in an arbitrary number of mutational events. Our test conditions on the number of samples with a mutation as well as per-event, per-sample mutation probabilities. We provide a recursive formula to compute P-values for the weighted test exactly as well as a highly accurate and efficient saddlepoint approximation of the test. We use our test to approximate a commonly used permutation test for exclusivity that conditions on per-event, per-sample mutation frequencies. However, our test is more efficient and it recovers more significant results than the permutation test. We use our Weighted Exclusivity Test (WExT) software to analyze hundreds of colorectal and endometrial samples from The Cancer Genome Atlas, which are two cancer types that often have extremely high mutation rates. On both cancer types, the weighted test identifies sets of mutually exclusive mutations in cancer genes with fewer false positives than earlier approaches. AVAILABILITY AND IMPLEMENTATION: See http://compbio.cs.brown.edu/projects/wext for software. CONTACT: braphael@cs.brown.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic response of RNA editing to temperature in Drosophila.

BACKGROUND: Adenosine-to-inosine RNA editing is a highly conserved process that post-transcriptionally modifies mRNA, generating proteomic diversity, particularly within the nervous system of metazoans. Transcripts encoding proteins involved in neurotransmission predominate as targets of such modifications. Previous reports suggest that RNA editing is responsive to environmental inputs in the form of temperature alterations. However, the molecular determinants underlying temperature-dependent RNA editing responses are not well understood. RESULTS: Using the poikilotherm Drosophila, we show that acute temperature alterations within a normal physiological range result in substantial changes in RNA editing levels. Our examination of particular sites reveals diversity in the patterns with which editing responds to temperature, and these patterns are conserved across five species of Drosophilidae representing over 10 million years of divergence. In addition, we show that expression of the editing enzyme, ADAR (adenosine deaminase acting on RNA), is dramatically decreased at elevated temperatures, partially, but not fully, explaining some target responses to temperature. Interestingly, this reduction in editing enzyme levels at elevated temperature is only partially reversed by a return to lower temperatures. Lastly, we show that engineered structural variants of the most temperature-sensitive editing site, in a sodium channel transcript, perturb thermal responsiveness in RNA editing profile for a particular RNA structure. CONCLUSIONS: Our results suggest that the RNA editing process responds to temperature alterations via two distinct molecular mechanisms: through intrinsic thermo-sensitivity of the RNA structures that direct editing, and due to temperature sensitive expression or stability of the RNA editing enzyme. Environmental cues, in this case temperature, rapidly reprogram the Drosophila transcriptome through RNA editing, presumably resulting in altered proteomic ratios of edited and unedited proteins.