A transfer function approach to measuring cell inheritance
© Rees et al; licensee BioMed Central Ltd. 2011
Received: 23 August 2010
Accepted: 22 February 2011
Published: 22 February 2011
The inheritance of cellular material between parent and daughter cells during mitosis is highly influential in defining the properties of the cell and therefore the population lineage. This is of particular relevance when studying cell population evolution to assess the impact of a disease or the perturbation due to a drug treatment. The usual technique to investigate inheritance is to use time-lapse microscopy with an appropriate biological marker, however, this is time consuming and the number of inheritance events captured are too low to be statistically meaningful.
Here we demonstrate the use of a high throughput fluorescence measurement technique e.g. flow cytometry, to measure the fluorescence from quantum dot markers which can be used to target particular cellular sites. By relating, the fluorescence intensity measured between two time intervals to a transfer function we are able to deconvolve the inheritance of cellular material during mitosis. To demonstrate our methodology we use CdTe/ZnS quantum dots to measure the ratio of endosomes inherited by the two daughter cells during mitosis in the U2-OS, human osteoscarcoma cell line. The ratio determined is in excellent agreement with results obtained previously using a more complex and computational intensive bespoke stochastic model.
The use of a transfer function approach allows us to utilise high throughput measurement of large cell populations to derive statistically relevant measurements of the inheritance of cellular material. This approach can be used to measure the inheritance of organelles, proteins etc. and also particles introduced to cells for drug delivery.
The function of an organism is determined by the evolution of a cell population all descended from a single progenitor cell. The behaviour of individual cells is thus determined not only by their environment but also by the composition of cellular material inherited from the parent cell during mitosis. Therefore, the segregation of cellular material during mitosis is critical in determining the fate of the daughter cells. Previously, it was assumed that mitotic partitioning of cellular material was equal between daughter cells , however more recent studies have shown that certain cell components partition asymmetrically at mitosis [2–5] and there is a growing realisation that division asymmetry is a fundamental property of biological cells. For example, it has been demonstrated that proteins destined for proteasomal degradation within aggresomes are preferentially inherited by one daughter and it has been suggested that the generation of different daughter cells promotes long term lineage survival . Stem cell differentiation has also been linked to asymmetric inheritance of endosome function in daughter cells [2, 7].
Whilst division asymmetry is now clearly recognised, it has been difficult to accurately quantify as most studies rely on time-lapse microscopy, which can provide measurements on a limited number of cells. There is a need therefore for statistically relevant measures of inheritance quantifying the inter-generational transformation across a cell population. In previous studies, we have demonstrated that through flow cytometric analysis of colloidal quantum dot (QD) fluorescence the asymmetric inheritance of endosomes can be accurately quantified, based on measurements from sets of > 104 cells . This technique employed a stochastic cell model of all individual cells within a population to predict the generational redistribution of the QD-loaded endosome fluorescence over time. The model parameters were optimised, via an evolutionary algorithm, to maximise correlation between the fluorescent histograms generated numerically to that experimentally measured. This process predicts an asymmetric redistribution of QDs across the population  as opposed to the intuitive guess of symmetric partitioning. While this approach proved a robust method for analysing cellular inheritance, the algorithms employed to predict the evolvution of the in-silico population are numerically bespoke requiring an in-depth knowledge of the biological principles and processes involved during the complex cell-cycle.
The purpose of this paper is to outline a method for assessing the inheritance of cellular material using high throughput fluorescence measurements combined with simple and more accessible models derived from systems theory. The focus of systems biology is usually on understanding the complex interactions which occur between the components of a biological system e.g. the protein interaction networks driving a specific cellular function or response. By identifying the relevant components within the complex system a host of mathematical tools can be used to attempt to identify the nature of the mechanisms between them. Here we take inspiration from traditional systems engineering to take an alternative viewpoint that concentrates on the dynamics of the whole system (cellular population evolution) through relation of initial and final states of that system by a transfer function. A transfer function is simply a mathematical expression which relates the input (function) and the output (function) of a (non-)linear system. It must the formulated to include the physical (or biological) processes involved in the system which modifies the variable or quantity that is being measured. We use the transfer function approach to encapsulate and transform the stochastic population dynamics of the stochastic cell cycle model into a continuous form. This entails a significant simplification of the complex processes associated with the biological system processes, as they are now described purely phenomenologically through their influence on the evolving cellular state. However we will demonstrate that appropriately designed experiments coupled with simple mathematical models of the transfer function can accurately elucidate biological function and provide a generalised analysis method that delivers relevant information on the system without requiring a-priori knowledge of sub-cellular processes. Casting the problem in terms of a transfer function between an initial and final state-function opens up the possibility of using system analysis such as convolution theory to deconvolve cell inheritance information from flow measurements.
Results and Discussion
The aim of the analysis is to effectively determine this distribution i.e determine α and Δα for the mitosis events within a cell population. The technique does not measure the inheritance for a single cell, but rather a distribution of inheritance based on all cells within the population. The assumption that we are dealing with large cell numbers allows us to use the continuous Gaussian distribution to represent the random variation of the inheritance around the mean value.
where f t is the fraction of cells, having undergone mitosis during t and δ is the Dirac delta function.
We note that the implementation of the convolution integral in equation 2, requires a summation over the fluorescence intensity and therefore it is worth discussing how the intensity data is prepared to provide the input function N i (I) for the transfer function T(I,I',t). The high-throughput fluorescence imaging technique would need to measure the total intensity per cell, however, in most flow cytometers the fluorescence intensity of the cell is measured using a logarithmic amplifier in order to maximise the dynamic range of measurement and subsequently the output is in the form of a histogram with logarithmic bin spacing . This can easily be dealt with by either converting to a linear intensity histogram or applying the transfer function with an integral summed over variable intensity width appropriate for the logarithm scale. While flow cytometry is an ideal technique for implementation of this technique, we also note instruments such as high-throughput microscopes, imaging cytometers etc. can also be used to provide the average fluorescence intensity of each cell.
The first term on the right hand side of equation 3 describes the probability of a change in fluorescence intensity from I' to I given that f t is the fraction of cells having undergone mitosis in the time t when the second fluorescence intensity measurement is performed. The second term represents the cells which have not undergone mitosis and therefore their fluorescent intensity remains unchanged. Now provided we can determine f t we are in a position to apply equation 3 to the first measurement and optimise parameters α and Δα to obtain the best fit with the fluorescent intensity measurement at time t. There are several approaches we can take to determine f t . The simplest approach is to use the variable as a fitting parameter which gives three unknown variables within the transfer function to optimise. This does not cause any problems for any fitting algorithm employed as the cells which have not undergone mitosis have the highest intensity (while the daughter cells have a lower intensity) and therefore the shape of the intensity profile at high values is highly sensitive to the value of f t while remains insensitive to the values of α and Δα.
where I 0 and I t are the mean intensity of the cell population measured at time 0 and t respectively, obtained from the experimental flow histograms. This is of course only valid if we measure the fluorescence intensity from the same number of cells in each measurement. Using the data from the fluorescence histograms in figure 3 we obtain a value f t = 0.969 which is excellent agreement with the value determined using the fitting procedure. While this gives us confidence in the fitting procedure used, it is also possible to use this experimental value in the transfer function and therefore just fit for values of α and Δα.
We note most of the cells within the population would have undergone mitosis during the time duration chosen between the measurements. When designing this type of experiment it is important to maximise the number of mitosis events which helps increase the accuracy of the splitting parameters determined by the fitting procedure. However if the time interval between the two fluorescence intensity measurements becomes too long then the daughter cells themselves would start to undergo mitosis which adds another layer of complexity to the analysis. The measure of the fraction of cells that have undergone mitosis determined using equation 4 can provide a useful tool (if the intermitotic time of the cell line is unknown) in deciding the time interval between the two fluorescence measurements. The relative number of cells in each generation can be calculated for example using branching theory  however this requires a knowledge of the cell intermitotic time. The issue of the time between measurements also raises a further consideration for the planning of such an experiment. The analysis can only measure the inheritance of labelled cellular material and therefore the inheritance of any new cellular material which develops during the cell cycle cannot be elucidated upon.
As discussed earlier the method describe so far relies on the time period of the experiment to be sufficiently short so that no cell in the population divides twice. This is simply due to the way we have defined our transfer function (equation 3) in terms of cells that have split once and those that have not undergone mitosis. If an experiment necessitated a time period where several generations of cells might be observed then the analysis is perfectly valid providing the transfer function in equation3 included an extra term for each generation. For example, if 3 generations were to be observed an extra term would be included to describe the fraction of cells which had 'undergone two mitosis events'. The probability distribution associated with the cells would be a modified form of equation 1 which would now include 4 intensity partitioning Gaussian distributions for the 4 cells (originating from the one progenitor cell) in the third generation.
While for this study we have used fluorescent quantum dot markers to label the organelles, this protocol can be used with any other fluorescent dye providing a couple of criteria are met. Firstly the dye must fluoresce with intensity proportional to the quantity of organelles to be measured. Secondly the intensity of the dye must remain stable for the duration of the experiment. Also if the dye fluorescence intensity does deteriorate during the experiment by a quantifiable amount then this can be compensated for in the transfer function as a correction to the magnitude of the Gaussian distributions in equation 1.
In summary, we have presented a transfer function method for determining the inheritance of cellular material; the partitioning of endosomes loaded with colloidal quantum dots between daughter cells has been measured using flow cytometry. The use of a transfer function methodology provides a generic framework that is suitable to a wide range of cytological analyses. In particular it relates transformative states of the cell and in doing so provides quantitative information on the system evolution without the need for a-priori knowledge of the system components such as molecular networks or gene expression profiles. This whole-system approach is grounded in experimental observation, based as it is on the mathematical linking of two observational cell descriptor functions. In this work we compare whole cell parameters over a large population but the approach is equally valid for analysis of the evolution of a single cell captured using fluorescence microscopy. Cellular inheritance provides an ideal test system for the use of transfer functions as bifurcations at cell division give clearly marked events in the system history. Furthermore, the studies reported here can easily be extended to measure the inheritance of other cellular organelles [11, 12] or to describe the differing fate of daughter cells . More subtle system evolution within a single cell generation could also be tracked for instance cell-cycle progression through cyclin-based fluorescent reporters  or cellular metabolism via measurements of fluorescence lifetime .
QD loading and preparation of cells
Human osteosarcoma cell populations (U-2 OS cell line, ATCC HTB-96) were maintained in McCoy's 5a medium supplemented with 10% foetal calf serum (FCS), 1 mM glutamine, and antibiotics and incubated at 37°C in an atmosphere of 5% CO2 in air. The cells were loaded with commercially available targeted nanocrystals using the Qtracker® 705 (QTracker705) Cell Labelling Kit (4 nM) (Invitrogen (Q25061MP). The QTracker system uses a proprietary peptide targeting molecule, attached to the dot surface to enable receptor mediated endocytosis. Compared to untargeted QDots this produces a rapid cellular uptake (loading rate time of ~ 15 minutes compared to 4 hours for untargeted, carboxyl coated QDs) . The internalization of the QDots into discrete vesicles has been confirmed by confocal microscopy. The nature of the vesicle i.e. endosome or lysosome has not been established and will change over time, however this does not affect the division analysis as this considers numbers of labeled vesicles and is independent of their biological identity. The reagents in the Qtracker® 705 Cell Labelling Kit use a custom targeting peptide (9-arginine peptide) to deliver near-infrared-fluorescent nanocrystals into the cytoplasm of live cells via the endosomal pathway. Qtracker reagent A and B were premixed and then incubated for 5 mins at room temperature. 1 ml of fresh full growth media was added to the tube and vortexed for 30 seconds. This labelling solution was then added to each well of the cells and incubated for 1 hour at 37°C after which they were washed twice with fresh media and split into two flasks. To prepare for flow analysis samples were incubated in FACS Buffer (PBS/0.2% BSA/0.05% Sodium Azide) for 30 minutes before re-suspending in 200 μl of PBS and placed in a refrigerator at 4°C until later data acquisition on a flow cytometer. One sample was analysed using the flow cytometer while the second flask was incubated at 37°C for 19 hours and then analysed.
Cell images were acquired using an Imagestream100 Cell Analyser (Amnis Corporation, Seattle). A 488 nm wavelength laser was used to excite QD fluorescence which was collected via the 660-735 nm spectral detection channel. 104 cells were imaged for each sample and analyzed using the manufacturer's software. Gating of the acquired data to ensure focused images of viable cells reduced the analysed population to 5,000 cells. The system and peak image analysis algorithms were used to identify intensity clusters and calculate their number assuming a discrimination level of intensity at twice the intensity of the background. Spot counting accuracy was confirmed by manual verification of confocal images. The number of cells per ml within the gated population was also determined using the Imagestream100, which continuously runs speed calibration beads to provide absolute calibration of fluid flow rates.
This work was supported by the Engineering and Physical Sciences Research Council, U.K. under grant ref: EP/H008683/1, 'Nanoparticle Cytometrics' and the European Regional Development Funded, Swansea Centre for Nanohealth.
- Bergeland T, Widerberg J, Bakke O, Nordeng TW: Mitotic partitioning of endosomes and lysosomes. Current Biology 2001, 11: 644-651. 10.1016/S0960-9822(01)00177-4View ArticlePubMedGoogle Scholar
- Singhvi A, Garriga G: Asymmetric divisions, aggresomes and apoptosis. Trends in Cell Biol 2008, 19: 1-7. 10.1016/j.tcb.2008.10.004View ArticleGoogle Scholar
- Warren G: Membrane partitioning during cell division. Annu Rev Biochem 1993, 62: 323-348. 10.1146/annurev.bi.62.070193.001543View ArticlePubMedGoogle Scholar
- Dunster K, Hock Toh B, Sentry JW: Early endosomes, late endosomes, and lysosomes display distinct partitioning strategies of inheritance with similarities to Golgi-derived membranes. Eur Journ of Cell Biol 2002, 81: 117-124. 10.1078/0171-9335-00232View ArticleGoogle Scholar
- Errington RJ, Brown MR, Silvestre O, Njoh KL, Chappell SC, Khan I, Rees P, Wilks SP, Smith PJ, Summers HD: Single cell nanoparticle tracking to model cell cycle dynamics and compartmental inheritance. Cell Cycle 2010, 9: 1-10. 10.4161/cc.9.1.10246View ArticleGoogle Scholar
- Rujano MA, Bosveld F, Salomons FA, Dijk F, van Waarde MAWH, van der Want JJL, de Vos RAI, Brunt ER, Sibon OCM, Kampinga HH: Polarised asymmetric inheritance of accumulated protein damage in higher eukaryotes. PLoS Biology 2006, 4: 2325-2335. 10.1371/journal.pbio.0040417View ArticleGoogle Scholar
- Macara IG, Mili S: Polarity and differential inheritance - universal attributes of life? Cell 2008, 135: 801-812. 10.1016/j.cell.2008.11.006PubMed CentralView ArticlePubMedGoogle Scholar
- Brown MR, Summers HD, Rees P, Silvestre O, Njoh KL, Chappell SC, Khan I, Rees P, Wilks SP, Smith PJ, Errington RJ: Long-term time series analysis of quantum dot encoded cells by deconvolution of the auto fluorescence signal. Cytometry A 2010.Google Scholar
- Shapiro HM: Practical flow cytometry. John Wiley, Hoboken; 2003.View ArticleGoogle Scholar
- Kimmel M, Axelrod DE: Branching Processes in Biology. Springer-Verlag, New York; 2002.View ArticleGoogle Scholar
- Nunnari J, Walter P: Regulation of Organelle Biogenesis. Cell 1996, 84: 389-394. 10.1016/S0092-8674(00)81283-0View ArticlePubMedGoogle Scholar
- Warren G, Wickner W: Organelle Inheritance. Cell 1996, 84: 395-400. 10.1016/S0092-8674(00)81284-2View ArticlePubMedGoogle Scholar
- Darzynkiewicz Z, Gong J, Juan G, Ardelt B, Traganos F: Cytometry of Cyclin Proteins. Cytometry 1996, 25: 1-13. 10.1002/(SICI)1097-0320(19960901)25:1<1::AID-CYTO1>3.0.CO;2-NView ArticlePubMedGoogle Scholar
- Borchardt RT, Smith PL, Wilson G: Models for Assessing Drug Absorption and Metabolism. Springer; 1996.View ArticleGoogle Scholar
- Summers HD, Rees P, Holton MD, Brown MR, Chappell SC, Smith PJ, Errington RJ: Statistical analysis of nanoparticle dosing in a dynamic cellular system. Nature Nanotechnology 2011.Google Scholar