Alignment of (dA)·(dT) Homopolymer Tracts in Gene Flanking Sequences Suggests Nucleosomal Periodicity in D. discoideum DNA

Abstract

It has been shown that the frequency versus size distribution of A and T overlapping and non- overlapping homopolymer tracts of N>5 in D. discoideum gene flanking and intron regions are significantly greater than in coding regions(1). In the present report we demonstrate, that a spatial periodicity exists in long A and T tracts (N> 10) in long flanking sequences by scored alignments of those tracts (N> 10) with the nucleosomal repeat. A tract spacing was found at 185–190 bp that corresponds to a maximum alignment score. This is exactly the average spacing of D. discoideum nucleosomes determined experimentally. A majority of A and T tracts in flanking sequences are often spaced by short DNA stretches and the total length of adjacent A and T tracts plus the interrupting short DNA stretch corresponds closely to the average experimentally measured nucleosomal linker DNA size in D. discoideum-42 bp. These data suggest a model which has A and T runs of N> 10 bp in flanking DNA of D. discoideum organized in a regular phase with nonhomopolymer sequences along the DNA. This model has functional implications for A and T tracts, suggesting that they are found in nucleosomal linker DNA regions of chromatin during some necessary portion(s) of the life of the cell.     

The binding of poly (rA) and poly (rU) to denatured DNA. II. Studies with natural DNAs.

We have studied the interaction of poly(rA) and poly(rU) with natural DNAs containing (dA.dT)n sequences. The results indicate that hybridization of poly(rA) to denatured DNA can be used to estimate the size and frequency of large (dA.dT)n tracts, whereas hybridization with poly(rU) does not give reliable information on these points. In 6.6 M CsCl, poly(rU) can form stable complexes with denatured DNA containing short (dA)n tracts (n less than or equal to 6), whereas binding of poly(rA) to denatured DNA under these conditions requires much larger (dT)n tracts (estimated n greater than 13). Moreover, binding of poly(rA) requires pre-hybridization in low salt, because free poly(rA) precipitates in 6.6 M CsCl.     

Homopolymer tract length dependent enrichments in functional regions of 27 eukaryotes and their novel dependence on the organism DNA (G+C)% composition

Background  

DNA homopolymer tracts, poly(dA).poly(dT) and poly(dG).poly(dC), are the simplest of simple sequence repeats. Homopolymer tracts have been systematically examined in the coding, intron and flanking regions of a limited number of eukaryotes. As the number of DNA sequences publicly available increases, the representation (over and under) of homopolymer tracts of different lengths in these regions of different genomes can be compared.     

Long poly(A) tracts in the human genome are associated with the Alu family of repeated elements

Long poly(dA).poly(dT) tracts (poly(A) tracts), regions of DNA containing at least 20 contiguous dA residues on one strand and dT residues on the complementary strand, are found in about 2 X 10(4) copies interspersed throughout the human genome. Using poly(dA).poly(dA) as a hybridization probe, we identified recombinant lambda phage that contained inserts of human DNA with poly(A) tracts. Three such tracts have been characterized by restriction mapping and sequence analysis. One major poly(A) tract is present within each insert and is composed of from 28 to 35 A residues. In each case, the poly(A) tract directly abuts the 3' end of the human Alu element, indicating that the major class of poly(A) tracts in the human genome is associated with this family of repeats. The poly(A) tracts are also adjacent to A-rich sequences and, in one case, to a polypurine tract, having the structure GA3-GA3-GA4-GA6-GA5-GA4. We suggest that repetitive cycles of unequal crossing over may give rise to both the long poly(A) and polypurine tracts observed in this study.     

Nucleosomes will not form on double-stranded RNa or over poly(dA).poly(dT) tracts in recombinant DNA.

We have been unable to "force" double-stranded RNA to fold into nucleosome-like structures using several different histone-RNA "reconstitution" procedures. Even if the histones are first stabilized in octameric form by dimethylsuberimidate cross-linking they are still unable to form specific complexes with the RNA. Moreover double-stranded RNA is unable to induce histones to assemble into octamers although we confirm that the non-nucleic acid homopolymer, polyglutamic acid, has this ability. We have also determined, using pyrimidine tract analysis, that nucleosomes will not form over a sufficiently long segment of poly(dA).poly(dT) in a recombinant DNA molecule. Thus nucleosomes cannot fold DNA containing an 80 base pair poly(dA).poly(dT) segment but a 20 base pair segment can be accommodated in nucleosomes fairly well. Segments of intermediate length can be accommodated but are clearly selected against. Poly(dA).poly(dT) differs only slightly from natural DNA in helix structure. Therefore either this homopolymer resists folding, or nucleosomes are very exacting in the nucleic acid steroid parameters they will tolerate. Such constraints may be relevant to nucleosome positioning in chromatin.     

DNA curvature does not require bifurcated hydrogen bonds or pyrimidine methyl groups.

Short tracts of the homopolymer dA.dT confer intrinsic curvature on the axis of the DNA double helix. This phenomenon is assumed to be a consequence of such tracts adopting a stable B'-DNA conformation that is distinct from B-form structure normally assumed by other DNA sequences. The more stable B' structure of dA.dT tracts has been attributed to several possible stabilizing factors: (1) optimal base stacking interactions consequent upon the high propeller twist, (2) bifurcated hydrogen bonds between adjacent dA.dT base-pairs, (3) stacking interactions involving the dT methyl groups, and finally (4) a putative spine of ordered water molecules in the minor groove. DNA oligodeoxynucleotides have been synthesized that enable these hypotheses to be tested; of particular interest is the combination of effects due to bifurcation (2) and methylation of the pyrimidines nucleotides (3). The data indicate that neither bifurcated hydrogen bonds nor pyrimidine methyl groups nor both are essential for DNA curvature. The data further suggest that the influence of the minor groove spine of hydration on the B'-formation is small. The experiments favor the hypothesis that base stacking interactions are the dominant force in stabilizing the B'-form structure.     

Structural repeat units of Chinese hamster ovary chromatin. Evidence for variations in repeat unit DNA size in higher eukaryotes.

DNA lengths in the structural repeat units of Chinese hamster ovary (CHO) and chicken erythrocyte chromatin were compared by analyzing the sizes of DNA fragments produced after treatment of nuclei with staphylococcal nuclease. The repeat length of CHO chromatin (173 +- 4 BP) is about 20 base pairs (BP) smaller than that of chicken erythrocyte chromatin (194 +- 8 BP). Repeat lengths of rat liver and calf thymus chromatin were found to be about 10 BP shorter than that of chicken erythrocyte chromatin. Thus significant variations occur in repeat units of chromatin of higher eukaryotes. These variations occur in the lengths of "spacer" (or "internucleosomal") DNA segments, not in "core particle" (or "nucleosomal") DNA lengths. The concept of spacer regions and the possible influence of H1 histones is discussed.     

Differential sequence dynamics of homopolymeric and alternating AT tracts in a small plasmid DNA

The location of OsO4 bispyridine hyper- and hyporeactivity in a small deletion derivative of plasmid ColE1 (PTC12, 1727 bp) has been determined for approximately 70% of the molecule. Thymine bases in homopolymeric (dA)n.(dT)n tracts (n greater than or equal to 4) were always found to be resistant toward OsO4 modification. DNA supercoiling did not destabilize these tracts. The extent of OsO4 bispyridine reactivity of homopolymeric (dA)n.(dT)n tracts, where n = 3, was found to be dependent on the rate of base unpairing of the sequence immediately 5' and 3' to the tract. Repressed OsO4 reactivity of thymine bases in (dA)3.(dT)3 tracts was observed if immediately both 5' and 3' to the tract were stable DNA sequences composed of GC base pairs and/or a homopolymeric (dA)n.(dT)n tract (n greater than or equal to 4). Homopolymeric tracts of n = 3 not having adjacent sequences with repressed unpairing rates did not show reduced levels of OsO4 bispyridine reactivity. Alternating d(TA)n tracts (n greater than or equal to 2) were found to exhibit hyperreactivity with OsO4. The extent of this hyperreactivity was dependent on the length of the tract and superhelical torsional stress. The distribution and frequency of homopolymeric (dA)n.(dT)n (n greater than or equal to 4) tracts in Escherichia coli promoter sequences were examined, and the possible implications of these tracts on promoter function are discussed.     

Inhibition of DNA replication by berenil in plasmids containing Poly(dA)poly(dT) sequences

The effect of berenil on plasmid DNA replication was studied on pBR322-derived plasmids containing poly(dA)poly(dT) sequences. In comparison to the parental plasmid pBR322, plasmid pKH47 harboring 100 bp of poly(dA)poly(dT) at the PvuII site showed a decrease in plasmid yield in the presence of berenil. This effect was also observed in pVL26, a related plasmid in which the location of the poly(dA)poly(dT) region had been shifted to the EcoRV site in pBR322. [(3)H]Thymidine incorporation experiments indicated that DNA synthesis may be affected in these plasmids in the presence of the drug. Bromodeoxyuridine incorporation experiments coupled to Cs(2)SO(4) equilibrium density gradient centrifugation indicated that the lower plasmid yield was due to an inhibition of DNA replication by berenil. We have also found that berenil induces DNA degradation in plasmids containing the homopolymer. Our studies strongly suggest that the effect of berenil on plasmid replication and DNA stability results from its binding to the poly(dA)poly(dT) region present in these plasmids. Moreover, we have found a correlation between the position of the poly(dA)poly(dT) region and this inhibitory effect. Thus, plasmid pKH47, containing the poly(dA)poly(dT) region most proximal to the origin of pBR322 replication, was most severely affected.     

Preferential, cooperative binding of DNA topoisomerase II to scaffold-associated regions.

DNA elements termed scaffold-associated regions (SARs) are AT-rich stretches of several hundred base pairs which are known to bind specifically to nuclear or metaphase scaffolds and are proposed to specify the base of chromatin loops. SARs contain sequences homologous to the DNA topoisomerase II cleavage consensus and this enzyme is known to be the major structural component of the mitotic chromosome scaffold. We find that purified topoisomerase II preferentially binds and aggregates SAR-containing DNA. This interaction is highly cooperative and, with increasing concentrations of topoisomerase II, the protein titrates quantitatively first SAR-containing DNA and then non-SAR DNA. About one topoisomerase II dimer is bound per 200 bp of DNA. SARs exhibit a Circe effect; they promote in cis topoisomerase II-mediated double-strand cleavage in SAR-containing DNA fragments. The AT-rich SARs contain several oligo(dA).oligo(dT) tracts which determine their protein-binding specificity. Distamycin, which is known to interact highly selectively with runs of A.T base pairs, abolishes the specific interaction of SARs with topoisomerase II, and the homopolymer oligo(dA).oligo(dT) is, above a critical length of 240 bp, a highly specific artificial SAR. These results support the notion of an involvement of SARs and topoisomerase II in chromosome structure.     

Homopolymer tract organization in the human malarial parasite Plasmodium falciparum and related Apicomplexan parasites

Background

Homopolymeric tracts, particularly poly dA.dT, are enriched within the intergenic sequences of eukaryotic genomes where they appear to act as intrinsic regulators of nucleosome positioning. A previous study of the incomplete genome of the human malarial parasite Plasmodium falciparum reports a higher than expected enrichment of poly dA.dT tracts, far above that anticipated even in this highly AT rich genome. Here we report an analysis of the relative frequency, length and spatial arrangement of homopolymer tracts for the complete P. falciparum genome, extending this analysis to twelve additional genomes of Apicomplexan parasites important to human and animal health. In addition, using nucleosome-positioning data available for P. falciparum, we explore the correlation of poly dA.dT tracts with nucleosome-positioning data over key expression landmarks within intergenic regions.

Results

We describe three apparent lineage-specific patterns of homopolymeric tract organization within the intergenic regions of these Apicomplexan parasites. Moreover, a striking pattern of enrichment of overly long poly dA.dT tracts in the intergenic regions of Plasmodium spp. uniquely extends into protein coding sequences. There is a conserved spatial arrangement of poly dA.dT immediately flanking open reading frames and over predicted core promoter sites. These key landmarks are all relatively depleted in nucleosomes in P. falciparum, as would be expected for poly dA.dT acting as nucleosome exclusion sequences.

Conclusions

Previous comparative studies of homopolymer tract organization emphasize evolutionary diversity; this is the first report of such an analysis within a single phylum. Our data provide insights into the evolution of homopolymeric tracts and the selective pressures at play in their maintenance and expansion.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-848) contains supplementary material, which is available to authorized users.     

Sequence context of oligomer tracts in eukaryotic DNA: biological and conformational implications

Recent studies of homooligomer tracts suggest different characteristics from random sequence DNA (dA).(dT) and (dG).(dC) tracts are frequent in upstream regions and in some cases have been shown to be essential for regulation. Here we examine homooligomer occurrences in non-coding and coding eukaryotic sequences, focusing on the context in which the homooligomers occur. This analysis of sequences in the junction areas yields distinct and consistent characteristics. In particular, the nucleotide interrupting a run is most frequently complementary to the run. The base next to it is most frequently identical to the one constituting the run. For A or T runs the least frequent nearest and next to nearest neighbors are G or C. For G or C tracts the least frequent are A or T. Complementary oligomers behave similarly. These and additional trends are strongest for run lengths greater than or equal to 3. The computations are carried out on the whole eukaryotic database of greater than 4 x 10(6) nucleotides, separately for coding and non-coding regions. These same trends are evident for both groups, but are somewhat stronger for the non-coding regions. The context in which the homooligomers occur may yield some clues to DNA conformation and its biological implications.