Detecting horizontally transferred and essential genes based on dinucleotide relative abundance.

Various methods have been developed to detect horizontal gene transfer in bacteria, based on anomalous nucleotide composition, assuming that compositional features undergo amelioration in the host genome. Evolutionary theory predicts the inevitability of false positives when essential sequences are strongly conserved. Foreign genes could become more detectable on the basis of their higher order compositions if such features ameliorate more rapidly and uniformly than lower order features. This possibility is tested by comparing the heterogeneities of bacterial genomes with respect to strand-independent first- and second-order features, (i) G + C content and (ii) dinucleotide relative abundance, in 1 kb segments. Although statistical analysis confirms that (ii) is less inhomogeneous than (i) in all 12 species examined, extreme anomalies with respect to (ii) in the Escherichia coli K12 genome are typically co-located with essential genes.

Methods for comparing sources of strand compositional asymmetry in microbial chromosomes.

Significant compositional biases in bacterial chromosomes have been explained by replication- and transcription-coupled repair mechanisms, the latter causing GC skew to indicate the direction of replication when gene polarity is correspondingly entrained. Correlations between indicators of replication direction, skew, and transcription polarity are computed for the complete nucleotide sequences of 20 microbial chromosomes and interpreted through statistical tests. A second quantitative method, previously applied to the first complete draft of the Escherichia coli K12 genome, characterizes the sequences by average skew and net skew due to replication. These methods generally agree in finding the coexistence of replication- and translation-coupled effects and in identifying atypical sequences in which one influence is clearly dominant. The replication-dominated class is exemplified by two chlamydial sequences and the transcription-dominated class by three archaea. The preference for leading-strand transcription in two mycoplasmas is stronger than the skew implies. These concordant methods provide an objective framework for comparing sources of strand compositional asymmetry and interpreting skew diagrams.

Pervasive properties of the genomic signature.

BACKGROUND: The dinucleotide relative abundance profile can be regarded as a genomic signature because, despite diversity between species, it varies little between 50 kilobase or longer windows on a given genome. Both the causes and the functional significance of this phenomenon could be illuminated by determining if it persists on smaller scales. The profile is computed from the base step "odds ratios" that compare dinucleotide frequencies to those expected under the assumption of stochastic equilibrium (thorough shuffling). Analysis is carried out on 22 sequences, representing 19 species and comprised of about 53 million bases all together, to assess stability of the signature in windows ranging in size from 50 kilobases down to 125 bases. RESULTS: Dinucleotide relative abundance distance from the global signature is computed locally for all non-overlapping windows on each sequence. These distances are log-normally distributed with nearly constant variance and with means that tend to zero slower than reciprocal square root of window size. The mean distance within genomes is larger for protist, plant, and human chromosomes, and smaller for archaea, bacteria, and yeast, for any window size. CONCLUSIONS: The imprint of the global signature is locally pervasive on all scales considered in the sequences (either genomes or chromosomes) that were scanned.