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.     

Effect of Vibrio cholerae neuraminidase on the generation of cell-mediated cytotoxicity in vitro.

Vibrio cholerae neuraminidase (VCN))12.5 units/2 X 10(6) cells/ml) continuously present for a standard 5-day MLC will significant (p less than 0.02) increase the cytotoxic activity generated by a given number of responding spleen cells without reducing the specificity. Heat-inactiviated VCN produced no such augmentation. This augmented cytotoxicity could be reproduced by preincubating (1 hr) the responding spleen cells with VCN (25 units/2.5 X 10(6) cells/ml) before addition of stimulating spleen cells. Preincubating the stimulator spleen cells with VCN had no effect. VCN preincubation of target cells or presensitized effector cells produced no augmentation. The addition of soluble VCN to the killing assay also did not increase cytotoxicity. Thus, VCN acts only during the generation of specifically sensitized cytotoxic T cells. When the effect of VCN on MLC reactivity, cell recovery and total cytotoxicity (lytic units/10(6) cells) were compared, it became apparent that VCN increases the proliferation of responder cells after stimulation resulting in both an increased number of cells and also an increase in the proportion of specifically sensitized cytotoxic cells in the culture. VCN treatment of responder cell membrane apparently permits a more ready response to allogenic antigens in culture facilitating both increased proliferation and the increased development of specific cytotoxic killers.     

Red cell aging. II. Anomalous electrophoretic properties of neuraminidase treated human erythrocytes.

Desialylation of human red blood cells (RBC) by Vibrio cholerae neuraminidase (VCN) was found to produce cells with electrophoretic properties which were inconsistent with the view of simple loss of N-acetylneuraminic acid (NANA) as the sole effect of VCN treatment. Modification of human RBC with 50--350 U VCN/10(10) RBC for one hour at 37 degrees C releases 90-100% of the NANA and produces a progressive decrease towards zero in their electrophoretic mobilities when measured in 0.15 M NaCl (pH 7.2) at 25 degrees C. The appearance of positive groups on the desialylated cells was indicated by the VCN-treated cells displaying positive mobilities below approximately pH 5.5 and increased negative mobilities at approximately pH 9 as well as substantial increases in their mobility at neutral pH following treatment with formaldehyde. Adsorption of about 95% of the VCN activity at 0 degrees C to the RBC did not produce any significant change in their electrophoretic mobilities thus indicating that the observed changes in the electrophoretic properties of the RBC following VCN treatment could not be attributable to adsorption of VCN. These studies indicate that the cationic charge groups which appear at the electrophoretic surface of the RBC after VCN treatment are probably of endogenous origin. It is suggested that this alteration rather than simple NANA release may operate to shorten the in vivo survival time of desialylated red cells.     

Statistical evaluation of cell kinetic data from DNA flow cytometry (FCM) by the EM algorithm

Flow cytometric DNA measurements yield the amount of DNA for each of a large number of cells. A DNA histogram normally consists of a mixture of one or more constellations of G0/G1-, S-, G2/M-phase cells, together with internal standards, debris, background noise, and one or more populations of clumped cells. We have modelled typical DNA histograms as a mixed distribution with Gaussian densities for the G0/G1 and G2/M phases, an S-phase density, assumed to be uniform between the G0/G1 and G2/M peaks, observed with a Gaussian error, and with Gaussian densities for standards of chicken and trout red blood cells. The debris is modelled as a truncated exponential distribution, and we also have included a uniform background noise distribution over the whole observation interval. We have explored a new approach for maximum-likelihood analyses of complex DNA histograms by the application of the EM algorithm. This algorithm was used for four observed DNA histograms of varying complexity. Our results show that the algorithm works very well, and it converges to reasonable values for all parameters. In simulations from the estimated models, we have investigated bias, variance, and correlations of the estimates.     

Use of the enzyme neuraminidase for immunotherapeutic treatment of chemically induced carcinogenesis]

The effect of vibrio cholerae neuraminidase (VCN) on the growth of dimethylbenzanthracene-induced sarcoma cells in the inbred CBA mice was investigated. The use of this preparation was started after the appearance of the tumour. Injection of 50 units of VCN twice a week for three months was effective at the early stages of carcinogenesis. An increase of the life-span of mice of comparison with control animals was also observed in the animals inoculated intraperitoneally with induced syngeneic sarcoma cells pretreated with VCN and simultaneously injected (into the developing tumour) with sensitized lymphocytes received from the syngeneic tumour-bearing mice. Lymphocytes were inoculated into the growing tumour. No positive effect ensued when the lymphocytes inoculated into the tumour region were pretreated with VCN. A simultaneous inoculation of neuraminidase into the growing tumour and of syngenous induced-sarcoma cells pretreated with this enzyme intraperitoneally was the most effective. Possibilities of application of neuraminidase under clinical conditions are discussed.     

The problem of cell debris in pulse-cytophotometry: A probability theoretic model

The problem of cell debris in pulse-cytophotometry is considered on the basis of a probability-theoretic model and the results are compared with the model of a pure exponential function for the portions of the DNA fragments in the channels of the flow cyclometer, which is generally used in practice. The model is based on the assumption that the possible decomposition of the DNA content of a single cell into some parts as a consequence of the necessary preparation of the cell material for measurement can be interpreted as a Poisson point process. Therefore the number of divisions of the DNA content of a cell is assumed as a Poisson random variable. It follows a well-defined distribution function of the DNA content in the common population of intact cells and fragments depending on the apriori distribution before pretreatment. This distribution determines a theoretical histogram which can be compared with the measurements. The results corroborate the assumption of a pure exponential course for the portions of DNA fragments only in the area of very small DNA contents. For greater DNA contents there are differences between the pure exponential course and the model described here depending on the a priori distribution and on the intensity of the cell disintegration process. As a byproduct, an estimation of the fraction of disintegrated cells is obtained with respect to the fraction of fragments. A corresponding computer program was applied to a number of experimental histograms. The results have shown that the theory describes the measured histograms satisfactory in most cases.     

Salicylate-Induced Hearing Loss Trigger Structural Synaptic Modifications in the Ventral Cochlear Nucleus of Rats via Medial Olivocochlear (MOC) Feedback Circuit

Lesion-induced cochlear damage can result in synaptic outgrowth in the ventral cochlear nucleus (VCN). Tinnitus may be associated with the synaptic outgrowth and hyperactivity in the VCN. However, it remains unclear how hearing loss triggers structural synaptic modifications in the VCN of rats subjected to salicylate-induced tinnitus. To address this issue, we evaluated tinnitus-like behavior in rats after salicylate treatment and compared the amplitude of the distortion product evoked otoacoustic emission (DPOAE) and auditory brainstem response (ABR) between control and treated rats. Moreover, we observed the changes in the synaptic ultrastructure and in the expression levels of growth-associated protein (GAP-43), brain-derived neurotrophic factor (BDNF), the microglial marker Iba-1 and glial fibrillary acidic protein (GFAP) in the VCN. After salicylate treatment (300 mg/kg/day for 4 and 8 days), analysis of the gap prepulse inhibition of the acoustic startle showed that the rats were experiencing tinnitus. The changes in the DPOAE and ABR amplitude indicated an improvement in cochlear sensitivity and a reduction in auditory input following salicylate treatment. The treated rats displayed more synaptic vesicles and longer postsynaptic density in the VCN than the control rats. We observed that the GAP-43 expression, predominantly from medial olivocochlear (MOC) neurons, was significantly up-regulated, and that BDNF- and Iba-1-immunoreactive cells were persistently decreased after salicylate administration. Furthermore, GFAP-immunoreactive astrocytes, which is associated with synaptic regrowth, was significantly increased in the treated groups. Our study revealed that reduced auditory nerve activity triggers synaptic outgrowth and hyperactivity in the VCN via a MOC neural feedback circuit. Structural synaptic modifications may be a reflexive process that compensates for the reduced auditory input after salicylate administration. However, massive increases in excitatory synapses in the VCN may represent a detrimental process that causes central hyperactivity, leading to tinnitus.     

Variance within homogeneous phytoplankton populations, II: Analysis of clonal cultures

In part I of this series of articles, a framework was presented for interpreting histograms of volume or fluorescence as measured by a flow cytometer on homogeneous phytoplankton populations. In this paper, the analytical framework is applied to flow cytometric histograms from laboratory experiments involving clonal phytoplankton cultures. The density function derived in part I was modified to include a third parameter representing a linear shift of the origin. This modified density function was fitted to chlorophyll fluorescence histograms for populations believed to be asynchronous (grown in continuous light) and also to histograms from populations grown on a 14:10 (h:h) light/dark cycle. Near-synchronous subpopulations sorted from an asynchronous population were also analyzed. In populations in which underlying assumptions (asynchronous divisions, constant growth) are valid, curve fits provide estimates of the inherent variability among cells at age 0. The implication of fitting the density function to populations in which these assumptions are not valid is discussed.     

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.     

Acrylonitrile interaction with testicular DNA in rats.

In the present study we report the in vivo interaction of acrylonitrile (VCN) with testicular tissue in rats. Covalent binding of radioactivity to testicular tissue DNA was examined for a period of 72 hr after a single oral dose (46.5 mg/kg) of [2,3-14C] VCN. Maximal covalent binding was observed at 0.5 hr (8.9 mumol VCN equivalent/mol nucleotide). Binding decreased gradually thereafter but was still detected (2.5 mumol VCN equivalent/mol nucleotide) at 72 hr following VCN administration. Further, we examined the effects of VCN on DNA synthesis and repair in the testes of rats following a single oral dose (46.5 mg/kg) of VCN to clarify the impact of the covalent binding observed on the testicular genetic material. A significant decrease in DNA synthesis (80% of control) was observed at 0.5 hr after treatment. At 24 hr following acrylonitrile administration, testicular DNA synthesis was severely inhibited (38% of control). Testicular DNA repair was increased 1.5-fold at 0.5 hr and more than 3.3-fold at 24 hr following treatment with VCN. These results suggest that VCN can act as a multipotent genotoxic agent by alkylating DNA in testicular tissue and may affect the male reproductive function by interfering with testicular DNA synthesis and repair processes.     

Partial phase synchronization of neural populations due to random Poisson inputs

We show that populations of identical uncoupled neurons exhibit partial phase synchronization when stimulated with independent, random unidirectional current spikes with interspike time intervals drawn from a Poisson distribution. We characterize this partial synchronization using the phase distribution of the population, and consider analytical approximations and numerical simulations of phase-reduced models and the corresponding conductance-based models of typical Type I (Hindmarsh-Rose) and Type II (Hodgkin-Huxley) neurons, showing quantitatively how the extent of the partial phase synchronization depends on the magnitude and mean interspike frequency of the stimulus. Furthermore, we present several simple examples that disprove the notion that phase synchrony must be strongly related to spike synchrony. Instead, the importance of partial phase synchrony is shown to lie in its influence on the response of the population to stimulation, which we illustrate using first spike time histograms.     

Molecular interaction of [2,3-14C] acrylonitrile with DNA in gastric tissue of rat

Acrylonitrile (VCN) is used extensively in polymer industries, and is known to induce gastric cancer following oral administration, A paucity of information exists regarding the mechanism(s) by which acrylonitrile induces gastric neoplasia. The time course for uptake of radioactivity by gastric tissue and covalent binding of [2,3-¹⁴C] VCN or its metabolites to gastric DNA were determined following a single oral dose of 46.5 mg/kg. The rates of DNA synthesis and repair, as measured by unscheduled DNA synthesis in the gastric tissue of VCN-treated rats, were also studied. Maximum tissue uptake and covalent binding of radioactivity to gastric DNA were observed at 15 minutes following [2,3-¹⁴C] VCN administration. At 6 hours following VCN administration, significant inhibition (37% of control) in gastric replicative DNA synthesis was observed. A rebound followed by an increase (211% of control) in replicative DNA synthesis was observed at 24 hours. A three-fold elevation in unscheduled DNA synthesis was observed at 24 hours following treatment with VCN. These results indicate that VCN or its metabolites irreversibly interact with gastric DNA, causing DNA damage. The results also indicate that the delayed VCN-induced DNA repair, determined as unscheduled DNA synthesis, is inefficient for the removal of the resulting DNA lesions.     

Acrylonitrile interaction with testicular DNA in rats

In the present study we report the in vivo interaction of acrylonitrile (VCN) with testicular tissue in rats. Covalent binding of radioactivity to testicular tissue DNA was examined for a period of 72 hr after a single oral dose (46.5 mg/kg) of [2, 3-¹⁴C] VCN. Maximal covalent binding was observed at 0.5 hr (8.9 μmol VCN equivalent/mol nucleotide). Binding decreased gradually thereafter but was still detected (2.5 μmol VCN equivalent/mol nucleotide) at 72 hr following VCN administration. Further, we examined the effects of VCN on DNA synthesis and repair in the testes of rats following a single oral dose (46.5 mg/kg) of VCN to clarify the impact of the covalent binding observed on the testicular genetic material. A significant decrease in DNA synthesis (80% of control) was observed at 0.5 hr after treatment. At 24 hr following acrylonitrile administration, testicular DNA synthesis was severely inhibited (38% of control). Testicular DNA repair was increased 1.5-fold at 0.5 hr and more than 3.3-fold at 24 hr following treatment with VCN. These results suggest that VCN can act as a multipotent genotoxic agent by alkylating DNA in testicular tissue and may affect the male reproductive function by interfering with testicular DNA synthesis and repair processes.     

Molecular interaction of acrylonitrile and potassium cyanide with rat blood

The interaction of acrylonitrile (VCN) with rat blood has been investigated at the molecular level in an attempt to understand the possible mechanism of its toxicity. The results obtained were compared to those with potassium cyanide (KCN), a compound known to liberate cyanide (CN^?) in biologic conditions. The radioactivity derived from K¹⁴CN was eliminated faster than that from [1-¹⁴C]VCN. Up to a maximum of 94% of ¹⁴C from VCN in erythrocytes was detected covalently bound to cytoplasmic and membrane proteins, whereas 90% of the radioactivity from KCN in erythrocytes was found in the heme fraction of hemoglobin. Determination of specific activity showed that binding occurred more in vivo than in vitro which indicated that the VCN molecule was bioactivated inside erythrocytes. These results indicate that KCN interacts mainly through CN^? liberation and binding to heme, whereas VCN, which binds to cytoplasmic and membrane proteins, may cause damage to red cells by mechanisms other than release of CN^?.     

Modeling modification of chemical mutagenesis in human cells. 1. Elimination of chromosomal aberrations]

Elimination of chromosome aberrations was studied in populations of dividing cells. For this purpose, on the basis of the corresponding theory of Carrano-Heddle assuming the Poisson distribution, a theory is advanced by the authors based on geometrical distribution, describing the distribution of lesions caused by the action of tioTEF. Parameters of elimination are obtained observed in case of addition of a protector, aminopropy-laminoethylthiophosphoric acid (APAETF) in a lymphocyte culture of human peripheral blood. It is shown that the addition of the protector diminishes the probability of the transmission of chromosome aberrations to daughter cells.     

Augmentation of immune response to syngeneic fibrosarcoma T241 with in vivo protection

Summary Lewis T241 fibrosarcoma, a syngeneic tumor in C57 BL/6J mice, was found to be poorly immunogenic. When tumor-bearing animals (TBA) were challenged with tumor cells either concomitantly or after excision of a growing tumor no protection was observed. In vivo (Winn) neutralization assays also showed a lack of tumor immunogenicity. However, in vitro studies showed that a significant proliferative response could be elicited from the spleen cells of TBA when these cells were cultured with either mitomycin-C-treated tumor cells or KCl tumor extract. Similarly, macrophage migration inhibition factor (MIF) was produced by TBA spleen cells upon incubation with KCl tumor extract, but no cell-mediated cytotoxicity to T241 target cells was observed with various lymphoid cell populations at any stage of tumor growth. Immunization of syngeneic animals with Vibrio cholerae neuraminidase(VCN)-treated, irradiated tumor cells alone or admixed with Freund's complete adjuvant (FCA) resulted in decreased tumor growth and fewer pulmonary metastases following challenge with 10⁶ tumor cells. No complete tumor rejection was observed. In contrast, 13 of 16 animals immunized with irradiated tumor cells admixed with FCA rejected 10⁵ tumor cells. Animals that grew tumors had significantly reduced tumor growths and pulmonary metastases. Lymph node and peritoneal exudate cells (PEC) of immunized animals showed significant cytotoxicity to T241 cells.     

THE APPLICATION OF AGE DISTRIBUTION THEORY IN THE ANALYSIS OF CYTOFLUORIMETRIC DNA HISTOGRAM DATA

Age distribution theory has been employed in a model to analyse a variety of histograms of the DNA content of single cells in samples from experimental tumours growing in tissue culture. The method has produced satisfactory correspondence with the experimental data in which there was a wide variation in the proportions of cells in the intermitotic phases, and generally good agreement between the ³H-thymidine labelling index and the computed proportion in S phase. The model has the capacity to analyse data from populations which contain a proportion of non-cycling cells. However, it is concluded that reliable results for the growth fraction and also for the relative durations of the intermitotic phase times cannot be obtained for the data reported here from the DNA histograms alone. To obtain reliable estimates of the growth fraction the relative durations of the phase times must be known, and conversely, reliable estimates of the relative phase durations can only be obtained if the growth fraction is known.     

Effect of drug efflux blockers on vital staining of cellular DNA with Hoechst 33342

The present study shows that staining of certain live cells, e.g., adriamycin-resistant P388 cells, by Hoechst 33342 is difficult because of the presence of a rapid efflux pump, which reduces intracellular dye concentration. Coincubation of these refractory cells in the presence of efflux blockers such as phenothiazines (trifluoperazine) or Ca++ channel blockers (verapamil) enhances dye retention and thus leads to generation of normal DNA distribution histograms. Laser flow cytometric data is confirmed by fluorometric assays, which show that P388/R cells retain one-third the amount of Hoechst 33342, and coincubation with efflux blockers increases Hoechst retention to values similar to those of drug-sensitive P388 cells. DNA histograms of mouse splenocytes incubated with Hoechst 33342 alone have a bimodal distribution possibly because of the presence of subpopulations that do not retain the fluorochrome owing to rapid efflux. Coincubation with an efflux blocker results in the generation of unimodal DNA histograms from these cells. These preliminary studies suggest that reduced retention of Hoechst 33342 in certain cell types (because of rapid efflux) can be blocked by efflux blockers, thus leading to generation of typical DNA distribution histograms.