CRISPR-Cas9 screens reveal common essential miRNAs in human cancer cell lines.

BACKGROUND: Genome-wide functional screening using the CRISPR-Cas9 system is a powerful tool to uncover tumor-specific and common genetic dependencies across cancer cell lines. Current CRISPR-Cas9 knockout libraries, however, primarily target protein-coding genes. This limits functional genomics-based investigations of miRNA function. METHODS: We designed a novel CRISPR-Cas9 knockout library (lentiG-miR) of 8107 distinct sgRNAs targeting a total of 1769 human miRNAs and benchmarked its single guide RNA (sgRNA) composition, predicted on- and off-target activity, and screening performance against previous libraries. Using a total of 45 human cancer cell lines, representing 16 different tumor entities, we performed negative selection screens to identify miRNA fitness genes. Fitness miRNAs in each cell line were scored using a combination of supervised and unsupervised essentiality classifiers. Common essential miRNAs across distinct cancer cell lines were determined using the 90th percentile method. For subsequent validation, we performed knockout experiments for selected common essential miRNAs in distinct cancer cell lines and gene expression profiling. RESULTS: We found significantly lower off-target activity for protein-coding genes and a higher miRNA gene coverage for lentiG-miR as compared to previously described miRNA-targeting libraries, while preserving high on-target activity. A minor fraction of miRNAs displayed robust depletion of targeting sgRNAs, and we observed a high level of consistency between redundant sgRNAs targeting the same miRNA gene. Across 45 human cancer cell lines, only 217 (12%) of all targeted human miRNAs scored as a fitness gene in at least one model, and fitness effects for most miRNAs were confined to small subsets of cell lines. In contrast, we identified 49 common essential miRNAs with a homogenous fitness profile across the vast majority of all cell lines. Transcriptional profiling verified highly consistent gene expression changes in response to knockout of individual common essential miRNAs across a diverse set of cancer cell lines. CONCLUSIONS: Our study presents a miRNA-targeting CRISPR-Cas9 knockout library with high gene coverage and optimized on- and off-target activities. Taking advantage of the lentiG-miR library, we define a catalogue of miRNA fitness genes in human cancer cell lines, providing the foundation for further investigation of miRNAs in human cancer.

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.

EZH2 Cooperates with BRD4-NUT to Drive NUT Carcinoma Growth by Silencing Key Tumor Suppressor Genes.

NUT carcinoma is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of progrowth genes. BET bromodomain inhibitors (BETi) are a promising treatment for NUT carcinoma that can impede BRD4-NUT's ability to activate genes, but the efficacy of BETi as monotherapy is limited. Here, we demonstrated that enhancer of zeste homolog 2 (EZH2), which silences genes through establishment of repressive chromatin, is a dependency in NUT carcinoma. Inhibition of EZH2 with the clinical compound tazemetostat potently blocked growth of NUT carcinoma cells. Epigenetic and transcriptomic analysis revealed that tazemetostat reversed the EZH2-specific H3K27me3 silencing mark and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4-NUT-regulated genes. Indeed, H3K27me3 and H3K27ac domains were found to be mutually exclusive in NUT carcinoma cells. CDKN2A was identified as the only gene among all tazemetostat-derepressed genes to confer resistance to tazemetostat in a CRISPR-Cas9 screen. Combined inhibition of EZH2 and BET synergized to downregulate cell proliferation genes, resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In preclinical models, combined tazemetostat and BETi synergistically blocked tumor growth and prolonged survival of NUT carcinoma-xenografted mice, with complete remission without relapse in one cohort. Identification of EZH2 as a dependency in NUT carcinoma substantiates the reliance of NUT carcinoma tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary, chromatin regulatory pathways to maintain NUT carcinoma growth. SIGNIFICANCE: Repression of tumor suppressor genes, including CDKN2A, by EZH2 provides a mechanistic rationale for combining EZH2 and BET inhibitors for the clinical treatment of NUT carcinoma. See related commentary by Kazansky and Kentsis, p. 3827.

EZH2 synergizes with BRD4-NUT to drive NUT carcinoma growth through silencing of key tumor suppressor genes.

NUT carcinoma (NC) is an aggressive carcinoma driven by the BRD4-NUT fusion oncoprotein, which activates chromatin to promote expression of pro-growth genes. BET bromodomain inhibitors (BETi) impede BRD4-NUT's ability to activate genes and are thus a promising treatment but limited as monotherapy. The role of gene repression in NC is unknown. Here, we demonstrate that EZH2, which silences genes through establishment of repressive chromatin, is a dependency in NC. Inhibition of EZH2 with the clinical compound tazemetostat (taz) potently blocked growth of NC cells. Epigenetic and transcriptomic analysis revealed that taz reversed the EZH2-specific H3K27me3 silencing mark, and restored expression of multiple tumor suppressor genes while having no effect on key oncogenic BRD4- NUT-regulated genes. CDKN2A was identified as the only gene amongst all taz-derepressed genes to confer resistance to taz in a CRISPR-Cas9 screen. Combined EZH2 inhibition and BET inhibition synergized to downregulate cell proliferation genes resulting in more pronounced growth arrest and differentiation than either inhibitor alone. In pre-clinical models, combined taz and BETi synergistically blocked growth and prolonged survival of NC-xenografted mice, with all mice cured in one cohort. STATEMENT OF SIGNIFICANCE: Identification of EZH2 as a dependency in NC substantiates the reliance of NC tumor cells on epigenetic dysregulation of functionally opposite, yet highly complementary chromatin regulatory pathways to maintain NC growth. In particular, repression of CDKN2A expression by EZH2 provides a mechanistic rationale for combining EZH2i with BETi for the clinical treatment of NC.

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.

Structure-function analysis of the SHOC2-MRAS-PP1C holophosphatase complex.

Receptor tyrosine kinase (RTK)-RAS signalling through the downstream mitogen-activated protein kinase (MAPK) cascade regulates cell proliferation and survival. The SHOC2-MRAS-PP1C holophosphatase complex functions as a key regulator of RTK-RAS signalling by removing an inhibitory phosphorylation event on the RAF family of proteins to potentiate MAPK signalling1. SHOC2 forms a ternary complex with MRAS and PP1C, and human germline gain-of-function mutations in this complex result in congenital RASopathy syndromes2-5. However, the structure and assembly of this complex are poorly understood. Here we use cryo-electron microscopy to resolve the structure of the SHOC2-MRAS-PP1C complex. We define the biophysical principles of holoenzyme interactions, elucidate the assembly order of the complex, and systematically interrogate the functional consequence of nearly all of the possible missense variants of SHOC2 through deep mutational scanning. We show that SHOC2 binds PP1C and MRAS through the concave surface of the leucine-rich repeat region and further engages PP1C through the N-terminal disordered region that contains a cryptic RVXF motif. Complex formation is initially mediated by interactions between SHOC2 and PP1C and is stabilized by the binding of GTP-loaded MRAS. These observations explain how mutant versions of SHOC2 in RASopathies and cancer stabilize the interactions of complex members to enhance holophosphatase activity. Together, this integrative structure-function model comprehensively defines key binding interactions within the SHOC2-MRAS-PP1C holophosphatase complex and will inform therapeutic development .

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.

Acquired Resistance to KRASG12C Inhibition in Cancer.

BACKGROUND: Clinical trials of the KRAS inhibitors adagrasib and sotorasib have shown promising activity in cancers harboring KRAS glycine-to-cysteine amino acid substitutions at codon 12 (KRASG12C). The mechanisms of acquired resistance to these therapies are currently unknown. METHODS: Among patients with KRASG12C -mutant cancers treated with adagrasib monotherapy, we performed genomic and histologic analyses that compared pretreatment samples with those obtained after the development of resistance. Cell-based experiments were conducted to study mutations that confer resistance to KRASG12C inhibitors. RESULTS: A total of 38 patients were included in this study: 27 with non-small-cell lung cancer, 10 with colorectal cancer, and 1 with appendiceal cancer. Putative mechanisms of resistance to adagrasib were detected in 17 patients (45% of the cohort), of whom 7 (18% of the cohort) had multiple coincident mechanisms. Acquired KRAS alterations included G12D/R/V/W, G13D, Q61H, R68S, H95D/Q/R, Y96C, and high-level amplification of the KRASG12C allele. Acquired bypass mechanisms of resistance included MET amplification; activating mutations in NRAS, BRAF, MAP2K1, and RET; oncogenic fusions involving ALK, RET, BRAF, RAF1, and FGFR3; and loss-of-function mutations in NF1 and PTEN. In two of nine patients with lung adenocarcinoma for whom paired tissue-biopsy samples were available, histologic transformation to squamous-cell carcinoma was observed without identification of any other resistance mechanisms. Using an in vitro deep mutational scanning screen, we systematically defined the landscape of KRAS mutations that confer resistance to KRASG12C inhibitors. CONCLUSIONS: Diverse genomic and histologic mechanisms impart resistance to covalent KRASG12C inhibitors, and new therapeutic strategies are required to delay and overcome this drug resistance in patients with cancer. (Funded by Mirati Therapeutics and others; ClinicalTrials.gov number, NCT03785249.).

Comprehensive Mutational Analysis of the BRCA1-Associated DNA Helicase and Tumor-Suppressor FANCJ/BACH1/BRIP1.

FANCJ (BRIP1/BACH1) is a hereditary breast and ovarian cancer (HBOC) gene encoding a DNA helicase. Similar to HBOC genes, BRCA1 and BRCA2, FANCJ is critical for processing DNA inter-strand crosslinks (ICL) induced by chemotherapeutics, such as cisplatin. Consequently, cells deficient in FANCJ or its catalytic activity are sensitive to ICL-inducing agents. Unfortunately, the majority of FANCJ clinical mutations remain uncharacterized, limiting therapeutic opportunities to effectively use cisplatin to treat tumors with mutated FANCJ. Here, we sought to perform a comprehensive screen to identify FANCJ loss-of-function (LOF) mutations. We developed a FANCJ lentivirus mutation library representing approximately 450 patient-derived FANCJ nonsense and missense mutations to introduce FANCJ mutants into FANCJ knockout (K/O) HeLa cells. We performed a high-throughput screen to identify FANCJ LOF mutants that, as compared with wild-type FANCJ, fail to robustly restore resistance to ICL-inducing agents, cisplatin or mitomycin C (MMC). On the basis of the failure to confer resistance to either cisplatin or MMC, we identified 26 missense and 25 nonsense LOF mutations. Nonsense mutations elucidated a relationship between location of truncation and ICL sensitivity, as the majority of nonsense mutations before amino acid 860 confer ICL sensitivity. Further validation of a subset of LOF mutations confirmed the ability of the screen to identify FANCJ mutations unable to confer ICL resistance. Finally, mapping the location of LOF mutations to a new homology model provides additional functional information. IMPLICATIONS: We identify 51 FANCJ LOF mutations, providing important classification of FANCJ mutations that will afford additional therapeutic strategies for affected patients.

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency <0.001) variants in 10/25 binding sites. In addition, we found significant signals of selection and accelerated evolution in the ACE2 coding sequence across all mammals, and specific to the bat lineage. Our results, if confirmed by additional experimental data, may lead to the identification of intermediate host species for SARS-CoV-2, guide the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Broad Host Range of SARS-CoV-2 Predicted by Comparative and Structural Analysis of ACE2 in Vertebrates.

The novel coronavirus SARS-CoV-2 is the cause of Coronavirus Disease-2019 (COVID-19). The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of 410 vertebrates, including 252 mammals, to study cross-species conservation of ACE2 and its likelihood to function as a SARS-CoV-2 receptor. We designed a five-category ranking score based on the conservation properties of 25 amino acids important for the binding between receptor and virus, classifying all species from very high to very low. Only mammals fell into the medium to very high categories, and only catarrhine primates in the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 binding, and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (<0.1%) variants in 10/25 binding sites. In addition, we observed evidence of positive selection in ACE2 in multiple species, including bats. Utilized appropriately, our results may lead to the identification of intermediate host species for SARS-CoV-2, justify the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Genetic screens in isogenic mammalian cell lines without single cell cloning.

Isogenic pairs of cell lines, which differ by a single genetic modification, are powerful tools for understanding gene function. Generating such pairs of mammalian cells, however, is labor-intensive, time-consuming, and, in some cell types, essentially impossible. Here, we present an approach to create isogenic pairs of cells that avoids single cell cloning, and screen these pairs with genome-wide CRISPR-Cas9 libraries to generate genetic interaction maps. We query the anti-apoptotic genes BCL2L1 and MCL1, and the DNA damage repair gene PARP1, identifying both expected and uncharacterized buffering and synthetic lethal interactions. Additionally, we compare acute CRISPR-based knockout, single cell clones, and small-molecule inhibition. We observe that, while the approaches provide largely overlapping information, differences emerge, highlighting an important consideration when employing genetic screens to identify and characterize potential drug targets. We anticipate that this methodology will be broadly useful to comprehensively study gene function across many contexts.

Acquired HER2 mutations in ER+ metastatic breast cancer confer resistance to estrogen receptor-directed therapies.

Seventy percent of breast cancers express the estrogen receptor (ER), and agents that target the ER are the mainstay of treatment. However, virtually all people with ER+ breast cancer develop resistance to ER-directed agents in the metastatic setting. Beyond mutations in the ER itself, which occur in 25-30% of people treated with aromatase inhibitors1-4, knowledge about clinical resistance mechanisms remains incomplete. We identified activating HER2 mutations in metastatic biopsies from eight patients with ER+ metastatic breast cancer who had developed resistance to aromatase inhibitors, tamoxifen or fulvestrant. Examination of treatment-naive primary tumors in five patients showed no evidence of pre-existing mutations in four of five patients, suggesting that these mutations were acquired under the selective pressure of ER-directed therapy. The HER2 mutations and ER mutations were mutually exclusive, suggesting a distinct mechanism of acquired resistance to ER-directed therapies. In vitro analysis confirmed that the HER2 mutations conferred estrogen independence as well as-in contrast to ER mutations-resistance to tamoxifen, fulvestrant and the CDK4 and CDK6 inhibitor palbociclib. Resistance was overcome by combining ER-directed therapy with the irreversible HER2 kinase inhibitor neratinib.

Phenotypic Characterization of a Comprehensive Set of MAPK1/ERK2 Missense Mutants.

Tumor-specific genomic information has the potential to guide therapeutic strategies and revolutionize patient treatment. Currently, this approach is limited by an abundance of disease-associated mutants whose biological functions and impacts on therapeutic response are uncharacterized. To begin to address this limitation, we functionally characterized nearly all (99.84%) missense mutants of MAPK1/ERK2, an essential effector of oncogenic RAS and RAF. Using this approach, we discovered rare gain- and loss-of-function ERK2 mutants found in human tumors, revealing that, in the context of this assay, mutational frequency alone cannot identify all functionally impactful mutants. Gain-of-function ERK2 mutants induced variable responses to RAF-, MEK-, and ERK-directed therapies, providing a reference for future treatment decisions. Tumor-associated mutations spatially clustered in two ERK2 effector-recruitment domains yet produced mutants with opposite phenotypes. This approach articulates an allele-characterization framework that can be scaled to meet the goals of genome-guided oncology.

DNA repair. Mechanism of DNA interstrand cross-link processing by repair nuclease FAN1.

DNA interstrand cross-links (ICLs) are highly toxic lesions associated with cancer and degenerative diseases. ICLs can be repaired by the Fanconi anemia (FA) pathway and through FA-independent processes involving the FAN1 nuclease. In this work, FAN1-DNA crystal structures and biochemical data reveal that human FAN1 cleaves DNA successively at every third nucleotide. In vitro, this exonuclease mechanism allows FAN1 to excise an ICL from one strand through flanking incisions. DNA access requires a 5'-terminal phosphate anchor at a nick or a 1- or 2-nucleotide flap and is augmented by a 3' flap, suggesting that FAN1 action is coupled to DNA synthesis or recombination. FAN1's mechanism of ICL excision is well suited for processing other localized DNA adducts as well.

Structure of the FANCI-FANCD2 complex: insights into the Fanconi anemia DNA repair pathway.

Fanconi anemia is a cancer predisposition syndrome caused by defects in the repair of DNA interstrand cross-links (ICLs). Central to this pathway is the Fanconi anemia I-Fanconi anemia D2 (FANCI-FANCD2) (ID) complex, which is activated by DNA damage-induced phosphorylation and monoubiquitination. The 3.4 angstrom crystal structure of the ~300 kilodalton ID complex reveals that monoubiquitination and regulatory phosphorylation sites map to the I-D interface, suggesting that they occur on monomeric proteins or an opened-up complex and that they may serve to stabilize I-D heterodimerization. The 7.8 angstrom electron-density map of FANCI-DNA crystals and in vitro data show that each protein has binding sites for both single- and double-stranded DNA, suggesting that the ID complex recognizes DNA structures that result from the encounter of replication forks with an ICL.

The ObgE/CgtA GTPase influences the stringent response to amino acid starvation in Escherichia coli.

The stringent response is important for bacterial survival under stressful conditions, such as amino acid starvation, and is characterized by the accumulation of ppGpp and pppGpp. ObgE (CgtA, YhbZ) is an essential conserved GTPase in Escherichia coli and several observations have implicated the protein in the control of the stringent response. However, consequences of the protein on specific responses to amino acid starvation have not been noted. We show that ObgE binds to ppGpp with biologically relevant affinity in vitro, implicating ppGpp as an in vivo ligand of ObgE. ObgE mutants increase the ratio of pppGpp to ppGpp within the cell during the stringent response. These changes are correlated with a delayed inhibition of DNA replication by the stringent response, delayed resumption of DNA replication after release, as well as a decreased survival after amino acid deprivation. With these data, we place ObgE as an active effector of the response to amino acid starvation in vivo. Our data correlate the pppGpp/ppGpp ratio with DNA replication control under bacterial starvation conditions, suggesting a possible role for the relative balance of these two nucleotides.

Mechanisms of recombination: lessons from E. coli.

The genetics and biochemistry of genetic recombination in E. coli has been studied for over four decades and provides a useful model system to understand recombination in other organisms. Here we provide an overview of the mechanisms of recombination and how such processes contribute to DNA repair. We describe the E. coli functions that are known to contribute to these mechanisms, step by step, and summarize their biochemical properties in relation to the role these proteins play in vivo. We feature areas of investigation that are newly emerging, as well as work that provides a historical perspective to the field. Finally, we highlight some of the questions that remain unanswered.

Chromosome segregation control by Escherichia coli ObgE GTPase.

Escherichia coli cells depleted of the conserved GTPase, ObgE, show early chromosome-partitioning defects and accumulate replicated chromosomes in which the terminus regions are colocalized. Cells lacking ObgE continue to initiate replication, with a normal ratio of the origin to terminus. Localization of the SeqA DNA binding protein, normally seen as punctate foci, however, was disturbed. Depletion of ObgE also results in cell filamentation, with polyploid DNA content. Depletion of ObgE did not cause lethality, and cells recovered fully after expression of ObgE was restored. We propose a model in which ObgE is required to license chromosome segregation and subsequent cell cycle events.