The Luria–Delbrück experiment (1943) (also called the Fluctuation Test) demonstrates that in bacteria, genetic mutations arise in the absence of selection, rather than being a response to selection. Therefore, Darwin's theory of natural selection acting on random mutations applies to bacteria as well as to more complex organisms. Max Delbrück and Salvador Luria won the 1969 Nobel Prize in Physiology or Medicine in part for this work.
By the 1940s the ideas of inheritance and mutation were generally accepted, though the role of DNA as the hereditary material had not yet been established. It was thought that bacteria were somehow different and could develop heritable genetic mutations depending on the circumstances they found themselves: in short, was the mutation in bacteria pre-adaptive (pre-existent) or post-adaptive (directed adaption)? Luria in particular was obsessed with this idea and was determined to test it. He conceived the experiment at a faculty dance at Indiana University while watching a slot machine.
In their experiment, Luria and Delbrück inoculated a small number of bacteria (Escherichia coli) into separate culture tubes. After a period of growth, they plated equal volumes of these separate cultures onto agar containing the T1 phage (virus). If resistance to the virus in bacteria were caused by an induced activation in bacteria i.e. if resistance were not due to heritable genetic components, then each plate should contain roughly the same number of resistant colonies. Assuming a constant rate of mutation, Luria hypothesized that if mutations occurred after and in response to exposure to the selective agent, the number of survivors would be distributed according to a Poisson distribution with the mean equal to the variance. This was not what Delbrück and Luria found: Instead the number of resistant colonies on each plate varied drastically: the variance was considerably greater than the mean.
Luria and Delbrück proposed that these results could be explained by the occurrence of a constant rate of random mutations in each generation of bacteria growing in the initial culture tubes. Based on these assumptions Delbrück derived a probability distribution (now called the Luria–Delbrück distribution) that gives a relationship between moments consistent with the experimentally obtained values. The distribution that follows from the directed adaptation hypothesis (the Poisson distribution) predicted moments inconsistent with the data. Therefore, the conclusion was that mutations in bacteria, as in other organisms, are random rather than directed.
The results of Luria and Delbrück were confirmed in more graphical, but less quantitative, way by Newcombe. Newcombe incubated bacteria in a Petri dish for a few hours, then replica plated it onto two new Petri dishes treated with phage. The first plate was left unspread, and the second plate was then respread, that is, bacterial cells were moved around allowing single cells in some colony to form their own new colonies. If colonies contained resistant bacterial cells before entering into contact with the phage virus, one would expect that some of these cells would form new resistant colonies on the respread dish and so to find a higher number of surviving bacteria there. When both plates were incubated for growth, there were actually as much as 50 times greater number of bacterial colonies on the respread dish. This showed that bacterial mutations to virus resistance had randomly occurred during the first incubation. Once again, the mutations occurred before selection was applied.
More recently, the results of Luria and Delbrück were questioned by Cairns and others, who studied mutations in sugar metabolism as a form of environmental stress. Some scientists suggest that this result may have been caused by selection for gene amplification and/or a higher mutation rate in cells unable to divide. Others have defended the research and propose mechanisms which account for the observed phenomena consistent with adaptive mutagenesis.
This distribution appears to have been first determined by Haldane. An unpublished manuscript was discovered in 1991 at University College London describing this distribution. The derivation is different but the results are difficult to compute without the use of a computer.
Description of the test
A small number of cells are used to inoculate parallel cultures in a non-selective medium. The cultures are grown to saturation to obtain equal cell densities. The cells are plated onto selective media to obtain the number of mutants (r). Dilutions are plated onto rich medium to calculate the total number of viable cells ( Nt ). The number of mutants that appear in the saturated culture is a measure of both the mutation rate and when the mutants arise during the growth of the culture: mutants appearing early in the growth of the culture will propagate many more mutants than those that arise later during growth. These factors cause the frequency ( r / Nt ) to vary greatly, even if the number of mutational events ( m ) is the same. Frequency is not a sufficiently accurate measure of mutation and the mutation rate (m / Nt) should always be calculated.
The estimation of the mutation rate (μ) is complex. Luria and Delbruck estimated this parameter from the mean of the distribution but this estimator was subsequently shown to be biased.
- r = median number of colonies on one plate containing the indicator (e.g. rifampicin, sodium chlorate, streptomycin)
- m = a variable which will be varied, corresponds to the mutations/culture
- The value of the variable m is adjusted until the total value of the equation is close to 0. Then the mutation rate (probablitity of a mutation/cell/division or generation) can be calculated as one of three formulae:
- where Nt is the median of the number of viable cells on a non-indicator plate (often LB agar with no additive)
- The choice of which formula to use depends on at which stage in the cell division that the mutations are expected to occur. 
This method has since been improved on but these more accurate methods are complex. The Ma-Sandri-Sarkar maximum likelihood estimator is currently the best known estimator. A number of additional methods and estimates from experimental data have been described.
Two web-applications for the calculation of the mutation rate are freely available: Falcor  and bz-rates. Bz-rates implements a generalized version of the Ma-Sandri-Sarkar maximum likelihood estimator that can take into account the relative differential growth rate between mutant and wild-type cells as well as a generating function estimator that can estimate both the mutation rate and the differential growth rate. A worked example is shown in this paper by Jones et al.
In all these models the mutation rate (μ) and growth rate (β) were assumed to be constant. The model can be easily generalized to relax these and other constraints. These rates are likely to differ in non experimental settings. The models also require that Nt μ >> 1 where Nt is the total number of organisms. This assumption is likely to hold in most realistic or experimental settings.
Luria and Delbrück estimated the mutation rate from the equation
where β is the cellular growth rate, n0 is the initial number of bacteria in each culture, t is the time, and
where Ns is the number of cultures without resistant bacteria and N is the total number of cultures.
Lea and Coulson's model differed from the original in that they considered a collection of independent Yule processes (a filtered Poisson process). Numerical comparisons of these two models with realistic values of the parameters has shown that they differ only slightly. The generating function for this model was found by Bartlett in 1978 and is
where μ is the mutation rate (assumed to be constant), φ = 1 − e−βt with β as the cellular growth rate (also assumed to be constant) and t as the time.
The determination of μ from this equation has proved difficult but a solution was discovered in 2005. Differentiation of the generating function with respect to μ allows the application of the Newton–Raphson method which together with the use of a score function allows one to obtain confidence intervals for μ.
The mechanism of resistance to the phage T1 appears to have been due to mutations in the fhuA gene - a membrane protein that acts as the T1 receptor. The tonB gene product is also required for infection by T1. The FhuA protein is actively involved in the transport of ferrichrome, albomycin and rifamycin. It also confers sensitivity to microcin J25 and colicin M and acts as a receptor for the phages T5 and phi80 as well as T1.
The FhuA protein has a beta-barrel domain (residues 161 to 714) that is closed by a globular cork domain (residues 1 to 160). Within the cork domain is the TonB binding region (residues 7 to 11). The large membrane spanning monomeric β-barrel domains have 22 β-strands of variable length, several of which extend significantly beyond the membrane hydrophobic core into the extracellular space. There are 11 extracellular loops numbered L1 to L11. The L4 loop is where the T1 phage binds.
- Luria SE (1984) A slot machine, a broken test tube: An autobiography. Harper & Row
- Zheng, Q. (1999). "Progress of a half century in the study of the Luria–Delbrück distribution". Mathematical Biosciences. 162 (1–2): 1–32. doi:10.1016/S0025-5564(99)00045-0. PMID 10616278.
- Zheng, Q. (2010). "The Luria-Delbrück distribution: early statistical thinking about evolution". Chance. 23: 15–18. doi:10.1007/s00144-010-0017-y.
- Luria, S. E.; Delbrück, M. (1943). "Mutations of Bacteria from Virus Sensitivity to Virus Resistance". Genetics. 28 (6): 491–511.
- Newcombe, H. B. (1949). "Origin of Bacterial Variants". Nature. 164 (4160): 150–151. Bibcode:1949Natur.164..150N. doi:10.1038/164150a0. PMID 18146850. S2CID 4119793.
- Cairns, J.; Overbaugh, J.; Miller, S. (1988). "The Origin of Mutants". Nature. 335 (6186): 142–145. Bibcode:1988Natur.335..142C. doi:10.1038/335142a0. PMID 3045565. S2CID 4304995.
- Slechta, E. S.; Liu, J.; Andersson, D. I.; Roth, J. R. (2002). "Evidence that selected amplification of a bacterial lac frameshift allele stimulates Lac(+) reversion (adaptive mutation) with or without general hypermutability". Genetics. 161 (3): 945–956. PMC 1462195. PMID 12136002.
- Foster, Patricia L. (2004). "Adaptive Mutation in Escherichia coli". Journal of Bacteriology. 186 (15): 4846–4852. doi:10.1128/jb.186.15.4846-4852.2004. PMC 451643. PMID 15262917.
- Sarkar, S (1991). "Haldane's solution of the Luria-Delbruck distribution". Genetics. 127 (2): 257–261. PMC 1204353. PMID 2004702.
- Hall, BM; Ma, CX; Liang, P; Singh, KK (2009). "Fluctuation analysis CalculatOR: a web tool for the determination of mutation rate using Luria-Delbruck fluctuation analysis". Bioinformatics. 25 (12): 1564–1565. doi:10.1093/bioinformatics/btp253. PMC 2687991. PMID 19369502.
- Lea, DE; Coulson, CA (1949). "The distribution of the numbers of mutants in bacterial populations". J Genet. 49 (3): 264–285. doi:10.1007/bf02986080. PMID 24536673. S2CID 30301620.
- Foster, Patricia L. (2006), Methods for Determining Spontaneous Mutation Rates, Methods in Enzymology, 409, Elsevier, pp. 195–213, doi:10.1016/s0076-6879(05)09012-9, ISBN 978-0-12-182814-1, PMC 2041832, PMID 16793403, retrieved 2020-10-19
- Zheng, Q (2000). "Statistical and algorithmic methods for fluctuation analysis with SALVADOR as an implementation". Math Biosci. 176 (2): 237–252. doi:10.1016/S0025-5564(02)00087-1. PMID 11916511.
- Rosche, WA; Foster, PL (2000). "Determining mutation rates in bacterial populations". Methods. 20 (1): 4–17. doi:10.1006/meth.1999.0901. PMC 2932672. PMID 10610800.
- Jones, ME; Thomas, SM; Rogers, A (1994). "Luria-Delbruk experiments: Design and analysis". Genetics. 136: 1209–1216.
- Houchmandzadeh, B. (2015). "General formulation of Luria-Delbrück distribution of the number of mutants". Phys. Rev. E. 92 (1): 012719. arXiv:1505.06108. Bibcode:2015PhRvE..92a2719H. doi:10.1103/PhysRevE.92.012719. PMID 26274214. S2CID 4834465.
- Zheng, Q (1999). "Progress of a half century in the study of the Luria–Delbrück distribution". Mathematical Biosciences. 162 (1–2): 1–32. doi:10.1016/s0025-5564(99)00045-0. PMID 10616278.
- Bartlett M. (1978) An introduction to stochastic processes. Cambridge University Press, Cambridge, 3rd edition
- Carvajal-Rodríguez, A. (2012). "Teaching the fluctuation test in silico by using mutate: A program to distinguish between the adaptive and spontaneous mutation hypotheses". Biochemistry and Molecular Biology Education. 40 (4): 277–283. doi:10.1002/bmb.20615. PMID 22807434. S2CID 22732741.
- Endriss, F; Braun, M; Killmann, H; Braun, V (2003). "Mutant analysis of the Escherichia coli FhuA protein reveals sites of FhuA activity". J Bacteriol. 185 (16): 4683–4692. doi:10.1128/jb.185.16.4683-4692.2003. PMC 166461. PMID 12896986.
- Killmann, H; Braun, M; Herrmann, C; Braun, V (2001). "FhuA barrel-cork hybrids are active transporters and receptors". J Bacteriol. 183 (11): 3476–3487. doi:10.1128/jb.183.11.3476-3487.2001. PMC 99646. PMID 11344156.