Stop Codon Usage Research Articles

Stop codons in bacteria are not selectively equivalent

BackgroundThe evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes.ResultsWe show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG.ConclusionsOur data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications.ReviewersThis article was reviewed by Michail Gelfand, Arcady Mushegian and Shamil Sunyaev. For the full reviews, please go to the Reviewers’ Comments section.

Evolution and Diversification of the Organellar Release Factor Family

Translation termination is accomplished by proteins of the Class I release factor family (RF) that recognize stop codons and catalyze the ribosomal release of the newly synthesized peptide. Bacteria have two canonical RFs: RF1 recognizes UAA and UAG, RF2 recognizes UAA and UGA. Despite that these two release factor proteins are sufficient for de facto translation termination, the eukaryotic organellar RF protein family, which has evolved from bacterial release factors, has expanded considerably, comprising multiple subfamilies, most of which have not been functionally characterized or formally classified. Here, we integrate multiple sources of information to analyze the remarkable differentiation of the RF family among organelles. We document the origin, phylogenetic distribution and sequence structure features of the mitochondrial and plastidial release factors: mtRF1a, mtRF1, mtRF2a, mtRF2b, mtRF2c, ICT1, C12orf65, pRF1, and pRF2, and review published relevant experimental data. The canonical release factors (mtRF1a, mtRF2a, pRF1, and pRF2) and ICT1 are derived from bacterial ancestors, whereas the others have resulted from gene duplications of another release factor. These new RF family members have all lost one or more specific motifs relevant for bona fide release factor function but are mostly targeted to the same organelle as their ancestor. We also characterize the subset of canonical release factor proteins that bear nonclassical PxT/SPF tripeptide motifs and provide a molecular-model-based rationale for their retained ability to recognize stop codons. Finally, we analyze the coevolution of canonical RFs with the organellar genetic code. Although the RF presence in an organelle and its stop codon usage tend to coevolve, we find three taxa that encode an RF2 without using UGA stop codons, and one reverse scenario, where mamiellales green algae use UGA stop codons in their mitochondria without having a mitochondrial type RF2. For the latter, we put forward a “stop-codon reinvention” hypothesis that involves the retargeting of the plastid release factor to the mitochondrion.

Coevolution ofaah: Adps-Like Gene with the Host Bacterium Revealed by Comparative Genomic Analysis

A protein named AAH was isolated from the bacterium Microbacterium arborescens SE14, a gut commensal of the lepidopteran larvae. It showed not only a high sequence similarity to Dps-like proteins (DNA-binding proteins from starved cell) but also reversible hydrolase activity. A comparative genomic analysis was performed to gain more insights into its evolution. The GC profile of the aah gene indicated that it was evolved from a low GC ancestor. Its stop codon usage was also different from the general pattern of Actinobacterial genomes. The phylogeny of dps-like proteins showed strong correlation with the phylogeny of host bacteria. A conserved genomic synteny was identified in some taxonomically related Actinobacteria, suggesting that the ancestor genes had incorporated into the genome before the divergence of Micrococcineae from other families. The aah gene had evolved new function but still retained the typical dodecameric structure.

Comparative Genomic and Transcriptomic Characterization of the Toxigenic Marine Dinoflagellate Alexandrium ostenfeldii

Many dinoflagellate species are notorious for the toxins they produce and ecological and human health consequences associated with harmful algal blooms (HABs). Dinoflagellates are particularly refractory to genomic analysis due to the enormous genome size, lack of knowledge about their DNA composition and structure, and peculiarities of gene regulation, such as spliced leader (SL) trans-splicing and mRNA transposition mechanisms. Alexandrium ostenfeldii is known to produce macrocyclic imine toxins, described as spirolides. We characterized the genome of A. ostenfeldii using a combination of transcriptomic data and random genomic clones for comparison with other dinoflagellates, particularly Alexandrium species. Examination of SL sequences revealed similar features as in other dinoflagellates, including Alexandrium species. SL sequences in decay indicate frequent retro-transposition of mRNA species. This probably contributes to overall genome complexity by generating additional gene copies. Sequencing of several thousand fosmid and bacterial artificial chromosome (BAC) ends yielded a wealth of simple repeats and tandemly repeated longer sequence stretches which we estimated to comprise more than half of the whole genome. Surprisingly, the repeats comprise a very limited set of 79–97 bp sequences; in part the genome is thus a relatively uniform sequence space interrupted by coding sequences. Our genomic sequence survey (GSS) represents the largest genomic data set of a dinoflagellate to date. Alexandrium ostenfeldii is a typical dinoflagellate with respect to its transcriptome and mRNA transposition but demonstrates Alexandrium-like stop codon usage. The large portion of repetitive sequences and the organization within the genome is in agreement with several other studies on dinoflagellates using different approaches. It remains to be determined whether this unusual composition is directly correlated to the exceptionally genome organization of dinoflagellates with a low amount of histones and histone-like proteins.

Adaptation of the short intergenic spacers between co-directional genes to the Shine-Dalgarno motif among prokaryote genomes

BackgroundIn prokaryote genomes most of the co-directional genes are in close proximity. Even the coding sequence or the stop codon of a gene can overlap with the Shine-Dalgarno (SD) sequence of the downstream co-directional gene. In this paper we analyze how the presence of SD may influence the stop codon usage or the spacing lengths between co-directional genes.ResultsThe SD sequences for 530 prokaryote genomes have been predicted using computer calculations of the base-pairing free energy between translation initiation regions and the 16S rRNA 3' tail. Genomes with a large number of genes with the SD sequence concentrate this regulatory motif from 4 to 11 bps before the start codon. However, not all genes seem to have the SD sequence. Genes separated from 1 to 4 bps from a co-directional upstream gene show a high SD presence, though this regulatory signal is located towards the 3' end of the coding sequence of the upstream gene. Genes separated from 9 to 15 bps show the highest SD presence as they accommodate the SD sequence within an intergenic region. However, genes separated from around 5 to 8 bps have a lower percentage of SD presence and when the SD is present, the stop codon usage of the upstream gene changes to accommodate the overlap between the SD sequence and the stop codon.ConclusionThe SD presence makes the intergenic lengths from 5 to 8 bps less frequent and causes an adaptation of the stop codon usage. Our results introduce new elements to the discussion of which factors affect the intergenic lengths, which cannot be totally explained by the pressure to compact the prokaryote genomes.

Distinct eRF3 Requirements Suggest Alternate eRF1 Conformations Mediate Peptide Release during Eukaryotic Translation Termination

Organisms that use the standard genetic code recognize UAA, UAG, and UGA as stop codons, whereas variant code species frequently alter this pattern of stop codon recognition. We previously demonstrated that a hybrid eRF1 carrying the Euplotes octocarinatus domain 1 fused to Saccharomyces cerevisiae domains 2 and 3 (Eo/Sc eRF1) recognized UAA and UAG, but not UGA, as stop codons. In the current study, we identified mutations in Eo/Sc eRF1 that restore UGA recognition and define distinct roles for the TASNIKS and YxCxxxF motifs in eRF1 function. Mutations in or near the YxCxxxF motif support the cavity model for stop codon recognition by eRF1. Mutations in the TASNIKS motif eliminated the eRF3 requirement for peptide release at UAA and UAG codons, but not UGA codons. These results suggest that the TASNIKS motif and eRF3 function together to trigger eRF1 conformational changes that couple stop codon recognition and peptide release during eukaryotic translation termination.

Ciliates use both variant and universal genetic codes: Evidence of omnipotent eRF1s in the class Litostomatea

Stop codon reassignments have occurred very frequently in ciliates. In some ciliate species, the universal stop codons UAA and UAG are translated into glutamine, while in some other species, the universal stop codon UGA appears to be translated into cysteine or tryptophan. The class Litostomatea has been hypothesized to be the only group of ciliates using the universal genetic code. However, the hypothesis was based on a statistical analysis of quite small sequence dataset which was insufficient to elucidate the codon usage of the class among such highly deviated phylum. In this study, together with the updated database sequence analysis for the class, we approached the problem of stop codon usage by examining the capacity of the translation termination factor eRF1 for recognizing stop codons. Using in vivo assay systems in budding yeast, we estimated the activity of eRF1 from two litostome species Didinium nasutum and Dileptus margaritifer. The results clearly showed that Didinium and Dileptus eRF1s efficiently recognize all three stop codons. This is the first experimental evidence that strongly supports the hypothesis that litostome ciliates use universal genetic code.

Distinct patterns of evolution between respiratory syncytial virus subgroups A and B From New Zealand isolates collected over thirty‐seven years

Respiratory syncytial virus (RSV) is the most important cause of viral lower respiratory tract infections in infants and children worldwide. In New Zealand, infants with RSV disease are hospitalized at a higher rate than other industrialized countries, without a proportionate increase in known risk factors. The molecular epidemiology of RSV in New Zealand has never been described. Therefore, we analyzed viral attachment glycoprotein (G) gene sequences from 106 RSV subgroup A isolates collected in New Zealand between 1967 and 2003, and 38 subgroup B viruses collected between 1984 and 2004. Subgroup A and B sequences were aligned separately, and compared to sequences of viruses isolated from other countries during a similar period. Genotyping and clustering analyses showed RSV in New Zealand is similar and temporally related to viruses found in other countries. By quantifying temporal clustering, we found subgroup B viruses clustered more strongly than subgroup A viruses. RSV B sequences displayed more variability in stop codon usage and predicted protein length, and had a higher degree of predicted O-glycosylation site changes than RSV A. The mutation rate calculated for the RSV B G gene was significantly higher than for RSV A. Together, these data reveal that RSV subgroups exhibit different patterns of evolution, with subgroup B viruses evolving faster than A.

Genetic heterogeneity of G and F protein genes from Argentinean human metapneumovirus strains

Human metapneumovirus (hMPV) is a newly identified paramixovirus, associated with respiratory illnesses in all age groups. Two genetic groups of hMPV have been described. The nucleotide sequences of the G and F genes from 11 Argentinean hMPV strains (1998-2003) were determined by RT-PCR and direct sequencing. Phylogenetic analysis showed that hMPV strains clustered into two main genetic lineages, A and B. Strains clustered into A group were split into two sublineages, A1 and A2. All strains belonging to group B clustered with representative strains from sublineage B1. No Argentinean strains belonged to sublineage B2. F sequences showed high percentage identities at nucleotide and amino acid levels. In contrast, G sequences showed high diversity between A and B groups. Most changes observed in the deduced G protein sequence were amino acid substitutions in the extracellular domain, and changes in stop codon usage leading to different lengths in the G proteins. High content of serine and threonine residues were also shown, suggesting that this protein would be highly glycosylated. The potential sites for N- and O-glycosylation seem to have a different conservation pattern between the two main groups. This is the first report on the genetic variability of the G and F protein genes of hMPV strains in South America. Two main genetic groups and at least three subgroups were revealed among Argentinean hMPV strains. The F protein seems to be highly conserved, whereas the G protein showed extensive diversity between groups A and B.

Are stop codons recognized by base triplets in the large ribosomal RNA subunit?

The precise mechanism of stop codon recognition in translation termination is still unclear. A previously published study by Ivanov and colleagues proposed a new model for stop codon recognition in which 3-nucleotide Ter-anticodons within the loops of hairpin helices 69 (domain IV) and 89 (domain V) in large ribosomal subunit (LSU) rRNA recognize stop codons to terminate protein translation in eubacteria and certain organelles. We evaluated this model by extensive bioinformatic analysis of stop codons and their putative corresponding Ter-anticodons across a much wider range of species, and found many cases for which it cannot explain the stop codon usage without requiring the involvement of one or more of the eight possible noncomplementary base pairs. Involvement of such base pairs may not be structurally or thermodynamically damaging to the model. However, if, according to the model, Ter-anticodon interaction with stop codons occurs within the ribosomal A-site, the structural stringency which that site imposes on sense codon.tRNA anticodon interaction should also extend to stop codon.Ter-anticodon interactions. Moreover, with Ter-tRNA in place of an aminoacyl-tRNA, for each of the various Ter-anticodons there is a sense codon that can interact with it preferentially by complementary and wobble base-pairing. Both these considerations considerably weaken the arguments put forth previously.

Relationships Among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

It is well known that stop codons play a critical role in the process of protein synthesis. However, little effort has been made to investigate whether stop codon usage exhibits biases, such as widely seen for synonymous codon usage. Here we systematically investigate stop codon usage bias in various eukaryotes as well as its relationships with its context, GC3 content, gene expression level, and secondary structure. The results show that there is a strong bias for stop codon usage in different eukaryotes, i.e., UAA is overrepresented in the lower eukaryotes, UGA is overrepresented in the higher eukaryotes, and UAG is least used in all eukaryotes. Different conserved patterns for each stop codon in different eukaryotic classes are found based on information content and logo analysis. GC3 contents increase with increasing complexity of organisms. Secondary structure prediction revealed that UAA is generally associated with loop structures, whereas UGA is more uniformly present in loop and stem structures, i.e., UGA is less biased toward having a particular structure. The stop codon usage bias, however, shows no significant relationship with GC3 content and gene expression level in individual eukaryotes. The results indicate that genomic complexity and GC3 content might contribute to stop codon usage bias in different eukaryotes. Our results indicate that stop codons, like synonymous codons, exhibit biases in usage. Additional work will be needed to understand the causes of these biases and their relationship to the mechanism of protein termination.

Genetic and Antigenic Variability of Human Respiratory Syncytial Virus (Groups A and B) Isolated over Seven Consecutive Seasons in Argentina (1995 to 2001)

The genetic and antigenic variability of human respiratory syncytial virus (HRSV) strains isolated in Buenos Aires from 1995 to 2001 was evaluated by partial nucleotide sequencing of the G gene and enzyme-linked immunosorbent assay analysis with anti-G monoclonal antibodies. Phylogenetic analyses showed that 37 group A strains clustered into five genotypes, whereas 20 group B strains clustered into three genotypes. Group A showed more genetic variability than group B. A close correlation between genotypes and antigenic patterns was observed. Changes detected in the G protein of viruses from both groups included (i) amino acid substitutions and(ii) differences in protein length due to either changes in stop codon usage or sequence duplications. Three B strains from 1999 exhibited a duplication of 20 amino acids, while one B strain from 2001 had 2 amino acids duplicated. The comparison among Argentinean HRSV strains and viruses isolated in other geographical areas during different epidemics is discussed.

Unusual tubulins with deviant genetic code in the hypotrich ciliate Euplotidium itoi

We characterized nine macronuclear tubulin gene minichromosomes in the hypotrich ciliate Euplotidium itoi: two α‐tubulin genes, five β‐tubulin genes, and two γ‐tubulin genes. In particular, one peculiar β‐tubulin gene was identified. It shows only 81.5% amino acid similarity to the other characterized β‐tubulin types and an UAG stop codon in frame. Based on this result we started to analyse stop codon usage in Euplotidium itoi macronuclear tubulin genes and some unrelated macronuclear minichromosomes. Most of the ciliates independently evolved diverse non‐canonical genetic codes. Differences to the standard code consist in the reassignment of one or two of the three standard stop codons to encode a particular amino acid. Class Spirotrichea represents a particular case comprising two groups of organisms that do not share any stop codon. In oxytrichids, UAR is translated as glutamine and UGA is the single stop codon. This is in contrast to the genus Euplotes, close relatives of Euplotidium, in which UAR encodes a stop signal and UGA is translated as cysteine. Fourteen complete genes were characterized in Euplotidium itoi; in all cases UAA revealed as the single triplet used to code a stop, UGA was never and UAG rarely observed in frame. The UAG stop codon was observed in the peculiar type 3 β‐tubulin gene, in two dynein partial sequences and in some other putative genes. We could confirm the active expression of three of these genes, the unusual β‐tubulin, a putative lipase and a possible steroid‐binding protein, by RT‐PCR. Finally, a detailed analysis on the structure of the macronuclear tubulin genes, including nucleotide and amino acid differences among tubulin gene paralogs, polyadenylation sites, putative chromosome fragmentation sites, ATGC‐content in coding and non‐coding regions, telomere sequence and length will be presented.

Newly sequenced eRF1s from ciliates: the diversity of stop codon usage and the molecular surfaces that are important for stop codon interactions

The genetic code of nuclear genes in some ciliates was found to differ from that of other organisms in the assignment of UGA, UAG, and UAA codons, which are normally assigned as stop codons. In some ciliate species, the universal stop codons UAA and UAG instead encode glutamine. In some other ciliates, the universal stop codon UGA appears to be translated as cysteine or tryptophan. Eukaryotic release factor 1 (eRF1) is a key protein in stop codon recognition, thus, the protein is believed to play an important role in the stop codon reassignment in ciliates. We have cloned, sequenced, and analyzed the cDNA of eRF1 from four ciliate species of three different classes: Karyorelictea ( Loxodes striatus), Heterotrichea ( Blepharisma musculus), and Litostomatea ( Didinium nasutum, Dileptus margaritifer). Phylogenetic analysis of these eRF1s supports the hypothesis that the genetic code in ciliates has deviated independently several times from the universal genetic code, and that different ciliate eRF1s may have undergone different processes to change the codon specificity. Using computational methods, we have also suggested areas on the surface of eRF1s that are important for stop codon recognition in ciliate eRF1s.

Different stop codon usage in two pseudohypotrich ciliates

Based on rRNA phylogeny, morphologic and morphogenetic characters, two major groups of hypotrich ciliates can be distinguished: euhypotrichs and pseudohypotrichs. Through the sequencing of actin genes, we show here that, interestingly, the pseudohypotrichs Dyophrys sp. and Euplotes vannus have a different stop codon usage. In fact, the stop codon usage of the former species resembles that of euhypotrichs. This unexpected result is used to discuss the origin and acquisition of genetic code deviations in ciliates.

Isolation, Cloning, and Complete Nucleotide Sequence of a Phenotypically Distinct Brazilian Isolate of Human T-Lymphotropic Virus Type II (HTLV-II)

Analysis of human T-lymphotropic virus type II (HTLV-II) isolates from North America and Europe have demonstrated the existence of two molecular subtypes of the virus, HTLV-IIa and HTLV-IIb. Recently, studies on HTLV-II infections in Brazil have revealed isolates that are related phylogenetically to the HTLV-IIa subtype but have a HTLV-IIb phenotype with respect to the transactivating protein, tax. To more clearly define this relationship, HTLV-II was isolated from peripheral blood of an IVDA from Sao Paulo, Brazil (SP-WV), and the complete provirus was cloned and sequenced. Comparison of HTLV-IISP-WV nucleotide sequences to other available complete HTLV-II proviral sequences revealed that HTLV-IISP-WV is most closely related to HTLV-IIMo, the prototypic HTLV-IIa subtype sequence. Phylogenetic analysis of LTR, env, and tax regions unequivocally demonstrated that HTLV-IISP-WV and all other Brazilian sequences examined are members of the IIa subtype. The predicted amino acid sequences of the major coding regions of HTLV-IISP-WV are also most closely related to HTLV-IIMo, with the important exception of tax. The tax protein encoded by HTLV-IISP-WV is 96–99% identical to the tax of IIb isolates and is similar in that it has an additional 25 amino acids at the carboxy-terminus compared to the HTLV-IIMo tax with which it shares 91% identity. Analysis of tax stop codon usage of a number of HTLV-IIa isolates from North American, Europe, and Brazil demonstrated that isolates from the last region appear to be unique in their extended tax phenotype. It could be demonstrated that the extended tax proteins in the HTLV-IIb and Brazilian isolates had equivalent ability to transactivate the viral LTR, and studies with deletion mutants indicated that the extended C-terminus is not essential for transactivation. In contrast, the HTLV-IIa tax was found to have a greatly diminished ability to transactivate the viral LTR, which appeared to be a consequence of reduced expression of the protein. The studies show that although the Brazilian strains do not represent an entirely new subtype based on nucleotide sequence analysis they are a phenotypically unique molecular variant within the HTLV-IIa subtype.

Translation in Bacillus subtilis: roles and trends of initiation and termination, insights from a genome analysis.

We analysed the Bacillus subtilis protein coding sequences termini, and compared it to other genomes. The analysis focused on signals, com-positional biases of nucleotides, oligonucleotides, codons and amino acids and mRNA secondary structure. AUG is the preferred start codon in all genomes, independent of their G+C content, and seems to induce less stable mRNA structures. However, it is not conserved between homologous genes neither is it preferred in highly expressed genes. In B.subtilis the ribosome binding site is very strong. We found that downstream boxes do not seem to exist either in Escherichia coli or in B.subtilis. UAA stop codon usage is correlated with the G+C content and is strongly selected in highly expressed genes. We found less stable mRNA structures at both termini, which we related to mRNA-ribosome and mRNA-release-factor interactions. This pattern seems to impose a peculiar A-rich nucleotide and codon usage bias in these regions. Finally the analysis of all proteins from B.subtilis revealed a similar amino acid bias near both termini of proteins consisting of over-representation of hydrophilic residues. This bias near the stop codon is partially release-factor specific.

Function of 3' non-coding sequences and stop codon usage in expression of the chloroplast psaB gene in Chlamydomonas reinhardtii.

The rate of mRNA decay is an important step in the control of gene expression in prokaryotes, eukaryotes and cellular organelles. Factors that determine the rate of mRNA decay in chloroplasts are not well understood. Chloroplast mRNAs typically contain an inverted repeat sequence within the 3' untranslated region that can potentially fold into a stem-loop structure. These stem-loop structures have been suggested to stabilize the mRNA by preventing degradation by exonuclease activity, although such a function in vivo has not been clearly established. Secondary structures within the translation reading frame may also determine the inherent stability of an mRNA. To test the function of the inverted repeat structures in chloroplast mRNA stability mutants were constructed in the psaB gene that eliminated the 3' flanking sequences of psaB or extended the open reading frame into the 3' inverted repeat. The mutant psaB genes were introduced into the chloroplast genome of Chlamydomonas reinhardtii. Mutants lacking the 3' stem-loop exhibited a 75% reduction in the level of psaB mRNA. The accumulation of photosystem I complexes was also decreased by a corresponding amount indicating that the mRNA level is limiting to PsaB protein synthesis. Pulse-chase labeling of the mRNA showed that the decay rate of the psaB mRNA was significantly increased demonstrating that the stem-loop structure is required for psaB mRNA stability. When the translation reading frame was extended into the 3' inverted repeat the mRNA level was reduced to only 2% of wild-type indicating that ribosome interaction with stem-loop structures destabilizes chloroplast mRNAs. The non-photosynthetic phenotype of the mutant with an extended reading frame allowed us to test whether infrequently used stop codons (UAG and UGA) can terminate translation in vivo. Both UAG and UGA are able to effectively terminate PsaB synthesis although UGA is never used in any of the Chlamydomonas chloroplast genes that have been sequenced.

The translational termination signal database (TransTerm) now also includes initiation contexts.

The TransTerm database of termination codon contexts has been extended to include sense codon usage, and initiation codon contexts. The database was constructed from 23,721 coding sequences from 93 organisms. The database contains: a) the sequence around the termination codon (-10, +10); b) the sequence around the initiation codon (-20, +10); c) the length, 'G+C%' of the third position of codons (GC3), the 'codon adaptation index' (CAI) and the 'effective number of codons' statistic (Nc); d) summary tables for each organism including total codon usage, stop codon and tetranucleotide stop-signal usage, and matrices tallying base frequencies at each position around the initiation and termination codons. The data are arranged to facilitate investigation of the relationships between the three phases of protein synthesis. The database is available electronically from EMBL.

The translational termination signal database.

The Translational Termination Database (TransTerm) consists of the immediate context sequences around the natural termination codons from 45 organisms, and summary tables. The influence of termination codon context on their effectivness as stop signals has been widely documented. The SPECIES--TRI.DAT table shows trinucleotide stop codon usage in each organism and for comparison the occurrence of these sequences in the noncoding region. The SPECIES--TETRA.DAT table contains is a similar table of tetranucleotide stop signal usage. The database is available from EMBL.