A practical graphical method for estimating the fraction of cells in S in DNA histograms from clinical tumor samples containing aneuploid cell populations

A graphical method for the analysis of unperturbed DNA histograms is presented in which the area of the normalized histogram subtended by the fraction of cells in S is represented by a trapezoid whose dimensions are dependent on features common to all such histograms. The technique takes measurement variability into account. This method was applied to a variety of synthetic DNA histograms. Overall, calculated values for the fraction of cells in S correlated well with actual values. This method was applied to 36 diploid cases of non-Hodgkin lymphoma; results correlated well with those obtained by a computer-based method. The results of the graphical-method were also highly reproducible between different observers. The graphical method can be used in the presence of aneuploid cell populations. Techniques for calculating S fractions in the presence of aneuploidy in clinical samples are described. These techniques were applied to synthetic histograms of mixed diploid and aneuploid populations. Calculated values correlated well with actual values.     

Different opinions on classification of DNA histograms produced from paraffin-embedded tissue

Flow cytometric DNA ploidy determination has been regarded as an objective prognostic parameter in several types of human cancer. To test whether DNA histograms are similarly interpreted, a series of flow cytometric DNA histograms was posted to six investigators working in the field for independent classification. The histograms were produced from paraffin-embedded adrenal adenomas or non-neoplastic tissue and had several different patterns. Only 44% of the histograms were similarly classified by all investigators, and 85% by five of the six participants, when DNA ploidy was evaluated. Different criteria for tetraploidy existed, and also some uncertainty in classifying peridiploid and small aneuploid peaks. It is concluded that lack of consensus on histogram classification may result in widely varying percentages of DNA aneuploid tumors found even if the data are similar. Until general agreement is reached on the definition of DNA aneuploidy and its subclasses, classification of DNA histograms is variable and subjective.     

Reliability of DNA cytometric S-phase analysis in surgical biopsies: assessment of systematic and sampling errors and comparison between results obtained by image and flow cytometry

Three major parameters in DNA histograms that contribute to the reliability of S-phase analysis were evaluated. These parameters are (1) the extent of background in relation to the amount of S-phase cells (and the validity of its subtraction), (2) the size of the "free" S-phase range (S(free)), and (3) the sampling error of cell counting. Tests in histograms obtained from surgical biopsies by flow cytometry (FCM) showed that the background subtraction is reliable if the found S-phase fraction is higher than the fraction of background events in the histogram range of the cell population. The size of S(free) was determined in computer-generated test histograms as a function of variables such as the coefficient of variation (CV) and the DNA index (DI). To calculate the sampling error of cell counting above background and in S(free), a model was developed that was validated by experimental data. This error can serve as an indicator of the uncertainty in S-phase analysis. The poor correlation found between %S values measured by image cytometry (ICM) and FCM in surgical biopsies was assigned to high uncertainty by low cell numbers in ICM histograms. A method is proposed to estimate quantitatively the reliability of S-phase analysis that can facilitate the interpretation of results.     

Flow cytometry and cell proliferation kinetics

Flow cytometric techniques are presented which allow to determine parameters of cell proliferation kinetics by means of histogram sequences after special manipulations of the cell culture under investigation: (a) In the stathmokinetic method metaphase blocking agents are applied which allow the cells of the population to continue progression through interphase and accumulate at 4C DNA content. The development of DNA specific histograms during this process is analysed as to the G1 phase duration and the fraction of nonproliferating cells. (b) In the BUdR/Hoechst method the suppression of Hoechst fluorescence after BUdR incorporation during S phase is taken as a means for inducing a temporal change of histogram shapes without perturbing the cell cycle progression of the population. This temporal development of histogram shapes is analysed as to phase duration, whole cycle time and fraction of nonproliferating cells. (c) By combining the BUdR/Hoechst technique with a simultanous DNA specific stain and analysing with a two-parametrical flow cytometer, more information is obtained from each histogram after BUdR incorporation: The location of cells in the cycle at the beginning of the experiment, the cycle stage at cell harvest, and from this the distance and velocity of progression through the cycle during drug incubation. By introduction of these dynamic methods flow cytometry has become a powerful tool for the study of cell proliferation kinetics in culture.     

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.     

Estimation of cell DNA synthesis rate from flow-cytometric histograms

The rate of DNA synthesis along S-phase can be estimated on the basis of a simple model of a cell population, from a single FCM histogram under conditions of balanced exponential growth, and from two histograms suitably spaced in time under conditions of locally exponential S-phase influx. An algorithm that constructs a piecewise constant approximation of the synthesis rate and estimates the phase fractions from a single histogram is presented. The results of the application of the algorithm to a set of computer generated histograms in balanced exponential growth are reported.     

Dependence of DNA-histograms on the sampling techniques in fine needle aspirates of the breast.

48 fine needle aspiration biopsy (FNAB) samples from 25 breast cancer cases, originally used for cytodiagnosis were subjected to DNA cytometry. There were air dried smears stained with the MGG method, and samples stained with HE or PAP stain after 50% ethanol fixation and cytocentrifugation. Different sampling strategies were applied. Four methods were tested: method 1: cell groups measured, method 2: all cells measured, method 3: free cells measured, and method 4: atypical free cells measured. Method 4 showed most often DNA aneuploid histogram patterns, sampling method 1 had the highest number of DNA diploid histogram patterns. Diagnostic approaches may benefit from a sampling method detecting the hiding aneuploid cell population. Grading of neoplasm could potentially benefit from other approaches.     

Optimizing flow cytometric DNA ploidy and S-phase fraction as independent prognostic markers for node-negative breast cancer specimens.

Developing a reliable and quantitative assessment of the potential virulence of a malignancy has been a long-standing goal in clinical cytometry. DNA histogram analysis provides valuable information on the cycling activity of a tumor population through S-phase estimates; it also identifies nondiploid populations, a possible indicator of genetic instability and subsequent predisposition to metastasis. Because of conflicting studies in the literature, the clinical relevance of both of these potential prognostic markers has been questioned for the management of breast cancer patients. The purposes of this study are to present a set of 10 adjustments derived from a single large study that optimizes the prognostic strength of both DNA ploidy and S-phase and to test the validity of this approach on two other large multicenter studies. Ten adjustments to both DNA ploidy and S-phase were developed from a single node-negative breast cancer database from Baylor College (n = 961 cases). Seven of the adjustments were used to reclassify histograms into low-risk and high-risk ploidy patterns based on aneuploid fraction and DNA index optimum thresholds resulting in prognostic P values changing from little (P < 0.02) or no significance to P < 0.000005. Other databases from Sweden (n = 210 cases) and France (n = 220 cases) demonstrated similar improvement of DNA ploidy prognostic significance, P < 0.02 to P < 0.0009 and P < 0.12 to P < 0.002, respectively. Three other adjustments were applied to diploid and aneuploid S-phases. These adjustments eliminated a spurious correlation between DNA ploidy and S-phase and enabled them to combine independently into a powerful prognostic model capable of stratifying patients into low, intermediate, and high-risk groups (P < 0.000005). When the Baylor prognostic model was applied to the Sweden and French databases, similar significant patient stratifications were observed (P < 0.0003 and P < 0.00001, respectively). The successful transference of the Baylor prognostic model to other studies suggests that the proposed adjustments may play an important role in standardizing this test and provide valuable prognostic information to those involved in the management of breast cancer patients.     

Flow cytometric S-phase fraction measurement in breast carcinoma: Influence of software and histogram resolution

BACKGROUND: S-phase fraction (SPF) measurement by flow cytometry is a clinically useful prognostic factor in patients with breast carcinoma. Standardized SPF determination is essential. As part of a multicenter study, we evaluated the influence of the choice of software and histogram resolution (256, 512, or 1,024 channels) on SPF quantification. METHODS: One hundred thirty-three DNA histograms were analyzed in three laboratories with Modfit 5.2, Modfit LT, and Multicycle AV software. Strict rules for histogram interpretation and software management were applied. The following five options were compared: MF 5.2 1024, MF 5.2 256, MF LT 256, MC AV 256, and MC AV 512. RESULTS: In the DNA diploid and aneuploid groups, SPF distributions were not statistically different among the five options. Excellent quantitative correlations were obtained between pairs of options. When using tertiles as cutpoints for SPF classification, concordance rates ranged from 79.7% to 93.2% for DNA diploid samples and from 87.8% to 95.9% for DNA aneuploid samples, the best results being obtained with software working with a similar histogram resolution. CONCLUSIONS: Standardized use of commercially available software, including the choice of histogram resolution, provides comparable SPF results.     

Prospective clinical study comparing DNA flow cytometry and HPV typing as predictive tests for persistence and progression of CIN I/II.

A cohort of 70 consecutive women at a university hospital colposcopy clinic with untreated CIN I and CIN II (CIN I/II) confirmed by cytology and histology was followed for 1 year in the setting of a prospective trial. In the lesions, the presence of DNA from HPV types was examined by restriction fragment length polymorphism (RFLP) analysis. Aneuploid cell lines were demonstrated by aneuploid histograms generated by high-resolution DNA flow cytometry. HPV type 16 infection and the existence of aneuploid cell lines proved to be significant risk factors for CIN I/II lesions to persist or progress to CIN III in the 1-year follow-up period in the same cohort of patients. The relative risks and 95% confidence intervals (CI) were 1.81 (1.44-2.76) for aneuploid cell lines and 1.74 (1.10-2.76) for HPV type 16 infection in CIN I/II lesions. As a predictive diagnostic test for CIN I/II lesions to persist or progress, the specificity and positive predictive value (PPV) for aneuploid histograms were 100% (CI, 73.5-100%) and 100% (CI, 86.8-100%), respectively. The low sensitivity of 27.3% (CI, 14.9-42.8%) restricted the clinical application of the test, leaving 32 of 44 women with persisting or progressing CINI/II with diploid histograms. HPV type 16 positivity by FRLP had a PPV of 68.4% (CI, 43.5-87.4%) as a prognostic test. Six of 19 HPV 16 infected women showed complete remission of their CIN lesion. A combination of the two tests did not provide any additional information.     

Two-color multiparametric method for flow cytometric DNA analysis of carcinomas using staining for cytokeratin and leukocyte-common antigen

Solid tumors contain heterogenous cell populations, resulting in flow cytometric (FCM) DNA quantitations of a mixture of tumor and host cells. Such mixed populations can result in dilution of the tumor cells by the host cells, in difficulty defining the diploid reference mean and in histogram peak overlap, precluding cell-cycle analysis. In this study, epithelial (tumor) cells and contaminating host cells in 100 consecutively accessioned human mammary and colorectal carcinomas were segregated in a multiparametric two-color FCM DNA analysis of intact, ethanol-fixed cells. These two carcinomas and bladder carcinomas contain a cytoskeleton of simple epithelium that is selectively stained with an FITC-labeled monoclonal antibody (MAb) to cytokeratin (CK: CAM 5.2-FITC). This MAb detects the CK 8, CK 18 and CK 19 consistently present in all layers of normal and neoplastic urothelium, colonic epithelium and mammary epithelium. Gating on CK in these tumors enables the nonstaining leukocytes, stromal fibroblasts and endothelial cells to be excluded from DNA analysis. A separate aliquot of each tumor evaluated was labeled with an MAb to leukocyte-common antigen (LCA-FITC) to serve as a patient-specific intrinsic diploid reference standard. Both the CK-labeled and LCA-labeled cells were then dual labeled for DNA with propidium iodide. This method (1) correctly identified the intrinsic diploid (LCA-positive) channel, allowing an accurate definition of normal cell DNA content for calculation of the DNA index; and (2) resulted in an increased sensitivity in the identification of both diploid and abnormal hyperdiploid tumor cell populations. It also (3) limited DNA cell cycle analysis to urothelial, colonic and mammary epithelial cells, the majority of which were neoplastic in carefully selected tumor samples. In addition, this method (4) clarified near-tetraploid populations that overlap the normal nonepithelial G2M region by diminishing the normal G2M peak and accentuating the aneuploid tetraploid G0G1 peak and (5) deconvoluted overlapping histograms composed of normal host and diploid-range or aneuploid tumor cells by gating on tissue-specific markers. This exclusion of host cells in both classes of tumors resulted in more accurate cell-cycle calculations in the former and allowed calculation of the S-phase fractions in the latter.     

Methodologic aspects of DNA assessment by means of image cytometry in tumors of the salivary glands. A comparison between the results obtained using sections and cytospin preparations from the same paraffin-embedded specimens

The image cytometric nuclear DNA assessments on paraffin-embedded tissue sections and on Cytospin preparations of disaggregated specimens from the same cases were compared in 98 salivary gland tumors, including 21 acinic cell carcinomas, 29 mucoepidermoid carcinomas, 21 adenocarcinomas and 27 adenoid cystic carcinomas. The histogram type (diploid, tetraploid or aneuploid) and the number of cells with DNA values greater than 2.5c (expressed in relative units) were considered as variables in the correlation. A high correlation between the results in different specimens was found in acinic cell carcinomas, mucoepidermoid carcinomas and adenocarcinomas; the histogram type and the number of cells with DNA values greater than 2.5c were essentially the same between specimen types in these three tumor entities. The cases of adenoid cystic carcinomas showed a considerably lower degree of correlation: in 8 of the 27 cases, the Cytospin preparations yielded diploid histograms, while the tissue sections yielded aneuploid histograms. The number of cells with DNA values greater than 2.5c was notably lower in the Cytospin preparations from adenoid cystic carcinoma; the reasons for this exceptional behavior of the cells of adenoid cystic carcinoma are discussed. These findings demonstrate that paraffin-embedded specimens of different tumor entities, even from the same organ, can be affected differently by disaggregation procedures. While retrospective studies on disaggregated paraffin-embedded specimens can yield reliable results, comparative assessments using both DNA analysis techniques, as in this study, should be performed before a large number of cases is evaluated.     

Measurement variability in DNA flow cytometry of replicate samples

A Bladder Cancer Flow Cytometry Network study has been carried out aimed at identification of the sources of inter- and intralaboratory variability. Replicate "cocktail" samples containing a mixture of peripheral blood lymphocytes and an aneuploid cell line and samples of peripheral blood lymphocytes serving as a DNA reference standard were distributed to five network laboratories. The samples were stained for DNA using propidium iodide, with each laboratory using its own staining protocol. Sets of these samples were analyzed by flow cytometry to obtain cellular DNA distributions. DNA index and hyperdiploid fraction were calculated for each histogram using an automated technique. Results were evaluated by analysis of variance to identify sources of variability. Three important sources of variation were found that affect flow cytometry in general and- the transportability of flow cytometry results to routine clinical use in particular. The significant variation among laboratories that is constant across time most probably represents stable differences in instrumentation, instrument set-up, and laboratory techniques. This variation can be compensated for, if it is known and stable, to develop transportable classification criteria. The second type of variation, termed the interaction component, represents differences among laboratories that are not constant across time. Sources of this variation include inconsistency in sample preparation, staining, and analysis. The elimination of this type of variation is required for meaningful comparison of data within and among laboratories and the creation of interlaboratory data-bases. The third type of variation represents pure measurement variability and affects the sensitivity of the technique.     

Kinetics of protein and DNA synthesis studied by mathematical modelling of flow cytometric protein and DNA histograms

Abstract. Mathematical models for histograms of cellular protein content as measured by flow cytometry were developed, based on theoretical protein distributions. These were derived from the age distribution of cells and the accumulation function for cellular protein content as a function of age within the cell cycle. A model assuming an exponential age distribution and an exponential protein. accumulation function was found to give the best representation of protein histograms of exponentially growing NHIK 3025 cells. This is in good agreement with the known kinetic behaviour of such cells. By the combined use of the protein histogram model and a similar model for DNA content, and assuming linear DNA accumulation during S, the fraction of cells in S, as a function of cellular protein content, was simulated. This function showed good agreement with values of the [³H]TdR labelling index scored in cells sorted by flow cytometry from 5-channel intervals of the protein histogram. The protein and DNA histogram models were combined into a two-dimensional model for correlated protein/DNA measurements. Comparison between simulated data and experimentally derived two-dimensional protein/DNA histograms gave further support to the cell kinetic assumptions underlying the models, but also identified some minor deviations which could not be recognized in the analysis of the one-dimensional histograms.     

A method of DNA histogram analysis by computer and its evaluation by simultaneous combination with 3H-thymidine autoradiography

A method for determining the fractions of cells in the G1, S, and G2 + M phases of the cell life cycle, by quantifying DNA histograms derived from static fluorescence cytophotometry, was evaluated by simultaneous combination with 3H-thymidine autoradiography. DNA histograms were obtained by cytofluorometry on the Feulgen-stained autoradiographs of HeLa cells, and mouse and rat hepatocytes, after DNA labelling with 3H-thymidine. The synthetic histogram determined by "sum of discrete normal curves" technique was fitted to the experimental data according to a weighted least-squares method by a desk-top computer (HP 85F). The mean relative percent deviations of estimated cell cycle phase fractions from the actual phase fractions determined directly on an autoradiograph was 6.6 +/- 3.3%.     

Small-cell-type poorly differentiated squamous cell carcinoma of the lung. Cytologic, immunohistochemical and nuclear DNA content analysis.

The relationship between the nuclear DNA content, the immunohistochemical findings, the clinical characteristics (tumor volume doubling time and survival) and the cytomorphologic features of small cell poorly differentiated squamous cell carcinoma of the lung was studied in ten cases. There were no significant correlations between the immunohistochemical stainings for neuron-specific enolase and keratin and the clinical characteristics in these cases. The DNA histogram patterns were classified as type I or type II, depending on the degree of dispersion of values. There was no relationship between the immunohistochemical findings and the DNA histogram patterns. Only the DNA histogram patterns were related to some of the clinical characteristics: patients with type II histograms had significantly shorter tumor volume doubling times than did patients with type I histograms. Such information may aid in distinguishing the small cell type of poorly differentiated squamous carcinoma from classic small cell carcinoma of the lung, with which it may be confused.     

Flow cytometric DNA content of fresh tumor specimens using keratin-antibody as second stain for two-parameter analysis.

Studies concerning flow cytometric assessed DNA content reveal problems in interpretating DNA histograms of tumor specimens. The main problems are histograms with a broad coefficient of variation in the G0/G1 fraction; a high G2M fraction and samples with a low percentage of tumor cells. Therefore, in the present study, 382 fresh tumor specimens of carcinomas were analysed routinely, double labeled with, on the one hand, propidium-iodide for assessing DNA content and, on the other, a monoclonal keratin-antibody for marking epithelial and tumor cells. Of the 311 tumor samples, using single parameter analysis 165 (54%) were classified as DNA aneuploid and 146 (46%) as DNA "euploid." By double parameter analysis, 224 (72%) samples were keratin positive and 87 (27%) keratin negative and, of the 224 keratin positive tumors, 175 (78%) were DNA aneuploid and 49 (22%) DNA euploid. The DNA histograms of single and double parameter analysis were compared and it was concluded that in 24 cases (11%) keratin labeling was necessary to recognize DNA aneuploidy. In another 23 (10%) cases, keratin labeling was helpful in assessing DNA aneuploidy. Finally when the results of the 311 samples were combined, 215 (68%) were scored as DNA aneuploid and 99 (32%) DNA euploid. Thus the overall gain in assessing DNA aneuploidy using the double labeling technique is 14%. In conclusion, it is shown that keratin labeling on fresh tumor cell suspensions of epithelial tumors is of additional value in establishing DNA content. Because single parameter DNA assessment is adequate in approximately 60% of the tested samples, the double labeling technique can be performed routinely, or after initial single parameter DNA assessment. Histograms having a broad CV and/or a high G2M are good candidates for the double labeling technique. Using this technique, DNA-content assessment becomes more reliable.     

Correlation of DNA distribution and cytological differentiation of human cervical carcinomas

The DNA distribution of biopsy specimens from 46 patients suffering from cervical carcinoma was analysed by flow cytometry and compared with the cytological differentiation. According to morphological criteria the carcinomas were classified as highly differentiated, moderately differentiated and poorly differentiated. The results demonstrate that highly differentiated tumours contain hyperploid cells predominantly with hyperdiploid DNA content. Hyperploid cell populations in the moderately differentiated tumours are mainly in the hyperdiploid and tetraploid regions. Poorly differentiated tumours contain hypertetraploid and aneuploid cell populations. The significance of these findings is discussed.     

Flow cytometric bivariate analysis of DNA and cytokeratin in colorectal cancer.

Different opinions about flow cytometric estimates of DNA aneuploidy and/or S-phase fraction (SPF) as supplementary prognostic markers in colorectal cancer are to some degree associated with methodology. Using univariate DNA analysis, we have previously investigated the DNA ploidy in colorectal cancer, its heterogeneity within and between tumors and its relation to survival. To improve detection of DNA aneuploid subpopulations and particularly estimation of their SPF's we investigated a method for bivariate DNA/cytokeratin analysis on fine-needle aspirates of 728 frozen biopsies from 157 colorectal tumors. Unfixed aspirates were stained with propidium iodide and FITC-conjugated anti-cytokeratin antibody in a saponin-buffer. A significant association between SPF and debris was observed. There were no substantial difference in DNA ploidy patterns between univariate and bivariate measurements (concordance was 92-95%). No new DNA aneuploid subpopulations were detected in cytokeratin-gated compared to ungated or univariate histograms. Debris-adjusted SPF's of cytokeratin-gated histograms were significantly higher than of ungated histograms, also for subpopulations with DI>1.4 (p<0.0001). There was no significant association between SPF and survival.