- Open Access
Analysis of the use of codon pairs in the HE gene of the ISA virus shows a correlation between bias in HPR codon-pair use and mortality rates caused by the virus
© Tello et al.; licensee BioMed Central Ltd. 2013
- Received: 28 January 2013
- Accepted: 28 May 2013
- Published: 6 June 2013
Segment 6 of the ISA virus codes for hemoagglutinin-esterase (HE). This segment is highly variable, with more than 26 variants identified. The major variation is observed in what is called the high polymorphism region (HPR). The role of the different HPR zones in the viral cycle or evolution remains unknown. However viruses that present the HPR0 are avirulent, while viruses with important deletions in this region have been responsible for outbreaks with high mortality rates. In this work, using bioinformatic tools, we examined the influence of different HPRs on the adaptation of HE genes to the host translational machinery and the relationship to observed virulence.
Translational efficiency of HE genes and their HPR were estimated analyzing codon-pair bias (CPB), adaptation to host codon use (codon adaptation index - CAI) and the adaptation to available tRNAs (tAI). These values were correlated with reported mortality for the respective ISA virus and the ΔG of RNA folding. tRNA abundance was inferred from tRNA gene numbers identified in the Salmo salar genome using tRNAScan-SE. Statistical correlation between data was performed using a non-parametric test.
We found that HPR0 contains zones with codon pairs of low frequency and low availability of tRNA with respect to salmon codon-pair usage, suggesting that HPR modifies HE translational efficiency. Although calculating tAI was impossible because one third of tRNAs (~60.000) were tRNA-ala, translational efficiency measured by CPB shows that as HPR size increases, the CPB value of the HE gene decreases (P = 2x10-7, ρ = −0.675, n = 63) and that these values correlate positively with the mortality rates caused by the virus (ρ = 0.829, P = 2x10-7, n = 11). The mortality associated with different virus isolates or their corresponding HPR sizes were not related with the ΔG of HPR RNA folding, suggesting that the secondary structure of HPR RNA does not modify virulence.
Our results suggest that HPR size affects the efficiency of gene translation, which modulates the virulence of the virus by a mechanism similar to that observed in production of live attenuated vaccines through deoptimization of codon-pair usage.
- Codon Usage
- Translation Efficiency
- Synonymous Codon
- Infectious Salmon Anemia Virus
- Amino Acid Pair
The infectious salmon anemia virus (ISAV) is an orthomyxovirus, composed of eight segments. . Two phylogenetic groups have been identified, the European and the North American, which probably diverged at the beginning of the last century [2, 3]. Within each group the viruses present minimal genetic differences with the exception of certain regions of high variability localized in genes encoding for glycoprotein F  and HE. The region of hyper-variability in HE genes (HPR) is located between the transmembrane region and the globular domain with hemagglutinin-esterase activity [5, 6]. The function of HPR within the viral cycle is unknown. However, the type of HPR in the virus has been associated with the degree of observed virulence. For example, the HPR0 is found in avirulent ISAV strains, while HPR7b is associated with ISAV outbreaks in Chile with high mortality rates [7–9]. The molecular ratios that explain the relationship between the HPR and observed virulence have not yet been identified. However, the HPR does not modify esterase or hemagglutininin activity .
Far from being silent, synonymous mutations affect diverse cellular processes, such as translation kinetics, mRNA structure and protein folding . These mutations are distributed differently in the genes according to their levels of expression. Highly expressed genes tend to use preferentially the most abundant synonymous codons [12, 13]. The relationship between bias in synonymous codon usage and gene expression seems to be the consequence of the synchrony observed between the frequency of codon usage and the abundance of tRNA that recognizes it [12, 14–19].
As well as the bias in codon usage, there is a bias in codon-pair usage, which is characteristic for each organism and independent of the frequency of codon usage . As with the synonymous codons, the codon pairs are distributed in the genes differently according to their level of gene expression .
The ribosome traffic jam in the mRNA has been suggested as a common mechanism that explains the reduction in translation efficiency by introducing codons with low tRNA availability, or tRNAs with little compatibility among themselves, as well as the negative relationship between the formation of stable structures in the mRNA and the level of protein expression [22–25].
There are few examples among viruses of the influence of codon use on the viral cycle. Studies about the adaptation to host codon usage indicate that viral genes that codify for critical proteins tend to use the synonymous codons most represented in the host genome . In Influenza A, the structure acquired by the 5’end of the codifying zone is fundamental for viral packaging [27, 28]. In the hepatitis A virus, the presence of underrepresented codons in the host genome seems to be a strategy to avoid competition with cellular mRNAs and to increase adequate folding of the capsid protein [29, 30]. In the influenza A and polio viruses, the introduction of synonymous codon pairs in the viral genome with a negative bias produces viable attenuated viral particles with reduced rates of protein synthesis and replication [31–33].
It has been suggested that the different HPR are probably the result of in-frame deletions from sectors of HPR0 [6, 34]. These deletions could modify translation efficiency by elimination of underrepresented codons in a similar but opposite mechanism to what is observed in the generation in vitro of live attenuated vaccines by deoptimization of their codon-pair usage. In this work we analyze this hypothesis of the effect of HPR on the adaptation to host codon usage, tRNA availability, host codon-pair usage bias and the folding energy of the HE gene.
The effect of the HPR on translation efficiency
Translation efficiency is governed by secondary structure of the mRNA  and by its adaptation to the tRNA pool . To evaluate if the HPR affects virulence through the possible formation of a stable secondary structure that reduces the advance of ribosome, thus causing a decrease in the HE translation rate, we determined if there is a relationship between the HPR folding energy and observed virulence. Our results indicate the absence of a statistically significant relationship between folding ΔG and the size of the HPR region or associated mortality. However, the HPR0 region showed the most favorable energy (ΔG = −24.2 Kcal/mol) (Additional file 1: Figure S1).
The effect of the presence of the different HPRs on translation efficiency can also be evaluated by analyzing the degree of adaptation to the tRNA pool through the use of the tAI value, which corresponds to the geometric mean of the relative abundance of tRNA in the cells that recognize the codons present in the gene. The relative abundance of the tRNA can be calculated based on the number of copies of each tRNA present in the genome . Using the tRNAScan-SE program, we identified 63857 genes in the Salmo salar genome that codify for tRNAs. Approximately 21085 of these genes codify for pseudogenes and another 20546 codify for alanine tRNA. Of the rest of the tRNAs, their abundances range from 2724 for methionine tRNA to 236 for tryptophan tRNA. As well, 353 untyped tRNA were identified and 20 genes for tRNA suppressors (Additional file 2: Table S1). As in the calculation of tAI, the relative abundance of the tRNAs that recognize the codons in the gene was normalized with respect to the codon with greater availability of tRNAs. The extreme disproportion between Ala-tRNA copies and the other tRNA for the other amino acids makes interpretation of the data problematic and could lead to erroneous conclusions. To avoid this problem, we continue the study using adaptation to host codon use and bias in codon-pair usage.
Adaptation to host codon usage and host codon pair bias of the ISAV hemagglutinin-esterase gene
The codon adaptation index (CAI) is a good indicator of the adaptation to the host codon usage and a good predictor of translation efficiency . Previous studies have indicated that ISAV is the least adapted orthomyxovirus to it host’s genome, sharing CAI values with less than 3% of the genes of Salmo salar (Tello, unpublished). Among these, the North American virus presents the lowest level of adaptation. If the HPR modifies translation efficiency through the accumulation of high or low frequency codons, it can be expected that the CAI values of hemagglutinin genes are proportional to the size of the HPR and the protein.
When viruses were compared according to their country of isolation, there were significant differences (t-test, P = 0.008) in the mean CPB values of the HE gene of the Norwegian (CPB = 0.0394±0.007) and Chilean viruses (CPB = 0.0446±0.006), with the Chilean isolates presenting the higher CPB values (Additional file 1: Figure S4). When the CPB values of the HE gene are classified according to their HPR type, we note that the lowest CPB values are found in the HE genes containing HPR21 and HPR0, while the highest values are found in the HE genes containing HPR1, HPR15 and HRP7 (Additional file 1: Figure S5).
Relationship among codon-pair bias, HPR type and virulence
Our analysis consistently shows that the genes containing the HPR0 and HPR7 regions have extreme CPB values. The two HPRs are associated with opposite phenotypes: the HPR0 region is found in avirulent strains of ISAV, while the HPR7 region has been systematically found in Chilean outbreaks that have been responsible for eliminating 90% of salmon in production centers. To determine if there is a relationship among virulence, HPR type and codon-pair bias, we correlated our CPB values to the mortality rates reported by Mjeeland .
Extrapolation of mortality based on the CPB value of the HPR shows that the Norwegian viruses with HPR15 (AF378180.2 and AF302799.1) and the Chilean viruses with HPR1C (GU830908.1 and FJ594285.1) have an estimated mortality > 80%, which is higher than values observed in ISAV2 with HPR01. When we analyzed the CPB of the complete HE gene, we observed that most of the Chilean viruses have values between those of ISAV1 (HPR14) and ISAV2 (HPR01). The CPB value of the HE gene from Chilean isolates with HRP7b has a value of 0.122, which by extrapolation becomes a virus with similar characteristics to HPR1 (CPB = 0.146) in the isolates of ISAV2. ISAV2 has the highest mortality rate by ISAV injection (75%) and the second highest rate by cohabitation (28.9%), followed by the HPR14 in the ISAV1 isolate (46.1%) (Additional file 3: Table S2).
The results presented in this work indicate the existence of a relationship between bias in the use of codon pairs and the size of the HPR region or the HE gene. In particular, our analysis shows that the different types of HPR tend to eliminate regions that present underrepresented codon pairs (negative CPS index). In Influenza A and polio, the introduction of codon pairs with low CPS values is an effective way to reduce translation of the viral gene, reducing the production of viral particles and virulence. In contrast, the introduction of codon pairs with positive CPS increases protein expression [31, 32]. The molecular explanation for this relationship is not clear, but it is probably because not all tRNA are compatible among themselves upon simultaneously occupying the A and P ribosome sites, as multivariate analysis suggests .
A universal codon distribution pattern has recently been described. In general the genes tend to accumulate codons poorly adapted to their tRNA pool among the first 50 codons, while the more adapted codons are found among the last 150 to 300 nucleotides . It has been suggested that the presence of this adaptation ramp to the tRNA pool serves to avoid the ribosome traffic jam in the mRNA, the translation efficiency being inverse to the quantity of ribosomes in the mRNA. Thus, the presence of underrepresented codons could reduce translational efficiency and the advance and turn-over of ribosomes [22–24].
The HPR region is located at approximately 80 codons from the carboxyl terminal with respect to the gene with HPR0. According to the ramp observed in eukaryotic , this region should contain codons more adapted to the tRNA pool. The presence of underrepresented codon pairs in the HPR zone could decrease the velocity of the advance of ribosomes by the incompatibility among tRNA, being functionally equivalent to the use of codons with low levels of available tRNA. On the other hand, the elimination of these codon pairs could increase translation efficiency by decreasing the ribosome traffic jam in this region. The Chilean isolates with HPR7b present a ΣCPS profile with high CPS values and a continuous ramp toward the carboxyl terminal (Figure 3), which is interrupted in the virus with HPR0. Consistent with this, viruses containing HPR0 are associated with avirulent ISAV strains. On the other hand, HPR regions where codon pairs with negative CPS have been eliminated are associated with outbreaks with high mortality rates, such as HPR1 or HPR7b found in 80% of the pathogenic outbreaks in Chile [38, 39]. The virulence values extrapolated for the Chilean isolates are consistent with the effect of these viruses on the Chilean salmon industry.
Current efforts to cultivate the virus in vitro with the HPR0 genotype have not been fruitful, which is why it remains unknown if these viruses have lower HE protein synthesis rates than viruses with the HPR7b genotype. Nevertheless, the correlations between the CPB values in the HE gene with the size of the region and the degree of pathogenicity of the ISA virus are consistent with the effect of HPR on translation efficiency.
Both immune escape due to antigenic variation produced by different HPRs  and modification of hemagglutinin-esterase activity can act as selection forces for the different HPR described. Recent experiments have determined that HPR does not modify hemagglutinin-esterase activity , suggesting that the biological activity of the protein does not explain the difference in observed virulence among ISAV isolates. There are variations in the protein sequence of the glycoproteins of influenza A that give rise to antigenic drift . However, in influenza A, the variations occur by specific mutations that produce variations in the epitopes, in contrast to what occurs in ISAV, where the greatest variation observed in the HPR occurs by partial deletions in the HPR0 . In ISAV, lesser variations have been observed on the surface of the HE protein . However, it is not known if these mutations have an important effect on antigen-antibody recognition. Although antigenic variation in the ISAV hemagglutinin-esterase could act as a selective force, in our opinion this concept makes more sense in populations whose individuals are long-lived and acquire immunity after surviving a viral infection or vaccination. Unlike humans, salmon are short-lived (four years) and suffer outbreaks of ISAV with high levels of mortality. Moreover, variants in the HPR region appeared prior to the development of vaccines. Determining the degree of adaptation to the tRNA pool could help in confirming that the codons in the HPR region reduce translation efficiency. However, the estimation of the number of tRNA copies in Salmo salar, using the tRNAScan-SE program (which recognizes close to 100% of the tRNA in a genome), showed that Salmo salar has a large number of tRNA genes (>65000), which is much higher than the approximately 500 copies in humans and is on the order of the number found in other teleosts, like zebra fish (~45,000) . Interestingly, close to a third of the identified tRNA is tRNA-Ala and another third is pseudogenes. This disproportionately high level of tRNA suggests the presence of a repeated SINE-type element, which presents homologous regions to those of the tRNA . Repeated SINE-type elements have been described in salmonoids and other teleosts, but with leucine and lysine tRNA [44–47]. This suggests the presence of a new SINE-type in Salmo salar, although further studies will be necessary once the S. salar genome has been completely annotated to characterize this type of element. Despite the significant disproportion of Ala-tRNA compared to the other tRNAs, the W values of the codons in the HPR0 region show three areas, close to the area with negative CPS, where codons with low tRNA availability accumulate (Additional file 4: Table S3, Additional file 5: Table S4 and Additional file 1: Figure S8). This reinforces the idea that HPR0 is a zone of reduced translation efficiency.
In summary, our analysis shows a correlation between the CPB value of the HE gene or the HPR region and the virulence of the ISA virus. To date, we have not determined if the CPB value in the ISAV is related to gene translation. However, in vitro data from influenza A and the poliovirus strongly suggest the existence of this relationship. Although it is an over-simplification of the pathogenesis to explain virulence based solely on the adaptation of viral genes to the host’s codon pair usage, to our knowledge the correlations described in this work constitute the only natural example of de-optimization of codon pairs as the element correlating to the level of pathogenicity of a virus and constitute a reliable explanation of the role of the HPR region in the ISA virus as a modulating factor of virulence. Without doubt, the development of a reverse genetic technique could verify this relationship, but until this technique becomes available, sequence analysis is a useful tool.
Calculation of tAI
The adaptation of a gene to the tRNA pool is calculated according to equation 2, in which wi is defined as the quotient between Wi and Wmax (Wi/Wmax), ikg is the codon defined by the kth triplet of gene g, and lg is the length of gene g in codifying codons.
Determination of bias in codon-pair usage
Calculation of the CAI
The codon adaptation index was calculated with the Emboss program  based on the frequency of codon usage of Salmo salar.
Prediction of free folding energy
Free folding energy was calculated with the Mfold program through the web server (http://mfold.rna.albany.edu/?q=mfold) using the standard parameters. The energy was calculated for the region between residuals 330 and 371 with respect to HE SK779 HPR0.
The virulence of 11 ISA viruses with known genomes was obtained from data published by Mjaaland et al. in 2005, in which they reported mortality rates 60 days after infection. This work analyzed the mortality rates of two fish populations from the same parents and the mortality rates of fish with genetic heterogeneity. The fish with the same parents were infected by cohabitation (TGI and TGII) with genetically heterogeneous fish that were previously infected by intraperitoneal injection (COH) .
HE genes of ISAV
The complete sequence of the codifying regions of the ISAV HE genes was downloaded from Genbank. All the genes that did not begin with ATG were filtered out, ending with a stop codon and ambiguous nucleotides. A total of 63 sequences passed through the filter. The GenBank codes of these sequences can be found in the supplementary material.
Graphs and statistical methods
The graphs and statistical tests were made using the program SigmaPlot 11. The significance of the differences or correlations among the data groups obtained was evaluated with non-parametric tests (Rank Sum Test for comparing two groups or Spearman Rank Order for correlations) using a value of p < 0.05 as a cutoff.
We would like to thank to the Vice Rectoría de Investigación y Desarrollo (VRID) of the Universidad de Santiago de Chile (USACH) for administrative support and travel assistance, Project Bicentenario PDA20 to MTR, CONICYT Fellowship to JMS, and CONICYT-COOP. INTERNACIONAL ACE-02 and INNOVA-CORFO 09MCSS-6698.
- Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus SE: Field's Virology. 2007, 5th ednGoogle Scholar
- Krossoy B, Nilsen F, Falk K, Endresen C, Nylund A: Phylogenetic analysis of infectious salmon anaemia virus isolates from Norway, Canada and Scotland. Dis Aquat Organ. 2001, 44: 1-6.PubMedView ArticleGoogle Scholar
- Blake S, Bouchard D, Keleher W, Opitz M, Nicholson BL: Genomic relationships of the North American isolate of infectious salmon anemia virus (ISAV) to the Norwegian strain of ISAV. Dis Aquat Organ. 1999, 35: 139-144.PubMedView ArticleGoogle Scholar
- Devold M, Karlsen M, Nylund A: Sequence analysis of the fusion protein gene from infectious salmon anemia virus isolates: evidence of recombination and reassortment. J Gen Virol. 2006, 87: 2031-2040. 10.1099/vir.0.81687-0.PubMedView ArticleGoogle Scholar
- Devold M, Falk K, Dale B, Krossoy B, Biering E, Aspehaug V, Nilsen F, Nylund A: Strain variation, based on the hemagglutinin gene, in Norwegian ISA virus isolates collected from 1987 to 2001: indications of recombination. Dis Aquat Organ. 2001, 47: 119-128.PubMedView ArticleGoogle Scholar
- Mjaaland S, Hungnes O, Teig A, Dannevig BH, Thorud K, Rimstad E: Polymorphism in the infectious salmon anemia virus hemagglutinin gene: importance and possible implications for evolution and ecology of infectious salmon anemia disease. Virology. 2002, 304: 379-391. 10.1006/viro.2002.1658.PubMedView ArticleGoogle Scholar
- Godoy MG, Aedo A, Kibenge MJ, Groman DB, Yason CV, Grothusen H, Lisperguer A, Calbucura M, Avendano F, Imilan M, et al: First detection, isolation and molecular characterization of infectious salmon anaemia virus associated with clinical disease in farmed Atlantic salmon (Salmo salar) in Chile. BMC Vet Res. 2008, 4: 28-10.1186/1746-6148-4-28.PubMedPubMed CentralView ArticleGoogle Scholar
- Christiansen DH, Ostergaard PS, Snow M, Dale OB, Falk K: A low-pathogenic variant of infectious salmon anemia virus (ISAV-HPR0) is highly prevalent and causes a non-clinical transient infection in farmed Atlantic salmon (Salmo salar L.) in the Faroe Islands. J Gen Virol. 2011, 92: 909-918. 10.1099/vir.0.027094-0.PubMedView ArticleGoogle Scholar
- McBeath AJ, Bain N, Snow M: Surveillance for infectious salmon anaemia virus HPR0 in marine Atlantic salmon farms across Scotland. Dis Aquat Organ. 2009, 87: 161-169.PubMedView ArticleGoogle Scholar
- McBeath A, Fourrier M, Munro E, Falk K, Snow M: Presence of a full-length highly polymorphic region (HPR) in the ISAV haemagglutinin-esterase does not affect the primary functions of receptor binding and esterase activity. Arch Virol. 2011, 156: 2285-2289. 10.1007/s00705-011-1106-9.PubMedView ArticleGoogle Scholar
- Plotkin JB, Kudla G: Synonymous but not the same: the causes and consequences of codon bias. Nat Rev Genet. 2011, 12: 32-42. 10.1038/nrg2899.PubMedPubMed CentralView ArticleGoogle Scholar
- dos Reis M, Savva R, Wernisch L: Solving the riddle of codon usage preferences: a test for translational selection. Nucleic Acids Res. 2004, 32: 5036-5044. 10.1093/nar/gkh834.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Ikemura T: Codon usage and tRNA content in unicellular and multicellular organisms. Mol Biol Evol. 1985, 2: 13-34.PubMedGoogle Scholar
- Sorensen MA, Kurland CG, Pedersen S: Codon usage determines translation rate in Escherichia coli. J Mol Biol. 1989, 207: 365-377. 10.1016/0022-2836(89)90260-X.PubMedView ArticleGoogle Scholar
- Moriyama EN, Powell JR: Codon usage bias and tRNA abundance in Drosophila. J Mol Evol. 1997, 45: 514-523. 10.1007/PL00006256.PubMedView ArticleGoogle Scholar
- Varenne S, Buc J, Lloubes R, Lazdunski C: Translation is a non-uniform process. Effect of tRNA availability on the rate of elongation of nascent polypeptide chains. J Mol Biol. 1984, 180: 549-576. 10.1016/0022-2836(84)90027-5.PubMedView ArticleGoogle Scholar
- Kanaya S, Yamada Y, Kinouchi M, Kudo Y, Ikemura T: Codon usage and tRNA genes in eukaryotes: correlation of codon usage diversity with translation efficiency and with CG-dinucleotide usage as assessed by multivariate analysis. J Mol Evol. 2001, 53: 290-298. 10.1007/s002390010219.PubMedView ArticleGoogle Scholar
- Gingold H, Pilpel Y: Determinants of translation efficiency and accuracy. Mol Syst Biol. 2011, 7: 481-PubMedPubMed CentralView ArticleGoogle Scholar
- Tats A, Tenson T, Remm M: Preferred and avoided codon pairs in three domains of life. BMC Genomics. 2008, 9: 463-10.1186/1471-2164-9-463.PubMedPubMed CentralView ArticleGoogle Scholar
- Gutman GA, Hatfield GW: Nonrandom utilization of codon pairs in Escherichia coli. Proc Natl Acad Sci U S A. 1989, 86: 3699-3703. 10.1073/pnas.86.10.3699.PubMedPubMed CentralView ArticleGoogle Scholar
- Tuller T, Carmi A, Vestsigian K, Navon S, Dorfan Y, Zaborske J, Pan T, Dahan O, Furman I, Pilpel Y: An evolutionarily conserved mechanism for controlling the efficiency of protein translation. Cell. 2010, 141: 344-354. 10.1016/j.cell.2010.03.031.PubMedView ArticleGoogle Scholar
- Tuller T: Codon bias, tRNA pools and horizontal gene transfer. Mob Genet Elements. 2011, 1: 75-77. 10.4161/mge.1.1.15400.PubMedPubMed CentralView ArticleGoogle Scholar
- Kudla G, Murray AW, Tollervey D, Plotkin JB: Coding-sequence determinants of gene expression in Escherichia coli. Science. 2009, 324: 255-258. 10.1126/science.1170160.PubMedPubMed CentralView ArticleGoogle Scholar
- Jia M, Li Y: The relationship among gene expression, folding free energy and codon usage bias in Escherichia coli. FEBS Lett. 2005, 579: 5333-5337. 10.1016/j.febslet.2005.08.059.PubMedView ArticleGoogle Scholar
- Bahir I, Fromer M, Prat Y, Linial M: Viral adaptation to host: a proteome-based analysis of codon usage and amino acid preferences. Mol Syst Biol. 2009, 5: 311-PubMedPubMed CentralView ArticleGoogle Scholar
- Gog JR: Afonso Edos S, Dalton RM, Leclercq I, Tiley L, Elton D, von Kirchbach JC, Naffakh N, Escriou N, Digard P: Codon conservation in the influenza A virus genome defines RNA packaging signals. Nucleic Acids Res. 2007, 35: 1897-1907. 10.1093/nar/gkm087.PubMedPubMed CentralView ArticleGoogle Scholar
- Marsh GA, Rabadan R, Levine AJ, Palese P: Highly conserved regions of influenza a virus polymerase gene segments are critical for efficient viral RNA packaging. J Virol. 2008, 82: 2295-2304. 10.1128/JVI.02267-07.PubMedPubMed CentralView ArticleGoogle Scholar
- Pinto RM, Aragones L, Costafreda MI, Ribes E, Bosch A: Codon usage and replicative strategies of hepatitis A virus. Virus Res. 2007, 127: 158-163. 10.1016/j.virusres.2007.04.010.PubMedView ArticleGoogle Scholar
- Aragones L, Guix S, Ribes E, Bosch A, Pinto RM: Fine-tuning translation kinetics selection as the driving force of codon usage bias in the hepatitis A virus capsid. PLoS Pathog. 2010, 6: e1000797-10.1371/journal.ppat.1000797.PubMedPubMed CentralView ArticleGoogle Scholar
- Coleman JR, Papamichail D, Skiena S, Futcher B, Wimmer E, Mueller S: Virus attenuation by genome-scale changes in codon pair bias. Science. 2008, 320: 1784-1787. 10.1126/science.1155761.PubMedPubMed CentralView ArticleGoogle Scholar
- Mueller S, Coleman JR, Papamichail D, Ward CB, Nimnual A, Futcher B, Skiena S, Wimmer E: Live attenuated influenza virus vaccines by computer-aided rational design. Nat Biotechnol. 2010, 28: 723-726. 10.1038/nbt.1636.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Nylund A, Plarre H, Karlsen M, Fridell F, Ottem KF, Bratland A, Saether PA: Transmission of infectious salmon anaemia virus (ISAV) in farmed populations of Atlantic salmon (Salmo salar). Arch Virol. 2007, 152: 151-179. 10.1007/s00705-006-0825-9.PubMedView ArticleGoogle Scholar
- Gray NK, Hentze MW: Regulation of protein synthesis by mRNA structure. Mol Biol Rep. 1994, 19: 195-200. 10.1007/BF00986961.PubMedView ArticleGoogle Scholar
- Mjaaland S, Markussen T, Sindre H, Kjoglum S, Dannevig BH, Larsen S, Grimholt U: Susceptibility and immune responses following experimental infection of MHC compatible Atlantic salmon (Salmo salar L.) with different infectious salmon anaemia virus isolates. Arch Virol. 2005, 150: 2195-2216. 10.1007/s00705-005-0588-8.PubMedView ArticleGoogle Scholar
- Buchan JR, Aucott LS, Stansfield I: tRNA properties help shape codon pair preferences in open reading frames. Nucleic Acids Res. 2006, 34: 1015-1027. 10.1093/nar/gkj488.PubMedPubMed CentralView ArticleGoogle Scholar
- Kibenge FS, Godoy MG, Wang Y, Kibenge MJ, Gherardelli V, Mansilla S, Lisperger A, Jarpa M, Larroquete G, Avendano F, et al: Infectious salmon anaemia virus (ISAV) isolated from the ISA disease outbreaks in Chile diverged from ISAV isolates from Norway around 1996 and was disseminated around 2005, based on surface glycoprotein gene sequences. Virol J. 2009, 6: 88-10.1186/1743-422X-6-88.PubMedPubMed CentralView ArticleGoogle Scholar
- Cottet L: Cortez-San Martin M, Tello M, Olivares E, Rivas-Aravena A, Vallejos E, Sandino AM, Spencer E: Bioinformatic analysis of the genome of infectious salmon anemia virus associated with outbreaks with high mortality in Chile. J Virol. 2010, 84: 11916-11928. 10.1128/JVI.01202-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Kibenge FS, Kibenge MJ, McKenna PK, Stothard P, Marshall R, Cusack RR, McGeachy S: Antigenic variation among isolates of infectious salmon anaemia virus correlates with genetic variation of the viral haemagglutinin gene. J Gen Virol. 2001, 82: 2869-2879.PubMedView ArticleGoogle Scholar
- Webster RG, Laver WG, Air GM, Schild GC: Molecular mechanisms of variation in influenza viruses. Nature. 1982, 296: 115-121. 10.1038/296115a0.PubMedView ArticleGoogle Scholar
- Bermudez-Santana C, Attolini CS, Kirsten T, Engelhardt J, Prohaska SJ, Steigele S, Stadler PF: Genomic organization of eukaryotic tRNAs. BMC Genomics. 2010, 11: 270-10.1186/1471-2164-11-270.PubMedPubMed CentralView ArticleGoogle Scholar
- Kramerov DA, Vassetzky NS: Origin and evolution of SINEs in eukaryotic genomes. Heredity (Edinb). 2011, 107: 487-495. 10.1038/hdy.2011.43.View ArticleGoogle Scholar
- Matveev V, Nishihara H, Okada N: Novel SINE families from salmons validate Parahucho (Salmonidae) as a distinct genus and give evidence that SINEs can incorporate LINE-related 3'-tails of other SINEs. Mol Biol Evol. 2007, 24: 1656-1666. 10.1093/molbev/msm083.PubMedView ArticleGoogle Scholar
- Morescalchi MA, Barucca M, Stingo V, Capriglione T: Polypteridae (Actinopterygii: Cladistia) and DANA-SINEs insertions. Mar Genomics. 2010, 3: 79-84. 10.1016/j.margen.2010.06.001.PubMedView ArticleGoogle Scholar
- Kawagoe-Takaki H, Nameki N, Kajikawa M, Okada N: Probing the secondary structure of salmon SmaI SINE RNA. Gene. 2006, 365: 67-73.PubMedView ArticleGoogle Scholar
- Izsvak Z, Ivics Z, Garcia-Estefania D, Fahrenkrug SC, Hackett PB: DANA elements: a family of composite, tRNA-derived short interspersed DNA elements associated with mutational activities in zebrafish (Danio rerio). Proc Natl Acad Sci U S A. 1996, 93: 1077-1081. 10.1073/pnas.93.3.1077.PubMedPubMed CentralView ArticleGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964.PubMedPubMed CentralView ArticleGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.PubMedView ArticleGoogle Scholar
This article is published under license to 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.