Analysis of DNA distributions from flow cytometry

A recently proposed automatic procedure for analyzing DNA distribution from flow cytometric data is extensively tested against simulated data. After a discussion of the procedure itself and of the simulation program, the results obtained are reported. They are evidence of the reliability of the procedure in extracting the proper underlying DNA distribution from sets of data obtained under various simulated instances. The different sources of error are then analyzed, along with their quantitative effects on the fit of the fluorescence histogram.     

Application of method of moments analysis to fluorescence decay lifetime distributions

In fluorescence decay work, distributions of exponential decay lifetimes are anticipated where complex systems are examined. We describe here methods of gaining information on such distributions using the method of moments analysis approach. The information obtained may be as simple as the average and deviation of the lifetime distribution, quantities which we show may be estimated directly from the results of a multiexponential analysis. An approximation to the actual distribution shape may also be obtained using a procedure we call the variable filter analysis (VFA) method without making any assumptions about the shape of the distribution. Tests of VFA using both simulated and experimental data are described. Limitations of this method and of distribution analysis methods in general are discussed. Results of analyses on experimental decays for ethidium intercalated in core particles and in free DNA are reported.     

Multimerization-cyclization of DNA fragments as a method of conformational analysis

Ligation of short DNA fragments results in the formation of linear and circular multimers of various lengths. The distribution of products in such a reaction is often used to evaluate fragment bending caused by specific chemical modification, by bound ligands or by the presence of irregular structural elements. We have developed a more rigorous quantitative approach to the analysis of such experimental data based on determination of j-factors for different multimers from the distribution of the reaction products. j-Factors define the effective concentration of one end of a linear chain in the vicinity of the other end. To extract j-factors we assumed that kinetics of the reaction is described by a system of differential equations where j-factors appear as coefficients. The assumption was confirmed by comparison with experimental data obtained here for DNA fragments containing A-tracts. At the second step of the analysis j-factors are used to determine conformational parameters of DNA fragments: the equilibrium bend angle, the bending rigidity of the fragment axis, and the total twist of the fragments. This procedure is based on empirical equations that connect the conformational parameters with the set of j-factors. To obtain the equations, we computed j-factors for a large array of conformational parameters that describe model fragments. The approach was tested on both simulated and actual experimental data for DNA fragments containing A-tracts. A-tract DNA bend angle determined here is in good agreement with previously published data. We have established a set of experimental conditions necessary for the data analysis to be successful.     

Test statistics for detecting aneuploidy and hyperdiploidy

Possible approaches to the analytical evaluation of ploidy patterns are discussed and two specific problems are considered: detection of early onset of aneuploidy and detection of moderate hyperdiploidy. A statistical model for a euploid DNA pattern is formulated in terms of a mixture distribution. A test statistic for detecting deviations from this pattern is defined, and its performance is evaluated for simulated data representing differing degrees of severity of aneuploidy. An analysis based on a discriminant function using order statistics of the sample cumulative distribution functions is proposed for detecting hyperdiploidy. This procedure has the advantage of being relatively distribution-free; its performance is evaluated for simulated data and is compared with that of its classical counterparts. Although the results reported are only preliminary, they indicate that tailor-made statistical analyses can provide early detection of aneuploidy and hyperdiploidy with known and acceptable error rates using clinically reasonable sample sizes.     

Atomic force microscopy imaging of double stranded DNA and RNA.

A procedure for imaging long DNA and double stranded RNA (dsRNA) molecules using Atomic Force Microscopy (AFM) is described. Stable binding of double stranded DNA molecules to the flat mica surface is achieved by chemical modification of freshly cleaved mica under mild conditions with 3-aminopropyltriethoxy silane. We have obtained striking images of intact lambda DNA, Hind III restriction fragments of lambda DNA and dsRNA from reovirus. These images are stable under repeated scanning and measured contour lengths are accurate to within a few percent. This procedure leads to strong DNA attachment, allowing imaging under water. The widths of the DNA images lie in the range of 20 to 80nm for data obtained in air with commercially available probes. The work demonstrates that AFM is now a routine tool for simple measurements such as a length distribution. Improvement of substrate and sample preparation methods are needed to achieve yet higher resolution.     

Atomic Force Microscopy Imaging of Double Stranded DNA and RNA

Abstract

A procedure for imaging long DNA and double stranded RNA (dsRNA) molecules using Atomic Force Microscopy (AFM) is described. Stable binding of double stranded DNA molecules to the flat mica surface is achieved by chemical modification of freshly cleaved mica under mild conditions with 3-aminopropyltriethoxy silane. We have obtained striking images of intact lambda DNA, Hind III restriction fragments of lambda DNA and dsRNA from reovirus. These images are stable under repeated scanning and measured contour lengths are accurate to within a few percent. This procedure leads to strong DNA attachment, allowing imaging under water. The widths of the DNA images lie in the range of 20 to 80nm for data obtained in air with commercially available probes. The work demonstrates that AFM is now a routine tool for simple measurements such as a length distribution. Improvement of substrate and sample preparation methods are needed to achieve yet higher resolution.     

Interpretation of Sucrose Gradient Sedimentation Pattern of Deoxyribonucleic Acid Fragments Resulting from Random Breaks

Mass distribution in a sucrose gradient of deoxyribonucleic acid (DNA) fragments arising as a result of random breaks is predicted by analytical means from which computer evaluations are plotted. The analytical results are compared with the results of verifying experiments: (i) a Monte Carlo computer experiment in which simulated molecules of DNA were individuals of unit length subjected to random "breaks" applied by a random number generator, and (ii) an in vitro experiment in which molecules of T4 DNA, highly labeled with (32)P, were stored in liquid nitrogen for variable periods of time during which a precisely known number of (32)P atoms decayed, causing single-stranded breaks. The distribution of sizes of the resulting fragments was measured in an alkaline sucrose gradient. The profiles obtained in this fashion were compared with the mathematical predictions. Both experiments agree with the analytical approach and thus permit the use of the graphs obtained from the latter as a means of determining the average number of random breaks in DNA from distributions obtained experimentally in a sucrose gradient. An example of the application of this procedure to a previously unresolved problem is provided in the case of DNA from ultraviolet-irradiated phage which undergoes a dose-dependent intracellular breakdown. The relationship between the number of lethal hits and the number of single-stranded breaks was not previously established. A comparison of the calculated number of nicks per strand of DNA with the known dose in phage-lethal hits reveals a relationship closely approximating one lethal hit to one single-stranded break.     

Simulation of NOESY spectra of DNA segments: A new scaling procedure for iterative comparison of calculated and experimental NOE intensities

Summary A new algorithm for simulation of two-dimensional NOESY spectra of DNA segments has been developed. For any given structure, NOE intensities are calculated using the relaxation matrix approach and a new realistic procedure is suggested for 1:1 comparison of calculated and experimental intensities. The procedure involves a novel method for scaling of calculated NOE intensities to represent volumes of digitised cross peaks in NOESY spectra. A data base of fine structures of all the relevant cross peaks with Lorentzian line shapes and in-phase components, is generated in a digitised manner by two-dimensional Fourier transformation of simulated time domain data, assuming a total intensity of 1.0 for each of the cross peaks. With this procedure, it is shown that the integrated volumes of these digitised cross peaks above any given threshold scale exactly as the total intensity of the respective peaks. This procedure eliminates the repetitive generation of digitised cross peaks by two-dimensional Fourier transformation during the iterative process of structure alteration and NOE intensity calculation and thus enhances the speed of DNA structure optimization. Illustrative fits of experimental and calculated spectra obtained using the new procedure are shown.[/p]     

Averaged gene expressions for regression

Although averaging is a simple technique, it plays an important role in reducing variance. We use this essential property of averaging in regression of the DNA microarray data, which poses the challenge of having far more features than samples. In this paper, we introduce a two-step procedure that combines (1) hierarchical clustering and (2) Lasso. By averaging the genes within the clusters obtained from hierarchical clustering, we define supergenes and use them to fit regression models, thereby attaining concise interpretation and accuracy. Our methods are supported with theoretical justifications and demonstrated on simulated and real data sets.     

Point mutations as an optimal search process in biological evolution

Point mutations are pictured as jumps in a phase space representing the sequences of amino acids or nucleotides as discrete points. It is shown that this space can be given a natural metric by quantifying common physical and chemical properties of amino acid constituents in terms of a natural measure. Evolution through point mutations is simulated by the search for points in the phase space representing amino acid sequences of high survival fitness. Due to the local compactness of the distribution of these functionally allowed points in phase space any successful search procedure has characteristics qualitatively different from those in the case of a random distribution. This is demonstrated by model calculations. A specified distribution of allowed points is generated with subsequent evaluation of the success of the retrieval process as a function of the jump probabilities between lattice sites. The results of such simulations are compared with data obtained from the analysis of the DNA or mRNA sequences coding related proteins. By counting silent and expressed nucleotide replacement frequencies one can draw conclusions as to the efficacy of the natural evolutionary search processes in the phase space of amino acid sequences. There are cases, where the highest possible information gain of one bit per accepted point mutation is achieved. In general the information gain is found to be somewhat sub-maximal due to functional requirements.     

DNA assessment in thin sections of lymphomas by densitometric and stereological methods: comparison with imprint and flow cytometric results

The nuclear DNA content was estimated in 2 microns sections of 18 lymphoma cases by two methods: (1) Feulgen densitometry using QTM 900 with correction by Bins' procedure which allows size-independent DNA distributions; (2) stereological unfolding as proposed by Cruz-Orive giving sphere-size distributions. A general correlation was found between results and DNA measurements obtained by imprint and flow cytometric techniques in the same specimens. When histologic DNA profiles were compared to cytologic histograms, a high correlation was found between the distribution of ploidy classes by correspondence analysis. However two highly proliferating lymphomas were erroneously classified as aneuploid. Conversely, sphere-size distributions allowed the identification of the majority of aneuploid lymphomas but failed to recognize proliferating ones. It appears that when cytologic specimens are not available, densitometric studies on sections may provide valuable information on DNA content, with complementary data obtained from stereological procedures.     

Analysis of DNA synthesis rate of cultured cells from flow cytometric data

The rate of DNA synthesis along S phase is estimated from flow cytometric histograms on the basis of a mathematical model of a cell population. In the absence of loss, the model expresses the population kinetics in terms of DNA synthesis rate, S-phase influx, and population size. A single histogram is sufficient to determine the DNA synthesis rate when the population is in balanced exponential growth. Two suitably chosen histograms are necessary if the S-phase influx is exponential in a time interval longer than the S-phase duration. The analysis procedure was tested on published autoradiographic data and applied to three cultured cell lines (CM-S, 3LL, and M14 cells) that show various patterns of DNA distribution. In each case the cell-cycle fractions, the DNA synthesis rate, and the S-phase duration were obtained.     

Turn of promotor DNA by cAMP receptor protein characterized by bead model simulation of rotational diffusion

The rotation diffusion of DNA double helices and their complexes with the cAMP receptor protein (CRP) has been simulated by bead models, in order to derive information on their structure in solution by comparison with results obtained from dichroism decay measurements. Straight DNA double helices are simulated by linear, rigid strings of overlapping beads. The radius of the beads and the length of the string are increased simultaneously by the same increments from initial outer dimensions derived from crystallographic data to final values, which are fitted to experimental rotation time constants observed for short DNA fragments (less than 100 bp). The final values reflect the solvated structure with the same 'solvation layer' added in all three dimensions. The protein is simulated by overlapping beads, which are assembled to a structure very similar to that found by x-ray crystallography. Complexes of the protein with DNA are formed with the centres of palindromic DNA sites at the centre of the two helix-turn-helix-motifs of the protein with some overlap of the two components. Simulation of the experimental data obtained for CRP complexes with specific DNA in the presence of cAMP requires strong bending of the double helices. According to our simulation the DNA is almost completely wrapped around the protein both in the complexes with a 62 bp fragment containing the standard CRP site and with a 80 bp fragment containing the second binding site of the lac operon. Simulations of the data obtained for a 203 bp fragment with both binding sites suggest that the two bound CRP proteins are in contact with each other and that the DNA is wrapped around the two protein dimers. A stereochemical model is suggested with a tetrahedral arrangement of the four protein subunits, which provides the advantage that two binding sites of the protein formed by two subunits each are located favorable for tight contacts to two binding sites on bent DNA, provided that the DNA sites are separated by an integer number of helix turns. In summary, the simulations demonstrate strong bending, which can be reflected by an arc radius in the range around 50 A. According to these data the overall bending angle of our longest DNA fragment is approximately 180 degrees, and thus the protruding ends are sufficiently close to each other such that RNA polymerase, for example, could contact both helical segments.     

Turn of Promotor DNA by cAMP Receptor Protein Characterized by Bead Model Simulation of Rotational Diffusion

Abstract

The rotation diffusion of DNA double helices and their complexes with the cAMP receptor protein (CRP) has been simulated by bead models, in order to derive information on their structure in solution by comparison with results obtained from dichroism decay measurements. Straight DNA double helices are simulated by linear, rigid strings of overlapping beads. The radius of the beads and the length of the string are increased simultaneously by the same increments from initial outer dimensions derived from crystallographic data to final values, which are fitted to experimental rotation time constants observed for short DNA fragments (< 100 bp). The final values reflect the solvated structure with the same ‘solvation layer’ added in all three dimensions. The protein is simulated by overlapping beads, which are assembled to a structure very similar to that found by x-ray crystallography. Complexes of the protein with DNA are formed with the centres of palindromic DNA sites at the centre of the two helix- turn-helix-motifs of the protein with some overlap of the two components. Simulation of the experimental data obtained for CRP complexes with specific DNA in the presence of cAMP requires strong bending of the double helices. According to our simulation the DNA is almost completely wrapped around the protein both in the complexes with a 62 bp fragment containing the standard CRP site and with a 80 bp fragment containing the second binding site of the lac operon. Simulations of the data obtained for a 203 bp fragment with both binding sites suggest that the two bound CRP proteins are in contact with each other and that the DNA is wrapped around the two protein dimers. A stereochemical model is suggested with a tetrahedral arrangement of the four protein subunits, which provides the advantage that two binding sites of the protein formed by two subunits each are located favorable for tight contacts to two binding sites on bent DNA provided that the DNA sites are separated by an integer number of helix turns. In summary, the simulations demonstrate strong bending, which can be reflected by an arc radius in the range around 50 Å. According to these data the overall bending angle of our longest DNA fragment is approximately 180°, and thus the protruding ends are sufficiently close to each other such that RNA polymerase, for example, could contact both helical segments.     

Photometric determination of the DNA distribution in the 24 human chromosomes

From five normal individuals the DNA content and the DNA arm ratios of the 24 metaphase chromosomes were determined by means of scanning densitometry of photographic negatives of Feulgen-stained metaphase preparations. The results showed high reproducibility of the measuring procedure. The obtained DNA values for the 24 chromosomes showed general correspondence between the individuals. No differences between males and females were found. The DNA arm ratios showed somewhat higher inter-individual variability, especially for the acrocentric chromosomes. Our data are in agreement with other data published so far, which were obtained with somewhat different techniques, indicating that the DNA content of the individual human chromosomes in general is highly constant. Attempts were made to distinguish chromosomes by their DNA content and DNA ratio. It appears that classification of chromosomes using these parameters cannot compete with classification according to the banding patterns. Determination of the total DNA content and DNA distribution along the metaphase chromosomes may, however, provide a frame of reference for cytochemical methods directed towards the localization and quantification of molecular properties of the chromosome.     

Computation of identity by descent probabilities conditional on DNA markers via a Monte Carlo Markov Chain method

The accurate estimation of the probability of identity by descent (IBD) at loci or genome positions of interest is paramount to the genetic study of quantitative and disease resistance traits. We present a Monte Carlo Markov Chain method to compute IBD probabilities between individuals conditional on DNA markers and on pedigree information. The IBDs can be obtained in a completely general pedigree at any genome position of interest, and all marker and pedigree information available is used. The method can be split into two steps at each iteration. First, phases are sampled using current genotypic configurations of relatives and second, crossover events are simulated conditional on phases. Internal track is kept of all founder origins and crossovers such that the IBD probabilities averaged over replicates are rapidly obtained. We illustrate the method with some examples. First, we show that all pedigree information should be used to obtain line origin probabilities in F2 crosses. Second, the distribution of genetic relationships between half and full sibs is analysed in both simulated data and in real data from an F2 cross in pigs.     

Estimation of Additive Genetic Variance Components in Aquaculture Populations Selectively Pedigreed by DNA Fingerprinting

A method to estimate genetic variance components in populations partially pedigreed by DNA fingerprinting is presented. The focus is on aquaculture, where breeding procedures may produce thousands of individuals. In aquaculture populations the individuals available for measurement will often be selected, i.e. will come from the upper tail of a size‐at‐age distribution, or the lower tail of an age‐at‐maturity distribution etc. Selection typically occurs by size grading during grow‐out and/or choice of superior fish as broodstock. The method presented in this paper enables us to estimate genetic variance components when only a small proportion of individuals, those with extreme phenotypes, have been identified by DNA fingerprinting. We replace the usual normal density by appropriate robust least favourable densities to ensure the robustness of our estimates. Standard analysis of variance or maximum likelihood estimation cannot be used when only the extreme progeny have been pedigreed because of the biased nature of the estimates. In our model‐based procedure a full robust likelihood function is defined, in which the missing information about non‐extreme progeny has been taken into account. This robust likelihood function is transformed into a computable function which is maximized to get the estimates. The estimates of sire and dam additive variance components are significantly and uniformly more accurate than those obtained by any of the standard methods when tested on simulated population data and have desirable robustness properties.     

A quantitative method for evaluating bivariate flow cytometric data obtained using monoclonal antibodies to bromodeoxyuridine.

A method is presented for analyzing data from bivariate analysis of cell populations exposed to bromodeoxyuridine and subsequently examined both for the presence of BrdUrd and for the cellular DNA content. It is shown that certain features may be defined in the bivariate data which are constant independent both of cell type and, within limits, experimental variability. These landmark features include the ratio of red, DNA, fluorescence of G2 + M cells to G1 cells, the ratio of green fluorescence corresponding to the non-specific binding of unlabeled G2 + M cells to unlabeled G1 cells, and the distribution of green fluorescence in unlabeled cells. The landmarks make it possible to standardize rules for establishing the separation line between-labeled and unlabeled cells as required in these experiments to obtain estimates of cytokinetic parameters. Values obtained for the DNA synthesis time and the potential doubling time which result from different decision rules for distinguishing labeled from unlabeled are compared in two murine tumor lines. The potential doubling time, but not the DNA synthesis time is shown to depend sensitively on the separation line. Suggestions are presented for analyzing clinical data with this procedure.     

A rapid single step staining technique for DNA analysis by flow microfluorimetry

A new procedure is reported for the staining of DNA, for flow microfluorimetry. It allows the production of stained cell nuclei in a single step by incorporating the DNA stain with a solution of the nonionic detergent Triton-X-100. This method has been found to be applicable to all DNA fluorochromes tested (ethidium bromide, propidium iodide, mithramycin, DAPI, Hoechst 33342). DNA histograms obtained in this way are comparable to those using conventional staining techniques, e.g., ethanol fixation followed by staining. Using this procedure the DNA content distribution of solid tissue or cells from suspension or monolayer cultures can be generated in less than 5 min.