Characterisation of a non-canonical genetic code in the oxymonad Streblomastix strix

The genetic code is one of the most highly conserved characters in living organisms. Only a small number of genomes have evolved slight variations on the code, and these non-canonical codes are instrumental in understanding the selective pressures maintaining the code. Here, we describe a new case of a non-canonical genetic code from the oxymonad flagellate Streblomastix strix. We have sequenced four protein-coding genes from S.strix and found that the canonical stop codons TAA and TAG encode the amino acid glutamine. These codons are retained in S.strix mRNAs, and the legitimate termination codons of all genes examined were found to be TGA, supporting the prediction that this should be the only true stop codon in this genome. Only four other lineages of eukaryotes are known to have evolved non-canonical nuclear genetic codes, and our phylogenetic analyses of alpha-tubulin, beta-tubulin, elongation factor-1 alpha (EF-1 alpha), heat-shock protein 90 (HSP90), and small subunit rRNA all confirm that the variant code in S.strix evolved independently of any other known variant. The independent origin of each of these codes is particularly interesting because the code found in S.strix, where TAA and TAG encode glutamine, has evolved in three of the four other nuclear lineages with variant codes, but this code has never evolved in a prokaryote or a prokaryote-derived organelle. The distribution of non-canonical codes is probably the result of a combination of differences in translation termination, tRNAs, and tRNA synthetases, such that the eukaryotic machinery preferentially allows changes involving TAA and TAG.     

Environmental PCR survey to determine the distribution of a non-canonical genetic code in uncultivable oxymonads

The universal genetic code is conserved throughout most living systems, but a non-canonical code where TAA and TAG encode glutamine has evolved in several eukaryotes, including oxymonad protists. Most oxymonads are uncultivable, so environmental RT-PCR and PCR was used to examine the distribution of this rare character. A total of 253 unique isolates of four protein-coding genes were sampled from the hindgut community of the cockroach, Cryptocercus punctulatus , an environment rich in diversity from two of the five subgroups of oxymonad, saccinobaculids and polymastigids. Four α-tubulins were found with non-canonical glutamine codons. Environmental RACE confirmed that these and related genes used only TGA as stop codons, as expected for the non-canonical code, whereas other genes used TAA or TAG as stop codons, as expected for the universal code. We characterized α-tubulin from manually isolated Saccinobaculus ambloaxostylus , confirming it uses the universal code and suggesting, by elimination, that the non-canonical code is used by a polymastigid. HSP90 and EF-1α phylogenies also showed environmental sequences falling into two distinct groups, and are generally consistent with previous hypotheses that polymastigids and Streblomastix are closely related. Overall, we propose that the non-canonical genetic code arose once in a common ancestor of Streblomastix and a subgroup of polymastigids.     

A non-canonical genetic code in an early diverging eukaryotic lineage.

The nearly invariant nature of the ''Universal Genetic Code'' attests to its early establishment in evolution and to the difficulty of altering it now, since so many molecules are required for, and depend upon, faithful translation. Nevertheless, variations on the universal code are known in a handful of genomes. We have found one such variant in diplomonads, an early-diverging eukaryotic lineage. Genes for alpha-tubulin, beta-tubulin and elongation factor 1 alpha (EF-1alpha) from two unclassified strains of Hexamitidae were found to contain TAA and TAG (TAR) triplets at positions suggesting a variant code in which TAR codes for glutamine. We found confirmation of this hypothesis by identifying genes encoding glutamine-tRNAs with CUA and UUA anticodons. The alpha-tubulin and EF-1alpha genes from two other diplomonads, Spironucleus muris and Hexamita inflata, were also sequenced and shown to contain no such non-canonical codons. However, tRNA genes with the anticodons UUA and CUA were found in H.inflata, suggesting that this diplomonad also uses these codons, albeit infrequently. The high GC content of these genomes and the presence of two isoaccepting tRNAs compound the difficulty of understanding how this variant code arose by strictly neutral means.     

Consequences of stop codon reassignment on protein evolution in ciliates with alternative genetic codes

Tetrahymena thermophila and Paramecium tetraurelia are ciliates that reassign TAA and TAG from stop codons to glutamine codons. Because of the lack of full genome sequences, few studies have concentrated on analyzing the effects of codon reassignment in protein evolution. We used the recently sequenced genome of these species to analyze the patterns of amino acid substitution in ciliates that reassign the code. We show that, as expected, the codon reassignment has a large impact on amino acid substitutions in closely related proteins; however, contrary to expectations, these effects also hold for very diverged proteins. Previous studies have used amino acid substitution data to calculate the minimization of the genetic code; our results show that because of the lasting influence of the code in the patterns of substitution, such studies are tautological. These different substitution patterns might affect alignment of ciliate proteins, as alignment programs use scoring matrices based on substitution patterns of organisms that use the standard code. We also show that glutamine is used more frequently in ciliates than in other species, as often as expected based on the presence of the 2 new reassigned codons, indicating that the frequencies of amino acids in proteomes is mostly determined by neutral processes based on their number of codons.     

Complex phylogenetic distribution of a non-canonical genetic code in green algae

Background  

A non-canonical nuclear genetic code, in which TAG and TAA have been reassigned from stop codons to glutamine, has evolved independently in several eukaryotic lineages, including the ulvophycean green algal orders Dasycladales and Cladophorales. To study the phylogenetic distribution of the standard and non-canonical genetic codes, we generated sequence data of a representative set of ulvophycean green algae and used a robust green algal phylogeny to evaluate different evolutionary scenarios that may account for the origin of the non-canonical code.     

Tandem Stop Codons in Ciliates That Reassign Stop Codons

Tandem stop codons are extra stop codons hypothesized to be present downstream of genes to act as a backup in case of read-through of the real stop codon. Although seemingly absent from Escherichia coli, recent studies have confirmed the presence of such codons in yeast. In this paper we will analyze the genomes of two ciliate species—Paramecium tetraurelia and Tetrahymena thermophila—that reassign the stop codons TAA and TAG to glutamine, for the presence of tandem stop codons. We show that there are more tandem stop codons downstream of both Paramecium and Tetrahymena genes than expected by chance given the base composition of the downstream regions. This excess of tandem stop codons is larger in Tetrahymena and Paramecium than in yeast. We propose that this might be caused by a higher frequency of stop codon read-through in these species than in yeast, possibly because of a leaky termination machinery resulting from stop codon reassignment.     

A relationship between GC content and coding-sequence length

Since base composition of translational stop codons (TAG, TAA, and TGA) is biased toward a low G+C content, a differential density for these termination signals is expected in random DNA sequences of different base compositions. The expected length of reading frames (DNA segments of sense codons flanked by in-phase stop codons) in random sequences is thus a function of GC content. The analysis of DNA sequences from several genome databases stratified according to GC content reveals that the longest coding sequences—exons in vertebrates and genes in prokaryotes—are GC-rich, while the shortest ones are GC-poor. Exon lengthening in GC-rich vertebrate regions does not result, however, in longer vertebrate proteins, perhaps because of the lower number of exons in the genes located in these regions. The effects on coding-sequence lengths constitute a new evolutionary meaning for compositional variations in DNA GC content. Correspondence to: J. L. Oliver     

The Distribution of Elongation Factor-1 Alpha (EF-1α), Elongation Factor-Like (EFL), and a Non-Canonical Genetic Code in the Ulvophyceae: Discrete Genetic Characters Support a Consistent Phylogenetic Framework

ABSTRACT. The systematics of the green algal class Ulvophyceae have been difficult to resolve with ultrastructural and molecular phylogenetic analyses. Therefore, we investigated relationships among ulvophycean orders by determining the distribution of two discrete genetic characters previously identified only in the order Dasycladales. First, Acetabularia acetabulum uses the core translation GTPase Elongation Factor 1α (EF-1α) while most Chlorophyta instead possess the related GTPase Elongation Factor-Like (EFL). Second, the nuclear genomes of dasycladaleans A. acetabulum and Batophora oerstedii use a rare non-canonical genetic code in which the canonical termination codons TAA and TAG instead encode glutamine. Representatives of Ulvales and Ulotrichales were found to encode EFL, while Caulerpales, Dasycladales, Siphonocladales, and Ignatius tetrasporus were found to encode EF-1α, in congruence with the two major lineages previously proposed for the Ulvophyceae. The EF-1α of I. tetrasporus supports its relationship with Caulerpales/Dasycladales/Siphonocladales, in agreement with ultrastructural evidence, but contrary to certain small subunit rRNA analyses that place it with Ulvales/Ulotrichales. The same non-canonical genetic code previously described in A. acetabulum was observed in EF-1α sequences from Parvocaulis pusillus (Dasycladales), Chaetomorpha coliformis , and Cladophora cf. crinalis (Siphonocladales), whereas Caulerpales use the universal code. This supports a sister relationship between Siphonocladales and Dasycladales and further refines our understanding of ulvophycean phylogeny.     

The signal for the termination of protein synthesis in procaryotes.

The sequences around the stop codons of 862 Escherichia coli genes have been analysed to identify any additional features which contribute to the signal for the termination of protein synthesis. Highly significant deviations from the expected nucleotide distribution were observed, both before and after the stop codon. Immediately prior to UAA stop codons in E. coli there is a preference for codons of the form NAR (any base, adenine, purine), and in particular those that code for glutamine or the basic amino acids. In contrast, codons for threonine or branched nonpolar amino acids were under-represented. Uridine was over-represented in the nucleotide position immediately following all three stop codons, whereas adenine and cytosine were under-represented. This pattern is accentuated in highly expressed genes, but is not as marked in either lowly expressed genes or those that terminate in UAG, the codon specifically recognised by polypeptide chain release factor-1. These observations suggest that for the efficient termination of protein synthesis in E. coli, the 'stop signal' may be a tetranucleotide, rather than simply a tri-nucleotide codon, and that polypeptide chain release factor-2 recognises this extended signal. The sequence following stop codons was analysed in genes from several other procaryotes and bacteriophages. Salmonella typhimurium, Bacillus subtilis, bacteriophages and the methanogenic archaebacteria showed a similar bias to E. coli.     

The complete mitogenome of the Chinese bush cricket, Gampsocleis gratiosa （Orthoptera： Tettigonioidea）

The complete mitochondrial genome (mitogenome) of Gampsocleis gratiosa was determined. The 15, 929 bp in the size of G. gratiosa mitogenome contains a typical gene content, base composition, and codon usage found in metazoan. All 13 protein coding genes (PCGs) of the G. gratiosa mitogenome start with a typical ATN codon. The usual termination codons (TAA and TAG) were found from 10 PCGs. However, the atp6, nad4, and nad5 had incomplete termination codon (T). The anticodons of all tRNAs are identical to those observed in Drosophila yakuba and Locusta migratoria, and can be folded in the form of a typical clover leaf structure except for trnS (AGN). The secondary structure of trnS (AGN) was drawn according with the Steinberg-Cedergren tertiary structure. The A T content (67.4%) of the A T-rich region is relatively lower among the mitogenome regions, in contrast, it usually contains the highest A T content for most insects. Two isolated sequence repeat regions (202 bp) were found in the A T-rich region with mapping and secondary structure information.     

Role of premature stop codons in bacterial evolution

When the stop codons TGA, TAA, and TAG are found in the second and third reading frames of a protein-encoding gene, they are considered premature stop codons (PSC). Deinococcus radiodurans disproportionately favored TGA more than the other two triplets as a PSC. The TGA triplet was also found more often in noncoding regions and as a stop codon, though the bias was less pronounced. We investigated this phenomenon in 72 bacterial species with widely differing chromosomal GC contents. Although TGA and TAG were compositionally similar, we found a great variation in use of TGA but a very limited range of use of TAG. The frequency of use of TGA in the gene sequences generally increased with the GC content of the chromosome, while the frequency of use of TAG, like that of TAA, was inversely proportional to the GC content of the chromosome. The patterns of use of TAA, TGA and TAG as real stop codons were less biased and less influenced by the GC content of the chromosome. Bacteria with higher chromosomal GC contents often contained fewer PSC trimers in their genes. Phylogenetically related bacteria often exhibited similar PSC ratios. In addition, metabolically versatile bacteria have significantly fewer PSC trimers in their genes. The bias toward TGA but against TAG as a PSC could not be explained either by the preferential usage of specific codons or by the GC contents of individual chromosomes. We proposed that the quantity and the quality of the PSC in the genome might be important in bacterial evolution.     

Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes

BACKGROUND: Lateral gene transfer (LGT) is an important evolutionary mechanism among prokaryotes. The situation in eukaryotes is less clear; the human genome sequence failed to give strong support for any recent transfers from prokaryotes to vertebrates, yet a number of LGTs from prokaryotes to protists (unicellular eukaryotes) have been documented. Here, we perform a systematic analysis to investigate the impact of LGT on the evolution of diplomonads, a group of anaerobic protists.RESULTS: Phylogenetic analyses of 15 genes present in the genome of the Atlantic Salmon parasite Spironucleus barkhanus and/or the intestinal parasite Giardia lamblia show that most of these genes originated via LGT. Half of the genes are putatively involved in processes related to an anaerobic lifestyle, and this finding suggests that a common ancestor, which most probably was aerobic, of Spironucleus and Giardia adapted to an anaerobic environment in part by acquiring genes via LGT from prokaryotes. The sources of the transferred diplomonad genes are found among all three domains of life, including other eukaryotes. Many of the phylogenetic reconstructions show eukaryotes emerging in several distinct regions of the tree, strongly suggesting that LGT not only involved diplomonads, but also involved other eukaryotic groups.CONCLUSIONS: Our study shows that LGT is a significant evolutionary mechanism among diplomonads in particular and protists in general. These findings provide insights into the evolution of biochemical pathways in early eukaryote evolution and have important implications for studies of eukaryotic genome evolution and organismal relationships. Furthermore, "fusion" hypotheses for the origin of eukaryotes need to be rigorously reexamined in the light of these results.     

The mitochondrial genome of Apodemus peninsulae (Rodentia, Muridae)

The complete mitochondrial (mt) genome of the Korean field mouse Apodemus peninsulae was sequenced and found to be 16,266 bp in length. The mt protein-coding genes of A. peninsulae had ATG, GTG, ATC, and ATA as initiation codons and TAA, TAG, TA, and T as termination codons. Two forms each of trnL and trnS and the three tRNA clusters, IQM, WANCY, and HSL were identified, as in the typical Rodentia mt genome. Among tRNAs, abnormal cloverleaf structure of trnS((AGY)) was identified in DHU arm. The l-strand replication origin has the potential to form a stable stem-loop structure and control region has several conserved sequence elements.     

High mannose-binding lectin with preference for the cluster of alpha1-2-mannose from the green alga Boodlea coacta is a potent entry inhibitor of HIV-1 and influenza viruses

The complete amino acid sequence of a lectin from the green alga Boodlea coacta (BCA), which was determined by a combination of Edman degradation of its peptide fragments and cDNA cloning, revealed the following: 1) B. coacta used a noncanonical genetic code (where TAA and TAG codons encode glutamine rather than a translation termination), and 2) BCA consisted of three internal tandem-repeated domains, each of which contains the sequence motif similar to the carbohydrate-binding site of Galanthus nivalis agglutinin-related lectins. Carbohydrate binding specificity of BCA was examined by a centrifugal ultrafiltration-HPLC assay using 42 pyridylaminated oligosaccharides. BCA bound to high mannose-type N-glycans but not to the complex-type, hybrid-type core structure of N-glycans or oligosaccharides from glycolipids. This lectin had exclusive specificity for α1-2-linked mannose at the nonreducing terminus. The binding activity was enhanced as the number of terminal α1-2-linked mannose substitutions increased. Mannobiose, mannotriose, and mannopentaose were incapable of binding to BCA. Thus, BCA preferentially recognized the nonreducing terminal α1-2-mannose cluster as a primary target. As predicted from carbohydrate-binding propensity, this lectin inhibited the HIV-1 entry into the host cells at a half-maximal effective concentration of 8.2 nm. A high association constant (3.71 × 10(8) M(-1)) of BCA with the HIV envelope glycoprotein gp120 was demonstrated by surface plasmon resonance analysis. Moreover, BCA showed the potent anti-influenza activity by directly binding to viral envelope hemagglutinin against various strains, including a clinical isolate of pandemic H1N1-2009 virus, revealing its potential as an antiviral reagent.     

Cloning and sequence comparison of the mouse, human, and chicken engrailed genes reveal potential functional domains and regulatory regions.

We have isolated and characterized genomic DNA clones for the human and chicken homologues of the mouse En-1 and En-2 genes and determined the genomic structure and predicted protein sequences of both En genes in all three species. Comparison of these vertebrate En sequences with the Xenopus En-2 [Hemmati-Brivanlou et al., 1991) and invertebrate engrailed-like genes showed that the two previously identified highly conserved regions within the En protein ]reviewed in Joyner and Hanks, 1991] can be divided into five distinct subregions, designated EH1 to EH5. Sequences 5' and 3' to the predicted coding regions of the vertebrate En genes were also analyzed in an attempt to identify cis-acting DNA sequences important for the regulation of En gene expression. Considerable sequence similarity was found between the mouse and human homologues both within the putative 5' and 3' untranslated as well as 5' flanking regions. Between the mouse and Xenopus En-2 genes, shorter stretches of sequence similarity were found within the 3' untranslated region. The 5' untranslated regions of the mouse, chicken and Xenopus En-2 genes, however, showed no similarly conserved stretches. In a preliminary analysis of the expression pattern of the human En genes, En-2 protein and RNA were detected in the embryonic and adult cerebellum respectively and not in other tissues tested. These patterns are analogous to those seen in other vertebrates. Taken together these results further strengthen the suggestion that En gene function and regulation has been conserved throughout vertebrate evolution and, along with the five highly conserved regions within the En protein, raise an interesting question about the presence of conserved genetic pathways.     

Variation in Release Factor Abundance Is Not Needed to Explain Trends in Bacterial Stop Codon Usage

In bacteria stop codons are recognized by one of two class I release factors (RF1) recognizing TAG, RF2 recognizing TGA, and TAA being recognized by both. Variation across bacteria in the relative abundance of RF1 and RF2 is thus hypothesized to select for different TGA/TAG usage. This has been supported by correlations between TAG:TGA ratios and RF1:RF2 ratios across multiple bacterial species, potentially also explaining why TAG usage is approximately constant despite extensive variation in GC content. It is, however, possible that stop codon trends are determined by other forces and that RF ratios adapt to stop codon usage, rather than vice versa. Here, we determine which direction of the causal arrow is the more parsimonious. Our results support the notion that RF1/RF2 ratios become adapted to stop codon usage as the same trends, notably the anomalous TAG behavior, are seen in contexts where RF1:RF2 ratios cannot be, or are unlikely to be, causative, that is, at 3′untranslated sites never used for translation termination, in intragenomic analyses, and across archaeal species (that possess only one RF1). We conclude that specifics of RF biology are unlikely to fully explain TGA/TAG relative usage. We discuss why the causal relationships for the evolution of synonymous stop codon usage might be different from those affecting synonymous sense codon usage, noting that transitions between TGA and TAG require two-point mutations one of which is likely to be deleterious.