## A clean calculation of the Luria-Delbruck mathematics

Last updated: May 27, 2006

All biology students are made to read Luria and Delbruck's paper "Mutations of bacteria from virus sensitivity to virus resistance," which begins with a series of ad hoc probability calculations. When the paper was written in 1943, they represented the best tools available in probability. Today they are quaint to the point of being incomprehensible. Here is the whole calculation from the paper in modern form:

Assume that we have a population of bacteria growing exponentially and deterministically with growth rate 1, so the size of the population at time is , where is the size at time 0. These bacteria mutate to a state resistant to bacteriophage infections randomly and uniformly in time, but assume that the resistant bacteria are otherwise identical to their progenitors and that their numbers are always negligible compared to the whole population.

Define as the number of bacteria that mutate in the interval of time . is Poisson with mean , where is the mutation rate per cell. The rest of the calculation depends on two facts about :

Let be the number of mutant cells at time which originated from cells mutating in the interval . Then

Finally, let be the total number of mutant cells at time . Obviously,

Now we calculate the means and variances of .

= { definition of , and commutes with integration }

= { definition of , and is linear }

= { definition of ; algebra }

And the variance,

= { definition of }

= { pull through integral and expectation }

= { write out square of integral, renaming one dummy variable }
= { pull one integral inside the other }
= { expectation commutes with integrals }

= { property 2 of }

= { collapse the delta function, evaluate the integral }

The ratio of variance to mean is

and we can ignore the if we are looking very far out in time. The point of the Luria-Delbruck paper is to test whether mutants arise according to the process just described, at random before the addition of phage, or at random in response to the addition of phage. If it were in response to the addition of phage, the incidence should be Poisson distributed, with variance equal to mean. So we have a clean way to measure the difference between the two.

In the original paper, Delbruck prefers to change the lower limit of integration from 0 to the expected time of the first mutation. This produces a statistic which is on average more narrowly distributed, but is biased.

Fred Ross
May 2006
New York, NC