Luria–Delbrück experiment

The two possibilities tested by the Luria–Delbrück experiment. (A) If mutations are induced by the media, roughly the same number of mutants are expected to appear on each plate. (B) If mutations arise spontaneously during cell divisions prior to plating, each plate will have a highly variable number of mutants.

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.

History

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 (James Watson's PhD advisor) 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.^[1]

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^[2][3]) 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.^[4]

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.^[5]

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.^[6] 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.^[7] Others have defended the research and propose mechanisms which account for the observed phenomena consistent with adaptive mutagenesis.^[8]

This distribution appears to have been first determined by Haldane.^[9] An unpublished manuscript was discovered in 1991 in 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.^[10] 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 ( N_t ). 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 / N_t ) 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 / N_t) should always be calculated.

The estimation of the mutation rate (m) is complex. Luria and Delbruck estimated this parameter from the mean of the distribution but this estimator was subsequently shown to be biased. The method of the median was introduced in 1949.^[11] This method is based on the equation

${\frac {r}{m}}-\ln(m)-1.24=0.$

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.^[12] A number of additional methods and estimates from experimental data have been described.^[13]

Two web-applications for the calculation of the mutation rate are freely available: Falcor ^[10] 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.^[14]

Distribution

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.^[15] These rates are likely to differ in non experimental settings. The models also require that N_t m >> 0 where N_t is the total number of organisms. This assumption is likely to hold in most realistic or experimental settings.

Luria and Delbrück^[4] estimated the mutation rate from the equation

${\displaystyle \mu ={\frac {\beta \ln(\rho )}{n_{0}[1-\exp(\beta t)]}}}$

where β is the cellular growth rate, n₀ is the initial number of bacteria in each culture, t is the time, and

${\displaystyle \rho ={\frac {N_{s}}{N}}}$

where N_s is the number of cultures without resistant bacteria and N is the total number of cultures.

Lea and Coulson's model^[11] 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.^[16] The generating function for this model was found by Bartlett in 1978^[17] and is

$G(z,\mu ,\phi )=\exp \left({\frac {\mu }{\phi }}\left({\frac {1}{z}}-1\right)\log \left(1-\phi z\right)\right)$

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 μ.

Molecular biology

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.^[18] 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.^[19] 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).^[20] 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.

References

1. Luria SE (1984) A slot machine, a broken test tube: An autobiography. Harper & Row
2. 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.
3. Zheng, Q. (2010). "The Luria-Delbrück distribution: early statistical thinking about evolution". Chance. 23: 15–18. doi:10.1007/s00144-010-0017-y.
4. Luria, S. E.; Delbrück, M. (1943). "Mutations of Bacteria from Virus Sensitivity to Virus Resistance". Genetics. 28 (6): 491–511.
5. Newcombe, H. B. (1949). "Origin of Bacterial Variants". Nature. 164 (4160): 150–151. Bibcode:1949Natur.164..150N. doi:10.1038/164150a0.
6. 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.
7. 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.
8. 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.
9. Sarkar, S (1991). "Haldane's solution of the Luria-Delbruck distribution". Genetics. 127 (2): 257–261. PMC 1204353. PMID 2004702.
10. 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.
11. Lea, DE; Coulson, CA (1949). "The distribution of the numbers of mutants in bacterial populations". J Genet. 49: 264–285. doi:10.1007/bf02986080.
12. Zheng, Q (2000). "Statistical and algorithmic methods for fluctuation analysis with SALVADOR as an implementation". Math Biosci. 176 (2): 237–252.
13. Rosche, WA; Foster, PL (2000). "Determining mutation rates in bacterial populations". Methods. 20 (1): 4–17. doi:10.1006/meth.1999.0901.
14. Jones, ME; Thomas, SM; Rogers, A (1994). "Luria-Delbruk experiments: Design and analysis". Genetics. 136: 1209–1216.
15. B. Houchmandzadeh (2015), General formulation of Luria-Delbrück distribution of the number of mutants, Phys. Rev. E 92:012719
16. Zheng Q (1999) Progress of a half century in the study of the Luria–Delbrück distribution. Mathematical Biosciences 162 (1–2) 1–32
17. Bartlett M. (1978) An introduction to stochastic processes. Cambridge University Press, Cambridge, 3rd edition
18. 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.
19. 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.
20. 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.

External links

• On Mutation lab
• Profiles in Science: The Salvador E. Luria Papers Information on Salvador Luria from the National Library of Medicine