# Isolation and Genetic Characterization of Puumala Orthohantavirus Strains from France

Johann Vulin, Séverine Murri, Sarah Madrières, Maxime Galan, Caroline

Tatard, Sylvain Piry, Gabriele Vaccari, Claudia D'agostino, Nathalie

Charbonnel, Guillaume Castel, et al.

► To cite this version:

Johann Vulin, Séverine Murri, Sarah Madrières, Maxime Galan, Caroline Tatard, et al. Isolation and Genetic Characterization of Puumala Orthohantavirus Strains from France. Pathogens, 2021, 10 (3), pp.349. 10.3390/pathogens10030349. [https://hal.inrae.fr/hal-03231544](https://hal.inrae.fr/hal-03231544)

HAL Id: hal-03231544

https://hal.inrae.fr/hal-03231544

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Isolation and Genetic Characterization of Puumala Orthohantavirus Strains from France

Johann Vulin<sup>1,*†</sup>, Séverine Murri<sup>1,†</sup>, Sarah Madrières<sup>1,2</sup>, Maxime Galan<sup>2</sup>, Caroline Tatard<sup>2</sup>, Sylvain Piry<sup>2</sup>, Gabriele Vaccari<sup>3</sup>, Claudia D'Agostino<sup>3</sup>, Nathalie Charbonnel<sup>2</sup>, Guillaume Castel<sup>2</sup> and Philippe Marianneau<sup>1</sup>

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Vulin, J.; Murri, S.; Madrières, S.; Galan, M.; Tatard, C.; Piry, S.; Vaccari, G.; D'Agostino, C.; Charbonnel, N.; Castel, G.; et al. Isolation and Genetic Characterization of Puumala Orthohantavirus Strains from France. *Pathogens* **2021**, *10*, 349. https://doi.org/10.3390/pathogens10030349

Academic Editor: Lawrence S. Young

Received: 2 February 2021

Accepted: 12 March 2021

Published: 16 March 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 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 (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Puumala orthohantavirus (PUUV) causes a mild form of haemorrhagic fever with renal syndrome (HFRS) called nephropathia epidemica (NE), regularly diagnosed in Europe. France represents the western frontier of the expansion of NE in Europe with two distinct areas: an endemic area (north-eastern France) where PUUV circulates in rodent populations, with the detection of many human NE cases, and a non-endemic area (south-western France) where the virus is not detected, with only a few human cases being reported. In this study, we describe the different stages of the isolation of two PUUV strains from two distinct French geographical areas: Ardennes (endemic area) and Loiret (non-endemic area). To isolate PUUV efficiently, we selected wild bank voles (*Myodes glareolus*, the specific reservoir of PUUV) captured in these areas and that were seronegative for anti-PUUV IgG (ELISA) but showed a non-negligible viral RNA load in their lung tissue (qRT-PCR). With this study design, we were able to cultivate and maintain these two strains in Vero E6 cells and also propagate both strains in immunologically neutral bank voles efficiently and rapidly. High-throughput and Sanger sequencing results provided a better assessment of the impact of isolation methods on viral diversity.

**Keywords:** Hantavirus; Puumala; isolation

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.

<footnote n="1">¹ Laboratoire de Lyon, ANSES, Unité de Virologie, 69007 Lyon, France; severine.murri@anses.fr (S.M.); sarah.madrieres@anses.fr (S.M.); philippe.marianneau@anses.fr (P.M.)</footnote>

<footnote n="2">² CBGP, INRAE, CIRAD, IRD, Institut Agro, Univ Montpellier, 34000 Montpellier, France; maxime.galan@inrae.fr (M.G.); caroline.tatard@inrae.fr (C.T.); sylvain.piry@inrae.fr (S.P.); nathalie.charbonnel@inrae.fr (N.C.); guillaume.castel@inrae.fr (G.C.)</footnote>

<footnote n="3">³ Department of Food Safety, Nutrition and Veterinary Public Health, Istituto Superiore di Sanità, 00161 Rome, Italy; gabriele.vaccari@iss.it (G.V.); claudia.dagostino@iss.it (C.D.)</footnote>

<footnote>* Correspondence: johann.vulin@anses.fr</footnote>

<footnote>† These authors contributed equally to the work.</footnote>
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⁻¹ 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 × 10⁶ 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 × 10⁵ 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⁴ focus forming units (FFU)/mL and 1.5 × 10⁴ 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-Harnies 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.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** PUUV isolation with Vero E6 cells. Kinetics of PUUV RNA quantity during isolation process (A). Circles represent supernatants sampling. Picture of viral titration results (B) for Harnies and Vouzon compared to Sotkamo referential strain.</figcaption></figure>

<table><caption>**Table 1.** Serological results of bank vole experimentally infected with lung homogenates from animals trapped from Harnies and Vouzon.</caption><thead><tr><th rowspan="3">Capture Area</th><th rowspan="3">Day Post-Infection</th><th rowspan="3">Animal Number</th><th colspan="3">Antibody IgG-NPuu</th></tr><tr><th colspan="2">IgG-ELISA</th><th rowspan="2">Titration<br/>Last + Dilution<br/>(OD)</th></tr><tr><th>OD</th><th>Conclusion</th></tr></thead><tbody><tr><td rowspan="4">Ardennes-Harnies</td><td>D7</td><td>811</td><td>0.090</td><td>-</td><td>ND</td></tr><tr><td>D7</td><td>812</td><td>0.058</td><td>-</td><td>ND</td></tr><tr><td>D9</td><td>811</td><td>0.140</td><td>+</td><td>200 (0.102)</td></tr><tr><td>D9</td><td>812</td><td>0.117</td><td>+</td><td>100 (0.104)</td></tr><tr><td rowspan="4">Loiret-Vouzon</td><td>D7</td><td>791</td><td>0.032</td><td>-</td><td>ND</td></tr><tr><td>D7</td><td>792</td><td>0.029</td><td>-</td><td>ND</td></tr><tr><td>D9</td><td>791</td><td>0.230</td><td>+</td><td>800 (0.105)</td></tr><tr><td>D9</td><td>792</td><td>0.227</td><td>+</td><td>800 (0.106)</td></tr></tbody></table>

OD: Optical Density. ND: Not Defined.
<table><caption>**Table 2.** RNA detection results of bank vole experimentally infected with lung homogenates from trapped animals at Hargnies and Vouzon.</caption><thead><tr><th rowspan="2">Strain</th><th rowspan="2">Day Post-Infection</th><th rowspan="2">Animal Number</th><th colspan="4">RT-PCR (Mean Ct)</th></tr><tr><th>Lung</th><th>Liver</th><th>Sera-Day 7</th><th>Sera-Day 9</th></tr></thead><tbody><tr><td rowspan="2">Ardennes-Hargnies</td><td rowspan="2">D9</td><td>811</td><td>16.5 ± 0.13</td><td>19.7 ± 0.12</td><td>29.2 ± 0.07</td><td>32.1 ± 0.05</td></tr><tr><td>812</td><td>19.1 ± 0.00</td><td>20.0 ± 0.08</td><td>30.6 ± 0.11</td><td>32.7 ± 0.11</td></tr><tr><td rowspan="2">Loiret-Vouzon</td><td rowspan="2">D9</td><td>791</td><td>19.0 ± 0.05</td><td>21.0 ± 0.17</td><td>30.7 ± 0.20</td><td>31.0 ± 0.04</td></tr><tr><td>792</td><td>19.1 ± 0.06</td><td>20.8 ± 0.07</td><td>29.9 ± 0.19</td><td>30.5 ± 0.00</td></tr></tbody></table>

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.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparison of S and M segment amino-acid sequences. Synthesis of Sanger sequencing of S and M segments for natural isolates (in natura) (sequence as reference), cell cultivated viruses (in vitro) and experimentally infected bank voles (in vivo). Bank vole and cell drawings were designed by brgfx/Freepik.</figcaption></figure>

We then performed HTS analyses to investigate these genetic variations further. For the method details, see Supplementary Information S3, S4, S5 and S6 [16–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 bio-
logical 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 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.

<figure><img src="image_6.png" /><figcaption>Figure 3. *Cont.*</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>(B) *In vivo*</figcaption></figure>

<figure><img src="image_8.png" /></figure>

**Figure 3.** Comparison of viral diversity after (A) in vitro and (B) in vivo isolation process. *p*-values were determined using the Kruskal-Wallis test (*, *p* < 0.05; **, *p* < 0.01), X² = 13.27, *p* = 0.0041, *n* = 40.

## 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-Harnies 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éditerrané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

1. 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. [CrossRef]

2. 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. [CrossRef]

3. 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 Hebд.* **2017**, *23*, 492–499.

4. 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. [CrossRef] [PubMed]

5. Reynes, J.-M.; Carli, D.; Thomas, D.; Castel, G. Puumala Hantavirus Genotypes in Humans, France, 2012–2016. *Emerg. Infect. Dis.* **2019**, *25*, 140–143. [CrossRef] [PubMed]

6. 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. [CrossRef]

7. 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. [CrossRef] [PubMed]

8. 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. [CrossRef]

9. 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. [CrossRef] [PubMed]

10. 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. [CrossRef]

11. 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. [CrossRef] [PubMed]

12. 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. [CrossRef]

13. 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. [CrossRef] [PubMed]

14. 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. [CrossRef]

15. Sironen, T.; Vaheri, A.; Plyusnin, A. Molecular Evolution of Puumala Hantavirus. *J. Virol.* **2001**, *75*, 11803–11810. [CrossRef]

16. 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. [CrossRef]

17. 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. [CrossRef] [PubMed]

18. Kircher, M.; Sawyer, S.; Meyer, M. Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. *Nucleic Acids Res.* **2011**, *40*, e3. [CrossRef]

19. Martin, J.-F. Creating error-proof indexes for high throughput sequencing. *Zenodo* **2019**. [CrossRef]
20. 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. [CrossRef]

21. 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. [CrossRef]

22. Gregori, J.; Perales, C.; Rodriguez-Frias, F.; Esteban, J.I.; Quer, J.; Domingo, E. Viral quasispecies complexity measures. *Virology* **2016**, *493*, 227–237. [CrossRef]

23. Zhao, L.; Illingworth, C.J.R. Measurements of intrahost viral diversity require an unbiased diversity metric. *Virus Evol.* **2019**, *5*, vey041. [CrossRef]

24. 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. [CrossRef] [PubMed]

25. 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. [CrossRef] [PubMed]

26. 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. [CrossRef]

27. 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. [CrossRef] [PubMed]