Orthohantaviruses represent an increasing threat to humans due to their worldwide distribution, the increase in the number of infections and the emergence or re-emergence of new viruses [
1]. Puumala virus (PUUV) is the main orthohantavirus circulating in France. This tri-segmented, enveloped RNA virus, hosted by a rodent, the bank vole (
Myodes glaerolus), causes a mild form of haemorrhagic fever with renal syndrome (HFRS) called nephropathia epidemica (NE) [
2]. About 100 cases of NE are reported annually in north-eastern France [
3].
In France, PUUV sequences have been studied in bank vole samples [
4] and more recently in human samples [
5], but these studied strains have never been cultured or maintained until now. Considering the adaptation capacity of PUUV to cell culture conditions and changes in infectivity induced by cell culture passages [
6], obtaining a “wild-type” virus is a crucial step towards a better understanding of the biology of this virus. Unfortunately, due to its slow growth in cell culture, PUUV is often very difficult to isolate [
7].
In this study, we applied several methods that enabled the identification of animals carrying the virus, to isolate and cultivate the corresponding strains and to study the effect of the isolation protocol on virus diversity. All of these methods are shown in detail in
Supplementary Information S1.
Firstly, we determined a method for selecting bank voles containing live virus from which to isolate PUUV strains. Seto et al. [
8] isolated PUUV from bank voles that harboured the orthohantavirus nucleocapsid protein (NP) in their lungs, but showed no antibodies against PUUV. This detection of infection coupled with the absence of antibodies probably corresponds to an active viraemic phase. We adapted the Seto et al. approach by screening our collections of bank vole samples (composed of voles captured in various French forests over the last 10 years [
4,
9]) for seronegative animals with high amounts of viral RNA in their lung tissue. Serum samples were assayed using IgG ELISA and viral RNA was extracted from lung homogenates. More specifically, we screened bank voles captured in two distinct French geographical areas: (i) an NE endemic area: Ardennes, where PUUV circulates in rodent populations and where many cases of human NE are detected; (ii) a non-endemic zone: Loiret, where PUUV circulates in rodent populations and where no or very few cases of NE are diagnosed [
4]. In our previous studies, we showed that the virus from these two areas represents two distinct phylogenetics clusters [
4,
10,
11]. We selected seven bank voles captured in Ardennes during autumn 2011 (out of a total of 201 animals captured during this trapping session) and one bank vole captured in Loiret in summer 2014 (out of 44 animals). For our first isolation assays, we used a sample from Ardennes named Hargnies/2011 with the smallest cycle threshold value (Ct; i.e., 18.90 cycles) and a lung tissue sample from Loiret named Vouzon/2014 that showed a Ct of 23.65 cycles.
From these two animals, we tried to isolate live viruses from lung homogenates (5% w/v). The conditions that led to the in vitro isolation of the two viral strains consisted of the use of the lung homogenate centrifuged (300 rpm for 10 min), not filtered, diluted to 10−1 and applied to low-passaged Vero E6 cells (16 passages) in DMEM (Dulbecco′s modified Eagle′s medium) supplemented with 5% foetal calf serum (Gibco) and cultured to hyper-confluence (1 × 106 cells for one well in P6). It was necessary to carry out three successive passages to reach a sufficient quantity of live viruses. These passages were carried out using “co-culture”: a quarter of the initially infected cells were transferred to a culture of non-infected cells (2 × 105 cells/well in P6) accompanied by 500 µL of the supernatant from the infected culture added to 2 mL of fresh media.
The viral load measured in the supernatant of Vero E6 cells increased substantially: −15.2 Ct and −12.8 Ct for Ardennes-Hargnies and Loiret-Vouzon, respectively (
Figure 1A). Moreover, at the end of the cell isolation, viral titre was 1.6 × 10
4 focus forming units (FFU)/mL and 1.5 × 10
4 FFU/mL for Ardennes-Hargnies and Loiret-Vouzon, respectively. Even if the final viral titers were similar, the viral titer and viral RNA quantities were slightly different during the two first passages. After titration, their focus sizes were different to that of the PUUV control (Sotkamo strain) serially passaged in Vero E6 cells (
Figure 1B). The focus size of cultivated viruses was smaller (in particular for Loiret-Vouzon) than Sotkamo foci. Phenotypic variations have already been described for PUUV titration assays [
12]. Altogether, these results show that we were able to select viral isolates from natural host tissue and cultivate them in cell culture.
Furthermore, we carried out an in vivo isolation process that consists of biological assays on bank voles that were maintained at the ISS Institute (Rome, Italy). We used the same lung homogenates already used for the in vitro isolation process. We performed a subcutaneous injection of lung homogenate (5%
w/
v) on two animals for each strain. None of the inoculated animals showed clinical evidence of viral infection. Seven days after infection, we collected blood samples via retro-orbital sinus puncture and carried out PUUV serological detection and viral RNA analyses. The animals’ sera (pre-diluted at 1/10) were screened using IgG ELISA with PUUV recombinant nucleocapsid (N-PUUV) protein and negative controls. Samples were considered positive if the optical density (OD) was greater than 0.100. We did not detect any anti-PUUV antibodies (
Table 1), but we found viral RNA in all sera analyzed. Viral loads were similar regardless of the PUUV strain (Ardennes-Hargnies and Loiret-Vouzon), indicating that our protocol led to the effective transmission of PUUV in these animals. Two days later, at nine days post-infection (D9), we euthanized the four bank voles by cervical dislocation. We found that all animals had seroconverted, with a higher antibody titre observed for the Loiret-Vouzon strain (1/800) than for the Ardennes-Hargnies strain (animal 1: 1/200 and animal 2: 1/100). The RNA viral loads in the sera had already decreased in all four animals analysed, but were still positive. We also analysed the RNA viral loads in the lungs and liver, i.e., organs in which PUUV antigens have previously been detected [
13,
14]. We found high viral RNA loads (
Table 2), which confirmed the effective infection of all animals by PUUV.
Molecular analyses were carried out as already described [
7] to assess the PUUV levels during bank vole infection. RNA was extracted from lung tissue samples of positive bank voles (natural populations and experimentally infected) and from the supernatant of Vero E6 cells, then sequenced to obtain the complete coding sequence of the S and M segments (for a list of sequences and accession numbers, see
Supplementary Table S2). In order to provide a better assessment of the impact of isolation methods on intra-host viral diversity, amino-acid sequences were compared between naturally infected bank vole (in natura), cell-cultivated virus (in vitro) and experimentally infected bank vole (in vivo) for each strain (Ardennes-Hargnies or Loiret-Vouzon). For the M segment (3444 bp for coding region), the PUUV nucleotide sequences from Ardennes-Hargnies or Loiret-Vouzon were 100% identical in the three conditions analyzed by strains. For the S segment (1304 bp for coding region), the amino-acid sequences from Ardennes-Hargnies or Loiret-Vouzon revealed one mismatch between the sequences from bank vole (in natura and in vivo) and in vitro (
Figure 2). For Ardennes-Hargnies, bank voles (in natura and in vivo) showed an arginine at codon 63 (R63), whereas the virus cultivated in cells (in vitro) presented a glutamine (Q63). For Loiret-Vouzon, a difference was found in position 28 with an alanine residue (A28) for bank voles (in natura and in vivo) and serine (S28) for in vitro. It is important to note that these two mutations are located in the major antigenic domain (MAD) previously described [
15]. Therefore, the in vitro isolation protocol seemed to have influenced the amino-acid sequences of both PUUV strains.
We then performed HTS analyses to investigate these genetic variations further. For the method details, see
Supplementary Information S3, S4, S5 and S6 [
16,
17,
18,
19,
20,
21]. More specifically, we assessed whether the different isolation protocols maintained the initial intra-host viral diversity detected in the original samples (from which the strains were isolated) presented as “in natura” results below. We thus characterized the viral variants pool for the three viral “sources”: in natura, in vitro, and in vivo. Two indices were used to compare the viral diversity between the two isolation protocols and the two strains: the number of polymorphic sites [
22,
23] and the mean percent complexity (the number of unique sequence reads/total reads × 100) calculated on the 10 amplicons [
24]. A Kruskal–Wallis test followed by a Dunn multiple comparison test was conducted to compare the mean percent complexity between strains of different origins (in natura, in vitro or in vivo). For both propagation processes—in vivo and in vitro—the RNA viral sequences discussed after came from endpoint processes (biological assay at D9 for in vivo and D14 of 3rd passage for in vitro).
In natura, the intra-host viral diversity of the Ardennes-Hargnies and Loiret-Vouzon strains were similar. We found no difference in the number of polymorphic sites or in the mean percent complexity between areas. After the in vitro isolation process, one variant of the Ardennes-Hargnies strain (R63), mostly found in natura, switched to another major variant, such as Q63. Nevertheless, R63 was still present at a high frequency (about 23%) instead of a viral diversity decreased after the cell culture passages, regardless of the diversity index considered. The Loiret-Vouzon strain, as observed with the Sanger method, showed a change in the major variant from A28 to S28. However, HTS results were contradictory: the number of polymorphic sites increased, but the mean percent complexity remained stable. The in vitro isolation process lead to a significant difference in viral diversity between PUUV strains, but not in the same way for both strains (
Figure 3A).
Regarding the in vivo isolation process, we only analyzed lungs at endpoint of biological assay (D9). We did not perform any kinetic or multi-organ analyses. We detected a lower number of polymorphic sites compared with the natural strains. However, we found no significant difference in the mean percent complexity between PUUV strains or in in natura and in vitro conditions (
Figure 3B). The sequence of the majority variant did not differ between these two conditions, regardless of the isolation protocol considered. Results were similar in bank voles infected with the same strain.
Conclusions
In this study, we were able, for the first time, to isolate and maintain in cell culture two PUUV strains from two distinct French areas. Molecular analyses of the S and M segments of PUUV originating from natural isolates, from experimentally infected bank voles and from cell cultures revealed only one amino-acid mismatch for the Ardennes-Hargnies and Loiret-Vouzon strains. Both mismatches were identified in cell culture. HTS confirmed this result, thereby leading to the definition of a “bank vole” and a “cell culture” molecular profile. Having these two types of PUUV wild strain cultures is an important asset. Cell-adapted strains provide well-characterized viruses that can be used for antiviral candidate studies or cell-level experiments to assess the role of apoptosis, PUUV propagation or control [
25]. Furthermore, wild strains maintained in their natural host can help contribute to improving our knowledge of the PUUV ecology and evolution [
26]. The bank voles used in the in vivo experiments came from a “neutral” origin, i.e., not from Ardennes or Loiret. For the moment, bank vole-PUUV interactions are poorly documented, but may be a key point in viral diversity.
Previous studies showed that there is significant inter-regional viral genetic diversity [
4,
5,
27], which may explain, at least in part, the regional differences in NE epidemiological status in France. Here, we showed that the viral diversity of PUUV circulating in the NE endemic area (Ardennes) and in the NE non-endemic area (Loiret) evolves differently when passaged on Vero E6 cells. Further research dedicated to the characterization of the now-available French Puumala strains will help further our knowledge on the NE epidemiological situation.
Supplementary Materials
The following are available online at
https://www.mdpi.com/2076-0817/10/3/349/s1, Supplementary information S1: methods description, Table S2: PUUV sequences obtained in this study, Table S3: Sequences of primer cocktails used in Miseq, Supplementary information S4: Library preparation adapted from Galan et al. 2018, Supplementary information S5: Definition of the 0.32% threshold, Supplementary information S6: Means number of read used polymorphic sites quantification.
Author Contributions
Conceptualization, J.V., S.M. (Séverine Murri), G.C. and P.M.; methodology, J.V., S.M. (Séverine Murri), S.M. (Sarah Madrières), M.G.; validation, G.C., P.M.; formal analysis, J.V., S.M. (Séverine Murri), S.M. (Sarah Madrières), C.T., M.G. and G.C.; resources, J.V., S.M. (Séverine Murri), S.M. (Sarah Madrières), M.G., C.T., S.P., G.V., C.D., N.C., G.C. and P.M.; writing—original draft preparation, J.V., S.M. (Séverine Murri), S.M. (Sarah Madrières), G.C. and P.M.; writing—review and editing, J.V., S.M. (Séverine Murri), S.M. (Sarah Madrières), M.G., C.T., G.V., C.D., N.C., G.C. and P.M.; supervision, P.M. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
All animal studies were conducted in accordance with French and European regulations on the care and protection of laboratory animals (French Decree No 2013-118 of 1 February 2013 and Directive 2010/63/EU of 22 September 2010).
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Acknowledgments
Data used in this work were partly produced through sequencing facilities at ISEM (Institut des Sciences et de l’Evolution-Montpellier) and the LabEx CEMEB (Centre Méditérranéen de l’Environnement et de la Biodiversité). Sarah Madrières was funded by an INRA-EFPA/ANSES fellowship.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Lee, H.W.; Vaheri, A.; Schmaljohn, S. Discovery of hantaviruses and of the Hantavirus genus: Personal and historical perspectives of the Presidents of the International Society of Hantaviruses. Virus Res. 2014, 187, 2–5. [Google Scholar] [CrossRef]
- Brummer-Korvenkontio, M.; Vaheri, A.; Hovi, T.; Von Bonsdorff, C.-H.; Vuorimies, J.; Manni, T.; Penttinen, K.; Oker-Blom, N.; Lahdevirta, J. Nephropathia Epidemica: Detection of Antigen in Bank Voles and Serologic Diagnosis of Human Infection. J. Infect. Dis. 1980, 141, 131–134. [Google Scholar] [CrossRef]
- Reynes, J.M.; Renaudin, B.; Fizet, A.; Bour, J.B.; Brodard, V.; Cart-Tanneur, E.; Dewilde, A.; El Hamri, M.; Fleury, H.; Hecquet, D.; et al. Surveillance of human hantavirus infections in metropolitan France, 2012–2016. Bull. Epidémiologique Hebd. 2017, 23, 492–499. [Google Scholar]
- Castel, G.; Couteaudier, M.; Sauvage, F.; Pons, J.-B.; Murri, S.; Plyusnina, A.; Pontier, D.; Cosson, J.-F.; Plyusnin, A.; Marianneau, P.; et al. Complete Genome and Phylogeny of Puumala Hantavirus Isolates Circulating in France. Viruses 2015, 7, 5476–5488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Reynes, J.-M.; Carli, D.; Thomas, D.; Castel, G. Puumala Hantavirus Genotypes in Humans, France, 2012–2016. Emerg. Infect. Dis. 2019, 25, 140–143. [Google Scholar] [CrossRef] [PubMed]
- Lundkvist, A.; Cheng, Y.; Sjolander, K.B.; Niklasson, B.; Vaheri, A.; Plyusnin, A. 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. 1997, 71, 9515–9523. [Google Scholar] [CrossRef] [Green Version]
- Binder, F.; Reiche, S.; Roman-Sosa, G.; Saathoff, M.; Ryll, R.; Trimpert, J.; Kunec, D.; Höper, D.; Ulrich, R.G. Isolation and characterization of new Puumala orthohantavirus strains from Germany. Virus Genes 2020, 56, 448–460. [Google Scholar] [CrossRef] [PubMed]
- Seto, T.; Tkachenko, E.A.; Morozov, V.G.; Tanikawa, Y.; Kolominov, S.I.; Belov, S.N.; Nakamura, I.; Hashimoto, N.; Kon, Y.; Balakiev, A.E.; et al. An efficient in vivo method for the isolation of Puumala virus in Syrian hamsters and the characterization of the isolates from Russia. J. Virol. Methods 2011, 173, 17–23. [Google Scholar] [CrossRef]
- Klingström, J.; Heyman, P.; Escutenaire, S.; Sjölander, K.B.; De Jaegere, F.; Henttonen, H.; Lundkvist, Å. Rodent host specificity of European hantaviruses: Evidence of Puumala virus interspecific spillover. J. Med. Virol. 2002, 68, 581–588. [Google Scholar] [CrossRef] [PubMed]
- Castel, G.; Chevenet, F.; Razzauti, M.; Murri, S.; Marianneau, P.; Cosson, J.-F.; Tordo, N.; Plyusnin, A. Phylogeography of Puumala orthohantavirus in Europe. Viruses 2019, 11, 679. [Google Scholar] [CrossRef] [Green Version]
- Murri, S.; Madrières, S.; Tatard, C.; Piry, S.; Benoit, L.; Loiseau, A.; Pradel, J.; Artige, E.; Audiot, P.; Leménager, N.; et al. Detection and Genetic Characterization of Puumala Orthohantavirus S-Segment in Areas of France Non-Endemic for Nephropathia Epidemica. Pathogens 2020, 9, 721. [Google Scholar] [CrossRef] [PubMed]
- Sundstrom, K.B.; Stoltz, M.; Lagerqvist, N.; Lundkvist, A.; Nemirov, K.; Klingstrom, J. Characterization of two substrains of Puumala virus that show phenotypes that are different from each other and from the original strain. J. Virol. 2011, 85, 1747–1756. [Google Scholar] [CrossRef] [Green Version]
- Gavrilovskaya, I.N.; Apekina, N.S.; Myasnikov, Y.A.; Bernshtein, A.D.; Ryltseva, E.V.; Gorbachkova, E.A.; Chumakov, M.P. Features of circula-tion of hemorrhagic fever with renal syndrome (HFRS) virus among small mammals in the European U.S.S.R. Arch. Virol. 1983, 75, 313–316. [Google Scholar] [CrossRef] [PubMed]
- Yanagihara, R.; Amyx, H.L.; Gajdusek, D.C. Experimental infection with Puumala virus, the etiologic agent of nephropathia epidemica, in bank voles (Clethrionomys glareolus). J. Virol. 1985, 55, 34–38. [Google Scholar] [CrossRef] [Green Version]
- Sironen, T.; Vaheri, A.; Plyusnin, A. Molecular Evolution of Puumala Hantavirus. J. Virol. 2001, 75, 11803–11810. [Google Scholar] [CrossRef] [Green Version]
- Gouy, M.; Guindon, S.; Gascuel, O. SeaView Version 4: A Multiplatform Graphical User Interface for Sequence Alignment and Phylogenetic Tree Building. Mol. Biol. Evol. 2009, 27, 221–224. [Google Scholar] [CrossRef] [Green Version]
- Galan, M.; Pons, J.-B.; Tournayre, O.; Pierre, É.; Leuchtmann, M.; Pontier, D.; Charbonnel, N. Metabarcoding for the parallel identification of several hundred predators and their prey: Application to bat species diet analysis. Mol. Ecol. Resour. 2018, 18, 474–489. [Google Scholar] [CrossRef] [PubMed]
- Kircher, M.; Sawyer, S.; Meyer, M. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res. 2011, 40, e3. [Google Scholar] [CrossRef] [Green Version]
- Martin, J.-F. Creating error-proof indexes for high throughput sequencing. Zenodo 2019. [Google Scholar] [CrossRef]
- Escudié, F.; Auer, L.; Bernard, M.; Mariadassou, M.; Cauquil, L.; Vidal, K.; Maman, S.; Hernandez-Raquet, G.; Combes, S.; Pascal, G. FROGS: Find, Rapidly, OTUs with Galaxy Solution. Bioinformatics 2018, 34, 1287–1294. [Google Scholar] [CrossRef]
- Callahan, B.J.; Mcmurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [Green Version]
- Gregori, J.; Perales, C.; Rodriguez-Frias, F.; Esteban, J.I.; Quer, J.; Domingo, E. Viral quasispecies complexity measures. Virology 2016, 493, 227–237. [Google Scholar] [CrossRef]
- Zhao, L.; Illingworth, C.J.R. Measurements of intrahost viral diversity require an unbiased diversity metric. Virus Evol. 2019, 5, vey041. [Google Scholar] [CrossRef] [Green Version]
- Cousins, M.M.; Ou, S.-S.; Wawer, M.J.; Munshaw, S.; Swan, D.; Magaret, C.A.; Mullis, C.E.; Serwadda, D.; Porcella, S.F.; Gray, R.H.; et al. Comparison of a High-Resolution Melting Assay to Next-Generation Sequencing for Analysis of HIV Diversity. J. Clin. Microbiol. 2012, 50, 3054–3059. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Sola-Riera, C.; Gupta, S.; Ljunggren, H.G.; Klingstrom, J. Orthohantaviruses belonging to three phylogroups all inhibit apoptosis in infected target cells. Sci. Rep. 2019, 9, 834. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Madrières, S.; Castel, G.; Murri, S.; Vulin, J.; Marianneau, P.; Charbonnel, N. The Needs for Developing Experiments on Reservoirs in Hantavirus Research: Accomplishments, Challenges and Promises for the Future. Viruses 2019, 11, 664. [Google Scholar] [CrossRef] [Green Version]
- Monchatre-Leroy, E.; Murri, S.; Castel, G.; Calavas, D.; Boué, F.; Hénaux, V.; Marianneau, P. First insights into Puumala orthohantavirus circulation in a rodent population in Alsace, France. Zoonoses Public Health 2018, 65, 540–551. [Google Scholar] [CrossRef] [PubMed]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).