Evidence for a bimodal distribution of Escherichia coli doubling times below a threshold initial cell concentration.  
Jump to Full Text  
MedLine Citation:

PMID: 20678197 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

BACKGROUND: In the process of developing a microplatebased growth assay, we discovered that our test organism, a native E. coli isolate, displayed very uniform doubling times (tau) only up to a certain threshold cell density. Below this cell concentration (<or= 100 1,000 CFU mL1 ; <or= 27270 CFU well1) we observed an obvious increase in the tau scatter. RESULTS: Working with a foodborne E. coli isolate we found that tau values derived from two different microtiter platereaderbased techniques (i.e., optical density with growth time {=OD[t]} fit to the sigmoidal Boltzmann equation or time to calculated 1/2maximal OD {=tm} as a function of initial cell density {=tm[CI]}) were in excellent agreement with the same parameter acquired from total aerobic plate counting. Thus, using either LuriaBertani (LB) or defined (MM) media at 37 degrees C, tau ranged between 1718 (LB) or 5154 (MM) min. Making use of such OD[t] data we collected many observations of tau as a function of manifold initial or starting cell concentrations (CI). We noticed that tau appeared to be distributed in two populations (bimodal) at low CI. When CI <or=100 CFU mL1 (stationary phase cells in LB), we found that about 48% of the observed tau values were normally distributed around a mean (mutau1) of 18 +/ 0.68 min (+/ sigmatau1) and 52% with mutau2 = 20 +/ 2.5 min (n = 479). However, at higher starting cell densities (CI>100 CFU mL1), the tau values were distributed unimodally (mutau = 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 tau distribution results were also observed using E. coli cells diluted from midlog phase cultures. Similar results were also obtained when using either an E. coli O157:H7 or a Citrobacter strain. When sterilefiltered LB supernatants, which formerly contained relatively low concentrations of bacteria(1,00010,000 CFU mL1), 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. 
Authors:

Peter L Irwin; LyHuong T Nguyen; George C Paoli; ChinYi Chen 
Related Documents
:

40957  Sucrose transport by the escherichia coli lactose carrier. 18944927  Production of a mating inhibitor by ustilago hordei. 9570127  Identification of the mgc1 mutation which affects matingpheromoneinduced morphogenesi... 348097  Susceptibility of escherichia coli to bactericidal action of lactoperoxidase, peroxide,... 22422377  Efficient intracellular delivery of camptothecin by silica/titania hollow nanoparticles. 7262207  Scanning electron microscope observations on surface changes occurring in murine granul... 
Publication Detail:

Type: Journal Article; Research Support, NonU.S. Gov't Date: 20100802 
Journal Detail:

Title: BMC microbiology Volume: 10 ISSN: 14712180 ISO Abbreviation: BMC Microbiol. Publication Date: 2010 
Date Detail:

Created Date: 20100818 Completed Date: 20101116 Revised Date: 20130529 
Medline Journal Info:

Nlm Unique ID: 100966981 Medline TA: BMC Microbiol Country: England 
Other Details:

Languages: eng Pagination: 207 Citation Subset: IM 
Affiliation:

Molecular Characterization of Foodborne Pathogens Research Unit, Agricultural Research Service, US Department of Agriculture, Eastern Regional Research Center, Wyndmoor, PA 19038, USA. peter.irwin@ars.usda.gov 
Export Citation:

APA/MLA Format Download EndNote Download BibTex 
MeSH Terms  
Descriptor/Qualifier:

Animals Chickens / microbiology Escherichia coli / chemistry, growth & development* Kinetics 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmta): BMC Microbiol ISSN: 14712180 Publisher: BioMed Central 
Article Information Copyright ©2010 Irwin et al; licensee BioMed Central Ltd. openaccess: Received Day: 28 Month: 1 Year: 2010 Accepted Day: 2 Month: 8 Year: 2010 collection publication date: Year: 2010 Electronic publication date: Day: 2 Month: 8 Year: 2010 Volume: 10First Page: 207 Last Page: 207 Publisher Id: 1471218010207 PubMed Id: 20678197 DOI: 10.1186/1471218010207 
Evidence for a bimodal distribution of Escherichia coli doubling times below a threshold initial cell concentration  
Peter L Irwin1  Email: peter.irwin@ars.usda.gov 
LyHuong T Nguyen1  Email: Ly.Nguyen@ars.usda.gov 
George C Paoli1  Email: george.paoli@ars.usda.gov 
ChinYi Chen1  Email: ChinYi.Chen@ars.usda.gov 
1Molecular Characterization of Foodborne Pathogens Research Unit, Eastern Regional Research Center, Agricultural Research Service, U. S. Department of Agriculture, 600 East Mermaid Lane, Wyndmoor, PA 19038 USA 
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, probabilitybased 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 coworkers studied the effects of various stress factors (temperature, starvation, salt concentration, etc.) on individual cellbased 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 populationbased information on specific growth rate is lacking.
The findings presented herein developed from work associated with the attachment of various Gramnegative bacteria to antiSalmonella and antiE. coli O157 immunomagnetic beads or IMBs [^{9}^{}^{11}]. For these IMB investigations microplate (ODbased) MPN methods were utilized because of the low limits of bacterial detection [^{12},^{13}] necessary to characterize the nonspecific attachment of background food organisms to various capture surfaces. Because of large interbacterial 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/2maximal 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 nonspecifically 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.
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 LuriaBertani (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 withinmedium 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 microplatebased 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 {46} × 10^{18 }moles O_{2 }sec^{1 }CFU^{1 }(Fig. 1B). Even the baffled flask culture showed a drop in [O_{2}] of 4057%. 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 areatovolume 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.
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 Gaussiandistributed (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 concentrationassociated 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. 56: all t_{m }& Tbased 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 timerelated 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 coworkers [^{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 Gammadistributed.
The experiments represented in Fig. 2 were repeated using midlog phaseassociated 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 unimodallydistributed with a maximum spread of ca. 17 to 22 min (159 observations; μ_{τ }± σ_{τ }= 17.9 ± 0.645 min). Similar to the stationary phasebased 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 phasebased τ 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 stochasticallybased explanation is the fact that a significant fraction of the scatter in τ (Figs. 2 and 4) occurs between C_{I }= 10100 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 binomialbased 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 midlog phase E. coli O157:H7 and Citrobacter in LB, E. coli in MM or LB with 75 mM ethyl acetate (EA; solvent for Nacyl homoserine lactones). The stationary or log phasebased 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 Nacyl homoserine lactones (AHL; in Gramnegative 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 sterilefiltered 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 midlog phase LB medium was centrifuged, sterilefiltered 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/heatkilled 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).
Working with a native, foodborne E. coli isolate grown in either LB or MM, we found that microplatebased 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 sterilefiltered 'conditioned' LB media (formerly contained relatively low concentrations of bacteria or sonicated/heatkilled 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 timecell concentration dependency shown in these results might indicate that bacteria exude a labile biochemical which controls τ, or a need for celltocell physical contact. The latter proposal seems unlikely inasmuch as the probability of random contact would be small at such low cell densities (C_{I }~ 1001,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.
Escherichia coli (nonpathogenic chicken isolate) [^{11}], E. coli O157:H7 (CDC isolate B1409), and Citrobacter freundii (nonpathogenic 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. 24 (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 2022°C (lab temperature) for 1618 hours. Using Microsoft Excel's formulaic protocol, the TAPCbased doubling time = 1/LINEST(LOG(TAPC_{1}:TAPC_{n},2),t_{1}:t_{n}) where the values TAPC_{1 }through TAPC_{n }are loglinear with respect to associated growth times t_{1 }to t_{n}; n was typically 68 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.
O_{2 }levels ([O_{2}], units of μM) were measured using a Clarktype oxygen electrode (Model 5300, Yellow Spring Instruments) connected to a Gilson waterjacketed chamber (1.42 mL; circulating water bath attached, 37°C) containing a magnetic stirring bar. Airsaturated 37°C water was used for calibration. To determine steadystate [O_{2}] in shaking/bubbled cultures, samples were withdrawn with a syringe from bacterial culture flasks at various time points during midto latelog 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 backcalculate 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 densitynormalized [O_{2}] (TAPC plating was performed simultaneously on LB) as a function of time (apparent K_{m }~ 15 ± 6 μM) [^{18}].
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 X100 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 earlyto midlog phase culture. When C_{0 }≥ 10^{9 }CFU mL^{1}, the cells were sampled from a stationary phase culture. Typically, each 96well 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 96well plate was then covered with the Tritontreated top, placed in a temperatureequilibrated PerkinElmer HTS 7000+ 96well microplate reader, and monitored for optical density (OD) under the following conditions: λ = 590 nm; the time between points (Δt) = 1025 min; total points = 50110; 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 CFUbased 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}].
After completion of any OD with time growth experiment, a tabdelimited 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 curvefitted (nonlinear regression) to the Boltzmann equation (Eq. 1), a wellknown sigmoidal function used in various physiological studies [^{19}]
(1)
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 curvefitting was performed using a GaussNewton algorithm on an Excel spreadsheet [^{20}]. In Eq. 1, OD_{I }is the estimated initial optical density (0.050.1), OD_{F }is the calculated final OD (0.50.7), k is the firstorder rate constant, and t_{m }is the time to OD_{F }÷ 2. The Boltzmann relationship appears to be generally useful with opticallybased 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 210 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 curvefitting. Of course, t_{m }can also be estimated from the xaxis value where the center of symmetry in ΔOD/Δt occurs (Fig. 8). We have tested two other microplate readers (BioTek 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 PerkinElmer instrument consistently gave the lowest τ values (τ = 18 ± 0.99 min) followed by the BioTek (τ = 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 readerassociated differences in τ are due to instrumentbased 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
(2)
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
(3)
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 nonlinear. For our wildtype 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
(4)
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 midlog 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 midlog phase culture; mediacorrected OD_{590600 }< 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 semilog 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 }(semilog 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 semilog plots (10^{2 }< C_{I }< 10^{7 }CFU mL^{1 }region only). T calculations (Eq. 6) from the growth of stationary phasediluted 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 phasediluted cells produced nearzero values (T = 11 ± 15 min; average of 8 experiments; C_{I }> 100 CFU mL^{1}). Thus, the total offset between log and stationary phasederived cells shown in Fig. 3 was about 52 min and implies that stationary phase cells require about an hour to revert to logphase. 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 curvefit [^{20}] frequency of occurrence of τ data to the individuallyobserved τ integers. The bimodal form consisted of the sum of two Gaussians (Eq. 7) whereupon α + β = 1
(7)
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 collegelevel statistics table of Fvalues.
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.
All funding was from ARS base funds associated with Current Research Information System (CRIS) Project Number 19354200005800 D (Integrated BiosensorBased Processes for Multipathogenic Analyte Detection).
References
Oscar T,Validation of Lag Time and Growth Rate Models for Salmonella Typhimurium: Acceptable Prediction Zone MethodJ Food SciYear: 200570M129M137  
Kutalik Z,Razaz M,Baranyi J,Connection between stochastic and deterministic modeling of microbial growthJ Theor BiolYear: 200523228529910.1016/j.jtbi.2004.08.01315530497  
Lopez S,Prieto M,Dijkstra J,Dhanoa M,France J,Statistical Evaluation of Mathematical Models for Microbial GrowthInt J Food MicrobiolYear: 20049628930010.1016/j.ijfoodmicro.2004.03.02615454319  
Elfwing A,LeMarc Y,Baranyi J,Ballagi A,Observing growth and division of large numbers of individual bacteria by image analysisAppl Environ MicrobiolYear: 20047067567810.1128/AEM.70.2.675678.200414766541  
Guillier L,Pardon P,Augustin JC,Influence of Stress on Individual Lag Time Distributions of Listeria monocytogenesAppl Environ MicrobiolYear: 2005712940294810.1128/AEM.71.6.29402948.200515932988  
Guillier L,Pardon P,Augustin JC,Automated image analysis of bacterial colony growth as a tool to study individual lag time distributions of immobilized cellsJ Microbiol MethodsYear: 20066532433410.1016/j.mimet.2005.08.00716185781  
Metris A,George S,Peck M,Baranyi J,Distribution of turbidity detection times produced by single cellgenerated bacterial populationsJ Microbiol MethodsYear: 20035582182710.1016/j.mimet.2003.08.00614607428  
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 analysisJ Microbiol MethodsYear: 20066531131710.1016/j.mimet.2005.08.00616169621  
Irwin P,Damert W,Brewster J,Gehring A,Tu SI,Immunomagnetic bead mass transport and capture efficiency at low target cell densities in phosphatebuffered salineJ Rapid Methods Autom MicrobiolYear: 20021012914710.1111/j.17454581.2002.tb00018.x  
Irwin P,Damert W,Immunomagnetic bead mass transport and capture efficiency at high target cell densities in phosphatebuffered salineJ Rapid Methods Autom MicrobiolYear: 20041126528410.1111/j.17454581.2003.tb00267.x  
Irwin P,Gehring A,Tu SI,Chen CY,Blocking nonspecific adsorption of native foodborne microorganisms by immunomagnetic beads withιcarrageenanCarbohydr ResYear: 200433961362110.1016/j.carres.2003.10.03315013398  
Brewster J,A simple microgrowth assay for enumerating bacteriaJ Microbiol MethodsYear: 200253778610.1016/S01677012(02)002269  
Irwin P,Gehring A,Tu SI,Brewster J,Fanelli F,Ehrenfeld E,Minimum Detectable Level of Salmonellae Using a BinomialBased Bacterial Ice Nucleation Detection AssayJ AOAC IntYear: 2000831087109511048849  
Balaban N,Merrin J,Chait R,Kowalik L,Leibler S,Bacterial Persistence as a Phenotypic SwitchScienceYear: 20043051622162510.1126/science.109939015308767  
Kussell E,Leibler S,Phenotypic Diversity, Population Growth, and Information in Fluctuating EnvironmentsScienceYear: 20053092075207810.1126/science.111438316123265  
Irwin P,Nguyen LH,Chen CY,Paoli G,Binding of nontarget microorganisms from food washes to antiSalmonella and antiE. coli O157 immunomagnetic beads: most probable composition of background EubacteriaAnal Bioanal ChemYear: 200839152553610.1007/s0021600819592  
Chen CY,Nace G,Irwin P,A 6× 6 drop plate method for simultaneous colony counting and MPN enumeration of Campylobacter jejuni, Listeria monocytogenes, and Escherichia coliJ Microbiol MethodsYear: 20035547547910.1016/S01677012(03)00194514529971  
Irwin P,Brouillette J,Germann M,Hicks K,Kurantz M,Damert W,Calculation of immobilized enzyme reaction progress curves from nested orderedsequential rate expressionsEnzyme Microb TechnolYear: 19992467568610.1016/S01410229(99)000046  
Valiunas V,Manthey D,Vogel R,Willecke K,Weingart R,Biophysical properties of mouse connexin30 gap junction channels studied in transfected human HeLa cellsJ PhysiolYear: 199951963164410.1111/j.14697793.1999.0631n.x10457079  
Irwin P,Damert W,Doner L,Curve fitting in nuclear magnetic resonance: illustrative examples using a spreadsheet and microcomputerConcepts Magn ResonYear: 19946576710.1002/cmr.1820060105  
Balagadde F,You L,Hansen C,Arnold F,Quake S,LongTerm Monitoring of Bacteria Undergoing Programmed Population Control in a MicrochemostatScienceYear: 200530913714010.1126/science.110917315994559 
Figures
Tables
Comparison of doubling time (= τ) observations based on total aerobic plate counting with growth time (= TAPC[t]), time to 1/2maximal OD with growth time (= t_{m}[t]), and optical density at 590 nm with time (= OD[t]) in either LB broth or MM at 37°C.
Method  τ (min)  LB  



Exp. 1  2  3  average  F_{2,4}  






TAPC[t]  18.6  17.3  18.1  18.0  3.43 
t_{m}[Φ_{i}]  17.1  17.4  16.8  17.1  P >0.1 
OD[t]  17.9  17.9  17.7  17.8  


Method  τ (min) MM  


Exp. 1  2  3  average  F2,4  






TAPC[t]  52.7  50.1  51.9  51.6  0.886 
t_{m}[Φ_{i}]  50.8  59.9  52.1  54.3  P >>0.1 
OD[t]  50.1  53.8  49.4  51.1 
Comparison of doubling time distribution parameters (Eq. 1) for E. coli, E. coli O157:H7, and Citrobacter in LB, LB + ethyl acetate (EA, 75 mM), or MM at 37°C; S = Stationary phase, L = Log Phase.
C_{I }≤ 100 CFU mL^{1 }  C_{I }≥ 1000 CFU mL^{1}  



Organism (phase)  Medium LB  α  μ_{τ 1 }± σ_{τ1}  β  μ_{τ2 }± σ_{τ2}  Δμ_{τ}  μ_{τ }± σ_{τ} 








E. coli (S)  LB  0.48  18.0 ± 0.678  0.52  19.9 ± 2.48  1.87  17.6 ± 0.708 
E. coli (L)  LB  0.35  18.2 ± 0.660  0.65  20.0 ± 2.11  1.79  17.9 ± 0.645 
E. coli (L)  LB+EA  0.15  19.6 ± 0.999  0.85  21.7 ± 2.25  2.13  21.4 ± 2.06 
E. coli (L)  MM  0.30  51.1 ± 1.75  0.70  56.9 ± 8.32  5.77  52.0 ± 2.09 
E. coli O157:H7 (L)  LB  0.40  18.5 ± 0.401  0.60  20.1 ± 2.01  1.60  18.1 ± 0.438 
Citrobacter (L)  LB  0.6  42.5 ± 3.75  0.40  50.7 ± 6.50  8.24  42.4 ± 3.72 
Comparison of doubling time distribution parameters (Eq. 1) for E. coli in LB, or in LB with sonicated and heatkilled cells at 37°C; S = stationary phase, L = Log phase.
C_{I }≤ 100 CFU mL^{1}  C_{I }≥ 1000 CFU mL^{1}  



Organism (phase)  α  μ_{τ1 }± σ_{τ1}  β  μ_{τ2 }± σ_{τ2}  Δμ_{τ}  μ_{τ }± σ_{τ} 







Control LB (S)  0.48  18.0 ± 0.678  0.52  19.9 ± 2.48  1.87  17.6 ± 0.708 
Conditioned LB (S)  0.21  17.8 ± 0.553  0.79  18.8 ± 1.99  1.03  17.5 ± 1.06 
Control LB (L)  0.35  18.2 ± 0.660  0.65  20.0 ± 2.11  1.79  17.9 ± 0.645 
Conditioned LB (L)  0.31  19.1 ± 0.627  0.69  20.1 ± 2.10  0.994  18.9 ± 0.700 
Sonicated, Heatkilled Cells in LB (L)  0.54  21.0 ± 0.690  0.46  21.3 ± 2.58  0.300  21.1 ± 0.646 
Article Categories:

Previous Document: Association of invasionpromoting tenascinC additional domains with breast cancers in young women.
Next Document: Reliability and validity of the Spanish version of the Child Health and Illness Profile (CHIP) Child...