A computer program for analyzing bivariate flow karyotypes

This article describes a computer program for analyzing bivariate flow karyotypes of human chromosomes stained with Hoechst 33258 (HO) and chromomycin A3 (CA). The karyotype first is divided into regions that contain chromosome peaks. The chromosomes that are associated with those areas are identified. The distributions in these areas then are fitted with mathematical functions of increasing complexity. The process starts by fitting a specified number of univariate Gauss functions to projections of the HO and CA distributions of each area. The final fit can include multiple bivariate Gauss functions, including a background function for debris subtraction. The results of one stage in the fitting process serve as seed values for the next, more complex step. Since the program autonomously estimates the starting values for the iterative fitting procedures, the fit results are insensitive to operator bias and the program will consistently converge to the same solutions. The resulting table of parameter values can be used to compare flow karyotypes to a reference data set.     

Semi-automated detection of aberrant chromosomes in bivariate flow karyotypes.

A method is described that is designed to compare, in a standardized procedure, bivariate flow karyotypes of Hoechst 33258 (HO)/Chromomycin A3 (CA) stained human chromosomes from cells with aberrations with a reference flow karyotype of normal chromosomes. In addition to uniform normalization of normal and abnormal flow karyotypes, the main purpose is detection of structurally abnormal chromosomes in often complex karyotypes of tumor cells. The method, which has been implemented in a computer program, consists of a comparison of individual chromosome peaks with the positions of peaks in the flow karyotype constituted by normal chromosomes and takes into account the natural variability in base composition of normal chromosomes among healthy individuals. Flow-karyotypes are normalized using an iterative fitting procedure, using corrections for (1) amplification of HO and CA fluorescence, (2) cross-talk between the fluorescence signals of HO and CA, and (3) offset of the HO and CA origin. Flow karyotypes of two cell lines, one with a simple deletion and the other with more complex karyotypic changes, were analyzed. The results of flow analysis were found to be in general agreement with the cytogenetic analysis of quinacrine banded karyotypes.     

High resolution chromosome analysis: one and two parameter flow cytometry

Isolated mammalian chromosomes have been quantitatively classified by high resolution flow cytometry. Chinese hamster chromosomes stained with 33258 Hoechst and excited in the UV showed a fluorescence distribution in which the 14 types of Chinese hamster chromosomes were resolved into 16 groups seen as distinct peaks in the distributions. Chinese hamster chromosomes were also stained with both 33258 Hoechst (HO) and chromomycin A3 (CA3); the two dye contents were measured by selective excitation in the UV and at 458 nm in a dual beam flow cytometer. The resulting two parameter distribution (HO versus CA3) showed 10 chromosome groups¹. Human strain LLL 761 chromosomes stained with HO and excited in the UV showed a fluorescence distribution in which the 23 types of human chromosomes were resolved into 12 groups. Human chromosomes stained with both HO and CA3 and measured in the dual beam flow cytometer produced two parameter fluorescence distributions which showed 20 groups. The chromosomes associated with each group were determined by quinacrine banding analysis of sorted chromosomes and by DNA cytophotometry of preidentified metaphase chromosomes. The relative HO and CA3 stain content and frequency of occurrence of chromosomes in each group were determined from the fluorescence distributions and compared to the results from DNA cytophotometry. The chromosome to chromosome variations in HO and CA3 staining are attributed to variations in chromosomal base composition.     

Different Hoechst 33342 and DAPI fluorescence of the human Y chromosome in bivariate flow karyotypes

Summary The staining properties of AT-specific dyes Hoechst 33342 and DAPI as revealed by Hoechst 33342/mithramycin and mithramycin/DAPI bivariate human flow karyotype patterns are different for chromosomes rich in heterochromatin. The peak corresponding to chromosome Y of a given cell line is higher on the A/T axis with mithramycin/DAPI staining than with mithramycin/Hoechst. The chromosome 1 was found slightly more fluorescent in mithramycin/DAPI than in mithramycin/Hoechst as judged by the slight displacement of its area on the Hoechst-DAPI axis. The other peaks did not show major differences. On the same flow karyotypes, chromosomal rearrangements were detected.     

Counterstaining human chromosomes for flow karyology

Isolated human metaphase chromosomes stained with the fluorochromes 4'-6-diamidino-2-phenylindole (DAPI) and chromomycin A3(CA3), and counterstained with nonfluorescent netropsin (NTR), have been analyzed by dual-laser flow cytometry. Counterstaining with NTR reduces DAPI fluorescence except at regions on chromosomes 1,9,15,16, and Y, corresponding to C-band heterochromatin. Bivariate flow karyology of human chromosomes treated with this triple-stain combination resolves chromosomes 1,9, and Y distinctly from the remaining chromosomes and resolves variations between chromosome homologues not detected by staining with propidium iodide (PI) or with the double stain combination Hoechst 33258(HO) and CA3.     

Chromosome heteromorphism quantified by high-resolution bivariate flow karyotyping.

Maternal and paternal homologues of many chromosome types can be differentiated on the basis of their peak position in Hoechst 33258 versus chromomycin A3 bivariate flow karyotypes. We demonstrate here the magnitude of DNA content differences among normal chromosomes of the same type. Significant peak-position differences between homologues were observed for an average of four chromosome types in each of the karyotypes of 98 different individuals. The frequency of individuals with differences in homologue peak positions varied among chromosome types: e.g., chromosome 15, 61%; chromosome 3, 4%. Flow karyotypes of 33 unrelated individuals were compared to determine the range of peak position among normal chromosomes. Chromosomes Y, 21, 22, 15, 16, 13, 14, and 19 were most heteromorphic, and chromosomes 2-8 and X were least heteromorphic. The largest chromosome 21 was 45% larger than the smallest 21 chromosome observed. The base composition of the variable regions differed among chromosome types. DNA contents of chromosome variants determined from flow karyotypes were closely correlated to measurements of DNA content made of gallocyanin chrome alum-stained metaphase chromosomes on slides. Fluorescence in situ hybridization with chromosome-specific repetitive sequences indicated that variability in their copy number is partly responsible for peak-position variability in some chromosomes. Heteromorphic chromosomes are identified for which parental flow karyotype information will be essential if de novo rearrangements resulting in small DNA content changes are to be detected with flow karyotyping.     

Application of flow karyotyping in prenatal detection of chromosome aberrations

This paper describes the application of bivariate flow karyotyping to (1) classification of chromosomes isolated from cultures of cells taken by amniocentesis and (2) detection of numerical and structural aberrations. Chromosomes were isolated from primary cultures 2-5 wk after amniocentesis, stained with Hoechst 33258 and chromomycin A3, and analyzed using dual beam flow cytometry. Information about chromosome DNA content and DNA base composition was derived from the locations of the peaks in the flow karyotypes, each peak being produced by one or more chromosome types with similar DNA content and DNA base composition. Information about the relative frequency of each chromosome type was determined on the basis of the relative volume of the peak for that chromosome type. Cytogenetic information determined on the basis of flow karyotypes was compared with that obtained by visual analysis following G-banding. Variability among the peak means and volumes in flow karyotypes was determined from analyses of 50 normal amniocyte cultures. Numerical aberrations involving chromosomes 21, 18, and Y were detected correctly in all of 28 analyses, including eight in a blind study. Structural aberrations involving chromosomes 1, 2, 3, 6, 9-12, 13, 14, 15, 21, and 22 were detected in all of seven cultures in a blind study. Flow karyotypes proved to be insensitive to small, normally occurring chromosome polymorphisms detected by banding analysis. In addition, a few samples were erroneously scored as having numerical aberrations.     

Background and peak evaluation of one parameter flow karyotypes on a PC/AT computer

Flow cytometry has become a fast, quantitative method for the classification of metaphase chromosomes in suspension (flow karyotyping) stained with fluorescent dyes. Such a flow karyotype (frequency distribution of the fluorescence signals) consists of several peaks. The peak pattern characterizes the analyzed chromosome complement. In many cases flow karyotypes contain a continuum of an unspecific background deriving from chromosome fragments or chromosome aggregates. For the quantitative evaluation of a flow karyotype this background has to be subtracted by a suitable background function. In this approach the application of chi 2-functions is described. The feasibility of this method to flow karyotypes has been concluded from a computer simulation of chromosome breaking under different conditions. In spite of the rather rough assumptions of the model compared to the complex reasons that influence chromosome breaking, the chi 2-function fits the background better than the exponential function in current use. The approximation of a Gaussian distribution function by the chi 2-function also makes it possible to use the same subtraction procedure for chromosome aggregates. The procedure was tested for isolated chromosomes of Chinese hamster cell lines under different states of breaking. For further evaluation of one parameter flow karyotypes a setup of computer routines has been developed for PC/AT and compatible computer systems. Different peak values of these flow karyotypes can be determined (e.g. peak mean, standard deviation, absolute and relative peak area etc.). The applied method is to fit Gaussian curves to each peak of an experimentally measured histogram by using an interactive program. Fluctuations depending on 'noise' may be suppressed by a 'k-nearest-neighbours' smoothing procedure.     

Flow cytometric analysis of the chromosomes and stability of a wheat cell-culture line

A rapidly growing, long-term suspension culture derived from Triticum aestivum L. (wheat) was synchronized using hydroxyurea and colchicine, and a chromosome suspension with chromosomes was made. After staining with the DNA-specific fluorochromes Hoechst 33258 and Chromomycin univariate and bivariate flow-cytometry histograms showed 15 clearly resolved peaks corresponding to individual chromosome types or groups of chromosomes with similar DNA contents. The flow karyotype was closely similar to a histogram of DNA content measurements of Feulgen-stained chromosomes made by microdensitometry. We were able to show the stability of the flow karyotype of the cell line over a year, while a parallel subculture had a slightly different, stable, karyotype following different growth conditions. The data indicate that flow cytometric analysis of plant karyotypes enables accurate, statistically precise chromosome classification and karyotyping of cereals. There was little overlap between individual flow-histogram peaks, so the method is useful for flow sorting and the construction of chromosome specific-recombinant DNA libraries. Using bivariate analysis, the AT:GC ratio of all the chromosomes was remarkably similar, in striking contrast to mammalian flow karyotypes. We speculate about a fundamental difference in organization and homogenization of DNA sequences between chromosomes within mammalian and plant genomes. Received: 24 April 1996 / Accepted: 24 May 1996     

Immunofluorescent labeling of centromeres for flow cytometric analysis.

A procedure to stain the centromeric region of chromosomes for dual beam flow cytometric analysis is described. Serum from a CREST (Scleroderma syndrome) patient presenting a high titer of anticentromeric antibodies was chosen on the basis of specificity of labeling of cells on slides. The high affinity of the antibodies to centromeres and low binding to chromosomal arms allowed the development of an indirect immunofluorescent labeling procedure using isolated and unfixed chromosomes stabilized by Mg++ ions. Discontinuous Ficoll gradients were used to separate chromosomes from unbound antibodies. With this procedure, chromosome clumping and degradation were minimal. The chromosomes were then stained with the DNA dyes Hoechst 33258 and chromomycin A3, before dual beam flow cytometric analysis. Flow karyotypes, with good chromosome peak resolution, were obtained for both human and hamster chromosomes subjected to the immunolabeling procedure. For quantification of FITC fluorescence due to bound antibody, chromosomes were counterstained with Hoechst only. The FITC intensity of antibody-labeled human and hamster chromosomes were 4-10 and 20 times greater than control chromosomes, respectively. These results suggest that the staining procedure may be suitable for immunolabeling of chromosomes with antibodies recognizing other nuclear proteins and their subsequent quantification by flow cytometry.     

Swine chromosomes: flow sorting and spot blot hybridization

Flow cytometry analysis was applied to swine chromosomes prepared from phytohemagglutinin (PHA) stimulated peripheral blood lymphocytes. Flow karyotypes from both sexes and from t(3;7) translocation carrier females were obtained. A certain number of chromosome pairs could be assigned to various peaks. In fact, 13 peaks were observed for 18 autosomal pairs plus X and Y. Moreover, abnormalities owing to the t(3;7) translocation were readily observable. The number of base pairs for chromosomes associated with the various peaks was estimated by comparison with human flow karyotypes. The following four peaks were thus sorted: the peak assumed to represent the translocated chromosome 7 plus the normals associated with it; the corresponding peak from a normal swine; the peak assumed to contain among others the normal chromosome 7; and finally the peak corresponding to swine chromosome 1. Chromosomes of each peak were collected on Pall Biodyne membrane. Following appropriate denaturation and prehybridization, the four samples were hybridized with a human leucocyte antigen (HLA) class I 32P-labelled cDNA probe, representing most of the coding sequence of the HLA B7 gene. The results confirmed previous data from other techniques that assigned the swine MHC(SLA) to chromosome 7. Subsequently, sorted samples were hybridized with a porcine genomic Interferon alpha probe in order to confirm the mapping of this gene family on porcine chromosome 1.     

Flow Cytometric Analysis and Chromosome Sorting of Barley (Hordeum vulgare L.)

Flow cytometric analysis was systematically performed to optimize the concentration and duration of hydroxyurea (DNA synthesis inhibitor) and trifluralin (metaphase blocking reagent) treatments for synchronizing the cell cycle and accumulating metaphase chromosomes in barley root tips. A high metaphase index (76.5% in the root tip meristematic area) was routinely achieved. Seedlings of about 1.0-cm length were treated with 1.25 mM hydroxyurea for 14 h to synchronize the root tip meristem cells at the S/G2 phase. After rinsing with hydroxyurea, the seedlings were incubated in a hydroxyurea-free solution for 2 h and were treated with 1 microM trifluralin for 4 h to accumulate mitotic cells in the metaphase. The consistent high metaphase index depended on the uniform germination of seeds prior to treatment. High-quality and high-quantity isolated metaphase chromosomes were suitable for flow cytometric analysis and sorting. Flow karyotypes of barley chromosomes were established via univariate and bivariate analysis. A variation of flow karyotypes was detected among barley lines. Two single chromosome types were identified and sorted. Bivariate analysis showed no variation among barley individual chromosomes in AT and GC content.     

Flow Cytometric Analysis and Chromosome Sorting of Barley (Hordeum vulgare L.)

Flow cytometric analysis was systematically performed to optimize the concentration and duration of hydroxyurea (DNA synthesis inhibitor) and trifluralin (metaphase blocking reagent) treatments for synchronizing the cell cycle and accumulating metaphase chromosomes in barley root tips. A high metaphase index (76.5% in the root tip meristematic area) was routinely achieved. Seedlings of about 1.0-cm length were treated with 1.25 mM hydroxyurea for 14 h to synchronize the root tip meristem cells at the S/G2 phase. After rinsing with hydroxyurea, the seedlings were incubated in a hydroxyurea-free solution for 2 h and were treated with 1 M trifluralin for 4 h to accumulate mitotic cells in the metaphase. The consistent high metaphase index depended on the uniform germination of seeds prior to treatment. High-quality and high-quantity isolated metaphase chromosomes were suitable for flow cytometric analysis and sorting. Flow karyotypes of barley chromosomes were established via univariate and bivariate analysis. A variation of flow karyotypes was detected among barley lines. Two single chromosome types were identified and sorted. Bivariate analysis showed no variation among barley individual chromosomes in AT and GC content.     

Bivariate flow karyotyping of acute myelocytic leukemia in the BNML rat model

Univariate as well as bivariate flow karyotyping has been performed on chromosome suspensions obtained from the Brown Norway myelocytic leukemia (BNML), a rat model for human acute myelocytic leukemia (AML). Flow karyograms were obtained from both the in vivo transplantable parent line and from an in vitro established cell line. Density gradient centrifugation performed on cells arrested in mitosis resulted in an enrichment of mitotic cells. Furthermore, with this procedure leukemic and nonleukemic cells could be separated. Univariate analysis with propididum iodide (PI) as a DNA stain revealed the position of the several tumor-specific marker chromosomes in the in vitro cell line. Estimations of the peak position of the various chromosomes was done by comparing the univariate flow karyogram with a computer-simulated karyogram from the BNML that was derived from the mean length of the individual chromosomes in conventionally prepared metaphase slides. By comparing the bivariate flow karyogram of the in vivo BNML cells with the flow karyogram of normal BN cells, it was clearly demonstrated which peaks are involved in the altered chromosomal pattern of the BNML. No differences were found between the flow karyograms of the in vitro- and the ex vivo-derived chromosome suspensions in this rat leukemia model.     

Swine chromosomal DNA quantification by bivariate flow karyotyping and karyotype interpretation.

Human and swine chromosomes were analyzed separately and as a mix to obtain bivariate flow karyotypes. They were normalized to each other in order to use the human chromosomal DNA content as standard. Our results led to the characterization of the "DNA line" in swine identical to the human "DNA line." Estimation of the DNA content in mega-base pairs of the swine chromosomes is proposed. Chromosomal assignment to the various resolved peaks on the bivariate swine flow karyotype is suggested from the relation between DNA content quantified by flow cytometry and chromosomal size. Swine chromosomes 1, 13, 6, 5, 10, 16, 11, 18, and Y were assigned to peaks A, B, C, K, L, N, O, Q, and Y, respectively. Peaks D and E were assumed to contain chromosomes 2 and 14, but without specific assignment. Similarly, P and M peaks were expected to correspond to chromosomes 12 and 17. Of the remaining chromosomes (3, 7, X, 8, 15, 9, and 4), chromosomes 3, 7, and X, which were assigned previously to peaks F, G, and H, respectively, led us to deduce that chromosomes 15 and 8 belonged to peaks I and J, and chromosomes 9, 4, and X to peak H.     

Complete homology maps of the rabbit (Oryctolagus cuniculus) and human by reciprocal chromosome painting

Fluorescence in situ hybridization (FISH) was used to construct a homology map to analyse the extent of evolutionary conservation of chromosome segments between human and rabbit (Oryctolagus cuniculus, 2n = 44). Chromosome-specific probes were established by bivariate fluorescence activated flow sorting followed by degenerate oligonucleotide-primed PCR (DOP-PCR). Painting of rabbit probes to human chromosomes and vice versa allowed a detailed analysis of the homology between these species. All rabbit chromosome paints, except for the Y paint, hybridized to human chromosomes. All human chromosome paints, except for the Y paint, hybridized to rabbit chromosomes. The results obtained revealed extensive genome conservation between the two species. Rabbit chromosomes 12, 19 and X were found to be completely homologous to human chromosomes 6, 17 and X, respectively. All other human chromosomes were homologous to two or sometimes three rabbit chromosomes. Many conserved chromosome segments found previously in other mammals (e.g. cat, pig, cattle, Indian muntjac) were also found to be conserved in rabbit chromosomes.     

Preparation and bivariate analysis of suspensions of human chromosomes

Chromosomes were isolated from a variety of human cell types using a HEPES-buffered hypotonic solution (pH 8.0) containing KCl, MgSO4, dithioerythritol, and RNase. The chromosomes isolated by this procedure could be stained with a variety of fluorescent stains including propidium iodide, chromomycin A3, and Hoechst 33258. Addition of sodium citrate to the stained chromosomes was found to improve the total fluorescence resolution. High-quality bivariate Hoechst vs. chromomycin fluorescence distributions were obtained for chromosomes isolated from a human fibroblast cell strain, a human colon carcinoma cell line, and human peripheral blood lymphocyte cultures. Good flow karyotypes were also obtained from primary amniotic cell cultures. The Hoechst vs. chromomycin flow karyotypes of a given cell line, made at different times and at dye concentrations varying over fourfold ranges, show little variation in the relative peak positions of the chromosomes. The size of the DNA in chromosomes isolated using this procedure ranges from 20 to over 50 kilobases. The described isolation procedure is simple, it yields high-quality flow karyotypes, and it can be used to prepare chromosomes from clinical samples.     

Anti-kinetochore staining for single laser, bivariate flow sorting of Indian muntjac chromosomes.

Flow cytometric chromosome sorting typically relies upon dual-laser, bivariate analysis after staining with two different base pair-specific dyes for resolution of chromosomes with similar DNA content. The availability of FITC-conjugated antibodies offers the possibility of single-laser bivariate analysis when combined with propidium iodide (PI) DNA staining, but requires exploitable antigenic differences between chromosomes of interest. A technique was developed for indirect immunofluorescent anti-kinetochore staining of Indian muntjac chromosomes in suspension. Primary antibody binding within permeabilized whole cells minimized centrifugation-induced loss of chromosomal integrity. Subsequent FITC-conjugated second antibody binding was not affected by concurrent PI-counterstaining. Anti-kinetochore staining facilitated resolution of chromosomes No. 2 and X, which formed a doublet peak upon univariate DNA content analysis, as well as recognition of the Y2 peak which was indistinguishable from debris by univariate analysis. The technique allowed greater than 90% purification of each Indian muntjac chromosome.     

Quantification of the DNA content of structurally abnormal X chromosomes and X chromosome aneuploidy using high resolution bivariate flow karyotyping

Quantification of the Hoechst and chromomycin A3 fluorescence intensities of mitotic human chromosomes isolated from karyotypically normal and abnormal cells was performed with a dual beam flow cytometer. The resultant flow karyotypes contain information about the relative DNA content and base composition of chromosomes and their relative frequencies in the mitotic cell sample. The relative copy number of X and Y chromosomes was determined for 38 normal males and females and 6 cell lines with X or Y chromosome aneuploidy. Flow karyotype diagnoses corresponded with conventional cytogenetic results in all cases. We show that chromosome DNA content can be derived from peak position in Hoechst vs. chromomycin flow karyotypes. These values are linearly related to propidium iodide staining intensity as measured with flow cytometry and to the binding of gallocyanin chrome alum to phosphate groups as measured with slide-based scanning photometry. Cell lines with deleted or dicentric X chromosomes ranging in length from 0.53 to 1.95 times normal were analyzed by using flow cytometry. The measured difference in DNA content between a normal X and each of the structurally abnormal chromosomes was linearly correlated to the difference predicted from cytogenetics and/or probe analyses. Deletions of 3-5 Mb, which were at and below the detection limits of conventional cytogenetics, could be quantified by flow karyotyping in individuals with X-linked diseases such as Duchenne muscular dystrophy, choroideremia, and ocular albinism/ichthyosis. The results show that the use of flow karyotyping to quantify the size of restricted regions of the genome can complement conventional cytogenetics and other physical mapping techniques in the study of genetic disorders.     

Preparation and flow cytometric analysis of metaphase chromosomes of tomato

Summary A procedure for the preparation of tomato chromosome suspensions suitable for flow cytometric analysis is described. Rapidly growing cell suspension cultures of Lycopersicon esculentum cv VFNT cherry and L. pennellii LA716 were treated with colchicine to enrich for metaphase chromosomes. Metaphase indices between 20 and 35% were routinely obtained when cultures were exposed to 0.1% colchicine for 15–18 h after 2 days of subculture. Mitotic cells were isolated by brief treatment with cell wall digesting enzymes in a medium with low osmolarity (325 mOsm/kg of H₅₂O). The low osmolarity medium was needed to avoid the chromosome clumping and decondensation seen in standard media. Suspensions of intact chromosomes were prepared by lysing swollen protoplasts in various buffers (MgSO₄, polyamines, hexylene glycol, or KCl-propidium iodide) similar in contents to the buffers used to isolate mammalian chromosomes. For univariate flow cytometric analysis, chromosome suspensions were stained with a fluorescent DNA-binding stain (propidium iodide, Hoechst 33258, mithramycin, or chromomycin A3) and analyzed using an EPICS flow cytometer (Profile Analyzer or 753). Peaks for the chromosomes, chromatids, clumps of chromosomes, nuclei, and fluorescent debris were seen on a histogram of log of fluorescence intensity, and were confirmed by microscopic examination of the objects collected by flow-sorting. Chromosome suspensions prepared in MgSO₄ buffer have the highest frequency of intact chromosomes and the least fluorescent cellular debris. Peaks similar to theoretical univariate flow karyotypes of tomato chromosomes were seen on the observed univariate flow karyotypes, but were not as well resolved. Bivariate flow analysis of tomato chromosome suspension using double-stain combination, Hoechst 33258 and chromomycin A3, and two laser beams showed better resolution of some chromosomes.