Previous Article | Next Article ![]()
Applied and Environmental Microbiology, June 2008, p. 3757-3763, Vol. 74, No. 12
0099-2240/08/$08.00+0 doi:10.1128/AEM.02551-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

School of Food Biosciences, The University of Reading, P.O. Box 226, Whiteknights, Reading RG6 6AP, United Kingdom,1 Dstl Porton Down, Salisbury, Wiltshire SP4 0JQ, United Kingdom2
Received 12 November 2007/ Accepted 11 April 2008
|
|
|---|
|
|
|---|
Lag time distributions can be determined directly by observing single cells under the microscope (5, 11, 22, 25, 32) or indirectly by measuring the time to produce a detectable optical density change in replicate cultures inoculated with approximately one cell using an automated growth analyzer such as the Bioscreen apparatus (7, 12, 19, 26, 29, 31, 32). An alternative indirect method consists of measuring the time to reach a certain colony size from single cells inoculated on an agar plate (8). Indirect methods are more convenient than microscopic methods and are essential when studying severely stressed populations containing only a few viable survivors. Distributions of detection times are a very good measure of the distribution of lag times (20), but they do not always give accurate absolute values for single-cell lag times because extrapolation from population to single-cell levels amplifies the effect of small measurement uncertainties. For this reason, microscopic methods based on direct observation are often regarded as the "gold standard" for studying single-cell behavior. Growth and division of single cells have been monitored in agar slide cultures viewed by light microscopy (10, 23, 24) and in a flow chamber in which daughter cells are flushed away from adherent mother cells after cell division (5). These methods are complementary, but each has certain limitations. The flow chamber method of Elfwing et al. (5) allows large numbers of cells to be monitored and gives a clear-cut point of cell division, but, under some conditions, problems can arise with spontaneous detachment of cells during the period of observation (19). The slide culture method is simple and can be used with a wide range of organisms but is subject to uncertainty when used to measure lag times because it is difficult to decide consistently and objectively exactly when cells have divided. However, a method was recently described that enables the time to first cell division to be determined objectively using image analysis (22). The time of first cell division is calculated from the "box area ratio" (BAR), which is the area of the smallest rectangle that can be drawn around a cell divided by the area of the cell itself. The box area ratio was found to increase abruptly during cell growth at the point of cell division as judged by eye (22). This allows the time of division to be measured objectively and opens up opportunities for systematic studies of the effects of previous history and current growth conditions on lag time distributions.
The aim of this work was to use the box area ratio method to determine the times to first division of individual Escherichia coli cells and thus examine the effect of heat injury to the inoculum and the salt content of the growth medium on lag time distributions of cells growing on an agar medium. We also determined the relative division time (RDT), which is a measure of the "work to be done" during lag, from direct observation of single cells.
|
|
|---|
Microscopic observation of growing cells.
Estimates of the times to first division of individual cells were made following the observation of cells growing on an agar surface, as described previously (22). Coverwell Press-Seal imaging chambers (Z 365866; Sigma-Aldrich Company, Gillingham, United Kingdom) were used to carry a thin layer of tryptone soya agar (TSA; Oxoid, Basingstoke, United Kingdom). These consisted of a 0.8-mm-thick silicone spacer gasket, incorporating a 20-mm circular well bonded to a polycarbonate microscope slide. TSA (27.5 µl) was added to the well and left to dry for 20 min under aseptic conditions. The agar was then inoculated with a drop (5 µl) of stationary-phase culture of E. coli, and a coverslip was placed over the well. Slides were incubated at room temperature (ca. 20°C) during observation of bacterial growth by phase contrast using a Nikon Microphot-SA microscope equipped with a x20 objective lens. Digital images were captured at 10- to 20-min intervals using a CoolSNAP-Pro cf monochrome camera incorporating a 1,392-by-1,040 pixel charge-coupled device and 12-bit 20-MHz digitization.
Determination of times to first division of individual cells.
The time from inoculation to the first cell division was taken as the measure of individual cell lag time. This was determined by automated processing of digital images to calculate the BAR, as described previously by Niven et al. (22). Digital images in TIFF format showing growing cells on agar were analyzed using Image-Pro Plus V4.5 image analysis software (Media Cybernetics). This identified each cell in the field of view and measured the parameters required to calculate the BAR (cell area, length, and width). These data, along with the x and y coordinates of each cell in the image, were downloaded into Microsoft Excel for further processing. A time course of BAR values was produced for each observed cell from the series of images. Since each image contained a large number of cells, a Microsoft Visual Basic program written in house was used to "track" individual cells through the sequence of images using the xy coordinates and to generate an Excel file containing time courses of area, length, and width for each cell (22). Another in-house Visual Basic program was used to calculate the BAR values as (length x width)/area and to identify the inflection points in plots of BAR against time, which had previously been shown to indicate the point of division (22). In order to increase the number of cells observed in each experiment, four fields of view were selected and the same fields were examined at each time point during an experiment. Data were thus generated for approximately 300 cells per experiment.
Where distributions of times to first cell division are shown, these are the results of a single experiment but are typical of the data produced by several replicate experiments. The Anderson-Darling test was used as a measure of how well data fitted each distribution using a maximum likelihood estimation. This was carried out using Minitab release 14.20 software (Minitab, Inc.).
Determination of RDT.
The time to first cell division divided by the cell size doubling time was taken as a measure of RDT. This was calculated for individual cells by image analysis, the time to first division being determined by the BAR method as described above, while the size doubling time was that of the two-dimensional cell area. The increase in cell area was observed to follow an exponential pattern after an initial period of accelerating growth rate, thus enabling the doubling time to be calculated from the gradient of the linear portion of plots of log cell area against time. The period of logarithmic growth was identified in the measured area data for individual cells, and the associated cell area doubling times were calculated using a Visual Basic program written for the purpose. The values quoted in the text refer to the average RDTs determined for all of the cells in the analyzed population.
|
|
|---|
![]() View larger version (10K): [in a new window] |
FIG. 1. Results of three replicate analyses showing distributions of times to first cell division following inoculation of cells grown to stationary phase at 37°C and subsequently grown at 20°C on agar. %N, percentage of analyzed cells.
|
|
View this table: [in a new window] |
TABLE 1. Anderson-Darling values for fits of various distribution functions to three replicate distributions of individual cell times to first division
|
![]() View larger version (14K): [in a new window] |
FIG. 2. Influence of preinoculation heating on distribution of times to first division. Stationary-phase cells grown at 37°C were heated at 50°C for the stated times prior to inoculation and analysis during growth at 25°C.
|
![]() View larger version (13K): [in a new window] |
FIG. 3. Effect of heating time at 50°C on mean times to first division (closed symbols) and the coefficient of variation of division times (open symbols). Stationary-phase cells grown at 37°C were heated prior to inoculation and growth at 20°C.
|
![]() View larger version (16K): [in a new window] |
FIG. 4. Influence of NaCl in the growth medium on the distribution of times to first division. Stationary-phase cells grown at 37°C were inoculated on agar medium containing the stated concentrations of NaCl (wt/vol), and the division times were determined during growth at 20°C.
|
![]() View larger version (12K): [in a new window] |
FIG. 5. Effect of NaCl in the growth medium on mean times to first division (closed symbols) and the coefficient of variation of division times (open symbols). Stationary-phase cells grown at 37°C were inoculated onto agar medium containing the stated concentrations of NaCl and incubated at 20°C.
|
The mean RDT for unstressed cultures was estimated as 1.8. In cultures that were heat stressed prior to inoculation, the RDT increased with the severity of the stress, rising to 4.3 after heating at 50°C for 80 min (Fig. 6a). This followed the pattern already described for the increase in the time to first division with the relative increase in time to first division and RDT being similar. The influence of postinoculation salt stress on RDT was comparatively slight, with no increase observed at salt concentrations less than 2.5% and an increase to 2.7 in the presence of 4% NaCl (Fig. 6b). In this case, the proportional increase in the RDT was lower than that of the time to first division. This indicated that the reduced growth rate of cells in the presence of elevated salt concentrations was a significant factor in increasing the time to first division, rather than the increased amount of work required to adapt to the new growth conditions. As with the time to first division, the standard deviations of all the RDT values for a given population were proportional to the mean value. The coefficients of variation for heat- and salt-stressed cells were 0.27 ± 0.04 and 0.27 ± 0.03, respectively. The coefficients of variation were thus relatively unaffected by the stress imposed on the cultures (Fig. 6).
![]() View larger version (17K): [in a new window] |
FIG. 6. Mean RDTs (closed symbols) and coefficients of variation (open symbols) for cultures subjected to heating at 50°C prior to inoculation (A) and to growth in the presence of NaCl (B). The RDTs were calculated as time to first division for individual cells divided by the doubling time of the cell two-dimensional area.
|
![]() View larger version (11K): [in a new window] |
FIG. 7. Mean cell two-dimensional areas at division for cultures subjected to heating at 50°C prior to inoculation (A) and to growth in the presence of NaCl (B).
|
|
|
|---|
There have been relatively few investigations of individual-cell lag using direct observation methods, and most of these have examined Listeria monocytogenes or Listeria innocua (5, 11, 25). The only other study that we are aware of in which lag times of stressed cells were measured directly using a slide culture method is that of Rasch et al. (25), who observed an increase in lag time and its variance in L. innocua cells growing on agar medium in the presence of the antimicrobial molecule reuterin (β-hydroxypropionaldehyde). Wu et al. (32) also determined lag time distributions in single cells of L. monocytogenes by microscopy but did not examine the effect of stress treatments.
The effects of sodium chloride on individual-cell lag times of E. coli have been examined using the microscope flow chamber and were found to increase with increasing sodium chloride concentration, but quantitative values were not presented (5). There are insufficient published data available to determine if lag times measured on a solid agar surface are similar to those in liquid under all conditions, but Wu et al. (32) found that lag times of unstressed E. coli cells in broth measured in the Bioscreen were longer than those of cells on agar. Interestingly, the mean lag times of unstressed E. coli measured using the flow chamber at room temperature were 5 h (21), somewhat longer than the 3 h observed with unstressed cells in the current study. To some extent, such differences in the absolute division times between studies may reflect the criteria used to objectively determine the point of cell division. The BAR method is based on the degree of flexibility at the point of attachment of daughter cells (22), while the flow chamber method requires physical separation of the cells with detachment of one daughter from the matrix (5). Liquid media are generally considered to provide more favorable conditions for recovery of injured cells than solid media, and growth rates on agar were shown to be slower than those in liquid (4). On the other hand, Guillier et al. (8) found no differences in the lag of stressed cells on agar and in broth media. Cells in foods are immobilized in a solid matrix rather than being present as planktonic cells in a liquid, so it seems important to determine whether the physical environment affects lag behavior of individual cells and under what conditions this may occur.
The E. coli distributions observed in the current study were slightly skewed toward higher division times, which is a common feature of lag distributions (6, 7, 9, 12, 20, 26, 31). It was not possible to assign any particular mathematical description to the distributions, since good fits to the data were found for a range of different functions and the histograms of division times lacked sufficient resolution to distinguish between them. A wide range of distributions have been examined in other studies, but no one distribution has proven to be best for all stress conditions and environments. These include, for example, gamma (1, 6, 19), exponential (6), Weibull (6), lognormal (12, 15), and extreme value (7, 9, 29) distributions. Given the diverse physiological effects of different stresses and the extreme sensitivity of stressed cells to small changes in their growth environment (13), it is not surprising that lag behavior fails to conform consistently to a mathematical ideal. In addition, it is likely that current methods for measuring individual cell lag times do not offer sufficient precision to enable different distribution functions to be adequately differentiated on the basis of experimental observations. However, this also implies that, for predictive modeling purposes, the exact function used to describe the distribution of lag times may not be critical as long as a reasonable basis exists for establishing values for the required variables. That being the case, it is perhaps more appropriate to have a method that is "good enough" and that has a satisfactory theoretical basis rather than attempting to use a range of different empirical fits to experimental data to accommodate different sets of conditions.
In our experiments, the mean single-cell lag time divided by the standard deviation was consistently around 0.2 for heat-shocked or osmotically stressed E. coli cells. Relatively mild stresses were examined of necessity because microscopic methods are not suitable for use after severe pretreatments when only a small proportion of the visible population is viable. Coefficients of variation for single-cell lag times reported in the literature tend to increase according to the severity of the stress, reaching values of 0.8 to >1.0 (6, 9). Guillier and Augustin (9) found that the logarithms of the means and standard deviations of individual cell lag times were linearly related for a wide range of conditions. These authors showed that the product of the growth rate and lag time at the single-cell level was constant and could be taken as representative of the physiological state of the cells.
If lag time is expressed as a multiple of the population doubling time during exponential growth under the same conditions, the resultant quantity, known as the relative lag time, is a means of comparing on the same basis lag times under different environmental conditions (18, 27, 28). Relative lag time can be regarded as a measure of the "work to be done" to prepare for cell division and is analogous to the physiological state parameter (h0) in the model of Baranyi and Roberts (2). To allow comparison of work to be done at the single-cell level, the time to first division was expressed as a multiple of the single-cell biomass doubling time, as estimated from the exponential increase in cell area. We have called this the RDT. Determinations of RDT indicated that to a large extent the extension of lag time of cells in the presence of NaCl was due to a reduction in the rate at which cells could grow, rather than additional time required to adapt to the conditions. This contrasted with heat-treated cells where the higher RTD values observed suggested greater amounts of work required to repair heat damage. For heat-injured cells, the cell size at division increased progressively with increasing degree of stress. For unheated cultures, the cells grew on average by a factor of 2.2 before first division, but after the most severe heat treatment this increased to 2.8, the area at division being 23% greater than for unstressed cells. Part of the work to be done may therefore be the need to achieve a larger cell size before division could take place. It is possible to speculate that the increase in size at division may be a result of a limiting step in the repair process that prevented cell division but not cell growth. The cells therefore grew to a larger size due to the additional time required complete this step. Kutalik et al. (11) also noticed that heat-shocked cells of L. monocytogenes needed to increase their length to 2.2. times their size at inoculation before division occurred. They also found that heat treatment caused a reduction in cell size in L. monocytogenes, which we did not observe with E. coli. So far as we are aware, there are no similar published data on relative lag based on direct measurements of single cells.
Despite the increasing number of studies of population heterogeneity among microbes, we remain ignorant of the underlying reasons for cell-to-cell variability. We propose that two classes of effect could give rise to such variability. The first is physiological or nongenetic heterogeneity within the population (3, 30). This arises from random differences in the distributions of molecules and macromolecular assemblies among cells that affect their capacity to withstand stress, to recover from stress, or to adapt to new growth conditions. The second class of effects is the stochastic variation in the interaction of cells with their environment both during and after exposure to stress. For example, in a heated suspension of cells the chances of a cell component being damaged or inactivated will be a probabilistic process (14). During recovery, differences in the microenvironment may also occur at longer distance scales, especially in nonhomogenous food material. The net result is that cells in a population may start out with inherent differences, and those differences may be compounded by variation in the degree of stress actually experienced by individual cells and by differences in their local recovery conditions.
Recent studies have begun to reveal how the scatter of properties within populations of cells is influenced by inoculum history and current growth conditions. This will lead to improvements in predictive models that can be used in quantitative risk assessments. This is particularly important when low numbers of pathogens in food pose a significant risk, as with highly infective pathogens such as E. coli O157 or with an organism such as Listeria monocytogenes that can grow to infective levels even under refrigeration. Further work is needed to understand the physiological basis of heterogeneity and the sequence of events that occur during lag. Direct observations of single-cell behavior, including measurements of physiological heterogeneity, under different conditions will undoubtedly contribute to that understanding.
We are grateful to Antony Vannard for the statistical analysis of distribution function fits and to Aline Métris for helpful discussions on lag time distributions.
Published ahead of print on 18 April 2008. ![]()
|
|
|---|
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»