A structured population model of clonal selection in acute leukemias with multiple maturation stages

Mathematics of neural stem cells: Linking data and processes

Adult stem cells are described as a discrete population of cells that stand at the top of a hierarchy of progressively differentiating cells. Through their unique ability to self-renew and differentiate, they regulate the number of end-differentiated cells that contribute to tissue physiology. The question of how discrete, continuous, or reversible the transitions through these hierarchies are and the precise parameters that determine the ultimate performance of stem cells in adulthood are the subject of intense research. In this review, we explain how mathematical modelling has improved the mechanistic understanding of stem cell dynamics in the adult brain. We also discuss how single-cell sequencing has influenced the understanding of cell states or cell types. Finally, we discuss how the combination of single-cell sequencing technologies and mathematical modelling provides a unique opportunity to answer some burning questions in the field of stem cell biology.

Local asymptotic stability of a system of integro-differential equations describing clonal evolution of a self-renewing cell population under mutation

In this paper we consider a system of non-linear integro-differential equations (IDEs) describing evolution of a clonally heterogeneous population of malignant white blood cells (leukemic cells) undergoing mutation and clonal selection. We prove existence and uniqueness of non-trivial steady states and study their asymptotic stability. The results are compared to those of the system without mutation. Existence of equilibria is proved by formulating the steady state problem as an eigenvalue problem and applying a version of the Krein-Rutmann theorem for Banach lattices. The stability at equilibrium is analysed using linearisation and the Weinstein-Aronszajn determinant which allows to conclude local asymptotic stability.

Computational Reconstruction of Clonal Hierarchies From Bulk Sequencing Data of Acute Myeloid Leukemia Samples

Acute myeloid leukemia is an aggressive cancer of the blood forming system. The malignant cell population is composed of multiple clones that evolve over time. Clonal data reflect the mechanisms governing treatment response and relapse. Single cell sequencing provides most direct insights into the clonal composition of the leukemic cells, however it is still not routinely available in clinical practice. In this work we develop a computational algorithm that allows identifying all clonal hierarchies that are compatible with bulk variant allele frequencies measured in a patient sample. The clonal hierarchies represent descendance relations between the different clones and reveal the order in which mutations have been acquired. The proposed computational approach is tested using single cell sequencing data that allow comparing the outcome of the algorithm with the true structure of the clonal hierarchy. We investigate which problems occur during reconstruction of clonal hierarchies from bulk sequencing data. Our results suggest that in many cases only a small number of possible hierarchies fits the bulk data. This implies that bulk sequencing data can be used to obtain insights in clonal evolution.

Dynamical properties of feedback signalling in B lymphopoiesis: A mathematical modelling approach

Haematopoiesis is the process of generation of blood cells. Lymphopoiesis generates lymphocytes, the cells in charge of the adaptive immune response. Disruptions of this process are associated with diseases like leukaemia, which is especially incident in children. The characteristics of self-regulation of this process make them suitable for a mathematical study. In this paper we develop mathematical models of lymphopoiesis using currently available data. We do this by drawing inspiration from existing structured models of cell lineage development and integrating them with paediatric bone marrow data, with special focus on regulatory mechanisms. A formal analysis of the models is carried out, giving steady states and their stability conditions. We use this analysis to obtain biologically relevant regions of the parameter space and to understand the dynamical behaviour of B-cell renovation. Finally, we use numerical simulations to obtain further insight into the influence of proliferation and maturation rates on the reconstitution of the cells in the B line. We conclude that a model including feedback regulation of cell proliferation represents a biologically plausible depiction for B-cell reconstitution in bone marrow. Research into haematological disorders could benefit from a precise dynamical description of B lymphopoiesis.

CELLULAR FEEDBACK NETWORKS AND THEIR RESILIENCE AGAINST MUTATIONS

Many tissues undergo a steady turnover, where cell divisions are on average balanced with cell deaths. Cell fate decisions such as stem cell (SC) differentiations, proliferations, or differentiated cell (DC) deaths, may be controlled by cell populations through cell-to-cell signaling. Here, we examine a class of mathematical models of turnover in SC lineages to understand engineering design principles of control (feedback) loops, that may operate in such systems. By using ordinary differential equations that describe the co-dynamics of SCs and DCs, we study the effect of different types of mutations that interfere with feedback present within cellular networks. For instance, we find that mutants that do not participate in feedback are less dangerous in the sense that they will not rise from low numbers, whereas mutants that do not respond to feedback signals could rise and replace the wild-type population. Additionally, we asked if different feedback networks can have different degrees of resilience against such mutations. We found that all minimal networks, that is networks consisting of exactly one feedback loop that is sufficient for homeostatic stability of the wild-type population, are equally vulnerable. Mutants with a weakened/eliminated feedback parameter might expand from lower numbers and either enter unlimited growth or reach an equilibrium with an increased number of SCs and DCs. Therefore, from an evolutionary viewpoint, it appears advantageous to combine feedback loops, creating redundant feedback networks. Interestingly, from an engineering prospective, not all such redundant systems are equally resilient. For some of them, any mutation that weakens/eliminates one of the loops will lead to a population growth of SCs. For others, the population of SCs can actually shrink as a result of “cutting” one of the loops, thus slowing down further unwanted transformations.

Using mathematical models to improve risk-scoring in acute myeloid leukemia

Acute myeloid leukemia (AML) is an aggressive cancer of the blood forming (hematopoietic) system. Due to the high patient variability of disease dynamics, risk-scoring is an important part of its clinical management. AML is characterized by impaired blood cell formation and the accumulation of so-called leukemic blasts in the bone marrow of patients. Recently, it has been proposed to use counts of blood-producing (hematopoietic) stem cells (HSCs) as a biomarker for patient prognosis. In this work, we use a non-linear mathematical model to provide mechanistic evidence for the suitability of HSC counts as a prognostic marker. Using model analysis and computer simulations, we compare different risk-scores involving HSC quantification. We propose and validate a simple approach to improve risk prediction based on HSC and blast counts measured at the time of diagnosis.