Increased CXCL10 (IP-10) associates with advanced myeloproliferative neoplasms and its loss dampens erythrocytosis in mouse models.

Key studies in pre-leukemic disorders have linked increases in pro-inflammatory cytokines with accelerated phases of disease, but the precise role of the cellular microenvironment in disease initiation and evolution remains poorly understood. In myeloproliferative neoplasms (MPNs), higher levels of specific cytokines have been previously correlated with increased disease severity (TNF-α, IP-10) and decreased survival (IL-8). Whereas TNF-α and IL-8 have been studied by numerous groups, there is a relative paucity of studies on IP-10 (CXCL10). Here we explore the relationship of IP-10 levels with detailed genomic and clinical data and undertake a complementary cytokine screen alongside functional assays in a wide range of MPN mouse models. Similar to patients, levels of IP-10 were increased in mice with more severe disease phenotypes (e.g., JAK2V617F/V617F TET2-/- double mutant mice) compared to those with less severe phenotypes (e.g., CALRdel52 or JAK2+/V617F mice) and WT littermate controls. While exposure to IP-10 did not directly alter proliferation or survival in single hematopoietic stem cells (HSCs) in vitro, IP-10-/- mice transplanted with disease initiating HSCs developed an MPN phenotype more slowly, suggesting that the effect of IP-10 loss was non-cell autonomous. To explore the broader effects of IP-10 loss, we crossed IP-10-/- mice into a series of MPN mouse models and show that its loss reduces the erythrocytosis observed in mice with the most severe phenotype. Together these data point to a potential role for blocking IP-10 activity in the management of MPNs.

Inherited polygenic effects on common hematological traits influence clonal selection on JAK2V617F and the development of myeloproliferative neoplasms.

Myeloproliferative neoplasms (MPNs) are chronic cancers characterized by overproduction of mature blood cells. Their causative somatic mutations, for example, JAK2V617F, are common in the population, yet only a minority of carriers develop MPN. Here we show that the inherited polygenic loci that underlie common hematological traits influence JAK2V617F clonal expansion. We identify polygenic risk scores (PGSs) for monocyte count and plateletcrit as new risk factors for JAK2V617F positivity. PGSs for several hematological traits influenced the risk of different MPN subtypes, with low PGSs for two platelet traits also showing protective effects in JAK2V617F carriers, making them two to three times less likely to have essential thrombocythemia than carriers with high PGSs. We observed that extreme hematological PGSs may contribute to an MPN diagnosis in the absence of somatic driver mutations. Our study showcases how polygenic backgrounds underlying common hematological traits influence both clonal selection on somatic mutations and the subsequent phenotype of cancer.

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.

Different niches for stem cells carrying the same oncogenic driver affect pathogenesis and therapy response in myeloproliferative neoplasms.

Aging facilitates the expansion of hematopoietic stem cells (HSCs) carrying clonal hematopoiesis-related somatic mutations and the development of myeloid malignancies, such as myeloproliferative neoplasms (MPNs). While cooperating mutations can cause transformation, it is unclear whether distinct bone marrow (BM) HSC-niches can influence the growth and therapy response of HSCs carrying the same oncogenic driver. Here we found different BM niches for HSCs in MPN subtypes. JAK-STAT signaling differentially regulates CDC42-dependent HSC polarity, niche interaction and mutant cell expansion. Asymmetric HSC distribution causes differential BM niche remodeling: sinusoidal dilation in polycythemia vera and endosteal niche expansion in essential thrombocythemia. MPN development accelerates in a prematurely aged BM microenvironment, suggesting that the specialized niche can modulate mutant cell expansion. Finally, dissimilar HSC-niche interactions underpin variable clinical response to JAK inhibitor. Therefore, HSC-niche interactions influence the expansion rate and therapy response of cells carrying the same clonal hematopoiesis oncogenic driver.

Iron homeostasis governs erythroid phenotype in polycythemia vera.

Polycythemia vera (PV) is a myeloproliferative neoplasm driven by activating mutations in JAK2 that result in unrestrained erythrocyte production, increasing patients' hematocrit and hemoglobin concentrations, placing them at risk of life-threatening thrombotic events. Our genome-wide association study of 440 PV cases and 403 351 controls using UK Biobank data showed that single nucleotide polymorphisms in HFE known to cause hemochromatosis are highly associated with PV diagnosis, linking iron regulation to PV. Analysis of the FinnGen dataset independently confirmed overrepresentation of homozygous HFE variants in patients with PV. HFE influences the expression of hepcidin, the master regulator of systemic iron homeostasis. Through genetic dissection of mouse models of PV, we show that the PV erythroid phenotype is directly linked to hepcidin expression: endogenous hepcidin upregulation alleviates erythroid disease whereas hepcidin ablation worsens it. Furthermore, we demonstrate that in PV, hepcidin is not regulated by expanded erythropoiesis but is likely governed by inflammatory cytokines signaling via GP130-coupled receptors. These findings have important implications for understanding the pathophysiology of PV and offer new therapeutic strategies for this disease.

Diverse mutational landscapes in human lymphocytes.

The lymphocyte genome is prone to many threats, including programmed mutation during differentiation1, antigen-driven proliferation and residency in diverse microenvironments. Here, after developing protocols for expansion of single-cell lymphocyte cultures, we sequenced whole genomes from 717 normal naive and memory B and T cells and haematopoietic stem cells. All lymphocyte subsets carried more point mutations and structural variants than haematopoietic stem cells, with higher burdens in memory cells than in naive cells, and with T cells accumulating mutations at a higher rate throughout life. Off-target effects of immunological diversification accounted for approximately half of the additional differentiation-associated mutations in lymphocytes. Memory B cells acquired, on average, 18 off-target mutations genome-wide for every on-target IGHV mutation during the germinal centre reaction. Structural variation was 16-fold higher in lymphocytes than in stem cells, with around 15% of deletions being attributable to off-target recombinase-activating gene activity. DNA damage from ultraviolet light exposure and other sporadic mutational processes generated hundreds to thousands of mutations in some memory cells. The mutation burden and signatures of normal B cells were broadly similar to those seen in many B-cell cancers, suggesting that malignant transformation of lymphocytes arises from the same mutational processes that are active across normal ontogeny. The mutational landscape of normal lymphocytes chronicles the off-target effects of programmed genome engineering during immunological diversification and the consequences of differentiation, proliferation and residency in diverse microenvironments.

Clonal dynamics of haematopoiesis across the human lifespan.

Age-related change in human haematopoiesis causes reduced regenerative capacity1, cytopenias2, immune dysfunction3 and increased risk of blood cancer4-6, but the reason for such abrupt functional decline after 70 years of age remains unclear. Here we sequenced 3,579 genomes from single cell-derived colonies of haematopoietic cells across 10 human subjects from 0 to 81 years of age. Haematopoietic stem cells or multipotent progenitors (HSC/MPPs) accumulated a mean of 17 mutations per year after birth and lost 30 base pairs per year of telomere length. Haematopoiesis in adults less than 65 years of age was massively polyclonal, with high clonal diversity and a stable population of 20,000-200,000 HSC/MPPs contributing evenly to blood production. By contrast, haematopoiesis in individuals aged over 75 showed profoundly decreased clonal diversity. In each of the older subjects, 30-60% of haematopoiesis was accounted for by 12-18 independent clones, each contributing 1-34% of blood production. Most clones had begun their expansion before the subject was 40 years old, but only 22% had known driver mutations. Genome-wide selection analysis estimated that between 1 in 34 and 1 in 12 non-synonymous mutations were drivers, accruing at constant rates throughout life, affecting more genes than identified in blood cancers. Loss of the Y chromosome conferred selective benefits in males. Simulations of haematopoiesis, with constant stem cell population size and constant acquisition of driver mutations conferring moderate fitness benefits, entirely explained the abrupt change in clonal structure in the elderly. Rapidly decreasing clonal diversity is a universal feature of haematopoiesis in aged humans, underpinned by pervasive positive selection acting on many more genes than currently identified.

STAT1 is essential for HSC function and maintains MHCIIhi stem cells that resist myeloablation and neoplastic expansion.

Adult hematopoietic stem cells (HSCs) are predominantly quiescent and can be activated in response to acute stress such as infection or cytotoxic insults. STAT1 is a pivotal downstream mediator of interferon (IFN) signaling and is required for IFN-induced HSC proliferation, but little is known about the role of STAT1 in regulating homeostatic hematopoietic stem/progenitor cells (HSPCs). Here, we show that loss of STAT1 altered the steady state HSPC landscape, impaired HSC function in transplantation assays, delayed blood cell regeneration following myeloablation, and disrupted molecular programs that protect HSCs, including control of quiescence. Our results also reveal STAT1-dependent functional HSC heterogeneity. A previously unrecognized subset of homeostatic HSCs with elevated major histocompatibility complex class II (MHCII) expression (MHCIIhi) displayed molecular features of reduced cycling and apoptosis and was refractory to 5-fluorouracil-induced myeloablation. Conversely, MHCIIlo HSCs displayed increased megakaryocytic potential and were preferentially expanded in CALR mutant mice with thrombocytosis. Similar to mice, high MHCII expression is a feature of human HSCs residing in a deeper quiescent state. Our results therefore position STAT1 at the interface of stem cell heterogeneity and the interplay between stem cells and the adaptive immune system, areas of broad interest in the wider stem cell field.

Life histories of myeloproliferative neoplasms inferred from phylogenies.

Mutations in cancer-associated genes drive tumour outgrowth, but our knowledge of the timing of driver mutations and subsequent clonal dynamics is limited1-3. Here, using whole-genome sequencing of 1,013 clonal haematopoietic colonies from 12 patients with myeloproliferative neoplasms, we identified 580,133 somatic mutations to reconstruct haematopoietic phylogenies and determine clonal histories. Driver mutations were estimated to occur early in life, including the in utero period. JAK2V617F was estimated to have been acquired by 33 weeks of gestation to 10.8 years of age in 5 patients in whom JAK2V617F was the first event. DNMT3A mutations were acquired by 8 weeks of gestation to 7.6 years of age in 4 patients, and a PPM1D mutation was acquired by 5.8 years of age. Additional genomic events occurred before or following JAK2V617F acquisition and as independent clonal expansions. Sequential driver mutation acquisition was separated by decades across life, often outcompeting ancestral clones. The mean latency between JAK2V617F acquisition and diagnosis was 30 years (range 11-54 years). Estimated historical rates of clonal expansion varied substantially (3% to 190% per year), increased with additional driver mutations, and predicted latency to diagnosis. Our study suggests that early driver mutation acquisition and life-long growth and evolution underlie adult myeloproliferative neoplasms, raising opportunities for earlier intervention and a new model for cancer development.

Zinc-dependent multimerization of mutant calreticulin is required for MPL binding and MPN pathogenesis.

Calreticulin (CALR) is mutated in the majority of JAK2/MPL-unmutated myeloproliferative neoplasms (MPNs). Mutant CALR (CALRdel52) exerts its effect by binding to the thrombopoietin receptor MPL to cause constitutive activation of JAK-STAT signaling. In this study, we performed an extensive mutagenesis screen of the CALR globular N-domain and revealed 2 motifs critical for CALRdel52 oncogenic activity: (1) the glycan-binding lectin motif and (2) the zinc-binding domain. Further analysis demonstrated that the zinc-binding domain was essential for formation of CALRdel52 multimers, which was a co-requisite for MPL binding. CALRdel52 variants incapable of binding zinc were unable to homomultimerize, form CALRdel52-MPL heteromeric complexes, or stimulate JAK-STAT signaling. Finally, treatment with zinc chelation disrupted CALRdel52-MPL complexes in hematopoietic cells in conjunction with preferential eradication of cells expressing CALRdel52 relative to cells expressing other MPN oncogenes. In addition, zinc chelators exhibited a therapeutic effect in preferentially impairing growth of CALRdel52-mutant erythroblasts relative to unmutated erythroblasts in primary cultures of MPN patients. Together, our data implicate zinc as an essential cofactor for CALRdel52 oncogenic activity by enabling CALRdel52 multimerization and interaction with MPL, and suggests that perturbation of intracellular zinc levels may represent a new approach to abrogate the oncogenic activity of CALRdel52 in the treatment of MPNs.

Hematopoietic Stem Cell Heterogeneity Is Linked to the Initiation and Therapeutic Response of Myeloproliferative Neoplasms.

The implications of stem cell heterogeneity for disease pathogenesis and therapy are poorly defined. JAK2V617F+ myeloproliferative neoplasms (MPNs), harboring the same mutation in hematopoietic stem cells (HSCs), display diverse phenotypes, including polycythemia vera (PV), essential thrombocythemia (ET), and primary myelofibrosis (PMF). These chronic malignant disorders are ideal models to analyze the pathological consequences of stem cell heterogeneity. Single-cell gene expression profiling with parallel mutation detection demonstrated that the megakaryocyte (Mk)-primed HSC subpopulation expanded significantly with enhanced potential in untreated individuals with JAK2V617F+ ET, driven primarily by the JAK2 mutation and elevated interferon signaling. During treatment, mutant HSCs were targeted preferentially in the Mk-primed HSC subpopulation. Interestingly, homozygous mutant HSCs were forced to re-enter quiescence, whereas their heterozygous counterparts underwent apoptosis. This study provides important evidence for the association of stem cell heterogeneity with the pathogenesis and therapeutic response of a malignant disease.

The stem/progenitor landscape is reshaped in a mouse model of essential thrombocythemia and causes excess megakaryocyte production.

Frameshift mutations in CALR (calreticulin) are associated with essential thrombocythemia (ET), but the stages at and mechanisms by which mutant CALR drives transformation remain incompletely defined. Here, we use single-cell approaches to examine the hematopoietic stem/progenitor cell landscape in a mouse model of mutant CALR-driven ET. We identify a trajectory linking hematopoietic stem cells (HSCs) with megakaryocytes and prospectively identify a previously unknown intermediate population that is overrepresented in the disease state. We also show that mutant CALR drives transformation primarily from the earliest stem cell compartment, with some contribution from megakaryocyte progenitors. Last, relative to wild-type HSCs, mutant CALR HSCs show increases in JAK-STAT signaling, the unfolded protein response, cell cycle, and a previously undescribed up-regulation of cholesterol biosynthesis. Overall, we have identified a novel megakaryocyte-biased cell population that is increased in a mouse model of ET and described transcriptomic changes linking CALR mutations to increased HSC proliferation and megakaryopoiesis.

Longitudinal Cytokine Profiling Identifies GRO-α and EGF as Potential Biomarkers of Disease Progression in Essential Thrombocythemia.

Myeloproliferative neoplasms (MPNs) are characterized by deregulation of mature blood cell production and increased risk of myelofibrosis (MF) and leukemic transformation. Numerous driver mutations have been identified but substantial disease heterogeneity remains unexplained, implying the involvement of additional as yet unidentified factors. The inflammatory microenvironment has recently attracted attention as a crucial factor in MPN biology, in particular whether inflammatory cytokines and chemokines contribute to disease establishment or progression. Here we present a large-scale study of serum cytokine profiles in more than 400 MPN patients and identify an essential thrombocythemia (ET)-specific inflammatory cytokine signature consisting of Eotaxin, GRO-α, and EGF. Levels of 2 of these markers (GRO-α and EGF) in ET patients were associated with disease transformation in initial sample collection (GRO-α) or longitudinal sampling (EGF). In ET patients with extensive genomic profiling data (n = 183) cytokine levels added significant prognostic value for predicting transformation from ET to MF. Furthermore, CD56+CD14+ pro-inflammatory monocytes were identified as a novel source of increased GRO-α levels. These data implicate the immune cell microenvironment as a significant player in ET disease evolution and illustrate the utility of cytokines as potential biomarkers for reaching beyond genomic classification for disease stratification and monitoring.

Mutant Calreticulin in the Myeloproliferative Neoplasms.

Mutations in the gene for calreticulin (CALR) were identified in the myeloproliferative neoplasms (MPNs) essential thrombocythaemia (ET) and primary myelofibrosis (MF) in 2013; in combination with previously described mutations in JAK2 and MPL, driver mutations have now been described for the majority of MPN patients. In subsequent years, researchers have begun to unravel the mechanisms by which mutant CALR drives transformation and to understand their clinical implications. Mutant CALR activates the thrombopoietin receptor (MPL), causing constitutive activation of Janus kinase 2 (JAK2) signaling and cytokine independent growth in vitro. Mouse models show increased numbers of hematopoietic stem cells (HSCs) and overproduction of megakaryocytic lineage cells with associated thrombocytosis. In the clinic, detection of CALR mutations has been embedded in World Health Organization and other international diagnostic guidelines. Distinct clinical and laboratory associations of CALR mutations have been identified together with their prognostic significance, with CALR mutant patients showing increased overall survival. The discovery and subsequent study of CALR mutations have illuminated novel aspects of megakaryopoiesis and raised the possibility of new therapeutic approaches.

Cohesin-dependent regulation of gene expression during differentiation is lost in cohesin-mutated myeloid malignancies.

Cohesin complex disruption alters gene expression, and cohesin mutations are common in myeloid neoplasia, suggesting a critical role in hematopoiesis. Here, we explore cohesin dynamics and regulation of hematopoietic stem cell homeostasis and differentiation. Cohesin binding increases at active regulatory elements only during erythroid differentiation. Prior binding of the repressive Ets transcription factor Etv6 predicts cohesin binding at these elements and Etv6 interacts with cohesin at chromatin. Depletion of cohesin severely impairs erythroid differentiation, particularly at Etv6-prebound loci, but augments self-renewal programs. Together with corroborative findings in acute myeloid leukemia and myelodysplastic syndrome patient samples, these data suggest cohesin-mediated alleviation of Etv6 repression is required for dynamic expression at critical erythroid genes during differentiation and how this may be perturbed in myeloid malignancies.

Notch Signaling Mediates Secondary Senescence.

Oncogene-induced senescence (OIS) is a tumor suppressive response to oncogene activation that can be transmitted to neighboring cells through secreted factors of the senescence-associated secretory phenotype (SASP). Currently, primary and secondary senescent cells are not considered functionally distinct endpoints. Using single-cell analysis, we observed two distinct transcriptional endpoints, a primary endpoint marked by Ras and a secondary endpoint marked by Notch activation. We find that secondary oncogene-induced senescence in vitro and in vivo requires Notch, rather than SASP alone, as previously thought. Moreover, Notch signaling weakens, but does not abolish, SASP in secondary senescence. Global transcriptomic differences, a blunted SASP response, and the induction of fibrillar collagens in secondary senescence point toward a functional diversification between secondary and primary senescence.

Clonal approaches to understanding the impact of mutations on hematologic disease development.

Interrogation of hematopoietic tissue at the clonal level has a rich history spanning over 50 years, and has provided critical insights into both normal and malignant hematopoiesis. Characterization of chromosomes identified some of the first genetic links to cancer with the discovery of chromosomal translocations in association with many hematological neoplasms. The unique accessibility of hematopoietic tissue and the ability to clonally expand hematopoietic progenitors in vitro has provided fundamental insights into the cellular hierarchy of normal hematopoiesis, as well as the functional impact of driver mutations in disease. Transplantation assays in murine models have enabled cellular assessment of the functional consequences of somatic mutations in vivo. Most recently, next-generation sequencing-based assays have shown great promise in allowing multi-"omic" characterization of single cells. Here, we review how clonal approaches have advanced our understanding of disease development, focusing on the acquisition of somatic mutations, clonal selection, driver mutation cooperation, and tumor evolution.

Classification and Personalized Prognosis in Myeloproliferative Neoplasms.

BACKGROUND: Myeloproliferative neoplasms, such as polycythemia vera, essential thrombocythemia, and myelofibrosis, are chronic hematologic cancers with varied progression rates. The genomic characterization of patients with myeloproliferative neoplasms offers the potential for personalized diagnosis, risk stratification, and treatment. METHODS: We sequenced coding exons from 69 myeloid cancer genes in patients with myeloproliferative neoplasms, comprehensively annotating driver mutations and copy-number changes. We developed a genomic classification for myeloproliferative neoplasms and multistage prognostic models for predicting outcomes in individual patients. Classification and prognostic models were validated in an external cohort. RESULTS: A total of 2035 patients were included in the analysis. A total of 33 genes had driver mutations in at least 5 patients, with mutations in JAK2, CALR, or MPL being the sole abnormality in 45% of the patients. The numbers of driver mutations increased with age and advanced disease. Driver mutations, germline polymorphisms, and demographic variables independently predicted whether patients received a diagnosis of essential thrombocythemia as compared with polycythemia vera or a diagnosis of chronic-phase disease as compared with myelofibrosis. We defined eight genomic subgroups that showed distinct clinical phenotypes, including blood counts, risk of leukemic transformation, and event-free survival. Integrating 63 clinical and genomic variables, we created prognostic models capable of generating personally tailored predictions of clinical outcomes in patients with chronic-phase myeloproliferative neoplasms and myelofibrosis. The predicted and observed outcomes correlated well in internal cross-validation of a training cohort and in an independent external cohort. Even within individual categories of existing prognostic schemas, our models substantially improved predictive accuracy. CONCLUSIONS: Comprehensive genomic characterization identified distinct genetic subgroups and provided a classification of myeloproliferative neoplasms on the basis of causal biologic mechanisms. Integration of genomic data with clinical variables enabled the personalized predictions of patients' outcomes and may support the treatment of patients with myeloproliferative neoplasms. (Funded by the Wellcome Trust and others.).