Email updates

Keep up to date with the latest news and content from Virology Journal and BioMed Central.

Open Access Research

Patterns and influencing factor of synonymous codon usage in porcine circovirus

Xin-sheng LIU, Yong-guang Zhang, Yu-zhen Fang and Yong-lu Wang*

Author Affiliations

State Key Laboratory of Veterinary Etiological Biology, National Foot and Mouth Disease Reference Laboratory, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou 730046, People's Republic of China

For all author emails, please log on.

Virology Journal 2012, 9:68  doi:10.1186/1743-422X-9-68

The electronic version of this article is the complete one and can be found online at: http://www.virologyj.com/content/9/1/68


Received:13 September 2011
Accepted:15 March 2012
Published:15 March 2012

© 2012 LIU et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Analysis of codon usage can reveal much about the molecular evolution of the viruses. Nevertheless, little information about synonymous codon usage pattern of porcine circovirus (PCV) genome in the process of its evolution is available. In this study, to give a new understanding on the evolutionary characteristics of PCV and the effects of natural selection from its host on the codon usage pattern of the virus, Patterns and the key determinants of codon usage in PCV were examined.

Methods

We carried out comprehensive analysis on codon usage pattern in the PCV genome, by calculating relative synonymous codon usage (RSCU), effective number of codons (ENC), dinucleotides and nucleic acid content of the PCV genome.

Results

PCV genomes have relatively much lower content of GC and codon preference, this result shows that nucleotide constraints have a major impact on its synonymous codon usage. The results of the correspondence analysis indicate codon usage patterns of PCV of various genotypes, various subgenotypes changed greatly, and significant differences in codon usage patterns of Each virus of Circoviridae.There is much comparability between PCV and its host in their synonymous codon usage, suggesting that the natural selection pressure from the host factor also affect the codon usage patterns of PCV. In particular, PCV genotype II is in synonymous codon usage more similar to pig than to PCV genotype I, which may be one of the most important molecular mechanisms of PCV genotype II to cause disease. The calculations results of the relative abundance of dinucleotides indicate that the composition of dinucleotides also plays a key role in the variation found in synonymous codon usage in PCV. Furthermore, geographic factors, the general average hydrophobicity and the aromaticity may be related to the formation of codon usage patterns of PCV.

Conclusion

The results of these studies suggest that synonymous codon usage pattern of PCV genome are the result of interaction between mutation pressure and natural selection from its host. The information from this study may not only have theoretical value in understanding the characteristics of synonymous codon usage in PCV genomes, but also have significant value for the molecular evolution of PCV.

Background

Genetic information is transmitted from mRNA to protein in a mode of triplet codon. Each amino acid matches with at least one codon, at most six codons. The codons encoding the same amino acid is called synonymous codon. During biosynthesis of protein, usage probability of those synonymous codons is different. Some species or some genes are usually prone to use one or several particular synonymous codons. These codons are called preferable codons, which is called as codon bias. Usage bias of codons from various species has been studied, and it is found that during protein biosynthesis synonymous codons encoding amino acid is not used randomly [1-3]. Many studies have indicated that obvious bias exists between different genes from different species or the same species [4-6]. Usage bias of codons is influenced mainly by mutation bias, translation selection, secondary protein structure, replication and selective transcription, hydrophobia and hydrophilia of protein, and external environment [7-13].

PCV belongs to genus of porcine circovirus, family of porcine circovirus. It has two genotype, PCV genotype I and PCV genotype II, and it is the smallest virus which has been discovered so far [14]. Among the different genotype, PCV genotype II infection and its related diseases have become one big problem across the globe for pig feeding, which threatens greatly to normal development of the industry of pig feeding. The PCV genome is a single-stranded negative circular DNA, and very small; full length of the PCV genotype I is only 1,759 bp, and PCV genotype II, 1,767 bp or 1,768 bp. The genome contains 11 open reading frames (ORF), among which, ORF1 encodes replication-associated proteins (Rep and Rep'); ORF2, structural proteins (viral capsid proteins, Cap); ORF3, toxicity-associated proteins, which can cause apoptosis [15,16]. By analyzing the whole sequence of PCV genome, it is found that ORF2 has smaller selective pressure than ORF1, and more mutation. Nucleotide sequences among various strains in the same genotype are very conservative, their homology is over 90%, while similarity between nucleotide sequences from various strains respectively from the two genotypes is less than 80% [17,18]. However, so far, studies have not related to usage of PCV codons. Explanation of codon usage pattern of PCV has significance on PCV evolution, gene prediction, gene classification, design of high expressed genes and viral vectors, and understanding of interaction between PCV and its host cells. Therefore, in this study, we first performed comprehensive analysis on codon usage pattern of PCV genome and the related factors affecting on codon usage. This study will play a major role in explanation of evolution process of PCV genome and further studies.

Results

The characteristics of synonymous codon usage in PCV

In order to investigate usage pattern of the PCV codons, we calculated various RSCU values of various codons from 28 different strains from different genotype. It can be seen from the three-dimension mesh plot of the analysis results of correspondence between 59 synonymous codons in the PCV genome (See Figure 1), range of Z axis (f'1) is between 15 to 2.5, which indicates that synonymous codons usage in the PCV genome is not balanced, that is to say, among all the 59 synonymous codons, a part of codon is rarely used, while others have higher usage frequency. Additionally, it can be seen from Table 1 that among 18 preferable codons, 12 ones have the end base of G or U, only 6 have the end base of A or C, and so those codons with the end base of G or U are prone to use in PCV genome. Nevertheless, compared with other vertebrate DNA viruses, PCV genome has lower GC%, from 48.35% to 49.12%, with an average content of 48.61% and SD value of 0.19 (Table 2). And hence, the phenomenon, that in PCV genome GC content is lower while the condons with the end base of G is used in a way of bias, suggests that content of G or C as the end base of codons has effect on usage pattern of synonymous codons. Apart from this, we can also see from Table 2, ENC values between PCV genomes has less fluctuation, with a range from 55.32 to 58.67, and an average value of 56.80 and SD value of 0.85, which indicates that codon bias of the PCV genome is stable.

thumbnailFigure 1. Three-dimension mesh plot of correspondence analysis value of 59 synonymous codons in PCV. It can be seen from the plot, change range on the Z axis(f'1)s from 15 to 2.5, which indicates PCV codon usage is not balanced.

Table 1. The relative synonymous codon usage frequency (RSCU) of PCV and swine

Table 2. Nucleotide content of 28 PCV genomes

Nucleotide content of all PCV genomes

Natural selection and mutation pressure has been considered to be two key factors which have effect on codon usage patterns of organisms [19]. In order to explore whether determinative factors for codon usage mutation in PCV is mutation pressure or natural selection, we compared correlation between A3%, U3%, G3%, C3%, GC3% and U3%, G3%, C3%, GC3% with correlation analysis (Table 3). The analysis results show, except GC%, GC3% has marked correlation with A%, U%, G% and C%. This indicates that GC3% can reflect interaction between natural selection and mutation pressure to some extent.

Table 3. The correlation analysis between the A, U, C, G contents and A3, U3, C3, G3 contents in all ORF of PCVa

In addition, the correlation between f'1,f'2 value and A%, U%, G%,C%, GC%, A3%, U3%, G3%, C3%, GC3% was analyzed (Table 4). The results showed that significant correlations exist between synonymous codon usage pattern and nucleotide composition in PCV. This result further verifies the conclusion that during the shaping of synonymous codon usage pattern of PCV, Composition constraints play and very important role.

Table 4. The correlation analysis between the first two axes in CA and the nucleotide contents of PCVa b

Genetic relationship based on synonymous codon usage

In order to compare synonymous codon usage patterns between different PCV genomes from different genotype, we carried out analysis on codon usage of different PCV genotype with correspondence analysis (CA). In correspondence analysis, the first dimension variable f'1 and the second dimension variable f'2 can reflect 39.87 and 23.83 percent of total mutation respectively. We can see from Figure 2, except for two strains of PCV genotype I in deviation from the cluster, other all the strains of PCV genotype II lie in the same cluster and overlap partly with each other. Obviously, PCV genotype I and PCV genotype II lay in two independent areas, which demonstrate that codon usage between the different PCV genotype is of great significance. Meanwhile, PCV genotype II-A and PCV genotype II-B almost lie in the same area, but they have tiny difference (Figure 2), which suggests that different sub-genotype from the same PCV genotype have difference in the aspect of codon usage.

thumbnailFigure 2. Correspondence analysis plot of relative synonymous codon usage in PCV. The first dimension variable f'1 and the second variable f'2 can reflect 39.87% and 23.83% of total mutation respectively.

Relationships between codon usage pattern of PCV and that of the host

From ENC values and corresponding relation distribution diagram of GC3% (Figure 3) we can see, most points are near or under the theoretical curve, which suggests that apart from mutation pressure which influences on codon usage pattern in PCV, the usage pattern is also influenced by other factors. As parasitic organisms, virus's codon usage pattern would be subject to its host to some extent [6]. In this study, patterns of codon usage are compared in PCV and its natural host, and then found that there are high similarities between them. In detail, the high frequently used codons in the swine were also the non-preferred codons of PCV, such as CUG, GUG, CAG, AAU, GAC, GAA, UGU, AGA and UUU. Further more, all preferentially used codons of the genome of PCV and swine were all G-ended or U-ended codons (Table 1). These results suggest that the selection pressure from the host affects codon usage pattern of PCV. It is worthwhile to note that PCV genotype II have high similarities with swine than PCV genotype I (Figure 4). In details, the values of RSCU in PCV genotype II and swine codon such as, AGA for Arg, GUG for Val, AGC for Ser, ACC for Thr were clearly different from that of PCV genotype I. It may be important one of Molecular mechanisms of infection and pathogenesis of PCV genotype II.

thumbnailFigure 3. The relationship between the effective number of codons (ENC) and the GC content of the third codon position (GC3%).

thumbnailFigure 4. Compare the codon the codon usage pattern among PCV genotype I, PCV genotype II and swine.

Relationship between dinucleotide biases and codon usage in PCV

Researches in recent years indicated that dinucleotide biases can affect codon bias [20]. To study the possible effect of dinucleotide composition on codon usage of the PCV genome, the relative abundances of the 16 dinucleotides in genome of the 28 PCV strains were calculated (Table 5). The result show that the occurrences of dinucleotides are not randomly distributed and no dinucleotides were present at the expected frequencies. The relative abundance of CpG showed a strong deviation from the "normal range" (mean ± SD = 0.622 ± 0.029) and was markedly under represented. Correspondingly, the synonymous codon containing the CpG were inhibited because of CpG was present at the expected frequencies. In detail, low RSCU values were present in 7 of all 8 codon containing CpG and they are all not preferentially used codons (such as CCG, UCG, ACG, CGC, CGG, CGA, CGU) (Table 1). These observations indicate that the composition of dinucleotides also plays a key role in the variation found in synonymous codon usage among PCV.

Table 5. Relative abundance of dinucleotides in PCV

Synonymous codon usage in different viruses of circoviridae is virus specific

The correspondence analysis has been performed in order to compare the synonymous codon usage pattern between the viruses of Circoviridae. From which we could detect one major trend in the f'1 which accounted for 20.27% of the total variation, and another major trend in the f'1 for 15.46%of the total variation. A plot of the f'1 and the f'1 of each virus of Circoviridae was shown in Figure 5. We can see from the plot that there were considerable differences for codon usage patterns among PCV, DCV, GCV, CoCV, CAV and BFDV. Concretely, PCV belong to the different genotype tends to come together. Moreover, DCV, GCV, CoCV and BFDV tend to come together and CAV was alone in separate area. It's clear that although different virus tends to come together, differences of synonymous codon usage pattern still exist between each virus. Therefore, the synonymous codon usage pattern of each virus of Circoviridae varies by the species of virus.

thumbnailFigure 5. Correspondence analysis plot of relative synonymous codon usage in Circoviridae . Each virus strain was divided by species.

Effect of other factors on codon usage

To investigate whether there is a correlation between the codon usage of PCV and geographic factor, 28 virus genes of PCV were divided into eleven groups according to obtained area, and correspondence analysis was also used. As can be seen from the plot, coordinate of virus isolates from different country is separated, and these relatively isolated spots tend to cluster into several groups according to the genotype (Figure 6). All above imply that these strains of PCV isolated from different places have different trend in codon usage variation. In addition, we performed another correlation analysis on f'1 in principal component analysis between GRAVY and the aromaticity score of each protein. From the result, we found that f'1 was high positive correlation with GRAVY (Spearman, r = 0.875, p < 0.001), and high negative correlation with aromaticity (Spearman, r = 0 905 p < 0.001).

thumbnailFigure 6. A plot of value of the first and second axis of each PCV strain for correspondence analysis. Each PCV strain was divided by geographical area.

Discussion

During protein biosynthesis synonymous codon encoding amino acids are not used randomly, and some species or some gene always prefers to use of one or several particular synonymous codons, which is called as codon usage bias. Precious studies reveal that different genes from different species or the same one have obvious codon usage bias [21,22]. Codon usage bias is influenced mainly by mutation bias [23,24], translation selection [25,26], secondary protein structure [20,27], replication and transcription selection [28], secondary mRNA structure [29], gene length [30], tRNA abundance [31], gene function and external environment [32]. However, most of these studies focus on some higher organism and many microorganisms with large genome and more genes, and there are few studies on virus with small genome and few genes or comparison between virus and host. Relatively, there are more reports on codon usage in genomes from viruses with great harm to mankind, such as SARS, human immunodeficiency virus, influenza virus A and hepatitis virus. PCV is a primary pathogen of postweaning multisystemic wasting syndrome (PMWS), which has threatened the development of pig feeding industry seriously because in recent year's occurrence of this disease has increased so as to bring about great economic loss in the world industry of pig feeding. Further studies on codon usage pattern in PCV have great significance on mutation pattern and molecular evolution of PCV. However, reports on codon usage pattern in PCV are rare, and this study is the first report.

By comparison with reported DNA viruses such as Duck plague virus, Duck enteritis virus, Iridovirus, Herpesvirus [33-36], synonymous codon usage bias in the PCV genome is low at large (average ENC values is 56.80, and minimum is 55.). This suggests that low codon bias may result from increase in itself replication efficiency in PCV in order to adapt to replication system of its hosts.

In this study, relation between main indices (f'1 and f'2)for the correspondence analysis on PCV usage cofon usage and its nucleotide composition (See Table 2) indicates, mutation pressure has a significant role in PCV codon usage. Other factors which can influence on PCV codon usage are also analysed and the initial results show that mutation pressure is the main factor to influence on PCV codon usage variation.

There were reports that natural selection can influence on synonymous codon usage pattern in viruses and the same conclusions are also obtained from this study. Three evidences support this conclusion. The first evidence is that PCV genome is GC3% -poor (average value = 47.08, SD = 2.88), but most of preferentially used codons are G/T-ended codons. Meanwhile, average of A3% is higher than that of T3%, but among the codons which PCV prefers to using, there are only three preferable codons with the end base of A3% while six those with the end base of T3%. The second evidence is that the high similarities exist between PCV and its natural host. The third evidence is that CpG and the synonymous codon including it were inhibited. The three above evidences both state that natural selection is involved in formation of PCV synonymous codon usage pattern.

At present, according to pathogenicity, antigenicity and nucleotide sequence difference, PCV is divided into two genotypes, PCV genotype I and PCV genotype II, of which PCV genotype II includes various subtypes. From significance of PCV codon usage between different genotypes in Figure 1, we can see that PCV codon bias may have association with genotypes. In addition to this, the results in this study also reveal that geological factor may almost have relation with codon usage in PCV. In some reports, gene length has certain correlation with codon usage [30]. Similarly, in some viruses, gene length has no effect on codon usage [22]. With correlation analysis we surveyed codon usage bias and gene length in PCV, and it is found that in these viral genes, codon usage bias has no notable correlation with gene length (Spearman, r = 0.075, p > 0.1). The results indicate that PCV gene length has no effect on synonymous codon usage. Other factors, including GRAVY and aromaticity may also significantly influence codon usage of PCV

Conclusions

Taken together, the codon usage patterns of PCV possibly result from interactions between natural selection and mutation pressure. These results not only provide an insight into the variation of codon usage pattern among the genomes of PCV, but also may help in understanding the processes governing the evolution of PCV.

Materials and methods

Sequence data

The information of 28 PCV genomes, including the genotype, length value, the isolated area and GenBank accession numbers of these strains was listed in the Table 6. In order to compare the differences between PCV and its host, twenty swine gene were gained and detailed information of these genes is listed in Table 7. In addition, to compare the codon usage patterns among different viruses, twenty-five viral genomes of Circoviridae were taken into account (Table 8). All of the sequences were downloaded from NCBI (http://www.ncbi.nlm.nih.gov/Genbank/ webcite). Each general nucleotide composition (T%,A%,C% and G%) and each nucleotide composition in the third site of codon (T3%,A3%,C3% and G3%) in PCV coding sequence were calculated by biosoftware DNAStar7.0 for windows.

Table 6. PCV genome sequences included in this study

Table 7. Swine genes used in this study

Table 8. Viral sequence of Circoviridae used in this study

Synonymous codon usage measures

In order to eliminate the influence of amino acid composition on codon usage and directly reflect the usage characteristics of codon, the study evaluates synonymous codon usage bias through statistical estimation on relative synonymous codon usage frequency (RSCU) [37]. RSCU value refers to the ratio between the usage frequency of one codon in gene sample and expected frequency in the synonymous codon family. If the synonymous codon usage of one amino acid has no preferences, that is, codon usage frequency is close to expected frequency, the RSCU values of codons are equal to 1; if a codon RSCU value is greater than 1, the codon use frequency is higher than expected frequency, whereas it is less than expected value.

The definition on a single gene codon bias is mainly based on effective number of codons (ENC) [38]. ENC values can reflect the preference degree of synonymous codon non-equilibrium use in codon family. The range of ENC values is from 20 (each amino acid only uses one codon) to 61 (all synonymous codons are equivalently used). ENC value is closer to 20, the degree of being used non-randomly is higher, and the bias is stronger. It is generally believed that the genes are provided with significant codon bias when ENC ≤ 35. The values of RSCU and ENC were obtained by codonW program.

A comparison of actual and expected dinucleotide frequencies of the 16 dinucleotides in coding region of PCV genomes was also undertaken using SPSS 17.0.

Correspondence analysis (CA)

Correspondence analysis is mainly used for detecting the changes of codon RSCU values in genes [39]. It is an effective multivariate statistical method of studying the internal relation between the variables and samples, and it is successfully applied to the study of codon. In correspondence analysis, all genes in samples are distributed in a 59-dimensional (59 justice codons, in addition to the stop codon, Met, and Trp) vector space, each gene is described with 59 (f'1, f'2,..., f'59) variables, the results can be applied for finding out the major factors affecting codon usage bias in genes [40,41]. This was done using the CodonW program.

Correlation analysis

Correlation analysis of PCV was used to identify the relationship between nucleotide composition and synonymous codon usage pattern [42]. All statistical processes were carried out by with statistical software SPSS17.0 for windows.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

XS L and YL W conceived of the study. XS L downloaded these sequences, calculated the data, analyzed the results and drafted the manuscript; YZ F assisted with data analysis; YG Z supervised the research and helped draft the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study was supported in parts by grants from National Pig Industrial System (CARS-36-06B) and Research and Demonstration on evaluation technologies of clinical immune responses of serious animal diseases vaccine (201203039)

References

  1. Dittmar KA, Goodenbour JM, Pan J: Tissue-specific differences in human transfer RNA expression.

    PLoS Genet 2006, 2:2107-2115. OpenURL

  2. Lloyd AT, Sharp PM: Evolution of codon usage patterns: the extent and nature of divergence between Candida albicans and Saccharomyces cerevisiae.

    Nucleic Acids Res 1992, 20:5289-5295. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Xie T, Ding D, Tao X, Dafu D: The relationship between synonymous codon usage and protein structure.

    FEBS Lett 1998, 434:93-96. PubMed Abstract | Publisher Full Text OpenURL

  4. Chiapello H, Lisacek F, Caboche M: Codon usage and gene function are related in sequences of Arabidopsisthaliana.

    Gene 1998, 209(1-2):GC1-GC38. PubMed Abstract | Publisher Full Text OpenURL

  5. Adams MJ, Antoniw JF: Codon usage bias amongst plant viruses.

    Arch Virol 2003, 149(1):113-135. PubMed Abstract | Publisher Full Text OpenURL

  6. Zhou H, Wang H, Huang LF: Heterogeneity in codon usages of sobemovirus genes.

    Arch Virol 2005, 150(8):1591-1605. PubMed Abstract | Publisher Full Text OpenURL

  7. Levin DB, Whittome B: Codon usage in nucleopolyhedroviruses.

    J Gen Virol 2000, 81:2313-2325. PubMed Abstract | Publisher Full Text OpenURL

  8. Gupta SK, Ghosh TC: Gene expressivity is the main factor in dictating the codon usage variation among the genes in Pseudomonas aeruginosa.

    Gene 2001, 273(1):63-70. PubMed Abstract | Publisher Full Text OpenURL

  9. D'Onofrio G, Ghosh TC, Bernardi G: The base composition of the genes is correlated with the secondary structures of the encoded Proteins.

    Gene 2002, 300(1):179-187. PubMed Abstract | Publisher Full Text OpenURL

  10. Gu W, Zhou T, Ma J, Sun X, Lu Z: The relationship between synonymous codon usage and Protein structure in Escherichia coli and Homo sapiens.

    Biosystems 2004, 73(2):89-97. PubMed Abstract | Publisher Full Text OpenURL

  11. Wang Meng, Zhang Jie, Zhou Jian-hua, Chen Hao-tai, Ma Li-na, Ding Yao-zhong, Liu Wen-qian, Liu Yong-sheng: Analysis of codon usage in bovine viral diarrhea virus.

    Arch Virol 2010, 156(1):153-160. PubMed Abstract | Publisher Full Text OpenURL

  12. Romero H, Zavala A, Musto H: Compositional pressure and translational selection determine codon usage in the extremely GC poor unicellular eukaryote Entamoeba histolytica.

    Gene 2000, 242(1-2):307-311. PubMed Abstract | Publisher Full Text OpenURL

  13. Van der Linden MG, de Farias ST: Correlation between codon usage and thermostability.

    Extremophiles 2006, 10(5):479-481. PubMed Abstract | Publisher Full Text OpenURL

  14. Mankertz A, Persson F, Mankertz J, Blaess G, Buhk HJ: Mapping and characterization of the origin of DNA replication of porcine circovirus.

    J Virol 1997, 71(3):2562-2566. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Karuppannan AK, Kwang J: ORF3 of porcine circovirus 2 enhances the in vitro and in vivo spread of the of the virus.

    Virology 2011, 410:248-256. PubMed Abstract | Publisher Full Text OpenURL

  16. Liu J, Chen I, Du Q: The ORF3 protein of porcine circovirus type 2 is involved in viral pathogenesis in vivo.

    Virol 2006, 80(10):5065-5073. Publisher Full Text OpenURL

  17. Andre L, Hamel , Lihua L, Lin , Gopi PS, Nayar : Nucleotide sequence of porcine circovirus associated with postweaning multisystemic wasting syndrome in pigs.

    J Virol 1998, 72:5262-5267. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Meehan BM, McNeilly F, Todd D, Kennedy S, Jewhurst VA, Ellis JA, Hassard LE, Clark EG, Haines DM, Allan GM: Characterization of novel circovirus DNAs associated with wasting syndromes in pigs.

    J Gen Virol 1998, 79:2171-2179. PubMed Abstract | Publisher Full Text OpenURL

  19. Zhou T, Wanjun Gu, Ma J, Sun X, Zuhong Lu: Analysis of synonymous codon usage in H5N1 virus and other influenza A viruses.

    Biosystems 2005, 81:77-86. PubMed Abstract | Publisher Full Text OpenURL

  20. Chiusano Maria Luisa, Alvarez-Valin Fernando, Giulio Massimo Di, D'Onofrio Giuseppe, Ammirato Gaetano, Colonna Giovanni, Bernardi Giorgio: Second codon positions of genes and the secondary structures of proteins.Relationships and implications for the origin of the genetic code.

    Gene 2000, 261:63-69. PubMed Abstract | Publisher Full Text OpenURL

  21. Sharp PM, Tuohy TM, Mosurski KR: Codon usage in yeast:cluster analysis clearly differentiates highly and lowly expressed genes.

    Nucleic Acids Res 1986, 14:5125-5143. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Gu WJ, Zhou T, Ma JM, Sun X, Lu ZH: Analysis of synonymous codon usage in SARS Coronavirus and other viruses in the Nidovirales.

    Virus Res 2004, 101(2):155-161. PubMed Abstract | Publisher Full Text OpenURL

  23. Liu Wen-qian, Zhang Jie, Zhang Yi-qiang, Zhou Jian-hua, Chen Hao-tai, Ma Li-na, Ding Yao-zhong, Liu Yongsheng: Compare the differences of synonymous codon usage between the two species within cardiovirus.

    Virol J 2011, 8:325. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Gareth M, Jenkins , Pagel Mark, Ernest A, Gould , Zanotto DeA, Paolo M, Edward C, Holmes : Evolution of base composition and codon usage bias in the genus Flavivirus.

    J Mol Evol 2001, 52:383-390. PubMed Abstract | Publisher Full Text OpenURL

  25. Peixoto L, Zavala A, Romero H, Musto H: The strength of translational selection for codon usage varies in the three replicons of Sinorhizobium meliloti.

    Gene 2003, 320(27):109-116. PubMed Abstract | Publisher Full Text OpenURL

  26. Romero H, Zavala A, Musto H, Bernardi G: The influence of translational selection on codon usage in fishes from the family cyprinidae.

    Gene 2003, 317(1-2):141-147. PubMed Abstract | Publisher Full Text OpenURL

  27. Gupta SK, Majumdar S, Bhattacharya TK, Ghosh TC: Studies on the relationships between the synonymous codon usage and protein secondary structural units.

    Biochem Biophys Res Commun 2000, 269:692-696. PubMed Abstract | Publisher Full Text OpenURL

  28. James O, McInerney : Replicational and transcriptional selection on codon usage in Borrelia burgdorferi.

    PNAS 1998, 95(18):10698-10703. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Zama M: Codon usage and secondary strueture of mRNA.

    Nucleic Acids Symp Ser 1990, 22:93-94. PubMed Abstract OpenURL

  30. Stoletzki N, Eyre-Walker A: Synonymous codon usage in Escherichia coli: selection for translational accuracy.

    Mol Biol Evol 2007, 24(2):374-381. PubMed Abstract | Publisher Full Text OpenURL

  31. Rocha PC, Eduardo : Codon usage bias from tRNA's point of view: redundancy, specialization, and efficient decoding for translation optimization.

    Genome Res 2004, 14(11):2279-2281. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Lynn DJ, GAC: Synonymous codon usage is subject to selection in thermophilic bacteria.

    Nucleic Acids Res 2002, 30(19):4272-4277. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Minghui Fu: Codon usage bias in herpesvirus.

    Arch Virol 2010, 155:391-396. PubMed Abstract | Publisher Full Text OpenURL

  34. Tsai C-T, Lin C-H, Chang C-Y: Analysis of codon usage bias and base compositional constraints in iridovirus genomes.

    Virus Res 2007, 126:196-206. PubMed Abstract | Publisher Full Text OpenURL

  35. Cai M-S, An-Chun Cheng, Wang M-S, Li-Chan Zhao: Characterization of Synonymous Codon Usage Bias in the Duck Plague Virus UL35 Gene.

    Intervirology 2009, 52:266-278. PubMed Abstract | Publisher Full Text OpenURL

  36. Jia R, Cheng A, Wang M: Analysis of synonymous codon usage in the UL24 gene of duck enteritis virus.

    Virus Genes 2009, 38:96-103. PubMed Abstract | Publisher Full Text OpenURL

  37. Wright F: The effective number of codons used in a gene.

    Gene 1990, 87:23-29. PubMed Abstract | Publisher Full Text OpenURL

  38. Tao P, Dai L, Luo M, Tien Fangqiang Tang Po, Pan Z: Analysis of synonymous codon usage in classical swine fever virus.

    Virus Genes 2009, 38:104-112. PubMed Abstract | Publisher Full Text OpenURL

  39. Hao S, Zhang Q, Liu X, Wang X, Zhang H, Wu Y, Jiang F: Analysis of synonymous codon usage in 11 Human Bocavirus isolates.

    Biosystems 2008, 92:207-214. PubMed Abstract | Publisher Full Text OpenURL

  40. Sau K, Gupta SK, Sau S, Mandal SC, Ghosh TC: Factors influencing synonymous codon and amino acid usage biases in Mimivirus.

    Biosystems 2006, 85:107-113. PubMed Abstract | Publisher Full Text OpenURL

  41. Ewens WJ, Grant GR: Statistical Methods in Bioinformatics. New York: Springer; 2001. OpenURL

  42. Drake JW, Holland JJ: Mutation rates among RNA viruses.

    Proc Natl Acad Sci USA 1999, 96:13910-13913. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL