An IL-1ß-driven neutrophil-stromal cell axis fosters a BAFF-rich protumor microenvironment in individuals with multiple myeloma.

Human bone marrow permanently harbors high numbers of neutrophils, and a tumor-supportive bias of these cells could significantly impact bone marrow-confined malignancies. In individuals with multiple myeloma, the bone marrow is characterized by inflammatory stromal cells with the potential to influence neutrophils. We investigated myeloma-associated alterations in human marrow neutrophils and the impact of stromal inflammation on neutrophil function. Mature neutrophils in myeloma marrow are activated and tumor supportive and transcribe increased levels of IL1B and myeloma cell survival factor TNFSF13B (BAFF). Interactions with inflammatory stromal cells induce neutrophil activation, including BAFF secretion, in a STAT3-dependent manner, and once activated, neutrophils gain the ability to reciprocally induce stromal activation. After first-line myeloid-depleting antimyeloma treatment, human bone marrow retains residual stromal inflammation, and newly formed neutrophils are reactivated. Combined, we identify a neutrophil-stromal cell feed-forward loop driving tumor-supportive inflammation that persists after treatment and warrants novel strategies to target both stromal and immune microenvironments in multiple myeloma.

The multiple myeloma microenvironment is defined by an inflammatory stromal cell landscape.

Progression and persistence of malignancies are influenced by the local tumor microenvironment, and future eradication of currently incurable tumors will, in part, hinge on our understanding of malignant cell biology in the context of their nourishing surroundings. Here, we generated paired single-cell transcriptomic datasets of tumor cells and the bone marrow immune and stromal microenvironment in multiple myeloma. These analyses identified myeloma-specific inflammatory mesenchymal stromal cells, which spatially colocalized with tumor cells and immune cells and transcribed genes involved in tumor survival and immune modulation. Inflammatory stromal cell signatures were driven by stimulation with proinflammatory cytokines, and analyses of immune cell subsets suggested interferon-responsive effector T cell and CD8+ stem cell memory T cell populations as potential sources of stromal cell-activating cytokines. Tracking stromal inflammation in individuals over time revealed that successful antitumor induction therapy is unable to revert bone marrow inflammation, predicting a role for mesenchymal stromal cells in disease persistence.

A single oncogenic enhancer rearrangement causes concomitant EVI1 and GATA2 deregulation in leukemia.

Chromosomal rearrangements without gene fusions have been implicated in leukemogenesis by causing deregulation of proto-oncogenes via relocation of cryptic regulatory DNA elements. AML with inv(3)/t(3;3) is associated with aberrant expression of the stem-cell regulator EVI1. Applying functional genomics and genome-engineering, we demonstrate that both 3q rearrangements reposition a distal GATA2 enhancer to ectopically activate EVI1 and simultaneously confer GATA2 functional haploinsufficiency, previously identified as the cause of sporadic familial AML/MDS and MonoMac/Emberger syndromes. Genomic excision of the ectopic enhancer restored EVI1 silencing and led to growth inhibition and differentiation of AML cells, which could be replicated by pharmacologic BET inhibition. Our data show that structural rearrangements involving the chromosomal repositioning of a single enhancer can cause deregulation of two unrelated distal genes, with cancer as the outcome.

Spatial multi-omic map of human myocardial infarction.

Myocardial infarction is a leading cause of death worldwide1. Although advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality2. Here we generate an integrative high-resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from patients with myocardial infarction and controls. Multi-modal data integration enabled us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell states of major cell types and analysed them in their spatial context, evaluating their dependency on other cell types. Our data elucidate the molecular principles of human myocardial tissue organization, recapitulating a gradual cardiomyocyte and myeloid continuum following ischaemic injury. In sum, our study provides an integrative molecular map of human myocardial infarction, represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Decoding myofibroblast origins in human kidney fibrosis.

Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist1-3. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood1,2,4. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.

HNRNPC haploinsufficiency affects alternative splicing of intellectual disability-associated genes and causes a neurodevelopmental disorder.

Heterogeneous nuclear ribonucleoprotein C (HNRNPC) is an essential, ubiquitously abundant protein involved in mRNA processing. Genetic variants in other members of the HNRNP family have been associated with neurodevelopmental disorders. Here, we describe 13 individuals with global developmental delay, intellectual disability, behavioral abnormalities, and subtle facial dysmorphology with heterozygous HNRNPC germline variants. Five of them bear an identical in-frame deletion of nine amino acids in the extreme C terminus. To study the effect of this recurrent variant as well as HNRNPC haploinsufficiency, we used induced pluripotent stem cells (iPSCs) and fibroblasts obtained from affected individuals. While protein localization and oligomerization were unaffected by the recurrent C-terminal deletion variant, total HNRNPC levels were decreased. Previously, reduced HNRNPC levels have been associated with changes in alternative splicing. Therefore, we performed a meta-analysis on published RNA-seq datasets of three different cell lines to identify a ubiquitous HNRNPC-dependent signature of alternative spliced exons. The identified signature was not only confirmed in fibroblasts obtained from an affected individual but also showed a significant enrichment for genes associated with intellectual disability. Hence, we assessed the effect of decreased and increased levels of HNRNPC on neuronal arborization and neuronal migration and found that either condition affects neuronal function. Taken together, our data indicate that HNRNPC haploinsufficiency affects alternative splicing of multiple intellectual disability-associated genes and that the developing brain is sensitive to aberrant levels of HNRNPC. Hence, our data strongly support the inclusion of HNRNPC to the family of HNRNP-related neurodevelopmental disorders.

GATA2 heterozygosity causes an epigenetic feedback mechanism resulting in myeloid and erythroid dysplasia.

The transcription factor GATA2 has a pivotal role in haematopoiesis. Heterozygous germline GATA2 mutations result in a syndrome characterized by immunodeficiency, bone marrow failure and predispositions to myelodysplastic syndrome (MDS) and acute myeloid leukaemia. Clinical symptoms in these patients are diverse and mechanisms driving GATA2-related phenotypes are largely unknown. To explore the impact of GATA2 haploinsufficiency on haematopoiesis, we generated a zebrafish model carrying a heterozygous mutation of gata2b (gata2b+/-), an orthologue of GATA2. Morphological analysis revealed myeloid and erythroid dysplasia in gata2b+/- kidney marrow. Because Gata2b could affect both transcription and chromatin accessibility during lineage differentiation, this was assessed by single-cell (sc) RNA-seq and single-nucleus (sn) ATAC-seq. Sn-ATAC-seq showed that the co-accessibility between the transcription start site (TSS) and a -3.5-4.1 kb putative enhancer was more robust in gata2b+/- zebrafish HSPCs compared to wild type, increasing gata2b expression and resulting in higher genome-wide Gata2b motif use in HSPCs. As a result of increased accessibility of the gata2b locus, gata2b+/- chromatin was also more accessible during lineage differentiation. scRNA-seq data revealed myeloid differentiation defects, that is, impaired cell cycle progression, reduced expression of cebpa and cebpb and increased signatures of ribosome biogenesis. These data also revealed a differentiation delay in erythroid progenitors, aberrant proliferative signatures and down-regulation of Gata1a, a master regulator of erythropoiesis, which worsened with age. These findings suggest that cell-intrinsic compensatory mechanisms, needed to obtain normal levels of Gata2b in heterozygous HSPCs to maintain their integrity, result in aberrant lineage differentiation, thereby representing a critical step in the predisposition to MDS.

Truncated CSF3 receptors induce pro-inflammatory responses in severe congenital neutropenia.

Severe congenital neutropenia (SCN) patients are prone to develop myelodysplastic syndrome (MDS) or acute myeloid leukaemia (AML). Leukaemic progression of SCN is associated with the early acquisition of CSF3R mutations in haematopoietic progenitor cells (HPCs), which truncate the colony-stimulating factor 3 receptor (CSF3R). These mutant clones may arise years before MDS/AML becomes overt. Introduction and activation of CSF3R truncation mutants in normal HPCs causes a clonally dominant myeloproliferative state in mice treated with CSF3. Paradoxically, in SCN patients receiving CSF3 therapy, clonal dominance of CSF3R mutant clones usually occurs only after the acquisition of additional mutations shortly before frank MDS or AML is diagnosed. To seek an explanation for this discrepancy, we introduced a patient-derived CSF3R-truncating mutation in ELANE-SCN and HAX1-SCN derived and control induced pluripotent stem cells and compared the CSF3 responses of HPCs generated from these lines. In contrast to CSF3R-mutant control HPCs, CSF3R-mutant HPCs from SCN patients do not show increased proliferation but display elevated levels of inflammatory signalling. Thus, activation of the truncated CSF3R in SCN-HPCs does not evoke clonal outgrowth but causes a sustained pro-inflammatory state, which has ramifications for how these CSF3R mutants contribute to the leukaemic transformation of SCN.

Rapid in vitro generation of bona fide exhausted CD8+ T cells is accompanied by Tcf7 promotor methylation.

Exhaustion is a dysfunctional state of cytotoxic CD8+ T cells (CTL) observed in chronic infection and cancer. Current in vivo models of CTL exhaustion using chronic viral infections or cancer yield very few exhausted CTL, limiting the analysis that can be done on these cells. Establishing an in vitro system that rapidly induces CTL exhaustion would therefore greatly facilitate the study of this phenotype, identify the truly exhaustion-associated changes and allow the testing of novel approaches to reverse or prevent exhaustion. Here we show that repeat stimulation of purified TCR transgenic OT-I CTL with their specific peptide induces all the functional (reduced cytokine production and polyfunctionality, decreased in vivo expansion capacity) and phenotypic (increased inhibitory receptors expression and transcription factor changes) characteristics of exhaustion. Importantly, in vitro exhausted cells shared the transcriptomic characteristics of the gold standard of exhaustion, CTL from LCMV cl13 infections. Gene expression of both in vitro and in vivo exhausted CTL was distinct from T cell anergy. Using this system, we show that Tcf7 promoter DNA methylation contributes to TCF1 downregulation in exhausted CTL. Thus this novel in vitro system can be used to identify genes and signaling pathways involved in exhaustion and will facilitate the screening of reagents that prevent/reverse CTL exhaustion.

Increased CXCL4 expression in hematopoietic cells links inflammation and progression of bone marrow fibrosis in MPN.

Primary myelofibrosis (PMF) is a myeloproliferative neoplasm (MPN) that leads to progressive bone marrow (BM) fibrosis. Although the cellular mutations involved in the pathogenesis of PMF have been extensively investigated, the sequential events that drive stromal activation and fibrosis by hematopoietic-stromal cross-talk remain elusive. Using an unbiased approach and validation in patients with MPN, we determined that the differential spatial expression of the chemokine CXCL4/platelet factor-4 marks the progression of fibrosis. We show that the absence of hematopoietic CXCL4 ameliorates the MPN phenotype, reduces stromal cell activation and BM fibrosis, and decreases the activation of profibrotic pathways in megakaryocytes, inflammation in fibrosis-driving cells, and JAK/STAT activation in both megakaryocytes and stromal cells in 3 murine PMF models. Our data indicate that higher CXCL4 expression in MPN has profibrotic effects and is a mediator of the characteristic inflammation. Therefore, targeting CXCL4 might be a promising strategy to reduce inflammation in PMF.

Identification of osteolineage cell-derived extracellular vesicle cargo implicated in hematopoietic support.

Osteolineage cell-derived extracellular vesicles (EVs) play a regulatory role in hematopoiesis and have been shown to promote the ex vivo expansion of human hematopoietic stem and progenitor cells (HSPCs). Here, we demonstrate that EVs from different human osteolineage sources do not have the same HSPC expansion promoting potential. Comparison of stimulatory and non-stimulatory osteolineage EVs by next-generation sequencing and mass spectrometry analyses revealed distinct microRNA and protein signatures identifying EV-derived candidate regulators of ex vivo HSPC expansion. Accordingly, the treatment of umbilical cord blood-derived CD34+ HSPCs with stimulatory EVs-altered HSPC transcriptome, including genes with known roles in cell proliferation. An integrative bioinformatics approach, which connects the HSPC gene expression data with the candidate cargo in stimulatory EVs, delineated the potentially targeted biological functions and pathways during hematopoietic cell expansion and development. In conclusion, our study gives novel insights into the complex biological role of EVs in osteolineage cell-HSPC crosstalk and promotes the utility of EVs and their cargo as therapeutic agents in regenerative medicine.

Expression of Plet1 controls interstitial migration of murine small intestinal dendritic cells.

Under homeostatic conditions, dendritic cells (DCs) continuously patrol the intestinal lamina propria. Upon antigen encounter, DCs initiate C-C motif chemokine receptor 7 (CCR7) expression and migrate into lymph nodes to direct T cell activation and differentiation. The mechanistic underpinnings of DC migration from the tissues to lymph nodes have been largely elucidated, contributing greatly to our understanding of DC functionality and intestinal immunity. In contrast, the molecular mechanisms allowing DCs to efficiently migrate through the complex extracellular matrix of the intestinal lamina propria prior to antigen encounter are still incompletely understood. Here we show that small intestinal murine CD11b+ CD103+ DCs express Placenta-expressed transcript 1 (Plet1), a glycophoshatidylinositol (GPI)-anchored surface protein involved in migration of keratinocytes during wound healing. In the absence of Plet1, CD11b+ CD103+ DCs display aberrant migratory behavior, and accumulate in the small intestine, independent of CCR7 responsiveness. RNA-sequencing indicated involvement of Plet1 in extracellular matrix-interactiveness, and subsequent in-vitro migration assays revealed that Plet1 augments the ability of DCs to migrate through extracellular matrix containing environments. In conclusion, our findings reveal that expression of Plet1 facilitates homeostatic interstitial migration of small intestinal DCs.

MBD4 guards against methylation damage and germ line deficiency predisposes to clonal hematopoiesis and early-onset AML.

The tendency of 5-methylcytosine (5mC) to undergo spontaneous deamination has had a major role in shaping the human genome, and this methylation damage remains the primary source of somatic mutations that accumulate with age. How 5mC deamination contributes to cancer risk in different tissues remains unclear. Genomic profiling of 3 early-onset acute myeloid leukemias (AMLs) identified germ line loss of MBD4 as an initiator of 5mC-dependent hypermutation. MBD4-deficient AMLs display a 33-fold higher mutation burden than AML generally, with >95% being C>T in the context of a CG dinucleotide. This distinctive signature was also observed in sporadic cancers that acquired biallelic mutations in MBD4 and in Mbd4 knockout mice. Sequential sampling of germ line cases demonstrated repeated expansion of blood cell progenitors with pathogenic mutations in DNMT3A, a key driver gene for both clonal hematopoiesis and AML. Our findings reveal genetic and epigenetic factors that shape the mutagenic influence of 5mC. Within blood cells, this links methylation damage to the driver landscape of clonal hematopoiesis and reveals a conserved path to leukemia. Germ line MBD4 deficiency enhances cancer susceptibility and predisposes to AML.

An autonomous CEBPA enhancer specific for myeloid-lineage priming and neutrophilic differentiation.

Neutrophilic differentiation is dependent on CCAAT enhancer-binding protein α (C/EBPα), a transcription factor expressed in multiple organs including the bone marrow. Using functional genomic technologies in combination with clustered regularly-interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 genome editing and in vivo mouse modeling, we show that CEBPA is located in a 170-kb topological-associated domain that contains 14 potential enhancers. Of these, 1 enhancer located +42 kb from CEBPA is active and engages with the CEBPA promoter in myeloid cells only. Germ line deletion of the homologous enhancer in mice in vivo reduces Cebpa levels exclusively in hematopoietic stem cells (HSCs) and myeloid-primed progenitor cells leading to severe defects in the granulocytic lineage, without affecting any other Cebpa-expressing organ studied. The enhancer-deleted progenitor cells lose their myeloid transcription program and are blocked in differentiation. Deletion of the enhancer also causes loss of HSC maintenance. We conclude that a single +42-kb enhancer is essential for CEBPA expression in myeloid cells only.

Mutational spectrum of myeloid malignancies with inv(3)/t(3;3) reveals a predominant involvement of RAS/RTK signaling pathways.

Myeloid malignancies bearing chromosomal inv(3)/t(3;3) abnormalities are among the most therapy-resistant leukemias. Deregulated expression of EVI1 is the molecular hallmark of this disease; however, the genome-wide spectrum of cooperating mutations in this disease subset has not been systematically elucidated. Here, we show that 98% of inv(3)/t(3;3) myeloid malignancies harbor mutations in genes activating RAS/receptor tyrosine kinase (RTK) signaling pathways. In addition, hemizygous mutations in GATA2, as well as heterozygous alterations in RUNX1, SF3B1, and genes encoding epigenetic modifiers, frequently co-occur with the inv(3)/t(3;3) aberration. Notably, neither mutational patterns nor gene expression profiles differ across inv(3)/t(3;3) acute myeloid leukemia, chronic myeloid leukemia, and myelodysplastic syndrome cases, suggesting recognition of inv(3)/t(3;3) myeloid malignancies as a single disease entity irrespective of blast count. The high incidence of activating RAS/RTK signaling mutations may provide a target for a rational treatment strategy in this high-risk patient group.

ICL-induced miR139-3p and miR199a-3p have opposite roles in hematopoietic cell expansion and leukemic transformation.

Interstrand crosslinks (ICLs) are toxic DNA lesions that cause severe genomic damage during replication, especially in Fanconi anemia pathway-deficient cells. This results in progressive bone marrow failure and predisposes to acute myeloid leukemia (AML). The molecular mechanisms responsible for these defects are largely unknown. Using Ercc1-deficient mice, we show that Trp53 is responsible for ICL-induced bone marrow failure and that loss of Trp53 is leukemogenic in this model. In addition, Ercc1-deficient myeloid progenitors gain elevated levels of miR-139-3p and miR-199a-3p with age. These microRNAs exert opposite effects on hematopoiesis. Ectopic expression of miR-139-3p strongly inhibited proliferation of myeloid progenitors, whereas inhibition of miR-139-3p activity restored defective proliferation of Ercc1-deficient progenitors. Conversely, the inhibition of miR-199a-3p functions aggravated the myeloid proliferation defect in the Ercc1-deficient model, whereas its enforced expression enhanced proliferation of progenitors. Importantly, miR-199a-3p caused AML in a pre-leukemic mouse model, supporting its role as an onco-microRNA. Target genes include HuR for miR-139-3p and Prdx6, Runx1, and Suz12 for miR-199a-3p. The latter genes have previously been implicated as tumor suppressors in de novo and secondary AML. These findings show that, in addition to TRP53-controlled mechanisms, miR-139-3p and miR-199a-3p are involved in the defective hematopoietic function of ICL-repair deficient myeloid progenitors.

Source tracking of Enterococcus moraviensis and E. haemoperoxidus.

Enterococci were detected occasionally in 100 L samples of water abstracted from a shallow aquifer in a natural dune infiltration area for drinking water production. Enterococcus moraviensis was the species most frequently identified in these samples. Because there are no existing reports of faecal sources of E. moraviensis and the closely related E. hemoperoxidus, this study aimed to find such sources of these two species in the dunes. Faecal samples from various animal species living in the vicinity of abstraction wells, were analysed for enterococci on Slanetz and Bartley Agar. From these samples, enterococci isolates (1,386 in total) were subsequently identified using matrix assisted laser desorption ionization-time of flight (MALDI-TOF) mass spectrometry. E. moraviensis was found in the faeces of geese, foxes and rabbits. Also, E. haemoperoxidus was isolated from goose faeces. Using hierarchical clustering, the species composition of Enterococcus spp. isolated from abstracted water formed one cluster with the species composition found in geese droppings. A sanitary survey supported the indication that feral geese may provide a substantial faecal load in particular parts of this dune infiltration area, close to the water abstraction system. This study confirms the faecal origin of E. moraviensis and E. haemoperoxidus from specific animals, which strengthens their significance as faecal indicators.

Two splice-factor mutant leukemia subgroups uncovered at the boundaries of MDS and AML using combined gene expression and DNA-methylation profiling.

Mutations in splice factor (SF) genes occur more frequently in myelodysplastic syndromes (MDS) than in acute myeloid leukemias (AML). We sequenced complementary DNA from bone marrow of 47 refractory anemia with excess blasts (RAEB) patients, 29 AML cases with low marrow blast cell count, and 325 other AML patients and determined the presence of SF-hotspot mutations in SF3B1, U2AF35, and SRSF2. SF mutations were found in 10 RAEB, 12 AML cases with low marrow blast cell count, and 25 other AML cases. Our study provides evidence that SF-mutant RAEB and SF-mutant AML are clinically, cytologically, and molecularly highly similar. An integrated analysis of genomewide messenger RNA (mRNA) expression profiling and DNA-methylation profiling data revealed 2 unique patient clusters highly enriched for SF-mutant RAEB/AML. The combined genomewide mRNA expression profiling/DNA-methylation profiling signatures revealed 1 SF-mutant patient cluster with an erythroid signature. The other SF-mutant patient cluster was enriched for NRAS/KRAS mutations and showed an inferior survival. We conclude that SF-mutant RAEB/AML constitutes a related disorder overriding the artificial separation between AML and MDS, and that SF-mutant RAEB/AML is composed of 2 molecularly and clinically distinct subgroups. We conclude that SF-mutant disorders should be considered as myeloid malignancies that transcend the boundaries of AML and MDS.

Cooperativity of RUNX1 and CSF3R mutations in severe congenital neutropenia: a unique pathway in myeloid leukemogenesis.

Severe congenital neutropenia (CN) is a preleukemic bone marrow failure syndrome with a 20% risk of evolving into leukemia or myelodysplastic syndrome (MDS). Patterns of acquisition of leukemia-associated mutations were investigated using next-generation deep-sequencing in 31 CN patients who developed leukemia or MDS. Twenty (64.5%) of the 31 patients had mutations in RUNX1. A majority of patients with RUNX1 mutations (80.5%) also had acquired CSF3R mutations. In contrast to their high frequency in CN patients who developed leukemia or MDS, RUNX1 mutations were found in only 9 of 307 (2.9%) patients with de novo pediatric acute myeloid leukemia. A sequential analysis at stages prior to overt leukemia revealed RUNX1 mutations to be late events in leukemic transformation. Single-cell analyses in 2 patients showed that RUNX1 and CSF3R mutations were present in the same malignant clone. Functional studies demonstrated elevated granulocyte colony-stimulating factor (G-CSF)-induced proliferation with diminished myeloid differentiation of hematopoietic CD34(+) cells coexpressing mutated forms of RUNX1 and CSF3R. The high frequency of cooperating RUNX1 and CSF3R mutations in CN patients suggests a novel molecular pathway of leukemogenesis: mutations in the hematopoietic cytokine receptor (G-CSFR) in combination with the second mutations in the downstream hematopoietic transcription fator (RUNX1). The detection of both RUNX1 and CSF3R mutations could be used as a marker for identifying CN patients with a high risk of progressing to leukemia or MDS.