Dysregulation of the phosphatidylinositol 3-kinase (PI3K) signaling pathway occurs frequently in human cancer. PTEN tumor suppressor or PIK3CA oncogene mutations both direct PI3K-dependent tumorigenesis largely through activation of the AKT/PKB kinase. However, here we show through phosphoprotein profiling and functional genomic studies that many PIK3CA mutant cancer cell lines and human breast tumors exhibit only minimal AKT activation and a diminished reliance on AKT for anchorage-independent growth. Instead, these cells retain robust PDK1 activation and membrane localization and exhibit dependency on the PDK1 substrate SGK3. SGK3 undergoes PI3K- and PDK1-dependent activation in PIK3CA mutant cancer cells. Thus, PI3K may promote cancer through both AKT-dependent and AKT-independent mechanisms. Knowledge of differential PI3K/PDK1 signaling could inform rational therapeutics in cancers harboring PIK3CA mutations.

Constitutive expression of telomerase in human cells prevents the onset of senescence and crisis by maintaining telomere homeostasis. However, accumulating evidence suggests that the human telomerase reverse transcriptase catalytic subunit (TERT) contributes to cell physiology independently of its ability to elongate telomeres. Here we show that TERT interacts with the RNA component of mitochondrial RNA processing endoribonuclease (RMRP), a gene that is mutated in the inherited pleiotropic syndrome cartilage-hair hypoplasia. Human TERT and RMRP form a distinct ribonucleoprotein complex that has RNA-dependent RNA polymerase (RdRP) activity and produces double-stranded RNAs that can be processed into small interfering RNA in a Dicer (also known as DICER1)-dependent manner. These observations identify a mammalian RdRP composed of TERT in complex with RMRP.

The phagocyte oxidative burst, mediated by Nox2 NADPH oxidase-derived reactive oxygen species, confers host defense against a broad spectrum of bacterial and fungal pathogens. Loss-of-function mutations that impair function of the Nox2 complex result in a life-threatening immunodeficiency, and genetic variants of Nox2 subunits have been implicated in pathogenesis of inflammatory bowel disease (IBD). Thus, alterations in the oxidative burst can profoundly impact host defense, yet little is known about regulatory mechanisms that fine-tune this response. Here we report the discovery of regulatory nodes controlling oxidative burst by functional screening of genes within loci linked to human inflammatory disease. Implementing a multi-omics approach, we define transcriptional, metabolic and ubiquitin-cycling nodes controlled by Rbpj, Pfkl and Rnf145, respectively. Furthermore, we implicate Rnf145 in proteostasis of the Nox2 complex by endoplasmic reticulum-associated degradation. Consequently, ablation of Rnf145 in murine macrophages enhances bacterial clearance, and rescues the oxidative burst defects associated with Ncf4 haploinsufficiency.

Functional genomics attempts to understand the genome by perturbing the flow of information from DNA to RNA to protein, in order to learn how gene dysfunction leads to disease. CRISPR/Cas9 technology is the newest tool in the geneticist's toolbox, allowing researchers to edit DNA with unprecedented ease, speed and accuracy, and representing a novel means to perform genome-wide genetic screens to discover gene function. In this review, we first summarize the discovery and characterization of CRISPR/Cas9, and then compare it to other genome engineering technologies. We discuss its initial use in screening applications, with a focus on optimizing on-target activity and minimizing off-target effects. Finally, we comment on future challenges and opportunities afforded by this technology.

Genomic structural variation (SV), a common hallmark of cancer, has important predictive and therapeutic implications. However, accurately detecting SV using high-throughput sequencing data remains challenging, especially for 'targeted' resequencing efforts. This is critically important in the clinical setting where targeted resequencing is frequently being applied to rapidly assess clinically actionable mutations in tumor biopsies in a cost-effective manner. We present BreaKmer, a novel approach that uses a 'kmer' strategy to assemble misaligned sequence reads for predicting insertions, deletions, inversions, tandem duplications and translocations at base-pair resolution in targeted resequencing data. Variants are predicted by realigning an assembled consensus sequence created from sequence reads that were abnormally aligned to the reference genome. Using targeted resequencing data from tumor specimens with orthogonally validated SV, non-tumor samples and whole-genome sequencing data, BreaKmer had a 97.4% overall sensitivity for known events and predicted 17 positively validated, novel variants. Relative to four publically available algorithms, BreaKmer detected SV with increased sensitivity and limited calls in non-tumor samples, key features for variant analysis of tumor specimens in both the clinical and research settings.

Metastasis remains the most common cause of death in most cancers, with limited therapies for combating disseminated disease. While the primary tumour microenvironment is an important regulator of cancer progression, it is less well understood how different tissue environments influence metastasis. We analysed tumour-stroma interactions that modulate organ tropism of brain, bone and lung metastasis in xenograft models. We identified a number of potential modulators of site-specific metastasis, including cathepsin S as a regulator of breast-to-brain metastasis. High cathepsin S expression at the primary site correlated with decreased brain metastasis-free survival in breast cancer patients. Both macrophages and tumour cells produce cathepsin S, and only the combined depletion significantly reduced brain metastasis in vivo. Cathepsin S specifically mediates blood-brain barrier transmigration through proteolytic processing of the junctional adhesion molecule, JAM-B. Pharmacological inhibition of cathepsin S significantly reduced experimental brain metastasis, supporting its consideration as a therapeutic target for this disease.

Cancer cells that express oncogenic alleles of RAS typically require sustained expression of the mutant allele for survival, but the molecular basis of this oncogene dependency remains incompletely understood. To identify genes that can functionally substitute for oncogenic RAS, we systematically expressed 15,294 open reading frames in a human KRAS-dependent colon cancer cell line engineered to express an inducible KRAS-specific shRNA. We found 147 genes that promoted survival upon KRAS suppression. In particular, the transcriptional coactivator YAP1 rescued cell viability in KRAS-dependent cells upon suppression of KRAS and was required for KRAS-induced cell transformation. Acquired resistance to Kras suppression in a Kras-driven murine lung cancer model also involved increased YAP1 signaling. KRAS and YAP1 converge on the transcription factor FOS and activate a transcriptional program involved in regulating the epithelial-mesenchymal transition (EMT). Together, these findings implicate transcriptional regulation of EMT by YAP1 as a significant component of oncogenic RAS signaling.

Inhibitors that target the receptor tyrosine kinase (RTK)/Ras/mitogen-activated protein kinase (MAPK) pathway have led to clinical responses in lung and other cancers, but some patients fail to respond and in those that do resistance inevitably occurs (Balak et al., 2006; Kosaka et al., 2006; Rudin et al., 2013; Wagle et al., 2011). To understand intrinsic and acquired resistance to inhibition of MAPK signaling, we performed CRISPR-Cas9 gene deletion screens in the setting of BRAF, MEK, EGFR, and ALK inhibition. Loss of KEAP1, a negative regulator of NFE2L2/NRF2, modulated the response to BRAF, MEK, EGFR, and ALK inhibition in BRAF-, NRAS-, KRAS-, EGFR-, and ALK-mutant lung cancer cells. Treatment with inhibitors targeting the RTK/MAPK pathway increased reactive oxygen species (ROS) in cells with intact KEAP1, and loss of KEAP1 abrogated this increase. In addition, loss of KEAP1 altered cell metabolism to allow cells to proliferate in the absence of MAPK signaling. These observations suggest that alterations in the KEAP1/NRF2 pathway may promote survival in the presence of multiple inhibitors targeting the RTK/Ras/MAPK pathway.

DNA damage activates a signalling network that blocks cell-cycle progression, recruits DNA repair factors and/or triggers senescence or programmed cell death. Alterations in chromatin structure are implicated in the initiation and propagation of the DNA damage response. Here we further investigate the role of chromatin structure in the DNA damage response by monitoring ionizing-radiation-induced signalling and response events with a high-content multiplex RNA-mediated interference screen of chromatin-modifying and -interacting genes. We discover that an isoform of Brd4, a bromodomain and extra-terminal (BET) family member, functions as an endogenous inhibitor of DNA damage response signalling by recruiting the condensin II chromatin remodelling complex to acetylated histones through bromodomain interactions. Loss of this isoform results in relaxed chromatin structure, rapid cell-cycle checkpoint recovery and enhanced survival after irradiation, whereas functional gain of this isoform compacted chromatin, attenuated DNA damage response signalling and enhanced radiation-induced lethality. These data implicate Brd4, previously known for its role in transcriptional control, as an insulator of chromatin that can modulate the signalling response to DNA damage.

The transcription factor HIF1A is a key mediator of the cellular response to hypoxia. Despite the importance of HIF1A in homeostasis and various pathologies, little is known about how it regulates RNA polymerase II (RNAPII). We report here that HIF1A employs a specific variant of the Mediator complex to stimulate RNAPII elongation. The Mediator-associated kinase CDK8, but not the paralog CDK19, is required for induction of many HIF1A target genes. HIF1A induces binding of CDK8-Mediator and the super elongation complex (SEC), containing AFF4 and CDK9, to alleviate RNAPII pausing. CDK8 is dispensable for HIF1A chromatin binding and histone acetylation, but it is essential for binding of SEC and RNAPII elongation. Global analysis of active RNAPII reveals that hypoxia-inducible genes are paused and active prior to their induction. Our results provide a mechanistic link between HIF1A and CDK8, two potent oncogenes, in the cellular response to hypoxia.

The application of RNA interference (RNAi) to mammalian cells has provided the means to perform phenotypic screens to determine the functions of genes. Although RNAi has revolutionized loss-of-function genetic experiments, it has been difficult to systematically assess the prevalence and consequences of off-target effects. The Connectivity Map (CMAP) represents an unprecedented resource to study the gene expression consequences of expressing short hairpin RNAs (shRNAs). Analysis of signatures for over 13,000 shRNAs applied in 9 cell lines revealed that microRNA (miRNA)-like off-target effects of RNAi are far stronger and more pervasive than generally appreciated. We show that mitigating off-target effects is feasible in these datasets via computational methodologies to produce a consensus gene signature (CGS). In addition, we compared RNAi technology to clustered regularly interspaced short palindromic repeat (CRISPR)-based knockout by analysis of 373 single guide RNAs (sgRNAs) in 6 cells lines and show that the on-target efficacies are comparable, but CRISPR technology is far less susceptible to systematic off-target effects. These results will help guide the proper use and analysis of loss-of-function reagents for the determination of gene function.

We previously piloted the concept of a Connectivity Map (CMap), whereby genes, drugs, and disease states are connected by virtue of common gene-expression signatures. Here, we report more than a 1,000-fold scale-up of the CMap as part of the NIH LINCS Consortium, made possible by a new, low-cost, high-throughput reduced representation expression profiling method that we term L1000. We show that L1000 is highly reproducible, comparable to RNA sequencing, and suitable for computational inference of the expression levels of 81% of non-measured transcripts. We further show that the expanded CMap can be used to discover mechanism of action of small molecules, functionally annotate genetic variants of disease genes, and inform clinical trials. The 1.3 million L1000 profiles described here, as well as tools for their analysis, are available at https://clue.io.

T- and NK-cell lymphomas (TCL) are a heterogenous group of lymphoid malignancies with poor prognosis. In contrast to B-cell and myeloid malignancies, there are few preclinical models of TCLs, which has hampered the development of effective therapeutics. Here we establish and characterize preclinical models of TCL. We identify multiple vulnerabilities that are targetable with currently available agents (e.g., inhibitors of JAK2 or IKZF1) and demonstrate proof-of-principle for biomarker-driven therapies using patient-derived xenografts (PDXs). We show that MDM2 and MDMX are targetable vulnerabilities within TP53-wild-type TCLs. ALRN-6924, a stapled peptide that blocks interactions between p53 and both MDM2 and MDMX has potent in vitro activity and superior in vivo activity across 8 different PDX models compared to the standard-of-care agent romidepsin. ALRN-6924 induced a complete remission in a patient with TP53-wild-type angioimmunoblastic T-cell lymphoma, demonstrating the potential for rapid translation of discoveries from subtype-specific preclinical models.

Genome-wide association studies (GWAS) have identified thousands of variants associated with human diseases and traits. However, the majority of GWAS-implicated variants are in non-coding regions of the genome and require in depth follow-up to identify target genes and decipher biological mechanisms. Here, rather than focusing on causal variants, we have undertaken a pooled loss-of-function screen in primary hematopoietic cells to interrogate 389 candidate genes contained in 75 loci associated with red blood cell traits. Using this approach, we identify 77 genes at 38 GWAS loci, with most loci harboring 1-2 candidate genes. Importantly, the hit set was strongly enriched for genes validated through orthogonal genetic approaches. Genes identified by this approach are enriched in specific and relevant biological pathways, allowing regulators of human erythropoiesis and modifiers of blood diseases to be defined. More generally, this functional screen provides a paradigm for gene-centric follow up of GWAS for a variety of human diseases and traits.

Clear-cell carcinomas (CCCs) are a histological group of highly aggressive malignancies commonly originating in the kidney and ovary. CCCs are distinguished by aberrant lipid and glycogen accumulation and are refractory to a broad range of anti-cancer therapies. Here we identify an intrinsic vulnerability to ferroptosis associated with the unique metabolic state in CCCs. This vulnerability transcends lineage and genetic landscape, and can be exploited by inhibiting glutathione peroxidase 4 (GPX4) with small-molecules. Using CRISPR screening and lipidomic profiling, we identify the hypoxia-inducible factor (HIF) pathway as a driver of this vulnerability. In renal CCCs, HIF-2α selectively enriches polyunsaturated lipids, the rate-limiting substrates for lipid peroxidation, by activating the expression of hypoxia-inducible, lipid droplet-associated protein (HILPDA). Our study suggests targeting GPX4 as a therapeutic opportunity in CCCs, and highlights that therapeutic approaches can be identified on the basis of cell states manifested by morphological and metabolic features in hard-to-treat cancers.

CRISPR-Cas9 provides the means to perform genome editing and facilitates loss-of-function screens. However, we and others demonstrated that expression of the Cas9 endonuclease induces a gene-independent response that correlates with the number of target sequences in the genome. An alternative approach to suppressing gene expression is to block transcription using a catalytically inactive Cas9 (dCas9). Here we directly compare genome editing by CRISPR-Cas9 (cutting, CRISPRc) and gene suppression using KRAB-dCas9 (CRISPRi) in loss-of-function screens to identify cell essential genes. CRISPRc identified 98% of previously defined cell essential genes. After optimizing library construction by analysing transcriptional start sites (TSS), CRISRPi identified 92% of core cell essential genes and did not show a bias to regions involved in copy number alterations. However, bidirectional promoters scored as false positives in CRISRPi. We conclude that CRISPRc and CRISPRi have different off-target effects and combining these approaches provides complementary information in loss-of-function genetic screens.

The CRISPR-Cas9 system has revolutionized gene editing both at single genes and in multiplexed loss-of-function screens, thus enabling precise genome-scale identification of genes essential for proliferation and survival of cancer cells. However, previous studies have reported that a gene-independent antiproliferative effect of Cas9-mediated DNA cleavage confounds such measurement of genetic dependency, thereby leading to false-positive results in copy number-amplified regions. We developed CERES, a computational method to estimate gene-dependency levels from CRISPR-Cas9 essentiality screens while accounting for the copy number-specific effect. In our efforts to define a cancer dependency map, we performed genome-scale CRISPR-Cas9 essentiality screens across 342 cancer cell lines and applied CERES to this data set. We found that CERES decreased false-positive results and estimated sgRNA activity for both this data set and previously published screens performed with different sgRNA libraries. We further demonstrate the utility of this collection of screens, after CERES correction, for identifying cancer-type-specific vulnerabilities.

Most human epithelial tumors harbor numerous alterations, making it difficult to predict which genes are required for tumor survival. To systematically identify cancer dependencies, we analyzed 501 genome-scale loss-of-function screens performed in diverse human cancer cell lines. We developed DEMETER, an analytical framework that segregates on- from off-target effects of RNAi. 769 genes were differentially required in subsets of these cell lines at a threshold of six SDs from the mean. We found predictive models for 426 dependencies (55%) by nonlinear regression modeling considering 66,646 molecular features. Many dependencies fall into a limited number of classes, and unexpectedly, in 82% of models, the top biomarkers were expression based. We demonstrated the basis behind one such predictive model linking hypermethylation of the UBB ubiquitin gene to a dependency on UBC. Together, these observations provide a foundation for a cancer dependency map that facilitates the prioritization of therapeutic targets.

Heterobifunctional proteolysis-targeting chimeric compounds leverage the activity of E3 ligases to induce degradation of target oncoproteins and exhibit potent preclinical antitumor activity. To dissect the mechanisms regulating tumor cell sensitivity to different classes of pharmacological "degraders" of oncoproteins, we performed genome-scale CRISPR-Cas9-based gene editing studies. We observed that myeloma cell resistance to degraders of different targets (BET bromodomain proteins, CDK9) and operating through CRBN (degronimids) or VHL is primarily mediated by prevention of, rather than adaptation to, breakdown of the target oncoprotein; and this involves loss of function of the cognate E3 ligase or interactors/regulators of the respective cullin-RING ligase (CRL) complex. The substantial gene-level differences for resistance mechanisms to CRBN- versus VHL-based degraders explains mechanistically the lack of cross-resistance with sequential administration of these two degrader classes. Development of degraders leveraging more diverse E3 ligases/CRLs may facilitate sequential/alternating versus combined uses of these agents toward potentially delaying or preventing resistance.

Genetic and chemical perturbations impact diverse cellular phenotypes, including multiple indicators of cell health. These readouts reveal toxicity and antitumorigenic effects relevant to drug discovery and personalized medicine. We developed two customized microscopy assays, one using four targeted reagents and the other three targeted reagents, to collectively measure 70 specific cell health phenotypes including proliferation, apoptosis, reactive oxygen species, DNA damage, and cell cycle stage. We then tested an approach to predict multiple cell health phenotypes using Cell Painting, an inexpensive and scalable image-based morphology assay. In matched CRISPR perturbations of three cancer cell lines, we collected both Cell Painting and cell health data. We found that simple machine learning algorithms can predict many cell health readouts directly from Cell Painting images, at less than half the cost. We hypothesized that these models can be applied to accurately predict cell health assay outcomes for any future or existing Cell Painting dataset. For Cell Painting images from a set of 1500+ compound perturbations across multiple doses, we validated predictions by orthogonal assay readouts. We provide a web app to browse predictions: http://broad.io/cell-health-app. Our approach can be used to add cell health annotations to Cell Painting datasets.