Fractional proliferation: a method to deconvolve cell population dynamics from singlecell data.  
Jump to Full Text  
MedLine Citation:

PMID: 22886092 Owner: NLM Status: MEDLINE 
Abstract/OtherAbstract:

We present an integrated method that uses extended timelapse automated imaging to quantify the dynamics of cell proliferation. Cell counts are fit with a quiescencegrowth 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 'oncogeneaddicted' 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 oncogenetargeted 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: 20120812 
Journal Detail:

Title: Nature methods Volume: 9 ISSN: 15487105 ISO Abbreviation: Nat. Methods Publication Date: 2012 Sep 
Date Detail:

Created Date: 20120904 Completed Date: 20121113 Revised Date: 20130712 
Medline Journal Info:

Nlm Unique ID: 101215604 Medline TA: Nat Methods Country: United States 
Other Details:

Languages: eng Pagination: 9238 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 ProteinTyrosine Kinases / antagonists & inhibitors SingleCell Analysis / methods* StructureActivity Relationship Time Factors Tumor Cells, Cultured 
Grant Support  
ID/Acronym/Agency:

2 UL1 TR00044506/TR/NCATS NIH HHS; 5U54 CA11300707/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 RR02497501/RR/NCRR NIH HHS; UL1 TR000445/TR/NCATS NIH HHS 
Chemical  
Reg. No./Substance:

0/Protein Kinase Inhibitors; EC 2.7.10.1/ProteinTyrosine Kinases 
Comments/Corrections 
Full Text  
Journal Information Journal ID (nlmjournalid): 101215604 Journal ID (pubmedjrid): 32338 Journal ID (nlmta): Nat Methods Journal ID (isoabbrev): Nat. Methods ISSN: 15487091 ISSN: 15487105 
Article Information Download PDF License: nihmssubmitted publication date: Day: 15 Month: 8 Year: 2012 Electronic publication date: Day: 12 Month: 8 Year: 2012 Print publication date: Month: 9 Year: 2012 pmcrelease 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 singlecell 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 
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 counts^{1}. 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 singlecell 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 timeresolved 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 singlecell levels and integrate them into a model of cell proliferation dynamics (Fig. 1). Briefly, we fit cell count data obtained from timelapse automated microscopy with a novel QuiescenceGrowth 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 (nondividing) cells.
We used our approach to monitor the proliferation dynamics of drugtreated cells. Cell counts extracted from the image data show that vehicleand cycloheximidetreated PC9 cells exhibit linear proliferation (in log scale). In contrast, treatment with erlotinib, lapatinib, PLX4720 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 (mAGgeminin)^{5} in drugtreated cells: erlotinib, lapatinib and PLX4720 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 erlotinibtreated 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 ^{8–10}. 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 “QuiescenceGrowth” model that includes two compartments, a dividing and a nondividing (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 QuiescenceGrowth model reflects arrest at any position in the cell cycle. We applied the QuiescenceGrowth model to fit the cell count data obtained from drugtreated 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 PLX4720 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).
Good fits of the cell count data with the QuiescenceGrowth 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 timelapse 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 nondividing (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 nonGaussian 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 wellfit model (such as an EMG), the rate of entry into quiescence could be significantly overestimated.
We experimentally derived the three QuiescenceGrowth 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 singlecell tracking data, EMG parameters and estimated quiescence fraction confirmed a decreased rate of division, a slight increase in death rate and a 10fold increase in the rate of entry into quiescence at either erlotinib concentration. These experimentallymeasured rates were in agreement with the rates estimated from fitting the QuiescenceGrowth model to the cell count data (Supplementary Fig. 3 and Supplementary Table 1), indicating that the QuiescenceGrowth 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).
We then used the parameters obtained from singlecell tracking data (Supplementary Video 2) in the QuiescenceGrowth 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%.
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 “oncogeneaddicted,” and hypersensitive (GI50 <20 nM) to EGFR tyrosine kinase inhibitors (TKI). These cells are representative of EGFRmutated tumor cells in lung cancer patients that respond favorably to treatment with erlotinib (or other EGFRspecific 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, singlecell 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 singlecell 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 singlecell tracking data and produced Fractional Proliferation graphs (Fig. 4d). As for CA1d cells, using these rates in the QuiescenceGrowth 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 >40fold from DMSO to 62.5 nM erlotinibtreated cells (Fig. 4d). These results indicate that in “oncogeneaddicted” PC9 cells, the antiproliferative response to erlotinib is primarily due to increased entry into a nondividing 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 Ki67 (Supplementary Fig. 7) and p27immunostained cells (not shown). We further confirmed these observations in another oncogeneaddicted model cell line (A375, representing BRaf V600Emutated melanoma) (Supplementary Note 3 and Supplementary Figs. 7,8). The nondividing 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 (Ki67positive) 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.
We assessed the performance of our approach on primary cells using a baculovirusbased transduction of a geneticallyencoded 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.
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 readouts, such as immunofluorescence on molecular markers or live cell reporters of molecular activity.
Timelapse video microscopy has been used for decades, but its highthroughput 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 timelapse movies were recently developed ^{16–18}. 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 signaltonoise 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 carbocyaninebased 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 oncogeneaddicted 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 ^{19–21}.
Erlotinib is of special interest because it is currently used clinically as adjuvant or firstline 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 EGFRmutated 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 QuiescenceGrowth 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 timedependent 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 populationlevel response.
We provide: i) an interactive Mathematica CDF Player file to explore the effects of altering QuiescenceGrowth model parameter values on cell proliferation plots (Supplementary Software 1); ii) an ImageJ^{28,29} macro for automated cell counting (enumeration of nuclei) (Supplementary Note 2); and iii) fracprolif, an extension of the freely available statistical software package R^{30} (http://www.Rproject.org) that incorporates all of the code (Supplementary Software 2) used to analyze singlecell tracking and cell count proliferation data and generate Fractional Proliferation Graphs (described in Supplementary Note 2).
The following cell lines were used: MCF10A, MCF10ACA1d (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 tumorderived cells were obtained from a patient with squamous cell carcinoma of the tongue and grown continuously in culture in keratinocyte (serumfree) growth medium (Invitrogen). Cell lines were engineered to express the histone H2B/monomeric red fluorescent protein (H2BmRFP) fusion protein using lentivirusmediated transduction as previously described ^{31}. Brightly fluorescent singlecell 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 H2BmRFP using recombinant baculoviral particles (CellLight® Nucleus, Invitrogen) at 20 particles per cell, according to the manufacturer’s instructions.
Imaging was performed on 96well plates (BD cat# 353219) using a BD Pathway 855 with a 20X (0.75 NA) objective in a CO_{2} and temperaturecontrolled 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 μm^{2}) 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 concentrationresponse curves on PC9 cells were performed more than five times.
Nuclei were counted from ETRAMgenerated 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.
An example of a region of a manually tracked ETRAMgenerated 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.
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 ETRAMgenerated 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.
The QuiescenceGrowth 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.
x′=(xqa)xy′=qxay 
x(t;x0,d,q,a)=x0e(dqa)ty0+x0(q(dq)(e(dq)t1)eat) 
On a log scale the model is nonlinear due to the quiescence compartment (y). Also of note, when d=q the following solution applies.
limd→qy(t;x0,y0,d,q,a)=eat(x0dt+y0) 
When d is greater than or equal to q + a, the following asymptotic behavior is observed:
x+y~ddqx0e(dqa)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 singlecell data with sufficient statistical power is not yet feasible, the alternative model reverts to the form used in this method.
To determine whether a model could describe the observed data with sufficient statistical accuracy, a onesided KolmogorovSmirnoff test was performed with a twosided 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 ShapiroWilk test was used to test for normality. The likelihood of one of two specific models (e.g. a linear model compared to the QuiescenceGrowth model) correctly fitting the data was determined using Akaike’s Information Criteria <http://en.wikipedia.org/wiki/Akaike_information_criterion>.
2 Click here for additional data file (NIHMS397921supplement2.txt)
3 Click here for additional data file (NIHMS397921supplement3.txt)
Video 1 Click here for additional data file (NIHMS397921supplementVideo_1.mov)
Video 2 Click here for additional data file (NIHMS397921supplementVideo_2.mov)
Video 3 Click here for additional data file (NIHMS397921supplementVideo_3.mov)
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. cowrote 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 mAGgeminin plasmid and Gideon Bollag (Plexxikon) for the generous gift of PLX4720. This work was supported by the US National Institutes of Health/National Cancer Institute Integrative Cancer Biology Program (5U54 CA11300707). 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 RR02497501) and is now at the National Center for Advancing Translational Sciences (2 UL1 TR00044506).
LITERATURE CITED
1.  Hughes M,et al. Sittampalam GS,et al.Early Drug Discovery and Development Guidelines: For Academic Researchers, Collaborators, and Startup 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 timeseries dataNat Protoc220572067Year: 200717853861 
5.  SakaueSawano A,et al. Visualizing spatiotemporal dynamics of multicellular cellcycle progressionCell132487498Year: 200818267078 
6.  d’Onofrio A,Fasano A,Monechi B. A generalization of Gompertz law compatible with the GyllenbergWebb theory for tumour growthMath Biosci2304554Year: 201121232543 
7.  Gyllenberg M,Webb GF. Agesize 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 EGFRdependent lung adenocarcinomasPLoS Med4e294Year: 200717927446 
15.  Sharma SV,et al. A chromatinmediated reversible drugtolerant 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: timeresolved phenotype annotation in highthroughput live cell imagingNat Methods7747754Year: 201020693996 
18.  Sigoillot FD,et al. A timeseries method for automated measurement of changes in mitotic and interphase duration from timelapse 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 numberdependent 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 upregulation and nuclear translocation in association with cell growth inhibition and G1/S phase arrest in human nonsmallcell 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 EGFRmutant nonsmallcell lung cancerNat Rev Cancer10760774Year: 201020966921 
24.  Riely GJ,et al. Clinical course of patients with nonsmall 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 highcontent automated microscopy of single cellsMethods Enzymol4672357Year: 200919897088 
Article Categories:

Previous Document: Insertion sequences shared by Bordetella species and implications for the biological diagnosis of pe...
Next Document: Accurate gene synthesis with tagdirected retrieval of sequenceverified DNA molecules.