Poly(dC) segments and cloning artifacts in databases

Nuclear ribosomal pseudogenes resolve a corroborated monophyly of the eucalypt genus Corymbia despite misleading hypotheses at functional ITS paralogs

Divergent paralogs can create both obstacles and opportunities for phylogenetic reconstruction. Phylogenetic relationships among eucalypt genera have been incongruent among datasets in previous studies, where morphological characters supported monophyly of the genus Corymbia, while intergenic spacers of the nuclear ribosomes (ITS) and chloroplast loci (trnL, trnH, psbA) showed Corymbia as either equivocal or paraphyletic. Ribosomal DNA occurs in multiple copies in a genome. We cloned and sequenced the nrITS to investigate if gene duplication was the cause of incongruence among trees in the eucalypts. Three ITS riboforms, two of them widespread, were recovered within some genomes. One of the ITS riboforms recovered a robust phylogeny showing Corymbia as a monophyletic genus, corroborating the evidence from morphology, fossil data, a recent ITS/ETS dataset and microsatellites (SSRs). Compelling evidence suggested that this divergent riboform is a pseudogene, i.e., non-functional paralog: comparatively lower GC content suggesting lower structural stability, deamination-like mutations at potential methylation sites, lack of conserved helices and hairpins and conspicuously lower thermodynamic stability in secondary structures. Phylogenies from the apparently functional riboform retained Corymbia as paraphyletic. We show here that pseudogenes can recover a well-corroborated phylogeny whereas their functional paralogs show misleading hypotheses. We explain that phylogenetic signals may be obscured when functional constraints in ITS necessitate compensatory mutations in the secondary structure helices involved in RNA transcription, whereas pseudogenes mutate under neutrality.

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.

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 × 10⁷ 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.

Why nuclear ribosomal DNA spacers (ITS) tell different stories in Quercus

The molecular systematics of Quercus (Fagaceae) was recently assessed by two teams using independently generated ITS sequences. Although the results disagreed in several remarkable features, the phylogenetic trees for either hypothesis were highly supported by bootstrap resampling. We have reanalyzed the ITS sequences used by both teams (eight taxa) to reveal the underlying patterns of this divergence. Within species, conspicuous length and G + C% divergence were evident in most sequence comparisons. In addition, a high rate of substitutions and deletions involving highly conserved motifs in both ITS spacers were present in a set of sequences. This was coupled with a less thermodynamic stability in the RNA structure, lacking conserved hairpins that are putatively involved in the processing of the RNA transcripts. Compelling evidence suggests that the divergent ITS alleles reported by one team are pseudogenes, i.e., nonfunctional paralogous loci. The hypothesis that the contrasting phylogenetic histories drawn from Quercus using ITS data are not strictly related to technical differences between laboratories, but that they have rather been generated from the analysis of paralogous sequences, best reconciles the available data. The risk of incorporating ITS paralogues in plant evolutionary studies which can distortion the phylogenetic signal should caution molecular systematists. Without a detailed inspection of some basic features of the sequence, including the integrity of the conserved motifs and the thermodynamic stability of the secondary structures of the RNA transcripts, errors in evolutionary inferences could be easily overlooked.

DNA sequencing and comparative sequence analysis reveal that the Escherichia coli genomic DNA may replace the target DNA during molecular cloning: evidence for the erroneous assembly of E. coli DNA into database sequences

DNA sequencing and similarity search of databases provide experimental evidence that portions of the host Escherichia coli genome may get ligated into the cloning vector, resulting in clones containing nontargeted inserts. Several lines of evidence suggest that this non-targeted ligation, as observed by us while subcloning troponin I cDNA, is presumably due to a recombination-mediated mechanism by which host DNA replaces the target DNA in the cloning vector. The E. coli genome mapping to 64-65 min and 92.8-00.1 min, the latter containing insertion sequences, appears to be the hotspot regions involved in this process. We examined the possibility that some sequences reported in the databases may also contain genomic sequences of E. coli. A search of current databases revealed that a rat hepatic glutathione transporter cDNA contains a 2.2-kb-long portion of the E. coli genome that has been wrongly assembled into its 5' untranslated and coding regions. In addition, about 30 sequences in databases, including a Yersinia pestis toxin gene, showed relatively high sequence identity with those portions of the E. coli genome that were present in the nonauthentic clones.