ABSTRACT
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.
Subject(s)
Leukemia, Myeloid, Acute , Myeloid-Lymphoid Leukemia Protein , Gene Rearrangement , Humans , Interferon Regulatory Factors/genetics , Interferon Regulatory Factors/metabolism , Leukemia, Myeloid, Acute/genetics , Myeloid-Lymphoid Leukemia Protein/genetics , Myeloid-Lymphoid Leukemia Protein/metabolism , Oncogenes/geneticsABSTRACT
BACKGROUND: Divarasib (GDC-6036) is a covalent KRAS G12C inhibitor that was designed to have high potency and selectivity. METHODS: In a phase 1 study, we evaluated divarasib administered orally once daily (at doses ranging from 50 to 400 mg) in patients who had advanced or metastatic solid tumors that harbor a KRAS G12C mutation. The primary objective was an assessment of safety; pharmacokinetics, investigator-evaluated antitumor activity, and biomarkers of response and resistance were also assessed. RESULTS: A total of 137 patients (60 with non-small-cell lung cancer [NSCLC], 55 with colorectal cancer, and 22 with other solid tumors) received divarasib. No dose-limiting toxic effects or treatment-related deaths were reported. Treatment-related adverse events occurred in 127 patients (93%); grade 3 events occurred in 15 patients (11%) and a grade 4 event in 1 patient (1%). Treatment-related adverse events resulted in a dose reduction in 19 patients (14%) and discontinuation of treatment in 4 patients (3%). Among patients with NSCLC, a confirmed response was observed in 53.4% of patients (95% confidence interval [CI], 39.9 to 66.7), and the median progression-free survival was 13.1 months (95% CI, 8.8 to could not be estimated). Among patients with colorectal cancer, a confirmed response was observed in 29.1% of patients (95% CI, 17.6 to 42.9), and the median progression-free survival was 5.6 months (95% CI, 4.1 to 8.2). Responses were also observed in patients with other solid tumors. Serial assessment of circulating tumor DNA showed declines in KRAS G12C variant allele frequency associated with response and identified genomic alterations that may confer resistance to divarasib. CONCLUSIONS: Treatment with divarasib resulted in durable clinical responses across KRAS G12C-positive tumors, with mostly low-grade adverse events. (Funded by Genentech; ClinicalTrials.gov number, NCT04449874.).
Subject(s)
Antineoplastic Agents , Carcinoma, Non-Small-Cell Lung , Colorectal Neoplasms , Enzyme Inhibitors , Lung Neoplasms , Humans , Antineoplastic Agents/administration & dosage , Antineoplastic Agents/adverse effects , Antineoplastic Agents/therapeutic use , Carcinoma, Non-Small-Cell Lung/drug therapy , Carcinoma, Non-Small-Cell Lung/genetics , Colorectal Neoplasms/drug therapy , Colorectal Neoplasms/genetics , Lung Neoplasms/drug therapy , Lung Neoplasms/genetics , Mutation , Proto-Oncogene Proteins p21(ras)/genetics , Administration, Oral , Enzyme Inhibitors/administration & dosage , Enzyme Inhibitors/adverse effects , Enzyme Inhibitors/therapeutic useABSTRACT
Malaria parasites elude eradication attempts both within the human host and across nations. At the individual level, parasites evade the host immune responses through antigenic variation. At the global level, parasites escape drug pressure through single nucleotide variants and gene copy amplification events conferring drug resistance. Despite their importance to global health, the rates at which these genomic alterations emerge have not been determined. We studied the complete genomes of different Plasmodium falciparum clones that had been propagated asexually over one year in the presence and absence of drug pressure. A combination of whole-genome microarray analysis and next-generation deep resequencing (totaling 14 terabases) revealed a stable core genome with only 38 novel single nucleotide variants appearing in seventeen evolved clones (avg. 5.4 per clone). In clones exposed to atovaquone, we found cytochrome b mutations as well as an amplification event encompassing the P. falciparum multidrug resistance associated protein (mrp1) on chromosome 1. We observed 18 large-scale (>1 kb on average) deletions of telomere-proximal regions encoding multigene families, involved in immune evasion (9.5×10(-6) structural variants per base pair per generation). Six of these deletions were associated with chromosomal crossovers generated during mitosis. We found only minor differences in rates between genetically distinct strains and between parasites cultured in the presence or absence of drug. Using these derived mutation rates for P. falciparum (1.0-9.7×10(-9) mutations per base pair per generation), we can now model the frequency at which drug or immune resistance alleles will emerge under a well-defined set of assumptions. Further, the detection of mitotic recombination events in var gene families illustrates how multigene families can arise and change over time in P. falciparum. These results will help improve our understanding of how P. falciparum evolves to evade control efforts within both the individual hosts and large populations.
Subject(s)
Antigens , Atovaquone/administration & dosage , Drug Resistance, Multiple , Host-Parasite Interactions , Plasmodium falciparum , Antigenic Variation/drug effects , Antigenic Variation/genetics , Antigens/drug effects , Antigens/genetics , Cytochromes b/genetics , Drug Resistance, Multiple/drug effects , Drug Resistance, Multiple/genetics , Evolution, Molecular , Genome, Protozoan/drug effects , High-Throughput Nucleotide Sequencing , Host-Parasite Interactions/genetics , Host-Parasite Interactions/immunology , Humans , Malaria, Falciparum/genetics , Malaria, Falciparum/immunology , Mitosis/genetics , Multidrug Resistance-Associated Proteins/genetics , Multidrug Resistance-Associated Proteins/immunology , Multidrug Resistance-Associated Proteins/metabolism , Plasmodium falciparum/drug effects , Plasmodium falciparum/genetics , Plasmodium falciparum/immunologyABSTRACT
Intraerythrocytic malaria parasites can obtain nearly their entire amino acid requirement by degrading host cell hemoglobin. The sole exception is isoleucine, which is not present in adult human hemoglobin and must be obtained exogenously. We evaluated two compounds for their potential to interfere with isoleucine utilization. Mupirocin, a clinically used antibacterial, kills Plasmodium falciparum parasites at nanomolar concentrations. Thiaisoleucine, an isoleucine analog, also has antimalarial activity. To identify targets of the two compounds, we selected parasites resistant to either mupirocin or thiaisoleucine. Mutants were analyzed by genome-wide high-density tiling microarrays, DNA sequencing, and copy number variation analysis. The genomes of three independent mupirocin-resistant parasite clones had all acquired either amplifications encompassing or SNPs within the chromosomally encoded organellar (apicoplast) isoleucyl-tRNA synthetase. Thiaisoleucine-resistant parasites had a mutation in the cytoplasmic isoleucyl-tRNA synthetase. The role of this mutation in thiaisoleucine resistance was confirmed by allelic replacement. This approach is generally useful for elucidation of new targets in P. falciparum. Our study shows that isoleucine utilization is an essential pathway that can be targeted for antimalarial drug development.
Subject(s)
Isoleucine-tRNA Ligase/metabolism , Isoleucine/metabolism , Plasmodium falciparum/metabolism , Protozoan Proteins/metabolism , Adult , Animals , Anti-Bacterial Agents/pharmacology , Dose-Response Relationship, Drug , Drug Resistance/genetics , Erythrocytes/parasitology , Genome, Protozoan/genetics , Green Fluorescent Proteins/genetics , Green Fluorescent Proteins/metabolism , Hemoglobins/metabolism , Humans , Isoleucine/analogs & derivatives , Isoleucine/pharmacology , Isoleucine-tRNA Ligase/genetics , Microscopy, Fluorescence , Mupirocin/pharmacology , Mutation , Plasmodium falciparum/drug effects , Plasmodium falciparum/genetics , Polymorphism, Single Nucleotide , Protozoan Proteins/geneticsABSTRACT
Tissue-specific differences in the expression of paralog genes, which are not essential in most cell types due to the buffering effect of the partner pair, can make for highly selective gene dependencies. To identify selective paralogous targets for acute myeloid leukemia (AML), we integrated the Cancer Dependency Map with numerous datasets characterizing protein-protein interactions, paralog relationships, and gene expression in cancer models. Here, we identified ATP1B3 as a context-specific, paralog-related dependency in AML. ATP1B3, the beta subunit of the sodium-potassium pump (Na/K-ATP pump), interacts with the alpha subunit ATP1A1 to form an essential complex for maintaining cellular homeostasis and membrane potential in all eukaryotic cells. When ATP1B3's paralog ATP1B1 is poorly expressed, elimination of ATP1B3 leads to the destabilization of the Na/K-ATP pump. ATP1B1 expression is regulated through epigenetic silencing in hematopoietic lineage cells both through histone and DNA methylation in the promoter region. Loss of ATP1B3 in AML cells induced cell death in vitro and reduced leukemia burden in vivo, which could be rescued by stabilizing ATP1A1 through overexpression of ATP1B1. ATP1B3 is thus a potential therapeutic target for AML and other hematologic malignancies with low expression of ATP1B1.
ABSTRACT
BACKGROUND: Paediatric hepatocellular carcinomas (HCC) traditionally arise in the context of a normal structural and functional liver and carry a dismal prognosis. While chemotherapy is the frontline standard, there is emerging interest in the study of immunotherapies for paediatric patients with relapsed/refractory disease. There is limited data to support whether immunotherapies will be of utility in this patient population. METHODS: Six paediatric patients (median age:16 years, range: 12-17 at the time of treatment) with advanced hepatocellular neosplams, either conventional hepatocellular or fibrolamellar carcinoma, were treated with immunotherapy. Patients were consented to institutional genomic profiling and biobanking protocols. Baseline samples and serial tissue samples, when available, were evaluated for somatic mutation rate, actionable gene mutations, and pan-immune bulk RNA expression profiling. Results were correlated with clinical course. FINDINGS: Three patients responded to checkpoint inhibition: one achieved a complete, durable response and the other two, prolonged stable disease. Three additional patients progressed. Diagnostic tissue from the complete responder demonstrated a higher relative mutational burden and robust immune infiltrate. Pre-treatment samples from the three responders demonstrated decreased expression of genes associated with T-cell dysfunction. INTERPRETATION: A subset of patients with primary paediatric hepatocellular tumours will respond to immunotherapy. Immunotherapies are currently under prospective study for relapsed/refractory liver tumours in paediatric patients. Results from this report support the prospective collection of serial serum and tissue samples which may further identify genomic and immunophenotypic patterns predictive of response. FUNDING: This work was supported by Philanthropic funds (Pan Mass Challenge, Team Angus and Team Perspective).
Subject(s)
Carcinoma, Hepatocellular , Immunophenotyping , Immunotherapy , Liver Neoplasms , Humans , Liver Neoplasms/therapy , Liver Neoplasms/immunology , Liver Neoplasms/pathology , Carcinoma, Hepatocellular/therapy , Carcinoma, Hepatocellular/immunology , Carcinoma, Hepatocellular/pathology , Carcinoma, Hepatocellular/genetics , Male , Female , Child , Adolescent , Immunotherapy/methods , Mutation , Treatment Outcome , Biomarkers, Tumor , Gene Expression ProfilingABSTRACT
PURPOSE: To inform prognosis, treatment response, disease biology, and KRAS G12C mutation heterogeneity, we conducted exploratory circulating tumor DNA (ctDNA) profiling on 134 patients with solid tumors harboring a KRAS G12C mutation treated with single-agent divarasib (GDC-6036) in a phase 1 study. EXPERIMENTAL DESIGN: Plasma samples were collected for serial ctDNA profiling at baseline (cycle 1 day 1 prior to treatment) and multiple on-treatment time points (cycle 1 day 15 and cycle 3 day 1). RESULTS: KRAS G12C ctDNA was detectable from plasma samples in 72.9% (43/59) and 92.6% (50/54) of patients with non-small cell lung cancer and colorectal cancer, respectively, the majority of whom were eligible for study participation based on a local test detecting the KRAS G12C mutation in tumor tissue. Baseline ctDNA tumor fraction was associated with tumor type, disease burden, and metastatic sites. A decline in ctDNA level was observed as early as cycle 1 day 15. Serial assessment showed a decline in ctDNA tumor fraction associated with response and progression-free survival. Except for a few cases of KRAS G12C sub-clonality, on-treatment changes in KRAS G12C variant allele frequency mirrored changes in the overall ctDNA tumor fraction. CONCLUSIONS: Across tumor types, the KRAS G12C mutation likely represents a truncal mutation in the majority of patients. Rapid and deep decline in ctDNA tumor fraction was observed in patients responding to divarasib treatment. Early on-treatment dynamics of ctDNA were associated with patient outcomes and tumor response to divarasib treatment.
Subject(s)
Biomarkers, Tumor , Circulating Tumor DNA , Mutation , Proto-Oncogene Proteins p21(ras) , Humans , Circulating Tumor DNA/genetics , Circulating Tumor DNA/blood , Proto-Oncogene Proteins p21(ras)/genetics , Female , Male , Middle Aged , Aged , Biomarkers, Tumor/genetics , Neoplasms/genetics , Neoplasms/drug therapy , Neoplasms/pathology , Neoplasms/blood , Prognosis , Adult , Genetic Heterogeneity , Treatment Outcome , Aged, 80 and over , Carcinoma, Non-Small-Cell Lung/genetics , Carcinoma, Non-Small-Cell Lung/drug therapy , Carcinoma, Non-Small-Cell Lung/pathology , Carcinoma, Non-Small-Cell Lung/bloodABSTRACT
Anaplastic large cell lymphoma (ALCL) is an aggressive, CD30+ T cell lymphoma of children and adults. ALK fusion transcripts or mutations in the JAK-STAT pathway are observed in most ALCL tumors, but the mechanisms underlying tumorigenesis are not fully understood. Here, we show that dysregulated STAT3 in ALCL cooccupies enhancers with master transcription factors BATF3, IRF4, and IKZF1 to form a core regulatory circuit that establishes and maintains the malignant cell state in ALCL. Critical downstream targets of this network in ALCL cells include the protooncogene MYC, which requires active STAT3 to facilitate high levels of MYC transcription. The core autoregulatory transcriptional circuitry activity is reinforced by MYC binding to the enhancer regions associated with STAT3 and each of the core regulatory transcription factors. Thus, activation of STAT3 provides the crucial link between aberrant tyrosine kinase signaling and the core transcriptional machinery that drives tumorigenesis and creates therapeutic vulnerabilities in ALCL.
Subject(s)
Lymphoma, Large-Cell, Anaplastic , Signal Transduction , Adult , Child , Humans , Signal Transduction/genetics , Anaplastic Lymphoma Kinase/genetics , Anaplastic Lymphoma Kinase/metabolism , Lymphoma, Large-Cell, Anaplastic/genetics , Lymphoma, Large-Cell, Anaplastic/metabolism , Lymphoma, Large-Cell, Anaplastic/pathology , Janus Kinases/metabolism , STAT Transcription Factors/metabolism , Cell Transformation, Neoplastic , Carcinogenesis/genetics , STAT3 Transcription Factor/geneticsABSTRACT
KRAS G12C mutation is prevalent in ~4% of colorectal cancer (CRC) and is associated with poor prognosis. Divarasib, a KRAS G12C inhibitor, has shown modest activity as a single agent in KRAS G12C-positive CRC at 400 mg. Epidermal growth factor receptor has been recognized as a major upstream activator of RAS-MAPK signaling, a proposed key mechanism of resistance to KRAS G12C inhibition in CRC. Here, we report on divarasib plus cetuximab (epidermal growth factor receptor inhibitor) in patients with KRAS G12C-positive CRC (n = 29) from arm C of an ongoing phase 1b trial. The primary objective was to evaluate safety. Secondary objectives included preliminary antitumor activity. The safety profile of this combination was consistent with those of single-agent divarasib and cetuximab. Treatment-related adverse events led to divarasib dose reductions in four patients (13.8%); there were no treatment withdrawals. The objective response rate was 62.5% (95% confidence interval: 40.6%, 81.2%) in KRAS G12C inhibitor-naive patients (n = 24). The median duration of response was 6.9 months. The median progression-free survival was 8.1 months (95% confidence interval: 5.5, 12.3). As an exploratory objective, we observed a decline in KRAS G12C variant allele frequency associated with response and identified acquired genomic alterations at disease progression that may be associated with resistance. The manageable safety profile and encouraging antitumor activity of divarasib plus cetuximab support the further investigation of this combination in KRAS G12C-positive CRC.ClinicalTrials.gov identifier: NCT04449874.
Subject(s)
Colorectal Neoplasms , Proto-Oncogene Proteins p21(ras) , Humans , Cetuximab/adverse effects , Cetuximab/genetics , Proto-Oncogene Proteins p21(ras)/genetics , Colorectal Neoplasms/drug therapy , Colorectal Neoplasms/genetics , Colorectal Neoplasms/pathology , ErbB Receptors/genetics , Progression-Free Survival , Mutation/geneticsABSTRACT
Here, we fully characterize the genomes of 14 Plasmodium falciparum patient isolates taken recently from the Iquitos region using genome scanning, a microarray-based technique that delineates the majority of single-base changes, indels, and copy number variants distinguishing the coding regions of two clones. We show that the parasite population in the Peruvian Amazon bears a limited number of genotypes and low recombination frequencies. Despite the essentially clonal nature of some isolates, we see high frequencies of mutations in subtelomeric highly variable genes and internal var genes, indicating mutations arising during self-mating or mitotic replication. The data also reveal that one or two meioses separate different isolates, showing that P. falciparum clones isolated from different individuals in defined geographical regions could be useful in linkage analyses or quantitative trait locus studies. Through pairwise comparisons of different isolates we discovered point mutations in the apicoplast genome that are close to known mutations that confer clindamycin resistance in other species, but which were hitherto unknown in malaria parasites. Subsequent drug sensitivity testing revealed over 100-fold increase of clindamycin EC(50) in strains harboring one of these mutations. This evidence of clindamycin-resistant parasites in the Amazon suggests that a shift should be made in health policy away from quinine + clindamycin therapy for malaria in pregnant women and infants, and that the development of new lincosamide antibiotics for malaria should be reconsidered.
Subject(s)
Chromosomal Instability , Chromosome Mapping , Clindamycin , Drug Resistance/genetics , Malaria, Falciparum/parasitology , Plasmodium falciparum/genetics , Antimalarials/therapeutic use , Base Sequence , Chromosomal Instability/genetics , Chromosome Mapping/methods , Clindamycin/therapeutic use , DNA Copy Number Variations , Female , Gene Frequency , Genome, Protozoan , Genotype , Humans , Infant , Malaria, Falciparum/diagnosis , Malaria, Falciparum/drug therapy , Malaria, Falciparum/genetics , Male , Models, Biological , Models, Molecular , Molecular Sequence Data , Pedigree , Peru , Pregnancy , Telomere/geneticsABSTRACT
Plasmodium vivax causes 25-40% of malaria cases worldwide, yet research on this human malaria parasite has been neglected. Nevertheless, the recent publication of the P. vivax reference genome now allows genomics and systems biology approaches to be applied to this pathogen. We show here that whole-genome analysis of the parasite can be achieved directly from ex vivo-isolated parasites, without the need for in vitro propagation. A single isolate of P. vivax obtained from a febrile patient with clinical malaria from Peru was subjected to whole-genome sequencing (30× coverage). This analysis revealed over 18,261 single-nucleotide polymorphisms (SNPs), 6,257 of which were further validated using a tiling microarray. Within core chromosomal genes we find that one SNP per every 985 bases of coding sequence distinguishes this recent Peruvian isolate, designated IQ07, from the reference Salvador I strain obtained in 1972. This full-genome sequence of an uncultured P. vivax isolate shows that the same regions with low numbers of aligned sequencing reads are also highly variable by genomic microarray analysis. Finally, we show that the genes containing the largest ratio of nonsynonymous-to-synonymous SNPs include two AP2 transcription factors and the P. vivax multidrug resistance-associated protein (PvMRP1), an ABC transporter shown to be associated with quinoline and antifolate tolerance in Plasmodium falciparum. This analysis provides a data set for comparative analysis with important potential for identifying markers for global parasite diversity and drug resistance mapping studies.
Subject(s)
Drug Resistance/genetics , Genes, Protozoan/genetics , Oligonucleotide Array Sequence Analysis/methods , Plasmodium vivax/genetics , Selection, Genetic , Sequence Analysis, DNA/methods , Erythrocytes/parasitology , Gene Expression Regulation , Humans , Leukocytes/parasitology , Malaria Vaccines/immunology , Multigene Family/genetics , Mutation/genetics , Peru , Plasmodium vivax/immunology , Plasmodium vivax/isolation & purification , Polymorphism, Genetic , Sequence Alignment , Transcription Factors/geneticsABSTRACT
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.
ABSTRACT
Transcription factors (TFs) are frequently mutated in cancer. Paediatric cancers exhibit few mutations genome-wide but frequently harbour sentinel mutations that affect TFs, which provides a context to precisely study the transcriptional circuits that support mutant TF-driven oncogenesis. A broadly relevant mechanism that has garnered intense focus involves the ability of mutant TFs to hijack wild-type lineage-specific TFs in self-reinforcing transcriptional circuits. However, it is not known whether this specific type of circuitry is equally crucial in all mutant TF-driven cancers. Here we describe an alternative yet central transcriptional mechanism that promotes Ewing sarcoma, wherein constraint, rather than reinforcement, of the activity of the fusion TF EWS-FLI supports cancer growth. We discover that ETV6 is a crucial TF dependency that is specific to this disease because it, counter-intuitively, represses the transcriptional output of EWS-FLI. This work discovers a previously undescribed transcriptional mechanism that promotes cancer.
Subject(s)
Sarcoma, Ewing , Child , Humans , Cell Line, Tumor , Cell Transformation, Neoplastic/genetics , Gene Expression Regulation, Neoplastic , Oncogene Proteins, Fusion/genetics , Oncogene Proteins, Fusion/metabolism , Proto-Oncogene Protein c-fli-1/genetics , Proto-Oncogene Protein c-fli-1/metabolism , Proto-Oncogene Proteins c-ets/genetics , RNA-Binding Protein EWS/genetics , RNA-Binding Protein EWS/metabolism , Sarcoma, Ewing/geneticsABSTRACT
Clinical progress in multiple myeloma (MM), an incurable plasma cell (PC) neoplasia, has been driven by therapies that have limited applications beyond MM/PC neoplasias and do not target specific oncogenic mutations in MM. Instead, these agents target pathways critical for PC biology yet largely dispensable for malignant or normal cells of most other lineages. Here we systematically characterized the lineage-preferential molecular dependencies of MM through genome-scale clustered regularly interspaced short palindromic repeats (CRISPR) studies in 19 MM versus hundreds of non-MM lines and identified 116 genes whose disruption more significantly affects MM cell fitness compared with other malignancies. These genes, some known, others not previously linked to MM, encode transcription factors, chromatin modifiers, endoplasmic reticulum components, metabolic regulators or signaling molecules. Most of these genes are not among the top amplified, overexpressed or mutated in MM. Functional genomics approaches thus define new therapeutic targets in MM not readily identifiable by standard genomic, transcriptional or epigenetic profiling analyses.
Subject(s)
Multiple Myeloma , Humans , Multiple Myeloma/genetics , Genomics , Genome , Clustered Regularly Interspaced Short Palindromic Repeats/geneticsABSTRACT
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.
Subject(s)
Leukemia, Myeloid, Acute , Histocompatibility Antigens Class II/genetics , Humans , Interferon Regulatory Factors , Leukemia, Myeloid, Acute/genetics , Recurrence , Transplantation, HomologousABSTRACT
Leukemic blasts are immune cells gone awry. We hypothesized that dysregulation of inflammatory pathways contributes to the maintenance of their leukemic state and can be exploited as cell-intrinsic, self-directed immunotherapy. To this end, we applied genome-wide screens to discover genetic vulnerabilities in acute myeloid leukemia (AML) cells implicated in inflammatory pathways. We identified the immune modulator IRF2BP2 as a selective AML dependency. We validated AML cell dependency on IRF2BP2 with genetic and protein degradation approaches in vitro and genetically in vivo. Chromatin and global gene-expression studies demonstrated that IRF2BP2 represses IL1ß/TNFα signaling via NFκB, and IRF2BP2 perturbation results in an acute inflammatory state leading to AML cell death. These findings elucidate a hitherto unexplored AML dependency, reveal cell-intrinsic inflammatory signaling as a mechanism priming leukemic blasts for regulated cell death, and establish IRF2BP2-mediated transcriptional repression as a mechanism for blast survival. SIGNIFICANCE: This study exploits inflammatory programs inherent to AML blasts to identify genetic vulnerabilities in this disease. In doing so, we determined that AML cells are dependent on the transcriptional repressive activity of IRF2BP2 for their survival, revealing cell-intrinsic inflammation as a mechanism priming leukemic blasts for regulated cell death. See related commentary by Puissant and Medyouf, p. 1617. This article is highlighted in the In This Issue feature, p. 1599.
Subject(s)
Leukemia, Myeloid, Acute , Humans , Inflammation/genetics , Leukemia, Myeloid, Acute/genetics , NF-kappa B/metabolism , Signal TransductionABSTRACT
CRISPR-Cas9-based genetic screens have successfully identified cell type-dependent liabilities in cancer, including acute myeloid leukemia (AML), a devastating hematologic malignancy with poor overall survival. Because most of these screens have been performed in vitro using established cell lines, evaluating the physiologic relevance of these targets is critical. We have established a CRISPR screening approach using orthotopic xenograft models to validate and prioritize AML-enriched dependencies in vivo, including in CRISPR-competent AML patient-derived xenograft (PDX) models tractable for genome editing. Our integrated pipeline has revealed several targets with translational value, including SLC5A3 as a metabolic vulnerability for AML addicted to exogenous myo-inositol and MARCH5 as a critical guardian to prevent apoptosis in AML. MARCH5 repression enhanced the efficacy of BCL2 inhibitors such as venetoclax, further highlighting the clinical potential of targeting MARCH5 in AML. Our study provides a valuable strategy for discovery and prioritization of new candidate AML therapeutic targets. SIGNIFICANCE: There is an unmet need to improve the clinical outcome of AML. We developed an integrated in vivo screening approach to prioritize and validate AML dependencies with high translational potential. We identified SLC5A3 as a metabolic vulnerability and MARCH5 as a critical apoptosis regulator in AML, both of which represent novel therapeutic opportunities.This article is highlighted in the In This Issue feature, p. 275.
Subject(s)
Antineoplastic Agents/therapeutic use , CRISPR-Cas Systems , Leukemia, Myeloid, Acute/drug therapy , Precision Medicine , Xenograft Model Antitumor Assays , Animals , Humans , Leukemia, Myeloid, Acute/geneticsABSTRACT
Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2ß, a transcription factor member of the lineage-defining transcriptional core regulatory circuitry (CRC) in NB. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed "JQAD1" that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at CRC enhancers and rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. Furthermore, JQAD1 activity is critically determined by cereblon (CRBN) expression across NB cells. SIGNIFICANCE: EP300, but not CBP, controls oncogenic CRC-driven transcription in high-risk NB by binding TFAP2ß. We developed JQAD1, a CRBN-dependent PROTAC degrader with preferential activity against EP300 and demonstrated its activity in NB. JQAD1 has limited toxicity to untransformed cells and is effective in vivo in a CRBN-dependent manner. This article is highlighted in the In This Issue feature, p. 587.
Subject(s)
Neuroblastoma , Regulatory Sequences, Nucleic Acid , Acetylation , Child , E1A-Associated p300 Protein/genetics , Humans , N-Myc Proto-Oncogene Protein/genetics , Neuroblastoma/drug therapy , Neuroblastoma/genetics , OncogenesABSTRACT
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.
Subject(s)
Epigenome , Glioma , Animals , Humans , Mutation , Glioma/genetics , Transcription Factors/genetics , Transcription Factors/metabolism , Neoplastic Stem Cells/metabolism , Mammals/genetics , Mammals/metabolism , DNA Helicases/genetics , Nuclear Proteins/geneticsABSTRACT
The combination of piperaquine and dihydroartemisinin has recently become the official first-line therapy in several Southeast Asian countries. The pharmacokinetic mismatching of these drugs, whose plasma half-lives are ~20 days and ~1 h, respectively, implies that recrudescent or new infections emerging shortly after treatment cessation will encounter piperaquine as a monotherapy agent. This creates substantial selection pressure for the emergence of resistance. To elucidate potential resistance determinants, we subjected cloned Plasmodium falciparum Dd2 parasites to continuous piperaquine pressure in vitro (47 nM; ~2-fold higher than the Dd2 50% inhibitory concentration [IC(50)]). The phenotype of outgrowth parasites was assayed in two clones, revealing an IC(50) against piperaquine of 2.1 µM and 1.7 µM, over 100-fold greater than that of the parent. To identify the genetic determinant of resistance, we employed comparative whole-genome hybridization analysis. Compared to the Dd2 parent, this analysis found (in both resistant clones) a novel single-nucleotide polymorphism in P. falciparum crt (pfcrt), deamplification of an 82-kb region of chromosome 5 (that includes pfmdr1), and amplification of an adjacent 63-kb region of chromosome 5. Continued propagation without piperaquine selection pressure resulted in "revertant" piperaquine-sensitive parasites. These retained the pfcrt polymorphism and further deamplified the chromosome 5 segment that encompasses pfmdr1; however, these two independently generated revertants both lost the neighboring 63-kb amplification. These results suggest that a copy number variation event on chromosome 5 (825600 to 888300) is associated with piperaquine resistance. Transgene expression studies are underway with individual genes in this segment to evaluate their contribution to piperaquine resistance.