TY - JOUR

T1 - Large population solution of the stochastic luria-delbruck evolution model

AU - Kessler, David A.

AU - Levine, Herbert

PY - 2013/7/16

Y1 - 2013/7/16

N2 - Luria and Delbruck introduced a very useful and subsequently widely adopted framework for quantitatively understanding the emergence of new cellular lineages. Here, we provide an analytical treatment of the fully stochastic version of the model, enabled by the fact that population sizes at the time of measurement are invariably very large and mutation rates are low. We show that the Lea-Coulson generating function describes the 'inner solution,' where the number of mutants is much smaller than the total population. We find that the corresponding distribution function interpolates between a monotonic 1=m m+1 decrease at relatively small populations, (compared with the inverse of the mutation probability), whereas it goes over to a Lévy a-stable distribution in the very large population limit. The moments are completely determined by the outer solution, and so are devoid of practical significance. The key to our solution is focusing on the fixed population size ensemble, which we show is very different from the fixed time ensemble due to the extreme variability in the evolutionary process.

AB - Luria and Delbruck introduced a very useful and subsequently widely adopted framework for quantitatively understanding the emergence of new cellular lineages. Here, we provide an analytical treatment of the fully stochastic version of the model, enabled by the fact that population sizes at the time of measurement are invariably very large and mutation rates are low. We show that the Lea-Coulson generating function describes the 'inner solution,' where the number of mutants is much smaller than the total population. We find that the corresponding distribution function interpolates between a monotonic 1=m m+1 decrease at relatively small populations, (compared with the inverse of the mutation probability), whereas it goes over to a Lévy a-stable distribution in the very large population limit. The moments are completely determined by the outer solution, and so are devoid of practical significance. The key to our solution is focusing on the fixed population size ensemble, which we show is very different from the fixed time ensemble due to the extreme variability in the evolutionary process.

UR - http://www.scopus.com/inward/record.url?scp=84880347707&partnerID=8YFLogxK

U2 - 10.1073/pnas.1309667110

DO - 10.1073/pnas.1309667110

M3 - ???researchoutput.researchoutputtypes.contributiontojournal.article???

C2 - 23818583

SN - 0027-8424

VL - 110

SP - 11682

EP - 11687

JO - Proceedings of the National Academy of Sciences of the United States of America

JF - Proceedings of the National Academy of Sciences of the United States of America

IS - 29

ER -