Germline variation contributes to false negatives in CRISPR-based experiments with varying burden across ancestries.

Reducing disparities is vital for equitable access to precision treatments in cancer. Socioenvironmental factors are a major driver of disparities, but differences in genetic variation likely also contribute. The impact of genetic ancestry on prioritization of cancer targets in drug discovery pipelines has not been systematically explored due to the absence of pre-clinical data at the appropriate scale. Here, we analyze data from 611 genome-scale CRISPR/Cas9 viability experiments in human cell line models to identify ancestry-associated genetic dependencies essential for cell survival. Surprisingly, we find that most putative associations between ancestry and dependency arise from artifacts related to germline variants. Our analysis suggests that for 1.2-2.5% of guides, germline variants in sgRNA targeting sequences reduce cutting by the CRISPR/Cas9 nuclease, disproportionately affecting cell models derived from individuals of recent African descent. We propose three approaches to mitigate this experimental bias, enabling the scientific community to address these disparities.

Rapid-kinetics degron benchmarking reveals off-target activities and mixed agonism-antagonism of MYB inhibitors.

Attenuating aberrant transcriptional circuits holds great promise for the treatment of numerous diseases, including cancer. However, development of transcriptional inhibitors is hampered by the lack of a generally accepted functional cellular readout to characterize their target specificity and on-target activity. We benchmarked the direct gene-regulatory signatures of six agents reported as inhibitors of the oncogenic transcription factor MYB against targeted MYB degradation in a nascent transcriptomics assay. The inhibitors demonstrated partial specificity for MYB target genes but displayed significant off-target activity. Unexpectedly, the inhibitors displayed bimodal on-target effects, acting as mixed agonists-antagonists. Our data uncover unforeseen agonist effects of small molecules originally developed as TF inhibitors and argue that rapid-kinetics benchmarking against degron models should be used for functional characterization of transcriptional modulators.

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.

Leukemia core transcriptional circuitry is a sparsely interconnected hierarchy stabilized by incoherent feed-forward loops.

Lineage-defining transcription factors form densely interconnected circuits in chromatin occupancy assays, but the functional significance of these networks remains underexplored. We reconstructed the functional topology of a leukemia cell transcription network from the direct gene-regulatory programs of eight core transcriptional regulators established in pre-steady state assays coupling targeted protein degradation with nascent transcriptomics. The core regulators displayed narrow, largely non-overlapping direct transcriptional programs, forming a sparsely interconnected functional hierarchy stabilized by incoherent feed-forward loops. BET bromodomain and CDK7 inhibitors disrupted the core regulators' direct programs, acting as mixed agonists/antagonists. The network is predictive of dynamic gene expression behaviors in time-resolved assays and clinically relevant pathway activity in patient populations.

Transcriptional Plasticity Drives Leukemia Immune Escape.

Relapse of acute myeloid leukemia (AML) after allogeneic bone marrow transplantation has been linked to immune evasion due to reduced expression of major histocompatibility complex class II (MHCII) genes through unknown mechanisms. In this work, we developed CORENODE, a computational algorithm for genome-wide transcription network decomposition that identified a transcription factor (TF) tetrad consisting of IRF8, MYB, MEF2C, and MEIS1, regulating MHCII expression in AML cells. We show that reduced MHCII expression at relapse is transcriptionally driven by combinatorial changes in the expression of these TFs, where MYB and IRF8 play major opposing roles, acting independently of the IFNγ/CIITA pathway. Beyond the MHCII genes, MYB and IRF8 antagonistically regulate a broad genetic program responsible for cytokine signaling and T-cell stimulation that displays reduced expression at relapse. A small number of cells with altered TF abundance and silenced MHCII expression are present at the time of initial leukemia diagnosis, likely contributing to eventual relapse. SIGNIFICANCE: Our findings point to an adaptive transcriptional mechanism of AML evolution after allogeneic transplantation whereby combinatorial fluctuations of TF expression under immune pressure result in the selection of cells with a silenced T-cell stimulation program. This article is highlighted in the In This Issue feature, p. 369.

Hypoxic, glycolytic metabolism is a vulnerability of B-acute lymphoblastic leukemia-initiating cells.

High-risk forms of B-acute lymphoblastic leukemia (B-ALL) remain a therapeutic challenge. Leukemia-initiating cells (LICs) self-renew and spark relapse and therefore have been the subject of intensive investigation; however, the properties of LICs in high-risk B-ALL are not well understood. Here, we use single-cell transcriptomics and quantitative xenotransplantation to understand LICs in MLL-rearranged (MLL-r) B-ALL. Compared with reported LIC frequencies in acute myeloid leukemia (AML), engraftable LICs in MLL-r B-ALL are abundant. Although we find that multipotent, self-renewing LICs are enriched among phenotypically undifferentiated B-ALL cells, LICs with the capacity to replenish the leukemic cellular diversity can emerge from more mature fractions. While inhibiting oxidative phosphorylation blunts blast proliferation, this intervention promotes LIC emergence. Conversely, inhibiting hypoxia and glycolysis impairs MLL-r B-ALL LICs, providing a therapeutic benefit in xenotransplantation systems. These findings provide insight into the aggressive nature of MLL-r B-ALL and provide a rationale for therapeutic targeting of hypoxia and glycolysis.

A distinct core regulatory module enforces oncogene expression in KMT2A-rearranged leukemia.

Acute myeloid leukemia with KMT2A (MLL) rearrangements is characterized by specific patterns of gene expression and enhancer architecture, implying unique core transcriptional regulatory circuitry. Here, we identified the transcription factors MEF2D and IRF8 as selective transcriptional dependencies of KMT2A-rearranged AML, where MEF2D displays partially redundant functions with its paralog, MEF2C. Rapid transcription factor degradation followed by measurements of genome-wide transcription rates and superresolution microscopy revealed that MEF2D and IRF8 form a distinct core regulatory module with a narrow direct transcriptional program that includes activation of the key oncogenes MYC, HOXA9, and BCL2. Our study illustrates a mechanism of context-specific transcriptional addiction whereby a specific AML subclass depends on a highly specialized core regulatory module to directly enforce expression of common leukemia oncogenes.

Synthetic Lethal Interaction between the ESCRT Paralog Enzymes VPS4A and VPS4B in Cancers Harboring Loss of Chromosome 18q or 16q.

Few therapies target the loss of tumor suppressor genes in cancer. We examine CRISPR-SpCas9 and RNA-interference loss-of-function screens to identify new therapeutic targets associated with genomic loss of tumor suppressor genes. The endosomal sorting complexes required for transport (ESCRT) ATPases VPS4A and VPS4B score as strong synthetic lethal dependencies. VPS4A is essential in cancers harboring loss of VPS4B adjacent to SMAD4 on chromosome 18q and VPS4B is required in tumors with co-deletion of VPS4A and CDH1 (E-cadherin) on chromosome 16q. We demonstrate that more than 30% of cancers selectively require VPS4A or VPS4B. VPS4A suppression in VPS4B-deficient cells selectively leads to ESCRT-III filament accumulation, cytokinesis defects, nuclear deformation, G2/M arrest, apoptosis, and potent tumor regression. CRISPR-SpCas9 screening and integrative genomic analysis reveal other ESCRT members, regulators of abscission, and interferon signaling as modifiers of VPS4A dependency. We describe a compendium of synthetic lethal vulnerabilities and nominate VPS4A and VPS4B as high-priority therapeutic targets for cancers with 18q or 16q loss.

Repeat expansions confer WRN dependence in microsatellite-unstable cancers.

The RecQ DNA helicase WRN is a synthetic lethal target for cancer cells with microsatellite instability (MSI), a form of genetic hypermutability that arises from impaired mismatch repair1-4. Depletion of WRN induces widespread DNA double-strand breaks in MSI cells, leading to cell cycle arrest and/or apoptosis. However, the mechanism by which WRN protects MSI-associated cancers from double-strand breaks remains unclear. Here we show that TA-dinucleotide repeats are highly unstable in MSI cells and undergo large-scale expansions, distinct from previously described insertion or deletion mutations of a few nucleotides5. Expanded TA repeats form non-B DNA secondary structures that stall replication forks, activate the ATR checkpoint kinase, and require unwinding by the WRN helicase. In the absence of WRN, the expanded TA-dinucleotide repeats are susceptible to cleavage by the MUS81 nuclease, leading to massive chromosome shattering. These findings identify a distinct biomarker that underlies the synthetic lethal dependence on WRN, and support the development of therapeutic agents that target WRN for MSI-associated cancers.