Epigenetic balance ensures mechanistic control of MLL amplification and rearrangement.

MLL/KMT2A amplifications and translocations are prevalent in infant, adult, and therapy-induced leukemia. However, the molecular contributor(s) to these alterations are unclear. Here, we demonstrate that histone H3 lysine 9 mono- and di-methylation (H3K9me1/2) balance at the MLL/KMT2A locus regulates these amplifications and rearrangements. This balance is controlled by the crosstalk between lysine demethylase KDM3B and methyltransferase G9a/EHMT2. KDM3B depletion increases H3K9me1/2 levels and reduces CTCF occupancy at the MLL/KMT2A locus, in turn promoting amplification and rearrangements. Depleting CTCF is also sufficient to generate these focal alterations. Furthermore, the chemotherapy doxorubicin (Dox), which associates with therapy-induced leukemia and promotes MLL/KMT2A amplifications and rearrangements, suppresses KDM3B and CTCF protein levels. KDM3B and CTCF overexpression rescues Dox-induced MLL/KMT2A alterations. G9a inhibition in human cells or mice also suppresses MLL/KMT2A events accompanying Dox treatment. Therefore, MLL/KMT2A amplifications and rearrangements are controlled by epigenetic regulators that are tractable drug targets, which has clinical implications.

Temporal dynamics and genomic programming of plasma cell fates.

Affinity-matured plasma cells (PCs) of varying lifespans are generated through a germinal center (GC) response. The developmental dynamics and genomic programs of antigen-specific PC precursors remain to be elucidated. Here, using a model antigen in mice, we demonstrate biphasic generation of PC precursors, with those generating long-lived bone marrow PCs preferentially produced in the late phase of GC response. Clonal tracing using single-cell RNA sequencing and B cell antigen receptor sequencing in spleen and bone marrow compartments, coupled with adoptive transfer experiments, reveals a new PC transition state that gives rise to functionally competent PC precursors. The latter undergo clonal expansion, dependent on inducible expression of TIGIT. We propose a model for the proliferation and programming of precursors of long-lived PCs, based on extended antigen encounters in the GC.

ThPOK is a critical multifaceted regulator of myeloid lineage development.

The transcription factor ThPOK (encoded by Zbtb7b) is well known for its role as a master regulator of CD4 lineage commitment in the thymus. Here, we report an unexpected and critical role of ThPOK as a multifaceted regulator of myeloid lineage commitment, differentiation and maturation. Using reporter and knockout mouse models combined with single-cell RNA-sequencing, progenitor transfer and colony assays, we show that ThPOK controls monocyte-dendritic cell versus granulocyte lineage production during homeostatic differentiation, and serves as a brake for neutrophil maturation in granulocyte lineage-specified cells through transcriptional regulation of lineage-specific transcription factors and RNA via altered messenger RNA splicing to reprogram intron retention.

Monocytes re-enter the bone marrow during fasting and alter the host response to infection.

Diet profoundly influences physiology. Whereas over-nutrition elevates risk for disease via its influence on immunity and metabolism, caloric restriction and fasting appear to be salutogenic. Despite multiple correlations observed between diet and health, the underlying biology remains unclear. Here, we identified a fasting-induced switch in leukocyte migration that prolongs monocyte lifespan and alters susceptibility to disease in mice. We show that fasting during the active phase induced the rapid return of monocytes from the blood to the bone marrow. Monocyte re-entry was orchestrated by hypothalamic-pituitary-adrenal (HPA) axis-dependent release of corticosterone, which augmented the CXCR4 chemokine receptor. Although the marrow is a safe haven for monocytes during nutrient scarcity, re-feeding prompted mobilization culminating in monocytosis of chronologically older and transcriptionally distinct monocytes. These shifts altered response to infection. Our study shows that diet-in particular, a diet's temporal dynamic balance-modulates monocyte lifespan with consequences for adaptation to external stressors.

Brain motor and fear circuits regulate leukocytes during acute stress.

The nervous and immune systems are intricately linked1. Although psychological stress is known to modulate immune function, mechanistic pathways linking stress networks in the brain to peripheral leukocytes remain poorly understood2. Here we show that distinct brain regions shape leukocyte distribution and function throughout the body during acute stress in mice. Using optogenetics and chemogenetics, we demonstrate that motor circuits induce rapid neutrophil mobilization from the bone marrow to peripheral tissues through skeletal-muscle-derived neutrophil-attracting chemokines. Conversely, the paraventricular hypothalamus controls monocyte and lymphocyte egress from secondary lymphoid organs and blood to the bone marrow through direct, cell-intrinsic glucocorticoid signalling. These stress-induced, counter-directional, population-wide leukocyte shifts are associated with altered disease susceptibility. On the one hand, acute stress changes innate immunity by reprogramming neutrophils and directing their recruitment to sites of injury. On the other hand, corticotropin-releasing hormone neuron-mediated leukocyte shifts protect against the acquisition of autoimmunity, but impair immunity to SARS-CoV-2 and influenza infection. Collectively, these data show that distinct brain regions differentially and rapidly tailor the leukocyte landscape during psychological stress, therefore calibrating the ability of the immune system to respond to physical threats.

Ontogeny and function of the circadian clock in intestinal organoids.

Circadian rhythms regulate diverse aspects of gastrointestinal physiology ranging from the composition of microbiota to motility. However, development of the intestinal circadian clock and detailed mechanisms regulating circadian physiology of the intestine remain largely unknown. In this report, we show that both pluripotent stem cell-derived human intestinal organoids engrafted into mice and patient-derived human intestinal enteroids possess circadian rhythms and demonstrate circadian phase-dependent necrotic cell death responses to Clostridium difficile toxin B (TcdB). Intriguingly, mouse and human enteroids demonstrate anti-phasic necrotic cell death responses to TcdB. RNA-Seq analysis shows that ~3-10% of the detectable transcripts are rhythmically expressed in mouse and human enteroids. Remarkably, we observe anti-phasic gene expression of Rac1, a small GTPase directly inactivated by TcdB, between mouse and human enteroids, and disruption of Rac1 abolishes clock-dependent necrotic cell death responses. Our findings uncover robust functions of circadian rhythms regulating clock-controlled genes in both mouse and human enteroids governing organism-specific, circadian phase-dependent necrotic cell death responses, and lay a foundation for human organ- and disease-specific investigation of clock functions using human organoids for translational applications.

Mouse models of neutropenia reveal progenitor-stage-specific defects.

Advances in genetics and sequencing have identified a plethora of disease-associated and disease-causing genetic alterations. To determine causality between genetics and disease, accurate models for molecular dissection are required; however, the rapid expansion of transcriptional populations identified through single-cell analyses presents a major challenge for accurate comparisons between mutant and wild-type cells. Here we generate mouse models of human severe congenital neutropenia (SCN) using patient-derived mutations in the GFI1 transcription factor. To determine the effects of SCN mutations, we generated single-cell references for granulopoietic genomic states with linked epitopes1, aligned mutant cells to their wild-type equivalents and identified differentially expressed genes and epigenetic loci. We find that GFI1-target genes are altered sequentially, as cells go through successive states of differentiation. These insights facilitated the genetic rescue of granulocytic specification but not post-commitment defects in innate immune effector function, and underscore the importance of evaluating the effects of mutations and therapy within each relevant cell state.

Stiffness Restricts the Stemness of the Intestinal Stem Cells and Skews Their Differentiation Toward Goblet Cells.

BACKGROUND & AIMS: Fibrosis and tissue stiffening are hallmarks of inflammatory bowel disease (IBD). We have hypothesized that the increased stiffness directly contributes to the dysregulation of the epithelial cell homeostasis in IBD. Here, we aim to determine the impact of tissue stiffening on the fate and function of the intestinal stem cells (ISCs). METHODS: We developed a long-term culture system consisting of 2.5-dimensional intestinal organoids grown on a hydrogel matrix with tunable stiffness. Single-cell RNA sequencing provided stiffness-regulated transcriptional signatures of the ISCs and their differentiated progeny. YAP-knockout and YAP-overexpression mice were used to manipulate YAP expression. In addition, we analyzed colon samples from murine colitis models and human IBD samples to assess the impact of stiffness on ISCs in vivo. RESULTS: We demonstrated that increasing the stiffness potently reduced the population of LGR5+ ISCs and KI-67+-proliferating cells. Conversely, cells expressing the stem cell marker, olfactomedin-4, became dominant in the crypt-like compartments and pervaded the villus-like regions. Concomitantly, stiffening prompted the ISCs to preferentially differentiate toward goblet cells. Mechanistically, stiffening increased the expression of cytosolic YAP, driving the extension of olfactomedin-4+ cells into the villus-like regions, while it induced the nuclear translocation of YAP, leading to preferential differentiation of ISCs toward goblet cells. Furthermore, analysis of colon samples from murine colitis models and patients with IBD demonstrated cellular and molecular remodeling reminiscent of those observed in vitro. CONCLUSIONS: Collectively, our findings highlight that matrix stiffness potently regulates the stemness of ISCs and their differentiation trajectory, supporting the hypothesis that fibrosis-induced gut stiffening plays a direct role in epithelial remodeling in IBD.

Erythroblastic islands foster granulopoiesis in parallel to terminal erythropoiesis.

The erythroblastic island (EBI), composed of a central macrophage surrounded by maturing erythroblasts, is the erythroid precursor niche. Despite numerous studies, its precise composition is still unclear. Using multispectral imaging flow cytometry, in vitro island reconstitution, and single-cell RNA sequencing of adult mouse bone marrow (BM) EBI-component cells enriched by gradient sedimentation, we present evidence that the CD11b+ cells present in the EBIs are neutrophil precursors specifically associated with BM EBI macrophages, indicating that erythro-(myelo)-blastic islands are a site for terminal granulopoiesis and erythropoiesis. We further demonstrate that the balance between these dominant and terminal differentiation programs is dynamically regulated within this BM niche by pathophysiological states that favor granulopoiesis during anemia of inflammation and favor erythropoiesis after erythropoietin stimulation. Finally, by molecular profiling, we reveal the heterogeneity of EBI macrophages by cellular indexing of transcriptome and epitope sequencing of mouse BM EBIs at baseline and after erythropoietin stimulation in vivo and provide a searchable online viewer of these data characterizing the macrophage subsets serving as hematopoietic niches. Taken together, our findings demonstrate that EBIs serve a dual role as niches for terminal erythropoiesis and granulopoiesis and the central macrophages adapt to optimize production of red blood cells or neutrophils.

Intragraft B cell differentiation during the development of tolerance to kidney allografts is associated with a regulatory B cell signature revealed by single cell transcriptomics.

Mouse kidney allografts are spontaneously accepted in select, fully mismatched donor-recipient strain combinations, like DBA/2J to C57BL/6 (B6), by natural tolerance. We previously showed accepted renal grafts form aggregates containing various immune cells within 2 weeks posttransplant, referred to as regulatory T cell-rich organized lymphoid structures, which are a novel regulatory tertiary lymphoid organ. To characterize the cells within T cell-rich organized lymphoid structures, we performed single-cell RNA sequencing on CD45+ sorted cells from accepted and rejected renal grafts from 1-week to 6-months posttransplant. Analysis of single-cell RNA sequencing data revealed a shifting from a T cell-dominant to a B cell-rich population by 6 months with an increased regulatory B cell signature. Furthermore, B cells were a greater proportion of the early infiltrating cells in accepted vs rejecting grafts. Flow cytometry of B cells at 20 weeks posttransplant revealed T cell, immunoglobulin domain and mucin domain-1+ B cells, potentially implicating a regulatory role in the maintenance of allograft tolerance. Lastly, B cell trajectory analysis revealed intragraft differentiation from precursor B cells to memory B cells in accepted allografts. In summary, we show a shifting T cell- to B cell-rich environment and a differential cellular pattern among accepted vs rejecting kidney allografts, possibly implicating B cells in the maintenance of kidney allograft acceptance.

Pbx4 limits heart size and fosters arch artery formation by partitioning second heart field progenitors and restricting proliferation.

Vertebrate heart development requires the integration of temporally distinct differentiating progenitors. However, few signals are understood that restrict the size of the later-differentiating outflow tract (OFT). We show that improper specification and proliferation of second heart field (SHF) progenitors in zebrafish lazarus (lzr) mutants, which lack the transcription factor Pbx4, produces enlarged hearts owing to an increase in ventricular and smooth muscle cells. Specifically, Pbx4 initially promotes the partitioning of the SHF into anterior progenitors, which contribute to the OFT, and adjacent endothelial cell progenitors, which contribute to posterior pharyngeal arches. Subsequently, Pbx4 limits SHF progenitor (SHFP) proliferation. Single cell RNA sequencing of nkx2.5+ cells revealed previously unappreciated distinct differentiation states and progenitor subpopulations that normally reside within the SHF and arterial pole of the heart. Specifically, the transcriptional profiles of Pbx4-deficient nkx2.5+ SHFPs are less distinct and display characteristics of normally discrete proliferative progenitor and anterior, differentiated cardiomyocyte populations. Therefore, our data indicate that the generation of proper OFT size and arch arteries requires Pbx-dependent stratification of unique differentiation states to facilitate both homeotic-like transformations and limit progenitor production within the SHF.

Maturation of heart valve cell populations during postnatal remodeling.

Heart valve cells mediate extracellular matrix (ECM) remodeling during postnatal valve leaflet stratification, but phenotypic and transcriptional diversity of valve cells in development is largely unknown. Single cell analysis of mouse heart valve cells was used to evaluate cell heterogeneity during postnatal ECM remodeling and leaflet morphogenesis. The transcriptomic analysis of single cells from postnatal day (P)7 and P30 murine aortic (AoV) and mitral (MV) heart valves uncovered distinct subsets of melanocytes, immune and endothelial cells present at P7 and P30. By contrast, interstitial cell populations are different from P7 to P30. P7 valve leaflets exhibit two distinct collagen- and glycosaminoglycan-expressing interstitial cell clusters, and prevalent ECM gene expression. At P30, four interstitial cell clusters are apparent with leaflet specificity and differential expression of complement factors, ECM proteins and osteogenic genes. This initial transcriptomic analysis of postnatal heart valves at single cell resolution demonstrates that subpopulations of endothelial and immune cells are relatively constant throughout postnatal development, but interstitial cell subpopulations undergo changes in gene expression and cellular functions in primordial and mature valves.

The Rhesus Macaque Serves As a Model for Human Lateral Branch Nephrogenesis.

BACKGROUND: Most nephrons are added in late gestation. Truncated extrauterine nephrogenesis in premature infants results in fewer nephrons and significantly increased risk for CKD in adulthood. To overcome the ethical and technical difficulties associated with studies of late-gestation human fetal kidney development, third-trimester rhesus macaques served as a model to understand lateral branch nephrogenesis (LBN) at the molecular level. METHODS: Immunostaining and 3D rendering assessed morphology. Single-cell (sc) and single-nucleus (sn) RNA-Seq were performed on four cortically enriched fetal rhesus kidneys of 129-131 days gestational age (GA). An integrative bioinformatics strategy was applied across single-cell modalities, species, and time. RNAScope validation studies were performed on human archival tissue. RESULTS: Third-trimester rhesus kidney undergoes human-like LBN. scRNA-Seq of 23,608 cells revealed 37 transcriptionally distinct cell populations, including naïve nephron progenitor cells (NPCs), with the prior noted marker genes CITED1, MEOX1, and EYA1 (c25). These same populations and markers were reflected in snRNA-Seq of 5972 nuclei. Late-gestation rhesus NPC markers resembled late-gestation murine NPC, whereas early second-trimester human NPC markers aligned to midgestation murine NPCs. New, age-specific rhesus NPCs (SHISA8) and ureteric buds (POU3F4 and TWIST) predicted markers were verified in late-gestation human archival samples. CONCLUSIONS: Rhesus macaque is the first model of bona fide LBN, enabling molecular studies of late gestation, human-like nephrogenesis. These molecular findings support the hypothesis that aging nephron progenitors have a distinct molecular signature and align to their earlier human counterparts, with unique markers highlighting LBN-specific progenitor maturation.

Resolving single-cell heterogeneity from hundreds of thousands of cells through sequential hybrid clustering and NMF.

MOTIVATION: The rapid proliferation of single-cell RNA-sequencing (scRNA-Seq) technologies has spurred the development of diverse computational approaches to detect transcriptionally coherent populations. While the complexity of the algorithms for detecting heterogeneity has increased, most require significant user-tuning, are heavily reliant on dimension reduction techniques and are not scalable to ultra-large datasets. We previously described a multi-step algorithm, Iterative Clustering and Guide-gene Selection (ICGS), which applies intra-gene correlation and hybrid clustering to uniquely resolve novel transcriptionally coherent cell populations from an intuitive graphical user interface. RESULTS: We describe a new iteration of ICGS that outperforms state-of-the-art scRNA-Seq detection workflows when applied to well-established benchmarks. This approach combines multiple complementary subtype detection methods (HOPACH, sparse non-negative matrix factorization, cluster 'fitness', support vector machine) to resolve rare and common cell-states, while minimizing differences due to donor or batch effects. Using data from multiple cell atlases, we show that the PageRank algorithm effectively downsamples ultra-large scRNA-Seq datasets, without losing extremely rare or transcriptionally similar yet distinct cell types and while recovering novel transcriptionally distinct cell populations. We believe this new approach holds tremendous promise in reproducibly resolving hidden cell populations in complex datasets. AVAILABILITY AND IMPLEMENTATION: ICGS2 is implemented in Python. The source code and documentation are available at http://altanalyze.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

cellHarmony: cell-level matching and holistic comparison of single-cell transcriptomes.

To understand the molecular pathogenesis of human disease, precision analyses to define alterations within and between disease-associated cell populations are desperately needed. Single-cell genomics represents an ideal platform to enable the identification and comparison of normal and diseased transcriptional cell populations. We created cellHarmony, an integrated solution for the unsupervised analysis, classification, and comparison of cell types from diverse single-cell RNA-Seq datasets. cellHarmony efficiently and accurately matches single-cell transcriptomes using a community-clustering and alignment strategy to compute differences in cell-type specific gene expression over potentially dozens of cell populations. Such transcriptional differences are used to automatically identify distinct and shared gene programs among cell-types and identify impacted pathways and transcriptional regulatory networks to understand the impact of perturbations at a systems level. cellHarmony is implemented as a python package and as an integrated workflow within the software AltAnalyze. We demonstrate that cellHarmony has improved or equivalent performance to alternative label projection methods, is able to identify the likely cellular origins of malignant states, stratify patients into clinical disease subtypes from identified gene programs, resolve discrete disease networks impacting specific cell-types, and illuminate therapeutic mechanisms. Thus, this approach holds tremendous promise in revealing the molecular and cellular origins of complex disease.

SKI controls MDS-associated chronic TGF-ß signaling, aberrant splicing, and stem cell fitness.

The transforming growth factor beta (TGF-ß) signaling pathway controls hematopoietic stem cell (HSC) behavior in the marrow niche; however, TGF-ß signaling becomes chronic in early-stage myelodysplastic syndrome (MDS). Although TGF-ß signaling normally induces negative feedback, in early-stage MDS, high levels of microRNA-21 (miR-21) contribute to chronic TGF-ß signaling. We found that a TGF-ß signal-correlated gene signature is sufficient to identify an MDS patient population with abnormal RNA splicing (eg, CSF3R) independent of splicing factor mutations and coincident with low HNRNPK activity. Levels of SKI messenger RNA (mRNA) encoding a TGF-ß antagonist are sufficient to identify these patients. However, MDS patients with high SKI mRNA and chronic TGF-ß signaling lack SKI protein because of miR-21 activity. To determine the impact of SKI loss, we examined murine Ski -/- HSC function. First, competitive HSC transplants revealed a profound defect in stem cell fitness (competitive disadvantage) but not specification, homing, or multilineage production. Aged recipients of Ski -/- HSCs exhibited mild phenotypes similar to phenotypes in those with macrocytic anemia. Second, blastocyst complementation revealed a dramatic block in Ski -/- hematopoiesis in the absence of transplantation. Similar to SKI-high MDS patient samples, Ski -/- HSCs strikingly upregulated TGF-ß signaling and deregulated expression of spliceosome genes (including Hnrnpk). Moreover, novel single-cell splicing analyses demonstrated that Ski -/- HSCs and high levels of SKI expression in MDS patient samples share abnormal alternative splicing of common genes (including those that encode splicing factors). We conclude that miR-21-mediated loss of SKI activates TGF-ß signaling and alternative splicing to impair the competitive advantage of normal HSCs (fitness), which could contribute to selection of early-stage MDS-genic clones.

Cross-platform single cell analysis of kidney development shows stromal cells express Gdnf.

The developing kidney provides a useful model for study of the principles of organogenesis. In this report we use three independent platforms, Drop-Seq, Chromium 10x Genomics and Fluidigm C1, to carry out single cell RNA-Seq (scRNA-Seq) analysis of the E14.5 mouse kidney. Using the software AltAnalyze, in conjunction with the unsupervised approach ICGS, we were unable to identify and confirm the presence of 16 distinct cell populations during this stage of active nephrogenesis. Using a novel integrative supervised computational strategy, we were able to successfully harmonize and compare the cell profiles across all three technological platforms. Analysis of possible cross compartment receptor/ligand interactions identified the nephrogenic zone stroma as a source of GDNF. This was unexpected because the cap mesenchyme nephron progenitors had been thought to be the sole source of GDNF, which is a key driver of branching morphogenesis of the collecting duct system. The expression of Gdnf by stromal cells was validated in several ways, including Gdnf in situ hybridization combined with immunohistochemistry for SIX2, and marker of nephron progenitors, and MEIS1, a marker of stromal cells. Finally, the single cell gene expression profiles generated in this study confirmed and extended previous work showing the presence of multilineage priming during kidney development. Nephron progenitors showed stochastic expression of genes associated with multiple potential differentiation lineages.

Caenorhabditis elegans SynMuv B gene activity is down-regulated during a viral infection to enhance RNA interference.

Small RNA pathways regulate eukaryotic antiviral defense. Many of the Caenorhabditis elegans mutations that were identified based on their enhanced RNAi, the synMuv B genes, also emerged from unrelated genetic screens for increased growth factor signaling. The dozen synMuv B genes encode homologues of the mammalian dREAM complex found in nearly all animals and plants, which includes the lin-35/retinoblastoma oncogene. We show that a set of highly induced mRNAs in synMuv B mutants is congruent with mRNAs induced by Orsay RNA virus infection of C. elegans. In wild type animals, a combination of a synMuv A mutation and a synMuv B mutation are required for the Muv phenotype of increased growth factor signaling. But we show that Orsay virus infection of a single synMuv A mutant can induce a Muv phenotype, unlike the uninfected single synMuv A mutant. This suggests that decreased synMuv B activity, which activates the antiviral RNAi pathway, is a defense response to viral infection. Small RNA deep sequencing analysis of various dREAM complex mutants uncovers distinct siRNA profiles indicative of such an siRNA response. We conclude that the synMuv B mutants maintain an antiviral readiness state even in the absence of actual infection. The enhanced RNAi and conservation of the dREAM complex mutants suggests new therapeutic avenues to boost antiviral defenses.

An Adipo-Pulmonary Axis Mediated by FABP4 Hormone Defines a Therapeutic Target Against Obesity-Induced Airway Disease.

Obesity-related airway disease is a clinical condition without a clear description and effective treatment. Here, we define this pathology and its unique properties, which differ from classic asthma phenotypes, and identify a novel adipo-pulmonary axis mediated by FABP4 hormone as a critical mediator of obesity-induced airway disease. Through detailed analysis of murine models and human samples, we elucidate the dysregulated lipid metabolism and immunometabolic responses within obese lungs, particularly highlighting the stress response activation and downregulation of surfactant-related genes, notably SftpC. We demonstrate that FABP4 deficiency mitigates these alterations, demonstrating a key role in obesity-induced airway disease pathogenesis. Importantly, we identify adipose tissue as the source of FABP4 hormone in the bronchoalveolar space and describe strong regulation in the context of human obesity, particularly among women. Finally, our exploration of antibody-mediated targeting of circulating FABP4 unveils a novel therapeutic avenue, addressing a pressing unmet need in managing obesity-related airway disease. These findings not only define the presence of a critical adipo-pulmonary endocrine link but also present FABP4 as a therapeutic target for managing this unique airway disease that we refer to as fatty lung disease associated with obesity. One Sentence Summary: Investigating FABP4's pivotal role in obesity-driven airway disease, this study unveils an adipo-pulmonary axis with potential therapeutic implications.

Leukemia aggressiveness is driven by chromatin remodeling and expression changes of core regulators.

Molecular mechanisms driving clonal aggressiveness in leukemia are not fully understood. We tracked and analyzed two mouse MLL-rearranged leukemic clones independently evolving towards higher aggressiveness. More aggressive subclones lost their growth differential ex vivo but restored it upon secondary transplantation, suggesting molecular memory of aggressiveness. Development of aggressiveness was associated with clone-specific gradual modulation of chromatin states and expression levels across the genome, with a surprising preferential trend of reversing the earlier changes between normal and leukemic progenitors. To focus on the core aggressiveness program, we identified genes with consistent changes of expression and chromatin marks that were maintained in vivo and ex vivo in both clones. Overexpressing selected core genes (Smad1 as aggressiveness driver, Irx5 and Plag1 as suppressors) affected leukemic progenitor growth in the predicted way and had convergent downstream effects on central transcription factors and repressive epigenetic modifiers, suggesting a broader regulatory network of leukemic aggressiveness.