Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, December 2005, p. 5081-5091, Vol. 49, No. 12
0066-4804/05/$08.00+0 doi:10.1128/AAC.49.12.5081-5091.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Patricia Geli,2,3,4,
Dan I. Andersson,5 and
Otto Cars1*
Antibiotic Research Unit, Department of Medical Sciences, Clinical Bacteriology and Infectious Diseases, Uppsala University, Uppsala, Sweden,1 Department of Epidemiology, Swedish Institute for Infectious Disease Control (SMI), Solna, Sweden,2 Division of Mathematical Statistics, Department of Mathematics, Stockholm University, Stockholm, Sweden,3 Stockholm Group for Epidemic Modeling (S-GEM), Stockholm, Sweden,4 Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden5
Received 7 January 2005/ Returned for modification 5 May 2005/ Accepted 1 October 2005
|
|
|---|
|
|
|---|
Among the most frequently used antibiotics are ß-lactams, such as penicillins and cephalosporins (26). ß-lactams interrupt the synthesis of the bacterial cell wall by forming a covalently bound complex with penicillin-binding proteins (PBPs), which are enzymes important in the final process of cell wall formation in bacteria (43, 44). The ability to produce TEM-ß-lactamases is the main mechanism for ß-lactam resistance in enteric gram-negative bacteria. The ß-lactamase enzymes inactivate penicillins and other ß-lactams by hydrolyzing the ß-lactam ring (24). The first plasmid-mediated ß-lactamase enzyme, TEM-1, was described shortly after the introduction of ampicillin for clinical use (6). Horizontal transfer of resistance genes led to a rapid interspecies spread of resistance, and today, TEM-1 is the most prevalent plasmid-mediated ß-lactamase found in gram-negative organisms (40, 41, 47). Antibiotic pressure has selected for over 130 TEM-1 ß-lactamase mutants with expanded hydrolytic capacities and activities against a variety of ß-lactam antibiotics, including monobactams, carbapenems, and extended-spectrum cephalosporins (25, 42). TEM-12 is a descendant of the TEM-1 enzyme and differs in a single substitution of arginine for serine at position 164 (22, 45). As a monomutated ß-lactamase, TEM-12 expresses an only slightly increased hydrolytic activity for cefotaxime. The most efficient TEM variants, which confer high-level resistance to cefotaxime, diverge from the native enzyme in several amino acids (4).
The growth of resistant subpopulations during treatment of a patient initially infected with susceptible bacteria presents an important problem. A number of in vitro studies have examined the effect of different dosing regimens in order to suppress the resistant subpopulations (1, 10, 23, 31, 35). A study by Negri et al. (31) revealed that low antibiotic concentrations can affect the selection of bacterial populations that show only small differences in susceptibility. Their work was based on a competition assay with Escherichia coli strains expressing different plasmid-borne variants of TEM-ß-lactamase enzymes. Negri detected a range of cefotaxime concentrations, a selective window, at which the selection of the strain with highest level of resistance was most intense. The experiments, however, were performed with static antibiotic concentrations in culture. Since antimicrobial therapy usually results in fluctuating drug concentrations in the patient, the selection process during treatment can be expected to differ from that in models with static antibiotic concentrations. Therefore, the outcome of the static model is difficult to apply on an individual patient level. In our study using a kinetic model, the selective window (SW) was defined as the concentration range between the MICs of two strains.
The purposes of this study were as follows: (i) use an in vitro kinetic model to study the selection of cefotaxime-resistant E. coli for different time periods within the SW, and (ii) construct a general mathematical model that describes the expected changes in the bacterial population as a function of pharmacokinetic parameters. Our hypothesis was that a longer time within the SW would increase the selection of the more resistant strain, when two strains were competing in the model. Unlike earlier studies examining the efficacy of various dosing regimens in preventing the emergence of resistance, this model incorporates the selection of both preexisting and newborn mutants and any potential post-MIC effect (PME). The PME is the period when regrowth is delayed even after antibiotic concentrations have fallen below the MIC (13, 21) and, like the in vivo postantibiotic effect, includes the effects of subinhibitory concentrations (9, 27). The model provides a convenient theoretical framework to understand experimental data and a theoretical basis for optimal dosing regimens, in order to maintain efficacy while simultaneously preventing the emergence of resistance.
|
|
|---|
Growth rates were determined for strains REL606(pBGTEM-1) and REL607(pBGTEM-12) separately in tubes. These strains will be referred to as TEM-1 and TEM-12 in this paper. Competition experiments were performed in the in vitro kinetic model (described below) with an initial 1:1 ratio of the competing strains. In addition, the competition experiments were performed with an inverse pair, REL607(pBGTEM-1) and REL606(pBGTEM-12), to confirm the neutrality of the plasmids and the arabinose marker of the host bacteria.
The plasmids have kanamycin resistance as a selective marker; hence, the strains were maintained on Columbia agar (Acumedia Manufacturers, Inc., Baltimore, MD) plates supplemented with 30 µg kanamycin/ml. The liquid medium used for bacterial growth was Mueller-Hinton broth (Difco Laboratories, Detroit, MI), and the solid medium in the assays was tetrazolium-arabinose indicator agar (18). The bacteria were grown at 35°C, and liquid cultures were incubated without shaking.
Antimicrobial agents. Cefotaxime powder was obtained from Aventis (Stockholm, Sweden) and was dissolved in 1 ml sterile distilled water to a concentration of 10 mg/ml. Fresh stock solutions were prepared on the day of use and diluted in Mueller-Hinton broth.
Susceptibility testing. The MICs of cefotaxime for the native strains were determined by a macrodilution technique according to CLSI (formerly NCCLS) standards (30) and were done in triplicate on separate occasions. The MICs for the strains containing TEM-1 and TEM-12 were 0.016 and 0.063 µg/ml, respectively, and these MIC values were used for the study design.
To detect the appearance of novel resistant mutants during exposure to cefotaxime, colonies were taken from the 24-h samples and analyzed with Etest on Columbia agar plates according to the instructions by the manufacturer (AB Biodisk, Solna, Sweden). The Etest method resulted in slightly lower MICs for the parental strains, 0.012 µg/ml for TEM-1 and 0.032 µg/ml for TEM-12, than with the macrodilution technique.
Determination of antibiotic concentrations.
The initial cefotaxime concentrations in the in vitro kinetic experiments were determined with a microbiological agar diffusion method. Plates with tryptone-glucose agar, pH 7.4, were seeded with a standardized inoculum of Escherichia coli MB3804. Antibiotic standards and samples from the experiments were applied to agar wells at a volume of 30 µl, and the plates were incubated overnight at 35°C. All assays were made in triplicate and the correlation coefficient for the standard curves was always
0.99.
In vitro kinetic model. The in vitro kinetic model used in this study has been described earlier (12, 21). It consists of a spinner flask (110 ml) with an open bottom that was placed on a holder with an outlet connected to a pump (P-500; Pharmacia Biotech, Uppsala, Sweden). A filter membrane with a pore size of 0.45 µm was supported by a metal rack between the flask and the holder, impeding the dilution of bacteria. A magnetic stirrer ensured a homogenous mixing of the culture and prevented membrane pore blockage. The spinner flask had two side arms: one with a silicone membrane inserted to enable repeated sampling and another connected to plastic tubing from a vessel containing fresh medium. The medium was drawn from the culture vessel, through the filter, at a given rate by the pump. Fresh medium was sucked into the flask at the same rate by the negative pressure built up inside. Antibiotic added to the flask was diluted according to the first-order kinetics according to equation 3 in the mathematical model. The apparatus was placed in a thermostatic room at 35°C during the experiment.
Study design: selective windows. Competition assays were performed with various times within the SW, i.e., time periods when the concentration of cefotaxime is below the MIC for TEM-12 but above the MIC for TEM-1. The flask was prepared with broth and the desired initial antibiotic concentration (Cmax, Table 1) and was installed in the thermostatic room (35°C). Bacteria from 6- to 7-h broth cultures were added to the flask to create a culture mixture of TEM-1 and TEM-12 at a proportion of 99:1. The initial bacterial concentrations of TEM-1 and TEM-12 ß-lactamase-producing strains were 105 CFU/ml and 103 CFU/ml, respectively. The time that the concentrations exceeded the MIC (T > MIC) for the TEM-12 strain was 2 h in all SWs, while T > MIC for the TEM-1-producing strain was varied. The elimination half-life (T1/2) in the kinetic model was adjusted accordingly and, if needed, changed during the experiments to obtain SWs of 1, 2, 4, 8, and 12 h (Table 1 and Fig. 1), and the experiments were run for 24 h. Samples of 200 to 400 µl were withdrawn at different time points and treated with penicillinase type IV (Sigma-Aldrich Chemie GmbH, Steinheim, Germany) for 20 min to prevent antibiotic carryover. Dilutions of the samples were then seeded on tetrazolium-arabinose indicator plates, and after 24 h at 35°C, the pink and red colonies were counted. The limit of detection for viable counts was 10 CFU/ml. The strains were easily discriminated in all experiments except for the 24-h sample in two SWs with increased Cmax where there was heavy growth of TEM-12 (see Results). The experiments were repeated five times except SW 2 h, for which 10 separate experiments were performed for estimation of parameters in the statistical model. MIC determinations were performed with Etest as described in Materials and Methods.
|
View this table: [in a new window] |
TABLE 1. Cmax and T1/2 values
|
![]() View larger version (16K): [in a new window] |
FIG. 1. Concentration profiles of cefotaxime for five selective windows (1, 2, 4, 8, and 12 h) with low Cmax. MICP1 and MICP12 indicate the MICs of the parental populations TEM-1 and TEM-12, respectively.
|
Characterization of high-level cefotaxime-resistant mutants.
The MICs of cefotaxime, chloramphenicol, and tetracycline were determined with Etest for parental and mutated strains TEM-1 and TEM-12, as well as for E. coli MG1655 and E. coli LM201 (ompF<>FRT; derived in Escherichia coli MG1655, the
ompF has been generated by homologous recombination technology). For PCR amplification and DNA sequencing of ompF, DNA was prepared from parental and mutated strains TEM-1 and TEM-12 using the Wizard genomic DNA purification kit (Promega, Madison, WI). The primer sequences used for PCR and sequencing were constructed from the ompF gene of E. coli K12: 1F (5'-CGTGAGATTGCTCTGGAAGG-3'), 3R (5'-CTCAACCTCTTGGCAACGGTA-3'), 2F (5'-TCGTACTTCAGACCAGTAGC-3'), 5R (5'-ACGGTGAAAACAGTTACGGT-3'), 4F (5'-ATTGATTTGAGTTTCCCCTTTA-3'), and 6R (5'-TGACGGTGTTCACAAAGTTCC-3'). PCR was carried out in 20-µl volumes containing 1 µM forward and reverse primers, 0.5 µl DNA sample, and 5 mM Mg2+ (3 mM for primers 4F and 6R). The reactions were run in a GeneAmp PCR system 9700 (Applied Biosystems, Foster City, CA) and the following temperature profile was used: initial denaturation at 95°C for 30 s; 30 cycles of 95°C for 30 s, 54°C for 30 s, and 72°C for 30 s; and a final extension at 72°C for 7 min. For primers 4F and 6R, the annealing temperature was 53°C. The PCR products were purified with a GFX-DNA purification kit (Amersham Biosciences, Uppsala, Sweden). A BigDye Terminator v 1.1 cycle sequencing kit (Applied Biosystems) was used for sequencing, and the analysis was performed with an ABI 3100 Genetic Analyzer, a multicolored-fluorescence-based DNA analyzing system. The parental and mutated E. coli strains TEM-1 and TEM-12 were also tested for organic solvent tolerance as previously described by Komp Lindgren et al. (16).
|
|
|---|
Selective windows. The mean initial cefotaxime concentrations were within 10% of expected values (coefficient of variation, 11%). When TEM-1 and TEM-12 were mixed at a proportion of 99:1 in the in vitro kinetic model and challenged with cefotaxime to obtain different times within the SW, an increase in the proportion of the TEM-12-producing strain was observed in the SWs for 1, 2, 4, and 8 h (Fig. 2A to D, left panels). Since regrowth of both strains was apparent already after 12 h, results only up to this time point are shown in the graphs. In the first three SWs (1, 2, and 4 h) there was a clear dominance of TEM-12 but, unexpectedly, the selection of TEM-12 appeared to be less effective in the SW of 8 h, and in the SW of 12 h, TEM-1 was selected (Fig. 2E, left panel).
![]() View larger version (16K): [in a new window] |
FIG.2. Competition assays with E. coli strains TEM-1 and TEM-12 exposed to cefotaxime in the in vitro kinetic model. Five selective windows were investigated: 1 h (A), 2 h (B), 4 h (C), 8 h (D), and 12 h (E). In the first series (low Cmax; left panels) T > MIC was varied for TEM-1 (3, 4, 6, 10, and 14 h), and fixed (2 h) for TEM-12. Each graph displays means of 5 experiments with the exception of SW 2 h, which shows means based on 10 experiments. The bars represent standard deviations. In the second series of SWs (high Cmax; right panels) the cefotaxime doses were four times higher to prevent the emergence of high-level resistant mutants. T > MIC was varied for strain TEM-1 (4, 5, 7, 11, and 15 h), and fixed (3 h) for strain TEM-12. Each graph displays the means of two experiments. Solid line, strain TEM-12; dashed line, strain TEM-1.
|
Selective windows with increased concentration of cefotaxime. To minimize the selection of high-level resistant mutants, experiments were performed with increased Cmax (Fig. 2A to E, right panels). In these experiments, growth of strain TEM-1 was reduced and possibly prevented in the SW of 8 and 12 h. Since the TEM-12-producing bacteria were in dominance, potential colonies of the TEM-1 strain could not be separated in the mixed population. Thus, they were scored as zero growth. No increase in MICs was seen for strain TEM-1 colonies except in one of the two 4-h SWs (MIC = 0.094). With a high antibiotic concentration, the growth of newly formed mutants of TEM-1-producing bacteria was prevented. As a result, selection of the TEM-12-producing strain was increased in the SW of 8 and of 12 h. Bacterial regrowth of TEM-12 was noted after 12 h in these two SWs (Fig. 2D to E, right panels), and at 24 h, TEM-12 had grown more than 7 log CFU, while TEM-1 was undetectable. With even higher concentrations of cefotaxime (4 µg/ml) both E. coli strains could be completely eliminated (data not shown).
Characterization of high-level cefotaxime-resistant mutants appearing in the competition experiments. The high-level cefotaxime-resistant mutants that appeared showed MICs of chloramphenicol and tetracycline that were four times higher than for the parental strains. This finding suggested that cefotaxime, chloramphenicol, and tetracycline resistance were caused by an inactivation of a transport function or activation of an efflux system. For example, an ompF mutation could cause the cefotaxime-resistant phenotype (31). To examine this possibility, the MICs of the high-level cefotaxime-resistant strains were compared with the MICs for two isogenic E. coli strains, one wild type and one with a deletion in ompF. However, the defined ompF mutation had a much smaller effect on the MICs for chloramphenicol and tetracycline than did those in our mutants, suggesting that the mutations were not in ompF. In addition, DNA sequencing of the ompF gene in parental and mutated strains TEM-1 and TEM-12 revealed no changes. To further investigate the high-level cefotaxime resistance, the organic solvent tolerance was measured, a phenotype associated with overexpression of the transmembrane AcrAB-TolC multidrug efflux pump (46). However, none of the tested bacteria were tolerant to cyclohexane.
Mathematical model. Since there was no simple mathematical relationship between the selection and the time within the SW, a mathematical model of pharmacokinetics and bacterial population dynamics was constructed, with the aim to predict how the time within the SW affects the selection and/or the emergence of resistance. A technical description of the model has previously been published as a Master's thesis from Stockholm University (11).
(i) Pharmacokinetics.
The elimination of the initial concentration of antibiotics, Cmax, follows first-order kinetics with a elimination rate k(t) that changes at some time points depending on the experimental setting. Thus
![]() | (1) |
![]() | (2) |
![]() | (3) |
(ii) The basic model of population dynamics.
The parental strains are denoted by P. Changes in these populations will, in a simple model, depend only on the net growth rate
(t) (can be negative or positive), i.e., the rate of cell division minus kill rate due to antibiotics, as follows:
![]() | (4) |
Extensions of the model. (i) Appearance of mutants.
It was assumed that mutations occurred with a constant rate
during the whole experimental period. Thus, the number of parental bacteria decreases at the same rate as mutants occur. Mutants are denoted by M, which adds the following term to the basic model:
![]() | (5) |
(t), depends on the concentration at each time point and will explain some of the PME observed in the experiments. Figure 3 shows the rates that determine the population dynamics. The sums of parental and mutant populations, S, are the numbers of bacteria that are observed in the experiments.
![]() View larger version (50K): [in a new window] |
FIG. 3. Illustration of the relative growth rates 1(t) and 12(t) and mutation rates 1 and 12. The sums of the parental population and the mutated population, S1 (P1 + M1) and S12 (P12 + M12), represent the numbers of bacteria that will be observed during the experiments.
|
and are synthesized with a rate ß. This can be expressed as:
![]() | (6) |
on the changes in PBP, meaning that:
![]() | (7) |
min, can be negative under antibiotic pressure and is assumed to be present when all PBPs are saturated. Conversely, in the absence of antibiotics when no PBPs are saturated, a maximal growth rate of bacteria,
max, is present. In this case Q(t) = 1, and
(t) reduces to
(t) =
max and
=
max
min.
![]() View larger version (11K): [in a new window] |
FIG. 4. Modeling of PME. Antibiotic saturation and synthesis of PBPs depends on the initial concentration of drug and the half-life time. The binding rate of antibiotics to PBPs is denoted by and the synthesis rate of new PBPs by ß. Open circles represent PBPs without bound antibiotic and filled circles represent PBPs to which antibiotic is bound.
|
(t) equals zero when the concentration at time point t is equal to the MIC. Furthermore, if the number of unbound PBPs at the moment when the concentration has reached the MIC (CMIC) is denoted by QMIC, equation 7 yields the following relationship:
![]() | (8) |
![]() | (9) |
![]() | (10) |
|
View this table: [in a new window] |
TABLE 2. Parameter estimates of the modela
|
Note that the MICs were estimated as unknown parameters for all strains. Therefore these estimates of MICs can be compared to those measured by Etest (Table 2). The difference between measured data and estimates from the model is small, which provides a validation of the model.
Prediction of the selection and the proportion of mutants. Predictions of the outcomes of the parental strains indicate that the selection of parental TEM-12 increases with the time within the SW (1, 2, or 4 h), as long as the level of antibiotics is low enough to allow regrowth of this strain (Fig. 5). Thus, our theory holds for the parental strains. Since the proportion of mutants appears to increase with the time within the SW (Fig. 6) it will no longer be possible to see a relationship between selection and time within SW. The proportion of mutant TEM-12 organisms becomes high later than mutant TEM-1, which is due to the initially smaller inoculum of parental TEM-12.
![]() View larger version (8K): [in a new window] |
FIG. 5. Predictions of the outcome of the parental strains P1 and P12 for competition assays with E. coli strains TEM-1 and TEM-12 with two of the selective windows that were investigated: 1 h (left) and 12 h (right). Dashed line, P1; solid line, P12.
|
![]() View larger version (10K): [in a new window] |
FIG. 6. The proportions of mutants, M1/(P1 + M1) (left), and M12/(P12 + M12) (right), estimated from the predicted values for the five selective windows that were investigated: 1, 2, 4, 8, and 12 h. Dashed line, strain TEM-1; solid line, strain TEM-12.
|
![]() View larger version (12K): [in a new window] |
FIG. 7. Predicted PME for two cases yielding the same AUC: constant Cmax of 0.1 µg/ml and T1/2 varying from 0.5 to 5 h; and constant T1/2 of 0.5 h and Cmax varying from 0.1 to 0.8 µg/ml.
|
![]() View larger version (18K): [in a new window] |
FIG. 8. Predicted and experimental data from competition assays with E. coli strains TEM-1 and TEM-12 in the in vitro kinetic model for two selective windows. Shown are data for 1 h (A) and 12 h (B), with low Cmax (left panel) and high Cmax (right panel). *, strain TEM-1 observed data; dashed line, TEM-1 predicted data; , strain TEM-12 observed data; solid line, TEM-12 predicted data. The bars correspond to 95% predictive intervals.
|
|
|
|---|
Here, we present a pharmacodynamic model that, in contrast to previous models, includes rates for the occurrence of mutants and the saturation and synthesis of PBPs. Thus, the model can be used to predict the selection of both preexisting and newborn mutants as well as the effect of any potential PME. By reestimating parameters, the model can be used for predictions of pathogens and antibiotics other than Escherichia coli and cefotaxime.
From an earlier study, it was expected that the selection of the more resistant parental strain (TEM-12) would increase with the time within the selective window (31). In concordance with this hypothesis, our experimental data showed a high dominance of the TEM-12 strain in the 1-, 2-, and 4-h SWs. However, when SWs of 8 and 12 h were tested, the selection of TEM-12 actually decreased and in the SW of 12 h, the low-resistance parental strain (TEM-1) was selected. The lack of correlation between the strength of selection and the time within the SW was a result of emergence of newborn high-level resistant mutants and the influence of PME. The experiments demonstrated that the TEM-1 strain repeatedly attained a high-level resistance in the SWs of 8 and 12 h, and occasionally in the SW of 4 h. The mutants that appeared from the TEM-1 strain had MICs about 12 times the original MIC, a resistance level higher than for the parental TEM-12 strain. This explains why selection of the TEM-12 strain was decreased for longer times within the SW. Increasing the initial cefotaxime concentration four times prevented the growth of new mutants from TEM-1 in almost every experiment. This led to an increased selection of TEM-12 in SW of 8 and 12 h, which better concurs with the hypothesis that longer time within the selective window increases selection of the more resistant parental strain (TEM-12).
To validate the model, we compared the predicted outcome with observed data. The predictions were found satisfactory regarding both the selection of preexisting and newborn mutants and PME, despite the fact that the following simplifications were made. First, the mutants were assumed to appear with a constant rate, and not randomly. Second, there was no fitness cost associated with the high-level mutants. The latter simplification does not alter the prediction for which strain is selected, but it influences the amount of the selection and may explain the increased deviance between observed data and predicted outcomes in experiments. Finally, a fundamental difference between our model and antimicrobial treatment in patients is the lack of a host immune response in the model. Thus, in vivo, the antimicrobial efficacy and potency of drugs are assisted by immune factors. To increase the predictive power of future refined pharmacodynamic models, relevant immunological parameters should be included.
In a situation where antibiotic concentrations are declining and newborn high-level resistant mutants are formed, the outcome becomes complex and will strongly depend on the concentration that prevents growth of the most resistant strain. Obviously, if drug concentrations are continuously maintained above this concentration, no resistant mutants will appear. Importantly, even shorter time periods above this concentration can effectively prevent appearance of newborn mutants (see Results and Fig. 2, right panels). The issue of suppression of resistant subpopulations has also been addressed by Jumbe et al., who used a mathematical model to calculate an AUC/MIC ratio that amplified a mutant subpopulation in vivo as well as a ratio that prevented the emergence of resistance (15). Here we showed that if drug concentrations are lower and are maintained in the selective window, selection of the more resistant parental strain (TEM-12) as well as mutants from both parental strains will occur and increase with longer time within the SW. Using a fixed AUC, selection will be minimized using a high-dose, short-elimination half-life regimen rather than a low-dose, long half-life regimen. With regard to the PME, it can vary with the pharmacokinetic profile (8, 9, 21). In our model, with a fixed AUC the PME is slightly more pronounced, with a long half-life rather than a short one. Thus, although a long PME would allow extended dosing intervals with preserved efficacy, it would also promote resistance.
In conclusion, our experimental data and mathematical modeling show that in a dynamic competition between strains with different levels of resistance, the appearance of newborn high-level resistant mutants from the parental strains and the post-MIC effect can strongly affect the outcome of the selection. Thus, it is important that pharmacodynamic models incorporate biologically relevant parameters to allow more realistic predictions of resistance development.
|
|
|---|
, ß, CMICP, and CMICM. Since the complete probability mechanism was too complicated to specify a full likelihood, a quasi-likelihood approach (36) to achieve robust inference was used. This means that only the mean and variance functions have to be specified, instead of the probability structure. The conditional mean and variance functions of S(tj), given S(tj1), were achieved by assuming that S(tj) followed a branching process (39). That means, let S(tj), for j = 0, 1, 2, 3, 4, or 5, denote the number of bacteria at the time points t0, t1, t2, t3, t4, or t5 (0, 2, 4, 6, 8, or 12) and for the experimental setting with five time points (0, 2, 4, 6, 10. If each single bacterium in generation zero, S0, produces new bacteria with a mean µ and variance
2, the total number of offspring will depend on the size of the previous generation. Thus, the size of the jth generation is
![]() | (A1) |
Furthermore, the variance for the size of the jth generation is
![]() | (A2) |
1. Now, assume that there are n generations between time point tj and tj 1. Since each bacterium produces offspring with a mean µ in each generation, the conditional mean of the number of bacteria at time point tj given the number at time point tj 1 is
![]() | (A3) |
![]() | (A4) |
1. Set
![]() | (A5) |
![]() | (A6) |
![]() | (A7) |
This work was supported by the EU Fifth Framework Program (O.C. and D.I.A.) and the Swedish Research Council (D.I.A.).
Authors contributed equally to the paper. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»