Evaluation of number average length analysis in quantifying double strand breaks in genomic DNAs

Double strand breaks in DNA can be quantified down to very low frequencies (a few per Gigabase pair) in nanogram quantities of nonradioactive, genomic DNA by dispersing the DNAs on electrophoretic gels, digitizing them by quantitative electronic imaging, and calculating the DNA lengths by number average length analysis. No specific distribution of damages is required for number average length analysis. To test the validity of this approach, we used DNA populations of known absolute lengths and break frequencies as experimental DNAs and calculated the number average lengths and double strand break levels. Experimental DNAs and length standards were dispersed using pulsed field electrophoretic modes (unidirectional pulsed field, contour clamped homogeneous field, or transverse alternating field) appropriate for their size range, stained with ethidium, destained, and a quantitative electronic image obtained. A dispersion curve was constructed from the migration-mobility relationships of the length standard DNAs, and the number average lengths of the experimental DNAs were calculated. The calculated DNA lengths agreed well with the actual lengths. Furthermore, the double strand break frequencies calculated through number average length analysis of DNAs dispersed by these pulsed field gel modes and digitized by quantitative electronic imaging were in excellent agreement with the actual values for populations of DNA over the size range of approximately 4 kbp to approximately 3 Mbp. The use of this approach in quantifying DNA damages is illustrated for double strand breaks and damage clusters (e.g., OxyPurine clusters recognized by Escherichia coli Fpg protein) induced in T7 DNA by ionizing radiation.     

Characterization of long and short repetitive sequences in the sea urchin genome

Long and short repetitive sequences were purified from the DNA of Paracentrotus lividus under conditions designed to optimize the yield of complete, end to end sequences. Double-stranded long repeat DNA prepared in this manner ranged in length from approximately 3000 to 15 000 nucleotide pairs with average sizes of approximately 6000 base pairs. In the electron microscope, long repeat DNA was observed to possess continuous sequences that often appeared to be terminated by one or more loops and/or fold backs. Long repeat DNA sequences, resheared to 300 base pairs, were found to have an average melting point identical to that for sheared native DNA. Thus, the reassociated duplexes of long repetitive DNA seem to possess very few mismatched base pairs. Reassociation kinetic analyses indicate that the majority of the long repeat sequences are reiterated only 4--7 times per haploid amount of DNA. Melt-reassociation analyses of short repetitive DNA, at several criteria, support the previously held concept that these sequences belong the sets or families of sequences which are inexact copies of one another. Our studies also support hypotheses suggesting that short repetitive sequences belong to families which may have arisen via distinct salttatory events. The relationships between long and short repetitive DNA sequences are considered with respect to widely held concepts of their sequence organization, evolution, and possible functions within eucaryotic genomes. A model for the possible organization of short repeats within long repetitive DNA sequences is also presented.     

Amplified DNA in Streptomyces fradiae.

A spontaneous mutant of Streptomyces fradiae contained an amplifiable unit of DNA with a sequence length of approximately 10.5 kilobases that was amplified to approximately 500 copies per chromosome. The amplified DNA appears to be cryptic. SalI fragments of the amplified DNA were cloned into Escherichia coli to construct a restriction map and characterize the amplified DNA. The amplified DNA contained tandem repeats of the amplifiable unit of DNA. The unit had an average base composition of 71% guanine plus cytosine, similar to the chromosomal DNA of Streptomyces species. At least a portion of the amplifiable unit of DNA was present at a low copy number in the wild-type strain. The phenotype of amplified DNA was designated Ads1SF for amplified DNA sequence 1 in S. fradiae.     

Formation of rings from segments of HeLa-cell nuclear deoxyribonucleic acid

Duplex segments of HeLa-cell nuclear DNA were generated by cleavage with DNA restriction endonuclease from Haemophilus influenzae. About 20-25% of the DNA segments produced, when partly degraded with exonuclease III and annealed, were found to form rings visible in the electron microscope. A further 5% of the DNA segments formed structures that were branched in configuration. Similar structures were generated from HeLa-cell DNA, without prior treatment with restriction endonuclease, when the complementary polynucleotide chains were exposed by exonuclease III action at single-chain nicks. After exposure of an average single-chain length of 1400 nucleotides per terminus at nicks in HeLa-cell DNA by exonuclease III, followed by annealing, the physical length of ring closures was estimated and found to be 0.02-0.1mum, or 50-300 base pairs. An almost identical distribution of lengths was recorded for the regions of complementary base sequence responsible for branch formation. It is proposed that most of the rings and branches are formed from classes of reiterated base sequence with an average length of 180 base pairs arranged intermittenly in HeLa-cell DNA. From the rate of formation of branched structures when HeLa-cell DNA segments were heat-denatured and annealed, it is estimated that the reiterated sequences are in families containing approximately 2400-24000 copies.     

A defined structure of the 30 nm chromatin fibre which accommodates different nucleosomal repeat lengths

Earlier work on the condensation of chromatins of different repeat lengths into the 30 nm fibre has been surveyed and it is shown that the external geometry of the fibre must be the same for all the chromatins. This can only be fitted by a helical coiling of nucleosomes into a solenoid with the linker DNA disposed internally. On this basis, various models were calculated and compared with published electric dichroism data. The only good fit is found with a 'reverse-loop' model, where the linker DNA forms a complete turn into the hole of the solenoid, of opposite hand to the nucleosomal DNA superhelix. This gives a topological linking number of one per nucleosome and would resolve the 'linking number paradox' if the DNA screw is the same in chromatin as in solution. The feasibility of a reverse-loop for short linkers (down to 15 base pairs) was investigated by model building and kinks of approximately 120 degrees into both DNA grooves are described, which will allow such packing. There will, however, be a 'forbidden' range for the linker DNA length, between approximately 1 and 14 bp, corresponding to nucleosomal repeats of 163 and 176 bp.     

The Estimation of the Number and the Length Distribution of Gene Conversion Tracts from Population DNA Sequence Data

DNA sequence variation studies report the transfer of small segments of DNA among different sequences caused by gene conversion events. Here, we provide an algorithm to detect gene conversion tracts and a statistical model to estimate the number and the length distribution of conversion tracts for population DNA sequence data. Two length distributions are defined in the model: (1) that of the observed tract lengths and (2) that of the true tract lengths. If the latter follows a geometric distribution, the relationship between both distributions depends on two basic parameters: ψ, which measures the probability of detecting a converted site, and , the parameter of the geometric distribution, from which the average true tract length, 1/(1 - ), can be estimated. Expressions are provided for estimating by the method of the moments and that of the maximum likelihood. The robustness of the model is examined by computer simulation. The present methods have been applied to the published rp49 sequences of Drosophila subobscura. Maximum likelihood estimate of for this data set is 0.9918, which represents an average conversion tract length of 122 bp. Only a small percentage of extant conversion events is detected.     

Resonant and localized breathing modes in terminal regions of the DNA double helix.

A Green's function approach is used in constructing a dynamic model of a semi-infinite length of the DNA homopolymer B poly(d) . poly(d). Considerable attention is focused on the hydrogen bond stretching close to the terminus. A melting (or breathing) coordinate (M) is defined as an average over the three linking hydrogen bond stretches in a unit cell. The thermal mean squared amplitude of (M) is enhanced at the chain end compared with the interior. Spectral branches at 69, 80 and 105 cm-1, as well as a local mode at 75 cm-1, are primary contributors to the enhancement. We suggest that this fact can affect the thermal melting of a DNA double helical homopolymer, enhancing the tendency to start from an end (if one is available). We show how certain infinite chain modes with small (M) amplitude can turn into breathing modes near the terminus, and suggest that the same phenomenon may occur near other specific base-pair sequences. There is also considerable attention paid to the low microwave region from approximately 0 to 1.75 cm-1. The thermally activated modes in this frequency region contribute approximately (0.02 A)2 to [M2(0)] at 40 K, approximately two orders of magnitude greater than for [M2(infinity)]. Most important however, is the existence of narrow resonant modes in this frequency region. Particularly pronounced resonances near 0.03 cm-1 and 0.08 cm-1 (approximately 0.9 and 2.4 GHz) amplify M2(0) at the terminus by about for orders of magnitude over the infinite chain value M2(infinity).     

Structure of the gene for human von Willebrand factor

von Willebrand factor is a large multimeric plasma protein composed of identical subunits which contain four types of repeated domains. von Willebrand factor is essential for normal hemostasis, and deficiency of von Willebrand factor is the most common inherited bleeding disorder of man. Four human genomic DNA cosmid libraries and one bacteriophage lambda library were screened with von Willebrand factor cDNA probes. Twenty positive overlapping clones were characterized that span the entire von Willebrand factor gene. A high-resolution restriction map was constructed for approximately 75% of the locus and a total of approximately 33.8 kilobases was sequenced on both strands including all intron-exon boundaries. The gene is approximately 178 kilobases in length and contains 52 exons. The exons vary from 40 to 1379 base pairs in length, and the introns vary from 97 base pairs to approximately 19.9 kilobases in length. The signal peptide and propeptide (von Willebrand antigen II) of von Willebrand factor are encoded by 17 exons in approximately 80 kilobases of DNA while the mature subunit of von Willebrand factor and 3' noncoding region are encoded by 35 exons in the remaining approximately 100 kilobases of the gene. A number of repetitive sequences were identified including 14 Alu repeats and a approximately 670-base pair TCTA simple repeat in intron 40 that is polymorphic. Regions of the gene that encode homologous domains have similar structures, supporting a model for their origin by gene segment duplication.     

Co-polymer tracts in eukaryotic, prokaryotic, and organellar DNA.

Large variations in DNA base composition and noticeable strand asymmetries are known to occur between different organisms and within different regions of the genomes of single organisms. Apparently such composition and sequence biases occur to fulfill structural rather than informational requirements. Here we report the wide occurrence of a more subtle biasing of DNA sequence that can have structural consequences: an increase or a suppression of the number of long tracts of two-base co-polymers. Strong biases were observed when the DNA sequences of the longest eukaryotic, prokaryotic, and organellar entries in the GenBank data base (totaling 773 kilobases) were analyzed for the number of occurrences of tracts of the two-base co-polymers (A,T)n, (G,C)n, and (A,C)n as a function of tract length. (The expression (A,T)n is used here to denote an uninterrupted tract, n nucleotides in length, of A and T bases in any proportion or order, terminated at each end by a G or C residue.) Characteristic differences are also observed in tract biases of eukaryotic vs. prokaryotic organisms.     

Magnetic ordering of DNA liquid crystals

Sonicated calf thymus DNA with an average length of approximately 100 base pairs has been found to form a cholesteric liquid crystal at a concentration of approximately 250 mg of DNA/mL of solution. Immediately after preparation, small ordered domains of a few micrometers are formed, resulting in an opaque solution. This liquid crystal can readily be oriented in the magnetic field of an NMR magnet, resulting in a clear birefringent phase. The DNA molecules align with their helix axes perpendicular to the field so that the cholesteric pitch axis was parallel with the field. A pitch length of approximately 2.5 microns for the cholesteric phase was determined both from optical measurements (optical light rotation) and from NMR measurements (solvent diffusion). The observation that DNA molecules can be magnetically oriented opens up new possibilities for studying the structure and dynamics of the aligned DNA molecules.     

Flexibility of duplex DNA on the submicrosecond timescale

Using a site-specific, Electron Paramagnetic Resonance (EPR)-active spin probe that is more rigidly locked to the DNA than any previously reported, the internal dynamics of duplex DNAs in solution were studied. EPR spectra of linear duplex DNAs containing 14-100 base pairs were acquired and simulated by the stochastic Liouville equation for anisotropic rotational diffusion using the diffusion tensor for a right circular cylinder. Internal motions have previously been assumed to be on a rapid enough time scale that they caused an averaging of the spin interactions. This assumption, however, was found to be inconsistent with the experimental data. The weakly bending rod model is modified to take into account the finite relaxation times of the internal modes and applied to analyze the EPR spectra. With this modification, the dependence of the oscillation amplitude of the probe on position along the DNA was in good agreement with the predictions of the weakly bending rod theory. From the length and position dependence of the internal flexibility of the DNA, a submicrosecond dynamic bending persistence length of around 1500 to 1700 A was found. Schellman and Harvey (Biophys. Chem. 55:95-114, 1995) have estimated that, out of the total persistence length of duplex DNA, believed to be about 500 A, approximately 1500 A is accounted for by static bends and 750 A by fluctuating bends. A measured dynamic persistence length of around 1500 A leads to the suggestion that there are additional conformations of the DNA that relax on a longer time scale than that accessible by linear CW-EPR. These measurements are the first direct determination of the dynamic flexibility of duplex DNA in 0.1 M salt.     

The subunit structure of chromatin from Physarum polycephalum.

Nucleosome DNA repeat lengths in Physarum chromatin, determined by nuclease digestion experiments, are shorter than those observed in most mammalian chromatin and longer than those reported for chromatin of certain other lower eukaryotes. After digestion with staphylococcal nuclease for short periods of time an average repeat length of 190 base pairs is measured. After more extensive digestion an average repeat length of 172 base pairs is measured. Upon prolonged digestion DNA is degraded to an average monomer subunit length of 160 base pairs, with only a small amount of DNA found in lengths of 130 base pairs or smaller. Mathematical analysis of the data suggests that the Physarum nucleosome DNA repeat comprises a protected DNA segment of about 159 base pairs with a nuclease-accessible interconnecting segment which ranges from 13 to 31 base pairs. The spacing data are compatible with measurements from electron micrographs of Physarum chromatin.     

A theory of thermal fluctuations in DNA miniplasmids.

A recent analysis of the normal modes of vibration of a ring formed by bringing together and sealing, with or without the addition of twist, the ends of rods that are straight when stress free is taken as the basis for a theory of the statistical thermodynamics of a canonical ensemble of DNA minicircles with specified linking number difference deltaLk and number N of base pairs. It is assumed that N corresponds to a circumference in the range of one or two persistence lengths. For such an ensemble, the theory yields an expression for the average writhe (Wr), which can be employed to calculate the free energy, entropy, and enthalpy of supercoiling, deltaGsc, deltaSsc, and deltaHsc. The results obtained for the dependence of deltaGsc on deltaLk and N are in accord with experimental observations of equilibrium distributions of topoisomers of plasmids with N approximately 200 bp.     

Length and sequence heterozygosity differentially affect HRAS1 minisatellite stability during meiosis in yeast

Minisatellites, one of the major classes of repetitive DNA sequences in eukaryotic genomes, are stable in somatic cells but destabilize during meiosis. We previously established a yeast model system by inserting the human Ha-ras/HRAS1 minisatellite into the HIS4 promoter and demonstrated that our system recapitulates all of the phenotypes associated with the human minisatellite. Here we demonstrate that meiotic minisatellite tract-length changes are half as frequent in diploid cells harboring heterozygous HRAS1 minisatellite tracts in which the two tracts differ by only two bases when compared to a strain with homozygous minisatellite tracts. Further, this decrease in alteration frequency is entirely dependent on DNA mismatch repair. In contrast, in a diploid strain containing heterozygous minisatellite tract alleles differing in length by three complete repeats, length alterations are observed at twice the frequency seen in a strain with homozygous tracts. Alterations consist of previously undetectable gene conversion events, plus nonparental length alteration events seen previously in strains with homozygous tracts. A strain containing tracts with both base and length heterozygosity exhibits the same level of alteration as a strain containing only length heterozygosity, indicating that base heterozygosity-dependent tract stabilization does not affect tract-length alterations occurring by gene conversion.     

The vasopressin mRNA poly(A) tract is unusually long and increases during stimulation of vasopressin gene expression in vivo.

We developed a method, termed an H-blot, by which the poly(A) tract of any specific mRNA may be detected by RNA filter hybridization after its removal from the body of the mRNA by a RNase H-catalyzed endonucleolytic cleavage in the 3' untranslated region. Using this method, we studied the modulation of the length of the poly(A) tract of rat vasopressin mRNA in vivo during changes in the levels of this mRNA resulting from a physiologic stimulus, osmotic stress. The poly(A) tract of hypothalamic vasopressin mRNA in hydrated rats was, quite remarkably, approximately 250 nucleotides in length, in contrast to that of somatostatin mRNA, which was approximately 30 nucleotides long. Vasopressin mRNA poly(A) tail length increased progressively from approximately 250 to approximately 400 nucleotides with the application of the hyperosmotic stimulus and declined to base line after its removal; somatostatin mRNA poly(A) tail length did not change during osmotic stress. The poly(A) tract length of total hypothalamic mRNA was between 35 and 140 nucleotides and also did not change with osmotic stress. Modulation of poly(A) tract length of specific mRNAs during stimulation of gene expression may provide an additional level of genetic regulation.     

Evidence for long poly(dA).poly(dT) tracts in D. discoideum DNA at high frequencies and their preferential avoidance of nucleosomal DNA core regions

The eukaryote, Dictyostelium discoideum, has one of the most (A+T) rich genomes studied to date. Isolated nuclear D. discoideum DNA (AX3 strain) was used to qualitatively determine the frequency and length distribution of long (dA).(dT) homopolymer tracts in this genome, in comparison to the less (A+T) rich calf thymus and Schistosoma mansoni DNAs that had few observable long tracts. These experimental data accurately reflect the significantly elevated frequencies of long tracts found computationally within the D. discoideum intron and flanking sequences, but not exons. PCR amplification of long (dA).(dT) homopolymer tract containing sequences was carried out. Then experimental biotinylated (dT)18 probe hybridization to the PCR amplified DNA showed that the long (dA).(dT) homopolymer tracts were enriched in D. discoideum sequences only hundreds of base pair in length, under conditions where no equivalent hybridization was observed to S. mansoni DNA or calf DNA sequences. Similar probe hybridization to DNA isolated following micrococcal nuclease digestion of D. discoideum chromatin demonstrated that long (dA).(dT) homopolymer tracts were more highly enriched in nucleosomal DNA lengths that included the internucleosomal linker as compared to shorter linker free mononucleosomal lengths. This observation is in agreement with the frequency of tract spacing results calculated from GenBank sequence data. These frequency data indicate that adjacent long tracts plus the intervening spacer DNA are found at peak lengths (average 42 bp), exactly characteristic of the internucleosomal spacer region of D. discoideum chromatin and are in sufficient number to be found in nearly half of all nucleosomes. Compared to shuffled tract sequence controls, these lengths of adjacent long tracts plus the intervening spacer DNA were found to be significantly enriched. Lesser enrichments are observed at lengths corresponding to adjacent tracts being separated by nucleosomal core length DNA sequences (145-185 bp). These data strongly suggest that adjacent long tracts occur spaced at selected lengths so as to avoid the central core regions of nucleosomes and instead are found localized within internucleosomal DNA linker and core edge regions in D. discoideum chromatin.     

Isolation and properties of the initiation of DNA replication in the sea urchin embryo

Using the instability of replication loops as a method for the isolation of double-stranded nascent DNA from embryos of the sea urchin Strongylocentrotus intermedius, extruded DNA enriched for replication origins was obtained. The average length of the fragments of the DNA of this fraction was estimated to be about 800 base pairs. The origin-enriched nascent DNA strands were assayed for the presence of inverted repeats. The results show that the origin-enriched DNA is also enriched in inverted repeats. The bulk of palindromes in the total nuclear DNA was estimated to be 200 base pairs in length and from the origin-enriched DNA-150 base pairs.     

Meiotic Gene Conversion Tract Length Distribution within the Rosy Locus of Drosophila Melanogaster

Employing extensive co-conversion data for selected and unselected sites of known molecular location in the rosy locus of Drosophila melanogaster, we determine the parameters of meiotic gene conversion tract length distribution. The tract length distribution for gene conversion events can be approximated by the equation P(L >/= n) = (n) where P is the probability that tract length (L) is greater than or equal to a specified number of nucleotides (n). From the co-conversion data, a maximum likelihood estimate with standard error for is 0.99717 +/- 0.00026, corresponding to a mean conversion tract length of 352 base pairs. (Thus, gene conversion tract lengths are sufficiently small to allow for extensive shuffling of DNA sequence polymorphisms within a gene.) For selected site conversions there is a bias towards recovery of longer tracts. The distribution of conversion tract lengths associated with selected sites can be approximated by the equation P(L >/= n| selected = (n)(1 - n + n/), where P is now the probability that a selected site tract length (L) is greater than or equal to a specified number of nucleotides (n). For the optimal value of determined from the co-conversion analysis, the mean conversion tract length for selected sites is 706 base pairs. We discuss, in the light of this and other studies, the relationship between meiotic gene conversion and P element excision induced gap repair and determine that they are distinct processes defined by different parameters and, possibly, mechanisms.