Detection and benchmarking of somatic mutations in cancer genomes using RNA-seq data

To detect functional somatic mutations in tumor samples, whole-exome sequencing (WES) is often used for its reliability and relative low cost. RNA-seq, while generally used to measure gene expression, can potentially also be used for identification of somatic mutations. However there has been little systematic evaluation of the utility of RNA-seq for identifying somatic mutations. Here, we develop and evaluate a pipeline for processing RNA-seq data from glioblastoma multiforme (GBM) tumors in order to identify somatic mutations. The pipeline entails the use of the STAR aligner 2-pass procedure jointly with MuTect2 from genome analysis toolkit (GATK) to detect somatic variants. Variants identified from RNA-seq data were evaluated by comparison against the COSMIC and dbSNP databases, and also compared to somatic variants identified by exome sequencing. We also estimated the putative functional impact of coding variants in the most frequently mutated genes in GBM. Interestingly, variants identified by RNA-seq alone showed better representation of GBM-related mutations cataloged by COSMIC. RNA-seq-only data substantially outperformed the ability of WES to reveal potentially new somatic mutations in known GBM-related pathways, and allowed us to build a high-quality set of somatic mutations common to exome and RNA-seq calls. Using RNA-seq data in parallel with WES data to detect somatic mutations in cancer genomes can thus broaden the scope of discoveries and lend additional support to somatic variants identified by exome sequencing alone

DUSP11-mediated control of 5′-triphosphate RNA regulates RIG-I sensitivity

Deciphering the mechanisms that regulate the sensitivity of pathogen recognition receptors is imperative to understanding infection and inflammation. Here we demonstrate that the RNA triphosphatase dual-specificity phosphatase 11 (DUSP11) acts on both host and virus-derived 5'-triphosphate RNAs rendering them less active in inducing a RIG-I-mediated immune response. Reducing DUSP11 levels alters host triphosphate RNA packaged in extracellular vesicles and induces enhanced RIG-I activation in cells exposed to extracellular vesicles. Virus infection of cells lacking DUSP11 results in a higher proportion of triphosphorylated viral transcripts and attenuated virus replication, which is rescued by reducing RIG-I expression. Consistent with the activity of DUSP11 in the cellular RIG-I response, mice lacking DUSP11 display lower viral loads, greater sensitivity to triphosphorylated RNA, and a signature of enhanced interferon activity in select tissues. Our results reveal the importance of controlling 5'-triphosphate RNA levels to prevent aberrant RIG-I signaling and demonstrate DUSP11 as a key effector of this mechanism

Cell-type specific and combinatorial usage of diverse transcription factors revealed by genome-wide binding studies in multiple human cells

Cell-type diversity is governed in part by differential gene expression programs mediated by transcription factor (TF) binding. However, there are few systematic studies of the genomic binding of different types of TFs across a wide range of human cell types, especially in relation to gene expression. In the ENCODE Project, we have identified the genomic binding locations across 11 different human cell types of CTCF, RNA Pol II (RNAPII), and MYC, three TFs with diverse roles. Our data and analysis revealed how these factors bind in relation to genomic features and shape gene expression and cell-type specificity. CTCF bound predominantly in intergenic regions while RNAPII and MYC preferentially bound to core promoter regions. CTCF sites were relatively invariant across diverse cell types, while MYC showed the greatest cell-type specificity. MYC and RNAPII co-localized at many of their binding sites and putative target genes. Cell-type specific binding sites, in particular for MYC and RNAPII, were associated with cell-type specific functions. Patterns of binding in relation to gene features were generally conserved across different cell types. RNAPII occupancy was higher over exons than adjacent introns, likely reflecting a link between transcriptional elongation and splicing. TF binding was positively correlated with the expression levels of their putative target genes, but combinatorial binding, in particular of MYC and RNAPII, was even more strongly associated with higher gene expression. These data illuminate how combinatorial binding of transcription factors in diverse cell types is associated with gene expression and cell-type specific biology