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.

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.

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.

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.

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.

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.

Tetracoding increases with body temperature in Lepidosauria.

Codons expanded by a silent position (quadruplet or tetracodons) may solve the conundrum that at life's origins, the weak tricodon-anticodon interactions could not promote translation in the absence of complex ribosomes. Modern genomes have isolated tetracodons resulting from insertion mutations. Some bioinformatic analyses suggest that tetracoding stretches overlap with regular mitochondrial protein coding genes. These tetragenes are probably decoded by (antisense) tRNAs with expanded anticodons. They are GC-rich, which produce stronger basepairs than A:T interactions, suggesting expression at high temperatures. The hypothesis that tetracoding is an adaptation to high temperatures is tested here by comparing predicted mitochondrial tetracoding in Lepidosauria (lizards, amphisbaenia, and Sphenodon), in relation to body temperature, expecting more tetracoding in species with high body temperature. The association between tRNAs with expanded anticodons and tetracoding previously described for mammals and Drosophila is confirmed for Lepidosauria. Independent evidence indicates that tetracoding increases with body temperature, supporting the hypothesis that tetracoding is an adaptation for efficient translation when conditions (temperature) make triplet codon-anticodons too unstable to allow efficient protein elongation.

Translational control of cell division by Elongator.

Elongator is required for the synthesis of the mcm(5)s(2) modification found on tRNAs recognizing AA-ending codons. In order to obtain a global picture of the role of Elongator in translation, we used reverse protein arrays to screen the fission yeast proteome for translation defects. Unexpectedly, this revealed that Elongator inactivation mainly affected three specific functional groups including proteins implicated in cell division. The absence of Elongator results in a delay in mitosis onset and cytokinesis defects. We demonstrate that the kinase Cdr2, which is a central regulator of mitosis and cytokinesis, is under translational control by Elongator due to the Lysine codon usage bias of the cdr2 coding sequence. These findings uncover a mechanism by which the codon usage, coupled to tRNA modifications, fundamentally contributes to gene expression and cellular functions.

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.

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.

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.

Clonorchis sinensis: codon usage in nuclear genes.

Codon usage in Clonorchis sinensis was analyzed using 12,515 codons from 38 coding sequences. Total GC content was 49.83%, and GC1, GC2 and GC3 contents were 56.32%, 43.15% and 50.00%, respectively. The effective number of codons converged at 51-53 codons. When plotted against total GC content or GC3, codon usage was distributed in relation to GC3 biases. Relative synonymous codon usage for each codon revealed a single major trend, which was highly correlated with GC content at the third position when codons began with A or U at the first two positions. In codons beginning with G or C base at the first two positions, the G or C base rarely occurred at the third position. These results suggest that codon usage is shaped by a bias towards G or C at the third base, and that this is affected by the first and second bases.

Genome sequence, full-length infectious cDNA clone, and mapping of viral double-stranded RNA accumulation determinant of hypovirus CHV1-EP721.

Cryphonectria parasitica strain EP721 is infected with a strain of hypovirus CHV1, CHV1-EP721, and exhibits typical hypovirulence-associated traits such as reduced pigmentation and reduced asexual sporulation. However, the accumulation of the viral double-stranded RNA (dsRNA) in this hypovirus-infected C. parasitica strain is atypically low. We now report the complete nucleotide sequence and construction of a full-length infectious cDNA clone for hypovirus CHV1-EP721. The genome sequence of CHV1-EP721 was determined to be 12,724 bp in length and to share extensive homology with two other hypovirus strains, CHV1-Euro7 and CHV1-EP713, with an average of 99% and 90% identities at the nucleotide level and 99% and 92% identities at the amino acid level, respectively. CHV1-EP721 was successfully introduced into virus-free fungal host strain EP721(-v) by transfection with transcripts derived from a full-length viral cDNA. The transfected strain had a phenotype indistinguishable from that of EP721, and the accumulation of CHV1-EP721 dsRNA in the transfectant was lower than those transfected by CHV1-Euro7 and CHV1-EP713 transcripts. Through the construction of chimeric viruses by domain swapping using infectious cDNA clones of CHV1-EP721, CHV1-EP713, and CHV1-Euro7 hypoviruses, the determinant for the low level of viral dsRNA accumulation in CHV1-EP721 was mapped to the second of two CHV1-EP721 open reading frames (ORFs), ORF B. Further refined swapping of domains within ORF B identified a 2.5-kb coding region between p48 and the polymerase domain of CHV1-EP721 as being responsible for the low viral dsRNA accumulation. Evidence is also provided that low rates of hypovirus transmission through conidial spores correlates with low viral dsRNA accumulation.

Analysis of codon usage patterns and predicted highly expressed genes for six phytopathogenic Xanthomonas genomes shows a high degree of conservation.

Members of the genus Xanthomonas are significant phytopathogens, which cause diseases in several economically important crops including rice, canola, tomato, citrus, etc. We have analyzed the genomes of six recently sequenced Xanthomonas strains for their synonymous codon usage patterns for all of protein coding genes and specific genes associated with pathogenesis, and determined the predicted highly expressed (PHX) genes by the use of the codon adaptation index (CAI). Our results show considerable heterogeneity among the genes of these moderately G+C rich genomes. Most of the genes were moderate to highly biased in their codon usage. However, unlike ribosomal protein genes, which were governed by translational selection, those genes associated with pathogenesis (GAP) were affected by mutational pressure and were predicted to have moderate to low expression levels. Only two out of 339 GAP genes were in the PHX category. PHX genes present in clusters of orthologous groups of proteins (COGs) were identified. Genes in the plasmids present in two strains showed moderate to low expression level and only a couple of genes featured in the PHX list. Common genes present in the top-20 PHX gene-list were identified and their possible functions are discussed. Correspondence analysis showed that genes are highly confined to a core in the plot.

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.

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.

Studies on codon usage in Thermoplasma acidophilum and its possible implications on the occurrences of lateral gene transfer.

Codon usage studies have been carried out on the coding sequences of Thermoplasma acidophilum, which is an archaeon and grows at very low pH and high temperature. Overall codon usage data analysis indicates that all the four bases are almost equifrequent at the third position of codons, which is expected (since genomic GC % of this genome is about 46%). However, multivariate statistical analysis indicates that there are two major trends in the codon usage variation among the genes in this organism. In the first major trend it is observed that genes having G and C ending codons are clustered at one end while, A and T ending ones are clustered at the other end. We have also found a significant positive correlation between the expressivities of genes and GC contents at the synonymous third codon positions. In the second major trend, it is seen that the genes are clustered into three distinct parts. A comparative analyses of codon usage data of T. acidophilum and Sulfolobus solfataricus reveals that one of the three clusters of genes of T. acidophilum is very similar to a considerable number of S. solfataricus genes, suggesting possible occurrences of lateral gene transfer between these two microorganisms as reported by earlier workers.