J Gen Virol Tips for Better Browsing
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Gen Virol 89 (2008), 1309-1313; DOI 10.1099/vir.0.83662-0

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.
Agricola
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.

Short Communication

Quasispecies dynamics and fixation of a synonymous mutation in hantavirus transmission

Tarja Sironen1, Eva R. Kallio2,3, Antti Vaheri1, Åke Lundkvist4 and Alexander Plyusnin1,4

1 Department of Virology, Haartman Institute, University of Helsinki, Finland
2 Finnish Forest Research Institute, Vantaa, Finland
3 Department of Biological and Environmental Science, University of Jyväskylä, Finland
4 Swedish Institute for Infectious Disease Control and MTC, Karolinska Institutet, Stockholm, Sweden

Correspondence
Tarja Sironen
Tarja.Sironen{at}helsinki.fi


   ABSTRACT
TOP
ABSTRACT
MAIN TEXT
REFERENCES
 
RNA-dependent RNA polymerases, the key enzymes in replication of RNA viruses, have a low fidelity; thus, these viruses replicate as a swarm of mutants termed viral quasispecies. Constant generation of new mutations allows RNA viruses to adapt swiftly to a novel environment through selection of both pre-existing and de novo-generated genetic variants. Here, quasispecies dynamics were studied in vivo in controlled hantavirus transmission from experimentally infected to naïve rodents through infested cage bedding. An elementary step of virus microevolution was apparent, as one synonymous mutation (A759G) repeatedly became fixed in the viral RNA quasispecies populations in the recipient animals.


   MAIN TEXT
TOP
ABSTRACT
MAIN TEXT
REFERENCES
 
Puumala virus (PUUV) is a hantavirus (family Bunyaviridae); it is carried by bank voles (Myodes glareolus) (Nichol et al., 2005Down) and causes a lifelong persistent and asymptomatic infection in these rodents (Plyusnin et al., 1996bDown). PUUV is pathogenic in humans, causing nephropathia epidemica (NE), a mild form of haemorrhagic fever with renal syndrome (Vapalahti et al., 2003Down). Hantavirus transmission occurs horizontally via direct contacts or through contaminated environment or aerosols in animals (Meyer & Schmaljohn, 2000Down). Humans acquire the infection by inhalation of contaminated rodent excreta.

For years, PUUV has been a model hantavirus in research on viral genetics and evolution. To our benefit, this virus can be studied both in cell culture and in the natural rodent host. The hantavirus genome consists of three segments of negative-polarity RNA (S, M and L) (Nichol et al., 2005Down) and, similarly to other RNA viruses, populations of PUUV in individual hosts are represented by quasispecies (Lundkvist et al., 1997Down; Plyusnin et al., 1995Down). Hantaviruses have adapted to their natural hosts during millions of years of co-evolution (Plyusnin & Morzunov, 2001Down) and thus the hantavirus quasispecies populations in the rodent hosts usually have a high fitness. Changes in the viral environment may, however, lead to rapid adaptive changes. For example, the adaptation of wild-type (wt) PUUV (passaged in colonized bank voles) to cultured primate Vero E6 cells is accompanied by increased genetic heterogeneity of the viral population and fixation of nucleotide changes in the S and L segments (Lundkvist et al., 1997Down; Nemirov et al., 2003Down). Even with minimal genetic differences between them, the two PUUV variants were clearly distinct phenotypically: the cell culture-adapted variant was non-infectious to bank voles. Interestingly, it is also non-pathogenic to cynomolgus macaques, which are used as a non-human primate model of NE (Groen et al., 1995Down). In contrast, the wt PUUV strain, passaged only in bank voles, caused in macaques a disease that mirrors NE of humans (Klingström et al., 2002Down).

Our aim was to investigate the quasispecies dynamics in bank voles when PUUV is transmitted from one animal to another in a controlled experiment. In vivo experiments on PUUV transmission have been described previously (Kallio et al., 2006Down) and are presented in Fig. 1Down. Briefly, colonized bank voles were inoculated subcutaneously with 100 bank vole ID50 of PUUV, strain Kazan-wt (Gavrilovskaya et al., 1983Down). Ten days later, they were placed into individual experimental cages for 7 days. The donor voles (labelled D) were removed and recipient voles (labelled R) were placed into the cages, where they were exposed to the excrement-contaminated bedding for 3 days. All precautions were taken to ensure that no aerosol transmission between the cages was possible (Kallio et al., 2006Down). After exposure, the recipient voles were moved to individual management cages for 15 days to develop the infection, after which the animals were sacrificed and samples collected. Five voles were subsequently exposed to each donor's bedding. PUUV infection was confirmed by using nested RT-PCR to detect PUUV S segment RNA in lung tissue. Bank voles positive for PUUV RNA in the RT-PCR test were used in this study.


Figure 1
View larger version (21K):
[in this window]
[in a new window]

 
Fig. 1. Scheme of PUUV infections. All voles were held in individual cages throughout the experiment; the cages were placed in HEPA-filtered isolators. Donor voles (D) were inoculated with PUUV and, after the infection had developed for 10 days, the donors were placed into the experimental cages (E). The donors excreted the virus into the cage bedding for 7 days, after which they were removed and sacrificed on day 17. On this day, the first recipient voles (R) were placed into the experimental cage, where they were exposed to the contaminated bedding; PUUV was transmitted into these voles. After 3 days exposure, the recipient voles were moved into the management (M) cages, in which the infection was developed for 15 days; after this, the recipient animals were sacrificed. Five voles were consequently exposed to each donor's bedding. Donors 1–5 (D1–D5) and all of their recipients were held in isolator 1, and donors D6–D10 and their recipients in isolator 2.

 
The region selected for analysis (nt 631–1630) included both the coding region (CR) and the 3' non-coding region of the S segment RNA. The selected part of the CR encoded both a well-conserved C terminus and a variable part of the N protein (Kaukinen et al., 2005Down). All three regions seem to evolve under different selective pressures. Altogether, approximately 60 kb was sequenced for analysis of the S segment quasispecies from two donor and four recipient bank voles. RNA was extracted from the lung tissue and the entire S segment was amplified by RT-PCR as described previously (Plyusnin et al., 1994Down). PCR products were cloned by using the pGEM-T cloning system (Promega). Ten cDNA clones from each of the two donors and from four recipients were sequenced (nt 631–1630), using primers PuuSF572 (5'-TATATGTATCCATGCCTACTGCC-3') and PuuSR1689 (5'-GATAGCTCAGTTTCACATTTTTG-3'). Sequencing with these two primers resulted in twofold coverage of nt 850–1350. Thus, regions 631–850 and 1350–1630 were initially sequenced in one reading, but if any ambiguous nucleotides were detected, the sequencing reaction was repeated.

Total mutation frequency in the donor sequences was 1.6x10–3 [31 mutations in 20 kb sequenced, mean±SD=(1.6±1.273)x10–3]; in the recipient sequences, it was almost twice as high [2.6x10–3, mean±SD=(2.63±1.234)x10–3] (Table 1Down), and the difference is statistically significant (Student's t-test, P=0.004). Nucleotide substitutions were distributed evenly across this region (data not shown) and also evenly across the different codon positions (Table 1Down), suggesting no stabilizing selection at the protein level for individual de novo-generated mutations. Mutations were predominantly transitions (92 %), and most mutations (78 %) were either A->G or U->C substitutions.


View this table:
[in this window]
[in a new window]

 
Table 1. Summary of mutations detected in the quasispecies swarm

The level of PUUV diversity seen in this study is seven to 15 times higher than the level of misincorporations produced by AmpliTaq polymerase (Lundkvist et al., 1997Down; Plyusnin et al., 1996aDown).

 
One synonymous mutation (A759G) appeared to become fixed in the viral quasispecies population during a single virus-transmission event. First, donor D1 and its two recipients, D1-R3 and D1-R4, were analysed and the mutation was seen in all 20 cDNA clones from the recipient voles, but in none of the D1-originating cDNA clones. Next, donor D2 and its two recipients, D2-R3 and D2-R4, were analysed. Again, the A759G mutation was observed in all of the recipient cDNA clones. Surprisingly, it was also seen in 60 % of donor D2 clones. An additional 30 clones of donor D1 were sequenced to see whether the mutation was pre-existing in its quasispecies population, but the mutation was not observed. Five additional donors were examined for the presence of this mutation by using direct sequencing of the partial S PCR product that was amplified by nested RT-PCR essentially as described previously (Plyusnin et al., 1997Down), but using the primers PUUSF688 [5'-ATGAGTCC(A/T)GTIATGGG(A/G)GTIATT-3'] and PUUSR1106 [5'-GCTGTICCTACAGT(C/T)TTIGATGCCAT-3'] in the second round. The mutation was found in only one of the donors. All PUUV-positive recipients of the five donors were also examined for the presence of this mutation (Fig. 2Down). In some recipient voles, the sequence with the A residue remaining at position 759 could also be seen in the sequencing file as a minor peak (data not shown). This suggested that, at least for one transmission event, the original nucleotide was preserved in the quasispecies swarm. Four additional donors and their five PUUV-positive recipients from the second trial of PUUV transmission (Kallio et al., 2006Down) were also studied, but the mutation A759G was not found in any of them.


Figure 2
View larger version (14K):
[in this window]
[in a new window]

 
Fig. 2. Occurrence of the mutation A759G. White voles carry the wild-type S sequence; black voles carry the mutation. In donor D2 (grey vole), 60 % of the cDNA clones are mutated. A dash represents either a PUUV-negative vole or a sample that was no longer available. Donors D1 and D2 and their recipient voles in groups R3 and R4 were used for quasispecies analysis via generation of cDNA clones. Partial S sequences from other voles were recovered by direct sequencing.

 
The wt S segment sequence (without the A759G mutation) was found in the quasispecies population, which has an overall mutation frequency of 1.4x10–3 at the nucleotide level. Notably, the mutation was found on a background of substantially higher genetic variability: the mutation frequency was 2.6x10–3 at the nucleotide level (t-test, P=0.002). The mutation is located at the third codon position of phenylalanine 239 and it changes the triplet in mRNA from UUU to UUC. Codon-usage analysis (Nakamura et al., 2000Down) showed that the UUC codon is more frequent in bank voles and humans, but in PUUV it is the rarer codon.

Quasineutral mutations have been detected in hantavirus genomes (Plyusnin et al., 1995Down, 1996aDown; Feuer et al., 1999Down), but in only one study (Lundkvist et al., 1997Down) has it been possible to follow whether they become dominant. In our study, we observed two alleles of the PUUV S gene in the population used to infect the donor voles. In the first trial (donors D1–D10) of the transmission experiment (Kallio et al., 2006Down), the mutation A759G was detected in all PCR-positive recipients (Fig. 2Up). In the second trial (donors D16–D21), however, it was not detected. The reason for this might be the somewhat lower transmission efficiency in the second trial (Kallio et al., 2006Down), which could result from as-yet-unknown factors, such as the humidity. Alternatively, it could have occurred because of sampling from the genetically heterogeneous virus stock used for inoculation of donor voles. When the variant with the A759G mutation nevertheless appeared, it seemed to have an advantage in natural transmission, as it was found in all cDNA clones from the recipients. This is, to our knowledge, the first time that it has been possible to follow the mutations appearing in a hantavirus quasispecies swarm in vivo, and to detect fixation of a synonymous mutation in transmission from one animal to another.

Viral RNA-dependent RNA polymerases are notorious for their high mutation rates (Domingo & Holland, 1997Down). New RNA virus mutants are continuously generated and assessed by positive (diversifying) and negative (stabilizing) selective forces (Domingo et al., 2001Down), which act at the population level rather than on individual variants (Vignuzzi et al., 2006Down). The evolution of PUUV is driven mainly by genetic drift through accumulation of point mutations or deletions/insertions (Plyusnin et al., 1995Down; Sironen et al., 2001Down). In theory, population bottlenecks – common events in virus transmission – should also affect the evolution of PUUV externally. A bottleneck reduces both the effective size and the genetic diversity of a viral population, as only a few particles are sampled randomly to found the novel population of the virus (Manrubia et al., 2005Down).

As mutation A759G is silent, the positive advantage brought by the mutation cannot arise on the level of protein structure. Changes in the predicted secondary structure of either the viral RNA or the mRNA appeared minor (not shown) and thus are unlikely to be the sole, if any, reason for the advantageous phenotypic effect of the mutation. Another possibility might be a general enhancement of viral replication, a phenomenon described recently for vesicular stomatitis virus and human immunodeficiency virus type 1 (Hamano et al., 2007Down; Novella et al., 2004Down).

In our experiments, the mutation A759G appears to be advantageous and selected for, as it becomes fixed in the PUUV master sequence in the recipient voles. As it was not feasible to sequence whole genomes in the frame of this project, the possibility remains that this substitution is a marker for an accompanying amino acid substitution located in a region not yet analysed. Further experiments with additional passages in animals would also reveal whether this mutation becomes fixed permanently in the virus genome. An important question is whether the mutation A759G is generated de novo in the process of viral adaptation or whether it is preserved, as a minority, in the quasispecies swarms within the donor voles. Notably, this particular mutation was seen earlier in two of 25 cDNA clones originating from the wt PUUV adapting to cell culture (Lundkvist et al., 1997Down). This suggested that the mutation can be maintained permanently, even as a minor variant, in the quasispecies swarm. This view is also supported by the finding of this mutation in two of the donors (D2 and D5) (Fig. 2Up). The adaptation to a completely different host (from voles to cell culture) was accomplished via only a few changes in the PUUV genome (Lundkvist et al., 1997Down; Nemirov et al., 2003Down). These small alterations, mainly on the RNA level, nevertheless changed the virus phenotype profoundly, and the virus was no longer able to infect the original host, the bank vole. Similarly, the silent mutation A759G may provide the virus with a strong advantage in transmission. Interestingly, both A and G nucleotides occur in the consensus S segment sequences of PUUV strains (all sequences from the Russian lineage, Kazan strain included, have A), suggesting that neither purine at position 759 is deleterious to the virus.

Genetic drift is the main mechanism of PUUV evolution (Plyusnin et al., 1995Down; Sironen et al., 2001Down). This study, for the first time, confirms the result in vivo. Furthermore, these data add to the increasing evidence (Hamano et al., 2007Down; Novella et al., 2004Down) that silent mutations in RNA viruses may be far from neutral and, on the contrary, can contribute significantly to adaptation of these agents to the ever-changing environment.


   ACKNOWLEDGEMENTS
 
Sonja Ahlfors is thanked for her assistance in sequencing. This work was supported by grants from the Academy of Finland, Helsinki Biomedical Graduate School, Sigrid Jusélius Foundation, the Swedish Research Council (project no. 12177), Oscar Öflund Foundation, Finnish Cultural Foundation and Nordforsk.


   REFERENCES
TOP
ABSTRACT
MAIN TEXT
REFERENCES
 
Domingo, E. & Holland, J. J. (1997). RNA virus mutations and fitness for survival. Annu Rev Microbiol 51, 151–178.[CrossRef][Medline]

Domingo, E., Biebricher, C. K., Eigen, M. & Holland, J. J. (2001). Quasispecies and RNA Virus Evolution: Principles and Consequences. Georgetown, TX: Landes Bioscience.

Feuer, R., Boone, J. D., Netski, D., Morzunov, S. P. & St Jeor, S. C. (1999). Temporal and spatial analysis of Sin Nombre virus quasispecies in naturally infected rodents. J Virol 73, 9544–9554.[Abstract/Free Full Text]

Gavrilovskaya, I. N., Chumakov, M. P., Apekina, N. S., Ryltseva, E. V., Martiyanova, L. I., Gorbachkova, E. A., Bernshtein, A. D., Zakharova, M. A. & Boiko, V. A. (1983). Adaptation to laboratory and wild animals of the haemorrhagic fever with renal syndrome virus present in the foci of European U.S.S.R. Brief report. Arch Virol 77, 87–90.[CrossRef][Medline]

Groen, J., Gerding, M., Koeman, J. P., Roholl, P. J., van Amerongen, G., Jordans, H. G., Niesters, H. G. & Osterhaus, A. D. (1995). A macaque model for hantavirus infection. J Infect Dis 172, 38–44.[Medline]

Hamano, T., Matsuo, K., Hibi, Y., Victoriano, A. F., Takahashi, N., Mabuchi, Y., Soji, T., Irie, S., Sawanpanyalert, P. & other authors (2007). A single-nucleotide synonymous mutation in the gag gene controlling human immunodeficiency virus type 1 virion production. J Virol 81, 1528–1533.[Abstract/Free Full Text]

Kallio, E. R., Klingström, J., Gustafsson, E., Manni, T., Vaheri, A., Henttonen, H., Vapalahti, O. & Lundkvist, Å. (2006). Prolonged survival of Puumala hantavirus outside the host: evidence for indirect transmission via the environment. J Gen Virol 87, 2127–2134.[Abstract/Free Full Text]

Kaukinen, P., Vaheri, A. & Plyusnin, A. (2005). Hantavirus nucleocapsid protein: a multifunctional molecule with both housekeeping and ambassadorial duties. Arch Virol 150, 1693–1713.[CrossRef][Medline]

Klingström, J., Plyusnin, A., Vaheri, A. & Lundkvist, Å. (2002). Wild-type Puumala hantavirus infection induces cytokines, C-reactive protein, creatinine, and nitric oxide in cynomolgus macaques. J Virol 76, 444–449.[Abstract/Free Full Text]

Lundkvist, Å., Cheng, Y., Sjölander, K. B., Niklasson, B., Vaheri, A. & Plyusnin, A. (1997). Cell culture adaptation of Puumala hantavirus changes the infectivity for its natural reservoir, Clethrionomys glareolus, and leads to accumulation of mutants with altered genomic RNA S segment. J Virol 71, 9515–9523.[Abstract]

Manrubia, S. C., Escarmis, C., Domingo, E. & Lazaro, E. (2005). High mutation rates, bottlenecks, and robustness of RNA viral quasispecies. Gene 347, 273–282.[CrossRef][Medline]

Meyer, B. J. & Schmaljohn, C. S. (2000). Persistent hantavirus infections: characteristics and mechanisms. Trends Microbiol 8, 61–67.[CrossRef][Medline]

Nakamura, Y., Gojobori, T. & Ikemura, T. (2000). Codon usage tabulated from international DNA sequence databases: status for the year 2000. Nucleic Acids Res 28, 292[Abstract/Free Full Text]

Nemirov, K., Lundkvist, Å., Vaheri, A. & Plyusnin, A. (2003). Adaptation of Puumala hantavirus to cell culture is associated with point mutations in the coding region of the L segment and in the noncoding regions of the S segment. J Virol 77, 8793–8800.[Abstract/Free Full Text]

Nichol, S. T., Beaty, B. J., Elliott, R. M., Goldbach, R., Plyusnin, A., Schmaljohn, C. S. & Tesh, R. B. (2005). Family Bunyaviridae. In Virus Taxonomy: Eighth Report of the International Committee on Taxonomy of Viruses, pp. 695–716. Edited by C. M. Fauquet, M. A. Mayo, J. Maniloff, U. Desselberger & L. A. Ball. San Diego, CA: Elsevier Academic Press.

Novella, I. S., Zarate, S., Metzgar, D. & Ebendick-Corpus, B. E. (2004). Positive selection of synonymous mutations in vesicular stomatitis virus. J Mol Biol 342, 1415–1421.[CrossRef][Medline]

Plyusnin, A. & Morzunov, S. P. (2001). Virus evolution and genetic diversity of hantaviruses and their rodent hosts. Curr Top Microbiol Immunol 256, 47–75.[Medline]

Plyusnin, A., Vapalahti, O., Ulfves, K., Lehväslaiho, H., Apekina, N., Gavrilovskaya, I., Blinov, V. & Vaheri, A. (1994). Sequences of wild Puumala virus genes show a correlation of genetic variation with geographic origin of the strains. J Gen Virol 75, 405–409.[Abstract/Free Full Text]

Plyusnin, A., Vapalahti, O., Lehväslaiho, H., Apekina, N., Mikhailova, T., Gavrilovskaya, I., Laakkonen, J., Niemimaa, J., Henttonen, H. & Brummer-Korvenkontio, M. (1995). Genetic variation of wild Puumala viruses within the serotype, local rodent populations and individual animal. Virus Res 38, 25–41.[CrossRef][Medline]

Plyusnin, A., Cheng, Y., Lehväslaiho, H. & Vaheri, A. (1996a). Quasispecies in wild-type Tula hantavirus populations. J Virol 70, 9060–9063.[Abstract]

Plyusnin, A., Vapalahti, O. & Vaheri, A. (1996b). Hantaviruses: genome structure, expression and evolution. J Gen Virol 77, 2677–2687.[Abstract/Free Full Text]

Plyusnin, A., Hörling, J., Kanerva, M., Mustonen, J., Cheng, Y., Partanen, J., Vapalahti, O., Kukkonen, S. K., Niemimaa, J. & other authors (1997). Puumala hantavirus genome in patients with nephropathia epidemica: correlation of PCR positivity with HLA haplotype and link to viral sequences in local rodents. J Clin Microbiol 35, 1090–1096.[Abstract]

Sironen, T., Vaheri, A. & Plyusnin, A. (2001). Molecular evolution of Puumala hantavirus. J Virol 75, 11803–11810.[Abstract/Free Full Text]

Vapalahti, O., Mustonen, J., Lundkvist, Å., Henttonen, H., Plyusnin, A. & Vaheri, A. (2003). Hantavirus infections in Europe. Lancet Infect Dis 3, 653–661.[CrossRef][Medline]

Vignuzzi, M., Stone, J. K., Arnold, J. J., Cameron, C. E. & Andino, R. (2006). Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature 439, 344–348.[CrossRef][Medline]

Received 4 December 2007; accepted 23 January 2008.


This article has been cited by other articles:


Home page
Mol Biol EvolHome page
C. Ramsden, E. C. Holmes, and M. A. Charleston
Hantavirus Evolution in Relation to Its Rodent and Insectivore Hosts: No Evidence for Codivergence
Mol. Biol. Evol., January 1, 2009; 26(1): 143 - 153.
[Abstract] [Full Text] [PDF]


Home page
J. Gen. Virol.Home page
M. Razzauti, A. Plyusnina, H. Henttonen, and A. Plyusnin
Accumulation of point mutations and reassortment of genomic RNA segments are involved in the microevolution of Puumala hantavirus in a bank vole (Myodes glareolus) population
J. Gen. Virol., July 1, 2008; 89(7): 1649 - 1660.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.
Agricola
Right arrow Articles by Sironen, T.
Right arrow Articles by Plyusnin, A.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS