Measures of genetic diversification in somatic tissues at bulk and single-cell resolution.

Intra-tissue genetic heterogeneity is universal to both healthy and cancerous tissues. It emerges from the stochastic accumulation of somatic mutations throughout development and homeostasis. By combining population genetics theory and genomic information, genetic heterogeneity can be exploited to infer tissue organization and dynamics in vivo. However, many basic quantities, for example the dynamics of tissue-specific stem cells remain difficult to quantify precisely. Here, we show that single-cell and bulk sequencing data inform on different aspects of the underlying stochastic processes. Bulk-derived variant allele frequency spectra (VAF) show transitions from growing to constant stem cell populations with age in samples of healthy esophagus epithelium. Single-cell mutational burden distributions allow a sample size independent measure of mutation and proliferation rates. Mutation rates in adult hematopietic stem cells are higher compared to inferences during development, suggesting additional proliferation-independent effects. Furthermore, single-cell derived VAF spectra contain information on the number of tissue-specific stem cells. In hematopiesis, we find approximately 2 × 105 HSCs, if all stem cells divide symmetrically. However, the single-cell mutational burden distribution is over-dispersed compared to a model of Poisson distributed random mutations. A time-associated model of mutation accumulation with a constant rate alone cannot generate such a pattern. At least one additional source of stochasticity would be needed. Possible candidates for these processes may be occasional bursts of stem cell divisions, potentially in response to injury, or non-constant mutation rates either through environmental exposures or cell-intrinsic variation.

Brownian motion in a growing population of ballistic particles.

We investigate the motility of a growing population of cells in a idealized setting: We consider a system of hard disks in which new particles are added according to prescribed growth kinetics, thereby dynamically changing the number density. As a result, the expected Brownian motion of the hard disks is modified. We compute the density-dependent friction of the hard disks and insert it in an effective Langevin equation to describe the system, assuming that the intercollision time is smaller than the timescale of the growth. We find that the effective Langevin description captures the changes in motility, in agreement with the simulation results. Our framework can be extended to other systems in which the transport coefficient varies with time.

Multistage feedback-driven compartmental dynamics of hematopoiesis.

Human hematopoiesis is surprisingly resilient to disruptions, providing suitable responses to severe bleeding, long-lasting immune activation, and even bone marrow transplants. Still, many blood disorders exist which push the system past its natural plasticity, resulting in abnormalities in the circulating blood. While proper treatment of such diseases can benefit from understanding the underlying cell dynamics, these are non-trivial to predict due to the hematopoietic system's hierarchical nature and complex feedback networks. To characterize the dynamics following different types of perturbations, we investigate a model representing hematopoiesis as a sequence of compartments covering all maturation stages-from stem to mature cells-where feedback regulates cell production to ongoing necessities. We find that a stable response to perturbations requires the simultaneous adaptation of cell differentiation and self-renewal rates, and show that under conditions of continuous disruption-as found in chronic hemolytic states-compartment cell numbers evolve to novel stable states.

Evolutionary dynamics of paroxysmal nocturnal hemoglobinuria.

Paroxysmal nocturnal hemoglobinuria (PNH) is an acquired clonal blood disorder characterized by hemolysis and a high risk of thrombosis, that is due to a deficiency in several cell surface proteins that prevent complement activation. Its origin has been traced to a somatic mutation in the PIG-A gene within hematopoietic stem cells (HSC). However, to date the question of how this mutant clone expands in size to contribute significantly to hematopoiesis remains under debate. One hypothesis posits the existence of a selective advantage of PIG-A mutated cells due to an immune mediated attack on normal HSC, but the evidence supporting this hypothesis is inconclusive. An alternative (and simpler) explanation attributes clonal expansion to neutral drift, in which case selection neither favours nor inhibits expansion of PIG-A mutated HSC. Here we examine the implications of the neutral drift model by numerically evolving a Markov chain for the probabilities of all possible outcomes, and investigate the possible occurrence and evolution, within this framework, of multiple independently arising clones within the HSC pool. Predictions of the model agree well with the known incidence of the disease and average age at diagnosis. Notwithstanding the slight difference in clonal expansion rates between our results and those reported in the literature, our model results lead to a relative stability of clone size when averaging multiple cases, in accord with what has been observed in human trials. The probability of a patient harbouring a second clone in the HSC pool was found to be extremely low ([Formula: see text]). Thus our results suggest that in clinical cases of PNH where two independent clones of mutant cells are observed, only one of those is likely to have originated in the HSC pool.