Bioinformatics analysis of rabbit haemorrhagic disease virus genome
Virology Journalvolume 8, Article number: 494 (2011)
Rabbit haemorrhagic disease virus (RHDV), as the pathogeny of Rabbit haemorrhagic disease, can cause a highly infectious and often fatal disease only affecting wild and domestic rabbits. Recent researches revealed that it, as one number of the Caliciviridae, has some specialties in its genome, its reproduction and so on.
In this report, we firstly analyzed its genome and two open reading frameworks (ORFs) from this aspect of codon usage bias. Our researches indicated that mutation pressure rather than natural is the most important determinant in RHDV with high codon bias, and the codon usage bias is nearly contrary between ORF1 and ORF2, which is maybe one of factors regulating the expression of VP60 (encoding by ORF1) and VP10 (encoding by ORF2). Furthermore, negative selective constraints on the RHDV whole genome implied that VP10 played an important role in RHDV lifecycle.
We conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. According to the results of the principal component analysis for ORF2 of RSCU, we firstly separated 30 RHDV into two genotypes, and the ENC values indicated ORF1 and ORF2 were independent among the evolution of RHDV.
Synonymous codons are not used randomly . The variation of codon usage among ORFs in different organisms is accounted by mutational pressure and translational selection as two main factors [2, 3]. Levels and causes of codon usage bias are available to understand viral evolution and the interplay between viruses and the immune response . Thus, many organisms such as bacteria, yeast, Drosophila, and mammals, have been studied in great detail up on codon usage bias and nucleotide composition . However, same researches in viruses, especially in animal viruses, have been less studied. It has been observed that codon usage bias in human RNA viruses is related to mutational pressure, G+C content, the segmented nature of the genome and the route of transmission of the virus . For some vertebrate DNA viruses, genome-wide mutational pressure is regarded as the main determinant of codon usage rather than natural selection for specific coding triplets . Analysis of the bovine papillomavirus type 1 (BPV1) late genes has revealed a relationship between codon usage and tRNA availability . In the mammalian papillomaviruses, it has been proposed that differences from the average codon usage frequencies in the host genome strongly influence both viral replication and gene expression . Codon usage may play a key role in regulating latent versus productive infection in Epstein-Barr virus . Recently, it was reported that codon usage is an important driving force in the evolution of astroviruses and small DNA viruses [10, 11]. Clearly, studies of synonymous codon usage in viruses can reveal much about the molecular evolution of viruses or individual genes. Such information would be relevant in understanding the regulation of viral gene expression.
Up to now, little codon usage analysis has been performed on Rabbit haemorrhagic disease virus (RHDV), which is the pathogen causing Rabbit haemorrhagic disease (RHD), also known as rabbit calicivirus disease (RCD) or viral haemorrhagic disease (VHD), a highly infectious and often fatal disease that affects wild and domestic rabbits. Although the virus infects only rabbits, RHD continues to cause serious problems in different parts of the world. RHDV is a single positive stranded RNA virus without envelope, which contains two open reading frames (ORFs) separately encoding a predicted polyprotein and a minor structural protein named VP10 . After the hydrolysis of self-coding 3C-like cysteinase, the polyprotein was finally hydrolyzed into 8 cleavage products including 7 nonstructural proteins and 1 structural protein named as VP60 [13, 14]. Studies on the phylogenetic relationship of RHDVs showed only one serotype had been isolated, and no genotyping for RHDV was reported. It reported that the VP10 was translated with an efficiency of 20% of the preceding ORF1 . In order to better understand the characteristics of the RHDV genome and to reveal more information about the viral genome, we have analyzed the codon usage and dinucleotide composition. In this report, we sought to address the following issues concerning codon usage in RHDV: (i) the extent and causes of codon bias in RHDV; (ii) A possible genotyping of RHDV; (iii) Codon usage bias as a factor reducing the expression of VP10 and (iiii) the evolution of the ORFs.
2. Materials and methods
The 30 available complete RNA sequences of RHDV were obtained from GenBank randomly in January 2011. The serial number (SN), collection dates, isolated areas and GenBank accession numbers are listed in Table 1.
2.2 The relative synonymous codon usage (RSCU) in RHDV
To investigate the characteristics of synonymous codon usage without the influence of amino acid composition, RSCU values of each codon in a ORF of RHDV were calculated according to previous reports (2 Sharp, Tuohy et al. 1986) as the followed formula:
Where gij is the observed number of the i th codon for j th amino acid which has ni type of synonymous codons. The codons with RSCU value higher than 1.0 have positive codon usage bias, while codons with value lower than 1.0 has relative negative codon usage bias. As RSCU values of some codons are nearly equal to 1.0, it means that these codons are chosen equally and randomly.
2.3 The content of each nucleotides and G+C at the synonymous third codon position (GC3s)
The index GC3s means the fraction of the nucleotides G+C at the synonymous third codon position, excluding Met, Trp, and the termination codons.
2.4 The effective number of codons (ENC)
The ENC, as the best estimator of absolute synonymous codon usage bias , was calculated for the quantification of the codon usage bias of each ORF . The predicted values of ENC were calculated as
where s represents the given (G+C)3% value. The values of ENC can also be obtained by EMBOSS CHIPS program .
2.5 Dn and ds of two ORFs
Analyses were conducted with the Nei-Gojobori model , involving 30 nucleotide sequences. All positions containing gaps and missing data were eliminated. The values of dn, ds and ω (dn/ds) were calculated in MEGA4.0 .
2.6 Correspondence analysis (COA)
Multivariate statistical analysis can be used to explore the relationships between variables and samples. In this study, correspondence analysis was used to investigate the major trend in codon usage variation among ORFs. In this study, the complete coding region of each ORF was represented as a 59 dimensional vector, and each dimension corresponds to the RSCU value of one sense codon (excluding Met, Trp, and the termination codons) .
2.7 Correlation analysis
Correlation analysis was used to identify the relationship between nucleotide composition and synonymous codon usage pattern . This analysis was implemented based on the Spearman's rank correlation analysis way.
All statistical processes were carried out by with statistical software SPSS 17.0 for windows.
3.1 Measures of relative synonymous codon usage
The values of nucleotide contents in complete coding region of all 30 RHDV genomes were analyzed and listed in Table 2 and Table 3. Evidently, (C+G)% content of the ORF1 fluctuated from 50.889 to 51.557 with a mean value of 51.14557, and (C+G)% content of the ORF2 were ranged from 35.593 to 40.113 with a mean value of 37.6624, which were indicating that nucleotides A and U were the major elements of ORF2 against ORF1. Comparing the values of A3%, U3%, C3% and G3%, it is clear that C3% was distinctly high and A3% was the lowest of all in ORF1 of RHDV, while U3% was distinctly high and C3% was the lowest of all in ORF2 of RHDV. The (C3+G3) % in ORF1 fluctuated from 57.014 to 58.977 with a mean value of 57.68287 and (C3+G3)% were range from 31.356 to 39.831 with a mean value of 34.8337. And the ENC values of ORF1 fluctuated from 54.192 to 55.491 with a mean value of 54.95 and ENC values of ORF2 displayed a far-ranging distribution from 39.771 to 51.964 with a mean value of 44.46. The ENC values of ORF1 were a little high indicating that there is a particular extent of codon preference in ORF1, but the codon usage is relatively randomly selected in ORF2 on the base of ENC values. The details of the overall relative synonymous codon usage (RSCU) values of 59 codons for each ORF in 30 RHDV genomes were listed in Table 4. Most preferentially used codons in ORF1 were C-ended or G-ended codons except Ala, Pro and Ser, however, A-ended or G-ended codons were preferred as the content of ORF2.
In addition, the dn, ds and ω(dN/dS) values of ORF1 were separately 0.014, 0.338 and 0.041, and the values of ORF2 were 0.034, 0.103 and 0.034, respectively. The ω values of two ORFs in RHDV genome are generally low, indicating that the RHDV whole genome is subject to relatively strong selective constraints.
3.2 Correspondence analysis
COA was used to investigate the major trend in codon usage variation between two ORFs of all 30 RHDV selected for this study. After COA for RHDV Genome, one major trend in the first axis (f'1) which accounted for 42.967% of the total variation, and another major trend in the second axis (f'2) which accounted for 3.632% of the total variation. The coordinate of the complete coding region of each ORF was plotted in Figure 1 defining by the first and second principal axes. It is clear that coordinate of each ORF is relatively isolated. Interestingly, we found that relatively isolated spots from ORF2 tend to cluster into two groups: the ordinate value of one group (marked as Group 1) is positive value and the other one (marked as Group 2) is negative value. Interestingly, all of those strains isolated before 2000 belonged to Group 2.
3.3 Correlation analysis
To estimate whether the evolution of RHDV genome on codon usage was regulated by mutation pressure or natural selection, the A%, U%, C%, G% and (C+G)% were compared with A3%, U3%, C3%, G3% and (C3+G3)%, respectively (Table 5). There is a complex correlation among nucleotide compositions. In detail, A3%, U3%, C3% and G3% have a significant negative correlation with G%, C%, U% and A% and positive correlation with A%, U%, C% and G%, respectively. It suggests that nucleotide constraint may influence synonymous codon usage patterns. However, A3% has non-correlation with U% and C%, and U3% has non-correlation with A% and G%, respectively, which haven't indicated any peculiarity about synonymous codon usage. Furthermore, C3% and G3% have non-correlation with A%, G% and U%, C%, respectively, indicating these data don't reflect the true feature of synonymous codon usage as well. Therefore, linear regression analysis was implemented to analyze the correlation between synonymous codon usage bias and nucleotide compositions. Details of correlation analysis between the first two principle axes (f'1 and f'2) of each RHDV genome in COA and nucleotide contents were listed in Table 6. In surprise, only f2 values are closely related to base nucleotide A and G content on the third codon position only, suggesting that nucleotide A and G is a factor influencing the synonymous codon usage pattern of RHDV genome. However, f'1 value has non-correlation with base nucleotide contents on the third codon position; it is observably suggest that codon usage patterns in RHDV were probably influenced by other factors, such as the second structure of viral genome and limits of host. In spite of that, compositional constraint is a factor shaping the pattern of synonymous codon usage in RHDV genome.
There have been more and more features that are unique to RHDV within the family Caliciviridae, including its single host tropism, its genome and its VP10 as a structural protein with unknown function. After we analyzed synonymous codon usage in RHDV (Table 2), we obtained several conclusions and conjectures as followed.
4.1 Mutational bias as a main factor leading to synonymous codon usage variation
ENC-plot, as a general strategy, was utilized to investigate patterns of synonymous codon usage. The ENC-plots of ORFs constrained only by a C3+G3 composition will lie on or just below the curve of the predicted values . ENC values of RHDV genomes were plotted against its corresponding (C3+G3) %. All of the spots lie below the curve of the predicted values, as shown in Figure 2, suggesting that the codon usage bias in all these 30 RHDV genomes is principally influenced by the mutational bias.
4.2 A proof for codon usage bias as a factor reducing the expression of VP10
As we know, the efficiency of gene expression is influenced by regulator sequences or elements and codon usage bias. It reported that the RNA sequence of the 3-terminal 84 nucleotides of ORF1were found to be crucial for VP10 expression instead of the encoded peptide. VP10 coding by ORF2 has been reported as a low expressive structural protein against VP60 coding by ORF1 . And its efficiency of translation is only 20% of VP60. According to results showed by Table 4, it revealed the differences in codon usage patterns of two ORFs, which is a possible factor reducing the expression of VP10.
4.3 Negative selective constraints on the RHDV whole genome
Although VP10 encoded by ORF2, as a minor structural protein with unknown functions, has been described by LIU as a nonessential protein for virus infectivity, the ω value of ORF2 suggests VP10 plays an important role in the certain stage of whole RHDV lifecycle. After combining with low expression and ω value of VP10, we conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. This mechanism has been confirmed in various positive-chain RNA viruses, including coxsackievirus, dengue virus, equine arterivirus, foot-and-mouth disease virus, hepatitis C virus, poliovirus, rhinovirus, and severe acute respiratory syndrome [23–29], although the details remain elusive.
4.4 Independent evolution of ORF1 and ORF2
As preceding description, ENC reflects the evolution of codon usage variation and nucleotide composition to some degree. After the correlation analysis of ENC values between ORF1 and ORF2 (Table 7), the related coefficient of ENC values of two ORFs is 0.230, and p value is 0.222 more than 0.05. These data revealed that no correlation existed in ENC values of two ORFs, indicating that codon usage patterns and evolution of two ORFs are separated each other. Further, this information maybe helps us well understand why RSCU and ENC between two ORFs are quite different.
4.5 A possible genotyping basis
Interestingly, we found that relatively isolated spots from ORF2 tend to cluster into two groups: the ordinate value of one group (marked as Group 1) is positive value and the other one (marked as Group 2) is negative value. And all of those strains isolated before 2000 belonged to Group 2, including Italy-90, RHDV-V351, RHDV-FRG, BS89, RHDV-SD and M67473.1. Although RHDV has been reported as only one type, this may be a reference on dividing into two genotypes.
In this report, we firstly analyzed its genome and two open reading frameworks (ORFs) from this aspect of codon usage bias. Our researches indicated that mutation pressure rather than natural is the most important determinant in RHDV with high codon bias, and the codon usage bias is nearly contrary between ORF1 and ORF2, which is maybe one of factors regulating the expression of VP60 (encoding by ORF1) and VP10 (encoding by ORF2). Furthermore, negative selective constraints on the RHDV whole genome implied that VP10 played an important role in RHDV lifecycle. We conjectured that VP10 might be beneficial for the replication, release or both of virus by inducing infected cell apoptosis initiate by RHDV. According to the results of the principal component analysis for ORF2 of RSCU, we firstly separated 30 RHDV into two genotypes, and the ENC values indicated ORF1 and ORF2 were independent among the evolution of RHDV. All the results will guide the next researches on the RHDV as a reference.
Karlin S, Mrazek J: What drives codon choices in human genes? J Mol Biol 1996,262(4):459-472. 10.1006/jmbi.1996.0528
Sharp PM, Tuohy TM, Mosurski KR: Codon usage in yeast: cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Res 1986,14(13):5125-5143. 10.1093/nar/14.13.5125
Lesnik T, Solomovici J, Deana A, Ehrlich R, Reiss C: Ribosome traffic in E-coli and regulation of gene expression. Journal of Theoretical Biology 2000,202(2):175-185. 10.1006/jtbi.1999.1047
Shackelton LA, Parrish CR, Holmes EC: Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses. J Mol Evol 2006,62(5):551-563. 10.1007/s00239-005-0221-1
Mooers AO, Holmes EC: The evolution of base composition and phylogenetic inference. Trends Ecol Evol 2000,15(9):365-369. 10.1016/S0169-5347(00)01934-0
Jenkins GM, Holmes EC: The extent of codon usage bias in human RNA viruses and its evolutionary origin. Virus Res 2003,92(1):1-7. 10.1016/S0168-1702(02)00309-X
Zhou J, Liu WJ, Peng SW, Sun XY, Frazer I: Papillomavirus capsid protein expression level depends on the match between codon usage and tRNA availability. J Virol 1999,73(6):4972-4982.
Zhao KN, Liu WJ, Frazer IH: Codon usage bias and A+T content variation in human papillomavirus genomes. Virus Res 2003,98(2):95-104. 10.1016/j.virusres.2003.08.019
Karlin S, Blaisdell BE, Schachtel GA: Contrasts in codon usage of latent versus productive genes of Epstein-Barr virus: data and hypotheses. J Virol 1990,64(9):4264-4273.
Sewatanon J, Srichatrapimuk S, Auewarakul P: Compositional bias and size of genomes of human DNA viruses. Intervirology 2007,50(2):123-132. 10.1159/000098238
Van Hemert FJ, Berkhout B, Lukashov VV: Host-related nucleotide composition and codon usage as driving forces in the recent evolution of the Astroviridae. Virology 2007,361(2):447-454. 10.1016/j.virol.2006.11.021
Wirblich C, Thiel HJ, Meyers G: Genetic map of the calicivirus rabbit hemorrhagic disease virus as deduced from in vitro translation studies. J Virol 1996,70(11):7974-7983.
Wirblich C, Sibilia M, Boniotti MB, Rossi C, Thiel HJ, Meyers G: 3C-like protease of rabbit hemorrhagic disease virus: identification of cleavage sites in the ORF1 polyprotein and analysis of cleavage specificity. J Virol 1995,69(11):7159-7168.
Meyers G, Wirblich C, Thiel HJ, Thumfart JO: Rabbit hemorrhagic disease virus: genome organization and polyprotein processing of a calicivirus studied after transient expression of cDNA constructs. Virology 2000,276(2):349-363. 10.1006/viro.2000.0545
Meyers G: Translation of the minor capsid protein of a calicivirus is initiated by a novel termination-dependent reinitiation mechanism. J Biol Chem 2003,278(36):34051-34060. 10.1074/jbc.M304874200
Comeron JM, Aguade M: An evaluation of measures of synonymous codon usage bias. J Mol Evol 1998, 47: 268-274. 10.1007/PL00006384
Wright F: The "effective number of codons" used in a gene. Gene 1990, 87: 23-29. 10.1016/0378-1119(90)90491-9
Lu H, Zhao WM, Zheng Y, Wang H, Qi M, Yu XP: Analysis of Synonymous Codon Usage Bias in Chlamydia . Acta Biochim Biophys Sin 2005, 37: 1-10. 10.1093/abbs/37.1.1
Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Molecular Biology and Evolution 1986, 3: 418-426.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution 2007, 24: 1596-1599. 10.1093/molbev/msm092
Mardia KV, Kent JT, Bibby JM: Multivariate analysis. NewYork. Academic press; 1979.
Ewens WJ, Grant GR: Statistical Methods in Bioinformatics. NewYork. Springer; 2001.
Egger D, Wolk B, Gosert R, Bianchi L, Blum HE, Moradpour D, Bienz K: Expression of Hepatitis C virus proteins induced distinct membrane alterations including a candidate viral replication complex. J Virol 2002,76(12):5974-5984. 10.1128/JVI.76.12.5974-5984.2002
Espert L, Codogno P, Biard-Piechaczyk K: Involvement of autophagy in viral infections: antiviral function and subversion by viruses. J Mol Med 2007, 85: 811-823. 10.1007/s00109-007-0173-6
Jackson WT, Giddings TH, Taylor MP, Mulinyawe S, Rabinovitch M, Ropito R, Kirkegaard K: Subversion of cellular autophagosomal machinery by RNA viruses. PLoS Biol 2005,3(5):e156. 10.1371/journal.pbio.0030156
Lee YR, Lei HY, Liu MT, Wang JR, Chen SH, Jiang-Shieh YF, Lin YS, Yeh TM, Liu CC, Liu HS: Autophagic machinery activated by dengue virus enhances virus replication. Virology 2008, 374: 240-248. 10.1016/j.virol.2008.02.016
Pedersen KW, van der Meer Y, Roos N, Snijder EJ: Open reading frame 1aencoded subunits of the arterivirus replicase induce endoplasmic reticulum derived double-membrane vesicles which carry the viral replication complex. J Virol 1999,73(3):2016-2026.
Wileman T: Aggresomes and autophagy generate sites for virus replication. Science 2006, 312: 875-878. 10.1126/science.1126766
Wong J, Zhang J, Xiaoning S, Gao G, Mao I, McManus BM, Luo H: Autophagosome supports Coxsackievirus B3 replication in host cells. J Virol 2008,82(18):9143-9153. 10.1128/JVI.00641-08
This work was supported by the fund of Special Social Commonweal Research Programs for Research Institutions (2005DIB4J041, China).
The authors declare that they have no competing interests.
XTT and BYL contributed equally to the original draft of the manuscript, and approved the final version. ZL and WQJ contributed to conception and design of the manuscript, and revised the manuscript. LJX is the corresponding author. All authors have read and approved the final manuscript.
Xiao-ting Tian, Bao-yu Li contributed equally to this work.