key: cord-0001039-imjiwslh authors: Goñi, Natalia; Iriarte, Andrés; Comas, Victoria; Soñora, Martín; Moreno, Pilar; Moratorio, Gonzalo; Musto, Héctor; Cristina, Juan title: Pandemic influenza A virus codon usage revisited: biases, adaptation and implications for vaccine strain development date: 2012-11-08 journal: Virol J DOI: 10.1186/1743-422x-9-263 sha: 9bd892682a0ebb3ab1361af5b375bbc11aa6d8b1 doc_id: 1039 cord_uid: imjiwslh BACKGROUND: Influenza A virus (IAV) is a member of the family Orthomyxoviridae and contains eight segments of a single-stranded RNA genome with negative polarity. The first influenza pandemic of this century was declared in April of 2009, with the emergence of a novel H1N1 IAV strain (H1N1pdm) in Mexico and USA. Understanding the extent and causes of biases in codon usage is essential to the understanding of viral evolution. A comprehensive study to investigate the effect of selection pressure imposed by the human host on the codon usage of an emerging, pandemic IAV strain and the trends in viral codon usage involved over the pandemic time period is much needed. RESULTS: We performed a comprehensive codon usage analysis of 310 IAV strains from the pandemic of 2009. Highly biased codon usage for Ala, Arg, Pro, Thr and Ser were found. Codon usage is strongly influenced by underlying biases in base composition. When correspondence analysis (COA) on relative synonymous codon usage (RSCU) is applied, the distribution of IAV ORFs in the plane defined by the first two major dimensional factors showed that different strains are located at different places, suggesting that IAV codon usage also reflects an evolutionary process. CONCLUSIONS: A general association between codon usage bias, base composition and poor adaptation of the virus to the respective host tRNA pool, suggests that mutational pressure is the main force shaping H1N1 pdm IAV codon usage. A dynamic process is observed in the variation of codon usage of the strains enrolled in these studies. These results suggest a balance of mutational bias and natural selection, which allow the virus to explore and re-adapt its codon usage to different environments. Recoding of IAV taking into account codon bias, base composition and adaptation to host tRNA may provide important clues to develop new and appropriate vaccines. Influenza A virus (IAV) is a member of the family Orthomyxoviridae and contains eight segments of a singlestranded RNA genome with negative polarity [1] . IAV is one of the most important infectious diseases in humans [2] . Unlike most pathogens where exposure leads to lasting immunity in the host, IAV presents a moving antigenic target [3] , evading specific immunity triggered by previous infections. This process, called antigenic drift, is the result of the selective fixation of mutations in the gene encoding the hemagglutinin (HA) protein and to a lesser extent in the neuraminidase (NA) protein [4] . Variants that best escape the host immune response are thought to have a significant reproductive advantage [5] . Another process, called reassortment, is also considered a major force in the evolution of IAV [4] . It occurs when the virus acquires an HA and/or NA of a different IAV subtype (via reassortation) of one or more gene segments. This process has been in the basis of the devastating influenza pandemics that occurred several times in the last century [6] . The first influenza pandemic of this century was declared in April of 2009, with the emergence of a novel H1N1 IAV strain (H1N1pdm) in Mexico and USA [7, 8] . By November of 2009, the virus was detected in about 207 countries, infecting more than 620,000 individuals worldwide and accounting for more than 7,800 deaths [7] . This strain was a multiple reassortant with genes derived from viruses that originally circulated in the swine, avian and human populations [9] . It has been observed that IAV is subjected to host immune selection pressure and undergoes rapid evolution, especially when the virus crosses the host species barrier [10] . The replication cycle of IAV depends on host machinery and the virus utilizes host cellular components for its protein synthesis. Therefore, the interplay of codon usage of virus and host could affect viral replication. For these reasons, a detailed understanding of IAV evolution and host adaptation is crucial. Due to the degeneracy of the genetic code, most amino acids are coded by more than one codon. Synonymous triplets are not used randomly. In several organisms, natural selection and mutational input seem to bias codon use toward a certain subset of codons [11] . Two major models have been proposed to explain codon usage: the translational selection and the mutational models [12] . Codon usage bias related to translation efficiency (at two different levels: speed and accuracy) seems to be linked to local cognate isoacceptors tRNAs abundances, which in turn determine the major codon preferences [13] . On the other hand, discrepancies on codon usage could be due to genome compositional constraints and mutational biases [14] . Nevertheless, these two models cannot be considered as mutually exclusive. Although previous studies have been performed on the general codon usage of IAV [2, 12, 15, 16] , a deep and comprehensive study to investigate the effect of selection pressure imposed by the human host on the codon usage of an emerging, pandemic IAV strain and the trends in viral codon usage involved over the pandemic time period is much needed. In order to gain insight into these matters, we performed a comprehensive codon usage analysis of 310 H1N1pdm IAV strains, isolated from April to September of 2009, for which the complete genome sequences are available. In order to study the extent of codon usage bias in H1N1pdm IAV strains in relation to seasonal H1N1 and H3N2 as well as human and swine host cells, the relative synonymous codon usage (RSCU) [14] values for each codon were calculated for the 310 H1N1pdm strains enrolled in these studies and compared with seasonal IAV strains and host organisms. The results of these studies are shown in Table 1 . All codons containing the dinucleotide CpG were underrepresented in all IAV viruses. Important differences were found between human and swine hosts and IAV strains. Particularly, high biased frequencies (Δ RSCU ≥ 0.30) were found for Leu, Ile, Val, Ser, Pro, Thr, Ala, His, Gln, Glu, Arg and Gly. Interestingly, the huge majority of preferred codons in the viruses are A-ended. In the case of Arg, there is a strong bias towards an increase in AGA and AGG, while the CGN codons are depleted (see Table 1 ). To observe if H1N1pdm IAV strain sequences display similar codon usage biases, the effective number of codons (ENC) [17] values were calculated for the 310 strains enrolled in this study (mean of 52.51 ± 0.05). ENC varies from 20 to 61, where the larger the extent of codon bias in a gene, the smaller the ENC value. Thus, a value of 52.5 strongly suggests that the overall codon usage among these strains is only slightly biased. Since codon usage by its very nature is multivariate, it is necessary to analyze the data using multivariate statistical techniques, like correspondence analysis (COA) [18] . The correlation between the position on the first dimensional factor generated by this analysis on RSCU (20.7% of the total variability) for each strain and the respective G + C content at synonymous variable third position (GC 3 s) values was significant (r = −0.47, p < 0.0001). Interestingly, this dimensional factor also significantly correlated with A content at synonymous variable third position (A 3 s, r = 0.68, p < 0.0001) and G content at the same position (G 3 s, r = −0.71, p < 0.0001) ( Figure 1 ). This means that the major factor shaping codon usage among these strains is an opposite trend between purines at third codon positions. Furthermore, this result is mainly due to the frequencies of the codons CGA (Arg) on one side of the distribution and GCG (Ala) and CGG (Arg) at the other side (see Additional file 1: Table S1 ). In other words, the differential usage of three low frequent codons (RSCU ≤ 0.63) is among the major factor shaping codon usage among these strains. It has been suggested that dinucleotide biases can affect codon bias [19] . To study the possible effect of dinucleotide composition on codon usage of the H1N1pdm IAV strains, the relative abundances of the 16 dinucleotides in the ORFs of the 310 strains enrolled in these studies were established. The results of these analyses are shown in Table 2 . As it can be seen in the table, the occurrences of dinucleotides are not randomly distributed and no dinucleotides were present at the expected frequencies ( Table 2 ). The relative abundance of CpG showed a strong deviation from the "normal range" (mean ± S.D. = 0.319 ± 0.0020) and were markedly underrepresented. Interestingly, when the second dimensional factor (11.1% of the total variability) was analyzed, we found that the position of each strain significantly correlated (r = 0.64, p < 0.0001) with the respective usage of the dinucleotide CpG. Besides, although the global usage of this dinucleotide is very low, we found that the correlation is due to the differential usage of CGU (Arg) and CCG (Pro) codons, since these triplets display the most extreme values on the second dimensional factor (see Additional file 1: Table S1 ). Importantly, we also found that the third and the fourth dimensional factors of COA (8.7% and 5.5% of the total variability, respectively), are again mainly linked to the low usage of codons containing the dinucleotide CpG, mainly at the positions 2 and 3. Moreover, among the 16 dinucleotides, 15 are highly correlated with the first dimensional factor value in COA (Table 2 ). These observations indicate that the composition of dinucleotides also plays a crucial role in the variation found in synonymous codon usage among H1N1pdm IAV ORFs. To study the possibility of codon usage variation in the H1N1pdm IAV genomes enrolled in this study, the distribution of the 310 strains in the plane defined by the first two axes of COA was established. The results of these studies are shown in Figure 2 . Interestingly, the distribution of the H1N1pdm IAV strains in the plane defined by the first two major axes showed that the principal dimensional factor splits the strains at least three major groups: two of them discriminated by the first dimensional factor, while the third is revealed by the extreme low values on the second dimensional factor ( Figure 2 ). As the translation process represents a key step in the viral infection cycle, it is important to explore the strategies employed by the virus to harness the translation machinery of the cell host. Since variation at the third codon position makes possible the wobble interaction between that base and the first one of the anticodon [20] , we wanted to gain further insight into the adaptation of H1N1pdm IAV strains to the respective host tRNA pool context. For this reason, the codon usage of virus (H1N1pdm IAV) was plotted against the codon usage of host (human cells) and the nucleotide that occupy the first anticodon position (wobble position) of the corresponding codon was identified. The results of these studies are shown in Figure 3 . As it can be seen in the figure, codon usage of virus and host is uncorrelated. The viral preference toward AT rich genomes and the T-headed anticodons is clear (Figure 3 ). This is also in agreement with the consequence of a differential usage of A3 s and G3 s (see also Figure 1 ). Comparison of these findings with the compilation of tRNAs species in the human genome [21] reveals that the virus highly preferred T-headed anticodons are not particularly adapted to the host transfer tRNA pool (Table 3) . Therefore, there is no obvious correlation between the number of human host isoacceptor tRNAs and codon usage of the IAV enrolled in these studies. As IAV relies on the host cell's machinery for its replication, codon usage bias could play a role in its adaptation to the host. The results of these studies revealed that codon usage in human IAV, including H1N1pdm, do not have the average codon usage pattern of their host's genes (see Table 1 ), in agreement with previous reports [12, 16] . Comparisons to previous results reported for other IAV such H5N1 (mean ENC = 50.91) [16, 22] [26] or Theilovirus (mean ENC = 51.08) [26] , revealed that the ENC values found in this study for H1N1pdm IAV strains (mean ENC value of 52.5) are roughly similar to these previous findings, indicating that the overall extent of codon usage in these viruses are only slightly biased. We have found a general link between codon usage bias and base composition, which is shown by the significant correlation of the position of each virus on the first dimensional factor of COA vs. the corresponding GC 3 s, together with the opposite trends in relation to purines at third codon position ( Figure 1A and B). Taken together, our results indicate that the mutational bias is a very important trend in the evolution of H1N1pdm IAV genomes. However, this does not per se discards a role of other natural selection mechanisms acting in the IAV strains enrolled in these studies. We have also found that CpG containing codons are sharply suppressed (see Table 1 ). This CpG deficiency was proposed to be related to the immunostimulatory properties of unmethylated CpG, which are recognized by the innate immune system of the host as a pathogen signature [24, 27] . This is triggered by the intracellular Pattern Recognition Receptor (PRR) Tool-like 9 (TLR9), which activates several immune response pathways [28] . It seems reasonable to suggest that exists among vertebrates a TLR9-like mechanism acting at the RNA level [29] . Interestingly, previous studies have shown that IAV strains originated from an avian reservoir and infecting human hosts since 1918 has been selected under strong pressure to reduce the frequency of CpG in its genome [30] . Marked CpG deficiency has been observed in several other RNA viruses [24, [31] [32] [33] [34] [35] , including H1N1pdm IAV [12, 30] . Then, escaping from the host antiviral response may act as another selective pressure contributing to codon usage in H1N1pdm IAV strains [36] . The distribution of the 310 H1N1 pdm IAV ORF's in the plane defined by the first two axes of COA shows the presence of at least three clusters of strains ( Figure 2 ). Since species with a close genetic relationship always present a similar codon usage pattern [37] (see also Table 1 ), the results of these studies suggests that a dynamic process occurred in the H1N1pdm strains enrolled in these studies. This is reflected in the variation of codon usage observed among them (see Figure 2) . These results suggest a balance of mutational bias and natural selection to shape codon usage in these strains, which allow the virus to explore and re-adapt its codon usage to different environments in a short period of time. From the classical point of view, the preferred codons are recognized by the most abundant isoacceptors tRNAs, which implies the action of natural selection [38] . The results shown in Table 3 strongly suggest that this is not the case for H1N1pdm IAV strains. In other words, codon usage of these viruses does not seem to be adapted to the tRNA pool of the human cells but probably reflects the influence of mutational biases. Interestingly, this has been observed for some other RNA viruses, like HIV [39] . Understanding the mechanisms used by IAV to properly express its genes could suggest a novel point of intervention and drug targets. Reduced translation efficiency, particularly of structural genes that are needed for the formation of new particles, could affect viral success [40] . The results of this work suggest that synthetic attenuated virus engineering (SAVE) could play a role in creating new vaccines for IAV. By deoptimization of codon usage (replacing wild-type codons with codons and codon combinations whose sequences impair replication and/or expression), it might be possible to attenuate a virus [41] . Moreover, as the codon changes do not alter the protein sequence, the antigenicity should not differ from the wild-type virus. Besides, codon changes tend to have individually small fitness effects, so many nucleotide changes will be required to restore wild-type fitness, itself requiring 100 s or more generations [42] [43] [44] [45] . This "death by a thousand cuts" strategy may provide an alternative method of attenuation [46] . Interestingly, it has been show that replacement of natural codons with synonymous triplets with increased frequencies of CpG gives rise to inactivation of Poliovirus infectivity [47] . Very recent studies revealed that this strategy can be applied to IAV [48] . Gly GGA UCC(9), GCC(15), CCC (7), ACC(0) 31 His CAU AUG(0), GUG (11) 11 Ile AUA UAU(5), AAU (14), GAU (8) 27 Leu CUA UAG(3), AAG(12), CAG(10), CAA (7), UAA (7), GAG(0) 39 Pro CCA UGG (7), AGG (10) Owing to known genome sequences, modern strategies of DNA synthesis have made it possible to recreate in principle all known viruses independent of natural templates [48] . Recoding of IAV to develop new vaccine candidates taking into account codon bias, base composition and adaptation to host tRNA by gene synthesis may provide important clues to elucidate virulence factors, identify targets for future drug intervention, and to develop new and appropriate vaccines [49] . Sequences from H1N1pdm IAV strains, isolated from April to December of 2009, were obtained from The Influenza Virus Resource at the National Center for Biotechnological Information [50] . The data set comprised the complete genome sequences (eight segments) of 310 strains. For each strain the ORFs were concatenated (PB2 + PB1 + PA + HA + NP + NA + MP + NS) and aligned using the MUSCLE program [51] . The alignment is available upon request. Codon usage, base dinucleotide composition, G + C at synonymous variable third position codons (GC 3 s), the relative synonymous codon usage (RSCU) [14] and the effective number of codons (ENC) [17] were calculated using the program CodonW (written by John Peden and available at http://sourceforge.net/projects/codonw/) as implemented in the Mobile server (http://mobyle. pasteur.fr). Codon usage data of influenza viral hosts, human (Homo sapiens) and domestic swine (Sus scrofa) were obtained from the codon usage database (available at: http://www.kazusa.or.jp/codon) [52] . The frequencies of tRNAs in human cells were retrieved from the GtRNAdb database [21] . COA is an ordination technique that identifies the major trends in the variation of the data and distributes genes along continuous axes in accordance with these trends. COA creates a series of orthogonal axes to identify trends that explain the data variation, with each subsequent dimensional factor explaining a decreasing amount of the variation [18] . Each ORF is represented as a 59-dimensional and each dimension is related to the RSCU value of each triplet (excluding AUG, UGG and stop codons). This was done using the CodonW program. Correlation analysis was carried out using Spearman's rank correlation analysis method [53] . Orthomyxovirus: replication, transcription, and polyadenylation Comparative study of the hemagglutinin and neuraminidase genes of Influenza A virus H3N2, H9N2 and H5N1 subtypes using bioinformatics techniques Long intervals of stasis punctuated by bursts of positive selection in the seasonal evolution of influenza A virus Realities and enigmas of human viral influenza: pathogenesis, epidemiology and control Influenza virus: a master of metamorphosis Ecological and immunological determinants of influenza evolution Influenza-like illness in the United States and Mexico Update: infections with a swine-origin influenza A (H1N1) virus -United States and other countries Origins and evolutionary genomics of the 2009 swine-origin H1N1 Influenza A epidemic Evolution of influenza A virus nucleocapsid genes: implications for the origins of H1N1 human and classical swine viruses Synonymous codon usage in Escherichia coli: selection for translational accuracy Codon usage bias and the evolution of Influenza A viruses. Codon usage biases of Influenza virus Correlation between the abundance of yeast transfer RNAs and the occurrence of the respective codons in protein genes. Differences in synonymous codon choice patterns of yeast and Escherichia coli with reference to the abundance of isoaccepting transfer RNAs An evolutionary perspective on synonymous codon usage in unicellular organisms Natural selection for nucleotide usage at synonymous and non-synonymous sites in the influenza A genes Analysis of synonymous codon usage in H5N1 virus and other influenza A viruses An evaluation of measures of synonymous codon usage bias Theory and applications of correspondence analysis. London: Academic Analysis of synonymous codon usage in classical swine fever virus Codon-anticodon pairing -Wobble hypothesis GtRNAdb: a database of transfer RNA genes detected in genomic sequence Analysis of synonymous codon usage bias in 09H1N1 Cytosine deamination and selection of CpG suppressed clones are the two major independent biological forces that shape codon usage bias in Coronaviruses Mutation pressures shapes codon usage in the GC-rich genome of foot-and-mouth disease virus Analysis of synonymous codon usage in the UL24 gene of duck enteritis virus Compare the differences of synonymous codon usage between the two species within cardiovirus Evolutionary basis of codon usage and nucleotide composition bias in vertebrate DNA viruses Clinical application of CpG-, non-CpG, and antisense oligodeoxynucleotides as immunomodulators Virus-host coevolution: common patterns of nucleotide motif usage in Flaviviridae and their hosts Patterns of evolution and host gene mimicry in influenza and other RNA viruses Comparison of avian and human influenza A viruses reveals a mutational bias on the viral genomes Mononucleotide and dinucleotide frequencies and codon usage in poliovirus RNA Why is CpG suppressed in the genomes of virtually all small eukaryotic viruses but not in those of large eukaryotic viruses Analysis of human P[4]G2 rotavirus strains isolated in Brazil reveals codon usage bias and strong compositional constraints A detailed comparative analysis on the overall codon usage patterns in hepatitis A virus Genome rhetoric and the emergence of compositional bias 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 Codon usage and tRNA content in unicellular and multicellular organisms HIV-1 modulates the tRNA pool to improve translation efficiency Quantitative effect of suboptimal codon usage on translational efficiency of mRNA encoding HIV-1 gag in intact T cells Virus attenuation by genome-scale changes in codon pair bias Modulation of poliovirus replicative fitness in HeLa cells by deoptimization of synonymous codon usage in the capsid region Reduction of the rate of poliovirus protein synthesis through large-scale codon deoptimization causes attenuation of viral virulence by lowering specific infectivity Virus attenuation by genome-scale changes in codon pair bias Slow fitness recovery in a codon-modified viral genome Live attenuated Influenza virus vaccines by computer-aided rational design Genetic inactivation of poliovirus infectivity by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous capsid region codons Synthetic poliovirus and other designer viruses: what have we learned from them? Synthetic viruses: a new opportunity to understand and prevent viral disease The Influenza Virus Resource at the National Center for Biotechnology Information MUSCLE: a multiple sequence alignment method with reduced time and space complexity Codon usage tabulated from international DNA sequence databases: status for the year 2000 Office for Research Development and Education, version 1.1.23-r7 Pandemic influenza A virus codon usage revisited: biases, adaptation and implications for vaccine strain development We acknowledge support by International Atomic Energy Agency, through Research Contract No. 15792 and Comisión Sectorial de Investigación Científica (CSIC), Universidad de la República, Uruguay, through I + D Project "Variabilidad genetica y evolución viral de virus Influenza A en Uruguay". Authors acknowledge support by Agencia Nacional de Investigación e Innovación (ANII) through project PE_ALI_2009_1_1603 and PEDECIBA, Uruguay. We also thank the support of Fondo Clemente Estable, 2007_ 722 to HM. Additional file 1: Table S1 . Each codon included in the correspondence analysis is represented by a row. Factor 1 and 2 columns contain the coordinate of the codon on the respective generated axis. The authors declare that they do not have competing interests.Authors' contributions JC conceived of the study, and participated in its design and coordination. NG and AI have made substantial contributions to the design of the study, acquisition of data and analysis. VC, MS, GM, GC, have been involved in revising the manuscript critically for important intellectual content. JC wrote the paper. HM helped to draft the manuscript and made substantial and fundamental contributions to the interpretation and discussion of the results found in this work. All authors read and approved the final manuscript.