- Research article
- Open Access
- Published:

# Evidence for a bimodal distribution of *Escherichia coli* doubling times below a threshold initial cell concentration

*BMC Microbiology*
**volume 10**, Article number: 207 (2010)

## Abstract

### Background

In the process of developing a microplate-based growth assay, we discovered that our test organism, a native *E. coli* isolate, displayed very uniform doubling times (τ) only up to a certain threshold cell density. Below this cell concentration (≤ 100 -1,000 CFU mL^{-1} ; ≤ 27-270 CFU well^{-1}) we observed an obvious increase in the τ scatter.

### Results

Working with a food-borne *E. coli* isolate we found that τ values derived from two different microtiter platereader-based techniques (i.e., optical density with growth time {=OD[t]} fit to the sigmoidal Boltzmann equation or time to calculated 1/2-maximal OD {=t_{m}} as a function of initial cell density {=t_{m}[C_{I}]}) were in excellent agreement with the same parameter acquired from total aerobic plate counting. Thus, using either Luria-Bertani (LB) or defined (MM) media at 37°C, τ ranged between 17-18 (LB) or 51-54 (MM) min. Making use of such OD[t] data we collected many observations of τ as a function of manifold initial or starting cell concentrations (C_{I}). We noticed that τ appeared to be distributed in two populations (bimodal) at low C_{I}. When C_{I} ≤100 CFU mL^{-1} (stationary phase cells in LB), we found that about 48% of the observed τ values were normally distributed around a mean (μ_{τ1}) of 18 ± 0.68 min (± σ_{τ1}) and 52% with μ_{τ2} = 20 ± 2.5 min (n = 479). However, at higher starting cell densities (C_{I}>100 CFU mL^{-1}), the τ values were distributed unimodally (μ_{τ} = 18 ± 0.71 min; n = 174). Inclusion of a small amount of ethyl acetate to the LB caused a collapse of the bimodal to a unimodal form. Comparable bimodal τ distribution results were also observed using *E. coli* cells diluted from mid-log phase cultures. Similar results were also obtained when using either an *E. coli* O157:H7 or a *Citrobacter* strain. When sterile-filtered LB supernatants, which formerly contained relatively low concentrations of bacteria(1,000-10,000 CFU mL^{-1}), were employed as a diluent, there was an evident shift of the two populations towards each other but the bimodal effect was still apparent using either stationary or log phase cells.

### Conclusion

These data argue that there is a dependence of growth rate on starting cell density.

## Background

Understanding the behavior of bacterial growth parameters (duration of lag phase, specific growth rate, and maximum cell density in stationary phase) under various environmental conditions is of some interest [1]. In particular, knowledge about growth parameter population distributions is needed in order to make better predictions about the growth of pathogens and spoilage organisms in food [1–3]. In fact, probability-based methods, such as microbial risk assessment [1], have to take into account the distribution of kinetic parameters in a population of cells [4]. There is a paucity of growth parameter distribution data because of the large number of data points required to obtain such results. The utilization of traditional microbiological enumeration methods (e.g., total aerobic plate count or TAPC) for such a body of work is daunting. For this reason various methods have been developed which enable more rapid observations related to one, or more, growth parameters. Recently, growth parameter distribution characterization has mainly focused on the duration of lag phase [4–8]. For instance, Guillier and co-workers studied the effects of various stress factors (temperature, starvation, salt concentration, etc.) on individual cell-based detection times in *Listeria monocytogenes*[5, 6]. Additionally, reporting on improved methods, various workers [4, 7, 8] have presented frequency distribution information concerning lag phase duration of individual bacterial cells (*Escherichia coli*, *L. monocytogenes*, and *Pseudomonas aeruginosa*) on solid media. However, similar population-based information on specific growth rate is lacking.

The findings presented herein developed from work associated with the attachment of various Gram-negative bacteria to anti-*Salmonella* and anti-*E. coli* O157 immunomagnetic beads or IMBs [9–11]. For these IMB investigations microplate (OD-based) MPN methods were utilized because of the low limits of bacterial detection [12, 13] necessary to characterize the non-specific attachment of background food organisms to various capture surfaces. Because of large inter-bacterial strain variability in the time requisite to reach a measurable level of turbidity, we found it necessary to characterize the growth rate and apparent lag time (time to 1/2-maximal OD or t_{m}) [12] of certain problematic organisms. Toward this end we began a routine investigation into the best microplate reader method to determine doubling time (τ). However, while performing this work we noticed that our test organism, a native *E. coli* isolate which non-specifically adheres to certain IMBs [11], seemed to display very uniform τ values only up to a certain threshold initial or starting cell density (C_{I}) beyond which we observed an obvious increase in the scatter. A larger number of observations were then made after various physiological perturbations (media used, growth phase, etc.) which have lead to the results discussed in this report.

## Results and Discussion

### Doubling Times from both TAPC and Microplate Observations

Table 1 shows analysis of variance data for τ calculated as described in the Methods Section from Optical Density with time (= OD[t]; **Eq. 1**) data, t_{m} as a function of C_{I} (= t_{m}[C_{I}]; **Eq. 6**), and total aerobic plate count with time (= TAPC[t]) on two different media at 37°C (C_{I} > 1,000 CFU mL^{-1}). These results indicate that doubling times derived from the aforementioned microplate techniques (i.e., OD[t] and t_{m}[C_{I}]) were in excellent agreement with τ values acquired from TAPC when using either Luria-Bertani (LB) or a defined minimal medium (MM) at 37°C. In these experiments τ varied 17 to 18 min (LB) or 51 to 54 min (MM) depending on media. The within-medium variation was not significant at even a 0.1 level (i.e., the probabilities of > 3.43 was 0.136 and >0.886 was 0.480). These results show that both microplate-based methods for measuring τ are equivalent to τ derived from TAPC. For low initial cell concentrations, the OD[t] method, as described in the Methods section, is obviously superior to t_{m}[C_{I}] since it makes no assumption about concentration dependence. However, for routine growth studies (e.g., antibiotic resistance) at a relatively high C_{I} the t_{m}[Φ_{I}] method (**Eq. 5**, Methods Section; Φ_{I} is the dilution factor used to make each C_{I}) for obtaining τ is preferable since t_{m} is easy to obtain without curve fitting albeit several dilutions need to be used.

The agreement between the *E. coli* τ from TAPC and microplate methods was somewhat unexpected inasmuch as solution agitation (i.e., oxygenation) of the media in each plate's wells would be less than that for solution agitation in either normal or baffled flasks which were used for the TAPC comparisons. However, we found (Fig. 1A, open symbols) that [O_{2}] levels in even highly agitated liquid *E. coli* cultures at 37°C dropped as much as 72% (LB, normal flask) with 200 RPM shaking while they were consuming approximately {4-6} × 10^{-18} moles O_{2} sec^{-1} CFU^{-1} (Fig. 1B). Even the baffled flask culture showed a drop in [O_{2}] of 40-57%. Simultaneously, no cultures (Fig. 1A, closed symbols) showed any perturbations in τ (~ 18 min); the 23 min τ seen with bubbling is probably greater due to evaporative cooling of the medium. Due to differences in both solution mixing and surface area-to-volume ratio, the [O_{2}] levels in microplate wells must be even lower than flask cultures at equivalent cell densities. Fig. 1 demonstrates that even at the lowest [O_{2}], the rates of growth were unaffected. Clearly, being a facultative anaerobe, *E. coli* is able to rapidly adjust to different levels of O_{2} with no apparent change in its specific growth rate, although the maximum cell density in stationary phase is usually greater in highly oxygenated samples by up to an order of magnitude.

### Effect of Initial or Starting CFU Concentration on τ

While performing studies related to comparing various assays for determining growth rate (Table 1), we noticed that our test organism, a nonpathogenic avian *E. coli* isolate, seemed to display uniform OD[t]-based τ values up to a threshold C_{I}, at which point there was an obvious increase in the observed τ scatter (Fig. 2). The main graph in Fig. 2 represents 653 measurements of τ derived from OD[t] data using **Eq. 1** (Methods Section) plotted as a function of C_{I} (diluted from stationary phase cells). When C_{I} > ca. 100 CFU mL^{-1}, τ was narrowly Gaussian-distributed (i.e., a unimodal distribution) with a total spread of ca. 16 to 20 min (n = 174 observations: μ_{τ} ± σ_{τ} = 17.6 ± 0.708 min using a single Gaussian distribution function: i.e., **Eq. 7** with α = 1 and β = 0; Methods Section). However, when C_{I} < ca. 100 CFU mL^{-1} there was a clear broadening in the range of observed τ values (ca. 10 to 34 min). At such low concentrations the CFUs per well should vary between 1 and 10 whereupon 44% of the wells should have 1 (± 1) CFU per well, 14% with 2 (± 1.4) CFUs per well, 8% with 3 (± 1.7) per well, 6% with 4 (± 2) per well, and 3% with between 5 (± 2.2) and10 (± 3.2) CFUs per well (assuming a Poisson distribution of CFU counts). The inset graph in Fig. 2 shows frequency of occurrence for all values of τ, which occur in the region of greatest scatter (C_{I}< 100 CFU mL^{-1}), with the best fit bimodal Gaussian distribution (**Eq. 7**) represented by the solid, black curve. The least squares bimodal distribution curve fit contains a narrow component (α ~0.48; μ_{τ1} ± σ_{τ1} = 18.0 ± 0.678 min) similar to the high cell concentration-associated unimodal distribution. Based upon area, there was also a nearly equivalent broad component (β ~ 0.52; μ_{τ2} ± σ_{τ2} = 19.9 ± 2.48 min). Each constituent of this bimodal distribution is shown as a solid, grey curve.

A similar increase in another growth parameter's scatter was also observed with the t_{m}[C_{I}]data at low C_{I} (Fig. 3) whereupon we saw that t_{m} values changed in a predictable way (e.g.,|∂t_{m}/∂Log_{2}C_{I}| = τ) up to CI ~ 100 - 1,000 CFU mL^{-1} at which point they began to show an obvious large deviation in t_{m} (between 6 and 11 hrs). These perturbations in t_{m} at low C_{I} confirm the τ observations because t_{m} is modulated, at least in part, by τ (**Eqs. 5-6**: all t_{m} & T-based equations are developed in the Methods Section) and therefore large deviations in τ (Fig. 2) should result in increased scatter in t_{m} as well. Working with stressed *Listeria monocytogenes*, Guillier and coworkers [5] observed numerous values of a lag time-related growth parameter with a similar asymmetric distribution pattern. Measuring the time of the first cell division in *E. coli* using a microscopic method, which should provide the true value of lag time, Niven and co-workers [8] were ableto make numerous observations (n = 434) which showed a very broad (μ_{T}~ 184 ± 45 min; our calculation assuming a unimodal distribution) asymmetric distribution. Asymmetry might be interpreted as weakly bimodal. Other workers [4] using a different method of observation showed that the distribution of individual times to the first cell division varied greatly based on salt concentration. In fact, at high salt concentrations, the distribution pattern appeared distinctly bimodal. However, in earlier work [7], such asymmetric population distributions were interpreted as being Gamma-distributed.

The experiments represented in Fig. 2 were repeated using mid-log phase-associated cells as described in the Experimental section and we saw qualitatively similar results (Fig. 4). The main graph in Fig. 4 represents 987 OD[t] observations with the calculated values of τ plotted as a function of C_{I}. At C_{I}s > ca. 1,000 CFU mL^{-1} the average τ was unimodally-distributed with a maximum spread of ca. 17 to 22 min (159 observations; μ_{τ} ± σ_{τ} = 17.9 ± 0.645 min). Similar to the stationary phase-based cells, we see that as C_{I} was decreased (C_{I} ≤ 200 CFU mL^{-1} or ≤ 54 ± 7.3 CFU/well), a striking increase occurred in the scatter of τ (spread between 12 and 36 min). The frequency of occurrence of all log phase-based τ values (C_{I} < 1,000 CFU mL^{-1}) are displayed in the inset graph of Fig. 4 (α ~ 0.35; μ_{τ1} ± σ_{τ1} = 18.2 ± 0.660 min; β ~ 0.65; μ_{τ2} ± σ_{τ1} = 20.0 ± 2.11 min).

It is important to keep in mind throughout this work that by the time we begin to observe an increase in OD (and therefore measure τ via **Eq. 1**), somewhere between 2 and 20 doublings will have occurred. This fact implies that the values we observe are somehow modulated based upon initial conditions. It should also be noted that low bacterial C_{I}s (i.e., ≤ 5 CFU mL^{-1}) would result in at least some single CFU occurrences per well (i.e., the average probability of observing 1 CFU per well should be about 32.0 ± 6.65%) at which point the first few events of cell division could modulate characteristics of both τ and true microbiological lag time (T). Thus, some of the increase in τ and T scatter we observe at low C_{I} could result from the random selection of isolates with particularly slow growth rates which would otherwise be masked by other isolates in the media with faster rates. However, arguing against such a stochastically-based explanation is the fact that a significant fraction of the scatter in τ (Figs. 2 and 4) occurs between C_{I} = 10-100 CFU mL^{-1} whereupon the probability of observing 1 CFU per well only ranges from 18.1 to ca. 0%. Under these conditions the random selection of one particular τ-component would be overwhelmed by the sheer number of other cells present. At slightly higher concentrations (e.g., 2 or 3 CFUs per well), any well which has 2 or 3 cells with τ values differing more than about 4 or 5 min would be obvious in the ∂OD[t]/∂t curves as additional peaks. Nevertheless, we just don't observe such behavior at these low C_{I}s. What we do observe are relatively uniform, monotypic growth curves (examples in Methods Section) indicative of one component (or, if more than one, the Δτs are small). The fact that we see much greater τ-based scatter at a relatively large threshold C_{I} argues that there is some other controlling factor in determining such binomial-based population growth rates.

In order to determine if the apparent C_{I} effect on τ was only associated with our native *E. coli* strain, we tested two other bacterial strains (*E. coli* O157:H7 and *Citrobacter*). Table 2 summarizes τ frequency distribution parameters (**Eq. 7**, Methods Section) from the experiments represented in Figs. 2 and 4 as well as results concerning mid-log phase *E. coli* O157:H7 and *Citrobacter* in LB, *E. coli* in MM or LB with 75 mM ethyl acetate (EA; solvent for N-acyl homoserine lactones). The stationary or log phase-based generic *E. coli* or *E. coli* O157:H7 growth data in LB gave similar results: for the narrower portion of the bimodal Gaussian distribution, the population mean τ values (μ_{τ1}) varied only 18.0 to 18.5 min (σ_{τ1} 0.401 to 0.678); the broader part of the distribution was also very similar (μ_{τ2} = 19.9 to 20.1 min; σ_{τ2} 2.01 to 2.48). Utilizing MM rather than LB with generic *E. coli* cells from log phase cultures, we saw that the τ distribution on initial cell concentration remained as apparent as the phenomenon in LB (μ_{τ1} ± σ_{τ1} = 51.1 ± 1.75 min; μ_{τ2} ± σ_{τ2} = 56.9 ± 8.32 min), which is consistent with other work (Table 1). The Gram negative bacterium *Citrobacter* (Table 2), which was also grown in LB with cells from log phase cultures, had relatively large doubling times but displayed a clear bimodal distribution in τ at low cell densities (α = 0.6, μ_{τ1} ± σ_{τ1} = 42.5 ± 3.75 min; β = 0.4, μ_{τ2} ± σ_{τ2} = 50.7 ± 6.5 min) similar to previous observations. However, the ethyl acetate set of experiments (LB with 75 mM EA) with *E. coli*, which were performed as a positive control for testing various N-acyl homoserine lactones (AHL; in Gram-negative bacteria AHL is one of two major types of quorum sensing compounds believed to regulate various aspects of bacterial physiology depending upon population size), showed that EA nearly collapsed the bimodal distribution (Fig. 5) to a unimodal form as a result. We observed that α dropped to 0.15 from an LB average of 0.41 (± 0.066), μ_{τ1} shifted upward 1.4 min, and σ_{τ1} broadened by 0.339 min. This result argues for a physiological basis for the increased τ scatter at C_{I} below 100 (stationary phase Fig. 2) to 1,000 (log phase Fig. 4) CFU mL^{-1}. Because of the relatively large effect of solvent alone, the AHL experiments were not performed.

Since there was an obvious dependence of τ on C_{I}, we were interested in determining if the bimodal effect could be reversed by growth in sterile-filtered LB media, which formerly contained the same bacterial isolate (i.e., 'conditioned' media), thus testing to see if an extracellular molecule modulated the bimodal distribution effect (i.e., related to quorum sensing). In one set of experiments (stationary phase inoculum) the LB diluent was made as follows: 37°C LB was inoculated with stationary phase *E. coli* cells and grown several hrs at 37°C (up to ca. 500 CFU mL^{-1}) followed by sterile filtering (2 μm) after centrifugation. These observations are plotted adjacent to control data (Fig. 2) in Fig. 6. A second (log phase cells) experiment was also performed (after harvesting an inoculum for the experiment, the mid-log phase LB medium was centrifuged, sterile-filtered and 20 μL added to each well for the growth experiment), with the results shown in Table 3. Both experiments showed that there was a shift in the low C_{I} bimodal populations (Δμ_{τ} from 1.8 to 1 min) but the bimodal effect was still apparent. The treatments depicted in Fig. 6 also clearly conceptualize the line broadening of the narrow distribution component, the relative decrease in α in the bimodal population, as well as the shift of the two bimodal components towards each other. Thus, some component exists in the media which somewhat modulates the growth process. Lastly, when approximately 2 × 10^{5} sonicated/heat-killed cells mL^{-1} in fresh LB were utilized as the diluent but with the starting innocula taken from a log phase culture, the effect was to induce the narrow component's average τ to shift to that of the broad component (e.g., μ_{τ1} ~ μ_{τ2}, Δ μ~ 0; Fig. 7A, left hand side of plots). Fig. 7B shows τ data plotted as a function of C_{I} and clearly shows the initial concentration effect of τ scatter below 100 CFU mL^{-1}. These results also argue for a physiological basis for the increased τ scatter at relatively low C_{I} (Figs. 2 and 4).

## Conclusion

Working with a native, food-borne *E. coli* isolate grown in either LB or MM, we found that microplate-based doubling times were bimodally distributed at low cell densities using either log or stationary phase cells as an initial inoculum. Qualitatively identical results were obtained for an *E. coli* O157:H7 and *Citrobacter* strain. When sterile-filtered 'conditioned' LB media (formerly contained relatively low concentrations of bacteria or sonicated/heat-killed cells) were employed as a diluent, there were apparent shifts in the two (narrow and broad) populations but the bimodal effect was still evident. However, the bimodal response was almost completely reversed when the growth media contained a small amount of ethyl acetate.

The clear doubling time-cell concentration dependency shown in these results might indicate that bacteria exude a labile biochemical which controls τ, or a need for cell-to-cell physical contact. The latter proposal seems unlikely inasmuch as the probability of random contact would be small at such low cell densities (C_{I} ~ 100-1,000 CFU mL^{-1}). Perhaps this anomalous bimodal distribution of doubling times is related to the recently proposed phenotypic switching [14, 15] which describes programmed variability in certain bacterial populations.

## Methods

### General

*Escherichia coli* (non-pathogenic chicken isolate) [11], *E. coli* O157:H7 (CDC isolate B1409), and *Citrobacter freundii* (non-pathogenic poultry isolate; identification based on 16 S rDNA analysis) [16] were cultured using LB (Difco) or MM (60 mM K_{2}HPO_{4}, 33 mM KH_{2}PO_{4}, 8 mM (NH_{4})_{2}SO_{4}, 2 mM C_{6}H_{5}O_{7}Na_{3} [Na Citrate], 550 μM MgSO_{4}, 14 μM C_{12}H_{18}C_{l2}Na_{4}OS [Thiamine•HCl], 12 mM C_{6}H_{12}O_{6} [glucose], pH 6.8). Liquid cultures were incubated with shaking (200 RPM) at 37°C for ca. 2-4 (for log phase cultures) or 18 hrs (stationary phase cultures) using either LB or MM. All total aerobic plate counts (TAPC) were performed using the 6 × 6 drop plate method [17] with LB followed by incubation at 20-22°C (lab temperature) for 16-18 hours. Using Microsoft Excel's formulaic protocol, the TAPC-based doubling time = 1/LINEST(LOG(TAPC_{1}:TAPC_{n},2),t_{1}:t_{n}) where the values TAPC_{1} through TAPC_{n} are log-linear with respect to associated growth times t_{1} to t_{n}; n was typically 6-8 points. All TAPC studies were performed using highly diluted stationary phase cells (initial colony forming unit [CFU] concentration or C_{I} ≥ 10^{3} CFU mL^{-1}) in either LB or MM.

### Steady State Oxygen

O_{2} levels ([O_{2}], units of μM) were measured using a Clark-type oxygen electrode (Model 5300, Yellow Spring Instruments) connected to a Gilson water-jacketed chamber (1.42 mL; circulating water bath attached, 37°C) containing a magnetic stirring bar. Air-saturated 37°C water was used for calibration. To determine steady-state [O_{2}] in shaking/bubbled cultures, samples were withdrawn with a syringe from bacterial culture flasks at various time points during mid-to late-log phase growth, and the oxygen consumption (e.g., [O_{2}] dropping with time) determined without vortexing. The time lapse between sample withdrawal and the first [O_{2}] data point was recorded and used to back-calculate the [O_{2}] at the time of sampling. These same samples were then vortexed ca. 15 sec and [O_{2}] measured again as a function of time. The rate of O_{2} consumption was calculated from the slope of cell density-normalized [O_{2}] (TAPC plating was performed simultaneously on LB) as a function of time (apparent K_{m} ~ 15 ± 6 μM) [18].

### 96-well Microplate Protocol

In order to avoid water condensation which might interfere with absorbance readings, the interior surface of microplate covers were rinsed with a solution of 0.05% Triton X-100 in 20% ethanol [12] and dried in a microbiological hood under UV light. About 270 μL of each bacterial cell concentration was pipetted into every well. Each initial concentration (C_{I}) is equal to C_{0} Φ_{I} where C_{0} is the cell density from liquid culture (either log or stationary phase). When C_{0} ≤ 10^{8} CFU mL^{-1}, the cells were sampled from an early-to mid-log phase culture. When C_{0} ≥ 10^{9} CFU mL^{-1}, the cells were sampled from a stationary phase culture. Typically, each 96-well microplate contained 2 replicates each of the 8 least dilute samples (Φ = 3×10^{-3} to 5×10^{-6}; 16 wells), 4 replicates of the next 4 highest dilutions (Φ = 2×10^{-6} to 5×10^{-7}; 16 wells), 8 replicates each of the following 2 dilutions (Φ = 2×10^{-7} or 1×10^{-7}; 16 wells), and, lastly, 24 replicates of the 2 most dilute samples (Φ = 6×10^{-8} or 3×10^{-8}; 48 wells). The 96-well plate was then covered with the Triton-treated top, placed in a temperature-equilibrated Perkin-Elmer HTS 7000+ 96-well microplate reader, and monitored for optical density (OD) under the following conditions: λ = 590 nm; the time between points (Δt) = 10-25 min; total points = 50-110; temperature = 37°C; 5 sec of moderate shaking before each reading (see Results section). These dilutions, listed above, produced at least some negative (no growth) readings mainly associated with the 4 most dilute sets of wells. This lack of growth in wells associated with these dilutions is evidence for single CFU-based growth occurrences at these low C_{I}. Thus, these low C_{I} have been diluted to such a degree that at least an occasional random sampling of 270 μL should contain no cells at all. Generally speaking, the most probable number (single dilution MPN) calculation for these dilutions agreed with the plate count estimate. The variability of growth parameters at such low concentrations (~ 1 CFU/well) has generated much recent interest [4, 6–8].

### Calculations

After completion of any OD with time growth experiment, a tab-delimited text file was generated and data pasted into a Microsoft Excel spreadsheet formatted to display the data arrays as individual well ODs associated with each time. Typical OD growth curves are presented in Fig. 8 which have been curve-fitted (non-linear regression) to the Boltzmann equation (**Eq. 1**), a well-known sigmoidal function used in various physiological studies [19]

While Eq. 1 is an empirical equation, it does rely on a first order rate constant (k) therefore the doubling time can be extracted as τ = k^{-1} Ln [2]. All curve-fitting was performed using a Gauss-Newton algorithm on an Excel spreadsheet [20]. In **Eq. 1**, OD_{I} is the estimated initial optical density (0.05-0.1), OD_{F} is the calculated final OD (0.5-0.7), k is the first-order rate constant, and t_{m} is the time to OD_{F} ÷ 2. The Boltzmann relationship appears to be generally useful with optically-based growth results since excellent fits were achieved (21°C growth in LB, τ = 3.26 ± 0.0292 hrs) when **Eq. 1** was utilized to fit previously published [21] bacterial growth data from a microchemostat.

As demonstrated previously [12], t_{m} can be used (for high C_{I}) as a method for estimating cell density. The inset plot in Fig. 8 shows both OD and first derivative (ΔOD/Δt) versus time data sets that were typically observed when growing our native *E. coli* isolate in MM. In order to achieve the best fit in the region which provides the most information (i.e., the exponential increase in OD), we have truncated these data and used only 2-10 points beyond the apparent t_{m} to fit to **Eq. 1**. Such data abbreviation had only minor effects on the growth parameters: e.g., if the OD[t] data points in the main plot of Fig. 8 were truncated to only 3 points past the calculated t_{m}, τ would change only from ~ 19.2 to 19.8 min and t_{m} only by 0.7 min. All values of τ and t_{m} reported herein are derived from such curve-fitting. Of course, t_{m} can also be estimated from the x-axis value where the center of symmetry in ΔOD/Δt occurs (Fig. 8). We have tested two other microplate readers (Bio-Tek EL 312e and Tecan Safire II) in order to determine the variability in τ (from OD[t] data; C_{I} > 1000 CFU mL^{-1}) due to the devices themselves. The Perkin-Elmer instrument consistently gave the lowest τ values (τ = 18 ± 0.99 min) followed by the Bio-Tek (τ = 19 ± 1.0 min) and Tecan (τ = 21 ± 1.2 min); {Error Mean Square ÷ n} ^{1/2}. = 0.42. It seems likely that the observed plate reader-associated differences in τ are due to instrument-based disparities in temperature.

During the log phase of growth [3], the rate of change in bacterial concentration with respect to time can be represented by the simple differential equation

in this relation, k is a first order rate constant, t is the growth time, and C is the bacterial concentration. Upon rearrangement, integration between initial (C_{I}) and final (C_{F}) values of C, expressing k in terms of a doubling or generation time (τ = k^{-1} Ln(2)) and solving for C_{F} we see that

where T is a time translation constant utilized to correct for the observed lag in cell growth. In our usage we assume that C_{F} is the cell density at which the relationship between OD and C becomes non-linear. For our wild-type *E. coli* isolate [11] C_{F} was typically about 5×10^{8} CFU mL^{-1}. Expressing **Eq. 3** in terms of the time it takes to reach C_{F} (OD ~ 0.6) we see that

Since it is facile to approximate the value of t when C = C_{F} ÷ 2 and t = t_{m} (Fig. 8), we have chosen to express **Eq. 4** in terms of t_{m}; making this alteration, substituting C_{0}Φ_{I} for C_{I} and rearrangement gives

In **Eq. 5** Φ_{I} is the dilution factor (e.g., for a C_{I} resulting from two 1:10 dilutions Φ_{I} = 0.1 × 0.1 = 10^{-2}) and C_{0} is the starting cell density (e.g., from either a mid-log or stationary phase suspension of cells) from which all dilutions are made. In this work C_{0} was either about 10^{8} (cells sampled from a mid-log phase culture; media-corrected OD_{590-600} < 0.1) or 10^{9} (stationary phase) CFU mL^{-1}. **Eq. 5** implies that τ can be determined by calculating the slope from a plot of t_{m} versus Log_{2} [Φ_{I}] (Excel τ = ABS (LINEST(t_{m,1} : t_{m,n}, LOG(Φ_{I,1} : Φ_{I,n,2})))). Fig. 9 displays both linear and semi-log plots of typical t_{m} data plotted as a function of Φ_{I}. Of course, identical results to the above are obtained if C_{I} replaces C_{0}Φ_{I} (i.e., **Eq. 5** with C_{0} deleted and C_{I} substituted for Φ_{I})

and a plot of t_{m} with Log_{2} [CI] is linear (Excel τ = ABS (LINEST(t_{m,1}:t_{m,n}, LOG(C_{I,1}:C_{I,n},2)))) with a slope equal to -τ and an intercept of (T + Log_{2} [C_{F}/2]).

**Eq. 6** implies that the time in lag phase (T) can be obtained knowing τ, C_{F}, and the intercept from a plot of t_{m} as a function of Log_{2} [C_{I}]. When numerous values of t_{m} are plotted against C_{I} (semi-log plot shown in Fig. 3) by diluting either log or stationary phase cells in LB one sees a significant perturbation in T (offsets in the intercept) of the semi-log plots (10^{2} < C_{I} < 10^{7} CFU mL^{-1} region only). T calculations (Eq. 6) from the growth of stationary phase-diluted cells (T = 41 ± 8.4 min; average of 10 experiments; C_{I} > 10^{2} CFU mL^{-1}) indicate that T was similar to lag times calculated from TAPC experiments (63 ± 9 min; average of 7 experiments). However, T values calculated in a similar fashion from log phase-diluted cells produced near-zero values (T = -11 ± 15 min; average of 8 experiments; C_{I} > 100 CFU mL^{-1}). Thus, the total offset between log and stationary phase-derived cells shown in Fig. 3 was about 52 min and implies that stationary phase cells require about an hour to revert to log-phase. However, because of the variability in the intercept and C_{F}, we believe that the value of T using **Eq. 6** has only a relative meaning. In other words, **Eqs. 5** & **6** show that variability in t_{m} can be due to either variability in T, τ or both.

In order to generate the frequency of occurrence of τ values (obtained using **Eq. 1**), we first created integers from the individual τ values, counted the number of occurrences of each τ then divided this by the total number counted. Thereafter a Gaussian or normal distribution function was used to curve-fit [20] frequency of occurrence of τ data to the individually-observed τ integers. The bimodal form consisted of the sum of two Gaussians (Eq. 7) whereupon α + β = 1

In **Eq. 7**, α is the fraction of the population associated with mean μ_{τ1} and standard deviation σ_{τ1}; a second Gaussian is characterized by β (= 1 - α), μ_{τ2}, and σ_{τ2}.

Regarding other statistical methods used in this work: analysis of variance tables were generated using Microsoft Excel and standard statistical formulae for a randomized complete block design. Values for F were taken from a college-level statistics table of F-values.

## References

- 1.
Oscar T: Validation of Lag Time and Growth Rate Models for

*Salmonella Typhimurium*: Acceptable Prediction Zone Method. J Food Sci. 2005, 70: M129-M137. 10.1111/j.1365-2621.2005.tb07103.x. - 2.
Kutalik Z, Razaz M, Baranyi J: Connection between stochastic and deterministic modeling of microbial growth. J Theor Biol. 2005, 232: 285-299. 10.1016/j.jtbi.2004.08.013.

- 3.
Lopez S, Prieto M, Dijkstra J, Dhanoa M, France J: Statistical Evaluation of Mathematical Models for Microbial Growth. Int J Food Microbiol. 2004, 96: 289-300. 10.1016/j.ijfoodmicro.2004.03.026.

- 4.
Elfwing A, LeMarc Y, Baranyi J, Ballagi A: Observing growth and division of large numbers of individual bacteria by image analysis. Appl Environ Microbiol. 2004, 70: 675-678. 10.1128/AEM.70.2.675-678.2004.

- 5.
Guillier L, Pardon P, Augustin J-C: Influence of Stress on Individual Lag Time Distributions of

*Listeria monocytogenes*. Appl Environ Microbiol. 2005, 71: 2940-2948. 10.1128/AEM.71.6.2940-2948.2005. - 6.
Guillier L, Pardon P, Augustin J-C: Automated image analysis of bacterial colony growth as a tool to study individual lag time distributions of immobilized cells. J Microbiol Methods. 2006, 65: 324-334. 10.1016/j.mimet.2005.08.007.

- 7.
Metris A, George S, Peck M, Baranyi J: Distribution of turbidity detection times produced by single cell-generated bacterial populations. J Microbiol Methods. 2003, 55: 821-827. 10.1016/j.mimet.2003.08.006.

- 8.
Niven G, Fuks T, Morton J, Rua S, Mackey B: A novel method for measuring lag times in division of individual bacterial cells using image analysis. J Microbiol Methods. 2006, 65: 311-317. 10.1016/j.mimet.2005.08.006.

- 9.
Irwin P, Damert W, Brewster J, Gehring A, Tu S-I: Immuno-magnetic bead mass transport and capture efficiency at low target cell densities in phosphate-buffered saline. J Rapid Methods Autom Microbiol. 2002, 10: 129-147. 10.1111/j.1745-4581.2002.tb00018.x.

- 10.
Irwin P, Damert W: Immuno-magnetic bead mass transport and capture efficiency at high target cell densities in phosphate-buffered saline. J Rapid Methods Autom Microbiol. 2004, 11: 265-284. 10.1111/j.1745-4581.2003.tb00267.x.

- 11.
Irwin P, Gehring A, Tu S-I, Chen C-Y: Blocking nonspecific adsorption of native foodborne microorganisms by immunomagnetic beads withι-carrageenan. Carbohydr Res. 2004, 339: 613-621. 10.1016/j.carres.2003.10.033.

- 12.
Brewster J: A simple micro-growth assay for enumerating bacteria. J Microbiol Methods. 2002, 53: 77-86. 10.1016/S0167-7012(02)00226-9.

- 13.
Irwin P, Gehring A, Tu S-I, Brewster J, Fanelli F, Ehrenfeld E: Minimum Detectable Level of Salmonellae Using a Binomial-Based Bacterial Ice Nucleation Detection Assay. J AOAC Int. 2000, 83: 1087-1095.

- 14.
Balaban N, Merrin J, Chait R, Kowalik L, Leibler S: Bacterial Persistence as a Phenotypic Switch. Science. 2004, 305: 1622-1625. 10.1126/science.1099390.

- 15.
Kussell E, Leibler S: Phenotypic Diversity, Population Growth, and Information in Fluctuating Environments. Science. 2005, 309: 2075-2078. 10.1126/science.1114383.

- 16.
Irwin P, Nguyen L-H, Chen C-Y, Paoli G: Binding of nontarget microorganisms from food washes to anti-Salmonella and anti-E. coli O157 immunomagnetic beads: most probable composition of background Eubacteria. Anal Bioanal Chem. 2008, 391: 525-536. 10.1007/s00216-008-1959-2.

- 17.
Chen C-Y, Nace G, Irwin P: A 6× 6 drop plate method for simultaneous colony counting and MPN enumeration of

*Campylobacter jejuni*,*Listeria monocytogenes*, and*Escherichia coli*. J Microbiol Methods. 2003, 55: 475-479. 10.1016/S0167-7012(03)00194-5. - 18.
Irwin P, Brouillette J, Germann M, Hicks K, Kurantz M, Damert W: Calculation of immobilized enzyme reaction progress curves from nested ordered-sequential rate expressions. Enzyme Microb Technol. 1999, 24: 675-686. 10.1016/S0141-0229(99)00004-6.

- 19.
Valiunas V, Manthey D, Vogel R, Willecke K, Weingart R: Biophysical properties of mouse connexin30 gap junction channels studied in transfected human HeLa cells. J Physiol. 1999, 519: 631-644. 10.1111/j.1469-7793.1999.0631n.x.

- 20.
Irwin P, Damert W, Doner L: Curve fitting in nuclear magnetic resonance: illustrative examples using a spreadsheet and microcomputer. Concepts Magn Reson. 1994, 6: 57-67. 10.1002/cmr.1820060105.

- 21.
Balagadde F, You L, Hansen C, Arnold F, Quake S: Long-Term Monitoring of Bacteria Undergoing Programmed Population Control in a Microchemostat. Science. 2005, 309: 137-140. 10.1126/science.1109173.

## Acknowledgements

All funding was from ARS base funds associated with Current Research Information System (CRIS) Project Number 1935-42000-058-00 D (Integrated Biosensor-Based Processes for Multipathogenic Analyte Detection).

## Author information

## Additional information

### Authors' contributions

PI designed all of the experiments, performed all calculations and statistical analyses, participated in running most of the experiments and drafting the manuscript. LN carried out all the TAPC and O_{2} electrode experiments and participated in drafting the manuscript. GP and CC assisted in the experiments using conditioned media, MM, and LB with disrupted cells and participated in O_{2} electrode experiments as well as drafting the manuscript. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Irwin, P.L., Nguyen, L.T., Paoli, G.C. *et al.* Evidence for a bimodal distribution of *Escherichia coli* doubling times below a threshold initial cell concentration.
*BMC Microbiol* **10, **207 (2010) doi:10.1186/1471-2180-10-207

#### Received

#### Accepted

#### Published

#### DOI

### Keywords

- Minimal Medium
- Stationary Phase Cell
- Initial Cell Concentration
- Coli Isolate
- Bimodal Effect