 Software
 Open Access
 Published:
Growthcurver: an R package for obtaining interpretable metrics from microbial growth curves
BMC Bioinformatics volume 17, Article number: 172 (2016)
Abstract
Background
Plate readers can measure the growth curves of many microbial strains in a highthroughput fashion. The hundreds of absorbance readings collected simultaneously for hundreds of samples create technical hurdles for data analysis.
Results
Growthcurver summarizes the growth characteristics of microbial growth curve experiments conducted in a plate reader. The data are fitted to a standard form of the logistic equation, and the parameters have clear interpretations on populationlevel characteristics, like doubling time, carrying capacity, and growth rate.
Conclusions
Growthcurver is an easytouse R package available for installation from the Comprehensive R Archive Network (CRAN). The source code is available under the GNU General Public License and can be obtained from Github (Sprouffske K, Growthcurver sourcecode, 2016).
Background
By tracking bacterial growth over time, important populationlevel 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 [1–3], including experimental evolution [4–7]. A variety of methods have been used to obtain metrics from such growth curves [1, 3, 4, 8–13]. Older methods relied on manually plotting the cell count or absorbance readings over time on semilog graph paper to obtain metrics like the maximum growth rate [8], an approach which has been computationally mirrored recently by GrowthRates [12]. Other recent methods fit the growth data to a variety of parametric growth models [11], but have a strong focus on generating dose response curves.
Here, we fit growth curve data to the standard form of the logistic equation common in ecology and evolution [14, 15] whose parameters (the growth rate, the initial population size, and the carrying capacity) provide meaningful populationlevel information with straightforward 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 workflow in the vignette.
Implementation
We developed an opensource R package, Growthcurver, to obtain a variety of easilyinterpretable metrics to summarize microbial growth curve data. Growthcurver is available from CRAN and the source code is available under the GNU General Public License.
Carrying capacity and growth rate
Growthcurver fits a basic form of the logistic equation common in ecology and evolution [14, 15] to experimentallyobtained growth curve data. The logistic equation gives the number of cells N _{ t } at time t.
Here, the population size at the beginning of the growth curve is given by N _{0}. 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 N _{0} for the growth curve data using the implementation of the nonlinear leastsquares LevenbergMarquardt algorithm [16] available in the minpack.lm R package [17]. The carrying capacity and growth rate values (K and r) are particularly useful for summarizing and comparing the growth dynamics of strains.
Area under the curve
Growthcurver computes the area under the logistic curve, which provides a metric (the logistic AUC) that integrates information from the logistic parameters (k, r, and N _{0}). We implemented this feature by evaluating the definite integral of the fitted logistic equation from time 0 to a userdefined 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 [18].
Doubling time
The 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 t _{ DT } possible for the population according to
[15], which occurs when the population has no restrictions on its growth (e.g., when it is far from the carrying capacity).
Results and discussion
Usage
Growthcurver 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.
In 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.
Example 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 [19].
Correlation between metrics
To 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 15758500Gf) supplemented with 1 μg/L glucose (‘DM1000’ media). For each measurement, we started a preculture in 2 mL DM1000 media and grew it at 37 °C with shaking for 24 h. We diluted the overnight culture 10000fold, transferred it into 96well 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.
The 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 nonparametric 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}).
Conclusions
Here, 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.
Growthcurver is available for installation from CRAN [20] or Github [21]. The accompanying vignette [19] provides stepbystep 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.
The highthroughput 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 [22]. Growthcurver allows for the straightforward analysis and interpretation of growth curve data collected in a highthroughput manner using plate readers.
Availability and requirements

Project name: Growthcurver

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

Operating system: Platform independent

Programming language: R (version 3.2.2)

Other requirements: The R Project for Statistical Computing

License: GNU General Public License

Any restrictions to use by nonacademics: According to GNU General Public License
Abbreviations
 AUC:

area under the curve
 CRAN:

comprehensive R archive network
References
 1
Buchanan RL, Whiting RC, Damert WC. When is simple good enough: a comparison of the Gompertz, Baranyi, and threephase linear models for fitting bacterial growth curves. Food Microbiol. 1997; 14:313–26.
 2
Toussaint M, Conconi A. Highthroughput and sensitive assay to measure yeast cell growth: a bench protocol for testing genotoxic agents. Nat Protoc. 2006; 1(4):1922–8.
 3
Koseki S, Nonaka J. Alternative approach to modeling bacterial lag time, using logistic regression as a function of time, temperature, pH, and sodium chloride concentration. Appl Environ Microbiol. 2012; 78(17):6103–112.
 4
Vasi F, Travisano M, Lenski RE. Longterm experimental evolution in Escherichia coli. II. Changes in lifehistory traits during adaptation to a seasonal environment. Am Nat. 1994; 144:432–56.
 5
Novak M, Pfeiffer T, Ackermann M, Bonhoeffer S. Bacterial growth properties at low optical densities. Anton Leeuw Int J Gen Mol Microbiol. 2009; 96(3):267–74.
 6
Lindsey HA, Gallie J, Taylor S, Kerr B. Evolutionary rescue from extinction is contingent on a lower rate of environmental change. Nature. 2013; 494(7438):463–7.
 7
Leiby N, Marx CJ. Metabolic erosion primarily through mutation accumulation, and not tradeoffs, drives limited evolution of substrate specificity in Escherichia coli. PLoS Biol. 2014; 12(2):1001789.
 8
Monod J. The growth of bacterial cultures. Annu Rev Microbiol. 1949; 3(1):371–94.
 9
Zwietering MH, Jongenburger I, Rombouts FM, van ’t Riet K. Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990; 56(6):1875–1881.
 10
Baranyi J, Roberts TA, McClure P. A nonautonomous differential equation to model bacterial growth. Food Microbiol. 1993; 10:43–59.
 11
Kahm M, Hasenbrink G, LichtenbergFraté H, Ludwig J, Kschischo M. Grofit: fitting biological growth curves with R. J Stat Softw. 2010; 33(7):1–21.
 12
Hall BG, Acar H, Nandipati A, Barlow M. Growth rates made easy. Mol Biol Evol. 2014; 31(1):232–8.
 13
Bukhman YV, DiPiazza NW, Piotrowski J, Shao J, Halstead AGW, Bui MD, Xie E, Sato TK. Modeling microbial growth curves with GCAT. BioEnergy Res. 2015; 8(3):1022–1030.
 14
Crow JF, Kimura M. An Introduction to Population Genetics Theory. Caldwell, NJ: Harper and Row, Publishers, Inc.; 1970.
 15
Rockwood LL. Introduction to Population Ecology, 2nd ed. West Sussex: Wiley; 2015.
 16
Moré JJ. The LevenbergMarquardt algorithm: implementation and theory. In: Numerical Analysis (Proc. 7th Biennial Conf., Univ. Dundee, Dundee, 1977). Springer: 1978. p. 105–16.
 17
Elzhov TV, Mullen KM, Spiess AN, Bolker B. minpack.lm: R Interface to the LevenbergMarquardt Nonlinear LeastSquares Algorithm Found in MINPACK, Plus Support for Bounds. https://cran.rproject.org/web/packages/minpack.lm. Accessed 24 Mar 2016.
 18
Hasenbrink G, Schwarzer S, Kolacna L, Ludwig J, Sychrova H, LichtenbergFraté H. Analysis of the mKir2.1 channel activity in potassium influx defective Saccharomyces cerevisiae strains determined as changes in growth characteristics. FEBS Lett. 2005; 579(7):1723–1731.
 19
Sprouffske K. Using Growthcurver. https://cran.rproject.org/web/packages/growthcurver/vignettes/Growthcurvervignette.html. Accessed 24 Mar 2016.
 20
Sprouffske K. Growthcurver: An R package for obtaining interpretable metrics from microbial growth curves. https://cran.rproject.org/web/packages/growthcurver/. Accessed 24 Mar 2016.
 21
Sprouffske K. Growthcurver sourcecode. https://github.com/sprouffske/growthcurver. Accessed 24 Mar 2016.
 22
ConcepciónAcevedo J, Weiss HN, Chaudhry WN, Levin BR. Malthusian parameters as estimators of the fitness of microbes: A cautionary tale about the low side of high throughput. PLOS ONE. 2015; 10(6):0126915.
Acknowledgements
The authors thank Sinisa Bratulic, Athena Chu, Rumen Kostadinov, and Macarena TollRiera for advice. KS acknowledges support through the Forschungskredit program of the University of Zurich, grant K74301030. AW acknowledges support through Swiss National Science Foundation grant 31003A_146137, as well as through the University Priority Research Program in Evolutionary Biology at the University of Zurich.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
KS conceived of the project, developed the package, did the growth curve experiments, and did the data analysis. KS wrote the paper with AW. Both authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Growth curve
 Logistic
 Experimental evolution