The described method provides a means to capture dynamic bacterial growth in a continuous, high-throughput manner by utilizing a 96-well format reader that takes multiple optical density measurements at various time intervals (from minutes to hours to days). The growth curves of an assortment of E. coli strains expressing various genomes can be precisely acquired in a single experiment (Figure 1A). In comparison to the described method, the traditional method (the CFU assay) generally requires longer sampling time intervals (Figure 1B) and is labor intensive if multiple culture is required. In addition, each culture sampling time for the CFU assay requires the use of a small volume of cell culture that cannot be used for repeated measurements. Additionally, limited numbers of time points for detection may miss the sampling points that are needed for the correct calculation of the growth rate and the maximal OD600. However, the CFU assay has a lower detection limit, 103 - 104 cells/mL, making it at least two orders of magnitude more sensitive than the optical turbidity assay (OD600= 0.001, which has a detection limit of approximately 105 - 106 cells/mL). The described method provides a practical and efficient experimental tool for systems-level studies on growth dynamics.
A key advantage of this protocol is that it takes repeated measurements from a common glycerol stock. Having multiple cell cultures at various dilution rates in early exponential growth phase is highly useful because it allows researchers to determine the culture time at any timescale and avoids having saturated cultures due to unpredictable events. Repeated culturing and measurements of the common stocks, which were prepared from one of the five cultures initiated using different dilution rates as described in Protocol section 2.2 (Figure 2A), presented highly similar outcomes of the growth dynamics analysis. Multiple growth curves of repeated and independent measurements (N= 6) overlapped well (Figure 2B, top panel), with little variance (Figure 2B, bottom panel), particularly during the exponential phase (Figure 2B, shadowed area). These results suggest that the described method delivers highly reproducible results.
Estimating the global bias of the 96-well reads is highly recommended. In the present method, the global bias is estimated by monitoring the growth of the wild-type E. coli strain W3110 as described in Protocol section 5. For example, the growth rates calculated from the growth curves showed significant variation among the 96 wells (Figure 3A), representing the global fluctuation derived from both the data manipulation and the device errors. These errors can be masked by loading the same sample into multiple wells at varied locations on the plate. By contrast, the maximal OD600 showed a clear location-dependent bias for the results, as gradual changes were observed along the vertical direction, that is, the maximal OD600 gradually increased from the wells of column 2 to those of column 11 on the 96-well plate (Figure 3B). To avoid the biased results, the same sample is recommended to be loaded into the wells located on the both sides of the 96-well plate for repeated measurements. The locational bias in maximal OD600 values is assumed to result from the variation in evaporation rates, which are determined by both heating and sealing efficiencies. Compared with the wells in the middle of the 96-well microplate, the wells located at the edges tended to show relatively higher values, reflecting the varied evaporation rates resulting from the sealing effect of the microplate. In addition, since other 96-well microplates subjected to the same bias test showed similar directional changes based on where the well was located on the plate, the extent of the bias might depend on the plate reader. Thus, it is important to evaluate the global bias of the 96-well reads in studies that determine maximal population density. Commercially available plate readers all have their specific global bias caused by the heating and shaking efficiencies, and the 96-well microplates are largely differentiated by sealing efficiency.
This method is designed for the systematic analysis of bacterial growth, and it can be used to analyze two major parameters, i.e., growth rate and cell density, that are commonly applied in a wide variety of fields, such as systems biology, evolution, and ecology16,17. The manual calculation using spreadsheet is introduced to show the processing of the raw reads to produce the evaluated parameters of growth rate and cell density. The hourly growth rates are first calculated from the raw reads of the growth curve (Figure 4A) according to Protocol section 4.6 as a data set of sequential changes in growth rates (Figure 4B). The means and standard deviations of every five consecutive growth rates are acquired sequentially (Figure 4C-D), and the largest mean growth rate with the smallest standard deviation is defined as the maximal growth rate of the growth curve (Figure 4C-D, highlighted). The maximal cell density (OD600) is analyzed similarly (Figure 4E-F). The calculated growth rate and cell density are subsequently subjected to further analysis. Note that the described calculations can be automatically performed using computational platforms (programming languages), e.g., R and Python.

Figure 1 Growth Curves Acquired with Diverse Methods. (A) Multiple acquisitions of growth curves. Growth curves of four different E. coli strains expressing various genomes were obtained in a high-throughput manner. The E. coli cells were grown in minimal medium M63, and the changes in OD600 reads were recorded by a microplate reader at 1-h intervals. The strains from left to right are the wild type E. coli W3110 (Number 0) and its reduced genomes numbers 1, 10 and 18 (described previously5), respectively. (B) Growth curve obtained by using the CFU assay. The wild-type E. coli cells (Number 0) were cultured and were sampled at the indicated time points over several hours. The growth changes were estimated using the CFU assay. Open circles indicate the recording or sampling points. Please click here to view a larger version of this figure.

Figure 2. Reproducible Repeated Measurements. (A) Overnight cultures of varied dilution rates. The E. coli cells were inoculated in tubes at multiple dilution rates, which varied from 1- to 10,000-fold, as indicated. The cells were cultured at 37 °C with shaking at 200 rpm for approximately 12 h. The final OD600 values of the five cell cultures (from left to right) were 1.51, 0.74, 0.05, 0.008 and 0.001, respectively. The culture with OD600 of 0.05 was used for the common glycerol stock. (B) Highly reproducible measurements. Six independent growth curves (top, gray lines) of the same E. coli strain were acquired and calculated from six positions on the microplate and two vials of the common glycerol stock. The black line represents the average of the six growth curves. The coefficients of variation (CV) of the six repeated measurements were calculated (bottom), and the steady low points are shadowed in pink. Please click here to view a larger version of this figure.

Figure 3 Global Bias of the 96-well Format. (A) Heat map of growth rates. The E. coli cell stock was suspended in minimal medium M63 and loaded equally in the 96 wells of the microplate (represented as a 12 × 8 table). The changes in OD600 were obtained using the microplate reader, as described in section 5 of the protocol. The gradation from white to red indicates the growth rates from 0.7 to 0.9 h-1. (B) Heat map of the maximal OD600. The gradation from white to red indicates the maximal OD600 values from 0.9 to 1.2. The calculations for the growth rate and maximal OD600 are as described in Figure 4. Please click here to view a larger version of this figure.

Figure 4 Analysis of the Growth Curves. (A) Raw data of growth curves. The hourly OD600 reads are plotted against the culture time. Open circles indicate the sampling times. (B) Temporal changes in growth rates. Hourly growth rates between two consecutive OD600 reads are calculated according to the equation described in 4.6. (C) Mean growth rates. The mean of five consecutive growth rates are calculated to give smooth growth rates. (D) Standard deviations of mean growth rates. The standard deviation of the same five consecutive growth rates is calculated. The maximal growth rate with the smallest standard deviation is highlighted in red. (E) Mean OD600. The mean of five consecutive OD600 values is calculated to determine the smooth OD600. (F) Standard deviations of mean OD600 values. Standard deviation of the same five consecutive OD600 values is calculated. The maximal OD600 with the smallest standard deviation is highlighted in red. Please click here to view a larger version of this figure.