Three million images and morphological profiles of cells treated with matched chemical and genetic perturbations.

The identification of genetic and chemical perturbations with similar impacts on cell morphology can elucidate compounds' mechanisms of action or novel regulators of genetic pathways. Research on methods for identifying such similarities has lagged due to a lack of carefully designed and well-annotated image sets of cells treated with chemical and genetic perturbations. Here we create such a Resource dataset, CPJUMP1, in which each perturbed gene's product is a known target of at least two chemical compounds in the dataset. We systematically explore the directionality of correlations among perturbations that target the same protein encoded by a given gene, and we find that identifying matches between chemical and genetic perturbations is a challenging task. Our dataset and baseline analyses provide a benchmark for evaluating methods that measure perturbation similarities and impact, and more generally, learn effective representations of cellular state from microscopy images. Such advancements would accelerate the applications of image-based profiling of cellular states, such as uncovering drug mode of action or probing functional genomics.

Comprehensive mutational scanning of EGFR reveals TKI sensitivities of extracellular domain mutants.

The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We perform saturation mutagenesis of EGFR and assess function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach reveals enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, not associated with approved targeted therapies, are sensitive to afatinib and dacomitinib in vitro. Two glioblastoma patients with somatic EGFR G598V dimerization domain mutations show responses to dacomitinib treatment followed by within-pathway resistance mutation in one case. In summary, this comprehensive screen expands the landscape of functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations.

Modular vector assembly enables rapid assessment of emerging CRISPR technologies.

The diversity of CRISPR systems, coupled with scientific ingenuity, has led to an explosion of applications; however, to test newly described innovations in their model systems, researchers typically embark on cumbersome, one-off cloning projects to generate custom reagents that are optimized for their biological questions. Here, we leverage Golden Gate cloning to create the Fragmid toolkit, a modular set of CRISPR cassettes and delivery technologies, along with a web portal, resulting in a combinatorial platform that enables scalable vector assembly within days. We further demonstrate that multiple CRISPR technologies can be assessed in parallel in a pooled screening format using this resource, enabling the rapid optimization of both novel technologies and cellular models. These results establish Fragmid as a robust system for the rapid design of CRISPR vectors, and we anticipate that this assembly approach will be broadly useful for systematic development, comparison, and dissemination of CRISPR technologies.

Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma.

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames (ORFs). To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a stepwise approach using multiple CRISPR-Cas9 screens to elucidate non-canonical ORFs and putative microproteins implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream ORFs (uORFs) exhibited selective functionality independent of main coding sequences. A microprotein encoded by one of these ORFs, ASNSD1-uORF or ASDURF, was upregulated, associated with MYC-family oncogenes, and promoted medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future studies seeking to define new cancer targets.

Modular vector assembly enables rapid assessment of emerging CRISPR technologies.

The diversity of CRISPR systems, coupled with scientific ingenuity, has led to an explosion of applications; however, to test newly-described innovations in their model systems, researchers typically embark on cumbersome, one-off cloning projects to generate custom reagents that are optimized for their biological questions. Here, we leverage Golden Gate cloning to create the Fragmid toolkit, a modular set of CRISPR cassettes and delivery technologies, along with a web portal, resulting in a combinatorial platform that enables scalable vector assembly within days. We further demonstrate that multiple CRISPR technologies can be assessed in parallel in a pooled screening format using this resource, enabling the rapid optimization of both novel technologies and cellular models. These results establish Fragmid as a robust system for the rapid design of CRISPR vectors, and we anticipate that this assembly approach will be broadly useful for systematic development, comparison, and dissemination of CRISPR technologies.

Optimization of Cas12a for multiplexed genome-scale transcriptional activation.

Cas12a CRISPR technology, unlike Cas9, allows for facile multiplexing of guide RNAs from a single transcript, simplifying combinatorial perturbations. While Cas12a has been implemented for multiplexed knockout genetic screens, it has yet to be optimized for CRISPR activation (CRISPRa) screens in human cells. Here, we develop a new Cas12a-based transactivation domain (TAD) recruitment system using the ALFA nanobody and demonstrate simultaneous activation of up to four genes. We screen a genome-wide library to identify modulators of growth and MEK inhibition, and we compare these results with those obtained with open reading frame (ORF) overexpression and Cas9-based CRISPRa. We find that the activity of multiplexed arrays is largely predictable from the best-performing guide and provide criteria for selecting active guides. We anticipate that these results will greatly accelerate the exploration of gene function and combinatorial phenotypes at scale.

Genome-wide pooled CRISPR screening in neurospheres.

Spheroid culture systems have allowed in vitro propagation of cells unable to grow in canonical cell culturing conditions, and may capture cellular contexts that model tumor growth better than current model systems. The insights gleaned from genome-wide clustered regularly interspaced short palindromic repeat (CRISPR) screening of thousands of cancer cell lines grown in conventional culture conditions illustrate the value of such CRISPR pooled screens. It is clear that similar genome-wide CRISPR screens of three-dimensional spheroid cultures will be important for future biological discovery. Here, we present a protocol for genome-wide CRISPR screening of three-dimensional neurospheres. While many in-depth protocols and discussions have been published for more typical cell lines, few detailed protocols are currently available in the literature for genome-wide screening in spheroidal cell lines. For those who want to screen such cell lines, and particularly neurospheres, we provide a step-by-step description of assay development tests to be performed before screening, as well as for the screen itself. We highlight considerations of variables that make these screens distinct from, or similar to, typical nonspheroid cell lines throughout. Finally, we illustrate typical outcomes of neurosphere genome-wide screens, and how neurosphere screens typically produce slightly more heterogeneous signal distributions than more canonical cancer cell lines. Completion of this entire protocol will take 8-12 weeks from the initial assay development tests to deconvolution of the sequencing data.

Optimizing the Cell Painting assay for image-based profiling.

In image-based profiling, software extracts thousands of morphological features of cells from multi-channel fluorescence microscopy images, yielding single-cell profiles that can be used for basic research and drug discovery. Powerful applications have been proven, including clustering chemical and genetic perturbations on the basis of their similar morphological impact, identifying disease phenotypes by observing differences in profiles between healthy and diseased cells and predicting assay outcomes by using machine learning, among many others. Here, we provide an updated protocol for the most popular assay for image-based profiling, Cell Painting. Introduced in 2013, it uses six stains imaged in five channels and labels eight diverse components of the cell: DNA, cytoplasmic RNA, nucleoli, actin, Golgi apparatus, plasma membrane, endoplasmic reticulum and mitochondria. The original protocol was updated in 2016 on the basis of several years' experience running it at two sites, after optimizing it by visual stain quality. Here, we describe the work of the Joint Undertaking for Morphological Profiling Cell Painting Consortium, to improve upon the assay via quantitative optimization by measuring the assay's ability to detect morphological phenotypes and group similar perturbations together. The assay gives very robust outputs despite various changes to the protocol, and two vendors' dyes work equivalently well. We present Cell Painting version 3, in which some steps are simplified and several stain concentrations can be reduced, saving costs. Cell culture and image acquisition take 1-2 weeks for typically sized batches of ≤20 plates; feature extraction and data analysis take an additional 1-2 weeks.This protocol is an update to Nat. Protoc. 11, 1757-1774 (2016): https://doi.org/10.1038/nprot.2016.105.

Translation of non-canonical open reading frames as a cancer cell survival mechanism in childhood medulloblastoma.

A hallmark of high-risk childhood medulloblastoma is the dysregulation of RNA translation. Currently, it is unknown whether medulloblastoma dysregulates the translation of putatively oncogenic non-canonical open reading frames. To address this question, we performed ribosome profiling of 32 medulloblastoma tissues and cell lines and observed widespread non-canonical ORF translation. We then developed a step-wise approach to employ multiple CRISPR-Cas9 screens to elucidate functional non-canonical ORFs implicated in medulloblastoma cell survival. We determined that multiple lncRNA-ORFs and upstream open reading frames (uORFs) exhibited selective functionality independent of the main coding sequence. One of these, ASNSD1-uORF or ASDURF, was upregulated, associated with the MYC family oncogenes, and was required for medulloblastoma cell survival through engagement with the prefoldin-like chaperone complex. Our findings underscore the fundamental importance of non-canonical ORF translation in medulloblastoma and provide a rationale to include these ORFs in future cancer genomics studies seeking to define new cancer targets.

Mapping the landscape of genetic dependencies in chordoma.

Identifying the spectrum of genes required for cancer cell survival can reveal essential cancer circuitry and therapeutic targets, but such a map remains incomplete for many cancer types. We apply genome-scale CRISPR-Cas9 loss-of-function screens to map the landscape of selectively essential genes in chordoma, a bone cancer with few validated targets. This approach confirms a known chordoma dependency, TBXT (T; brachyury), and identifies a range of additional dependencies, including PTPN11, ADAR, PRKRA, LUC7L2, SRRM2, SLC2A1, SLC7A5, FANCM, and THAP1. CDK6, SOX9, and EGFR, genes previously implicated in chordoma biology, are also recovered. We find genomic and transcriptomic features that predict specific dependencies, including interferon-stimulated gene expression, which correlates with ADAR dependence and is elevated in chordoma. Validating the therapeutic relevance of dependencies, small-molecule inhibitors of SHP2, encoded by PTPN11, have potent preclinical efficacy against chordoma. Our results generate an emerging map of chordoma dependencies to enable biological and therapeutic hypotheses.

Transcriptional Antagonism by CDK8 Inhibition Improves Therapeutic Efficacy of MEK Inhibitors.

Aberrant RAS/MAPK signaling is a common driver of oncogenesis that can be therapeutically targeted with clinically approved MEK inhibitors. Disease progression on single-agent MEK inhibitors is common, however, and combination therapies are typically required to achieve significant clinical benefit in advanced cancers. Here we focused on identifying MEK inhibitor-based combination therapies in neuroblastoma with mutations that activate the RAS/MAPK signaling pathway, which are rare at diagnosis but frequent in relapsed neuroblastoma. A genome-scale CRISPR-Cas9 functional genomic screen was deployed to identify genes that when knocked out sensitize RAS-mutant neuroblastoma to MEK inhibition. Loss of either CCNC or CDK8, two members of the mediator kinase module, sensitized neuroblastoma to MEK inhibition. Furthermore, small-molecule kinase inhibitors of CDK8 improved response to MEK inhibitors in vitro and in vivo in RAS-mutant neuroblastoma and other adult solid tumors. Transcriptional profiling revealed that loss of CDK8 or CCNC antagonized the transcriptional signature induced by MEK inhibition. When combined, loss of CDK8 or CCNC prevented the compensatory upregulation of progrowth gene expression induced by MEK inhibition. These findings propose a new therapeutic combination for RAS-mutant neuroblastoma and may have clinical relevance for other RAS-driven malignancies. SIGNIFICANCE: Transcriptional adaptation to MEK inhibition is mediated by CDK8 and can be blocked by the addition of CDK8 inhibitors to improve response to MEK inhibitors in RAS-mutant neuroblastoma, a clinically challenging disease.

BAF Complex Maintains Glioma Stem Cells in Pediatric H3K27M Glioma.

Diffuse midline gliomas are uniformly fatal pediatric central nervous system cancers that are refractory to standard-of-care therapeutic modalities. The primary genetic drivers are a set of recurrent amino acid substitutions in genes encoding histone H3 (H3K27M), which are currently undruggable. These H3K27M oncohistones perturb normal chromatin architecture, resulting in an aberrant epigenetic landscape. To interrogate for epigenetic dependencies, we performed a CRISPR screen and show that patient-derived H3K27M-glioma neurospheres are dependent on core components of the mammalian BAF (SWI/SNF) chromatin remodeling complex. The BAF complex maintains glioma stem cells in a cycling, oligodendrocyte precursor cell-like state, in which genetic perturbation of the BAF catalytic subunit SMARCA4 (BRG1), as well as pharmacologic suppression, opposes proliferation, promotes progression of differentiation along the astrocytic lineage, and improves overall survival of patient-derived xenograft models. In summary, we demonstrate that therapeutic inhibition of the BAF complex has translational potential for children with H3K27M gliomas. SIGNIFICANCE: Epigenetic dysregulation is at the core of H3K27M-glioma tumorigenesis. Here, we identify the BRG1-BAF complex as a critical regulator of enhancer and transcription factor landscapes, which maintain H3K27M glioma in their progenitor state, precluding glial differentiation, and establish pharmacologic targeting of the BAF complex as a novel treatment strategy for pediatric H3K27M glioma. See related commentary by Beytagh and Weiss, p. 2730. See related article by Mo et al., p. 2906.

PPM1D mutations are oncogenic drivers of de novo diffuse midline glioma formation.

The role of PPM1D mutations in de novo gliomagenesis has not been systematically explored. Here we analyze whole genome sequences of 170 pediatric high-grade gliomas and find that truncating mutations in PPM1D that increase the stability of its phosphatase are clonal driver events in 11% of Diffuse Midline Gliomas (DMGs) and are enriched in primary pontine tumors. Through the development of DMG mouse models, we show that PPM1D mutations potentiate gliomagenesis and that PPM1D phosphatase activity is required for in vivo oncogenesis. Finally, we apply integrative phosphoproteomic and functional genomics assays and find that oncogenic effects of PPM1D truncation converge on regulators of cell cycle, DNA damage response, and p53 pathways, revealing therapeutic vulnerabilities including MDM2 inhibition.

Noncanonical open reading frames encode functional proteins essential for cancer cell survival.

Although genomic analyses predict many noncanonical open reading frames (ORFs) in the human genome, it is unclear whether they encode biologically active proteins. Here we experimentally interrogated 553 candidates selected from noncanonical ORF datasets. Of these, 57 induced viability defects when knocked out in human cancer cell lines. Following ectopic expression, 257 showed evidence of protein expression and 401 induced gene expression changes. Clustered regularly interspaced short palindromic repeat (CRISPR) tiling and start codon mutagenesis indicated that their biological effects required translation as opposed to RNA-mediated effects. We found that one of these ORFs, G029442-renamed glycine-rich extracellular protein-1 (GREP1)-encodes a secreted protein highly expressed in breast cancer, and its knockout in 263 cancer cell lines showed preferential essentiality in breast cancer-derived lines. The secretome of GREP1-expressing cells has an increased abundance of the oncogenic cytokine GDF15, and GDF15 supplementation mitigated the growth-inhibitory effect of GREP1 knockout. Our experiments suggest that noncanonical ORFs can express biologically active proteins that are potential therapeutic targets.

Genotype-Fitness Maps of EGFR-Mutant Lung Adenocarcinoma Chart the Evolutionary Landscape of Resistance for Combination Therapy Optimization.

Cancer evolution poses a central obstacle to cure, as resistant clones expand under therapeutic selection pressures. Genome sequencing of relapsed disease can nominate genomic alterations conferring resistance but sample collection lags behind, limiting therapeutic innovation. Genome-wide screens offer a complementary approach to chart the compendium of escape genotypes, anticipating clinical resistance. We report genome-wide open reading frame (ORF) resistance screens for first- and third-generation epidermal growth factor receptor (EGFR) inhibitors and a MEK inhibitor. Using serial sampling, dose gradients, and mathematical modeling, we generate genotype-fitness maps across therapeutic contexts and identify alterations that escape therapy. Our data expose varying dose-fitness relationship across genotypes, ranging from complete dose invariance to paradoxical dose dependency where fitness increases in higher doses. We predict fitness with combination therapy and compare these estimates to genome-wide fitness maps of drug combinations, identifying genotypes where combination therapy results in unexpected inferior effectiveness. These data are applied to nominate combination optimization strategies to forestall resistant disease.

Resistance Mechanisms to SYK Inhibition in Acute Myeloid Leukemia.

Spleen tyrosine kinase (SYK) is a nonmutated therapeutic target in acute myeloid leukemia (AML). Attempts to exploit SYK therapeutically in AML have shown promising results in combination with chemotherapy, likely reflecting induced mechanisms of resistance to single-agent treatment in vivo. We conducted a genome-scale open reading frame (ORF) resistance screen and identified activation of the RAS-MAPK-ERK pathway as one major mechanism of resistance to SYK inhibitors. This finding was validated in AML cell lines with innate and acquired resistance to SYK inhibitors. Furthermore, patients with AML with select mutations activating these pathways displayed early resistance to SYK inhibition. To circumvent SYK inhibitor therapy resistance in AML, we demonstrate that a MEK and SYK inhibitor combination is synergistic in vitro and in vivo. Our data provide justification for use of ORF screening to identify resistance mechanisms to kinase inhibitor therapy in AML lacking distinct mutations and to direct novel combination-based strategies to abrogate these. SIGNIFICANCE: The integration of functional genomic screening with the study of mechanisms of intrinsic and acquired resistance in model systems and human patients identified resistance to SYK inhibitors through MAPK signaling in AML. The dual targeting of SYK and the MAPK pathway offers a combinatorial strategy to overcome this resistance.This article is highlighted in the In This Issue feature, p. 161.

A Compendium of Genetic Modifiers of Mitochondrial Dysfunction Reveals Intra-organelle Buffering.

Mitochondrial dysfunction is associated with a spectrum of human conditions, ranging from rare, inborn errors of metabolism to the aging process. To identify pathways that modify mitochondrial dysfunction, we performed genome-wide CRISPR screens in the presence of small-molecule mitochondrial inhibitors. We report a compendium of chemical-genetic interactions involving 191 distinct genetic modifiers, including 38 that are synthetic sick/lethal and 63 that are suppressors. Genes involved in glycolysis (PFKP), pentose phosphate pathway (G6PD), and defense against lipid peroxidation (GPX4) scored high as synthetic sick/lethal. A surprisingly large fraction of suppressors are pathway intrinsic and encode mitochondrial proteins. A striking example of such "intra-organelle" buffering is the alleviation of a chemical defect in complex V by simultaneous inhibition of complex I, which benefits cells by rebalancing redox cofactors, increasing reductive carboxylation, and promoting glycolysis. Perhaps paradoxically, certain forms of mitochondrial dysfunction may best be buffered with "second site" inhibitors to the organelle.

Synthetic Lethal Interaction of SHOC2 Depletion with MEK Inhibition in RAS-Driven Cancers.

The mitogen-activated protein kinase (MAPK) pathway is a critical effector of oncogenic RAS signaling, and MAPK pathway inhibition may be an effective combination treatment strategy. We performed genome-scale loss-of-function CRISPR-Cas9 screens in the presence of a MEK1/2 inhibitor (MEKi) in KRAS-mutant pancreatic and lung cancer cell lines and identified genes that cooperate with MEK inhibition. While we observed heterogeneity in genetic modifiers of MEKi sensitivity across cell lines, several recurrent classes of synthetic lethal vulnerabilities emerged at the pathway level. Multiple members of receptor tyrosine kinase (RTK)-RAS-MAPK pathways scored as sensitizers to MEKi. In particular, we demonstrate that knockout, suppression, or degradation of SHOC2, a positive regulator of MAPK signaling, specifically cooperated with MEK inhibition to impair proliferation in RAS-driven cancer cells. The depletion of SHOC2 disrupted survival pathways triggered by feedback RTK signaling in response to MEK inhibition. Thus, these findings nominate SHOC2 as a potential target for combination therapy.

Pooled Genomic Screens Identify Anti-apoptotic Genes as Targetable Mediators of Chemotherapy Resistance in Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) is often sensitive to initial treatment with platinum and taxane combination chemotherapy, but most patients relapse with chemotherapy-resistant disease. To systematically identify genes modulating chemotherapy response, we performed pooled functional genomic screens in HGSOC cell lines treated with cisplatin, paclitaxel, or cisplatin plus paclitaxel. Genes in the intrinsic pathway of apoptosis were among the top candidate resistance genes in both gain-of-function and loss-of-function screens. In an open reading frame overexpression screen, followed by a mini-pool secondary screen, anti-apoptotic genes including BCL2L1 (BCL-XL) and BCL2L2 (BCL-W) were associated with chemotherapy resistance. In a CRISPR-Cas9 knockout screen, loss of BCL2L1 decreased cell survival whereas loss of proapoptotic genes promoted resistance. To dissect the role of individual anti-apoptotic proteins in HGSOC chemotherapy response, we evaluated overexpression or inhibition of BCL-2, BCL-XL, BCL-W, and MCL1 in HGSOC cell lines. Overexpression of anti-apoptotic proteins decreased apoptosis and modestly increased cell viability upon cisplatin or paclitaxel treatment. Conversely, specific inhibitors of BCL-XL, MCL1, or BCL-XL/BCL-2, but not BCL-2 alone, enhanced cell death when combined with cisplatin or paclitaxel. Anti-apoptotic protein inhibitors also sensitized HGSOC cells to the poly (ADP-ribose) polymerase inhibitor olaparib. These unbiased screens highlight anti-apoptotic proteins as mediators of chemotherapy resistance in HGSOC, and support inhibition of BCL-XL and MCL1, alone or combined with chemotherapy or targeted agents, in treatment of primary and recurrent HGSOC. IMPLICATIONS: Anti-apoptotic proteins modulate drug resistance in ovarian cancer, and inhibitors of BCL-XL or MCL1 promote cell death in combination with chemotherapy.