Computational tools for assessing gene therapy under branching process models of mutation

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.