Document Detail

Fractional proliferation: a method to deconvolve cell population dynamics from single-cell data.
Jump to Full Text
MedLine Citation:
PMID:  22886092     Owner:  NLM     Status:  MEDLINE    
Abstract/OtherAbstract:
We present an integrated method that uses extended time-lapse automated imaging to quantify the dynamics of cell proliferation. Cell counts are fit with a quiescence-growth model that estimates rates of cell division, entry into quiescence and death. The model is constrained with rates extracted experimentally from the behavior of tracked single cells over time. We visualize the output of the analysis in fractional proliferation graphs, which deconvolve dynamic proliferative responses to perturbations into the relative contributions of dividing, quiescent (nondividing) and dead cells. The method reveals that the response of 'oncogene-addicted' human cancer cells to tyrosine kinase inhibitors is a composite of altered rates of division, death and entry into quiescence, a finding that challenges the notion that such cells simply die in response to oncogene-targeted therapy.
Authors:
Darren R Tyson; Shawn P Garbett; Peter L Frick; Vito Quaranta
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Extramural     Date:  2012-08-12
Journal Detail:
Title:  Nature methods     Volume:  9     ISSN:  1548-7105     ISO Abbreviation:  Nat. Methods     Publication Date:  2012 Sep 
Date Detail:
Created Date:  2012-09-04     Completed Date:  2012-11-13     Revised Date:  2013-07-12    
Medline Journal Info:
Nlm Unique ID:  101215604     Medline TA:  Nat Methods     Country:  United States    
Other Details:
Languages:  eng     Pagination:  923-8     Citation Subset:  IM    
Affiliation:
Department of Cancer Biology, Vanderbilt University School of Medicine, Nashville, Tennessee, USA. darren.tyson@vanderbilt.edu
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
Descriptor/Qualifier:
Carcinoma, Squamous Cell / drug therapy,  pathology*
Cell Count
Cell Proliferation
Humans
Lung Neoplasms / drug therapy,  pathology*
Microscopy, Video / methods*
Protein Kinase Inhibitors / pharmacology
Protein-Tyrosine Kinases / antagonists & inhibitors
Single-Cell Analysis / methods*
Structure-Activity Relationship
Time Factors
Tumor Cells, Cultured
Grant Support
ID/Acronym/Agency:
2 UL1 TR000445-06/TR/NCATS NIH HHS; 5U54 CA113007-07/CA/NCI NIH HHS; DK05840/DK/NIDDK NIH HHS; KL2 RR024977/RR/NCRR NIH HHS; KL2 TR000446/TR/NCATS NIH HHS; P30 CA68485/CA/NCI NIH HHS; TL1 RR024978/RR/NCRR NIH HHS; TL1 TR000447/TR/NCATS NIH HHS; U54 CA113007/CA/NCI NIH HHS; UL1 RR024975/RR/NCRR NIH HHS; UL1 RR024975-01/RR/NCRR NIH HHS; UL1 TR000445/TR/NCATS NIH HHS
Chemical
Reg. No./Substance:
0/Protein Kinase Inhibitors; EC 2.7.10.1/Protein-Tyrosine Kinases
Comments/Corrections

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
Journal ID (nlm-journal-id): 101215604
Journal ID (pubmed-jr-id): 32338
Journal ID (nlm-ta): Nat Methods
Journal ID (iso-abbrev): Nat. Methods
ISSN: 1548-7091
ISSN: 1548-7105
Article Information
Download PDF

License:
nihms-submitted publication date: Day: 15 Month: 8 Year: 2012
Electronic publication date: Day: 12 Month: 8 Year: 2012
Print publication date: Month: 9 Year: 2012
pmc-release publication date: Day: 01 Month: 3 Year: 2013
Volume: 9 Issue: 9
First Page: 923 Last Page: 928
PubMed Id: 22886092
ID: 3459330
DOI: 10.1038/nmeth.2138
ID: NIHMS397921

Fractional Proliferation: A method to deconvolve cell population dynamics from single-cell data
Darren R. Tyson
Shawn P. Garbett
Peter L. Frick
Vito Quaranta
Department of Cancer Biology, Vanderbilt University, 2220 Pierce Ave., Nashville, TN 37232
Correspondence: TO WHOM CORRESPONDENCE SHOULD BE ADDRESSED: Darren R. Tyson, darren.tyson@vanderbilt.edu

INTRODUCTION

Proliferation is a fundamental property of living cells. Altering or controlling cell—proliferation often by chemical compounds—is a major goal of several disciplines, including oncology, tissue engineering, and developmental biology. Many current proliferation assays rely on surrogate measurements of cell number (e.g. total ATP or DNA content) rather than direct cell counts1. Moreover, since cells are not directly visualized in most proliferation assays, their states (dividing, quiescent, apoptotic) are unknown; a given proliferation curve could be the consequence of one of many combination of states. Proliferation assays are also typically performed with very few time points and rarely account for proliferation dynamics. Flow cytometry provides measurements at the single-cell level, including: cell cycle position 2, expression of proliferative markers, fractions of dead or dying cells 3, and division tracking by dye dilution 4. However, adherent cells must be detached for end point analysis and single cells cannot be tracked or individual cell histories recorded.

Here we describe a methodology to quantify dynamic changes in a proliferating cell population and to deconvolve fractional cell fates over time (Fig. 1). We label cells with a nuclear fluorescent protein, carry out extended time-resolved automated microscopy (typically for ~96 h and as long as 10 days), segment nuclei in image stacks to directly count cells and track single cells for determining cell lifespans and fates (Methods). The large image datasets require analysis with computational tools and newly developed mathematical models (Supplementary Notes 1 and 2) that utilize information at both the population and single-cell levels and integrate them into a model of cell proliferation dynamics (Fig. 1). Briefly, we fit cell count data obtained from time-lapse automated microscopy with a novel Quiescence-Growth model, which incorporates rates of cell division, entry into quiescence and death. The model is constrained with experimental rates derived from tracked single cells and its output is visualized in Fractional Proliferation graphs. These graphs resolve the dynamic change in total cell numbers into fractions of dividing and quiescent (non-dividing) cells.


RESULTS
Quiescence-Growth model of cell population dynamics

We used our approach to monitor the proliferation dynamics of drug-treated cells. Cell counts extracted from the image data show that vehicle-and cycloheximide-treated PC9 cells exhibit linear proliferation (in log scale). In contrast, treatment with erlotinib, lapatinib, PLX-4720 or doxorubicin resulted in nonlinear effects on proliferation of PC9 cells (Fig. 2a). We observed similar effects, in other cell lines (Supplementary Fig. 1). We analyzed the expression of a marker of S or G2 phase (mAG-geminin)5 in drug-treated cells: erlotinib, lapatinib and PLX-4720 induced arrest in G1 (or G0) and doxorubicin in S or G2 (Fig. 2b), suggesting a correlation between cell cycle arrest and nonlinear proliferation. We also observed this correlation in erlotinib-treated PC9 cells analyzed using flow cytometry (Supplementary Fig. 2a, b) and in MCF10A cells deprived of serum (Supplementary Figure 2c). Nonlinear population dynamics are traditionally modeled using logistic or Gompertz equations 6,7, which have been applied to cultured cells as well as tumors 810. A major assumption of these models is the idea of a “carrying capacity” representing the maximum size a population can reach in a given environment. There is no clear relationship between carrying capacity and the biological processes that affect cell population size, such as cell division, quiescence and death.

To explain the nonlinear proliferation dynamics using biologically relevant parameters, we constructed a “Quiescence-Growth” model that includes two compartments, a dividing and a non-dividing (which we term quiescent) population, with death occurring from both compartments (Fig. 2c, Methods, and Supplementary Notes 1,2). The equation derived from this model incorporates three parameters (Fig. 2c), rates of division (d), death (a) and entry into quiescence (q). The term “quiescence” in the Quiescence-Growth model reflects arrest at any position in the cell cycle. We applied the Quiescence-Growth model to fit the cell count data obtained from drug-treated PC9 cells (Fig 2a). PC9 response to CHX was linear and best explained by decreasing the rate of division (d). In contrast, proliferation in response to erlotinib, lapatinib and PLX-4720 is nonlinear and can only be explained by varying the q parameter, rate of entry into quiescence (Fig. 2a). Varying the q parameter also best fits the response to doxorubicin (Fig. 2a), known to cause arrest in G2 (Fig. 2b).

Data-derived Quiescence-Growth model parameters

Good fits of the cell count data with the Quiescence-Growth model may be achieved by different combinations of parameter values (Supplementary Fig. 3, Supplementary Table 1 and Supplementary Software 1), particularly death and division rates, which oppose each other. We therefore constrained the model with experimentally measured rates from single cell tracking of time-lapse images.

We first tracked single cells across time lapse image stacks to quantify observed cell lifespans, which are defined as the time between an initial mitotic event and: 1) a death event, 2) another mitotic event (defining an intermitotic time; IMT), or 3) the end of the experiment (EoE) (Fig. 3a, Methods and Supplementary Video 1). Cell lifespans demarcated by an initial mitotic event and the EoE may belong to either dividing or non-dividing (quiescent) fractions.

For death rates, we identified death events by shrinkage/disintegration of nuclei (Supplementary Fig. 4a), tallied these events and converted them directly to rates (Methods).

For division rates, we first determined the variability of IMT in the population by examining IMT distributions, which are non-Gaussian by visual inspection (Fig. 3b, c and Supplementary Figs. 5,6), fail several statistical tests of normality (Supplementary Tables 2,3) and cause bias in the estimation of division rate if a simple mean value is used. We searched for model distributions that adequately capture the variability (Fig. 3b). Among several possible models we examined (Supplementary Table 2), an Exponentially Modified Gaussian (EMG) model 11 fits well the observed IMT distributions from a wide range of cells and conditions (Fig. 3, Supplementary Figs. 5,6), while minimizing the number of parameters to three. An additional benefit is that EMG model parameters are mathematically and biologically separable, i.e. their values are differentially affected by drugs (Supplementary Fig. 6).

We adapted a method of calculating division rates in bacterial cultures 12 to utilize all EMG parameters fit to the observed IMT distribution (Supplementary Notes 1,2). This method accounts for the dispersion of individual IMT, especially the slowly dividing cells in rightward skewed tails. In addition, this calculation of the division rate takes into account the age structure of an asynchronously dividing cell population under the assumption of exponential growth.

For quiescence rates, a crucial question is whether cells with lifespans demarcated by an initial mitotic event and the EoE are quiescent or whether they would have divided if the experiment had continued (Fig. 3a). We estimated the probability that a cell is quiescent by implementing a Likelihood of Quiescence Model based on a statistical survival formula 13 (Methods and Supplementary Notes 1,2). The older the undivided cell relative to the distribution, the lower the likelihood that it would have divided after the EoE, and the higher the likelihood that it is quiescent. The rate of entry into quiescence is then calculated using the fraction of quiescent cells in the population and is relative to the calculated rate of division. Importantly, without properly accounting for the slowly dividing cells in the population with a well-fit model (such as an EMG), the rate of entry into quiescence could be significantly overestimated.

We experimentally derived the three Quiescence-Growth model parameters from single cell tracking data of CA1d cells (Fig. 3c) treated with erlotinib (1 and 8 μM). Compared to control (DMSO), a large number of cells reached the EoE without dividing and were estimated to be quiescent by the Likelihood of Quiescence model (Table 1). The IMT distribution of cells treated with 8 μM erlotinib was flattened and elongated rightward (Fig. 3c). Rates calculated from the single-cell tracking data, EMG parameters and estimated quiescence fraction confirmed a decreased rate of division, a slight increase in death rate and a 10-fold increase in the rate of entry into quiescence at either erlotinib concentration. These experimentally-measured rates were in agreement with the rates estimated from fitting the Quiescence-Growth model to the cell count data (Supplementary Fig. 3 and Supplementary Table 1), indicating that the Quiescence-Growth model provides an accurate description of cell proliferation and validating its prediction that the nonlinear proliferation response to erlotinib (Fig. 2a) is explained by alteration of the rate of entry into quiescence (Fig. 3c). We observed that very few cells die during the experiment and cells considered quiescent did not exhibit preapoptotic nuclear condensation (not shown).

Fractional Proliferation graphs

We then used the parameters obtained from single-cell tracking data (Supplementary Video 2) in the Quiescence-Growth model to produce graphs of the total population deconvolved into fractions of dividing and quiescent cells (Fig. 3d and Supplementary Software 2). At ~24 h the fractions of cells in the two compartments are equal (Fig. 3d); in contrast, at 80 h the quiescent fraction is 82% and the dividing fraction is 18%.

Drug response of oncogene-addicted cells

We applied this approach to investigate the dynamic response of PC9 cells to erlotinib (Fig. 4 and Supplementary Fig. 1c). PC9 cells harbor activating mutations of the Epidermal Growth Factor Receptor (EGFR), are “oncogene-addicted,” and hypersensitive (GI50 <20 nM) to EGFR tyrosine kinase inhibitors (TKI). These cells are representative of EGFR-mutated tumor cells in lung cancer patients that respond favorably to treatment with erlotinib (or other EGFR-specific TKI). The prevalent view is that PC9 cells have a strong apoptotic response to EGFR TKI 14,15 but, to our knowledge, the relative contributions of cell death versus quiescence in the PC9 response to erlotinib remain uncharacterized.

From cell count data derived from image stacks, PC9 cells exhibit exponential proliferation with a rate = 0.0183 (doubling time = 38 h) in medium containing vehicle (DMSO); however, single-cell tracking data show that this value underestimates the rate of division and does not consider the contribution of quiescence and death (Fig. 4d). Proliferation decreases in a nonlinear fashion upon erlotinib treatment (Fig. 4a) and cells accumulate at the EoE (Fig. 3a and 4b). From single-cell tracking data, 53% and 73% of cells with initial mitotic events within the first 25 h entered quiescence in response to 15.6 nM and 62.5 nM erlotinib, respectively, compared to 4% in the control cells. The IMT distributions showed a rightward skew in response to erlotinib characterized by an increased k parameter value, while all other parameter values remain essentially unchanged (Fig. 4c). We calculated the rates of division, death and entry into quiescence from the single-cell tracking data and produced Fractional Proliferation graphs (Fig. 4d). As for CA1d cells, using these rates in the Quiescence-Growth model correctly predicted the experimentally determined cell counts (Fig. 4d). The rate of cell death increased from 0.0015 in DMSO to 0.0048 in 62.5 nM erlotinib, but this cannot by itself explain the decreased proliferation seen after 30 h. The q parameter increased >40-fold from DMSO to 62.5 nM erlotinib-treated cells (Fig. 4d). These results indicate that in “oncogene-addicted” PC9 cells, the antiproliferative response to erlotinib is primarily due to increased entry into a non-dividing or quiescent state, and not to apoptosis as commonly assumed.

We verified that erlotinib treatment increased the fraction of quiescent PC9 cells by flow cytometry on Ki-67 (Supplementary Fig. 7) and p27-immunostained cells (not shown). We further confirmed these observations in another oncogene-addicted model cell line (A375, representing B-Raf V600E-mutated melanoma) (Supplementary Note 3 and Supplementary Figs. 7,8). The non-dividing quiescent state we observe is not equivalent to a preapoptotic state; after 96 h in 1 μM erlotinib, a drug washout experiment showed that proliferation resumed to pretreatment levels (Supplementary Video 1 and Supplementary Fig. 9). Also, a large fraction of cells remained viable (Ki-67-positive) after more than 90 h of erlotinib treatment (Supplementary Fig. 7). Determining the eventual fate of these cells (division, death or extended quiescence) would require longer duration studies.

Application to primary cells

We assessed the performance of our approach on primary cells using a baculovirus-based transduction of a genetically-encoded fluorescent nuclear probe that is easy to use and commercially available (CellLight® Nucleus, Invitrogen). Labeled nuclei were detectable within 12 h, persisted for 5–6 days, and were amenable to automated counting and manual tracking. As an example, the proliferation of primary human squamous cell carcinoma cells is shown in Supplementary Figure 10 and Supplementary Video 3. It is worth noting that an EMG distribution also describes the distribution of IMT in primary cultured cells. Additional experimentation will be required to determine whether this distribution is applicable to all cell types.


DISCUSSION

The Fractional Proliferation method provides quantitative insight into cell proliferation in response to perturbations and permits deconvolution of the relative contribution of multiple cell fates to cell population dynamics. Using this approach, it is possible to capture the behavior of minor subpopulations of cells (e.g. stem or progenitor cells, clonal variants, drug resistant phenotypes) within perturbed populations. The importance of these subpopulations is increasingly appreciated, and methods to study them should be broadly applicable.

The integrated mathematical models we developed describe the emergence of population behavior from experimental data measuring single cell fates. We note that these models are not definitive; for example, they will need to be adjusted to more accurately accommodate age structure in asynchronously dividing populations; to apply them to stem cells, a rate of entry into one or more differentiated state would have to be included. However, image data sets such as those described here can be analyzed by future models. Furthermore, the approach could be enhanced by the incorporation of other read-outs, such as immunofluorescence on molecular markers or live cell reporters of molecular activity.

Time-lapse video microscopy has been used for decades, but its high-throughput implementation has only recently been enabled by technological advances in automated microscopy instrumentation coupled to computation. Methods to automate the extraction of mitotic phases and duration from time-lapse movies were recently developed 1618. These approaches mainly focus on extracting image features that correlate with mitotic events.

The labeling approach we employed involves the expression of a fluorescent nuclear protein (H2BmRFP) introduced by recombinant viral particles. Although other labeling techniques are feasible, and even phase contrast imaging can provide data amenable to this approach, the use of H2BmRFP overcomes limitations due to photobleaching, toxicity, low signal-to-noise ratio in images, and fluctuating subcellular localization. In our hands, nuclear dyes (e.g. Hoechst 33342) demonstrate significant toxicity and cytoplasmic and membrane dyes (e.g. carboxyfluorescein diacetate succinimidyl ester and carbocyanine-based dyes) are prone to photobleaching, tend to aggregate in vesicles and pose challenges to identification of individual cells. The development of nontoxic fluorescent nuclear dyes would enable rapid cell labeling without genetic manipulation.

Observation of cells for several days permits study of the cell response to perturbations as they approach population steady state, which usually happens over many days of treatment. Our discovery that a major component of the PC9 response to TKI erlotinib is entry into a quiescent state challenges current views that the major response to erlotinib in oncogene-addicted cells is cell death 15, but is consistent with previous reports that erlotinib and other EGFR TKI induce cell cycle arrest in G1 in many cancer cell lines 1921.

Erlotinib is of special interest because it is currently used clinically as adjuvant or first-line therapy in subsets of lung cancer patients with sensitizing EGFR mutationsfor reviews on the subject see references 22,23. Questions remain as to the basis for erlotinib hypersensitivity (“oncogene addiction”) in EGFR-mutated cancer cells, and the inevitable rise of resistance 24. As we have demonstrated, application of the Fractional Proliferation method should help answer these questions.

Describing IMT distributions with an EMG model provides value beyond calculating the rates of division and quiescence in the Quiescence-Growth model. The EMG model is mathematically separable into two components, Exponential and Gaussian. Intriguingly, erlotinib primarily affects the exponential (k) while CHX primarily affects the mean of the Gaussian (μ) (Supplementary Fig. 6). The possibility that these separable components reflect distinct biochemical mechanisms is currently under investigation. The EMG can be interpreted biologically as a combination of a Gaussian process (several random variables that are additive such as protein accumulation) and an exponential process (e.g. a checkpoint process such as the G1/S transition with a chance of passing at each check). Other descriptive distributions offer no simple interpretation biologically 25. The EMG may be more common in time-dependent cellular processes due to the mathematical properties of the chemical master equation which can have an exponential distribution of halting times after a well mixing time has passed 26,27.

It has long been known that antiproliferative responses are not limited to a single fate but the lack of methods to deconvolve behaviors has forced assumptions of linearity. The Fractional Proliferation method provides a means to accommodate nonlinear antiproliferative responses and to separate the underlying cellular fates that shape the population-level response.


METHODS
Software

We provide: i) an interactive Mathematica CDF Player file to explore the effects of altering Quiescence-Growth model parameter values on cell proliferation plots (Supplementary Software 1); ii) an ImageJ28,29 macro for automated cell counting (enumeration of nuclei) (Supplementary Note 2); and iii) fracprolif, an extension of the freely available statistical software package R30 (http://www.R-project.org) that incorporates all of the code (Supplementary Software 2) used to analyze single-cell tracking and cell count proliferation data and generate Fractional Proliferation Graphs (described in Supplementary Note 2).

Cell culture and labeling

The following cell lines were used: MCF10A, MCF10A-CA1d (abbreviated as CA1d), SQ20B, and PC9. MCF10A and CA1d cells were cultured in DMEM/F12 medium containing 10% equine serum, 5% fetal bovine serum, 20 ng/ml epidermal growth factor, 10 μg/ml insulin, 500 ng/ml hydrocortisone, and 100 ng/ml cholera toxin. SQ20B cells were cultured in DMEM containing 20% fetal bovine serum and 400 ng/ml hydrocortisone, and PC9 cells were cultured in RPMI 1640 medium containing 10% fetal bovine serum. Primary human tumor-derived cells were obtained from a patient with squamous cell carcinoma of the tongue and grown continuously in culture in keratinocyte (serum-free) growth medium (Invitrogen). Cell lines were engineered to express the histone H2B/monomeric red fluorescent protein (H2B-mRFP) fusion protein using lentivirus-mediated transduction as previously described 31. Brightly fluorescent single-cell clones were expanded and compared to parental populations using traditional proliferation measurements to ensure they were representative of the initial population. Primary cells were induced to express H2B-mRFP using recombinant baculoviral particles (CellLight® Nucleus, Invitrogen) at 20 particles per cell, according to the manufacturer’s instructions.

Imaging with extended temporally-resolved automated microscopy (ETRAM)

Imaging was performed on 96-well plates (BD cat# 353219) using a BD Pathway 855 with a 20X (0.75 NA) objective in a CO2- and temperature-controlled environment. Images were acquired every 6–30 min using BD Attovision 1.6.2 software with the instrument in confocal mode (spinning disk). Nine adjacent images were captured at 0.4 s exposure and 2X2 binning to comprise a single 3X3 montage (approximately 800 μm2) from each of approximately 40–60 wells per experiment, and images were acquired for at least 72 h. Example montaged images are shown in Supplementary Figure 1a. Assays minimally included duplicate or triplicate wells and complete experiments were performed at least twice. Cells were seeded at 2,500–5,000 cells per well and allowed to grow overnight yielding approximately 200–600 cells at the onset of imaging. Cells were imaged until confluence (approximately 3,000 cells in an image) was achieved in control wells. Erlotinib concentration-response curves on PC9 cells were performed more than five times.

Enumerating nuclei

Nuclei were counted from ETRAM-generated image stacks sampled at approximately 1 h intervals. Images were imported into the freely available ImageJ (http://rsb.info.nih.gov/ij/) program and subjected to a macro optimized for images obtained from a BD Pathway 855. The macro: i) corrects for uneven illumination using a 50 pixel diameter rolling ball filter; ii) converts images to binary using a predefined threshold intensity value; iii) segments individual nuclei using a watershed algorithm; and iv) quantifies objects within specified range of areas and circularity. Example images and resultant cell population plots are shown in Supplementary Figure 1a, b. The ImageJ macro is provided in Supplementary Note 2, and an example image sequence showing nuclei enumeration by ImageJ is shown in Supplementary Figure 4.

Quantifying rates of cell death

An example of a region of a manually tracked ETRAM-generated image stack is shown in Supplementary Video 1. Cell death was identified by nuclear shrinkage to less than 50% of average nuclear area and subsequent nuclear dispersion or detachment (no longer detectable) (Supplementary Fig. 4). Rates of cell death were determined by dividing the number of events (cell deaths) detected across all frames by the total cell observation time (number of cell nuclei in each frame multiplied by the time interval between frames); this rate does not require individual cells to be tracked over complete lifespans.

Quantifying intermitotic times and rate of division

For most cell lines, ETRAM was performed at 12 min intervals to minimize light exposure and phototoxicity, and for highly motile cells at 6 min intervals to eliminate bias against faster moving cells.

Intermitotic times were measured by manually tracking individual nuclei through the series of ETRAM-generated images, identifying mitotic events (metaphase chromosomes), and determining the number of frames between mitotic events for each individual cell lifespan. Lifespans of cell nuclei that divided once during the experiment but reached the EoE were also determined. 100 nuclei were tracked for each condition, and more if too many cells reached the EoE without having divided, since a minimum of 50 individual IMT represented a distribution.

Model of proliferation kinetics (Quiescence-Growth Model)

The Quiescence-Growth model is formulated as a pair of coupled ordinary differential equations where x represents the dividing cell number, y the quiescent cell number and x + y equals the total cell population. The rates are described by the three parameters: d for division or birth rate, q for quiescence rate and a for death rate.

[Formula ID: FD1]
x′=(x-q-a)xy′=qx-ay
which has an analytical solution of the following form:
[Formula ID: FD2]
x(t;x0,d,q,a)=x0e(d-q-a)ty0+x0(q(d-q)(e(d-q)t-1)e-at)

On a log scale the model is nonlinear due to the quiescence compartment (y). Also of note, when d=q the following solution applies.

[Formula ID: FD3]
limd→qy(t;x0,y0,d,q,a)=e-at(x0dt+y0)

When d is greater than or equal to q + a, the following asymptotic behavior is observed:

[Formula ID: FD4]
x+y~dd-qx0e(d-q-a)t,d>(q+a)

This shows that the model approaches exponential growth as long as the division rate exceeds the sum of the rates of death and entry into quiescence and mathematically proves that the observation of exponential proliferation of a population does not exclude the possibility of quiescence and death. Thus, only in the absence of any death or quiescence does the rate of cell division reflect the rate of proliferation of the population. On the other hand, note that the population could grow exponentially even if nearly half of the cells enter quiescence or die. Furthermore, entry into the quiescent compartment (with rate q) provides the only mechanism by which nonlinear proliferation curves can be achieved. An alternative of the model has been produced in which a different rate of death from each compartment is provided (Supplementary Note 1). However, since measuring these different rates from the single-cell data with sufficient statistical power is not yet feasible, the alternative model reverts to the form used in this method.

Statistical methods

To determine whether a model could describe the observed data with sufficient statistical accuracy, a one-sided Kolmogorov-Smirnoff test was performed with a two-sided test being the null hypothesis. This test determines the probability that the data and the model are sampled from the same distribution, and a P value of less than 0.05 was assumed to be statistically significant evidence that the data and the model represent different distributions.

The Shapiro-Wilk test was used to test for normality. The likelihood of one of two specific models (e.g. a linear model compared to the Quiescence-Growth model) correctly fitting the data was determined using Akaike’s Information Criteria <http://en.wikipedia.org/wiki/Akaike_information_criterion>.



Notes

FN2AUTHOR CONTRIBUTIONS

D.R.T. conceived of the approach, D.R.T. and P.L.F. cultured, treated and imaged cells, D.R.T. analyzed images; S.P.G. and D.R.T. developed the mathematical models, D.R.T. and S.P.G. fit model parameters to data; D.R.T. and V.Q. co-wrote the paper.

We would like to thank J. Hao and A. Udyavar for technical assistance, E. Pham and R. Feroze for assistance in manual image analysis, G. Ostheimer for assistance with flow cytometry and A. Weaver, G. Webb, B. Rexer, K. Dahlman, W. Yarbrough and L. Estrada for reviewing the manuscript and providing stimulating discussions. We acknowledge the generous gift of primary squamous cell carcinoma and SQ20B cells from W. Yarbrough and PC9 cells from W. Pao (both at Vanderbilt University School of Medicine). We also acknowledge A. Miyawaki (RIKEN Brain Science Institute) for the mAG-geminin plasmid and Gideon Bollag (Plexxikon) for the generous gift of PLX-4720. This work was supported by the US National Institutes of Health/National Cancer Institute Integrative Cancer Biology Program (5U54 CA113007-07). Flow cytometry experiments were performed in the Vanderbilt Medical Center Flow Cytometry Shared Resource, which is supported by the Vanderbilt Ingram Cancer Center (P30 CA68485) and the Vanderbilt Digestive Disease Research Center (DK058404). In addition, the project described was partially supported by the National Center for Research Resources (UL1 RR024975-01) and is now at the National Center for Advancing Translational Sciences (2 UL1 TR000445-06).


LITERATURE CITED
1. Hughes M,et al. Sittampalam GS,et al.Early Drug Discovery and Development Guidelines: For Academic Researchers, Collaborators, and Start-up CompaniesAssay Guidance ManualBethesda (MD)Year: 2004
2. Terry NH,White RA. Flow cytometry after bromodeoxyuridine labeling to measure S and G2+M phase durations plus doubling times in vitro and in vivoNat Protoc1859869Year: 200617406318
3. van Genderen H,et al. In vitro measurement of cell death with the annexin A5 affinity assayNat Protoc1363367Year: 200617406257
4. Hawkins ED,et al. Measuring lymphocyte proliferation, survival and differentiation using CFSE time-series dataNat Protoc220572067Year: 200717853861
5. Sakaue-Sawano A,et al. Visualizing spatiotemporal dynamics of multicellular cell-cycle progressionCell132487498Year: 200818267078
6. d’Onofrio A,Fasano A,Monechi B. A generalization of Gompertz law compatible with the Gyllenberg-Webb theory for tumour growthMath Biosci2304554Year: 201121232543
7. Gyllenberg M,Webb GF. Age-size structure in populations with quiescenceMath Biosci866795Year: 1987
8. Florian JA Jr,Eiseman JL,Parker RS. Accounting for quiescent cells in tumour growth and cancer treatmentSyst Biol (Stevenage)152185192Year: 200516986259
9. Kozusko F,Bajzer Z. Combining Gompertzian growth and cell population dynamicsMath Biosci185153167Year: 200312941534
10. Turner ME Jr,Bradley EL,Kirk KA,Pruitt KM. A theory of growthMath Biosci29367373Year: 1976
11. Grushka E. Characterization of exponentially modified Gaussian peaks in chromatographyAnal Chem4417331738Year: 197222324584
12. Powell EO. Growth rate and generation time of bacteria, with special reference to continuous cultureJ Gen Microbiol15492511Year: 195613385433
13. McLachlan GJ,Peel D. Finite Mixture ModelsJohn Wiley & Sons, IncYear: 2000
14. Gong Y,et al. Induction of BIM is essential for apoptosis triggered by EGFR kinase inhibitors in mutant EGFR-dependent lung adenocarcinomasPLoS Med4e294Year: 200717927446
15. Sharma SV,et al. A chromatin-mediated reversible drug-tolerant state in cancer cell subpopulationsCell1416980Year: 201020371346
16. Harder N,et al. Automatic analysis of dividing cells in live cell movies to detect mitotic delays and correlate phenotypes in timeGenome Res1921132124Year: 200919797680
17. Held M,et al. CellCognition: time-resolved phenotype annotation in high-throughput live cell imagingNat Methods7747754Year: 201020693996
18. Sigoillot FD,et al. A time-series method for automated measurement of changes in mitotic and interphase duration from time-lapse moviesPLoS One6e25511Year: 201121966537
19. Bos M,et al. PD153035, a tyrosine kinase inhibitor, prevents epidermal growth factor receptor activation and inhibits growth of cancer cells in a receptor number-dependent mannerClin Cancer Res320992106Year: 19979815602
20. Ciardiello F,Tortora G. A novel approach in the treatment of cancer: targeting the epidermal growth factor receptorClin Cancer Res729582970Year: 200111595683
21. Ling YH,et al. Erlotinib, an effective epidermal growth factor receptor tyrosine kinase inhibitor, induces p27KIP1 up-regulation and nuclear translocation in association with cell growth inhibition and G1/S phase arrest in human non-small-cell lung cancer cell linesMol Pharmacol72248258Year: 200717456787
22. Grunwald V,Hidalgo M. Developing inhibitors of the epidermal growth factor receptor for cancer treatmentJ Natl Cancer Inst95851867Year: 200312813169
23. Pao W,Chmielecki J. Rational, biologically based treatment of EGFR-mutant non-small-cell lung cancerNat Rev Cancer10760774Year: 201020966921
24. Riely GJ,et al. Clinical course of patients with non-small cell lung cancer and epidermal growth factor receptor exon 19 and exon 21 mutations treated with gefitinib or erlotinibClin Cancer Res12839844Year: 200616467097
25. Golubev A. Exponentially modified Gaussian (EMG) relevance to distributions related to cell proliferation and differentiationJ Theor Biol262257266Year: 201019825376
26. Sbano L,Kirkilionis M. Multiscale analysis of reaction networksTheory Biosci127107123Year: 200818446398
27. Chichagov VV. Asymptotic behavior of the first arrival time in Markovian random walksJ Math Sci29442948Year: 1987
28. Abramoff MD,Magalhaes PJ,Ram SJ. Image Processing with ImageJBiophotonics International113642Year: 2004
29. Rasband WS. ImageJU. S. National Institutes of HealthBethesda, MarylandYear: 1997–2012
30. R Development Core TeamR: A language and environment for statistical computingR Foundation for Statistical ComputingVienna, AustriaYear: 2010
31. Quaranta V,et al. Trait variability of cancer cells quantified by high-content automated microscopy of single cellsMethods Enzymol4672357Year: 200919897088

Article Categories:
  • Article


Previous Document:  Insertion sequences shared by Bordetella species and implications for the biological diagnosis of pe...
Next Document:  Accurate gene synthesis with tag-directed retrieval of sequence-verified DNA molecules.