 Software
 Open access
 Published:
Inferring timedependent population growth rates in cell cultures undergoing adaptation
BMC Bioinformatics volume 21, Article number: 583 (2020)
Abstract
Background
The population growth rate is an important characteristic of any cell culture. During sustained experiments, the growth rate may vary due to competition or adaptation. For instance, in presence of a toxin or a drug, an increasing growth rate indicates that the cells adapt and become resistant. Consequently, timedependent growth rates are fundamental to follow on the adaptation of cells to a changing evolutionary landscape. However, as there are no tools to calculate the timedependent growth rate directly by cell counting, it is common to use only end point measurements of growth rather than tracking the growth rate continuously.
Results
We present a computer program for inferring the growth rate over time in suspension cells using nothing but cell counts, which can be measured nondestructively. The program was tested on simulated and experimental data. Changes were observed in the initial and absolute growth rates, betraying resistance and adaptation.
Conclusions
For experiments where adaptation is expected to occur over a longer time, our method provides a means of tracking growth rates using data that is normally collected anyhow for monitoring purposes. The program and its documentation are freely available at https://github.com/Sandalmoth/ratrack under the permissive zlib license.
Background
Population growth rates are a core property of cell lines, and can be influenced by many factors. As such, in controlled experiments, cell growth rates, i.e. how quickly the population size changes over time, may correlate with the presence of drugs or toxins [1], temperature [2], or particular genetic changes in the cells [3], and other factors. The cell growth rate is an especially useful variable to feed into models of how a culture will progress [4,5,6]. Whereas well established mathematical definitions of the growth rate exist [7], it is not obvious how to estimate the growth rate from measurements of cell cultures. Suspension cells follow logistic growth to a decent approximation, and in that case growth rate could be derived via some form of logistic fitting.
The most common measurement for cell culture growth rate is the socalled population doubling time (PDT), i.e. the time it takes for a population to double its size [8]. The doubling time can be estimated from the population size at two points [8]. For cells growing exponentially this value is welldefined. However, the more a cell culture strays from exponential growth, the more inaccurate the PDT becomes a measure of population cell growth. It is no longer consistently applicable starting from any time. Exponential growth is a consequence of growth laws in individual cells causing regular division [9, 10]; the addition of factors such as cooperation or competition for resources creates a nonexponential population growth.
Under a different model of population growth, such as logistic growth, it is possible to separate the idea of growth rate from the decrease as the population approaches its carrying capacity. There are other approaches towards finding these growth parameters based on various methods, e.g. steepest descent optimization [11], and different means of Bayesian inference [12, 13]. Some models use highly detailed representations of the cell cycle [14]. However, these models generally do not consider a changing growth rate, and focus on inferring a static growth rate and carrying capacity. In reality, the growth rate sometimes changes over time, either because of densitydependent effects or because the cells are evolving. Examples of processes that result in modifications to the growth rate include adaptation to the presence of drugs [15, 16] and to lack of nutrients [17]. In general, any change in the evolutionary landscape might lead to a modification the growth rate.
Almost all available cell counting methods require fairly large sample sizes (\(> 10^5\) cells or \(> 1\) mL) [18]. Haemocytometer counting with some vital stain stands out as an exception, as it requires only a very minimal sample (5–10 \(\upmu \,\hbox {L}\), typically \(< 10^4\) cells). For monitoring a culture, especially in small volumes, requiring a large sample can be a limiting factor. Owing to its low cost and availability, haemocytometer counting has become the de facto standard in many applications. That, coupled with being measurable without significantly affecting an ongoing culture, makes it a prime choice for estimating the population growth rate.
The population size, i.e., the number of living cells, is inextricably tied to the growth rates, as the growth rate describes the positive component of the rate of change in population size. There is typically some basal death rate as well, although it is often difficult to separate apparent growth rate into the true growth rate and death rate. Unlike the population size, growth rate is not an observable variable and has to be inferred. By assuming a model of how cells grow, the growth rate can be deduced from the population size, which can be measured in ways that hardly affect the cell culture. Thus, it should be possible to track the (time dependent) growth rate in a cell culture over longer periods of time.
We developed a method to calculate the growth rates in a suspended cell lines. To use our approach, it is necessary to know the carrying capacity beforehand, but in return the growth rate is allowed to change over time. Our program is in principle agnostic to the type of measurement chosen, so long as it is proportional to the number of viable cells. However, it is structured primarily to work with cell counts (collected manually or by machine). Given the generality of the approach, it is possible to work with other measures as well, including optical densities, measurements of metabolites, etc.
Implementation
Cells in suspension culture typically follow logistic growth to a reasonable approximation. Logistic growth can be broken down into two elements: a given rate of division per cell, and a death rate that depends linearly on the population size [19]. Such a deathrate is equivalent to pairwise negative interactions between all cells in the population. Any genetic change or response in the cells might modify either of these parameters. We limit our model to changes in the growth rate, since death rates and the carrying capacity are coupled to nutrient availability, and are hence assumed to stay rather constant in any cell line.
The growth rate r(t) is represented as a piecewise linear function with evenly distributed segment lengths. This simple representation was chosen to reduce the risk of overfitting, as it requires as few as two (or even one, for a constant growth rate) parameters. The number of cells at any particular time, N(t), with carrying capacity K, can be obtained from
which we arrive at by substituting a timedependent growth rate into a standard logistic growth model [7]. Solving this Bernoulli differential equation [20] yields
Alternatively, the population curve can be obtained by simulating a logistic branching process [19], which is slower but likely more realistic, especially when the population size is low. In such a process, cells are modelled as individual particles that divide after some randomly distributed length of time. It is a type of model that very obviously connects to the reality of the situation. So, by extending the model in [19] with a timedependent growth rate, in the stochastic simulations, the growth rate is given by
and the death rate by
which is then simulated using the next reaction method [21]. Both simulation methods have been implemented; the stochastic simulator is recommended only for population sizes below \(2 \cdot 10^6\), as it is significantly slower for large populations since division and death events for all cells are simulated individually.
Experimental cell counting produces noisy and sparse timelines of the population size over time (e.g. Additional file 1: Figure S1). The piecewise linear growth rate cannot easily be obtained analytically from such data. To overcome this problem, we complete the population curve model above with imitation experimental noise. Approximate Bayesian computation (ABC) can thereafter be used to numerically estimate the growth rate.
The experimental noise consists of two parts:

1
Sampling noise—the error caused by counting cells in a small sample of the cell culture.

2
Counting noise—the errors induced by inaccurate counting.
The sampling noise is modelled by asserting that for well dispersed cells a sample contains a Poisson–distributed number of cells, thus sampling multiple times involves a nested series of Poisson random variables. This alone is known to underestimate the actual noise [22] as the counting noise is not included. Thus, to estimate the counting noise, we implemented a customizable set of filters, as disparate experimental setups may warrant different forms of noise. An easily motivated choice is a normally distributed counting noise, which is (approximately) consistent with each individual count having a certain chance of being incorrect. Dividing the noise into two parts is a pragmatic choice to make changing between counting methods easier, as the manual error rate is approximately known [23], whereas manufacturers of automatic counters commonly provide estimates.
This setup enables simulation from a piecewise linear growth rate over time using samples taken at discrete time points (Fig. 1). This forward model is inverted numerically using pyABC [24] and a sequential MonteCarlo (SMC) approach. Uniform priors over a userdefined range are used for the growth rate control points. Importantly, the deterministic simulator is wellbehaved also for negative growth rates. The output from ABC is samples from the posterior distribution, visualized using connected raincloud plots. We use raincloud plots as they combine a kernel density approximation, a boxplot and raw data plotting into a single figure [25], making it easy to interpret the results.
To summarize, the timedependent growth rate can be calculated numerically from the following:

1
a series of cell counts.

2
the dilution steps taken for each count.

3
the time when each of the counts was done.

4
an estimate of the carrying capacity.
1–3 are used to estimate the number of cells at a certain time taking measurement error into account, whereas 4 is a fundamental parameter of a certain cell type that has to be provided. This data is processed in a number of steps, which are automatically managed by a workflow system. See the “Usage and practical details” section for a more in depth explanation.
Results and discussion
Simulated data example
For an initial demonstration of the software capabilities, where the correct result is known, simulated data was created using the forward model assuming three different growth rate curves: a constant rate, an increasing rate, and a decreasing one. The growth rates for these simulated example data were then inferred using the program. All three were reproduced with sufficient accuracy to be qualitatively easily identifiable (Fig. 2). Quantitatively, some degree of inaccuracy is induced by sampling noise. Note for instance how all of the simulated samples in the constant growth rate case had a greater number of cells than they would on average, leading to a slightly overestimated growth rate. These sources of noise are however to a large degree captured by the width of the posterior distribution. Moreover, this type of error is expected from experiments as well.
Cell line tests
Two different experiments with cell lines were carried out to test viability in practice. First, KCL22 cells were cultured at two different concentrations of imatinib. KCL22 is a chronic myeloid leukaemia cell line reliant on the oncoprotein BcrAbl1 for survival, while imatinib is a BcrAbl1 inhibitor that effectively inhibits their growth [16]. With that in mind, the imatinib concentrations were selected such that the cells would grow, albeit slowly. The difference in growth rate can be clearly distinguished between the two concentrations (Fig. 3). Under normal conditions KCL22 cells have a growth rate of about 1 division per day (i.e. population size doubling per day). Imatinib has an inhibitory coefficient, IC50, of 240 nM (as measured by us). Thus, we expect to see growth rates of about 0.8 and 0.5 divisions per day, which is reasonably accurately represented given that IC50 measurements are typically not perfectly reproducible [26].
For the second cell line test, A batch of K562 cells was allowed to overgrow significantly. Each day, a subculture was taken and reseeded down to 1e5 cells/mL (a concentration at which they would normally grow well). These reseeded cultures were then counted every day. Allowing cells to overgrow typically damages their ability to reproduce, thus we expected that the reseeded cultures taken later on would grow slower initially. This was indeed observed; it could be seen that the longer cells had overgrown for, the worse their initial growth had become (Fig. 4). Additionally, the best growth occurred in the second subculture, which is consistent with the standard culturing procedure of splitting every other day. The observed decrease in growth rate in culture number two may have several explanations. First, in practice, as nutrient availability decreases with the increased number of living cells, we may expect a natural decrease in growth rate not accounted for by a simple logistic model of growth. Second, a slight overestimation of the carrying capacity parameter may cause such an artefact as the growing culture does not reach the expected population on time. During the four days for which they were tracked, all the subcultures recovered to grow at similar rates.
The endpoint posterior distributions are generally wider than the starting point. This occurs since when the population size approaches the carrying capacity, changes in the growth rate have an ever smaller effect on the relative change in population size. This is most evident in the extreme example of the final steady state at the carrying capacity, where changes in growth rate have no effects whatsoever. There is essentially less information about the growth rate as the population size approaches the carrying capacity. Note however that measurements of small population sizes are more noisy, which somewhat balances the qualitative effect on the size of the posterior distribution.
Given the width of the posterior distributions in Figs. 3 and 4 the magnitude of changes required to be clearly distinguishable is on the order of a 20% change in growth rate. This should be taken into account when designing experiments intended to use this method. For instance, to observe drug resistance, the drug dose used should inhibit growth of nonresistant cells by at least that much, so as to create an observable difference.
Conclusions
There are many tools which infer cellgrowth parameters from data, but we are not aware of any other tools for inferring a changing growth rate. The methodology is capable of distinguishing both absolute changes in growth rate, and differences specific to initial growth. The magnitude of the effects which could be identified is about a 20% change in growth rate or greater, though this depends on the specifics of the experiment. We did not test a decreasing growth rate on cells, as adaptation is normally positive for evolutionary reasons. Yet, since tracking decreasing rates worked for simulated data, this might be possible as well.
This work enables experiments wherein intentional or spontaneous changes to growth rate are of interest. For instance, it is possible to produce drug resistant cancer cells by culturing them with an inhibitor [16]. In cases where the adaptation happens quickly [15], this method could be used to follow the adaptation process. Another experiment made possible could be to examine the breakdown of some inhibitor, by observing how the growth changes as the concentration drops. This may be of use if the inhibitor is potent enough that the low concentrations are hard to measure directly.
Materials and methods
KCL22 and K562 cells were kindly provided by Prof. Leif Stenke (Karolinska Institutet, Solna, Sweden). They were cultured in \(37^{\circ }\hbox {C}\), 5% CO\(_2\) in RPMI1640 (Gibco, with Glutamax) with 10% HIFBS (Gibco) and 1% PenStrep (Gibco). Cell counting was performed using trypan blue exclusion on a LUNA2 cell counter (Logos Biosystems).
KCL22 with inhibitors
KCL22 cells were grown in under either 240 or 160 nM of imatinib, each in quadruplicate, in a 24 well plate. Cultures were initiated at 1e5 cells/mL in 1 mL medium, and were counted every day for 6 days.
K562 overgrowth
K562 cells were seeded to 1e5 cells/mL. Subcultures, diluted to 1e5 cells/mL, were taken daily for 8 days. All cultures were counted daily, as above.
Usage and practical details
The program runs on linux and has a simple command line usage powered by a snakemake [27] workflow that automatically processes input data, runs simulations and compiles an output report. It is also possible to manually run all intermediate steps if desired. As ABC is a computationally heavy approach, the deterministic simulator works on a regular laptop, but the stochastic simulator may require a more powerful machine to finish in a timely manner. Data input is prepared as two files, the first of which is a .csv featuring names, live cell counts, data collection times, and optionally information about the steps undertaken during sampling and possible subsequent dilution of the culture. These sampling and dilutions steps are used for simulating sampling noise. The second input file is a .toml file (a common human–readable config–file format) detailing simulation parameters. This file typically only requires a few parameters tuned from defaults. Most importantly: (1) The carrying capacity, i.e. the greatest population size the culture will reach if left alone. (2) How the named data should be grouped; grouping can be relevant if either several measurements are carried out on the same population, or if several populations are expected to behave similarly and one intends to elucidate the overall trend. Detailed annotated examples using the data from the tests below are available in the github repository. Figures 2, 3 and 4 were produced by the program, with minor additions made for clarity.
Usage example
Consider the following minimal example of cells growing logistically. Suppose we have cells growing in 1 mL of medium, which, if left alone, would reach a total of 3e6 cells (i.e. the carrying capacity). Suppose further we then recorded the following observations (Table 1) which, while presented here as a table for ease of reading, would be provided to the program as a regular .csv file. The measurements were done by first taking a \(10~\upmu \hbox {L}\) sample, mixing it with an equal volume of trypan blue dye, transferring the mixture to a haemocytometer, and finally counting the number of living cells in a \(0.4~\upmu \hbox {L}\) volume. In Table 1 the name column indicates that we are following one particular cell culture. The time is the time, in days, when the count was made (relative to the first one). count is selfexplanatory, whereas sample1 is the fraction of the total population taken in the first sampling step (i.e. the \(10~\upmu \hbox {L}\) sample), and sample2 is the fraction of the first sample counted in the second step, ignoring the dilution (as it does not change the number of cells in the sample). Note that \(0.4~\upmu \hbox {L}\) in a 1:1 dilution is \(0.2~\upmu \hbox {L}\) without dilution, and \(10~\upmu \hbox {L} \cdot 0.02 = 0.2~\upmu \hbox {L}\)
We run the program using the .toml file with parameters in Listing 1. Note that anything following a “#” are comments with no effect. We are assuming that the growth rate is somewhere between 0.01 and 3.0 divisions per cell and day, and that the growth rate function can be described as a piecewise linear function with either one (constant), two (linearly changing with time) or three (linearly changing with a central breakpoint) control points. Furthermore, we are running four simulations in parallel using the deterministic bernoulli equation simulator (Eq. 2). Finally, each individual count is assumed to be wrong 5% of the time (under a normal distribution approximation of the underlying binomial distribution).
Using data/minimal.csv and data/minimal.toml as input, the snakemake workflow automatically is used to run all simulation steps by calling snakemake results/minimal.pdf results/minimal.fit.csv which yields plots and a table of the inferred growth rates.
Availability and requirements

Project name: ratrack

Project home page: https://github.com/Sandalmoth/ratrack

Operating system: Linux

Programming language: Python3 and C++

Other requirements: Python 3.6 or higher, pyABC 9.14 or higher, snakemake 4.4 or higher. A complete list is provided as a conda environment file available on the project homepage.

License: zlib (permissive open source)

Any restrictions to use by nonacademics: Not applicable.
Availability of data and materials
The datasets supporting the conclusions of this article are included within the article (and its additional files). Simulation code (and datasets) are available at github.com/Sandalmoth/ratrack.
Abbreviations
 PDT:

Population doubling time
 HPDI:

Highest posterior density interval
 HIFBS:

Heatinactivated fetal bovine serum
 ABC:

Approximate bayesian computation
 SMC:

Sequential MonteCarlo
 IC50:

50% inhibitory coefficient
References
Zaharevitz DW, Holbeck SL, Bowerman C, Svetlik PA. COMPARE: a web accessible tool for investigating mechanisms of cell growth inhibition. J Mol Graph Model. 2002;20(4):297–303. https://doi.org/10.1016/s10933263(01)001267.
Watanabe I, Okada S. Effects of temperature on growth rate of cultured mammalian cells (L5178Y). J Cell Biol. 1967;32(2):309–23. https://doi.org/10.1083/jcb.32.2.309.
Skaggs BJ, Gorre ME, Ryvkin A, Burgess MR, Xie Y, Han Y, Komisopoulou E, Brown LM, Loo JA, Landaw EM, Sawyers CL, Graeber TG. Phosphorylation of the ATPbinding loop directs oncogenicity of drugresistant BCRABL mutants. Proc Natl Acad Sci. 2006;103(51):19466–71. https://doi.org/10.1073/pnas.0609239103.
Leder K, Foo J, Skaggs B, Gorre M, Sawyers CL, Michor F. Fitness conferred by BCRABL kinase domain mutations determines the risk of preexisting resistance in chronic myeloid leukemia. PLoS One. 2011;6(11):27682. https://doi.org/10.1371/journal.pone.0027682.
Mode CJ, Sleeman CK, Raj T. On the inclusion of self regulating branching processes in the working paradigm of evolutionary and population genetics. Front Genet. 2013;4:11. https://doi.org/10.3389/fgene.2013.00011.
Lindström HJG, de Wijn AS, Friedman R. Stochastic modelling of tyrosine kinase inhibitor rotation therapy in chronic myeloid leukaemia. BMC Cancer. 2019;19(1):508. https://doi.org/10.1186/s1288501956905.
Charlebois DA, Balázsi G. Modeling cell population dynamics. ISB. 2019;13(1–2):21–39. https://doi.org/10.3233/isb180470.
ATCC Animal Cell Culture Guide. https://www.lgcstandardsatcc.org/Documents/Marketing_Literature/Animal_Cell_Culture_Guide.aspx. Accessed 22 Oct 2020
Bell GI. Cell growth and division. Biophys J. 1968;8(4):431–44. https://doi.org/10.1016/s00063495(68)864987.
da Veiga J, Lafitte O, Schwartz L. A simple mathematical model for the growth and division of cells. Math Action. 2017;8(1):1–8. https://doi.org/10.5802/msia.10.
Juarez EF, Lau R, Friedman SH, Ghaffarizadeh A, Jonckheere E, Agus DB, Mumenthaler SM, Macklin P. Quantifying differences in cell line population dynamics using CellPD. BMC Syst Biol. 2016;10(1):92. https://doi.org/10.1186/s1291801603375.
Heydari J, Lawless C, Lydall DA, Wilkinson DJ. Fast Bayesian parameter estimation for stochastic logistic growth models. Biosystems. 2014;122:55–72. https://doi.org/10.1016/j.biosystems.2014.05.002.
Harris EA, Koh EJ, Moffat J, McMillen DR. Automated inference procedure for the determination of cell growth parameters. Phys Rev E. 2016;93(1):012402. https://doi.org/10.1103/physreve.93.012402.
Daukste L, Basse B, Baguley BC, Wall DJN. Mathematical determination of cell population doubling times for multiple cell lines. Bull Math Biol. 2012;74(10):2510–34. https://doi.org/10.1007/s1153801297647.
Yuan H, Wang Z, Gao C, Chen W, Huang Q, Yee JK, Bhatia R, Chen W. BCRABLGene expression is required for its mutations in a novel KCL22 cell culture model for acquired resistance of chronic myelogenous leukemia. J Biol Chem. 2009;285(7):5085–96. https://doi.org/10.1074/jbc.m109.039206.
Tang C, Schafranek L, Watkins DB, Parker WT, Moore S, Prime JA, White DL, Hughes TP. Tyrosine kinase inhibitor resistance in chronic myeloid leukemia cell lines: Investigating resistance pathways. Leuk Lymphoma. 2011;52(11):2139–47. https://doi.org/10.3109/10428194.2011.591013.
Gray DA, Dugar G, Gamba P, Strahl H, Jonker MJ, Hamoen LW. Extreme slow growth as alternative strategy to survive deep starvation in bacteria. Nat Commun. 2019;10(1):890. https://doi.org/10.1038/s41467019087198.
Butler M, Spearman M, Braasch K. Monitoring cell growth, viability, and apoptosis. Methods Mol Biol (Clifton, N.J.). 2014;1104:169–92. https://doi.org/10.1007/9781627037334_12.
Lambert A. The branching process with logistic growth. Ann Appl Probab. 2005;15(2):1506–35. https://doi.org/10.1214/105051605000000098.
Weisstein EW. Bernoulli differential equation. https://mathworld.wolfram.com/BernoulliDifferentialEquation.html. Accessed 20 Oct 2020
Anderson DF. A modified next reaction method for simulating chemical systems with time dependent propensities and delays. J Chem Phys. 2007;127(21):214107. https://doi.org/10.1063/1.2799998.
Nielson L, Smyth G, Greenfield P. Hemacytometer cell count distributions: implications of nonPoisson behavior. Biotechnol Prog. 1991;7(6):560–3. https://doi.org/10.1021/bp00012a600.
CadenaHerrera D, EsparzaDe Lara JE, RamírezIbañez ND, LópezMorales CA, Pérez NO, FloresOrtiz LF, MedinaRivero E. Validation of three viablecell counting methods: manual, semiautomated, and automated. Biotechnol Rep. 2015;7:9–16. https://doi.org/10.1016/j.btre.2015.04.004.
Klinger E, Rickert D, Hasenauer J. pyABC: distributed, likelihoodfree inference. Bioinformatics. 2018;34(20):3591–3. https://doi.org/10.1093/bioinformatics/bty361.
Allen M, Poggiali D, Whitaker K, Marshall TR, Kievit RA. Raincloud plots: a multiplatform tool for robust data visualization. Wellcome Open Res. 2019;4(63):63. https://doi.org/10.12688/wellcomeopenres.15191.1.
Soverini S, Rosti G, Iacobucci I, Baccarani M, Martinelli G. Choosing the best secondline tyrosine kinase inhibitor in imatinibresistant chronic myeloid leukemia patients harboring bcrabl kinase domain mutations: How reliable is the IC 50 ? The Oncologist. 2011;16(6):868–76. https://doi.org/10.1634/theoncologist.20100388.
Koster J, Rahmann S. Snakemake—a scalable bioinformatics workflow engine. Bioinformatics. 2012;28(19):2520–2. https://doi.org/10.1093/bioinformatics/bts480.
Acknowledgements
We thank Prof. Leif Stanke, Karolinksa Institutet, for the donation of the K562 and KCL22 cells and Dr. Elham Yekaterini Karin, Karolinska Institutet, for preparing and sending the cells.
Funding
Open Access funding provided by Linnaeus University. This work was supported by The Swedish Cancer Society (Cancerfonden) [CAN 2015/387 and CAN 2018/362 to RF]. The funding body had no role in the design and execution of the study or in writing the manuscript.
Author information
Authors and Affiliations
Contributions
The study was directed by RF. HJGL wrote software and performed the tests. HJGL wrote the manuscript with assistance from RF. Both authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Figure S1.
Cell count timelines collected from K562 cells.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Lindström, H.J.G., Friedman, R. Inferring timedependent population growth rates in cell cultures undergoing adaptation. BMC Bioinformatics 21, 583 (2020). https://doi.org/10.1186/s12859020038877
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12859020038877