Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves

1By tracking bacterial growth over time, important population-level information can be assessed, including doubling time and carrying capacity. Typical experiments entail measuring bacterial cell density at a series of time intervals, and then fitting these observations to an exponential growth model. Such measures can be made in parallel in a plate reader, and may result in hundreds or thousands of absorbance measurements over the course of 24 h. The resulting growth curves are commonly used in a variety of microbiological experiments, including experimental evolution. A variety of methods have been used to obtain metrics from such growth curves. Older methods relied on manually plotting the cell count or absorbance readings over time on semi-log graph paper to obtain metrics like the maximum growth rate, an approach which has been computationally mirrored recently by GrowthRates. Other recent methods fit the growth data to a variety of parametric growth models, but have a strong focus on generating dose response curves.

2Here, we fit growth curve data to the standard form of the logistic equation common in ecology and evolution whose parameters (the growth rate, the initial population size, and the carrying capacity) provide meaningful population-level information with straight-forward biological interpretation. We implemented this as the R package Growthcurver, available for download at Comprehensive R Archive Network (CRAN), and provide a simple data analysis work-flow in the vignette.

3We developed an open-source R package, Growthcurver, to obtain a variety of easily-interpretable metrics to summarize microbial growth curve data. Growthcurver is available from CRAN and the source code is available under the GNU General Public License.

4Growthcurver fits a basic form of the logistic equation common in ecology and evolution to experimentally-obtained growth curve data. The logistic equation gives the number of cells Nt at time t.

5Here, the population size at the beginning of the growth curve is given by N0. The maximum possible population size in a particular environment, or the carrying capacity, is given by K. The intrinsic growth rate of the population, r, is the growth rate that would occur if there were no restrictions imposed on total population size. Growthcurver finds the best values of K, r, and N0 for the growth curve data using the implementation of the non-linear least-squares Levenberg-Marquardt algorithm available in the minpack.lm R package. The carrying capacity and growth rate values (K and r) are particularly useful for summarizing and comparing the growth dynamics of strains.

6Growthcurver computes the area under the logistic curve, which provides a metric (the logistic AUC) that integrates information from the logistic parameters (k, r, and N0). We implemented this feature by evaluating the definite integral of the fitted logistic equation from time 0 to a user-defined time t. Growthcurver also computes an empirical AUC by summing the areas of the trapezoids made up by connecting consecutive data points of cell counts (or absorbance measurements) from time 0 to time t, similar to what was done by.

7The doubling time, also called the generation time, of a population is the time it takes for the number of individuals (or the absorbance reading) to double. Growthcurver computes the fastest doubling time tDT possible for the population according to, which occurs when the population has no restrictions on its growth (e.g., when it is far from the carrying capacity).

8Growthcurver can be used to compute metrics for growthcurves in two modes: individually for a single sample, or in batch mode for an entire plate of samples.

9In both cases, Growthcurver requires just two vectors of data for each growth curve: time measurements and absorbance readings taken at those times. The time measurements’ unit determines the unit of the metrics returned by Growthcurver (e.g., if the input data are in minutes, then the doubling time returned by Growthcurver is in minutes and the growth rate is in minutes −1). The absorbance of the media should be subtracted from the absorbance readings, and Growthcurver provides an option to do that automatically.

10Example calls to Growthcurver can be found in the R documentation for the functions SummarizeGrowth and SummarizeGrowthByPlate, and extensive sample code for obtaining, checking, and interpreting growth curves is provided in the accompanying vignette.

11To compare different metrics computed by Growthcurver, we measured 937 growth curves of 33 different Escherichia coli K12 strains in 200 μL of Davis Minimal broth (Fluka 15758-500G-f) supplemented with 1 μg/L glucose (‘DM1000’ media). For each measurement, we started a pre-culture in 2 mL DM1000 media and grew it at 37 °C with shaking for 24 h. We diluted the overnight culture 10000-fold, transferred it into 96-well plates (TPP 92096), and placed it in a plate reader (Tecan Infinite 200 Pro). We grew the resulting 96 populations at 37 °C with shaking for 24 h, and took optical density (OD) readings at a wavelength of 600 nm every 10 min. We removed the contribution of the media from the optical density by subtracting the minimum observed OD value for each growth curve from the rest of the values. Finally, we used Growthcurver to obtain the growth curve metrics for these populations.

12The growth rate is often used to summarize growth curve data, and we wondered to what extent the growth rate correlates with other metrics that summarize growth (Fig. 1). The growth rate correlates perfectly with the doubling time (Spearman, ρ=−1,p<2.2×10−16), which is unsurprising since the definition of doubling time relies on growth rate (see Implementation). The growth rate also correlates with the initial population size (Spearman, ρ=−0.88,p<2.2×10−16), the area under the curve (Spearman, ρ=0.57,p<2.2×10−16), and the carrying capacity (Spearman, ρ=0.32,p<2.2×10−16). The area under the curve is a promising metric to summarize a growth curve because it integrates the contributions of the initial population size, growth rate, and carrying capacity into a single value, and emphasizes growth rate (Spearman correlation between area under the curve and growth rate, ρ=0.57,p<2.2×10−16) and carrying capacity (Spearman correlation between area under the curve and carrying capacity, ρ=0.81,p<2.2×10−16). A viable non-parametric approach is to use the empirical area under the curve, which is correlated with the area under the curve (Spearman, ρ=1,p<2.2×10−16).

13Here, we have described the R package Growthcurver that provides several growth curve metrics with intuitive biological interpretation, including the growth rate, the doubling time, the carrying capacity, and the area under the logistic growth curve which integrates the contributions of the other metrics into a single value.

14Growthcurver is available for installation from CRAN or Github. The accompanying vignette provides step-by-step examples for analyzing a single sample and a plate of samples, as well as information on preparing the input data for analysis. Novice R users should be able to use the example code provided with few or no changes to analyze their own growth curve data.

15The high-throughput nature of obtaining replicate growth curves in plates allows us to easily characterize hundreds of growth curves, even though caution is needed to interpret these metrics in terms of their effects on competitive fitness. Growthcurver allows for the straightforward analysis and interpretation of growth curve data collected in a high-throughput manner using plate readers.

16Project name: GrowthcurverProject home page: http://github.com/sprouffske/growthcurverOperating system: Platform independentProgramming language: R (version 3.2.2)Other requirements: The R Project for Statistical ComputingLicense: GNU General Public LicenseAny restrictions to use by non-academics: According to GNU General Public License

17Project name: Growthcurver

18Project home page: http://github.com/sprouffske/growthcurver

19Operating system: Platform independent

20Programming language: R (version 3.2.2)

21Other requirements: The R Project for Statistical Computing

22License: GNU General Public License

23Any restrictions to use by non-academics: According to GNU General Public License