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.