Studies on the variations in generation times of rat hepatoma cells in culture

The proliferation kinetics of a cultured hepatoma cell line, HTC, was studied by time-lapse photography and autoradiography. Variability of G1 phase, S phase and intermitotic times of sister cells were compared with data on the variability of G1 phase, G2 phase, S phase and intermitotic times of unrelated cells. Evidence is presented that the variation in times spent in intermitosis by a pair of sisters is determined by the variation of G2 phase durations, while the variation in intermitotic times of unrelated cells is mainly determined by the variation of G1 phase times.     

PROLIFERATION KINETICS OF AN EXPERIMENTAL ASCITES TUMOUR OF THE MOUSE

The cell cycles of an experimental ascitic tumour of the C3H mouse (NCTC 2472) were determined at various times after the intraperitoneal injection of 10⁶ cells. It was found that, contrary to results in solid NCTC 2472 tumours, obtained with the same NCTC cells, the duration of the cell cycle and its phases lengthened with the age of the tumour while the growth fraction remained relatively constant. G₁ was the first phase to lengthen, while later T_s and T_G2 increased also. The amount of DNA per cell was determined by cytospectrophotometry. This method provides data on the evolution during growth of the relative number of cells in each phase of the cell cycle.     

THE DISCRETE–TIME KINETIC MODEL ANALYSIS OF DNA CONTENT DISTRIBUTIONS IN EXPERIMENTAL TUMOUR CELLS

A method was developed to analyse and characterize FMF measurements of DNA content distribution, utilizing the discrete time kinetic (DTK) model for cell kinetics analysis. The DTK model determines the time sequence of the cell age distribution during the proliferation of a tumor cell population and simulates the distribution pattern of the DNA content of cells in each age compartment of the cell cycle. The cells in one age compartment are distributed and spread into several compartments of the DNA content distribution to allow for different rates of DNA synthesis and instrument dispersion effects. It is assumed that the DNA content of cells in each age compartment has a Gaussian distribution. Thus, for a given cell age distribution the DNA content distribution depends on two parameters of the cells in each age compartment: the average DNA content and its coefficient of variation. As the DTK model generates the best fit DNA content distribution to the FMF measurement data, it enables one to estimate specific values of these two parameters in each stage of the cell cycle and to determine the fraction of cells in each cycle phase. The method was utilized to fit FMF measurements of DNA content distributions and to analyse their relationship to the cell kinetic parameters, namely cell loss rate, cell cycle times and growth fraction of exponentially growing Chinese hamster ovary cells in vitro and, also, with a wide range of coefficients of variation, of the L1210 ascites tumour during the growth period.     

Estimating cell death in G(2)M using bivariate BrdUrd/DNA flow cytometry.

BACKGROUND: In an accompanying paper (Asmuth et al.) it was found necessary to include cell death explicitly to estimate parameters of cell proliferation. The use of bivariate flow cytometry to estimate the phase durations and the doubling times of cells labeled with thymidine analogues is well established. However, these methods of analysis do not consider the possibility of cell death. This report demonstrates that estimating cell death in G(2)/M is possible. METHODS: Mathematical models for the experimental quantities, the fraction of labeled undivided cells, the fraction of labeled divided cells, and the relative movement were developed. These models include the possibility that, of the cells with G(2)/M DNA content, only a certain fraction will divide, with the remainder dying after some time T(R). Simulation studies were conducted to test the possibility of using simple methods to estimate phase durations and cell death rates. RESULTS: Cell death alters the estimates of phase transit times in a rather complex manner that depends on the lifetime of the doomed cells. However, it is still possible to obtain estimates of the phase durations of cells in S and G(2)/M and the death rates of cells in G(2)/M. CONCLUSIONS: The methods presented herein provide a new way to characterize cell populations that includes cell death rates and common measurements of cell proliferation.     

PERTURBED CELL PROLIFERATION KINETICS OF THE EMT6/M/AC TUMOUR IN MICE TREATED WITH ANTI-MOUSE LYMPHOCYTE SERUM

The cell proliferation kinetics of the EMT6/M/AC tumour were determined at two volumes, 6-3 mm³ and 180 mm³, in mice treated with anti-mouse lymphocyte serum, AMLS. Comparison of the growth curve with that obtained in non-AMLS treated animals showed a marked increase in the growth rate at all volumes in the treated group. In contrast, the cell cycle time and the intermitotic phase times were not significantly different in the treated and untreated groups at comparable volumes. The increase in the growth rate in AMLS treated mice was obtained in spite of decreases in both the rate constant for cell production and the growth fraction, and was due to a marked decrease in the rate constant for cell loss.     

THE CELL CYCLE IN TUMOURS: AN EXAMINATION OF DATA GAINED BY THE TECHNIQUE OF LABELLED MITOSES

Published data from a large series of investigations by the technique of labelled mitoses have been analysed using a process of computer simulation. The results of the analyses of tumour data are discussed, first from the point of view of the quality of simulation that can be achieved and secondly as a review of the extent of our knowledge of the intermitotic time of cells in tumours. Particularly in primary tumors the data are often consistent with a broad distribution of intermitotic times within a tumour and the use of the term 'cell cycle time'may be misleading.     

CORRELATION BETWEEN CELL SIZE AND POSITION WITHIN THE DIVISION CYCLE IN SUSPENSION CULTURES OF CHANG LIVER CELLS

Chang liver cells from exponentially growing suspension cultures have been separated by sedimentation at unit gravity. Determinations of the protein content per cell showed that the fractionation procedure resulted in good separation of cells of different size. On the other hand, the DNA content of individual cells from the fractions, as determined cytofluorimetrically, indicated considerable heterogeneity in the size of cells from the same stage of the division cycle. On the basis of earlier results on intermitotic growth and the variation in the length of the cell cycle in homogeneous cell populations, a mathematical model has been constructed and tested using a computer program. The present results on the size distribution of cells from the different stages of the mitotic cycle are consistent with a regeneration of size heterogeneity in each cell generation, as a result of the dispersion of intermitotic times. The variation in cell cycle times may be related to a probabilistic event in the G₁ period. In the mathematical model it was necessary to include a mechanism by which the regeneration of abnormally large cells is prevented. The experimental data are compatible with a gradually increasing inhibition of growth in cells larger than a certain size (circa 400 pg protein per cell).     

II. USE OF THE METHOD TO MONITOR THE IN VIVO KINETICS OF CELL POPULATIONS PERTURBED BY HYDROXYUREA

In a previous report, we described a new method called FP_i analysis to analyze time sequences of DNA histograms taken from a perturbed population of cells. In this paper we utilize the method to analyze the in vivo kinetic response of bone marrow and of lung metastases of the B16 tumor to various chemotherapeutic agents. We show that the technique allows useful kinetic data to be obtained with minimal processing of the raw histograms, thus allowing fast analysis of the data. We also show that, in order to monitor the kinetic response of living tissues, it is essential to collect multiparameter distributions; to monitor only the one dimensional fluorescence histogram can give rise to misleading results. Using these multiparameter histograms, we are also able to monitor the growth fraction of the lung metastases during treatment, allowing discrimination between cell synchrony and cell recruitment from the resting compartment.     

Sensitivity of flow cytometric data to variations in cell cycle parameters

Abstract We investigated to what extent flow cytometric DNA histograms are informative of cell cycle parameters. We created a computer program to simulate cell cycle progression in a generic and flexible way. Various scenarios, characterized by different models and distributions of cell cycle phase transit times, have been analysed in order to obtain the percentages of cells in the different cell cycle phases during exponential growth and their time course after mitotic block.
Cell percentages during exponential growth were insensitive to intercell variability in phase transit times and thus can be employed to estimate the relative mean phase transit times, even in the presence of non-cycling cells. However, this information is ambiguous if re-entry of such cells into the cycling status is permitted. The stathmokinetic outline gives the mean phase transit times, but also provides information about the spread, but not the form, of the phase transit time distributions, being particularly sensitive to the spread of G1 phase duration. The stathmokinetic outline also helps distinguish between scenarios considering only cycling cells, those forecasting a fraction of definitively non-cycling cells and those admitting a Go status with first-order output kinetics.     

Sensitivity of flow cytometric data to variations in cell cycle parameters

We investigated to what extent flow cytometric DNA histograms are informative of cell cycle parameters. We created a computer program to simulate cell cycle progression in a generic and flexible way. Various scenarios, characterized by different models and distributions of cell cycle phase transit times, have been analysed in order to obtain the percentages of cells in the different cell cycle phases during exponential growth and their time course after mitotic block. Cell percentages during exponential growth were insensitive to intercell variability in phase transit times and thus can be employed to estimate the relative mean phase transit times, even in the presence of non-cycling cells. However, this information is ambiguous if re-entry of such cells into the cycling status is permitted. The stathmokinetic outline gives the mean phase transit times, but also provides information about the spread, but not the form, of the phase transit time distributions, being particularly sensitive to the spread of G1 phase duration. The stathmokinetic outline also helps distinguish between scenarios considering only cycling cells, those forecasting a fraction of definitively non-cycling cells and those admitting a G0 status with first-order output kinetics.     

On the measurement of cytokinetics by continuous labeling with bromodeoxyuridine with applications to chick wing buds.

The cytokinetic properties, specifically the phase-transit times, TG1, TS, and TG2+M, of chick wing bud cells were estimated using data obtained from continuous labeling of stage 20 embryos with bromodeoxyuridine (BrdUrd). The presence of BrdUrd was detected with monoclonal antibodies, and the amount of DNA in the cells was determined with propidium iodide. The fraction of cells in each cell cycle phase, the fraction of labeled cells, and the relative movement, a measure of the mean DNA content, of all labeled cells were evaluated using bivariate flow cytometry at successive times following introduction of the label. Equations are presented to describe the fraction of unlabeled cells in G2 + M, which gives a direct estimate of TG2+M; the fraction of all labeled cells, which can then be used to estimate TG1; and, finally, the relative movement, which provides an estimate of TS. Thus, the data measured in these experiments together provide estimates of the progression through the cell cycle of limb mesoderm cells.     

Chromatin-bound PCNA as S-phase marker in mononuclear blood cells of patients with acute lymphoblastic leukaemia or multiple myeloma

Objectives: Proliferating cell nuclear antigen (PCNA) has often been used as a marker to aid assessment of tumour growth fraction. This paper addresses the question of whether it can be used as an S‐phase marker, when the non‐chromatin‐bound form of the protein is removed by pepsin treatment. Materials and methods: Cytofluorometric measurements were carried out after immunofluorescence staining of PCNA and counterstaining of DNA. S‐phase fraction was determined with the help of windows on PCNA versus DNA scattergrams, or mathematically from DNA histograms. Results: S‐phase fractions obtained using the two methods correlated well, but did not always agree, exact discrepancies depending on the mathematical model used for histogram analysis. Conclusions: Determination of S‐phase fractions with the help of PCNA immunofluorescence staining is possible, and probably more reliable than calculation of S‐fractions from DNA histograms. It thus offers an alternative to assays involving BrdU labelling in vivo.     

Analysis of flow cytometric data of irradiated cell populations

Summary Flow cytometry (FCM) and autoradiography have been applied to determine changes in the cell kinetics of irradiated cells. Synchronized L-929 cells were irradiated with 10 Gy of X-rays when progressing from G₁-to S-phase of the cell cycle. In this study three methods to analyse DNA histograms were tested for applicability on FCM data obtained from cell populations blocked or retarded in the cycle: A) the Gaussian integral method, B) the peak-half-reflection method, and C) the rectangle method. Since histograms from synchronized cells are heavily distorted as compared to those obtained from exponentially growing cells and are quite similar to histograms from irradiated cells, they were used to test the suitability of the evaluation methods. Comparing the evaluated FCM data with the autoradiographic results from the same experimental series, the Gaussian integral method proved to be superior to the two other relatively simple approximation methods. The FCM histograms of irradiated cells were therefore analyzed only by the Gaussian integral method. It was shown that a considerable fraction of cells is still in the S-phase 25 h post irradiation, the DNA synthesis of which has ceased, as shown by autoradiography. This indicated that parallel measurements using FCM and autoradiography yield additional information on cell kinetic changes that cannot be obtained from applying one of the two methods used.     

PROLIFERATION KINETICS OF SHAY CHLOROLEUKAEMIA CELLS GROWN IN DIFFUSION CHAMBERS IN VIVO

The proliferation kinetics of Shay chloroleukaemia cells was studied by the labelled mitoses technique and some other methods. Detailed estimates of the phase durations of the cell cycle for proliferating cells were obtained using young, exponentially growing diffusion chamber cultures; estimates of kinetic parameters were also obtained for 'old' cultures and for local subcutaneous chloroma tumours. The results enabled determination of the growth fraction and the cell loss factor.
Cell loss was found to be the dominant factor determining the growth rate of chloroleukaemia cell populations. The mechanism of cell loss was cell death. Some implications of these findings are discussed.     

A comparison of mathematical methods for the analysis of DNA histograms obtained by flow cytometry

Abstract. Twelve methods for analysing FCM-histograms were compared using the same set of data. Some of the histograms that were analysed were simulated by computer and some were taken from experiments. Simulated data were generated assuming asynchronously growing cell populations and (i) measurement coefficients of variation ( CV ) from 2 to 16%; (ii) constant measurement CV or CV 's increasing from G₁ to G₂ phase, and (iii) varying fractions of cells in each phase. Simulated data were also generated assuming synchronous cell populations in which a block in early S phase was applied and released. DNA histograms were measured for L-929 cells at various times after mitotic selection. Labelling indices were also measured for these cells at the same time.
The fractions of cells in the G₁, S, and (G₂+ M) phases were calculated by each analytical method and compared with the actual fractions used for simulation, or in case of experimental data, with autoradiographic results. Generally, all methods yielded reasonably accurate fractions of cells in each phase with relative errors in the range of 10–20%. However, most methods tended to overestimate G₁ fractions and underestimate S fractions. In addition, variations in the shape of the S phase distribution often caused considerable errors. Phase fractions were also calculated for histograms of kinetically perturbed populations, simulated as well as experimental The errors were only slightly larger than for histograms from asynchronously growing cell populations.     

A CHARACTERIZATION OF THE BOUNDARY BETWEEN THE CYCLING AND RESTING STATES IN ASCITES TUMOR CELLS

Growth deceleration of an Ehrlich ascites tumor with increasing mass is associated with a prolongation of the cell cycle and a decline in the growth fraction. These effects are reversed upon transfer of cells from an older tumor into a new host. Studies were made to locate the stages at which a cell cycle could be suspended or resumed. Transplantation caused a prompt rise in both mitotic and flash H³TdR labeling indices. When all the cells in cycle including mitoses were prelabeled with H³TdR in older tumors, the fraction of labeled mitoses did not decline for a considerable period after transplantation into new hosts. This suggests that the early rise in mitoses is not due to a flow of resting (G_o) cells from a G² store (G²-G^o transition). It appears rather to be a reflection of a lag of the mitotic process relative to other stages during the initial readjustment of the cycle. A prompt rise in flash H₃TdR indices in the transplants suggested cell entry into S from either a suspended G_I (G₁-G_o transition) or a suspended S (S-G_o transition). These possibilities were examined by relating micro-spectrophotometric estimates of DNA to the cell cycle stage as revealed by H³TdR autoradiography. Since G_o cells had DNA values corresponding to G_I, it was concluded that decycling or recycling could occur only after mitosis and before DNA synthesis.     

CELL KINETICS OF A CHONDROSARCOMA

The cell kinetics of the transplantable DC-II mouse chondrosarcoma have been studied by the pulse labelled mitoses method. The analysis gave the following estimates for the phases of the cell cycle: G, 10-5 hr; S, 9-5 hr; G₂, 4 hr with an intermitotic time of 23-5 hr. Consideration of the overall growth of the tumour indicated that the growth fraction and cell loss factor both had values of about 0–5. The results are compared with cell kinetic data from sarcomas and other cartilage tissues.     

Cell cycle analysis of tumour-derived cultures ofCrepis capillaris: A kinetic analogy with proliferation of animal tumour cells

Summary Analysis of the cell cycle by three methods has revealed unusual kinetics of proliferation in tumour derived suspensions ofCrepis capillaris. The different methods of analysis yield different estimates of cycle phase durations, and such discrepancies have been explained in terms of low growth fractions with rapid total cycle traverse. Specifically, confidence in the estimation of G₂ duration by the fraction of labelled mitosis analysis, and comparison with shorter G₂ estimates obtained by the two other methods, suggests that cells drop out in G₁. However, cells which do not drop out of the proliferative compartment traverse G₁ extremely rapidly. Extremely short cell cycle durations in which the G₁ phase is virtually non-existent are uncharacteristic of plant cell suspension cultures, in which the G₁ phase has previously been shown to be extended as compared with meristematic root tip cells. A model has been proposed in which a central core of rapidly dividing cells continuously loses cells into a subpopulation of resting or G₀ cells with the G₁ DNA content. Similarities between plant and animal tumours with respect to cell growth and division are discussed.     

CHARACTERISTICS OF THE ISOPRENALINE STIMULATED PROLIFERATIVE RESPONSE OF RAT SUBMAXILLARY GLAND

A simple stochastic model has been developed to determine the cell cycle kinetics of the isoprenaline stimulated proliferative response in rat acinar cells. The response was measured experimentally, using ³H-TdR labelling of interphase cells and cumulative collections of mitotic cells with vincristine. The rise and fall of the fraction of labelled interphase cells and of metaphase cells is expressed by the product of the proliferative fraction and a difference of probability distributions. The probability statements of the model were formulated and then compared by an iterative fitting procedure to experimental data to obtain estimates of the model parameters. The model when fitted to the combined fraction labelled interphase (FLIW) and fraction metaphase (FMW,) waves gave a mean G_is transit time of 21-2 hr, mean G_is+ S transit time of 270 hr, and mean G_is+ S + G₂ transit time of 35-8 hr for a single injection of isoprenaline, where G_is is the initiation to S phase time. When successive injections of isoprenaline were given at intervals of 24 and 28 hr the corresponding values after the third injection were 12-4 hr, 20-8 hr and 25-7 hr respectively. The variance of the G_is phase dropped from 18-1 to 1–3 while the other variances remained unchanged. The estimated proliferative fraction was 0–24 after a single injection of isoprenaline, and 0–31 after three injections of the drug. Independently determined values of the proliferative fraction, obtained from repeated ³H-TdR injections, were 0–21 and 0–36 respectively.     

THE KINETICS OF GRANULOSA CELLS IN DEVELOPING FOLLICLES IN THE MOUSE OVARY

This investigation describes the kinetics of the granulosa cells in medium-sized follicles type 3b, 4 and 5a in ovaries of 28-day-old Bagg mice. the method of labelling with ³H-thymidine followed by high resolution autoradiography is used in the experimental work, which consist of determining percentage labelled mitosis (PLM-) and continuous labelling (CL-) curves. In order to analyse the data by computer two alternative hypotheses A and B are set up. Both include the assumptions of no cell loss, exponential growth and a resting compartment Q. In hypothesis A cells from Q re-enter the mitotic cycle via the normal DNA-synthesis compartment S_p. Hypothesis B includes beside compartment S_p a special DNA-synthesis compartment S_q where only cells from Q are synthesizing DNA, and these cells re-enter the mitotic cycle via the G₂ compartment. the mean transit time in S_q is considered to be longer than the mean transit time in S_q. On the basis of the hypothesis mathematical expressions for the PLM- and CL-curves are obtained, and by means of a computer the theoretical curves are fitted to the experimental values: thereby all relevant cell kinetical parameters are estimated. Hypothesis B seems to give the best fit between the theoretical and experimental curves. the estimated parameters are: mean cycle times, μ_c= (56.1 hr, 56.1 hr and 22.3 hr for type 3b, 4 and 5a respectively), doubling times, T _D= (96.4 hr, 118.6 hr and 59.1 hr) and the proportion of cells in Q, p _Q = (0.60, 0.71 and 0.69).