Codon usage <Emphasis Type="BoldItalic">vis-a-vis</Emphasis> start and stop codon context analysis of three dicot species

To understand the variation in genomic composition and its effect on codon usage, we performed the comparative analysis of codon usage and nucleotide usage in the genes of three dicots, Glycine max, Arabidopsis thaliana and Medicago truncatula. The dicot genes were found to be A/T rich and have predominantly A-ending and/or T-ending codons. GC3s directly mimic the usage pattern of global GC content. Relative synonymous codon usage analysis suggests that the high usage frequency of A/T over G/C mononucleotide containing codons in AT-rich dicot genome is due to compositional constraint as a factor of codon usage bias. Odds ratio analysis identified the dinucleotides TpG, TpC, GpA, CpA and CpT as over-represented, where, CpG and TpA as under-represented dinucleotides. The results of (NcExp?NcObs)/NcExp plot suggests that selection pressure other than mutation played a significant role in influencing the pattern of codon usage in these dicots. PR2 analysis revealed the significant role of selection pressure on codon usage. Analysis of varience on codon usage at start and stop site showed variation in codon selection in these sites. This study provides evidence that the dicot genes were subjected to compositional selection pressure.     

Virus-Host Coevolution: Common Patterns of Nucleotide Motif Usage in Flaviviridae and Their Hosts

Virus-host biological interaction is a continuous coevolutionary process involving both host immune system and viral escape mechanisms. Flaviviridae family is composed of fast evolving RNA viruses that infects vertebrate (mammals and birds) and/or invertebrate (ticks and mosquitoes) organisms. These host groups are very distinct life forms separated by a long evolutionary time, so lineage-specific anti-viral mechanisms are likely to have evolved. Flaviviridae viruses which infect a single host lineage would be subjected to specific host-induced pressures and, therefore, selected by them. In this work we compare the genomic evolutionary patterns of Flaviviridae viruses and their hosts in an attempt to uncover coevolutionary processes inducing common features in such disparate groups. Especially, we have analyzed dinucleotide and codon usage patterns in the coding regions of vertebrate and invertebrate organisms as well as in Flaviviridae viruses which specifically infect one or both host types. The two host groups possess very distinctive dinucleotide and codon usage patterns. A pronounced CpG under-representation was found in the vertebrate group, possibly induced by the methylation-deamination process, as well as a prominent TpA decrease. The invertebrate group displayed only a TpA frequency reduction bias. Flaviviridae viruses mimicked host nucleotide motif usage in a host-specific manner. Vertebrate-infecting viruses possessed under-representation of CpG and TpA, and insect-only viruses displayed only a TpA under-representation bias. Single-host Flaviviridae members which persistently infect mammals or insect hosts (Hepacivirus and insect-only Flavivirus, respectively) were found to posses a codon usage profile more similar to that of their hosts than to related Flaviviridae. We demonstrated that vertebrates and mosquitoes genomes are under very distinct lineage-specific constraints, and Flaviviridae viruses which specifically infect these lineages appear to be subject to the same evolutionary pressures that shaped their host coding regions, evidencing the lineage-specific coevolutionary processes between the viral and host groups.     

Regularities of context-dependent codon bias in eukaryotic genes

Nucleotides surrounding a codon influence the choice of this particular codon from among the group of possible synonymous codons. The strongest influence on codon usage arises from the nucleotide immediately following the codon and is known as the N₁ context. We studied the relative abundance of codons with N₁ contexts in genes from four eukaryotes for which the entire genomes have been sequenced: Homo sapiens, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana. For all the studied organisms it was found that 90% of the codons have a statistically significant N₁ context-dependent codon bias. The relative abundance of each codon with an N₁ context was compared with the relative abundance of the same 4mer oligonucleotide in the whole genome. This comparison showed that in about half of all cases the context-dependent codon bias could not be explained by the sequence composition of the genome. Ranking statistics were applied to compare context-dependent codon biases for codons from different synonymous groups. We found regularities in N₁ context-dependent codon bias with respect to the codon nucleotide composition. Codons with the same nucleotides in the second and third positions and the same N₁ context have a statistically significant correlation of their relative abundances.     

Genome landscapes and bacteriophage codon usage

Across all kingdoms of biological life, protein-coding genes exhibit unequal usage of synonymous codons. Although alternative theories abound, translational selection has been accepted as an important mechanism that shapes the patterns of codon usage in prokaryotes and simple eukaryotes. Here we analyze patterns of codon usage across 74 diverse bacteriophages that infect E. coli, P. aeruginosa, and L. lactis as their primary host. We use the concept of a “genome landscape,” which helps reveal non-trivial, long-range patterns in codon usage across a genome. We develop a series of randomization tests that allow us to interrogate the significance of one aspect of codon usage, such as GC content, while controlling for another aspect, such as adaptation to host-preferred codons. We find that 33 phage genomes exhibit highly non-random patterns in their GC3-content, use of host-preferred codons, or both. We show that the head and tail proteins of these phages exhibit significant bias towards host-preferred codons, relative to the non-structural phage proteins. Our results support the hypothesis of translational selection on viral genes for host-preferred codons, over a broad range of bacteriophages.     

The Distribution of Synonymous Codon Choice in the Translation Initiation Region of Dengue Virus

Dengue is the most common arthropod-borne viral (Arboviral) illness in humans. The genetic features concerning the codon usage of dengue virus (DENV) were analyzed by the relative synonymous codon usage, the effective number of codons and the codon adaptation index. The evolutionary distance between DENV and the natural hosts (Homo sapiens, Pan troglodytes, Aedes albopictus and Aedes aegypti) was estimated by a novel formula. Finally, the synonymous codon usage preference for the translation initiation region of this virus was also analyzed. The result indicates that the general trend of the 59 synonymous codon usage of the four genotypes of DENV are similar to each other, and this pattern has no link with the geographic distribution of the virus. The effect of codon usage pattern of Aedes albopictus and Aedes aegypti on the formation of codon usage of DENV is stronger than that of the two primates. Turning to the codon usage preference of the translation initiation region of this virus, some codons pairing to low tRNA copy numbers in the two primates have a stronger tendency to exist in the translation initiation region than those in the open reading frame of DENV. Although DENV, like other RNA viruses, has a high mutation to adapt its hosts, the regulatory features about the synonymous codon usage have been ‘branded’ on the translation initiation region of this virus in order to hijack the translational mechanisms of the hosts.     

A general model of codon bias due to GC mutational bias

Background

In spite of extensive research on the effect of mutation and selection on codon usage, a general model of codon usage bias due to mutational bias has been lacking. Because most amino acids allow synonymous GC content changing substitutions in the third codon position, the overall GC bias of a genome or genomic region is highly correlated with GC3, a measure of third position GC content. For individual amino acids as well, G/C ending codons usage generally increases with increasing GC bias and decreases with increasing AT bias. Arginine and leucine, amino acids that allow GC-changing synonymous substitutions in the first and third codon positions, have codons which may be expected to show different usage patterns.

Principal Findings

In analyzing codon usage bias in hundreds of prokaryotic and plant genomes and in human genes, we find that two G-ending codons, AGG (arginine) and TTG (leucine), unlike all other G/C-ending codons, show overall usage that decreases with increasing GC bias, contrary to the usual expectation that G/C-ending codon usage should increase with increasing genomic GC bias. Moreover, the usage of some codons appears nonlinear, even nonmonotone, as a function of GC bias. To explain these observations, we propose a continuous-time Markov chain model of GC-biased synonymous substitution. This model correctly predicts the qualitative usage patterns of all codons, including nonlinear codon usage in isoleucine, arginine and leucine. The model accounts for 72%, 64% and 52% of the observed variability of codon usage in prokaryotes, plants and human respectively. When codons are grouped based on common GC content, 87%, 80% and 68% of the variation in usage is explained for prokaryotes, plants and human respectively.

Conclusions

The model clarifies the sometimes-counterintuitive effects that GC mutational bias can have on codon usage, quantifies the influence of GC mutational bias and provides a natural null model relative to which other influences on codon bias may be measured.     

Insights into the riddles of codon usage patterns and codon context signatures in fungal genus Puccinia,a persistent threat to global agriculture

The fungal genus Puccinia, comprising of several menacing pathogens, has been a persistent peril to global agriculture. Genome sequencing of various members of Puccinia offers a scope to excavate their genomic riddles. The present study has been addressed at exploring the complex niceties of codon and amino acid usage patterns and subsequent elucidation of the determinants that drive such behavior. Multivariate statistical analysis revealed a complex interplay of natural selection for translation and compositional bias to be operational on the codon usage patterns. Gene expression level was observed to be the most competent factor governing codon usage behavior of the genus. In spite of subtle AT richness of the genus, potential highly expressed gene sets were found to preferentially employ GC rich optimal codons. Estimation of relative dinucleotide abundance revealed preference toward the employment of GpA, CpA, TpC, and TpG dinucleotides and restraint from using TpA dinucleotide among the members of the genus. Extensive codon context analysis revealed that codon pairs with GpA, CpA, TpC, and TpG dinucleotides were over-represented and codon pairs with TpA dinucleotide were extensively avoided at the codon–codon (cP3–cA1) junctions. Amino acid usage signatures of the genus were found to be influenced considerably by several imperative factors like aromatic and hydrophobic character of the encoded gene products, genomic compositional constraint, and gene expressivity. Detailed know-how of the potential highly expressed gene sets and associated optimal codons in the genus promise to be informative for the scientific community engaged in combating Puccinia pathogenesis.     

Comparative genome sequence analysis of Sulfolobus acidocaldarius and 9 other isolates of its genus for factors influencing codon and amino acid usage

In the present study, major constraints for codon and amino acid usage of Sulfolobus acidocaldarius, Sulfolobus solfataricus, Sulfolobus tokodali, Sulfolobus islandis and 6 other isolates from islandicus species of genus Sulfolobus were investigated. Correspondence analysis revealed high significant correlation between the major trend of synonymous codon usage and gene expression level, as assessed by the “Codon Adaptation Index” (CAI). There is a significant negative correlation between Nc (Effective number of codons) and CAI demonstrating role of codon bias as an important determinant of codon usage. The significant correlation between major trend of synonymous codon usage and GC3s (G + C at third synonymous position) indicated dominant role of mutational bias in codon usage pattern. The result was further supported from SCUO (synonymous codon usage order) analysis. The amino acid usage was found to be significantly influenced by aromaticity and hydrophobicity of proteins. However, translational selection which causes a preference for codons that are most rapidly translated by current tRNA with multiple copy numbers was not found to be highly dominating for all studied isolates. Notably, 26 codons that were found to be optimally used by genes of S. acidocaldarius at higher expression level and its comparative analysis with 9 other isolates may provide some useful clues for further in vivo genetic studies on this genus.     

Selection for minimization of translational frameshifting errors as a factor in the evolution of codon usage

In a wide range of genomes, it was observed that the usage of synonymous codons is biased toward specific codons and codon patterns. Factors that are implicated in the selection for codon usage include facilitation of fast and accurate translation. There are two types of translational errors: missense errors and processivity errors. There is considerable evidence in support of the hypothesis that codon usage is optimized to minimize missense errors. In contrast, little is known about the relationship between codon usage and frameshifting errors, an important form of processivity errors, which appear to occur at frequencies comparable to the frequencies of missense errors. Based on the recently proposed pause-and-slip model of frameshifting, we developed Frameshifting Robustness Score (FRS). We used this measure to test if the pattern of codon usage indicates optimization against frameshifting errors. We found that the FRS values of protein-coding sequences from four analyzed genomes (the bacteria Bacillus subtilis and Escherichia coli, and the yeasts Saccharomyces cerevisiae and Schizosaccharomyce pombe) were typically higher than expected by chance. Other properties of FRS patterns observed in B. subtilis, S. cerevisiae and S. pombe, such as the tendency of FRS to increase from the 5′- to 3′-end of protein-coding sequences, were also consistent with the hypothesis of optimization against frameshifting errors in translation. For E. coli, the results of different tests were less consistent, suggestive of a much weaker optimization, if any. Collectively, the results fit the concept of selection against mistranslation-induced protein misfolding being one of the factors shaping the evolution of both coding and non-coding sequences.     

The Effects of Codon Context on In Vivo Translation Speed

We developed a bacterial genetic system based on translation of the his operon leader peptide gene to determine the relative speed at which the ribosome reads single or multiple codons in vivo. Low frequency effects of so-called “silent” codon changes and codon neighbor (context) effects could be measured using this assay. An advantage of this system is that translation speed is unaffected by the primary sequence of the His leader peptide. We show that the apparent speed at which ribosomes translate synonymous codons can vary substantially even for synonymous codons read by the same tRNA species. Assaying translation through codon pairs for the 5′- and 3′- side positioning of the 64 codons relative to a specific codon revealed that the codon-pair orientation significantly affected in vivo translation speed. Codon pairs with rare arginine codons and successive proline codons were among the slowest codon pairs translated in vivo. This system allowed us to determine the effects of different factors on in vivo translation speed including Shine-Dalgarno sequence, rate of dipeptide bond formation, codon context, and charged tRNA levels.     

Impact of the six nucleotides downstream of the stop codon on translation termination

The efficiency of translation termination is influenced by local contexts surrounding stop codons. In Saccharomyces cerevisiae, upstream and downstream sequences act synergistically to influence the translation termination efficiency. By analysing derivatives of a leaky stop codon context, we initially demonstrated that at least six nucleotides after the stop codon are a key determinant of readthrough efficiency in S. cerevisiae. We then developed a combinatorial-based strategy to identify poor 3′ termination contexts. By screening a degenerate oligonucleotide library, we identified a consensus sequence –CA(A/G)N(U/C/G)A–, which promotes >5% readthrough efficiency when located downstream of a UAG stop codon. Potential base pairing between this stimulatory motif and regions close to helix 18 and 44 of the 18S rRNA provides a model for the effect of the 3′ stop codon context on translation termination.     

Genetic Inactivation of Poliovirus Infectivity by Increasing the Frequencies of CpG and UpA Dinucleotides within and across Synonymous Capsid Region Codons

Replicative fitness of poliovirus can be modulated systematically by replacement of preferred capsid region codons with synonymous unpreferred codons. To determine the key genetic contributors to fitness reduction, we introduced different sets of synonymous codons into the capsid coding region of an infectious clone derived from the type 2 prototype strain MEF-1. Replicative fitness in HeLa cells, measured by plaque areas and virus yields in single-step growth experiments, decreased sharply with increased frequencies of the dinucleotides CpG (suppressed in higher eukaryotes and most RNA viruses) and UpA (suppressed nearly universally). Replacement of MEF-1 capsid codons with the corresponding codons from another type 2 prototype strain (Lansing), a randomization of MEF-1 synonymous codons, increased the %G+C without increasing CpG, and reductions in the effective number of codons used had much smaller individual effects on fitness. Poliovirus fitness was reduced to the threshold of viability when CpG and UpA dinucleotides were saturated within and across synonymous codons of a capsid region interval representing only ∼9% of the total genome. Codon replacements were associated with moderate decreases in total virion production but large decreases in the specific infectivities of intact poliovirions and viral RNAs. Replication of codon replacement viruses, but not MEF-1, was temperature sensitive at 39.5°C. Synthesis and processing of viral intracellular proteins were largely unaltered in most codon replacement constructs. Replacement of natural codons with synonymous codons with increased frequencies of CpG and UpA dinucleotides may offer a general approach to the development of attenuated vaccines with well-defined antigenicities and very high genetic stabilities.Diversification of genomic sequences is constrained in all biological systems. At the level of primary sequences, the range of variability in coding regions is restricted by the codon usage bias (CUB), whereby a subset of synonymous codons are preferentially used in translation (24, 53, 69). The intensity of the CUB and the specific set of preferred codons vary widely across biological systems (39). Intertwined with the CUB is the suppression of the dinucleotides CpG and TpA (or UpA in RNA viruses) in the genomes of higher eukaryotes (4, 7, 26, 61) and many of their RNA viruses and small DNA viruses (28, 49). Variation in the primary sequences of RNA virus genomes is further constrained by requirements to maintain essential secondary and higher-order structures (42, 54, 68).We previously described the modulation of the replicative fitness of the Sabin type 2 oral poliovirus vaccine (OPV) strain (Sabin 2) by systematically changing the CUB in the capsid region, replacing the naturally occurring preferred codons with an unpreferred synonymous codon (isocodon) for each of nine amino acids (8). We called our approach “codon deoptimization” to contrast with the process of codon optimization, which is frequently used to maximize expression of foreign proteins in heterologous host systems (1, 27, 70). Apart from its potential application to development of improved poliovirus vaccines (8, 13, 38), experimental investigations of codon deoptimization directly test the relationships between replicative fitness, the extent of CUB, and the intensity of CpG and UpA suppression. As a model system for such studies, polioviruses offer several favorable properties, including (i) intrinsically high error rates for the poliovirus RNA-dependent RNA polymerase (2, 14, 16, 65), (ii) very high evolution rates (25), (iii) short generation times (8 to 10 h) and large progeny yields of prototype polioviruses, and (iv) well-developed reverse genetics (9).In this report, we extend our codon deoptimization strategy to the type 2 wild poliovirus prototype strain MEF-1. As before, we restricted our replacement of synonymous codons to the capsid coding region, which encodes two of the defining properties of polioviruses, namely, (i) the capacity to bind the CD155 poliovirus receptor (PVR) (23) and (ii) the poliovirus type-specific neutralizing antigenic sites (35). No changes were made to the flanking 5′-untranslated region and noncapsid region sequences, as they contain essential secondary structural elements (42, 54, 68) and are frequently exchanged out by recombination during circulation of poliovirus in human populations (20, 30, 32). MEF-1 was selected because of its high fitness level (hence, its use as the type 2 component of the inactivated poliovirus vaccine [IPV]) and because of its neurovirulence for humans (15), for nontransgenic mice (52), and for transgenic mice expressing the PVR (71). Type 2 polioviruses were selected first for study because the Sabin 2 OPV strain is most frequently associated with vaccine-associated paralytic poliomyelitis in contacts of OPV recipients (57, 59), with prolonged excretion of immunodeficiency-associated vaccine-derived polioviruses (VDPVs) (10, 31, 60), and with the emergence of circulating VDPVs in areas of low OPV coverage (10, 31).Consistent with our previous findings, the fitness of MEF-1 decreased in proportion to the total number of synonymous replacement codons. Fitness was reduced most efficiently by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous codons. Saturation of CpG and UpA in a small capsid interval (representing only ∼9% of the genome) reduced fitness to the threshold of viability, even though the MEF-1 amino acid sequence was unaltered. The most prominent biological effect of deoptimization of codon usage and the large-scale incorporation of CpG and UpA was a sharp reduction in virus specific infectivities. In contrast, translation and processing of viral proteins and yields of intact virus particles with native antigenicities were reduced only moderately by increased CpG and UpA frequencies. Codon deoptimization with concurrent increases in the frequencies of CpG and/or UpA dinucleotides in RNA virus genomes may provide a novel general approach to the rational design of improved attenuated vaccines with predictable and stable genetic properties.     

Optimizing heterologous expression in dictyostelium: importance of 5' codon adaptation

Expression of heterologous proteins in Dictyostelium discoideum presents unique research opportunities, such as the functional analysis of complex human glycoproteins after random mutagenesis. In one study, human chorionic gonadotropin (hCG) and human follicle stimulating hormone were expressed in Dictyostelium. During the course of these experiments, we also investigated the role of codon usage and of the DNA sequence upstream of the ATG start codon. The Dictyostelium genome has a higher AT content than the human, resulting in a different codon preference. The hCG-β gene contains three clusters with infrequently used codons that were changed to codons that are preferred by Dictyostelium. The results reported here show that optimizing the first 5–17 codons of the hCG gene contributes to 4- to 5-fold increased expression levels, but that further optimization has no significant effect. These observations suggest that optimal codon usage contributes to ribosome stabilization, but does not play an important role during the elongation phase of translation. Furthermore, adapting the 5′-sequence of the hCG gene to the Dictyostelium ‘Kozak’-like sequence increased expression levels ~1.5-fold. Thus, using both codon optimization and ‘Kozak’ adaptation, a 6- to 8-fold increase in expression levels could be obtained for hCG.     

The Highly Conserved Codon following the Slippery Sequence Supports ?1 Frameshift Efficiency at the HIV-1 Frameshift Site

HIV-1 utilises −1 programmed ribosomal frameshifting to translate structural and enzymatic domains in a defined proportion required for replication. A slippery sequence, U UUU UUA, and a stem-loop are well-defined RNA features modulating −1 frameshifting in HIV-1. The GGG glycine codon immediately following the slippery sequence (the ‘intercodon’) contributes structurally to the start of the stem-loop but has no defined role in current models of the frameshift mechanism, as slippage is inferred to occur before the intercodon has reached the ribosomal decoding site. This GGG codon is highly conserved in natural isolates of HIV. When the natural intercodon was replaced with a stop codon two different decoding molecules—eRF1 protein or a cognate suppressor tRNA—were able to access and decode the intercodon prior to −1 frameshifting. This implies significant slippage occurs when the intercodon is in the (perhaps distorted) ribosomal A site. We accommodate the influence of the intercodon in a model of frame maintenance versus frameshifting in HIV-1.     

The Selective Advantage of Synonymous Codon Usage Bias in Salmonella

The genetic code in mRNA is redundant, with 61 sense codons translated into 20 different amino acids. Individual amino acids are encoded by up to six different codons but within codon families some are used more frequently than others. This phenomenon is referred to as synonymous codon usage bias. The genomes of free-living unicellular organisms such as bacteria have an extreme codon usage bias and the degree of bias differs between genes within the same genome. The strong positive correlation between codon usage bias and gene expression levels in many microorganisms is attributed to selection for translational efficiency. However, this putative selective advantage has never been measured in bacteria and theoretical estimates vary widely. By systematically exchanging optimal codons for synonymous codons in the tuf genes we quantified the selective advantage of biased codon usage in highly expressed genes to be in the range 0.2–4.2 x 10⁻⁴ per codon per generation. These data quantify for the first time the potential for selection on synonymous codon choice to drive genome-wide sequence evolution in bacteria, and in particular to optimize the sequences of highly expressed genes. This quantification may have predictive applications in the design of synthetic genes and for heterologous gene expression in biotechnology.     

Synonymous codon usage in chloroplast genome of Coffea arabica

Synonymous codon usage of 53 protein coding genes in chloroplast genome of Coffea arabica was analyzed for the first time to find out the possible factors contributing codon bias. All preferred synonymous codons were found to use A/T ending codons as chloroplast genomes are rich in AT. No difference in preference for preferred codons was observed in any of the two strands, viz., leading and lagging strands. Complex correlations between total base compositions (A, T, G, C, GC) and silent base contents (A₃, T₃, G₃, C₃, GC₃) revealed that compositional constraints played crucial role in shaping the codon usage pattern of C. arabica chloroplast genome. ENC Vs GC₃ plot grouped majority of the analyzed genes on or just below the left side of the expected GC₃ curve indicating the influence of base compositional constraints in regulating codon usage. But some of the genes lie distantly below the continuous curve confirmed the influence of some other factors on the codon usage across those genes. Influence of compositional constraints was further confirmed by correspondence analysis as axis 1 and 3 had significant correlations with silent base contents. Correlation of ENC with axis 1, 4 and CAI with 1, 2 prognosticated the minor influence of selection in nature but exact separation of highly and lowly expressed genes could not be seen. From the present study, we concluded that mutational pressure combined with weak selection influenced the pattern of synonymous codon usage across the genes in the chloroplast genomes of C. arabica.     

On finding poorly translated codons based on their usage frequency

Long stretches of “rare” codons are known to severely inhibit the efficiency of translation. Understanding the distribution of such rare codons is of critical importance in improving the efficiency of heterologous gene expression systems. Accurate estimates of codon usage take the abundance of each protein into consideration. In this paper, we analyze the correlation between approximate measures of codon usage and the availability of tRNA at various growth rates in E coli. We show that the computationally derived estimates of tRNA isoacceptor concentration enable the finding of poorly translated codons.     

The Effect of Multiple Evolutionary Selections on Synonymous Codon Usage of Genes in the Mycoplasma bovis Genome

Mycoplasma bovis is a major pathogen causing arthritis, respiratory disease and mastitis in cattle. A better understanding of its genetic features and evolution might represent evidences of surviving host environments. In this study, multiple factors influencing synonymous codon usage patterns in M. bovis (three strains’ genomes) were analyzed. The overall nucleotide content of genes in the M. bovis genome is AT-rich. Although the G and C contents at the third codon position of genes in the leading strand differ from those in the lagging strand (p<0.05), the 59 synonymous codon usage patterns of genes in the leading strand are highly similar to those in the lagging strand. The over-represented codons and the under-represented codons were identified. A comparison of the synonymous codon usage pattern of M. bovis and cattle (susceptible host) indicated the independent formation of synonymous codon usage of M. bovis. Principal component analysis revealed that (i) strand-specific mutational bias fails to affect the synonymous codon usage pattern in the leading and lagging strands, (ii) mutation pressure from nucleotide content plays a role in shaping the overall codon usage, and (iii) the major trend of synonymous codon usage has a significant correlation with the gene expression level that is estimated by the codon adaptation index. The plot of the effective number of codons against the G+C content at the third codon position also reveals that mutation pressure undoubtedly contributes to the synonymous codon usage pattern of M. bovis. Additionally, the formation of the overall codon usage is determined by certain evolutionary selections for gene function classification (30S protein, 50S protein, transposase, membrane protein, and lipoprotein) and translation elongation region of genes in M. bovis. The information could be helpful in further investigations of evolutionary mechanisms of the Mycoplasma family and heterologous expression of its functionally important proteins.     

Analysis of synonymous codon usage patterns in the genus Rhizobium

The codon usage patterns of rhizobia have received increasing attention. However, little information is available regarding the conserved features of the codon usage patterns in a typical rhizobial genus. The codon usage patterns of six completely sequenced strains belonging to the genus Rhizobium were analysed as model rhizobia in the present study. The relative neutrality plot showed that selection pressure played a role in codon usage in the genus Rhizobium. Spearman’s rank correlation analysis combined with correspondence analysis (COA) showed that the codon adaptation index and the effective number of codons (ENC) had strong correlation with the first axis of the COA, which indicated the important role of gene expression level and the ENC in the codon usage patterns in this genus. The relative synonymous codon usage of Cys codons had the strongest correlation with the second axis of the COA. Accordingly, the usage of Cys codons was another important factor that shaped the codon usage patterns in Rhizobium genomes and was a conserved feature of the genus. Moreover, the comparison of codon usage between highly and lowly expressed genes showed that 20 unique preferred codons were shared among Rhizobium genomes, revealing another conserved feature of the genus. This is the first report of the codon usage patterns in the genus Rhizobium.