Functional Genomic Screens Research Articles

Functional genomic screens with death rate analyses reveal mechanisms of drug action.

A common approach for understanding how drugs induce their therapeutic effects is to identify the genetic determinants of drug sensitivity. Because 'chemo-genetic profiles' are performed in a pooled format, inference of gene function is subject to several confounding influences related to variation in growth rates between clones. In this study, we developed Method for Evaluating Death Using a Simulation-assisted Approach (MEDUSA), which uses time-resolved measurements, along with model-driven constraints, to reveal the combination of growth and death rates that generated the observed drug response. MEDUSA is uniquely effective at identifying death regulatory genes. We apply MEDUSA to characterize DNA damage-induced lethality in the presence and absence of p53. Loss of p53 switches the mechanism of DNA damage-induced death from apoptosis to a non-apoptotic death that requires high respiration. These findings demonstrate the utility of MEDUSA both for determining the genetic dependencies of lethality and for revealing opportunities to potentiate chemo-efficacy in a cancer-specific manner.

Functionally-instructed modifiers of response to ATR inhibition in experimental glioma

BackgroundThe DNA damage response (DDR) is a physiological network preventing malignant transformation, e.g. by halting cell cycle progression upon DNA damage detection and promoting DNA repair. Glioblastoma are incurable primary tumors of the nervous system and DDR dysregulation contributes to acquired treatment resistance. Therefore, DDR targeting is a promising therapeutic anti-glioma strategy. Here, we investigated Ataxia telangiectasia and Rad3 related (ATR) inhibition (ATRi) and functionally-instructed combination therapies involving ATRi in experimental glioma.MethodsWe used acute cytotoxicity to identify treatment efficacy as well as RNAseq and DigiWest protein profiling to characterize ATRi-induced modulations within the molecular network in glioma cells. Genome-wide CRISPR/Cas9 functional genomic screens and subsequent validation with functionally-instructed compounds and selected shRNA-based silencing were employed to discover and investigate molecular targets modifying response to ATRi in glioma cell lines in vitro, in primary cultures ex vivo and in zebrafish and murine models in vivo.ResultsATRi monotherapy displays anti-glioma efficacy in vitro and ex vivo and modulates the molecular network. We discovered molecular targets by genome-wide CRISPR/Cas9 loss-of-function and activation screens that enhance therapeutic ATRi effects. We validated selected druggable targets by a customized drug library and functional assays in vitro, ex vivo and in vivo.ConclusionIn conclusion, our study leads to the identification of novel combination therapies involving ATRi that could inform future preclinical studies and early phase clinical trials.

Functional interactions between neurofibromatosis tumor suppressors underlie Schwann cell tumor de-differentiation and treatment resistance

Schwann cell tumors are the most common cancers of the peripheral nervous system and can arise in patients with neurofibromatosis type-1 (NF-1) or neurofibromatosis type-2 (NF-2). Functional interactions between NF1 and NF2 and broader mechanisms underlying malignant transformation of the Schwann lineage are unclear. Here we integrate bulk and single-cell genomics, biochemistry, and pharmacology across human samples, cell lines, and mouse allografts to identify cellular de-differentiation mechanisms driving malignant transformation and treatment resistance. We find DNA methylation groups of Schwann cell tumors can be distinguished by differentiation programs that correlate with response to the MEK inhibitor selumetinib. Functional genomic screening in NF1-mutant tumor cells reveals NF2 loss and PAK activation underlie selumetinib resistance, and we find that concurrent MEK and PAK inhibition is effective in vivo. These data support a de-differentiation paradigm underlying malignant transformation and treatment resistance of Schwann cell tumors and elucidate a functional link between NF1 and NF2.

Genome Editing, Transcriptional Regulation, and Forward Genetic Screening Using CRISPR-Cas12a Systems in Yarrowia lipolytica.

Class II Type V endonucleases have increasingly been adapted to develop sophisticated and easily accessible synthetic biology tools for genome editing, transcriptional regulation, and functional genomic screening in a wide range of organisms. One such endonuclease, Cas12a, presents itself as an attractive alternative to Cas9-based systems. The ability to mature its own guide RNAs (gRNAs) from a single transcript has been leveraged for easy multiplexing, and its lack of requirement of a tracrRNA element, also allows for short gRNA expression cassettes. To extend these functionalities into the industrially relevant oleaginous yeast Yarrowia lipolytica, we developed a set of CRISPR-Cas12a vectors for easy multiplexed gene knockout, repression, and activation. We further extended the utility of this CRISPR-Cas12a system to functional genomic screening by constructing a genome-wide guide library targeting every gene with an eightfold coverage. Pooled CRISPR screens conducted with this library were used to profile Cas12a guide activities and develop a machine learning algorithm that could accurately predict highly efficient Cas12a gRNA. In this protocols chapter, we first present a method by which protein coding genes may be functionally disrupted via indel formation with CRISPR-Cas12a systems. Further, we describe how Cas12a fused to a transcriptional regulator can be used in conjunction with shortened gRNA to achieve transcriptional repression or activation. Finally, we describe the design, cloning, and validation of a genome-wide library as well as a protocol for the execution of a pooled CRISPR screen, to determine guide activity profiles in a genome-wide context in Y. lipolytica. The tools and strategies discussed here expand the list of available synthetic biology tools for facile genome engineering in this industrially important host.

High-Throughput RNA-HCR-FISH Detection of Endogenous Pre-mRNA Splice Variants.

RNA-fluorescence in situ hybridization (RNA-FISH) is an essential and widely used tool for visualizing RNA molecules in intact cells. Recent advances have increased RNA-FISH sensitivity, signal detection efficiency, and throughput. However, detection of endogenous mRNA splice variants has been challenging due to the limits of visualization of RNA-FISH fluorescence signals and due to the limited number of RNA-FISH probes per target. HiFENS (high-throughput FISH detection of endogenous pre-mRNA splicing isoforms) is a method that enables visualization and relative quantification of mRNA splice variants at single-cell resolution in an automated high-throughput manner. HiFENS incorporates HCR (hybridization chain reaction) signal amplification strategies to enhance the fluorescence signal generated by low abundance transcripts or a small number of FISH probes targeting short stretches of RNA, such as single exons. The technique offers a significant advance in high-throughput FISH-based RNA detection and provides a powerful tool that can be used as a readout in functional genomics screens to discover and dissect cellular pathways regulating gene expression and alternative pre-mRNA splicing events.

Abstract A060: Discovery of potent and selective EP300 degraders with anti-cancer activity

Abstract E1A binding protein (EP300) and CREB binding protein (CBP) are paralog histone acetyltransferases involved in many cellular processes via their activity as transcriptional co-activators. Dysregulation of one or both proteins has been implicated in various cancers, and functional genomic screens have demonstrated a bidirectional synthetic lethal relationship between the two genes in tumor cells. Due to the high homology between EP300 and CBP, identifying chemical matter that selectively targets EP300 or CBP has proven challenging. Here, we describe a potent, highly selective heterobifunctional degrader of EP300 with biological activity in CBP-deficient and EP300-dependent tumor cells. This compound achieves robust degradation of EP300 in vivo at doses that sustain EP300 degradation with minimal effects on CBP protein level. Targeted degradation of EP300 protein resulted in a stronger suppression of cell growth and survival than targeting the bromodomain or HAT activity of EP300/CBP with small molecule inhibitors. Strong anti-proliferative effects have been demonstrated in multiple cancer types, including lymphomas and prostate tumors, highlighting the essential role of EP300 in the malignant phenotype of these tumors. Citation Format: Mark Zimmerman, Darshan Sappal, Ammar Adam, Benjamin Adams, Ketaki Adhikari, Wesley Austin, Breanna Bullock, Julie DiBernardo, Thomas Dixon, GiNell Elliot, Brian Ethell, Anais Gervais, Md Imran Hossain, David Huang, Dave Lahr, Mei Yun Lin, David Mayhew, Karolina Mizeracka, Solymar Negretti, Tyler Nguyen, Olga Prifti, Shawn Schiller, Brenna Sherbanee, David Terry, Nihan Ucisik, Elizabeth Wittenborn, Laura LaBonte, Danette Daniels. Discovery of potent and selective EP300 degraders with anti-cancer activity [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2023 Oct 11-15; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2023;22(12 Suppl):Abstract nr A060.

Abstract B143: HRO761, a first-in-class, clinical stage WRN inhibitor with potent preclinical anti-tumor activity in MSIhigh models

Abstract Colorectal carcinoma (CRC) as well as other indications that have lost competence in mismatch repair and have a phenotype known as microsatellite instability (MSI), remain a significant unmet medical need. Large-scale functional genomics screens across cell line panels have identified the Werner Syndrome RecQ helicase (WRN) as being synthetic lethal with MSIhigh cancers. Thus, WRN inhibitors may offer a new therapeutic avenue in MSIhigh cancers. We report the preclinical characteristics of HRO761, a first-in-class, highly potent and selective inhibitor of WRN helicase, with optimized drug-like properties. We evaluated the pharmacokinetics and pharmacodynamics of HRO761 when administered orally once daily in MSIhigh and also assessed the efficacy as single agent across cell (CDX) and patient-derived (PDX) xenografts as well as in combination with irinotecan, a prodrug derivative of camptothecin which acts as topoisomerase I inhibitor. Daily oral administration of HRO761 was highly efficacious across the CDX and PDX models evaluated with a disease control rate of 70% across all indications. In the SW48 CDX model, tumor regression was sustained for 2 months at the highest dose followed by slow relapse. HRO761 also induced WRN degradation, activated the DNA damage response and induced p53 target genes in a dose dependent manner. In vitro combination with irinotecan deepened sensitivity to HRO761 with a more sustained response in MSIhigh CRC cell lines. Combining administration of HRO761 daily with irinotecan weekly was highly efficacious and beneficial over single agent in CDX and PDX models. In CDX models, the combination showed a robust benefit, translating into homogeneous, sustained regression leading to a dose-dependent tumor growth delay even after treatment discontinuation. Overall, we demonstrate that the novel, selective, first-in-class WRN inhibitor HRO761 is a promising therapeutic approach for the treatment of MSIhigh cancer patients either as single agent or in combination with irinotecan. Clinical development is currently ongoing to assess the safety, tolerability and preliminary anti-tumor activity in patients with MSIhigh colorectal cancer and other MSIhigh solid tumors [clinicaltrials.gov NCT05838768]. Citation Format: Stephane Ferretti, Isabel Jaco, Andrea Decker, Christelle Hemmerlin, Clemens Scheufler, Cornelia Quadt, Dario Sterker, Elena Gavioli, Eloisa Jimenez Nunez, Ernesta Dammassa, Fanny Schaeffer, Genevieve Albrecht, Giorgia Clementi, Hansjoerg Martus, Jacques Hamon, Juergen Hinrichs, Laurent Laborde, Marion Dourdoigne, Michele Moschetta, Ramona Stump, Rita Andraos-Rey, Sarah Welly, Stephanie Barbe, Vincent Romanet, Markus Reschke, Ruben De Kanter, Michael Jensen, Henrik Moebitz, Marta Cortes Cros. HRO761, a first-in-class, clinical stage WRN inhibitor with potent preclinical anti-tumor activity in MSIhigh models [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2023 Oct 11-15; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2023;22(12 Suppl):Abstract nr B143.

Alzheimer’s disease risk genes modulate interferon‐responsive microglia

AbstractBackgroundHuman genetic studies and quantitative trait loci (eQTL) analyses have implicated microglia, the innate immune cells of the brain, as critical players in Alzheimer’s disease (AD). In response, the field has focused on defining microglial transcriptional states across many models of neurodegeneration. Integrating these studies, we find several microglial activation states commonly enriched in AD. These include Disease Associated Microglia, Antigen‐Presenting Microglia, and Interferon‐Responsive Microglia. To truly make use of this data, we now must move towards understanding the functional consequence of these microglial activation states and uncovering the mechanisms that drive changes in cell states during disease.MethodTo determine regulators of the interferon‐responsive microglial activation state, we performed a functional genomic screen using CRISPR interference (CRISPRi). To avoid bottlenecking of CRISPR guide RNAs and accommodate the large numbers of cells needed for screening, we implemented a novel transcription‐factor driven approach to differentiate iPSC‐derived microglia. After differentiation, microglia harboring individual gene knockdowns were primed with type I interferon and screened on expression of IFIT1, a key marker of the interferon‐responsive microglial state.ResultThese experiments have uncovered both canonical regulators of interferon signaling as well as novel regulators of the interferon‐responsive microglial state. Furthermore, several known AD risk genes were shown to inhibit microglial entry into this state, suggesting a potential mechanism whereby these nodes may influence disease risk.ConclusionRecent transcriptomic studies have uncovered significant heterogeneity of microglial states including the discovery of several disease associated states. However, it remains unclear whether these states represent beneficial, disease‐fighting functions or detrimental, disease‐promoting functions. Gaining the ability to selectively control entry into specific microglial states will allow novel investigations into how these specific microglial states interact with and guide disease trajectories.

CRISPR Activation Screen to Optimize Chimeric Antigen Receptor (CAR) T Cell Immunophenotype

Background Adoptive cellular therapies (ACT) offer potential cures for a wide range of human diseases. Immune cells genetically engineered with chimeric antigen receptors (CAR) have a clear role in B cell malignancies, and broader applications to other malignancies are under investigation. Significant obstacles to the long-term durability and potency of these therapies remain. As CAR T-cells remain persistently activated and exposed to the tumor microenvironment, they may become functionally impaired. This exhaustion can result in reduced production of cytokines, cytotoxic molecules, and other factors necessary for effective tumor cell killing. While CAR T-cell exhaustion remains a challenge, ongoing research and clinical trials aim to further understand the underlying mechanisms and develop strategies to overcome this limitation. Clinical studies have demonstrated that sustained in vivo persistence and function of CAR T-cells is correlated with infusion product phenotype. Persistence is associated with a less differentiated, memory-like phenotype. Naive (T N), Central Memory (T CM) and Stem Cell Memory (T SCM) are related to good proliferative potential and durability. CRISPR technologies offer powerful, functional genomic screens that permit the inference and modification of gene regulatory networks to optimize CAR T-cell persistence and function. Genome-wide CRISPR knock out screens have identified the CD2-CD58 axis as critical to T-cell potency and the knock-out of the mediator kinase module as beneficial to potency. Here we report on in vitro CRISPR activation screening of human CAR T-cells for gene expression that enhances T cell phenotypes associated with durable function. Method We designed 9,983 single guide RNA (sgRNA) to target the transcription start sites (TSS) of all known human transcription factors as well as a select set of genes with critical functions in T-cells. The vast majority of genes was targeted with 5 different sgRNA. To this we added a pool of 1,000 non-targeting sgRNA. We lentivirally transduced the sgRNA in tandem to a second-generation anti-CD20 CAR into primary T cells from healthy human donors. In a second lentiviral transduction we integrated catalytically dead Cas9 (dCas9) fused to the strong transcription activator, VP64 and in tandem to a puromycin resistance gene. After puromycin selection, we cultured the CRISPR activated CAR T-cells either alone or together with Raji or K562 cells for nine days before immunophenotypically sorting the CAR T-cells into CD8 + Naive, CD8 + Transitional Memory, CD8 + Central Memory, CD8 + Stem Cell Memory, CD8 + Transitional Effector, CD8 + Effector Memory, CD4 + Naive and CD4 + Memory T-cell populations. We extracted DNA from the sorted cells and prepared sequencing libraries targeting the sgRNA. We carried out two-way ANOVA analysis on the experimental conditions (trivial, Raji or K562 co-culture versus addition of CD2/CD28 activation beads) and T-cell subset to determine which genes targeted by CRISPRa resulted in favorable phenotypes. Result Across the different culture conditions, sgRNA activating IL23A and KLF4 were enriched in Transitional Memory versus Stem Cell Memory T-cells. Across the different CAR T-cell immunophenotypes, sgRNA targeting the TGFß (including KLF10), TNFα (including NFATC4) and FOXO signaling pathways were enriched in CAR T-cells cultured alone compared to CAR T-cells activated with CD2/CD28. Conclusion In this study we have investigated genes whose activation results in CAR T-cell immunophenotypes that are associated with more potent and durable CAR T-cell therapies. Our experiment suggests that the activation of genes, such as IL23 and KLF4, enriches for a Transitional Memory over Stem Cell Memory phenotype in a manner that is largely independent of environment (culture conditions). The activation of other genes, such as KLF10 and NFATC4, results in enrichment of CAR T-cells across all T-cell subsets, but in a manner that is dependent on environment. Gene expression requires coordinated control of gene-proximal and -distal cis-regulatory elements (CREs), whose state varies as cells differentiate or respond to environmental signals. Abrogating or adding CRE may result in more durable and potent adoptive cellular therapies independent of environment.

Selinexor Disrupts Epigenetic Programing and Modulates Immunogenicity in Multiple Myeloma

Introduction: Multiple myeloma (MM), remains an incurable, but treatable, malignancy that has benefitted from the integration of novel immunologic agents in clinical use, such as monoclonal antibody, bispecific T cell engager (TCE) and chimeric antigen receptor T (CAR-T) cell and antibody-drug conjugate based therapies. While these agents have shown great promise in treating relapsed/refractory MM patients, there is heterogeneity in response to therapy due to differences in patients' immune tumor microenvironments (iTME) and immunogenicity of the tumor cells. Previously (Sudalagunta et. al., ASH 2021), we used functional genomic screening of 53 MM samples to characterize the transcriptomic profile associated with sensitivity to the selective nuclear export inhibitor, selinexor (SELI), which was shown to be inversely correlated with the transcriptomic profile associated with sensitivity to an anti-CD38 monoclonal antibody in 106 MM samples. In this study, we investigate the cellular mechanisms of SELI resistance leading to an immunogenic cell state. Methods: Primary cells from an additional set of MM patients (n=65) were used for paired bulk RNAseq and an organotypic ex vivo drug sensitivity assay. Cellular viability was assessed using digital analysis of brightfield longitudinal images, and drug sensitivity was assessed as area under the dose-time-response curve (AUC) for each sample. We conducted gene set enrichment analysis (GSEA) on paired RNAseq data and ex vivo AUC of SELI. Individual genes and their respective correlations with ex vivo AUC of SELI were ranked from highest to lowest. We used RNAseq to examine the role of acute exposure to SELI on CD138-selected cells from paired pre/post treatment bone marrow specimens of five MM patients enrolled in the NCT02186834 clinical trial, taken prior to treatment with SELI and dexamethasone (DEX) and one/two weeks post-induction. Paired t-tests were carried out to quantify the changes. Results: In our previous analysis of 844 MM samples, a subset of the genes whose expression is positively correlated with SELI AUC are controlled by H3K27me3, which is catalyzed by nuclear accumulation of EZH2. Tripathi et. al. (AACR Annual Meeting 2022 & 2023) have demonstrated that XPO1 and EZH2 form a complex, and inhibition of XPO1 activity blocks export of EZH2. GSEA of a new cohort of 65 MM patient samples demonstrated increased ex vivo SELI sensitivity (lower AUC, hence negative correlation) in samples with higher expression of known epigenetic regulators of the polycomb repressive complex 2 (PRC2), such as EZH2, SUZ12, and EED; as shown in Figure a, which led to lowering of expression of genes bound by these proteins due to gene silencing via H3K27me3, shown in red in Figure a, according to GSEA of ChIPseq databases. Furthermore, gene expression profiles of paired pre/post SELI/DEX treatment samples showed reduced expression of genes associated with H3K27me3 histone modification (paired t-test p=2.06e-83, 1.1061e-81, 2.56e-80, 5.15e-12, 6.22e-4) in all five patients. These observations suggest that XPO1 inhibition leads to nuclear accumulation of epigenetic modulators like EZH2, SUZ12, and EED, which catalyze H3K27me3 activity leading to silencing of H3K27me3-regulated genes. GSEA of paired ex vivo AUC and RNAseq data for an EZH2 inhibitor (EPZ-6438) resulted in negative correlation between expression of genes in PRC2 and ex vivo sensitivity, similar to SELI. Interestingly, we observed that this epigenetic mechanism of resistance to SELI is inversely associated with key immunological targets in MM such as CD74, SLAMF7, CD38, CD46, TNFRSF13B, and LY9 as shown in Figure a. In Figure b, we present a graphical visualization of intracellular changes due to exposure to SELI in a naïve and a SELI-resistant MM cell. Conclusions: Here we propose that SELI disrupts EZH2 subcellular localization, dysregulating H3K27me3-gene suppression. Sequential BM biopsies from SELI treated patients demonstrate decreased expression of H3K27me3-regulated genes, suggesting increased EZH2 activity by nuclear accumulation. Finally, SELI-resistant primary cells have increased expression of EZH2 regulated immunotherapy target genes (e.g. CD38, SLAMF7), suggesting a possible therapeutic strategy of sequential treatment of SELI followed by immunotherapy. Using multiparameter flow cytometry of samples pre/post SELI/DEX induction, we will study SELI's effect on iTME.

A Functional Genomic Screen to Identify Novel Genes Involved in Multiple Myeloma Tumour Development

Multiple myeloma (MM) is an incurable malignancy characterised by uncontrolled proliferation of plasma cells (PCs) in the bone marrow (BM). MM is a genetically heterogeneous disease with each patient's PCs harbouring unique genetic mutations, however the development of MM tumours is not only dependent on the underlying genetics but also on the selective pressures applied by the BM microenvironment. Hence, we hypothesise that identifying the dependencies which promote MM cell outgrowth in the BM microenvironment will allow for the identification of new drug targets. This project aims to use an in vivo functional screen, combining CRISPR-Cas9 gene editing, using a single guide RNA (sgRNA) library targeting thousands of genes, with a murine model of MM, to identify novel genes involved in MM tumour development. The Bassik human apoptosis and cancer CRISPR knockout library (Addgene #101926) was used to transduce MM cells with 31,324 unique sgRNAs targeting 3,015 genes and 1,500 control regions. The human MM cell line OPM2 was transduced with the Cas9 transgene (FuCas9GFP), followed by the sgRNA expression vector (pMCB320) generating OPM2-Cas9-sgRNA cells. These cells were then injected (5x10 5 cells/mouse, n=8 mice + 1 tumour-naive control) into the tibiae of immunodeficient NOD.Cg-Prkdc scid Il2rg tm1Wjl/SzJ (NSG) mice. Four weeks post-injection, the primary tumour within the injected tibia was isolated and sgRNA frequencies were assessed via next generation sequencing and compared with that of the initial library and the injected cells using the MAGeCK algorithm. To identify genes that play roles in MM tumour formation in patients, we assessed the expression of our identified genes in CD138-selected BM PCs from newly diagnosed MM patients (n=155) when compared with normal PCs from healthy donors (n=5) in microarray dataset E-MTAB-363 (ArrayExpress). Furthermore, to identify potential pathway involvement of our identified genes we performed gene set enrichment analysis (GSEA) using the University of California San Diego and Broad Institute developed GSEA software, comparing to Hallmark Gene Sets. Comparison of the in vitro cultured cells with the initial library identified 34 genes with undetectable (p<0.05) sgRNAs post in vitro culture, suggesting these genes play key roles in the survival of OPM2 cells. This list contained a number of genes known to control MM cell survival, such as IRF4 and CCND2. Additionally, sgRNAs targeting a further 115 genes were significantly depleted (p<0.05) in vitro and were completely undetectable in vivo, suggesting that these are implicated in cell proliferation and tumour formation in vivo. This gene list contained some already known MM drivers such as MYC, KRAS and DIS3. Furthermore, we identified 28 genes where CRISPR knockout had no significant effect in vitro but were critical dependencies for tumour development, with sgRNAs being undetectable (p<0.05) in vivo. This list contained genes known to be important for MM and BM stromal cell interaction such as ALCAM and SPP1. Notably, 25% of the genes identified as critical dependencies for OPM2 cells in vitro and/or in vivo are significantly upregulated >2-fold (p<0.05; limma) in MM patient PCs when compared with normal PCs (E-MTAB-363), suggesting these may play a role in tumour development in patients. Gene set enrichment analysis revealed that the genes that were essential in vitro were predominantly involved in cell survival related pathways such as cell cycle regulation (G2M checkpoint) and DNA replication/repair; however, the uniquely in vivo genes are predominantly involved in cell metabolism pathways such as MTORC1 signalling, oxidative phosphorylation and glycolysis. Overall, this study has successfully used a functional genetic screen to identify novel gene dependencies required for in vitro MM cell growth and in vivo MM tumour development. A secondary genetic screen will be performed in vivo to validate and refine this list of candidate genes. These validated genes will be further investigated using patient derived data sets to determine if they represent common vulnerabilities in patients and in multiple MM cell lines to investigate their specific mechanistic and pathway involvement. Identification and understanding the driving dependencies which promote MM cell outgrowth in the BM microenvironment will allow for the identification and development of new drug targets to improve patient outcomes.

Thermococcus kodakarensis TK0353 is a novel AP lyase with a new fold

Hyperthermophilic organisms thrive in extreme environments prone to high levels of DNA damage. Growth at high temperature stimulates DNA base hydrolysis resulting in apurinic/apyrimidinic (AP) sites that destabilize the genome. Organisms across all domains have evolved enzymes to recognize and repair AP sites to maintain genome stability. The hyperthermophilic archaeon Thermococcus kodakarensis encodes several enzymes to repair AP site damage including the essential AP endonuclease TK Endonuclease IV. Recently, using functional genomic screening, we discovered a new family of AP lyases typified by TK0353. Here, using biochemistry, structural analysis and genetic deletion, we have characterized the TK0353 structure and function. TK0353 lacks glycosylase activity on a variety of damaged bases and is therefore either a monofunctional AP lyase or may be a glycosylase-lyase on a yet unidentified substrate. The crystal structure of TK0353 revealed a novel fold, which does not resemble other known DNA repair enzymes. The TK0353 gene is not essential for Thermococcus kodakarensis viability presumably due to redundant base excision repair enzymes involved in AP site processing. In summary, TK0353 is a novel AP lyase unique to hyperthermophiles that provides redundant repair activity necessary for genome maintenance.

MODL-47. A NEW MOUSE MODEL OF IDH MUTATED GLIOMAS IDENTIFIES TUMOR CELLS OF ORIGIN AND DETERMINANTS OF SENSITIVITY TO IDH INHIBITORS

Abstract Driver mutations in IDH1 and IDH2 characterize a substantial proportion of diffuse gliomas. These tumors are typically of lower grade at diagnosis, but many eventually transform to more aggressive disease. Furthermore, there is a lack of effective curative treatments. Mutated IDH molecules acquire neomorphic enzymatic activity, favoring the synthesis of D-2-hydroxyglutarate (D-2-HG). Accumulation of this metabolite perturbs many cellular functions and drives gliomagenesis through incompletely understood mechanisms. IDH mutations are early events in gliomagenesis, but the identity of the initiating cell type has been debated. Here, using genetically engineered mice, we found that combining Idh1R132H and Trp53 loss in oligodendrocyte progenitors leads to fully penetrant development of diffuse gliomas. The tumors recapitulate the cardinal features of the corresponding human disease. While heterogeneous, the mouse Idh1R132H and Idh1WT tumors show distinct patterns of transcriptional alterations and oncogenic copy number variants. Mutant IDH itself is an attractive drug target, as it is clonally expressed and usually retained during tumor evolution, even though its relevance for disease progression has been questioned. Nevertheless, recent clinical trials have demonstrated clear, but variable, anti-tumor effects of mutant IDH inhibitors in patients. To better understand the mechanistic basis for heterogeneous responses to mutant IDH inhibition, we performed CRISPR/Cas9 functional genomic screens in cell lines derived from the mouse Idh1R132H;Trp53MUT tumors. These revealed interactions between the mutant IDH inhibitor vorasidenib and molecules related to astrocytic differentiation, Notch signaling, and cell growth and metabolism. The translational relevance of these findings was supported by molecular data from human tumors, patients treated with IDH inhibitors, and human-derived cell models. Overall, these studies nominate oligodendrocyte progenitors as candidate cells of origin for IDH mutated gliomas, and point to strategies that may enhance the efficacy of IDH inhibitors.

EPCO-28. MULTIPLATFORM MOLECULAR AND FUNCTIONAL GENOMIC ANALYSES OF NF1 MUTANT GLIOBLASTOMA IDENTIFY DISTINCT CO-MUTATION PATTERNS AND DRUGGABLE DEPENDENCIES UNDERLYING SELUMETINIB RESPONSE

Abstract The tumor suppressor NF1 is recurrently mutated in sporadic IDH-wildtype glioblastoma and activates Ras/Raf/MEK signaling, motivating the use of targeted MEK inhibition. Here, we combine multiplatform molecular analysis of NF1 mutant glioblastoma human resection specimens with functional genomic CRISPRi screens to the MEK inhibitor selumetinib. METHODS: We identified 123 consecutive patients treated at the University of California San Francisco between 2016-2022 with NF1 mutant, IDH-wildtype glioblastoma. DNA methylation data was processed using the minfi package, and hierarchical clustering, principal component analysis (PCA), and survival analysis was performed in R. Human NF1 mutant GBM43 and mouse NRAS mutant SB28 GBM cells expressing CRISPRi machinery were used for genome wide selumetinib-resistance screens. Immunoblotting was performed to validate CRISPRi screen hits and selumetinib responses in cell lines. RESULTS: DNA sequencing of human glioblastomas identified 12 recurrent alterations in greater than 10% of cases including cell cycle regulators (CDKN2A/B, TP53, RB1), PI3K signaling (PTEN, PIK3R1, PIK3CA), epigenetic regulators (ATRX, SETD2), and Ras signaling (EGFR, PDGFRA, PTPN11). Methylation data (n=55 tumors) identified three groups marked by differences in MGMT methylation (p=0.04, chi-square test) and CDKN2A/B loss (p=0.01, chi-square test). CRISPRi screens in GBM43 and SB28 cells identified driver phenotypes in cell cycle regulators (TP53, CDKN2A, CDK4), PI3K signaling (PTEN, MTOR), and Ras signaling (RASA2, PTPN11) consistent with co-occurring mutations in the human NF1 mutant glioblastoma cohort. In contrast, mediators of selumetinib response converged on regulators of Ras/Raf/MEK signaling including the Ras effector SHOC2. Biochemical analysis validated SHOC2 repression is sufficient for selumetinib sensitization. CONCLUSIONS: Using a combination of human tumor genomic analysis, functional genetic screens, and biochemical validation, we identify three epigenetic groups of NF1 mutant glioblastoma with distinct co-mutation patterns and validate compensatory Ras/Raf/MEK activation as a conserved mechanism of selumetinib resistance, suggesting combination molecular therapy as a rational therapeutic combination.

EXTH-78. CRISPROFF EPIGENOME EDITING AND FUNCTIONAL GENOMICS DEFINE MEDIATORS OF CHEMOTHERAPY RESISTANCE IN GLIOBLASTOMA

Abstract Alkylating chemotherapies exhibit survival benefit for patients with glioblastoma (GBM). CRISPRoff is a programmable epigenetic memory writer that stably and heritably silences genes through DNA methylation. Epigenetic silencing of MGMT via promoter methylation predicts response to chemotherapy and is prognostic for overall survival. Here, we performed epigenome editing using CRISPRoff to stably silence MGMT through induced promoter methylation as a therapeutically tractable approach for potentiating GBM to chemotherapy. We then used genome-wide CRISPR interference (CRISPRi) screens to broadly define novel chemosensitizing targets of GBM cells. Epigenome editing was performed through electroporation of CRISPRoff mRNA and sgRNAs into MGMT hypomethylated GBM cell lines (LN18 and T98G) and primary GBM cells. Whole genome bisulfite sequencing, RT-qPCR, and western blot assessed gene silencing. Cell viability assays measured drug sensitivity. GBM xenografts with CRISPRoff-induced silencing of MGMT were intracranially transplanted in mice and treated with temozolomide. Genome-wide functional genomics screens were performed using CRISPR interference. CRISPRoff against MGMT reduced MGMT mRNA and protein levels by 99.97%, generating up to 150-fold sensitization to temozolomide and 9-fold sensitization to lomustine in GBM cells. Silencing remained stable and heritable in clonally isolated GBM cells at 145 days after transient delivery of CRISPRoff, and gene repression and drug sensitization were equivalent to Cas9-mediated homozygous deletion of MGMT in vitro and in vivo. CRISPRoff of MGMT in intracranial xenografts exhibited significant reduction in tumor growth with temozolomide treatment compared to temozolomide treatment with control CRISPRoff (p = 0.0047). CRISPRi screens validated MGMT and defined 266 and 238 additional modifiers of temozolomide and lomustine response, respectively, including the ATR (BRCA2), DNA repair (REV1), cell cycle (PSMD13), and Fanconi anemia pathways (FANCI, FANCD2). In summary, we integrate targeted epigenome editing with unbiased genome-wide approaches to establish a novel discovery and therapeutic platform against glioblastoma.

STEM-24. WDR5 REPRESENTS A THERAPEUTICALLY EXPLOITABLE TARGET FOR CANCER STEM CELLS IN GLIOBLASTOMA

Abstract Glioblastomas (GBMs) are heterogeneous, treatment-resistant tumors that are driven by populations of cancer stem cells (CSCs). In this study, we perform an epigenetic-focused functional genomics screen in GBM organoids and identify WDR5 as an essential epigenetic regulator in the SOX2-enriched, therapy resistant cancer stem cell niche. Despite their importance for tumor growth, few molecular mechanisms critical for CSC population maintenance have been exploited for therapeutic development. We developed a spatially resolved loss-of-function screen in GBM patient-derived organoids to identify essential epigenetic regulators in the SOX2-enriched, therapy resistant niche. Our niche-specific screens identified WDR5, an H3K4 histone methyltransferase responsible for activating specific gene expression, as indispensable for GBM CSC growth and survival. In GBM CSC models, WDR5 inhibitors blocked WRAD complex assembly and reduced H3K4 trimethylation and expression of genes involved in CSC-relevant oncogenic pathways. H3K4me3 peaks lost with WDR5 inhibitor treatment occurred disproportionally on POU transcription factor motifs, required for stem cell maintenance and including the POU5F1(OCT4)::SOX2 motif. We incorporated a SOX2/OCT4 motif driven GFP reporter system into our CSC cell models and found that WDR5 inhibitor treatment resulted in dose-dependent silencing of stem cell reporter activity. Further, WDR5 inhibitor treatment altered the stem cell state, disrupting CSC in vitro growth and self-renewal as well as in vivo tumor growth. Our results unveiled the role of WDR5 in maintaining the CSC state in GBM and provide a rationale for therapeutic development of WDR5 inhibitors for GBM and other advanced cancers. This conceptual and experimental framework can be applied to many cancers, and can unmask unique microenvironmental biology and rationally designed combination therapies.

Practical Characterization Strategies for Comparison, Qualification, and Selection of Cell Viability Detection Methods for Cellular Therapeutic Product Development and Manufacturing.

Cellular therapy development and manufacturing has focused on providing novel therapeutic cell-based products for various diseases. The International Organization for Standardization (ISO) has provided guidance on critical quality attributes (CQAs) that shall be considered when testing and releasing cellular therapeutic products. Cell count and viability measurements are two of the CQAs that are determined during development, manufacturing, testing, and product release. The ISO Cell Counting Standard Part 1 and 2 addressed the needs for improving the quality of cell counting results. However, there is currently no guidance on the qualification and selection of a fit-for-purpose cell viability detection method. In this work, we present strategies for the characterization and comparison of AO/PI and AO/DAPI staining methods using the heat-killed (HK) and low temperature/nutrient-deprived (LT/ND) cell death models to evaluate the comparability of cell viability measurements and identify potential causes of differences. We compared the AO/PI and AO/DAPI staining methods using HK and LT/ND-generated dead cells, investigated the staining time effects on cell viability measurements, and determined their viability linearity with different mixtures of live and dead cells. Furthermore, we validated AO/PI and AO/DAPI cell viability measurement with a long-term cell proliferation assay. Finally, we demonstrate a practical example of cell viability measurement comparison using AO/PI and AO/DAPI on antibiotic-selected transduced Jurkat and THP-1 cells to select a fit-for-purpose method for functional genomics screening. The proposed strategies may potentially enable scientists to properly characterize, compare, and select cell viability detection methods that are critical for cellular therapeutic product development and manufacturing.

Understanding cancer drug resistance with Sleeping Beauty functional genomic screens: Application to MAPK inhibition in cutaneous melanoma

Combined BRAF and MEK inhibition is an effective treatment for BRAF-mutant cutaneous melanoma. However, most patients progress on this treatment due to drug resistance. Here, we applied the Sleeping Beauty transposon system to understand how melanoma evades MAPK inhibition. We found that the specific drug resistance mechanisms differed across melanomas in our genetic screens of five cutaneous melanoma cell lines. While drivers that reactivated MAPK were highly conserved, many others were cell-line specific. One such driver, VAV1, activated a de-differentiated transcriptional program like that of hyperactive RAC1, RAC1P29S. To target this mechanism, we showed that an inhibitor of SRC, saracatinib, blunts the VAV1-induced transcriptional reprogramming. Overall, we highlighted the importance of accounting for melanoma heterogeneity in treating cutaneous melanoma with MAPK inhibitors. Moreover, we demonstrated the utility of the Sleeping Beauty transposon system in understanding cancer drug resistance.

Genetic Profiling and Genome-Scale Dropout Screening to Identify Therapeutic Targets in Mouse Models of Malignant Peripheral Nerve Sheath Tumor.

Malignant Peripheral Nerve Sheath Tumors (MPNSTs) are derived from Schwann cells or their precursors. In patients with the tumor susceptibility syndrome neurofibromatosis type 1 (NF1), MPNSTs are the most common malignancy and the leading cause of death. These rare and aggressive soft-tissue sarcomas offer a stark future, with 5-year disease-free survival rates of 34-60%. Treatment options for individuals with MPNSTs are disappointingly limited, with disfiguring surgery being the foremost treatment option. Many once-promising therapies such as tipifarnib, an inhibitor of Ras signaling, have failed clinically. Likewise, phase II clinical trials with erlotinib, which targets the epidermal growth factor (EFGR), and sorafenib, which targets the vascular endothelial growth factor receptor (VEGF), platelet-derived growth factor receptor (PDGF), and Raf, in combination with standard chemotherapy, have also failed to produce a response in patients. In recent years, functional genomic screening methods combined with genetic profiling of cancer cell lines have proven useful for identifying essential cytoplasmic signaling pathways and the development of target-specific therapies. In the case of rare tumor types, a variation of this approach known as cross-species comparative oncogenomics is increasingly being used to identify novel therapeutic targets. In cross-species comparative oncogenomics, genetic profiling and functional genomics are performed in genetically engineered mouse (GEM) models and the results are then validated in the rare human specimens and cell lines that are available. This paper describes how to identify candidate driver gene mutations in human and mouse MPNST cells using whole exome sequencing (WES). We then describe how to perform genome-scale shRNA screens to identify and compare critical signaling pathways in mouse and human MPNST cells and identify druggable targets in these pathways. These methodologies provide an effective approach to identifying new therapeutic targets in a variety of human cancer types.

Partial gene suppression improves identification of cancer vulnerabilities when CRISPR-Cas9 knockout is pan-lethal

BackgroundHundreds of functional genomic screens have been performed across a diverse set of cancer contexts, as part of efforts such as the Cancer Dependency Map, to identify gene dependencies—genes whose loss of function reduces cell viability or fitness. Recently, large-scale screening efforts have shifted from RNAi to CRISPR-Cas9, due to superior efficacy and specificity. However, many effective oncology drugs only partially inhibit their protein targets, leading us to question whether partial suppression of genes using RNAi could reveal cancer vulnerabilities that are missed by complete knockout using CRISPR-Cas9. Here, we compare CRISPR-Cas9 and RNAi dependency profiles of genes across approximately 400 matched cancer cell lines.ResultsWe find that CRISPR screens accurately identify more gene dependencies per cell line, but the majority of each cell line’s dependencies are part of a set of 1867 genes that are shared dependencies across the entire collection (pan-lethals). While RNAi knockdown of about 30% of these genes is also pan-lethal, approximately 50% have selective dependency patterns across cell lines, suggesting they could still be cancer vulnerabilities. The accuracy of the unique RNAi selectivity is supported by associations to multi-omics profiles, drug sensitivity, and other expected co-dependencies.ConclusionsIncorporating RNAi data for genes that are pan-lethal knockouts facilitates the discovery of a wider range of gene targets than could be detected using the CRISPR dataset alone. This can aid in the interpretation of contrasting results obtained from CRISPR and RNAi screens and reinforce the importance of partial gene suppression methods in building a cancer dependency map.