Search

[Abstract] Hypothetical proteins (HP) are those that are not characterized in the laboratory and so remain “orphaned” in genomic databases. In recent times there has been a lot of progress in characterizing HPs in the laboratory. Various methods, such as sequence capture and Next Generation Sequencing (NGS), have been used to rapidly identify HP functions and their encoded genes. Applications and methods, such as the isolation of single genes, are greatly facilitated by pull-down assays to characterize proteins. Furthermore, there are methods to extract proteins from either the whole cell or a subcellular fraction. But the weakness is that some assays are fairly expensive and laborious, and characterizing HP function is always imperfect. In the recent past, statistical interpretations of the in silico selection strategies have improved the identification of the most promising candidates, including those from various annotation methods, such as protein interaction networks (PIN). Given the improvements in technology that have permitted a substantial increase in computational annotation, we ask if the prediction of HP function in silico (validation of models through algorithms and data subsets) could likewise be improved. In this work, we apply a bioinformatics analogy to each step of a wet lab experiment performed to predict aspects confirming protein function. Although it may be a less bona fide approach, assigning a putative function from conservation observed in homologous protein sequences might be worthwhile to consider prior to a wet lab experiment.

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.

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.

There has recently been an increase in the number of open access journals showcasing the results of research, free of charge, in an affordable and easy to access online publication. In fact, it's a paradigm shift in publishing, and it has gained so much momentum and has become so favored institutionally that perhaps we can say that the one who shares, wins. We do wish to acknowledge the good that open access has achieved through journal readership, but we also want to mention here some of the problems and challenges brought about by these changes. These are issues we think authors should be aware of before submitting to open access journals.

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.

High resolution melting curves of total nuclearDictyostelium discoideumDNA (A×3 strain) are compared to theoretical melting calculated from GENBANK sequences(1.74 % of total) by the statistical thermodynamics program MELTSIM, parameterized for long DNA sequences(1,2). The lower and upper limits of simulated melting agree quantitatively with the experimental melting of total DNA. Calculated melting of coding, intron and flanking regions indicate that intron and flanking DNAs are extremely (A+T)-rich and account for the earliest melting DNA. There is no temperature overlap of these regions with coding DNA. A theoretical denaturation map of DNA containing the ribosomal DNA genes showed excellent agreement with subtransition positions of these genes in experimental curves. Agreement between these calculated and experimental melting data demonstrates our ability to accurately simulate DNA melting in complex eukaryotic genomes. This has important consequences for the understanding of sequence dependent energetic properties of nucleic acids and their potential use as biomaterials.