Exact site frequency spectra of neutrally evolving tumors: A transition between power laws reveals a signature of cell viability
Abstract
The site frequency spectrum (SFS) is a popular summary statistic of genomic data. While the SFS of a constant-sized population undergoing neutral mutations has been extensively studied in population genetics, the rapidly growing amount of cancer genomic data has attracted interest in the spectrum of an exponentially growing population. Recent theoretical results have generally dealt with special or limiting cases, such as considering only cells with an infinite line of descent, assuming deterministic tumor growth, or taking large-time or large-population limits. In this work, we derive exact expressions for the expected SFS of a cell population that evolves according to a stochastic branching process, first for cells with an infinite line of descent and then for the total population, evaluated either at a fixed time (fixed-time spectrum) or at the stochastic time at which the population reaches a certain size (fixed-size spectrum). We find that while the rate of mutation scales the SFS of the total population linearly, the rates of cell birth and cell death change the shape of the spectrum at the small-frequency end, inducing a transition between a
- Publication:
-
Theoretical Population Biology
- Pub Date:
- December 2021
- DOI:
- 10.1016/j.tpb.2021.09.004
- arXiv:
- arXiv:2102.11959
- Bibcode:
- 2021TPBio.142...67G
- Keywords:
-
- Mathematical modeling;
- Branching processes;
- Exponentially growing populations;
- Site frequency spectrum;
- Infinite sites model;
- Cancer evolution;
- Quantitative Biology - Populations and Evolution;
- Mathematics - Probability;
- 92D25;
- 92B05;
- 60J85
- E-Print:
- 47 pages, 15 figures, accepted for publication in Theoretical Population Biology