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.     

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.     

Local enrichment with homopolymeric (dA/dT) DNA in genomes of some lower dipterans and Drosophila melanogaster

An investigation into the chromosomal localization of homopolymeric dA/dT was carried out with species of the genera Rhynchosciara, Chironomus, Drosophila and several other taxa. In situ hybridisation probing mitotic and polytene chromosomes with RNA homopolymers was performed, followed by immunological detection of the DNA/RNA hybrid. Use of this method allowed us to assess specific regions of some dipteran genomes, where the signal was generally, but not always, located in heterochromatic regions. Human and Drosophila chromosome regions known to contain dA/dT runs of up to 153 bp were devoid of consistent labelling. The stability of the rA/dT hybrid formed in situ was in agreement with the T(m) for long rA/dT hybrid complexes, suggesting that the method used in this work is able to identify unusually long homopolymeric dA/dT tracts.     

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.     

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.     

kDNA minicircles of the major sequence class of C. fasciculata contain a single region of bent helix widely separated from the two origins of replication

The major sequence class of Crithidia fasciculata minicircles is shown to have a single region of bent helical DNA widely separated from the two replication origins located 180 degrees apart on the minicircle map. The position of the bend in the DNA has been mapped both by gel electrophoretic methods and by direct electron microscopic observation of the DNA. This sequence directed bending is apparently the result of homopolymeric dA X dT tracts 4-6 base pairs long repeated in phase with the helix screw. The region of the bend contains nineteen such homopolymeric tracts in a region of about 200 base pairs with sixteen of the tracts oriented in the same direction.     

The preferential binding of histone H1 to DNA scaffold-associated regions is determined by its C-terminal domain

Histone H1 preferentially binds and aggregates scaffold-associated regions (SARs) via the numerous homopolymeric oligo(dA).oligo(dT) tracts present within these sequences. Here we show that the mammalian somatic subtypes H1a,b,c,d,e and H1° and the male germline-specific subtype H1t, all preferentially bind to the Drosophila histone SAR. Experiments with the isolated domains show that whilst the C-terminal domain maintains strong and preferential binding, the N-terminal and globular domains show weak binding and poor specificity for the SAR. The preferential binding of SAR by the H1 molecule thus appears to be determined by its highly basic C-terminal domain. Salmine, a typical fish protamine, which could have its evolutionary origin in histone H1, also shows preferential binding to the SAR. The interaction of distamycin, a minor groove binder with high affinity for homopolymeric oligo(dA).oligo(dT) tracts, abolishes preferential binding of the C-terminal domain of histone H1 and protamine to the SAR, suggesting the involvement of the DNA minor groove in the interaction.     

Simple sequences are ubiquitous repetitive components of eukaryotic genomes

Simple sequences are stretches of DNA which consist of only one, or a few tandemly repeated nucleotides, for example poly (dA) X poly (dT) or poly (dG-dT) X poly (dC-dA). These two types of simple sequence have been shown to be repetitive and interspersed in many eukaryotic genomes. Several other types have been found by sequencing eukaryotic DNA. In this report we have undertaken a systematical survey for simple sequences. We hybridized synthetical simple sequence DNA to genome blots of phylogenetically different organisms. We found that many, probably even all possible types of simple sequence are repetitive components of eukaryotic genomes. We propose therefore that they arise by common mechanisms namely slippage replication and unequal crossover and that they might have no general function with regards to gene expression. This latter inference is supported by the fact that we have detected simple sequences only in the metabolically inactive micronucleus of the protozoan Stylonychia, but not in the metabolically active macronucleus which is derived from the micronucleus by chromosome diminution.     

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.     

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.     

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.     

Simple repetitive sequences in the genomes of archaebacteria

Stretches of simple sequences poly(dG-dT).poly(dC-dA), poly(dG-dA).poly(dC-dT), poly(dG).poly(dC) and poly(dA).poly(dT), the occurrence of which is a characteristic feature of eukaryotic genomes, are found in the genomes of archaebacteria Halobacterium halobium and Sulfolobus acidocaldarius. In S. acidocaldarius these sequences constitute a considerable portion of the genome; they belong to a class of repetitive sequences dispersed throughout the genome, being transcribed and found in RNAs of different lengths.     

Anomalous structure and properties of poly (dA).poly(dT). Computer simulation of the polynucleotide structure with the spine of hydration in the minor groove.

The results of the search for low-energy conformations of poly(dA).poly(dT) and of the poly(dA).poly(dT) "complex" with the spine of hydration similar to that found by Dickerson and co-workers (Kopka, M.L., Fratini, A.V., Drew, H.R. and Dickerson, R.E. (1983) J. Mol. Biol. 163, 129-146) in the minor groove of the CGCGAATTCGCG crystals are described. It is shown that the existence of such a spine in the minor groove of poly(dA).poly(dT) is energetically favourable. Moreover, the spine of hydration makes the polynucleotide conformation similar to the poly(dA).poly(dT) structure in fibers and to the conformation of the central part of CGCGAATTCGCG in crystals; it also acquires features characteristic of the structure of poly(dA).poly(dT) and DNA oligo(dA)-tracts in solution. It is shown that the existence of the TpA step in conformations characteristic of the poly(dA).poly(dT) complex with the spine of hydration is energetically unfavourable (in contrast to the ApT step) and therefore this step should result in destabilization of the spine of hydration in the DNA minor groove. Thus, it appears that the spine of hydration as described by Dickerson and co-workers is unlikely to exist in the poly d(A-T).poly d(A-T) structure. The data obtained permit us to interpret a large body of experimental facts concerning the unusual structure and properties of poly(dA).poly(dT) and oligo(dA)-tracts in DNA both in fibers and in solution. The results provide evidence of the existence of the minor groove spine of hydration both in fibers and in solution on A/T tracts of DNA which do not contain the TpA step. The spine plays an active role in the formation of the anomalous conformation of these tracts.     

Distamycin inhibition of topoisomerase I-DNA interaction: a mechanistic analysis.

Inhibition of eukaryotic DNA topoisomerase I by the minor groove binding ligand, distamycin A, was investigated. Low concentrations of the ligand selectively prevented catalytic action at a high affinity topoisomerase I binding sequence. A restriction enzyme protection assay indicated that the catalytic cycle was blocked at the binding step. Distamycin binding sites on DNA were localized by hydroxyl radical footprinting. A strongly preferred site mapped to a homopolymeric (dA).(dT)-tract partially included in the essential topoisomerase I binding region. Mutational elimination of the stable helix curvature associated with this ligand binding site demonstrated that (i) the intrinsic bend was unessential for efficient binding of topoisomerase I, and (ii) distamycin inhibition did not occur by deformation of a stable band. Alternative modes of inhibition are discussed.