Convergence of oncogenic cooperation at single-cell and single-gene levels drives leukemic transformation

Preleukemic single-cell landscapes reveal mutation-specific mechanisms and gene programs predictive of AML patient outcomes

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies.

Single-cell transcriptome analyses reveal critical roles of RNA splicing during leukemia progression

Leukemogenesis is proposed to be a multistep process by which normal hematopoietic stem and progenitor cells are transformed into full-blown leukemic cells, the details of which are not fully understood. Here, we performed serial single-cell transcriptome analyses of preleukemic and leukemic cells (PLCs) and constructed the cellular and molecular transformation trajectory in a Myc-driven acute myeloid leukemia (AML) model in mice, which represented the transformation course in patients. We found that the Myc targets were gradually up-regulated along the trajectory. Among them were splicing factors, which showed stage-specific prognosis for AML patients. Furthermore, we dissected the detailed gene network of a tipping point for hematopoietic stem and progenitor cells (HSPCs) to generate initiating PLCs, which was characterized by dramatically increased splicing factors and unusual RNA velocity. In the late stage, PLCs acquired explosive heterogeneity through RNA alternative splicing. Among them, the Hsp90aa1hi subpopulation was conserved in both human and mouse AML and associated with poor prognosis. Exon 4 skipping of Tmem134 was identified in these cells. While the exon skipping product Tmem134β promoted the cell cycle, full-length Tmem134α delayed tumorigenesis. Our study emphasized the critical roles of RNA splicing in the full process of leukemogenesis.

Layered immunity and layered leukemogenicity: Developmentally restricted mechanisms of pediatric leukemia initiation

Hematopoietic stem cells (HSCs) and multipotent progenitor cells (MPPs) arise in successive waves during ontogeny, and their properties change significantly throughout life. Ontological changes in HSCs/MPPs underlie corresponding changes in mechanisms of pediatric leukemia initiation. As HSCs and MPPs progress from fetal to neonatal, juvenile and adult stages of life, they undergo transcriptional and epigenetic reprogramming that modifies immune output to meet age-specific pathogenic challenges. Some immune cells arise exclusively from fetal HSCs/MPPs. We propose that this layered immunity instructs cell fates that underlie a parallel layered leukemogenicity. Indeed, some pediatric leukemias, such as juvenile myelomonocytic leukemia, myeloid leukemia of Down syndrome, and infant pre-B-cell acute lymphoblastic leukemia, are age-restricted. They only present during infancy or early childhood. These leukemias likely arise from fetal progenitors that lose competence for transformation as they age. Other childhood leukemias, such as non-infant pre-B-cell acute lymphoblastic leukemia and acute myeloid leukemia, have mutation profiles that are common in childhood but rare in morphologically similar adult leukemias. These differences could reflect temporal changes in mechanisms of mutagenesis or changes in how progenitors respond to a given mutation at different ages. Interactions between leukemogenic mutations and normal developmental switches offer potential targets for therapy.

Loss of Ezh2 function remodels the DNA replication initiation landscape

In metazoan cells, DNA replication initiates from thousands of genomic loci scattered throughout the genome called DNA replication origins. Origins are strongly associated with euchromatin, particularly open genomic regions such as promoters and enhancers. However, over a third of transcriptionally silent genes are associated with DNA replication initiation. Most of these genes are bound and repressed by the Polycomb repressive complex-2 (PRC2) through the repressive H3K27me3 mark. This is the strongest overlap observed for a chromatin regulator with replication origin activity. Here, we asked whether Polycomb-mediated gene repression is functionally involved in recruiting DNA replication origins to transcriptionally silent genes. We show that the absence of EZH2, the catalytic subunit of PRC2, results in increased DNA replication initiation, specifically in the vicinity of EZH2 binding sites. The increase in DNA replication initiation does not correlate with transcriptional de-repression or the acquisition of activating histone marks but does correlate with loss of H3K27me3 from bivalent promoters.

FLT3ITD drives context-specific changes in cell identity and variable interferon dependence during AML initiation

Acute myeloid leukemia (AML) initiation requires multiple rate-limiting mutations to cooperatively reprogram progenitor cell identity. For example, FLT3 internal tandem duplication (FLT3ITD) mutations cooperate with a variety of different initiating mutations to reprogram myeloid progenitor fate. These initiating mutations often skew toward either pediatric or adult AML patient populations, though FLT3ITD itself occurs at similar frequencies in both age groups. This raises the question of whether FLT3ITD might induce distinct transcriptional programs and unmask distinct therapeutic vulnerabilities when paired with pediatric, as opposed to adult AML-initiating mutations. To explore this possibility, we compared AML evolution in mice that carried Flt3ITD/NUP98-HOXD13 (NHD13) or Flt3ITD/Runx1DEL mutation pairs, which are respectively most common in pediatric and adult AML. Single-cell analyses and epigenome profiling revealed distinct interactions between Flt3ITD and its cooperating mutations. Whereas Flt3ITD and Flt3ITD/Runx1DEL caused aberrant expansion of myeloid progenitors, Flt3ITD/NHD13 drove the emergence of a pre-AML population that did not resemble normal hematopoietic progenitors. Differences between Flt3ITD/Runx1DEL and Flt3ITD/NHD13 cooperative target gene expression extended to fully transformed AML as well. Flt3ITD/NHD13 cooperative target genes were enriched in human NUP98-translocated AML. Flt3ITD/NHD13 selectively hijacked type I interferon signaling to drive expansion of the pre-AML population. Blocking interferon signaling delayed AML initiation and extended survival. Thus, common AML driver mutations, such as FLT3ITD, can coopt different mechanisms of transformation in different genetic contexts. Furthermore, pediatric-biased NUP98 fusions convey actionable interferon dependence.

Single-cell multi-omics of human clonal hematopoiesis reveals that DNMT3A R882 mutations perturb early progenitor states through selective hypomethylation

Somatic mutations in cancer genes have been detected in clonal expansions across healthy human tissue, including in clonal hematopoiesis. However, because mutated and wild-type cells are admixed, we have limited ability to link genotypes with phenotypes. To overcome this limitation, we leveraged multi-modality single-cell sequencing, capturing genotype, transcriptomes and methylomes in progenitors from individuals with DNMT3A R882 mutated clonal hematopoiesis. DNMT3A mutations result in myeloid over lymphoid bias, and an expansion of immature myeloid progenitors primed toward megakaryocytic-erythroid fate, with dysregulated expression of lineage and leukemia stem cell markers. Mutated DNMT3A leads to preferential hypomethylation of polycomb repressive complex 2 targets and a specific CpG flanking motif. Notably, the hypomethylation motif is enriched in binding motifs of key hematopoietic transcription factors, serving as a potential mechanistic link between DNMT3A mutations and aberrant transcriptional phenotypes. Thus, single-cell multi-omics paves the road to defining the downstream consequences of mutations that drive clonal mosaicism.