On the recombinational origin of protein-sequence-subunit structure

Escherichia coli genome is composed of two distinct types of nucleotide sequences

We calculated correlations of the nucleotide distributions along the E. coli genome. Subsequent cluster analysis of the correlation distributions showed that the genome was composed of two qualitatively different types of nucleotide sequences. The first type exhibited strong correlations of the genomic distributions of A with T and G with C, and high anticorrelations of A with C and G with T. In contrast, the second type was characterized by weak or negligible correlations typical of randomized sequences. Both types of sequences were almost equally abundant in the E. coli genome and their length varied from several hundred nucleotides to about 70 kilobases. They were not disjunct with respect to their (G + C) content but the high correlations and anticorrelations were rather characteristic for (A + T)-rich genomic segments. We offer possible explanations of the mosaic structure of the E. coli genome.

Variations of the mononucleotide and short oligonucleotide distributions in the genomes of various organisms

We calculated the variation coefficients of the mononucleotide and short oligonucleotide distributions in over 1700 long genomic sequences originating from six organisms to demonstrate that the human and Escherichia coli genomic sequences were the least and the most uniform, respectively. The most non-random genomic distributions were exhibited by the four canonical nucleotides, followed by the strong and weak nucleotides, while the distributions of purine or pyrimidine nucleotides and especially the distributions of (A+C) and (G+T) were significantly more uniform even in the human genome. In the human and mouse genomes, the highest coefficients of variation were further observed with the oligonucleotides where CG was combined with the strong nucleotides while its combination with the weak nucleotides significantly decreased the variation which, however, was still very high. High variation was also exhibited by the remaining oligonucleotides composed exclusively of the strong nucleotides or those containing only weak nucleotides. On the other hand, the distributions of oligonucleotides containing similar and especially the same numbers of the strong and weak nucleotides, but no CG or TA dinucleotide, were the most uniform. The information following from the present analysis will be useful not only in the identification of important genomic regions but also in computer simulations of the genomic nucleotide sequences in order to trace and reproduce the pathways of genome evolution.

Segmented structure of separate and transposable DNA and RNA elements as suggested by their size distributions

A collection of about 1000 different eukaryotic and prokaryotic DNA mobile and separate elements is compiled from literature-transposons, plasmids, extrachromosomal circular DNA, insertion sequences, as well as viral genomes and separate genome segments. Only small elements are collected, upto 2000 base pairs. Analysis of the sequence length distributions of the elements reveals that certain sizes are clearly preferred, namely those which correspond to multiples of about 345 bp in eukaryotes and multiples of about 210 bp in prokaryotes. This provides additional evidence in support of the theory (1) that segmented structure is characteristic of not only protein-coding sequences (2) but rather of genomes in general. In particular, it confirms the prediction (1) that mobile and separate elements would also be segmented.

Segmented structure of protein sequences and early evolution of genome by combinatorial fusion of DNA elements

A theory of an early stage of genome evolution by combinatorial fusion of circular DNA units is suggested, based on protein sequence "fossil" evidence. The evidence includes preference of protein sequence lengths for certain sizes--multiples of 123 aa for eukaryotes and multiples of 152 aa for prokaryotes. At the DNA level these sizes correspond to 350-450 base pairs--the known optimal range for DNA ring closure. The methionine residues repeatedly appear along the sequences with the same period of about 120 aa (in eukaryotes), presumably marking the sites of insertion of the early genes--rings of protein-coding DNA. No torsional constraint in this DNA results in very sharp estimate of the helical periodicity of the early DNA, indistinguishable from the experimental mean value for extant DNA. According to the combinatorial fusion theory, based on the above evidence, in the pregenomic, prerecombinational stage the genes and the noncoding sequences existed in form of autonomously replicating DNA rings of close to standard size, randomly segregating between dividing cells, like modern plasmids do. In the recombinational early genomic stage the rings started to fuse, forming larger DNA molecules consisting of several unit genes connected in various combinations and forming long protein-coding sequences (combinatorial fusion). This process, which involved, perhaps, noncoding sequences as well, eventually resulted in the formation of large genomes. The dispersed circular DNA--or, rather, evolutionarily advanced derivatives thereof--may still exist in the form of various mobile DNA elements.

Periodic recurrence of methionines: fossil of gene fusion?

As we have recently shown, approximately 20% of proteins are made of uniform size units of approximately 123 aa for eukaryotes and approximately 152 aa for prokaryotes. Such regularity may reflect certain past events in protein evolution by fusion (molecular recombination) of a spectrum of standard-size protein-coding DNA segments--the early genes. Consequently, methionines, as start residues, would mark those locations in proteins that correspond to the DNA recombination sites--the borders between the fused genes. This positional preference of the methionines may still survive as a fossil of the early protein sequence organization. In this study we address the question how methionines are distributed in modern protein sequences. This analysis of eukaryotic sequences shows that methionine residues do preferentially appear at the positions corresponding to the multiples of the unit size, as predicted.