1
|
Hermange G, Vainchenker W, Plo I, Cournède PH. Mathematical modelling, selection and hierarchical inference to determine the minimal dose in IFNα therapy against myeloproliferative neoplasms. MATHEMATICAL MEDICINE AND BIOLOGY : A JOURNAL OF THE IMA 2024; 41:110-134. [PMID: 38875109 DOI: 10.1093/imammb/dqae006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 01/30/2024] [Accepted: 05/02/2024] [Indexed: 06/16/2024]
Abstract
Myeloproliferative neoplasms (MPN) are blood cancers that appear after acquiring a driver mutation in a hematopoietic stem cell. These hematological malignancies result in the overproduction of mature blood cells and, if not treated, induce a risk of cardiovascular events and thrombosis. Pegylated IFN$\alpha $ is commonly used to treat MPN, but no clear guidelines exist concerning the dose prescribed to patients. We applied a model selection procedure and ran a hierarchical Bayesian inference method to decipher how dose variations impact the response to the therapy. We inferred that IFN$\alpha $ acts on mutated stem cells by inducing their differentiation into progenitor cells; the higher the dose, the higher the effect. We found that the treatment can induce long-term remission when a sufficient (patient-dependent) dose is reached. We determined this minimal dose for individuals in a cohort of patients and estimated the most suitable starting dose to give to a new patient to increase the chances of being cured.
Collapse
Affiliation(s)
- Gurvan Hermange
- Université Paris-Saclay, CentraleSupélec, Laboratory of Mathematics and Informatics (MICS), Gif-sur-Yvette, France
| | - William Vainchenker
- INSERM U1287 (INSERM, Gustave Roussy, Université Paris-Saclay), Villejuif, France
- Gustave Roussy, Villejuif, France
- Université Paris-Saclay, Villejuif, France
| | - Isabelle Plo
- INSERM U1287 (INSERM, Gustave Roussy, Université Paris-Saclay), Villejuif, France
- Gustave Roussy, Villejuif, France
- Université Paris-Saclay, Villejuif, France
| | - Paul-Henry Cournède
- Université Paris-Saclay, CentraleSupélec, Laboratory of Mathematics and Informatics (MICS), Gif-sur-Yvette, France
| |
Collapse
|
2
|
Del Core L, Pellin D, Wit EC, Grzegorczyk MA. Scalable inference of cell differentiation networks in gene therapy clonal tracking studies of haematopoiesis. Bioinformatics 2023; 39:btad605. [PMID: 37774002 PMCID: PMC10585354 DOI: 10.1093/bioinformatics/btad605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2022] [Revised: 07/15/2023] [Accepted: 09/28/2023] [Indexed: 10/01/2023] Open
Abstract
MOTIVATION Investigating cell differentiation under a genetic disorder offers the potential for improving current gene therapy strategies. Clonal tracking provides a basis for mathematical modelling of population stem cell dynamics that sustain the blood cell formation, a process known as haematopoiesis. However, many clonal tracking protocols rely on a subset of cell types for the characterization of the stem cell output, and the data generated are subject to measurement errors and noise. RESULTS We propose a stochastic framework to infer dynamic models of cell differentiation from clonal tracking data. A state-space formulation combines a stochastic quasi-reaction network, describing cell differentiation, with a Gaussian measurement model accounting for data errors and noise. We developed an inference algorithm based on an extended Kalman filter, a nonlinear optimization, and a Rauch-Tung-Striebel smoother. Simulations show that our proposed method outperforms the state-of-the-art and scales to complex structures of cell differentiations in terms of nodes size and network depth. The application of our method to five in vivo gene therapy studies reveals different dynamics of cell differentiation. Our tool can provide statistical support to biologists and clinicians to better understand cell differentiation and haematopoietic reconstitution after a gene therapy treatment. The equations of the state-space model can be modified to infer other dynamics besides cell differentiation. AVAILABILITY AND IMPLEMENTATION The stochastic framework is implemented in the R package Karen which is available for download at https://cran.r-project.org/package=Karen. The code that supports the findings of this study is openly available at https://github.com/delcore-luca/CellDifferentiationNetworks.
Collapse
Affiliation(s)
- Luca Del Core
- University of Groningen – Bernoulli Institute, 9747AG Groningen, The Netherlands
- University of Nottingham – School of Mathematical Sciences, Nottingham NG72RD, United Kingdom
| | - Danilo Pellin
- Harvard Medical School, Boston, MA 02115, United States
| | - Ernst C Wit
- University of Groningen – Bernoulli Institute, 9747AG Groningen, The Netherlands
- Università della Svizzera italiana – Institute of Computing, 6962 Lugano, Switzerland
| | - Marco A Grzegorczyk
- University of Groningen – Bernoulli Institute, 9747AG Groningen, The Netherlands
| |
Collapse
|
3
|
Stutz TC, Sinsheimer JS, Sehl M, Xu J. Computational tools for assessing gene therapy under branching process models of mutation. Bull Math Biol 2021; 84:15. [PMID: 34870755 DOI: 10.1007/s11538-021-00969-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2021] [Accepted: 11/15/2021] [Indexed: 11/28/2022]
Abstract
Multitype branching processes are ideal for studying the population dynamics of stem cell populations undergoing mutation accumulation over the years following transplant. In such stochastic models, several quantities are of clinical interest as insertional mutagenesis carries the potential threat of leukemogenesis following gene therapy with autologous stem cell transplantation. In this paper, we develop a three-type branching process model describing accumulations of mutations in a population of stem cells distinguished by their ability for long-term self-renewal. Our outcome of interest is the appearance of a double-mutant cell, which carries a high potential for leukemic transformation. In our model, a single-hit mutation carries a slight proliferative advantage over a wild-type stem cells. We compute marginalized transition probabilities that allow us to capture important quantitative aspects of our model, including the probability of observing a double-hit mutant and relevant moments of a single-hit mutation population over time. We thoroughly explore the model behavior numerically, varying birth rates across the initial sizes and populations of wild type stem cells and single-hit mutants, and compare the probability of observing a double-hit mutant under these conditions. We find that increasing the number of single-mutants over wild-type particles initially present has a large effect on the occurrence of a double-mutant, and that it is relatively safe for single-mutants to be quite proliferative, provided the lentiviral gene addition avoids creating single mutants in the original insertion process. Our approach is broadly applicable to an important set of questions in cancer modeling and other population processes involving multiple stages, compartments, or types.
Collapse
Affiliation(s)
- Timothy C Stutz
- Department of Computational Medicine, University of California, Los Angeles, CA, USA
| | - Janet S Sinsheimer
- Departments of Biostatistics, Computational Medicine, Human Genetics, University of California, Los Angeles, CA, USA
| | - Mary Sehl
- Department of Computational Medicine and Division of Hematology-Oncology, Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, CA, USA
| | - Jason Xu
- Department of Statistical Science, Duke University, Durham, NC, USA.
| |
Collapse
|
4
|
Abstract
PURPOSE OF REVIEW The hematopoietic compartment is tasked with the establishment and maintenance of the entire blood program in steady-state and in response to stress. Key to this process are hematopoietic stem cells (HSCs), which possess the unique ability to self-renew and differentiate to replenish blood cells throughout an organism's lifetime. Though tightly regulated, the hematopoietic system is vulnerable to both intrinsic and extrinsic factors that influence hematopoietic stem and progenitor cell (HSPC) fate. Here, we review recent advances in our understanding of hematopoietic regulation under stress conditions such as inflammation, aging, mitochondrial defects, and damage to DNA or endoplasmic reticulum. RECENT FINDINGS Recent studies have illustrated the vast mechanisms involved in regulating stress-induced hematopoiesis, including cytokine-mediated lineage bias, gene signature changes in aged HSCs associated with chronic inflammation, the impact of clonal hematopoiesis and stress tolerance, characterization of the HSPC response to endoplasmic reticulum stress and of several epigenetic regulators that influence HSPC response to cell cycle stress. SUMMARY Several key recent findings have deepened our understanding of stress hematopoiesis. These studies will advance our abilities to reduce the impact of stress in disease and aging through clinical interventions to treat stress-related outcomes.
Collapse
|
5
|
Lomeli LM, Iniguez A, Tata P, Jena N, Liu ZY, Van Etten R, Lander AD, Shahbaba B, Lowengrub JS, Minin VN. Optimal experimental design for mathematical models of haematopoiesis. J R Soc Interface 2021; 18:20200729. [PMID: 33499768 PMCID: PMC7879761 DOI: 10.1098/rsif.2020.0729] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Accepted: 01/04/2021] [Indexed: 11/12/2022] Open
Abstract
The haematopoietic system has a highly regulated and complex structure in which cells are organized to successfully create and maintain new blood cells. It is known that feedback regulation is crucial to tightly control this system, but the specific mechanisms by which control is exerted are not completely understood. In this work, we aim to uncover the underlying mechanisms in haematopoiesis by conducting perturbation experiments, where animal subjects are exposed to an external agent in order to observe the system response and evolution. We have developed a novel Bayesian hierarchical framework for optimal design of perturbation experiments and proper analysis of the data collected. We use a deterministic model that accounts for feedback and feedforward regulation on cell division rates and self-renewal probabilities. A significant obstacle is that the experimental data are not longitudinal, rather each data point corresponds to a different animal. We overcome this difficulty by modelling the unobserved cellular levels as latent variables. We then use principles of Bayesian experimental design to optimally distribute time points at which the haematopoietic cells are quantified. We evaluate our approach using synthetic and real experimental data and show that an optimal design can lead to better estimates of model parameters.
Collapse
Affiliation(s)
- Luis Martinez Lomeli
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
| | - Abdon Iniguez
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
| | - Prasanthi Tata
- Division of Hematology/Oncology, University of California Irvine, Irvine, CA, USA
| | - Nilamani Jena
- Division of Hematology/Oncology, University of California Irvine, Irvine, CA, USA
| | - Zhong-Ying Liu
- Division of Hematology/Oncology, University of California Irvine, Irvine, CA, USA
| | - Richard Van Etten
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
- Division of Hematology/Oncology, University of California Irvine, Irvine, CA, USA
- Department of Biological Chemistry, University of California Irvine, Irvine, CA, USA
- Center for Cancer Systems Biology, University of California Irvine, Irvine, CA, USA
- Chao Family Comprehensive Cancer Center, University of California Irvine, Irvine, CA, USA
| | - Arthur D. Lander
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
- Center for Cancer Systems Biology, University of California Irvine, Irvine, CA, USA
- Chao Family Comprehensive Cancer Center, University of California Irvine, Irvine, CA, USA
- Department of Developmental and Cell Biology, University of California Irvine, Irvine, CA, USA
- Department of Biomedical Engineering, University of California Irvine, Irvine, CA, USA
| | - Babak Shahbaba
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
- Center for Cancer Systems Biology, University of California Irvine, Irvine, CA, USA
- Department of Statistics, University of California Irvine, Irvine, CA, USA
| | - John S. Lowengrub
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
- Center for Cancer Systems Biology, University of California Irvine, Irvine, CA, USA
- Chao Family Comprehensive Cancer Center, University of California Irvine, Irvine, CA, USA
- Department of Biomedical Engineering, University of California Irvine, Irvine, CA, USA
- Department of Mathematics, University of California Irvine, Irvine, CA, USA
| | - Vladimir N. Minin
- Center for Complex Biological Systems, University of California Irvine, Irvine, CA, USA
- Center for Cancer Systems Biology, University of California Irvine, Irvine, CA, USA
- Department of Statistics, University of California Irvine, Irvine, CA, USA
| |
Collapse
|