Search

The slime mold, Dictyostelium discoideum, possesses an (A+T) rich eukaryotic genome that is being sequenced in the Human Genome Project. High resolution melting curves of isolated total and fractionated nuclear D. discoideum DNA(AX3 strain) were determined experimentally and are compared to melting curves calculated from GENBANK sequences (1.59% of genome) by the statistical thermodynamics program MELTSIM (1), parameterized for long DNA sequences (2,3). The lower and upper temperature limits of calculated melting agree well with the observed melting of total DNA. The experimental curve is unusual in that it contains a number of sharp peaks. MELTSIM allowed us to calculate positional denaturation maps of D. discoideum GENBANK sequence documents containing the 26S, 5.8S and 17S rDNA gene sequences, a major satellite DNA and repetitive sequence family present in 100-200 copies/nucleus. These denaturation maps contain subtransitions that correspond with a number of the experimentally observed peaks, some of which we show to correspond with rDNA gene enriched CsCl gradient fractions of D. discoideum DNA. MELTSIM calculated curves of coding, intron and flanking sequences indicate that both intron and flanking sequences are extremely (A+T) rich and account for most of the low temperature melting. There is no temperature overlap between thermal stabilities of these sequence domains and those of coding DNA. The latter must satisfy triplet codon constraints of higher (G+C) content. These large stability property differences enable a denaturation mapping feature of MELTSIM to clearly distinguish exon positions from those of introns and flanking DNA in long D. discoideum gene containing sequences.

1 International Institute of Tropical Agriculture, Nairobi, Kenya, 2 Faculty of Life Sciences, The University of Manchester, Manchester, United Kingdom, 3 Department of Computer and Information Sciences, Covenant University, Ota, Nigeria, 4 Institute of Bioinformatics, Johannes Kepler University, Linz, Austria, 5 Moroccan Society for Bioinformatics Institute, Morocco, 6 South African National Bioinformatics Institute, University of the Western Cape, Bellville, South Africa, 7 University of Cape Town, Cape Town, South Africa, 8 University of Notre Dame, South Bend, Indiana, United States of America, 9 Biotechnology Unit, University of Buea, Buea, South West Region, Cameroon, 10 International Livestock Research Institute, Nairobi, Kenya, 11 Biosciences Eastern and Central Africa, Nairobi, Kenya, 12 International Center of Insect Physiology and Ecology, Nairobi, Kenya, 13 Bioinformatics Organization, Hudson, Massachusetts, United States of America, 14 Bioinformatics Team, Center for Development of Advanced Computing, Pune University Campus, Pune, India, 15 Harvard School of Public Health, Boston, Massachusetts, United States of America

Long poly dA.dT tracts are enriched within the intergenic sequences of eukaryotic genomes where they appear to act as intrinsic regulators of nucleosome positioning. Previous studies of the incompl...

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. Results We carried out a survey of the extent of homopolymer tract over-representation (enrichment) and over-proportional length distribution (above expected length) primarily in the single gene documents, but including some whole chromosomes of 27 eukaryotics across the (G+C)% composition range from 20 – 60%. A total of 5.2 × 107 bases from 15,560 cleaned (redundancy removed) sequence documents were analyzed. Calculated frequencies of non-overlapping long homopolymer tracts were found over-represented in non-coding sequences of eukaryotes. Long poly(dA).poly(dT) tracts demonstrated an exponential increase with tract length compared to predicted frequencies. A novel negative slope was observed for all eukaryotes between their (G+C)% composition and the threshold length N where poly(dA).poly(dT) tracts exhibited over-representation and a corresponding positive slope was observed for poly(dG).poly(dC) tracts. Tract size thresholds where over-representation of tracts in different eukaryotes began to occur was between 4 – 11 bp depending upon the organism (G+C)% composition. The higher the GC%, the lower the threshold N value was for poly(dA).poly(dT) tracts, meaning that the over-representation happens at relatively lower tract length in more GC-rich surrounding sequence. We also observed a novel relationship between the highest over-representations, as well as lengths of homopolymer tracts in excess of their random occurrence expected maximum lengths. Conclusions We discuss how our novel tract over-representation observations can be accounted for by a few models. A likely model for poly(dA).poly(dT) tract over-representation involves the known insertion into genomes of DNA synthesized from retroviral mRNAs containing 3' polyA tails. A proposed model that can account for a number of our observed results, concerns the origin of the isochore nature of eukaryotic genomes via a non-equilibrium GC% dependent mutation rate mechanism. Our data also suggest that tract lengthening via slip strand replication is not governed by a simple thermodynamic loop energy model.

The yeast S. cereviseae represents the first eukaryotic organism whose genome has been entirely sequenced as a result of the Human Genome Project(1). In this report we demonstrate the good agreement between an experimental high resolution melting curve of total nuclear S. cereviseae DNA and the theoretical melting calculated for the complete yeast DNA genome (12,067,277 bp: Saccharomyces Genome Database) by the statistical thermodynamics program MELTSIM, parameterized for long DNA sequences(2,3). The experimental and theoretical melting curves are both fairly symmetrical and possess nearly identical Tmvalues. Calculated melting of coding and flanking DNA regions indicates that flanking DNAs are more (A+T)-rich than coding sequences and account for the earliest melting DNA. Calculated melting curves of the 16 individual yeast chromosomes are very similar and with few exceptions exhibit symmetric melting curves. MELTSIM was also used to calculate a theoretical denaturation map of Chromosome III DNA. The agreement between MELTSIM calculated and experimental melting data demonstrates our ability to accurately simulate long DNA sequence melting in complex eukaryotic genomes, whose sequences are becoming increasingly available for study in public databases. This has important consequences for the understanding of sequence dependent energetic properties of DNA in their biological sequence context and also for their potential use in biomaterials applications.