Despite their known transforming properties, the effects of leukemogenic FLT3-ITD mutations on hematopoietic stem and multipotent progenitor cells and on hematopoietic differentiation are not well understood. We report a mouse model harboring an ITD in the murine Flt3 locus that develops myeloproliferative disease resembling CMML and further identified FLT3-ITD mutations in a subset of human CMML. These findings correlated with an increase in number, cell cycling, and survival of multipotent stem and progenitor cells in an ITD dose-dependent manner in animals that exhibited alterations within their myeloid progenitor compartments and a block in normal B cell development. This model provides insights into the consequences of constitutive signaling by an oncogenic tyrosine kinase on hematopoietic progenitor quiescence, function, and cell fate.

Approximately 18% of acute myeloid leukemia (AML) cases express a fusion transcript. However, few fusions are recurrent across AML and the identification of these rare chimeras is of interest to characterize AML patients. Here, we studied the transcriptome of 8 adult AML patients with poorly described chromosomal translocation(s), with the aim of identifying novel and rare fusion transcripts. We integrated RNA-sequencing data with multiple approaches including computational analysis, Sanger sequencing, fluorescence in situ hybridization and in vitro studies to assess the oncogenic potential of the ZEB2-BCL11B chimera. We detected 7 different fusions with partner genes involving transcription factors (OAZ-MAFK, ZEB2-BCL11B), tumor suppressors (SAV1-GYPB, PUF60-TYW1, CNOT2-WT1) and rearrangements associated with the loss of NF1 (CPD-PXT1, UTP6-CRLF3). Notably, ZEB2-BCL11B rearrangements co-occurred with FLT3 mutations and were associated with a poorly differentiated or mixed phenotype leukemia. Although the fusion alone did not transform murine c-Kit+ bone marrow cells, 45.4% of 14q32 non-rearranged AML cases were also BCL11B-positive, suggesting a more general and complex mechanism of leukemogenesis associated with BCL11B expression. Overall, by combining different approaches, we described rare fusion events contributing to the complexity of AML and we linked the expression of some chimeras to genomic alterations hitting known genes in AML.

Genome-wide maps of histone modifications have been obtained for several plant species. However, most studies focus on model systems and do not enforce FAIR data management principles. Here we study the H3K27me3 epigenome and associated transcriptome of Brassica rapa, an important vegetable cultivated worldwide.

Peripheral arterial disease (PAD) is associated with a high risk of cardiovascular events and death and is postulated to be a critical socioeconomic cost in the future. Extracellular vesicles (EVs) have emerged as potential candidates for new biomarker discovery related to their protein and nucleic acid cargo. In search of new prognostic and therapeutic targets in PAD, we determined the prothrombotic activity, the cellular origin and the transcriptomic profile of circulating EVs. This prospective study included control and PAD patients. Coagulation time (Procoag-PPL kit), EVs cellular origin and phosphatidylserine exposure were determined by flow cytometry in platelet-free plasma (n = 45 PAD). Transcriptomic profiles of medium/large EVs were generated using the MARS-Seq RNA-Seq protocol (n = 12/group). The serum concentration of the differentially expressed gene S100A9, in serum calprotectin (S100A8/A9), was validated by ELISA in control (n = 100) and PAD patients (n = 317). S100A9 was also determined in EVs and tissues of human atherosclerotic plaques (n = 3). Circulating EVs of PAD patients were mainly of platelet origin, predominantly Annexin V positive and were associated with the procoagulant activity of platelet-free plasma. Transcriptomic analysis of EVs identified 15 differentially expressed genes. Among them, serum calprotectin was elevated in PAD patients (p < 0.05) and associated with increased amputation risk before and after covariate adjustment (mean follow-up 3.6 years, p < 0.01). The combination of calprotectin with hs-CRP in the multivariate analysis further improved risk stratification (p < 0.01). Furthermore, S100A9 was also expressed in femoral plaque derived EVs and tissues. In summary, we found that PAD patients release EVs, mainly of platelet origin, highly positive for AnnexinV and rich in transcripts related to platelet biology and immune responses. Amputation risk prediction improved with calprotectin and was significantly higher when combined with hs-CRP. Our results suggest that EVs can be a promising component of liquid biopsy to identify the molecular signature of PAD patients.

CRISPR/Cas9 technology offers novel approaches for the development of new therapies for many unmet clinical needs, including a significant number of inherited monogenic diseases. However, in vivo correction of disease-causing genes is still inefficient, especially for those diseases without selective advantage for corrected cells. We reasoned that substrate reduction therapies (SRT) targeting non-essential enzymes could provide an attractive alternative. Here we evaluate the therapeutic efficacy of an in vivo CRISPR/Cas9-mediated SRT to treat primary hyperoxaluria type I (PH1), a rare inborn dysfunction in glyoxylate metabolism that results in excessive hepatic oxalate production causing end-stage renal disease. A single systemic administration of an AAV8-CRISPR/Cas9 vector targeting glycolate oxidase, prevents oxalate overproduction and kidney damage, with no signs of toxicity in Agxt1-/- mice. Our results reveal that CRISPR/Cas9-mediated SRT represents a promising therapeutic option for PH1 that can be potentially applied to other metabolic diseases caused by the accumulation of toxic metabolites.

Altered transcription is a cardinal feature of acute myeloid leukemia (AML); however, exactly how mutations synergize to remodel the epigenetic landscape and rewire three-dimensional DNA topology is unknown. Here, we apply an integrated genomic approach to a murine allelic series that models the two most common mutations in AML: Flt3-ITD and Npm1c. We then deconvolute the contribution of each mutation to alterations of the epigenetic landscape and genome organization, and infer how mutations synergize in the induction of AML. Our studies demonstrate that Flt3-ITD signals to chromatin to alter the epigenetic environment and synergizes with mutations in Npm1c to alter gene expression and drive leukemia induction. These analyses also allow the identification of long-range cis-regulatory circuits, including a previously unknown superenhancer of Hoxa locus, as well as larger and more detailed gene-regulatory networks, driven by transcription factors including PU.1 and IRF8, whose importance we demonstrate through perturbation of network members.

Recurrent gene mutations often cooperate in a predefined stepwise and synergistic manner to alter global transcription, through directly or indirectly remodeling epigenetic landscape on linear and three-dimensional (3D) scales. Here, we present a multiomics data integration approach to investigate the impact of gene mutational synergy on transcription, chromatin states, and 3D chromatin organization in a murine leukemia model. This protocol provides an executable framework to study epigenetic remodeling induced by cooperating gene mutations and to identify the critical regulatory network involved. For complete details on the use and execution of this protocol, please refer to Yun et al. (2021).

Leukaemia stem cells (LSCs) are the critical seed for the growth of haematological malignancies, driving the clonal expansion that enables disease initiation, relapse and often resistance. Specifically, they display inherent phenotypic and epigenetic plasticity resulting in complex heterogenic diseases. In this review, we discuss the key principles of deregulation of epigenetic processes that shape this disease evolution. We consider measures to define and quantify clonal heterogeneity, combining information from recent studies assessing mutational, transcriptional and epigenetic landscapes at single cell resolution in myeloid neoplasms (MN). We highlight the importance of integrating epigenetic and genetic information to better understand inter- and intra-patient heterogeneity and discuss how this understanding further informs evolution and progression trajectories and subsequent clinical response in MN. Under this topic, we also discuss efforts to identify mechanisms of resistance, by longitudinal analyses of patient samples. Finally, we highlight how we might target these aberrant epigenetic processes for better therapeutic outcomes and to potentially eradicate LSCs.

Haematopoietic stem cells (HSC) reside in the bone marrow microenvironment (BMM), where they respond to extracellular calcium [eCa2+] via the G-protein coupled calcium-sensing receptor (CaSR). Here we show that a calcium gradient exists in this BMM, and that [eCa2+] and response to [eCa2+] differ between leukaemias. CaSR influences the location of MLL-AF9+ acute myeloid leukaemia (AML) cells within this niche and differentially impacts MLL-AF9+ AML versus BCR-ABL1+ leukaemias. Deficiency of CaSR reduces AML leukaemic stem cells (LSC) 6.5-fold. CaSR interacts with filamin A, a crosslinker of actin filaments, affects stemness-associated factors and modulates pERK, β-catenin and c-MYC signaling and intracellular levels of [Ca2+] in MLL-AF9+ AML cells. Combination treatment of cytarabine plus CaSR-inhibition in various models may be superior to cytarabine alone. Our studies suggest CaSR to be a differential and targetable factor in leukaemia progression influencing self-renewal of AML LSC via [eCa2+] cues from the BMM.

Epigenetic regulators, such as EZH2, are frequently mutated in cancer, and loss-of-function EZH2 mutations are common in myeloid malignancies. We have examined the importance of cellular context for Ezh2 loss during the evolution of acute myeloid leukemia (AML), where we observed stage-specific and diametrically opposite functions for Ezh2 at the early and late stages of disease. During disease maintenance, WT Ezh2 exerts an oncogenic function that may be therapeutically targeted. In contrast, Ezh2 acts as a tumor suppressor during AML induction. Transcriptional analysis explains this apparent paradox, demonstrating that loss of Ezh2 derepresses different expression programs during disease induction and maintenance. During disease induction, Ezh2 loss derepresses a subset of bivalent promoters that resolve toward gene activation, inducing a feto-oncogenic program that includes genes such as Plag1, whose overexpression phenocopies Ezh2 loss to accelerate AML induction in mouse models. Our data highlight the importance of cellular context and disease phase for the function of Ezh2 and its potential therapeutic implications.

Loss-of-function mutations of cyclic-AMP response element binding protein, binding protein (CREBBP) are prevalent in lymphoid malignancies. However, the tumour suppressor functions of CREBBP remain unclear. We demonstrate that loss of Crebbp in murine haematopoietic stem and progenitor cells (HSPCs) leads to increased development of B-cell lymphomas. This is preceded by accumulation of hyperproliferative lymphoid progenitors with a defective DNA damage response (DDR) due to a failure to acetylate p53. We identify a premalignant lymphoma stem cell population with decreased H3K27ac, which undergoes transcriptional and genetic evolution due to the altered DDR, resulting in lymphomagenesis. Importantly, when Crebbp is lost later in lymphopoiesis, cellular abnormalities are lost and tumour generation is attenuated. We also document that CREBBP mutations may occur in HSPCs from patients with CREBBP-mutated lymphoma. These data suggest that earlier loss of Crebbp is advantageous for lymphoid transformation and inform the cellular origins and subsequent evolution of lymphoid malignancies.

Acute myeloid leukemia (AML) is the most common acute leukemia in adults and has an unacceptably low cure rate. In recent years, a number of new treatment strategies and compounds were developed for the treatment of AML. There were several randomized controlled clinical trials with the objective to improve patients' management and patients' outcome in AML. Unfortunately, these trials are not always directly comparable since they do not measure the same outcomes, and currently there are no core outcome sets that can be used to guide outcome selection and harmonization in this disease area. The HARMONY (Healthcare Alliance for Resourceful Medicine Offensive against Neoplasms in Hematology) Alliance is a public-private European network established in 2017 and currently includes 53 partners and 32 associated members from 22 countries. Amongst many other goals of the HARMONY Alliance, Work Package 2 focuses on defining outcomes that are relevant to each hematological malignancy. Accordingly, this pilot study will be performed to define a core outcome set in AML.

Multiple myeloma (MM) is a plasma cell neoplasm associated with a broad variety of genetic lesions. In spite of this genetic heterogeneity, MMs share a characteristic malignant phenotype whose underlying molecular basis remains poorly characterized. In the present study, we examined plasma cells from MM using a multi-epigenomics approach and demonstrated that, when compared to normal B cells, malignant plasma cells showed an extensive activation of regulatory elements, in part affecting coregulated adjacent genes. Among target genes up-regulated by this process, we found members of the NOTCH, NF-kB, MTOR signaling, and TP53 signaling pathways. Other activated genes included sets involved in osteoblast differentiation and response to oxidative stress, all of which have been shown to be associated with the MM phenotype and clinical behavior. We functionally characterized MM-specific active distant enhancers controlling the expression of thioredoxin (TXN), a major regulator of cellular redox status and, in addition, identified PRDM5 as a novel essential gene for MM. Collectively, our data indicate that aberrant chromatin activation is a unifying feature underlying the malignant plasma cell phenotype.

Clinical recommendations for Acute Myeloid Leukemia (AML) classification and risk-stratification remain heavily reliant on cytogenetic findings at diagnosis, which are present in <50% of patients. Using comprehensive molecular profiling data from 3,653 patients we characterize and validate 16 molecular classes describing 100% of AML patients. Each class represents diverse biological AML subgroups, and is associated with distinct clinical presentation, likelihood of response to induction chemotherapy, risk of relapse and death over time. Secondary AML-2, emerges as the second largest class (24%), associates with high-risk disease, poor prognosis irrespective of flow Minimal Residual Disease (MRD) negativity, and derives significant benefit from transplantation. Guided by class membership we derive a 3-tier risk-stratification score that re-stratifies 26% of patients as compared to standard of care. This results in a unified framework for disease classification and risk-stratification in AML that relies on information from cytogenetics and 32 genes. Last, we develop an open-access patient-tailored clinical decision support tool.

Innate lymphoid cells (ILCs) are critical modulators of mucosal immunity, inflammation, and tissue homeostasis, but their full spectrum of cellular states and regulatory landscapes remains elusive. Here, we combine genome-wide RNA-seq, ChIP-seq, and ATAC-seq to compare the transcriptional and epigenetic identity of small intestinal ILCs, identifying thousands of distinct gene profiles and regulatory elements. Single-cell RNA-seq and flow and mass cytometry analyses reveal compartmentalization of cytokine expression and metabolic activity within the three classical ILC subtypes and highlight transcriptional states beyond the current canonical classification. In addition, using antibiotic intervention and germ-free mice, we characterize the effect of the microbiome on the ILC regulatory landscape and determine the response of ILCs to microbial colonization at the single-cell level. Together, our work characterizes the spectrum of transcriptional identities of small intestinal ILCs and describes how ILCs differentially integrate signals from the microbial microenvironment to generate phenotypic and functional plasticity.

Histone acetyltransferases (HATs) catalyze the transfer of an acetyl group from acetyl-CoA to lysine residues of histones and play a central role in transcriptional regulation in diverse biological processes. Dysregulation of HAT activity can lead to human diseases including developmental disorders and cancer. Through genome-wide CRISPR-Cas9 screens, we identified several HATs of the MYST family as fitness genes for acute myeloid leukemia (AML). Here we investigate the essentiality of lysine acetyltransferase KAT7 in AMLs driven by the MLL-X gene fusions. We found that KAT7 loss leads to a rapid and complete loss of both H3K14ac and H4K12ac marks, in association with reduced proliferation, increased apoptosis, and differentiation of AML cells. Acetyltransferase activity of KAT7 is essential for the proliferation of these cells. Mechanistically, our data propose that acetylated histones provide a platform for the recruitment of MLL-fusion-associated adaptor proteins such as BRD4 and AF4 to gene promoters. Upon KAT7 loss, these factors together with RNA polymerase II rapidly dissociate from several MLL-fusion target genes that are essential for AML cell proliferation, including MEIS1, PBX3, and SENP6. Our findings reveal that KAT7 is a plausible therapeutic target for this poor prognosis AML subtype.

Drug combinations are key to circumvent resistance mechanisms compromising response to single anti-cancer targeted therapies. The implementation of combinatorial approaches involving MEK1/2 or KRASG12C inhibitors in the context of KRAS-mutated lung cancers focuses fundamentally on targeting KRAS proximal activators or effectors. However, the antitumor effect is highly determined by compensatory mechanisms arising in defined cell types or tumor subgroups. A potential strategy to find drug combinations targeting a larger fraction of KRAS-mutated lung cancers may capitalize on the common, distal gene expression output elicited by oncogenic KRAS. By integrating a signature-driven drug repurposing approach with a pairwise pharmacological screen, here we show synergistic drug combinations consisting of multi-tyrosine kinase PKC inhibitors together with MEK1/2 or KRASG12C inhibitors. Such combinations elicit a cytotoxic response in both in vitro and in vivo models, which in part involves inhibition of the PKC inhibitor target AURKB. Proteome profiling links dysregulation of MYC expression to the effect of both PKC inhibitor-based drug combinations. Furthermore, MYC overexpression appears as a resistance mechanism to MEK1/2 and KRASG12C inhibitors. Our study provides a rational framework for selecting drugs entering combinatorial strategies and unveils MEK1/2- and KRASG12C-based therapies for lung cancer.

The introduction of highly selective ABL-tyrosine kinase inhibitors (TKIs) has revolutionized therapy for chronic myeloid leukemia (CML). However, TKIs are only efficacious in the chronic phase of the disease and effective therapies for TKI-refractory CML, or after progression to blast crisis (BC), are lacking. Whereas the chronic phase of CML is dependent on BCR-ABL, additional mutations are required for progression to BC. However, the identity of these mutations and the pathways they affect are poorly understood, hampering our ability to identify therapeutic targets and improve outcomes. Here, we describe a novel mouse model that allows identification of mechanisms of BC progression in an unbiased and tractable manner, using transposon-based insertional mutagenesis on the background of chronic phase CML. Our BC model is the first to faithfully recapitulate the phenotype, cellular and molecular biology of human CML progression. We report a heterogeneous and unique pattern of insertions identifying known and novel candidate genes and demonstrate that these pathways drive disease progression and provide potential targets for novel therapeutic strategies. Our model greatly informs the biology of CML progression and provides a potent resource for the development of candidate therapies to improve the dismal outcomes in this highly aggressive disease.

Chronic myeloid leukemia (CML) is initiated and maintained by the tyrosine kinase BCR-ABL which activates a number of signal transduction pathways, including PI3K/AKT signaling and consequently inactivates FOXO transcription factors. ABL-specific tyrosine kinase inhibitors (TKIs) induce minimal apoptosis in CML progenitor cells, yet exert potent antiproliferative effects, through as yet poorly understood mechanisms. Here, we demonstrate that in CD34+ CML cells, FOXO1 and 3a are inactivated and relocalized to the cytoplasm by BCR-ABL activity. TKIs caused a decrease in phosphorylation of FOXOs, leading to their relocalization from cytoplasm (inactive) to nucleus (active), where they modulated the expression of key FOXO target genes, such as Cyclin D1, ATM, CDKN1C, and BCL6 and induced G1 arrest. Activation of FOXO1 and 3a and a decreased expression of their target gene Cyclin D1 were also observed after 6 days of in vivo treatment with dasatinib in a CML transgenic mouse model. The over-expression of FOXO3a in CML cells combined with TKIs to reduce proliferation, with similar results seen for inhibitors of PI3K/AKT/mTOR signaling. While stable expression of an active FOXO3a mutant induced a similar level of quiescence to TKIs alone, shRNA-mediated knockdown of FOXO3a drove CML cells into cell cycle and potentiated TKI-induced apoptosis. These data demonstrate that TKI-induced G1 arrest in CML cells is mediated through inhibition of the PI3K/AKT pathway and reactivation of FOXOs. This enhanced understanding of TKI activity and induced progenitor cell quiescence suggests that new therapeutic strategies for CML should focus on manipulation of this signaling network.

In multicellular organisms, dedicated regulatory circuits control cell type diversity and responses. The crosstalk and redundancies within these circuits and substantial cellular heterogeneity pose a major research challenge. Here, we present CRISP-seq, an integrated method for massively parallel single-cell RNA sequencing (RNA-seq) and clustered regularly interspaced short palindromic repeats (CRISPR)-pooled screens. We show that profiling the genomic perturbation and transcriptome in the same cell enables us to simultaneously elucidate the function of multiple factors and their interactions. We applied CRISP-seq to probe regulatory circuits of innate immunity. By sampling tens of thousands of perturbed cells in vitro and in mice, we identified interactions and redundancies between developmental and signaling-dependent factors. These include opposing effects of Cebpb and Irf8 in regulating the monocyte/macrophage versus dendritic cell lineages and differential functions for Rela and Stat1/2 in monocyte versus dendritic cell responses to pathogens. This study establishes CRISP-seq as a broadly applicable, comprehensive, and unbiased approach for elucidating mammalian regulatory circuits.