We employed the Codonusage database to analyze the codon usage pattern in 31 organisms from all the main biological taxa. We tested the similarity in codon usage pattern between each organism and all the others by the Pearson linear correlation coefficient. The 465 values obtained were located in a 31 x 31 triangular matrix from which a correlation distance matrix was calculated. An evolutionary dendrogram was then constructed from these distances. The results showed a fair correlation between codon usage patterns and phylogenetic relationships, at least for organisms which diverged in rather recent times (end of Jurassic--beginning of Cretaceous).

We have recently reported that the statistical analysis of the frequency distribution of short oligonucleotides within mammalian and viral genomes allows the production of sets of DNA sequences enriched in signals for transcription factors. Such statistical approaches could facilitate the identification of new promoter regions playing a role in the transcriptional regulation of gene expression. In the case of mammalian oligonucleotides, we found that the published set of frequent decamers enriched in transcriptional motifs is not suitable for studies on genes of Homo sapiens and evolutionarily related genomes, because it contains decameric sequences belonging to genomic repeats. We report here that most of the decameric sequences of DNA repeats belong to Alu repeats. Accordingly, we produced a subset of Alu-free frequent decamers. In addition, we eliminated from the subset of Alu-free frequent decamers those that are frequently present within other common human repeats, including (GT)n, (AT)n, (CA)n, (ATT)n, (CAA)n and (GTT)n. The Alu-free (repeats-free) subset of frequent mammalian decamers is enriched in signals for transcription factors and allows the identification of putative signals in genes, such as those coding for plasminogen activator, adenosine deaminase and p53, that contain a large number of Alu-like repeats interspersed within our genomic sequences. The newly generated compilation of frequent decamers described here might be used to locate genomic regions playing functional roles in the expression of genes of Homo sapiens and related primates.

We studied the frequency distribution of 1,048,576 oligonucleotides 10 bp long in a sample of 1.961 Mbase of genes from plants, made of 635 sequences extracted from GenBank 71.0, with the aim of detecting transcription control signals. Among all decamers, 3255, or 0.3%, had a frequency 10 times higher than the mean and were subjected to further statistical analysis. For each of the 3255 decamers (parents), we counted the individual frequencies of the 30 decamers (progeny) differing from the parent by one base mutation, and calculated two variance/mean chi-squares for the progeny, with and without the parent decamer. By studying the distribution of the ratio between the two chi-squares we observed that out of 3255 decamers > 10 times frequent than average, 432 had a chi-square ratio > 1.9. In this residual set, which corresponds to < 0.04 per cent of all possible decamers, only 15 known eukaryotic transcription control elements were found; on the other hand, it included 29 decanucleotides that matched with decanucleotides of a set of Drosophila, 24 with a set from mammals, 13 with a set from yeast and four with a set of viruses--all sets identified with the statistical procedures here described. These decanucloetides are highly repetitive and seem to be present throughout all higher organisms, whereas they are uncommon in mammalian viruses.

We studied the frequency distribution of 1,048,576 oligonucleotides 10 bp long in a sample of 1.072 x 10(6) bases of genes from non-mammalian vertebrates, made of 322 sequences extracted from EMBL(R) 29.0, with the aim of detecting transcription control signals. Among all decamers, 2097 (0.2%) had a frequency 10 times higher than the mean and were subjected to further statistical analysis. For each of the 2097 decamers (parents), we counted the individual frequencies of the 30 decamers differing from the parent by one base mutation (progeny) and we calculated two variance/mean chi squares for the progeny, with and without the parent decamer. By studying the distribution of the ratio between the two chi squares we observed that out of 2097 decamers that occurred > 10 times more frequently than average, 1017 had a chi square ratio of between 1 and 1.5; in this final set, which corresponds to < 0.097% of all possible decamers, 75 decamers were found to contain 100 transcription control elements, like CCAAT and others. The final set contains a high excess of signals when compared to 100 random sets of 1017 decamers. Some of the decamers selected with the procedure are members of consensus sequences rather than unique sequences.

The large body of nucleic acid sequence data now available offers a unique opportunity for the characterization of individual oligonucleotides which may be specific to sequence functional domains. We have prepared algorithms for the study of the frequency distribution of all oligonucleotides of length 2-6 in DNA sequences. We have implemented them in the study of 634 mammalian DNA sequences spanning 1.782 Mb, and have obtained the distribution of the ratio between the observed frequency of oligonucleotides and their expected frequency based on independent nucleotide probabilities. We then studied the distribution of oligonucleotides (or k-tuples) of each length in a subset of 129 complete mammalian genes spanning 0.607 Mb. Eight distinct genomic regions, namely 5'-non-transcribed, first exon, first intron, intermediate exons, intermediate introns, last intron, last exon and 3'-non-transcribed, were considered. We observed that some oligonucleotides show a statistical behaviour and a regional distribution similar to that of known signal sequences. Moreover the frequency distribution of oligonucleotides of length 5 and 6 tends to become bimodal, indicating the existence of a population of very frequent oligonucleotides.