Universal Asymptotic Clone Size Distribution for General Population Growth

Deterministically growing (wild-type) populations which seed stochastically developing mutant clones have found an expanding number of applications from microbial populations to cancer. The special case of exponential wild-type population growth, usually termed the Luria–Delbrück or Lea–Coulson model, is often assumed but seldom realistic. In this article, we generalise this model to different types of wild-type population growth, with mutants evolving as a birth–death branching process. Our focus is on the size distribution of clones—that is the number of progeny of a founder mutant—which can be mapped to the total number of mutants. Exact expressions are derived for exponential, power-law and logistic population growth. Additionally, for a large class of population growth, we prove that the long-time limit of the clone size distribution has a general two-parameter form, whose tail decays as a power-law. Considering metastases in cancer as the mutant clones, upon analysing a data-set of their size distribution, we indeed find that a power-law tail is more likely than an exponential one.

[1]  E. Luebeck,et al.  Evaluation of screening strategies for pre-malignant lesions using a biomathematical approach. , 2008, Mathematical biosciences.

[2]  Lurias,et al.  MUTATIONS OF BACTERIA FROM VIRUS SENSITIVITY TO VIRUS RESISTANCE’-’ , 2003 .

[3]  H. M. Taylor,et al.  An introduction to stochastic modeling , 1985 .

[4]  W P Angerer,et al.  An explicit representation of the Luria–Delbrück distribution , 2001, Journal of mathematical biology.

[5]  Krishnendu Chatterjee,et al.  Evolutionary dynamics of cancer in response to targeted combination therapy , 2013, eLife.

[6]  P. A. P. Moran,et al.  An introduction to probability theory , 1968 .

[7]  L. Breuer Introduction to Stochastic Processes , 2022, Statistical Methods for Climate Scientists.

[8]  Daniel E. Geer,et al.  Power. Law , 2012, IEEE Secur. Priv..

[9]  D. Hudson Interval Estimation from the Likelihood Function , 1971 .

[10]  R. Durrett Probability: Theory and Examples , 1993 .

[11]  M. Nowak,et al.  Distant Metastasis Occurs Late during the Genetic Evolution of Pancreatic Cancer , 2010, Nature.

[12]  C. Tomasetti On the Probability of Random Genetic Mutations for Various Types of Tumor Growth , 2012, Bulletin of mathematical biology.

[13]  Daniel W. Lozier,et al.  NIST Digital Library of Mathematical Functions , 2003, Annals of Mathematics and Artificial Intelligence.

[14]  H. Simon,et al.  ON A CLASS OF SKEW DISTRIBUTION FUNCTIONS , 1955 .

[15]  S. Lang Complex Analysis , 1977 .

[16]  Ronald F. Boisvert,et al.  NIST Handbook of Mathematical Functions , 2010 .

[17]  S. Karlin,et al.  A second course in stochastic processes , 1981 .

[18]  E G Luebeck,et al.  A generalized Luria-Delbrück model. , 2005, Mathematical biosciences.

[19]  Herbert Levine,et al.  Scaling Solution in the Large Population Limit of the General Asymmetric Stochastic Luria–Delbrück Evolution Process , 2014, Journal of statistical physics.

[20]  T. Antal,et al.  Mutant number distribution in an exponentially growing population , 2014, 1410.3307.

[21]  Joel L. Schiff,et al.  The Laplace Transform , 1999 .

[22]  Anup Dewanji,et al.  Number and Size Distribution of Colorectal Adenomas under the Multistage Clonal Expansion Model of Cancer , 2011, PLoS Comput. Biol..

[23]  J. Schiff The Laplace Transform: Theory and Applications , 1999 .

[24]  D. Kendall,et al.  On some modes of population growth leading to R. A. Fisher's logarithmic series distribution. , 1948, Biometrika.

[25]  M. E. J. Newman,et al.  Power laws, Pareto distributions and Zipf's law , 2005 .

[26]  S. Redner,et al.  Organization of growing random networks. , 2000, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  M. Meerschaert Regular Variation in R k , 1988 .

[28]  James D. Murray Mathematical Biology: I. An Introduction , 2007 .

[29]  Stefano Longhi Spiral waves in optical parametric oscillators , 2001 .

[30]  Philippe Flajolet,et al.  Analytic Combinatorics , 2009 .

[31]  S. Tavaré The birth process with immigration, and the genealogical structure of large populations , 1987, Journal of mathematical biology.

[32]  Pierre Baldi,et al.  The fixed-size Luria-Delbruck model with a nonzero death rate. , 2007, Mathematical biosciences.

[33]  Johannes G. Reiter,et al.  The molecular evolution of acquired resistance to targeted EGFR blockade in colorectal cancers , 2012, Nature.

[34]  B. Houchmandzadeh General formulation of Luria-Delbrück distribution of the number of mutants , 2015, bioRxiv.

[35]  P. Flajolet,et al.  Analytic Combinatorics: RANDOM STRUCTURES , 2009 .

[36]  J. Kiefer,et al.  An Introduction to Stochastic Processes. , 1956 .

[37]  N. Rashevsky,et al.  Mathematical biology , 1961, Connecticut medicine.

[38]  William Feller,et al.  An Introduction to Probability Theory and Its Applications , 1951 .

[39]  Kenneth Dixon,et al.  Introduction to Stochastic Modeling , 2011 .

[40]  Peter Nijkamp,et al.  Accessibility of Cities in the Digital Economy , 2004, cond-mat/0412004.

[41]  David R. Anderson,et al.  Model Selection and Inference: A Practical Information-Theoretic Approach , 2001 .

[42]  Marc J. Williams,et al.  Identification of neutral tumor evolution across cancer types , 2016, Nature Genetics.

[43]  David A Kessler,et al.  Resistance to chemotherapy: patient variability and cellular heterogeneity. , 2014, Cancer research.

[44]  Mikko Alava,et al.  Branching Processes , 2009, Encyclopedia of Complexity and Systems Science.

[45]  C. A. Coulson,et al.  The distribution of the numbers of mutants in bacterial populations , 1949, Journal of Genetics.

[46]  M. Newman Power laws, Pareto distributions and Zipf's law , 2005 .

[47]  P. L. Krapivsky,et al.  Exact solution of a two-type branching process: clone size distribution in cell division kinetics , 2009, 0908.0484.

[48]  F. Michor,et al.  Evolution of acquired resistance to anti-cancer therapy. , 2014, Journal of theoretical biology.

[49]  Michael R. Frey,et al.  An Introduction to Stochastic Modeling (2nd Ed.) , 1994 .

[50]  C. Maley The Evolutionary Dynamics of Cancer , 2000 .

[51]  D. Kendall Birth-and-death processes, and the theory of carcinogenesis , 1960 .

[52]  Mark E. J. Newman,et al.  Power-Law Distributions in Empirical Data , 2007, SIAM Rev..

[53]  P. Krapivsky,et al.  Exact solution of a two-type branching process: models of tumor progression , 2011, 1105.1157.

[54]  Martin A Nowak,et al.  Evolution of Resistance During Clonal Expansion , 2006, Genetics.

[55]  Richard Durrett,et al.  Evolution of Resistance and Progression to Disease during Clonal Expansion of Cancer , 2009 .

[56]  Q Zheng,et al.  Progress of a half century in the study of the Luria-Delbrück distribution. , 1999, Mathematical biosciences.

[57]  Martin A Nowak,et al.  Timing and heterogeneity of mutations associated with drug resistance in metastatic cancers , 2014, Proceedings of the National Academy of Sciences.

[58]  Marco Zaider,et al.  A stochastic model for the sizes of detectable metastases. , 2006, Journal of theoretical biology.

[59]  A. R. Whiting X-ray induced visible mutations in Habrobracon oocytes. , 1963, Genetics.