Analysis of codon usage and nucleotide composition bias in polioviruses
Virology Journal volume 8, Article number: 146 (2011)
Poliovirus, the causative agent of poliomyelitis, is a human enterovirus and a member of the family of Picornaviridae and among the most rapidly evolving viruses known. Analysis of codon usage can reveal much about the molecular evolution of the viruses. However, little information about synonymous codon usage pattern of polioviruses genome has been acquired to date.
The relative synonymous codon usage (RSCU) values, effective number of codon (ENC) values, nucleotide contents and dinucleotides were investigated and a comparative analysis of codon usage pattern for open reading frames (ORFs) among 48 polioviruses isolates including 31 of genotype 1, 13 of genotype 2 and 4 of genotype 3.
The result shows that the overall extent of codon usage bias in poliovirus samples is low (mean ENC = 53.754 > 40). The general correlation between base composition and codon usage bias suggests that mutational pressure rather than natural selection is the main factor that determines the codon usage bias in those polioviruses. Depending on the RSCU data, it was found that there was a significant variation in bias of codon usage among three genotypes. Geographic factor also has some effect on the codon usage pattern (exists in the genotype-1 of polioviruses). No significant effect in gene length or vaccine derived polioviruses (DVPVs), wild viruses and live attenuated virus was observed on the variations of synonymous codon usage in the virus genes. The relative abundance of dinucleotide (CpG) in the ORFs of polioviruses are far below expected values especially in DVPVs and attenuated virus of polioviruses genotype 1.
The information from this study may not only have theoretical value in understanding poliovirus evolution, especially for DVPVs genotype 1, but also have potential value for the development of poliovirus vaccines.
When molecular sequence data started to be accumulated nearly 20 years ago, it was noted that synonymous codons are not used equally in different genomes, even in different genes of the same genome[1–3]. As an important evolutionary phenomenon, it is well known that synonymous codon usage bias exists in a wide range of biological systems from prokaryotes to eukaryotes [4, 5]. Codon usage analysis has been applied to prokaryote and eukaryote, such as Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Caenorhabditis elegans and human beings [6–8]. These observed patterns in synonymous codon usage varied among genes within a genome, and among genomes. The codon usage is attributable to the equilibrium between natural selection and mutation pressure [9, 10]. Recent studies of viral codon usage has shown that mutation bias may be a more important factor than natural selection in determining codon usage bias of some viruses, such as Picornaviridae, Pestivirus, plant viruses, and vertebrate DNA viruses [9, 11–13]. Meanwhile, recent report also showed that the G+C compositional constraint is the main factor that determines the codon usage bias in iridovirus genomes[11, 14]. Analysis of codon usage can reveal much about the molecular evolution or individual genes of the viruses.
Polioviruses belong to the family Picornaviridae and are classified as human enterovirus C (HEV-C) species in the genus Enterovirus according to the current taxonomy [15, 16]. Polioviruses can be divided into three different genotypes: 1, 2 and 3. The genome of each genotypes contains a single positive-stranded RNA with a size of approximately 6 kb consisting of a single large open reading frame (ORF) flanked by 5' and 3' untranslated region .
As we known, the Sabin oral poliovaccine (OPV) was among the best known viral vaccines . It has saved the lives and health of innumerable people, in particular children. However, poliovirus is highly genetically variable. OPV viruses may undergo transformation into circulating highly diverged VDPV, exhibiting properties hardly distinguishable from those of wild polioviruses . So far, little information about synonymous codon usage pattern of polioviruses genome has been acquired to date. To our knowledge, this is the first report of the codon usage analysis on polioviruses (including wild strains, attenuated live vaccine strains and VDPV strains). In this study, we analyzed the codon usage data and base composition of 48 available representative complete ORFs of poliovirus to obtain some clues to the features of genetic evolution of the virus.
A total of 48 poliovirus genomes were used in this study (Table 1). The serial number (SN), genotype, length value, isolated region, GenBank accession numbers, and other detail information about these strains were listed in Table 1. All of the sequences were downloaded from NCBI http://www.ncbi.nlm.nih.gov/Genbank/, and 48 poliovirus genomes were selected in the study. The other sequences with >98% sequence identities were excluded.
The actual and predicted values of the effective number of codon (ENC)
The ENC is used to measure the degree of departure from the equal use of synonymous codons of coding regions of polioviruses. The values of the effective number of codon (ENC) range from 20 to 61. In an extremely biased gene where only one codon is used for each amino acid, this value would be 20; if all codons are used equally, it would be 61; and if the value of ENC is greater than 40, the codon usage bias was regarded as low. The values of ENC were obtained by EMBOSS CHIPS program .
Genes, whose codon choice is constrained only by a mutation bias, will lie on or just below the curve of the predicted values. The predicted values of ENC were calculated as
where s represents the given (G+C)3% value .
The calculation of the relative synonymous codon usage (RSCU)
To investigate the pattern of relative synonymous codon usage (RSCU) without the influence of amino acid composition among all polioviruses samples, the RSCU values of codons in each ORF of polioviruses were calculated according to the formula of previous reports [22, 23].
where g ij is the observed number of the i th codon for j th amino acid which has n i type of synonymous codons. The codon with RSCU value more than 1.0 has positive codon usage bias, while the value <1.0 has relative negative codon usage bias. When RSCU value is equal to 1.0, it means that this codon is chosen equally and randomly.
Relative dinucleotide abundance in polioviruses
Because dinucleotide biases can affect codon bias, the relative abundance of dinucleotides in the coding regions of polioviruse genomes was assessed using the method described by Karlin and Burge . A comparison of actual and expected dinucleotide frequencies of the 16 dinucleotides in coding region of the 48 polioviruses genomes was also undertaken. The odds ratio ρxy=ƒxy/ƒyƒx, where ƒx denotes the frequency of the nucleotide X, ƒy denotes the frequency of the nucleotide Y, ƒyƒx the expected frequency of the dinucleotide XY and ƒxy the frequency of the dinucleotide XY, etc., for each dinucleotide were calculated. As a conservative criterion, for ρxy > 1.23 (or < 0.78), the XY pair is considered to be of over-represented (or under-represented) relative abundance compared with a random association of mononucleotides.
Principal component analysis (PCA) was carried out to analyze the major trend in codon usage pattern in different genomes of polioviruses (excluding non-coding regions). It is a statistical method that performs linear mapping to extract optimal features from an input distribution in the mean squared error sense and can be used by self-organizing neural networks to form unsupervised neural preprocessing modules for classification problems . In order to minimize the effect of amino acid composition on codon usage, each ORF is represented as a 59-dimensional vector. Each dimension corresponds to the RSCU value of one sense codon excluding Met, Trp and three stop codons.
Linear regression analysis was used to find the correlation between codon usage bias and gene length. Correlation analysis is used to identify the relationship between codon usage bias and synonymous codon usage pattern. This analysis is implemented based on the Spearman's rank correlation analysis way.
All statistical analyses were carried out using the statistical analysis software SPSS Version 17.0.
The characteristics of synonymous codon usage in polioviruses
In order to investigate the extent of codon usage bias in polioviruses, all RSCU values of different codon in 48 polioviruses strains were calculated. There is only two preferred codons UUG (Leu) and GUG (Val), choosing G at the third position, and most of preferred codons are ended with A (Table 2). Moreover, polioviruses genome is A redundant with A content ranging from 29.739 to 30.826.11, with the mean value of 30.367 and S.D. of 0.234; in contrast, low content of G ranging from 21.723 to 22.401 (mean = 22.118, S.D. of 0.147), suggesting that nucleotide contents influence the patterns of synonymous codon usage (Table 3). The values of ENC among these polioviruses ORFs are similar, which vary from 52.609 to 55.105 with a mean value of 53.754 and S.D. of 0.545. The data showed that the extent of codon preference in polioviruses genes was kept basically stable.
Compositional properties of ORFs of 48 polioviruses genomes
The values of A, U, C, G and C+G were compared with the values of A3, C3, G3, U3, (G+C) 3, respectively. An interesting and complex correlation was observed. In detail, the (C+G)3 have highly significant correlations with A, U, C, G and C+G, respectively, indicating C+G may reflect interaction between mutation pressure and natural selection. However, the A have no correlation with A3, G3 and C3, and U have no correlation with A3 (Table 4). Both cases suggested that the nucleotide constraint possibly influence synonymous codon usage of polioviruses. In addition, the correlation between the Axis 1 (calculated by PCA) and the values of A, C, G, U, A3, C3, G3, U3, (G+C), (G+C)3 of each strain was also analyzed. The significant correlation was found between nucleotide compositions and synonymous codon usage to some extent excluding Axis 1 and the value of A (Table 4). The analysis revealed that most of the codon usage bias among ORFs of polioviruses strains was directly related to the base composition. Finally, the ENC-plot [ENC plotted against (G+C)3%] was used as a part of general strategy to investigate patterns of synonymous codon usage and all of the spots lie below the expected curve (Figure 1). These imply that the codon bias can be explained mainly by an uneven base composition, in other words, by mutation pressure rather than natural selection.
Effect of other potential factors on codon usage
Principal component analysis was carried out to identify the codon usage bias among ORFs. From which we could detect one major trend in the Axis 1 which accounted for 20.815% of the total variation, and another major trend in the Axis 2 for 16.273% of the total variation. A plot of the Axis 1 and the Axis 2 of each gene was shown in Additional file 1, Figure S1. Obviously, those polioviruses belong to the same genotype tends to come together (except strain 48, isolated from Finland). Compared with the scattered groups of polioviruses genotype 1, genotype 2 and 3 strains aggregated more tightly to some degree. Although this graph is a little complex, it seems that there is a clear geographical demarcation in the polioviruses genotype 1 such as the VDPV strains isolated from USA, Dominica, China mainland and Taiwan. These may indicate that geographic is another factor on codon usage bias.
The frequencies of occurrence for dinucleotides were not randomly distributed and no dinucleotides were present at the expected frequencies. And the frequency of CpG and TpA was significantly low at all codon positions for coding region of 48 polioviruses genomes (mean ± S.D. = 0.490 ± 0.012; and mean ± S.D. = 0.748 ± 0.034. both < 0.78). The relative abundance of CpA and TpG also showed slight deviation from the ''normal range'' (mean ± S.D. = 1.253 ± 0.032 and 1.423 ± 0.023, respectively) (Table 5). In addition, the RSCU values of the eight codons containing CpG (CCG, GCG, UCG, ACG, CGC, CGG, CGU, and CGA) were analyzed, to reveal the possible effects of CpG under-represented on codon usage bias. All of these eight codons were not preferential codons and were markedly suppressed. The six codons containing TpA (UUA, CUA, GUA, UAU, UAC and AUA) were suppressed too. Conversely, the RSCU values of the eight codons containing CpA (UCA, CCA, ACA, GCA, CAA, CAG, CAU, CAC) and five codons containing UpG (UUG, GUG, UGU, UGC, CUG) are high, and most of them (8 out of 13) were preferential codons (Table 2 and Table 5). In addition, compared with DVPVs and live attenuated strain of polioviruses genotype 1, the wild viruses has higher frequencies of dinucleotides including CpG (Figure 2 and Table 5).
Furthermore, we also performed a linear regression analysis on ENC value and gene length of ORFs of 48 polioviruses genomes. However, there was no significant correlation between codon usage and gene length in these virus genes (Spearman P > 0.05).
Studies of synonymous codon usage in viruses can reveal much about viral genomes . The overall codon usage among 48 ORFs of polioviruses was analyzed in this study. First, the ENC values of all the poliovirus samples were analyzed, and the results showed that the majority of polioviruses do not have a strong codon bias (mean ENC = 53.754 > 40). In addition, together with published data on codon usage bias among some RNA viruses, such as BVDV, H5N1 influenza virus and SARS-covs with mean values of 51.43, 50.91 and 48.99, respectively, one possible explanation for this is that the weak codon bias of RNA virus is advantageous to replicate efficiently in vertebrate host cells, with potentially distinct codon preferences [26–28].
Natural selection and mutation pressure are thought to be the main factors that account for codon usage variation in different organisms [29–31]. In this study, the general association between codon usage bias and base composition suggests that mutational pressure, rather than natural selection is the mainly factors on codon usage pattern of polioviruses.
Codon usage can also be strongly influenced by underlying biases in dinucleotide frequency, which differs greatly among organisms. Specifically, after accounting for dinucleotide biases, the proportion of codon usage bias explained by mutation pressure often increases, as seen in human RNA viruses . Our study revealed that CpG and the eight CpG-containing codons are notably deficient in ORFs of 48 poliovirus genomes. The explanation for the CpG deficiency is immunologic escape. A high CpG content may be detrimental to small DNA (or RNA) viruses, as unmethylated CpGs are recognized by the host's innate immune system (Toll-like receptor 9) as a pathogen signature . As with vertebrate genomes, methylated viral genomes would face a high chance of mutation at CpGs, that would result in a reduction of this dinucleotide [9, 33]. We found that DVPVs and live attenuated virus of genotype 1 have lower frequencies of CpG dinucleotide compare with wild viruses of polioviruses genotype 1. The most popular explanation for lower frequencies of CpG in ORFs of DVPV genomes is that when OPV viruses turning into VDPV genotype 1, a lower frequencies of CpG dinucleotide maybe help VDPV out of the host immunity.
Although it seems speculative and complex, some researchers have found that reduction of the rate of poliovirus protein synthesis through large-scale utilization of codons that are not optimal has caused attenuation of viral virulence by lowering specific infectivity . Therefore, the information from this study may not only have theoretical value in understanding poliovirus evolution (especially for DVPVs genotype 1), but also have practical value for the development the poliovirus vaccine. However, a more comprehensive analysis is needed to reveal more information about codon usage bias variation within poliovirus and other responsible factors.
The information from this study may not only help to understand the evolution of the poliovirus, especially for DVPVs genotype 1, but also have potential value for the development of poliovirus vaccines.
Grantham R, Gautier C, Gouy M, Mercier R, Pave A: Codon catalog usage and the genome hypothesis. Nucleic Acids Res 1980, 8: r49-r62.
Ingvarsson PK: Molecular evolution of synonymous codon usage in Populus. BMC Evol Biol 2008, 8: 307. 10.1186/1471-2148-8-307
Supek F, Vlahovicek K: Comparison of codon usage measures and their applicability in prediction of microbial gene expressivity. BMC Bioinformatics 2005, 6: 182. 10.1186/1471-2105-6-182
Archetti M: Codon usage bias and mutation constraints reduce the level of error minimization of the genetic code. J Mol Evol 2004, 59: 258-266. 10.1007/s00239-004-2620-0
Liu X, Wu C, Chen AY: Codon usage bias and recombination events for neuraminidase and hemagglutinin genes in Chinese isolates of influenza A virus subtype H9N2. Arch Virol 2010, 155: 685-693. 10.1007/s00705-010-0631-2
Kanaya S, Kinouchi M, Abe T, Kudo Y, Yamada Y, Nishi T, Mori H, Ikemura T: Analysis of codon usage diversity of bacterial genes with a self-organizing map (SOM): characterization of horizontally transferred genes with emphasis on the E. coli O157 genome. Gene 2001, 276: 89-99. 10.1016/S0378-1119(01)00673-4
Karlin S, Mrazek J: What drives codon choices in human genes? J Mol Biol 1996, 262: 459-472. 10.1006/jmbi.1996.0528
Sharp PM, Cowe E, Higgins DG, Shields DC, Wolfe KH, Wright F: Codon usage patterns in Escherichia coli, Bacillus subtilis, Saccharomyces cerevisiae, Schizosaccharomyces pombe, Drosophila melanogaster and Homo sapiens; a review of the considerable within-species diversity. Nucleic Acids Res 1988, 16: 8207-8211. 10.1093/nar/16.17.8207
Shackelton LA, Parrish CR, Holmes EC: Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses. J Mol Evol 2006, 62: 551-563. 10.1007/s00239-005-0221-1
Sharp PM, Stenico M, Peden JF, Lloyd AT: Codon usage: mutational bias, translational selection, or both? Biochem Soc Trans 1993, 21: 835-841.
Zhong J, Li Y, Zhao S, Liu S, Zhang Z: Mutation pressure shapes codon usage in the GC-Rich genome of foot-and-mouth disease virus. Virus Genes 2007, 35: 767-776. 10.1007/s11262-007-0159-z
Tao P, Dai L, Luo M, Tang F, Tien P, Pan Z: Analysis of synonymous codon usage in classical swine fever virus. Virus Genes 2009, 38: 104-112. 10.1007/s11262-008-0296-z
Fu M: Codon usage bias in herpesvirus. Arch Virol 2010, 155: 391-396. 10.1007/s00705-010-0597-0
Tsai CT, Lin CH, Chang CY: Analysis of codon usage bias and base compositional constraints in iridovirus genomes. Virus Res 2007, 126: 196-206. 10.1016/j.virusres.2007.03.001
Savolainen-Kopra C, Blomqvist S: Mechanisms of genetic variation in polioviruses. Rev Med Virol 2010, 20: 358-371. 10.1002/rmv.663
Wimmer E, Hellen CU, Cao X: Genetics of poliovirus. Annu Rev Genet 1993, 27: 353-436. 10.1146/annurev.ge.27.120193.002033
Kew OM, Sutter RW, de Gourville EM, Dowdle WR, Pallansch MA: Vaccine-derived polioviruses and the endgame strategy for global polio eradication. Annu Rev Microbiol 2005, 59: 587-635. 10.1146/annurev.micro.58.030603.123625
Minor PD, Almond J: Poliovirus vaccines: molecular biology and immune response. Molecular biology of picornaviruses 2002, 381-390.
Agol VI: Vaccine-derived polioviruses. Biologicals 2006, 34: 103-108. 10.1016/j.biologicals.2006.02.007
Lu H, Zhao WM, Zheng Y, Wang H, Qi M, Yu XP: Analysis of synonymous codon usage bias in Chlamydia. Acta Biochim Biophys Sin (Shanghai) 2005, 37: 1-10. 10.1093/abbs/37.1.1
Wright F: The 'effective number of codons' used in a gene. Gene 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9
Sharp PM, Li WH: Codon usage in regulatory genes in Escherichia coli does not reflect selection for 'rare' codons. Nucleic Acids Res 1986, 14: 7737-7749. 10.1093/nar/14.19.7737
Zhou JH, Zhang J, Chen HT, Ma LN, Liu YS: Analysis of synonymous codon usage in foot-and-mouth disease virus. Vet Res Commun 2010, 34: 393-404. 10.1007/s11259-010-9359-4
Karlin S, Burge C: Dinucleotide relative abundance extremes: a genomic signature. Trends Genet 1995, 11: 283-290. 10.1016/S0168-9525(00)89076-9
Jenkins GM, Holmes EC: The extent of codon usage bias in human RNA viruses and its evolutionary origin. Virus Res 2003, 92: 1-7. 10.1016/S0168-1702(02)00309-X
Wang M, Zhang J, Zhou JH, Chen HT, Ma LN, Ding YZ, Liu WQ, Liu YS: Analysis of codon usage in bovine viral diarrhea virus. Arch Virol 2011, 156: 153-160. 10.1007/s00705-010-0848-0
Gu W, Zhou T, Ma J, Sun X, Lu Z: Analysis of synonymous codon usage in SARS Coronavirus and other viruses in the Nidovirales. Virus Res 2004, 101: 155-161. 10.1016/j.virusres.2004.01.006
Zhou T, Gu W, Ma J, Sun X, Lu Z: Analysis of synonymous codon usage in H5N1 virus and other influenza A viruses. Biosystems 2005, 81: 77-86. 10.1016/j.biosystems.2005.03.002
Schultz ST, Willis JH: Individual variation in inbreeding depression: the roles of inbreeding history and mutation. Genetics 1995, 141: 1209-1223.
Sharp PM, Li WH: The codon Adaptation Index--a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res 1987, 15: 1281-1295. 10.1093/nar/15.3.1281
Sniegowski PD, Gerrish PJ, Johnson T, Shaver A: The evolution of mutation rates: separating causes from consequences. Bioessays 2000, 22: 1057-1066. 10.1002/1521-1878(200012)22:12<1057::AID-BIES3>3.0.CO;2-W
Krieg AM: CpG DNA: trigger of sepsis, mediator of protection, or both? Scand J Infect Dis 2003, 35: 653-659. 10.1080/00365540310015999
Ambinder RF, Robertson KD, Tao Q: DNA methylation and the Epstein-Barr virus. Semin Cancer Biol 1999, 9: 369-375. 10.1006/scbi.1999.0137
Mueller S, Papamichail D, Coleman JR, Skiena S, Wimmer E: Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity. J Virol 2006, 80: 9687-9696. 10.1128/JVI.00738-06
This work was supported in parts by grants from National Key Technologies R&D Program (2009ZX08007-006B), International Science & Technology Cooperation Program of China (No.2010DFA32640) and Science and Technology Key Project of Gansu Province (No.0801NKDA034). This study was also supported by National Natural Science foundation of China (No.30700597 and No. 31072143).
The authors declare that they have no competing interests.
JZ, MW and YL designed the study and drafted the manuscript. WL, JZ, HC and LM collected the data and participated in the sequence alignment. YD and YG performed the statistical analysis. All authors read and approved the final manuscript.
Jie Zhang, Meng Wang contributed equally to this work.
Electronic supplementary material
About this article
Cite this article
Zhang, J., Wang, M., Liu, Wq. et al. Analysis of codon usage and nucleotide composition bias in polioviruses. Virol J 8, 146 (2011). https://doi.org/10.1186/1743-422X-8-146