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.

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.

Adaptation to ex vivo culture reduces human hematopoietic stem cell activity independently of cell cycle.

Loss of long-term hematopoietic stem cell (LT-HSC) function ex vivo hampers the success of clinical protocols reliant on culture. However, the kinetics and mechanisms by which this occurs remain incompletely characterized. Here, through time-resolved scRNA-Seq, matched in vivo functional analysis and the use of a reversible in vitro system of early G1 arrest, we define the sequence of transcriptional and functional events occurring during the first ex vivo division of human LT-HSCs. We demonstrate that the sharpest loss of LT-HSC repopulation capacity happens early on, between 6 and 24 hours of culture, before LT-HSCs commit to cell cycle progression. During this time window, LT-HSCs adapt to the culture environment, limiting global variability in gene expression and transiently upregulating gene networks involved in signaling and stress responses. From 24 hours, LT-HSC progression past early G1 contributes to the establishment of differentiation programmes in culture. However, contrary to current assumptions, we demonstrate that loss of HSC function ex vivo is independent of cell cycle progression. Finally, we show that targeting LT-HSC adaptation to culture by inhibiting early activation of JAK/STAT signaling improves HSC long-term repopulating function ex vivo. Collectively, our study demonstrates that controlling early LT-HSC adaptation to ex vivo culture, for example via JAK inhibition, is of critical importance to improve HSC gene therapy and expansion protocols.

Massive genomic rearrangement acquired in a single catastrophic event during cancer development.

Cancer is driven by somatically acquired point mutations and chromosomal rearrangements, conventionally thought to accumulate gradually over time. Using next-generation sequencing, we characterize a phenomenon, which we term chromothripsis, whereby tens to hundreds of genomic rearrangements occur in a one-off cellular crisis. Rearrangements involving one or a few chromosomes crisscross back and forth across involved regions, generating frequent oscillations between two copy number states. These genomic hallmarks are highly improbable if rearrangements accumulate over time and instead imply that nearly all occur during a single cellular catastrophe. The stamp of chromothripsis can be seen in at least 2%-3% of all cancers, across many subtypes, and is present in ∼25% of bone cancers. We find that one, or indeed more than one, cancer-causing lesion can emerge out of the genomic crisis. This phenomenon has important implications for the origins of genomic remodeling and temporal emergence of cancer.

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.

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.

Population dynamics of normal human blood inferred from somatic mutations.

Haematopoietic stem cells drive blood production, but their population size and lifetime dynamics have not been quantified directly in humans. Here we identified 129,582 spontaneous, genome-wide somatic mutations in 140 single-cell-derived haematopoietic stem and progenitor colonies from a healthy 59-year-old man and applied population-genetics approaches to reconstruct clonal dynamics. Cell divisions from early embryogenesis were evident in the phylogenetic tree; all blood cells were derived from a common ancestor that preceded gastrulation. The size of the stem cell population grew steadily in early life, reaching a stable plateau by adolescence. We estimate the numbers of haematopoietic stem cells that are actively making white blood cells at any one time to be in the range of 50,000-200,000. We observed adult haematopoietic stem cell clones that generate multilineage outputs, including granulocytes and B lymphocytes. Harnessing naturally occurring mutations to report the clonal architecture of an organ enables the high-resolution reconstruction of somatic cell dynamics in humans.

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.).

Thrombopoietin signaling to chromatin elicits rapid and pervasive epigenome remodeling within poised chromatin architectures.

Thrombopoietin (TPO) is a critical cytokine regulating hematopoietic stem cell maintenance and differentiation into the megakaryocytic lineage. However, the transcriptional and chromatin dynamics elicited by TPO signaling are poorly understood. Here, we study the immediate early transcriptional and cis-regulatory responses to TPO in hematopoietic stem/progenitor cells (HSPCs) and use this paradigm of cytokine signaling to chromatin to dissect the relation between cis- regulatory activity and chromatin architecture. We show that TPO profoundly alters the transcriptome of HSPCs, with key hematopoietic regulators being transcriptionally repressed within 30 minutes of TPO. By examining cis-regulatory dynamics and chromatin architectures, we demonstrate that these changes are accompanied by rapid and extensive epigenome remodeling of cis-regulatory landscapes that is spatially coordinated within topologically associating domains (TADs). Moreover, TPO-responsive enhancers are spatially clustered and engage in preferential homotypic intra- and inter-TAD interactions that are largely refractory to TPO signaling. By further examining the link between cis-regulatory dynamics and chromatin looping, we show that rapid modulation of cis-regulatory activity is largely independent of chromatin looping dynamics. Finally, we show that, although activated and repressed cis-regulatory elements share remarkably similar DNA sequence compositions, transcription factor binding patterns accurately predict rapid cis-regulatory responses to TPO.

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.

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.

Janus kinase deregulation in leukemia and lymphoma.

Genetic alterations affecting members of the Janus kinase (JAK) family have been discovered in a wide array of cancers and are particularly prominent in hematological malignancies. In this review, we focus on the role of such lesions in both myeloid and lymphoid tumors. Oncogenic JAK molecules can activate a myriad of canonical downstream signaling pathways as well as directly interact with chromatin in noncanonical processes, the interplay of which results in a plethora of diverse biological consequences. Deciphering these complexities is shedding unexpected light on fundamental cellular mechanisms and will also be important for improved diagnosis, identification of new therapeutic targets, and the development of stratified approaches to therapy.

Cytokine-induced megakaryocytic differentiation is regulated by genome-wide loss of a uSTAT transcriptional program.

Metazoan development is regulated by transcriptional networks, which must respond to extracellular cues including cytokines. The JAK/STAT pathway is a highly conserved regulatory module, activated by many cytokines, in which tyrosine-phosphorylated STATs (pSTATs) function as transcription factors. However, the mechanisms by which STAT activation modulates lineage-affiliated transcriptional programs are unclear. We demonstrate that in the absence of thrombopoietin (TPO), tyrosine-unphosphorylated STAT5 (uSTAT5) is present in the nucleus where it colocalizes with CTCF and represses a megakaryocytic transcriptional program. TPO-mediated phosphorylation of STAT5 triggers its genome-wide relocation to STAT consensus sites with two distinct transcriptional consequences, loss of a uSTAT5 program that restrains megakaryocytic differentiation and activation of a canonical pSTAT5-driven program which includes regulators of apoptosis and proliferation. Transcriptional repression by uSTAT5 reflects restricted access of the megakaryocytic transcription factor ERG to target genes. These results identify a previously unrecognized mechanism of cytokine-mediated differentiation.

SC3: consensus clustering of single-cell RNA-seq data.

Single-cell RNA-seq enables the quantitative characterization of cell types based on global transcriptome profiles. We present single-cell consensus clustering (SC3), a user-friendly tool for unsupervised clustering, which achieves high accuracy and robustness by combining multiple clustering solutions through a consensus approach (http://bioconductor.org/packages/SC3). We demonstrate that SC3 is capable of identifying subclones from the transcriptomes of neoplastic cells collected from patients.

Mutant calreticulin knockin mice develop thrombocytosis and myelofibrosis without a stem cell self-renewal advantage.

Somatic mutations in the endoplasmic reticulum chaperone calreticulin (CALR) are detected in approximately 40% of patients with essential thrombocythemia (ET) and primary myelofibrosis (PMF). Multiple different mutations have been reported, but all result in a +1-bp frameshift and generate a novel protein C terminus. In this study, we generated a conditional mouse knockin model of the most common CALR mutation, a 52-bp deletion. The mutant novel human C-terminal sequence is integrated into the otherwise intact mouse CALR gene and results in mutant CALR expression under the control of the endogenous mouse locus. CALRdel/+ mice develop a transplantable ET-like disease with marked thrombocytosis, which is associated with increased and morphologically abnormal megakaryocytes and increased numbers of phenotypically defined hematopoietic stem cells (HSCs). Homozygous CALRdel/del mice developed extreme thrombocytosis accompanied by features of MF, including leukocytosis, reduced hematocrit, splenomegaly, and increased bone marrow reticulin. CALRdel/+ HSCs were more proliferative in vitro, but neither CALRdel/+ nor CALRdel/del displayed a competitive transplantation advantage in primary or secondary recipient mice. These results demonstrate the consequences of heterozygous and homozygous CALR mutations and provide a powerful model for dissecting the pathogenesis of CALR-mutant ET and PMF.

Single-cell approaches identify the molecular network driving malignant hematopoietic stem cell self-renewal.

Recent advances in single-cell technologies have permitted the investigation of heterogeneous cell populations at previously unattainable resolution. Here we apply such approaches to resolve the molecular mechanisms driving disease in mouse hematopoietic stem cells (HSCs), using JAK2V617F mutant myeloproliferative neoplasms (MPNs) as a model. Single-cell gene expression and functional assays identified a subset of JAK2V617F mutant HSCs that display defective self-renewal. This defect is rescued at the single HSC level by crossing JAK2V617F mice with mice lacking TET2, the most commonly comutated gene in patients with MPN. Single-cell gene expression profiling of JAK2V617F-mutant HSCs revealed a loss of specific regulator genes, some of which were restored to normal levels in single TET2/JAK2 mutant HSCs. Of these, Bmi1 and, to a lesser extent, Pbx1 and Meis1 overexpression in JAK2-mutant HSCs could drive a disease phenotype and retain durable stem cell self-renewal in functional assays. Together, these single-cell approaches refine the molecules involved in clonal expansion of MPNs and have broad implications for deconstructing the molecular network of normal and malignant stem cells.

The unfolded protein response governs integrity of the haematopoietic stem-cell pool during stress.

The blood system is sustained by a pool of haematopoietic stem cells (HSCs) that are long-lived due to their capacity for self-renewal. A consequence of longevity is exposure to stress stimuli including reactive oxygen species (ROS), nutrient fluctuation and DNA damage. Damage that occurs within stressed HSCs must be tightly controlled to prevent either loss of function or the clonal persistence of oncogenic mutations that increase the risk of leukaemogenesis. Despite the importance of maintaining cell integrity throughout life, how the HSC pool achieves this and how individual HSCs respond to stress remain poorly understood. Many sources of stress cause misfolded protein accumulation in the endoplasmic reticulum (ER), and subsequent activation of the unfolded protein response (UPR) enables the cell to either resolve stress or initiate apoptosis. Here we show that human HSCs are predisposed to apoptosis through strong activation of the PERK branch of the UPR after ER stress, whereas closely related progenitors exhibit an adaptive response leading to their survival. Enhanced ER protein folding by overexpression of the co-chaperone ERDJ4 (also called DNAJB9) increases HSC repopulation capacity in xenograft assays, linking the UPR to HSC function. Because the UPR is a focal point where different sources of stress converge, our study provides a framework for understanding how stress signalling is coordinated within tissue hierarchies and integrated with stemness. Broadly, these findings reveal that the HSC pool maintains clonal integrity by clearance of individual HSCs after stress to prevent propagation of damaged stem cells.

Myeloproliferative neoplasms: from origins to outcomes.

Substantial progress has been made in our understanding of the pathogenetic basis of myeloproliferative neoplasms. The discovery of mutations in JAK2 over a decade ago heralded a new age for patient care as a consequence of improved diagnosis and the development of therapeutic JAK inhibitors. The more recent identification of mutations in calreticulin brought with it a sense of completeness, with most patients with myeloproliferative neoplasm now having a biological basis for their excessive myeloproliferation. We are also beginning to understand the processes that lead to acquisition of somatic mutations and the factors that influence subsequent clonal expansion and emergence of disease. Extended genomic profiling has established a multitude of additional acquired mutations, particularly prevalent in myelofibrosis, where their presence carries prognostic implications. A major goal is to integrate genetic, clinical, and laboratory features to identify patients who share disease biology and clinical outcome, such that therapies, both existing and novel, can be better targeted.

Effect of mutation order on myeloproliferative neoplasms.

BACKGROUND: Cancers result from the accumulation of somatic mutations, and their properties are thought to reflect the sum of these mutations. However, little is known about the effect of the order in which mutations are acquired. METHODS: We determined mutation order in patients with myeloproliferative neoplasms by genotyping hematopoietic colonies or by means of next-generation sequencing. Stem cells and progenitor cells were isolated to study the effect of mutation order on mature and immature hematopoietic cells. RESULTS: The age at which a patient presented with a myeloproliferative neoplasm, acquisition of JAK2 V617F homozygosity, and the balance of immature progenitors were all influenced by mutation order. As compared with patients in whom the TET2 mutation was acquired first (hereafter referred to as "TET2-first patients"), patients in whom the Janus kinase 2 (JAK2) mutation was acquired first ("JAK2-first patients") had a greater likelihood of presenting with polycythemia vera than with essential thrombocythemia, an increased risk of thrombosis, and an increased sensitivity of JAK2-mutant progenitors to ruxolitinib in vitro. Mutation order influenced the proliferative response to JAK2 V617F and the capacity of double-mutant hematopoietic cells and progenitor cells to generate colony-forming cells. Moreover, the hematopoietic stem-and-progenitor-cell compartment was dominated by TET2 single-mutant cells in TET2-first patients but by JAK2-TET2 double-mutant cells in JAK2-first patients. Prior mutation of TET2 altered the transcriptional consequences of JAK2 V617F in a cell-intrinsic manner and prevented JAK2 V617F from up-regulating genes associated with proliferation. CONCLUSIONS: The order in which JAK2 and TET2 mutations were acquired influenced clinical features, the response to targeted therapy, the biology of stem and progenitor cells, and clonal evolution in patients with myeloproliferative neoplasms. (Funded by Leukemia and Lymphoma Research and others.).