Augmented genetic decoding: global, local and temporal alterations of decoding processes and codon meaning.

The non-universality of the genetic code is now widely appreciated. Codes differ between organisms, and certain genes are known to alter the decoding rules in a site-specific manner. Recently discovered examples of decoding plasticity are particularly spectacular. These examples include organisms and organelles with disruptions of triplet continuity during the translation of many genes, viruses that alter the entire genetic code of their hosts and organisms that adjust their genetic code in response to changing environments. In this Review, we outline various modes of alternative genetic decoding and expand existing terminology to accommodate recently discovered manifestations of this seemingly sophisticated phenomenon.

Triplet-Based Codon Organization Optimizes the Impact of Synonymous Mutation on Nucleic Acid Molecular Dynamics.

Since the elucidation of the genetic code almost 50 years ago, many nonrandom aspects of its codon organization remain only partly resolved. Here, we investigate the recent hypothesis of 'dual-use' codons which proposes that in addition to allowing adjustment of codon optimization to tRNA abundance, the degeneracy in the triplet-based genetic code also multiplexes information regarding DNA's helical shape and protein-binding dynamics while avoiding interference with other protein-level characteristics determined by amino acid properties. How such structural optimization of the code within eukaryotic chromatin could have arisen from an RNA world is a mystery, but would imply some preadaptation in an RNA context. We analyzed synonymous (protein-silent) and nonsynonymous (protein-altering) mutational impacts on molecular dynamics in 13823 identically degenerate alternative codon reorganizations, defined by codon transitions in 7680 GPU-accelerated molecular dynamic simulations of implicitly and explicitly solvated double-stranded aRNA and bDNA structures. When compared to all possible alternative codon assignments, the standard genetic code minimized the impact of synonymous mutations on the random atomic fluctuations and correlations of carbon backbone vector trajectories while facilitating the specific movements that contribute to DNA polymer flexibility. This trend was notably stronger in the context of RNA supporting the idea that dual-use codon optimization and informational multiplexing in DNA resulted from the preadaptation of the RNA duplex to resist changes to thermostability. The nonrandom and divergent molecular dynamics of synonymous mutations also imply that the triplet-based code may have resulted from adaptive functional expansion enabling a primordial doublet code to multiplex gene regulatory information via the shape and charge of the minor groove.

Red fluorescent protein (DsRFP) optimization for Entamoeba histolytica expression.

Entamoeba histolytica genetic organization and genome structure is complex and under intense research. The genome is fully sequenced, and several tools have been developed for the molecular study of this organism. Nevertheless, good protein tracking tags that are easy to measure and image, like the fluorescent proteins are lacking. In this report, we codon-optimized the red fluorescent protein from the coral Discosoma striata (DsRFP) for its use in E. histolytica and demonstrated functionality in vivo. We envision that this protein can be widely used for the development of transcriptional reporter systems and protein-tagging applications.

A detailed analysis of codon usage patterns and influencing factors in Zika virus.

Recent outbreaks of Zika virus (ZIKV) in Africa, Latin America, Europe, and Southeast Asia have resulted in serious health concerns. To understand more about evolution and transmission of ZIKV, detailed codon usage analysis was performed for all available strains. A high effective number of codons (ENC) value indicated the presence of low codon usage bias in ZIKV. The effect of mutational pressure on codon usage bias was confirmed by significant correlations between nucleotide compositions at third codon positions and ENCs. Correlation analysis between Gravy values, Aroma values and nucleotide compositions at third codon positions also indicated some influence of natural selection. However, the low codon adaptation index (CAI) value of ZIKV with reference to human and mosquito indicated poor adaptation of ZIKV codon usage towards its hosts, signifying that natural selection has a weaker influence than mutational pressure. Additionally, relative dinucleotide frequencies, geographical distribution, and evolutionary processes also influenced the codon usage pattern to some extent.

Different patterns of codon usage in the overlapping polymerase and surface genes of hepatitis B virus suggest a de novo origin by modular evolution.

The polymerase (P) and surface (S) genes of hepatitis B virus (HBV) show the longest gene overlap in animal viruses. Gene overlaps originate by the overprinting of a novel frame onto an ancestral pre-existing frame. Identifying which frame is ancestral and which frame is de novo (the genealogy of the overlap) is an appealing topic. However, the P/S overlap of HBV is an intriguing paradox, because both genes are indispensable for virus survival. Thus, the hypothesis of a primordial virus without the surface protein or without the polymerase makes no biological sense. With the aim to determine the genealogy of the overlap, the codon usage of the overlapping frames P and S was compared to that of the non-overlapping region. It was found that the overlap of human HBV had two patterns of codon usage. One was localized in the 59 one-third of the overlap and the other in the 39 two-thirds. By extending the analysis to non-human HBVs, it was found that this feature occurred in all hepadnaviruses. Under the assumption that the ancestral frame has a codon usage significantly closer to that of the non-overlapping region than the de novo frame, the ancestral frames in the 59 and 39 region of the overlap could be predicted. They were, respectively, frame S and frame P. These results suggest that the spacer domain of the polymerase and the S domain of the surface protein originated de novo by overprinting. They support a modular evolution hypothesis for the origin of the overlap.

A genetic scale of reading frame coding.

The reading frame coding (RFC) of codes (sets) of trinucleotides is a genetic concept which has been largely ignored during the last 50 years. A first objective is the definition of a new and simple statistical parameter PrRFC for analysing the probability (efficiency) of reading frame coding (RFC) of any trinucleotide code. A second objective is to reveal different classes and subclasses of trinucleotide codes involved in reading frame coding: the circular codes of 20 trinucleotides and the bijective genetic codes of 20 trinucleotides coding the 20 amino acids. This approach allows us to propose a genetic scale of reading frame coding which ranges from 1/3 with the random codes (RFC probability identical in the three frames) to 1 with the comma-free circular codes (RFC probability maximal in the reading frame and null in the two shifted frames). This genetic scale shows, in particular, the reading frame coding probabilities of the 12,964,440 circular codes (PrRFC=83.2% in average), the 216 C(3) self-complementary circular codes (PrRFC=84.1% in average) including the code X identified in eukaryotic and prokaryotic genes (PrRFC=81.3%) and the 339,738,624 bijective genetic codes (PrRFC=61.5% in average) including the 52 codes without permuted trinucleotides (PrRFC=66.0% in average). Otherwise, the reading frame coding probabilities of each trinucleotide code coding an amino acid with the universal genetic code are also determined. The four amino acids Gly, Lys, Phe and Pro are coded by codes (not circular) with RFC probabilities equal to 2/3, 1/2, 1/2 and 2/3, respectively. The amino acid Leu is coded by a circular code (not comma-free) with a RFC probability equal to 18/19. The 15 other amino acids are coded by comma-free circular codes, i.e. with RFC probabilities equal to 1. The identification of coding properties in some classes of trinucleotide codes studied here may bring new insights in the origin and evolution of the genetic code.

Gene functionality's influence on the second codon: A large-scale survey of second codon composition in three domains.

The second codon of a transcript, besides encoding for an amino acid, is now known to also have multiple molecular functions and is involved in translation efficiency and protein turn-over and maturation processing. These multiple purposes therefore make the selection constraints on this codon's composition more complex. To examine the biological significance of various permutations of the second codon, we conducted a systematic survey of second codon composition from 442 selected genomes across three domains. The amino acid bias of the second codon is associated with specific protein functions. The most common amino acids (S, A, K and T) are significantly avoided in Cell Envelope-related genes but preferred in Translation or Energy Metabolism-related genes, suggesting that the function of a gene product is a significant factor influencing the composition of the second codon.

Plasmodium: mammalian codon optimization of malaria plasmid DNA vaccines enhances antibody responses but not T cell responses nor protective immunity.

We have evaluated the effect of mammalian codon optimization on the immunogenicity and protective efficacy of plasmid DNA vaccines encoding pre-erythrocytic stage Plasmodium falciparum and Plasmodium yoelii antigens in mice. Codon optimization significantly enhanced in vitro expression and in vivo antibody responses for P. falciparum circumsporozoite protein (PfCSP) and P. yoelii hepatocyte erythrocyte protein 17 kDa (PyHEP17) but not for P. yoelii circumsporozoite protein (PyCSP). Unexpectedly, more robust CD4+ and CD8+ T cell responses as measured by IFN-gamma ELIspot, lymphoproliferation, and cytotoxic T lymphocyte assays were noted with native as compared with codon optimization constructs. Codon optimization also failed to enhance CD8+ T cell dependent protection against P. yoelii sporozoite challenge as measured by liver-stage parasite burden. These data demonstrate that the effect of mammalian codon optimization is antigen-dependent and may not be beneficial for vaccines designed to induce T cell dependent protective immunity in this malaria model.

Codon optimization can improve expression of human genes in Escherichia coli: A multi-gene study.

The efficiency of heterologous protein production in Escherichia coli (E. coli) can be diminished by biased codon usage. Approaches normally used to overcome this problem include targeted mutagenesis to remove rare codons or the addition of rare codon tRNAs in specific cell lines. Recently, improvements in technology have enabled cost-effective production of synthetic genes, making this a feasible alternative. To explore this option, the expression patterns in E. coli of 30 human short-chain dehydrogenase/reductase genes (SDRs) were analyzed in three independent experiments, comparing the native and synthetic (codon-optimized) versions of each gene. The constructs were prepared in a pET-derived vector that appends an N-terminal polyhistidine tag to the protein; expression was induced using IPTG and soluble proteins were isolated by Ni-NTA metal-affinity chromatography. Expression of the native and synthetic gene constructs was compared in two isogenic bacterial strains, one of which contained a plasmid (pRARE2) that carries seven tRNAs recognizing rare codons. Although we found some degree of variability between experiments, in normal E. coli synthetic genes could be expressed and purified more readily than the native version. In only one case was native gene expression better. Importantly, in most but not all cases, expression of the native genes in combination with rare codon tRNAs mimicked the behavior of the synthetic genes in the native strain. The trend is that heterologous expression of some proteins in bacteria can be improved by altering codon preference, but that this effect can be generally recapitulated by introducing rare codon tRNAs into the host cell.

The use of the rare UUA codon to define "expression space" for genes involved in secondary metabolism, development and environmental adaptation in streptomyces.

In Streptomyces coelicolor, bldA encodes the only tRNA for a rare leucine codon, UUA. This tRNA is unnecessary for growth, but is required for some aspects of secondary metabolism and morphological development, as revealed by the phenotypes of bldA mutants in diverse streptomycetes. This article is a comprehensive review of out understanding of this unusual situation. Based on information from four sequenced genomes it now appears that, typically, about 2 approximately 3% of genes in any one streptomycete contain a TTA codon, most having been acquired through species-specific horizontal gene transfer. Among the few widely conserved TTA-containing genes, mutations in just one, the pleiotropic regulatory gene adpA, give an obvious phenotype: such mutants are defective in aerial growth and sporulation, but vary in the extent of their impairment in secondary metabolism in different streptomycetes. The TTA codon in adpA is largely responsible for the morphological phenotype of a bldA mutant of S. coelicolor. AdpA-dependent targets include several genes involved in the integrated action of extracellular proteases that, at least in some species, are involved in the conversion of primary biomass into spores. The effects of bldA mutations on secondary metabolism are mostly attributable to the presence of TTA codons in pathway-specific genes, particularly in transcriptional activator genes. This is not confined to S. coelicolor-it is true for about half of all known antibiotic biosynthetic gene sets from streptomycetes. Combined microarray and proteomic analysis of liquid (and therefore non-sporulating) S. coelicolor bldA mutant cultures revealed effects of the mutation during rapid growth, during transition phase, and in stationary phase. Some of these effects may be secondary consequences of changes in the pattern of ppGpp accumulation. It is argued that the preferential accumulation of the bldA tRNA under conditions in which growth is significantly constrained has evolved to favour the expression of genes that confer adaptive benefits in intermittently encountered sub-optimal environments. The evolution of this system may have been a secondary consequence of the selective pressure exerted by bacteriophage attack. Some biotechnological implications of bldA phenomenology are considered.

Analysis of codon usage pattern of infectious laryngotracheitis virus immunogenic glycoproteins and its biological implications.

Infectious laryngotracheitis virus (ILTV) is a highly contagious acute respiratory poultry pathogen. Modified live ILTV vaccines are the only control against ILT infections. Reversions and establishment of latent infections are the major concerns imparting the need to develop safer vaccines against ILTV infection. ILTV glycoprotein B and D (gB and gD) are major protective immunogens. The factors shaping synonymous codon usage bias and nucleotide composition in ILTV glycoprotein genes have not yet been reported. In the present study, we have analyzed the synonymous codon usage indices of ILTV gB and gD genes. Variation in the codon usage was seen in both the glycoproteins majorly by mutational pressure. The pattern was determined using the correspondence analysis, effective number of codon (Nc), GC3 plot and correlation analyses among different indices. The study is a comprehensive analysis of the codon usage patterns of ILTV glycoprotein genes. This will be helpful in understanding the codon usage bias of ILTV and related DNA viruses which could further explore its biology.

The repertoire of transfer RNA genes is tuned to codon usage bias in the genomes of Phytophthora sojae and Phytophthora ramorum.

In all, 238 and 155 transfer (t)RNA genes were predicted from the genomes of Phytophthora sojae and P. ramorum, respectively. After omitting pseudogenes and undetermined types of tRNA genes, there remained 208 P. sojae tRNA genes and 140 P. ramorum tRNA genes. There were 45 types of tRNA genes, with distinct anticodons, in each species. Fourteen common anticodon types of tRNAs are missing altogether from the genome in the two species; however, these appear to be compensated by wobbling of other tRNA anticodons in a manner which is tied to the codon bias in Phytophthora genes. The most abundant tRNA class was arginine in both P. sojae and P. ramorum. A codon usage table was generated for these two organisms from a total of 9,803,525 codons in P. sojae and 7,496,598 codons in P. ramorum. The most abundant codon type detected from the codon usage tables was GAG (encoding glutamic acid), whereas the most numerous tRNA gene had a methionine anticodon (CAT). The correlation between the frequencies of tRNA genes and the codon frequencies in protein-coding genes was very low (0.12 in P. sojae and 0.19 in P. ramorum); however, the correlation between amino acid tRNA gene frequency and the corresponding amino acid codon frequency in P. sojae and P. ramorum was substantially higher (0.53 in P. sojae and 0.77 in P. ramorum). The codon usage frequencies of P. sojae and P ramorum were very strongly correlated (0.99), as were tRNA gene frequencies (0.77). Approximately 60% of orthologous tRNA gene pairs in P sojae and P. ramorum are located in regions that have conserved synteny in the two species.

Codons 12 and 13 of H-ras protooncogene interrupt the progression of DNA synthesis catalyzed by DNA polymerase alpha.

Mutagenesis of protooncogenes has been postulated to contribute to the initiation and progression of human cancer. Activating mutations in the H-ras gene are predominantly single-base substitutions and are most frequently identified at codons 12, 13, and 61. We have analyzed the effects of DNA sequence context at specific codons that are hot spots for ras mutation with respect to abnormalities in copying by purified DNA polymerase alpha, a major eucaryotic replication enzyme. Exon 1 of H-ras gene was inserted into M13 mp19, single-stranded DNA constructs were isolated, and the progression of synthesis by polymerase alpha was measured. Strong termination sites were found in codons 12 and 13. Pausing at these codons is abolished when the template is mutated at the middle base of codon 12, the same alteration that converts H-ras into an activated oncogene. Resistance of codon 12 in double-stranded constructs to digestion with restriction enzymes and computer investigation of the ras sequence suggest that these termination sites are in a region of secondary structure. The frequency of sequence alterations within DNA chains that have been extended past codons 12 and 13 was found to be < 0.01. We consider a variety of mechanisms by which the potential secondary structure involving codons 12 and 13 may contribute to the pausing of DNA polymerase alpha and to the generation of clustered mutations at this site.

The 273rd codon mutants of p53 show growth modulation activities not correlated with p53-specific transactivation activity.

Cancer-related mutations of the p53 tumor suppressor gene are clustered in the four so-called 'hot spots', codons 175, 248, 273 and 281/282. By using recombination PCR in vitro mutagenesis, we introduced point mutations into the codon 273 of wild-type (wt) p53 (pC53-SN3) from Arg to His (pC53-273H [273H]), Asp (273D), Pro (273P), Lys (273K), Leu (273L) or Thr (273T), and compared their biological and biochemical activities with wt p53 and cancer-derived 175H, 248W and 273H/309S. Among them, 273H/309S, 273H and 273D as well as wt p53 transactivated the chloramphenicol acetyltransferase (CAT) gene placed downstream of the p53 binding consensus, while none of the other mutants including 273L did. Transcriptions from human c-fos and rat PCNA promoters were suppressed by wt p53 and 273D, while they were enhanced variously by all other mutants in Saos-2 and/or NIH3T3 cells. On the other hand, growth of human squamous carcinoma cell lines measured by the plating efficiency of G418-resistant colonies was enhanced by transfection of 175H, 248W, 273H/309S and 273P, while suppressed by not only wt p53, 273D and 273H but also 273L. Thus, 273H/309S enhanced cell growth in spite of its p53-specific transactivation activity, while 273L suppressed cell growth in spite of its complete loss of the p53-specific transactivation. We concluded that the sequence-specific transactivation of p53 is not always correlated with its growth inhibitory activity.

Selective inhibition of the polypeptide chain elongation in eukaryotic cells.

The effect of Cephalotaxus alkaloids--homoharringtonine and cephalotaxine--on translation in a cell-free system from rabbit reticulocytes and on phenylalanine polymerisation by human ribosomes was studied. The effect of the alkaloids on the nonenzymatic and the eEF-1-dependent Phe-tRNA(Phe) binding to poly(U)-programmed 80S ribosomes, diphenylalanine synthesis accompanying nonenzymatic Phe-tRNA(Phe) binding and acetylphenylalanyl-puromycin formation was examined. Homoharringtonine was shown to inhibit the formation of diphenylalanine and acetylphenylalanyl-puromycin catalysed by human and rat liver ribosomes, but was inactive as an inhibitor on the E. coli elongation system. Neither nonenzymatic nor enzymatic Phe-tRNA(Phe) binding was noticeably affected by the alkaloid. It has been proposed that the site of homoharringtonine binding to 80S ribosomes should overlap or coincide with the acceptor site of the ribosomal peptidyl transferase centre. The association constant of homoharringtonine for 80S human ribosomes was estimated to be (2.57 +/- 0.33).10(7) M-1 in the presence of puromycin. Cephalotaxine did not exert a significant influence on the polypeptide chain elongation.

Promoter cloning and characterisation of the transcriptional regulation of the human thyrostimulin A2 subunit.

The novel heterodimeric glycoprotein hormone thyrostimulin consists of two unique subunits, A2 and B5. To understand its yet unknown transcriptional regulation, we characterised the 3.1-kb immediate 5'-flanking region of the human A2 gene localised on chromosome 11q13. In transient transfection assays this sequence exhibited promoter activity, which could be confined to nucleotides -506 to -347 relative to the ATG start codon. Interestingly, this minimal promoter appeared to be non-tissue-specific. Deletional, mutational and gel shift analyses revealed regulatory elements that are essential for the regulation of the A2 gene expression. Another noteworthy feature of this gene is the presence of silencer elements upstream and downstream of the promoter. To surmise, our results provide an initial step toward a detailed analysis of the underlying molecular mechanisms of the human thyrostimulin gene expression.

Cloning and sequence analysis of a polygalacturonase-encoding gene from the phytopathogenic fungus Sclerotinia sclerotiorum.

The phytopathogenic fungus Sclerotinia sclerotiorum produces a number of extra-cellular pectin-degrading enzymes. We have cloned and determined the complete sequence of a gene (pg1) encoding an endopolygalacturonase (PG1). The coding region consists of a non-interrupted 1143-bp open reading frame. S. sclerotiorum pg1 was compared to other fungal PG-encoding genes. Basic transcription control sequences were identified in the 5' non-coding region. The deduced amino acid (aa) sequence (380 aa) of the enzyme is compared to seven fungal PG sequences and shows a high level of identity (41.5 to 59.8%). Predicted secondary structures were compared, revealing a similar protein organization most probably in antiparallel beta sheets. Hybridization analysis using a pg1 0.65-kb BamHI fragment as a probe allowed the identification of seven different recombinant phages from a genomic library. Analysis of the hybridizing restriction fragments suggests that PG-encoding genes are organized as a family.

A functional significance for codon third bases.

Most amino acids are specified by more than one trinucleotide codon. Here we show that amino acids of differing functional importance may be distinguished by the pattern of synonymous codon usage. GC-rich genes tend to be of a greater transcriptional (p<0.01) and mitogenic (p<0.0001) significance than AT-rich genes, consistent with GC-->AT mutational drift in methylated genomic regions. Third-base GC retention also identifies critical amino acids within individual proteins, as indicated by non-random patterns of codon variation between gene homologs and also by differential sequelae of site-directed mutagenesis. Sequence analysis of human receptor tyrosine kinase genes confirms that functionally important transmembrane hydrophobic amino acids are specified by codons containing GC third bases more often than are transmembrane neutral amino acids (chi(2)=134.2). Amino acids encoded by GC third bases thus appear more tightly linked to cell function and survival than are those encoded by AT third bases.

The psbC start codon in Synechocystis sp. PCC 6803.

The translation start codon for psbC, the gene encoding CP43, a chlorophyll-binding protein of photosystem II, has been identified for the cyanobacterium Synechosystis sp. PCC 6803 using site-directed mutagenesis. An AUG codon, about 50 bases upstream from the end of psbD-I had previously been assumed to be the translation start site of psbC. However, the fact that the AUG codon is not present in psbC from several other organisms, whereas a GUG codon 14 bases upstream from the end of psbD-I is strictly conserved suggests that CP43 translation starts at the latter codon. Mutation of GUG, but not of AUG, led to a loss of CP43 and photoautotrophic growth, indicating that the GUG codon is the sole initiation site for translation of the CP43 protein in Synechocystis sp. PCC 6803.

Mutagenesis using trinucleotide beta-cyanoethyl phosphoramidites.

There is no easy way to selectively introduce mixtures of codon triplets into mutagenesis libraries. Solid-phase-supported DNA synthesis using successive coupling of mixtures of mononucleotides can be made to supply 32 codons, which gives redundancies in coding for 20 natural amino acids, as well as an often unwanted stop codon. Resin-splitting methods have been described, but the representation of all permutations is limited by mechanical factors for a large library, and the method is experimentally cumbersome. To demonstrate a third, improved method, the 3'-cyanoethyl phosphoramidite codon triplets dATA, dCTT, dATC, dATG and dAGC were made by solution-phase methods, with protecting groups fully compatible with modern automated phosphoramidite DNA synthesis chemistry. The reagents were then used to synthesize a 54-mer DNA fragment, wherein 15 internal base pairs were randomized by coupling a mixture of the five codons five times. The fragment was amplified as a cDNA pool, which was subcloned into a phagemid vector, and 16 randomly selected recombinants from this mini-library were sequenced. These clones showed random incorporation of the proper transcribed codon sequences at the correct location. Other functional tests involving the trinucleotide phosphoramidites showed modest (ca. 70%) coupling efficiencies and structural integrity of the DNA produced.