Modelling acute myeloid leukaemia in a continuum of differentiation states

Detection of continuous hierarchical heterogeneity by single-cell surface antigen analysis in the prognosis evaluation of acute myeloid leukaemia

BACKGROUND

Acute myeloid leukaemia (AML) is characterised by the malignant accumulation of myeloid progenitors with a high recurrence rate after chemotherapy. Blasts (leukaemia cells) exhibit a complete myeloid differentiation hierarchy hiding a wide range of temporal information from initial to mature clones, including genesis, phenotypic transformation, and cell fate decisions, which might contribute to relapse in AML patients.

METHODS

Based on the landscape of AML surface antigens generated by mass cytometry (CyTOF), we combined manifold analysis and principal curve-based trajectory inference algorithm to align myelocytes on a single-linear evolution axis by considering their phenotype continuum that correlated with differentiation order. Backtracking the trajectory from mature clusters located automatically at the terminal, we recurred the molecular dynamics during AML progression and confirmed the evolution stage of single cells. We also designed a 'dispersive antigens in neighbouring clusters exhibition (DANCE)' feature selection method to simplify and unify trajectories, which enabled the exploration and comparison of relapse-related traits among 43 paediatric AML bone marrow specimens.

RESULTS

The feasibility of the proposed trajectory analysis method was verified with public datasets. After aligning single cells on the pseudotime axis, primitive clones were recognized precisely from AML blasts, and the expression of the inner molecules before and after drug stimulation was accurately plotted on the trajectory. Applying DANCE to 43 clinical samples with different responses for chemotherapy, we selected 12 antigens as a general panel for myeloblast differentiation performance, and obtain trajectories to those patients. For the trajectories with unified molecular dynamics, CD11c overexpression in the primitive stage indicated a good chemotherapy outcome. Moreover, a later initial peak of stemness heterogeneity tended to be associated with a higher risk of relapse compared with complete remission.

CONCLUSIONS

In this study, pseudotime was generated as a new single-cell feature. Minute differences in temporal traits among samples could be exhibited on a trajectory, thus providing a new strategy for predicting AML relapse and monitoring drug responses over time scale.

Comparison of cell state models derived from single-cell RNA sequencing data: graph versus multi-dimensional space

Single-cell sequencing technologies have revolutionized molecular and cellular biology and stimulated the development of computational tools to analyze the data generated from these technology platforms. However, despite the recent explosion of computational analysis tools, relatively few mathematical models have been developed to utilize these data. Here we compare and contrast two cell state geometries for building mathematical models of cell state-transitions with single-cell RNA-sequencing data with hematopoeisis as a model system; (i) by using partial differential equations on a graph representing intermediate cell states between known cell types, and (ii) by using the equations on a multi-dimensional continuous cell state-space. As an application of our approach, we demonstrate how the calibrated models may be used to mathematically perturb normal hematopoeisis to simulate, predict, and study the emergence of novel cell states during the pathogenesis of acute myeloid leukemia. We particularly focus on comparing the strength and weakness of the graph model and multi-dimensional model.

Differentiation of leukemic blasts is not completely blocked in acute myeloid leukemia

Hematopoiesis, the formation of blood cells, involves the hierarchical differentiation of immature blast cells into mature, functional cell types and lineages of the immune system. Hematopoietic stem cells precisely regulate self-renewal versus differentiation to balance the production of blood cells and maintenance of the stem cell pool. The canonical view of acute myeloid leukemia (AML) is that it results from a combination of molecular events in a hematopoietic stem cell that block differentiation and drive proliferation. These events result in the accumulation of primitive hematopoietic blast cells in the blood and bone marrow. We used mathematical modeling to determine the impact of varying differentiation rates on myeloblastic accumulation. Our model shows that, instead of the commonly held belief that AML results from a complete block of differentiation of the hematopoietic stem cell, even a slight skewing of the fraction of cells that differentiate would produce an accumulation of blasts. We confirmed this model by interphase fluorescent in situ hybridization (FISH) and sequencing of purified cell populations from patients with AML, which showed that different leukemia-causing molecular abnormalities typically thought to block differentiation were consistently present in mature myeloid cells such as neutrophils and monocytes at similar levels to those in immature myeloid cells. These findings suggest reduced or skewed, rather than blocked, differentiation is responsible for the development of AML. Approaches that restore normal regulation of hematopoiesis could be effective treatment strategies.

Opportunities for improving cancer treatment using systems biology

Current cancer therapies target a limited set of tumor features, rather than considering the tumor as a whole. Systems biology aims to reveal therapeutic targets associated with a variety of facets in an individual's tumor, such as genetic heterogeneity and its evolution, cancer cell-autonomous phenotypes, and microenvironmental signaling. These disparate characteristics can be reconciled using mathematical modeling that incorporates concepts from ecology and evolution. This provides an opportunity to predict tumor growth and response to therapy, to tailor patient-specific approaches in real time or even prospectively. Importantly, as data regarding patient tumors is often available from only limited time points during treatment, systems-based approaches can address this limitation by interpolating longitudinal events within a principled framework. This review outlines areas in medicine that could benefit from systems biology approaches to deconvolve the complexity of cancer.

Inferring population dynamics from single-cell RNA-sequencing time series data

Recent single-cell RNA-sequencing studies have suggested that cells follow continuous transcriptomic trajectories in an asynchronous fashion during development. However, observations of cell flux along trajectories are confounded with population size effects in snapshot experiments and are therefore hard to interpret. In particular, changes in proliferation and death rates can be mistaken for cell flux. Here we present pseudodynamics, a mathematical framework that reconciles population dynamics with the concepts underlying developmental trajectories inferred from time-series single-cell data. Pseudodynamics models population distribution shifts across trajectories to quantify selection pressure, population expansion, and developmental potentials. Applying this model to time-resolved single-cell RNA-sequencing of T-cell and pancreatic beta cell maturation, we characterize proliferation and apoptosis rates and identify key developmental checkpoints, data inaccessible to existing approaches.