Genome-scale RNAi libraries enable the systematic interrogation of gene function. However, the interpretation of RNAi screens is complicated by the observation that RNAi reagents designed to suppress the mRNA transcripts of the same gene often produce a spectrum of phenotypic outcomes due to differential on-target gene suppression or perturbation of off-target transcripts. Here we present a computational method, Analytic Technique for Assessment of RNAi by Similarity (ATARiS), that takes advantage of patterns in RNAi data across multiple samples in order to enrich for RNAi reagents whose phenotypic effects relate to suppression of their intended targets. By summarizing only such reagent effects for each gene, ATARiS produces quantitative, gene-level phenotype values, which provide an intuitive measure of the effect of gene suppression in each sample. This method is robust for data sets that contain as few as 10 samples and can be used to analyze screens of any number of targeted genes. We used this analytic approach to interrogate RNAi data derived from screening more than 100 human cancer cell lines and identified HNF1B as a transforming oncogene required for the survival of cancer cells that harbor HNF1B amplifications. ATARiS is publicly available at http://broadinstitute.org/ataris.

Glioblastoma (GBM) harbors subpopulations of therapy-resistant tumor-initiating cells (TICs) that are self-renewing and multipotent. To understand the regulation of the TIC state, we performed an image-based screen for genes regulating GBM TIC maintenance and identified ZFHX4, a 397 kDa transcription factor. ZFHX4 is required to maintain TIC-associated and normal human neural precursor cell phenotypes in vitro, suggesting that ZFHX4 regulates differentiation, and its suppression increases glioma-free survival in intracranial xenografts. ZFHX4 interacts with CHD4, a core member of the nucleosome remodeling and deacetylase (NuRD) complex. ZFHX4 and CHD4 bind to overlapping sets of genomic loci and control similar gene expression programs. Using expression data derived from GBM patients, we found that ZFHX4 significantly affects CHD4-mediated gene expression perturbations, which defines ZFHX4 as a master regulator of CHD4. These observations define ZFHX4 as a regulatory factor that links the chromatin-remodeling NuRD complex and the GBM TIC state.

Genome-wide screening in human and mouse cells using RNA interference and open reading frame over-expression libraries is rapidly becoming a viable experimental approach for many research labs. There are a variety of gene expression modulation libraries commercially available, however, detailed and validated protocols as well as the reagents necessary for deconvolving genome-scale gene screens using these libraries are lacking. As a solution, we designed a comprehensive platform for highly multiplexed functional genetic screens in human, mouse and yeast cells using popular, commercially available gene modulation libraries. The Gene Modulation Array Platform (GMAP) is a single microarray-based detection solution for deconvolution of loss and gain-of-function pooled screens.

Trans-splicing is an unusual process in which two separate RNA strands are spliced together to yield a mature mRNA. We present a novel computational approach which has an overall accuracy of 82% and can predict 92% of known trans-splicing sites. We have applied our method to chromosomes 1 and 3 of Leishmania major, with high-confidence predictions for 85% and 88% of annotated genes respectively. We suggest some extensions of our method to other systems.

Despite remarkable success of immune checkpoint inhibitors, the majority of cancer patients have yet to receive durable benefits. Here, in order to investigate the metabolic alterations in response to immune checkpoint blockade, we comprehensively profile serum metabolites in advanced melanoma and renal cell carcinoma patients treated with nivolumab, an antibody against programmed cell death protein 1 (PD1). We identify serum kynurenine/tryptophan ratio increases as an adaptive resistance mechanism associated with worse overall survival. This advocates for patient stratification and metabolic monitoring in immunotherapy clinical trials including those combining PD1 blockade with indoleamine 2,3-dioxygenase/tryptophan 2,3-dioxygenase   (IDO/TDO) inhibitors.

Using a genome-scale, lentivirally delivered shRNA library, we performed massively parallel pooled shRNA screens in 216 cancer cell lines to identify genes that are required for cell proliferation and/or viability. Cell line dependencies on 11,000 genes were interrogated by 5 shRNAs per gene. The proliferation effect of each shRNA in each cell line was assessed by transducing a population of 11M cells with one shRNA-virus per cell and determining the relative enrichment or depletion of each of the 54,000 shRNAs after 16 population doublings using Next Generation Sequencing. All the cell lines were screened using standardized conditions to best assess differential genetic dependencies across cell lines. When combined with genomic characterization of these cell lines, this dataset facilitates the linkage of genetic dependencies with specific cellular contexts (e.g., gene mutations or cell lineage). To enable such comparisons, we developed and provided a bioinformatics tool to identify linear and nonlinear correlations between these features.