Next Article in Journal
HAS-Flow May Be an Adequate Method for Evaluating Human T-Cell Leukemia Virus Type 1 Infected Cells in Human T-Cell Leukemia Virus Type 1-Positive Rheumatoid Arthritis Patients Receiving Antirheumatic Therapies: A Retrospective Cross-Sectional Observation Study
Next Article in Special Issue
Characterization of a Panel of Cross-Reactive Hantavirus Nucleocapsid Protein-Specific Monoclonal Antibodies
Previous Article in Journal
Viral Hemorrhagic Septicemia Virus Activates Integrated Stress Response Pathway and Induces Stress Granules to Regulate Virus Replication
Previous Article in Special Issue
Molecular Characterisation and Phylogeny of Tula Virus in Kazakhstan
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection

by
Elisa Heuser
1,2,
Stephan Drewes
1,
Jakob Trimpert
3,
Dusan Kunec
3,
Calvin Mehl
1,2,
Marieke P. de Cock
4,
Ankje de Vries
4,
Christiane Klier
5,
Martin Oskamp
6,
Peter Tenhaken
7,
Fatima Hashemi
8,9,
Daniela Heinz
8,
Mariana Nascimento
3,
Marc Boelhauve
10,
Rasa Petraityte-Burneikiene
11,
Dina Raafat
12,13,
Miriam Maas
4,
Detlev H. Krüger
14,
Andreas Latz
8,
Jörg Hofmann
14,
Gerald Heckel
15,
Johannes Dreesman
5 and
Rainer G. Ulrich
1,2,*
add Show full author list remove Hide full author list
1
Institute of Novel and Emerging Infectious Diseases, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, 17493 Greifswald, Germany
2
German Center for Infection Research (DZIF), Partner Site Hamburg-Lübeck-Borstel-Riems, 17493 Greifswald, Germany
3
Institute of Virology, Department of Veterinary Medicine, Freie Universität Berlin, 14163 Berlin, Germany
4
Centre for Infectious Disease Control, National Institute for Public Health and the Environment (RIVM), 3720 BA Bilthoven, The Netherlands
5
Department of Infectious Diseases, Public Health Agency of Lower Saxony, 30449 Hannover, Germany
6
Health Department of Grafschaft Bentheim County, 48527 Nordhorn, Germany
7
Health service for the district and city of Osnabrück, 49082 Osnabrück, Germany
8
Veterinary Department, NovaTec Immundiagnostica GmbH, 64859 Dietzenbach, Germany
9
University of Applied Sciences, 60318 Frankfurt am Main, Germany
10
Department of Agriculture, South Westphalia University of Applied Sciences, 59494 Soest, Germany
11
Institute of Biotechnology, Life Sciences Center, Vilnius University, LT-10257 Vilnius, Lithuania
12
Institute of Immunology, University Medicine Greifswald, 17475 Greifswald, Germany
13
Department of Microbiology and Immunology, Faculty of Pharmacy, Alexandria University, Alexandria 21521, Egypt
14
Institute of Virology, Charité–University Medicine Berlin, 10117 Berlin, Germany
15
Institute of Ecology and Evolution, University of Bern, 3012 Bern, Switzerland
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(2), 467; https://doi.org/10.3390/v15020467
Submission received: 21 December 2022 / Revised: 24 January 2023 / Accepted: 27 January 2023 / Published: 7 February 2023
(This article belongs to the Special Issue Rodent-Borne Viruses 2.0)

Abstract

:
Seoul orthohantavirus (SEOV) is a rat-associated zoonotic pathogen with an almost worldwide distribution. In 2019, the first autochthonous human case of SEOV-induced hemorrhagic fever with renal syndrome was reported in Germany, and a pet rat was identified as the source of the zoonotic infection. To further investigate the SEOV reservoir, additional rats from the patient and another owner, all of which were purchased from the same vendor, were tested. SEOV RNA and anti-SEOV antibodies were found in both of the patient’s rats and in two of the three rats belonging to the other owner. The complete coding sequences of the small (S), medium (M), and large (L) segments obtained from one rat per owner exhibited a high sequence similarity to SEOV strains of breeder rat or human origin from the Netherlands, France, the USA, and Great Britain. Serological screening of 490 rats from breeding facilities and 563 wild rats from Germany (2007–2020) as well as 594 wild rats from the Netherlands (2013–2021) revealed 1 and 6 seropositive individuals, respectively. However, SEOV RNA was not detected in any of these animals. Increased surveillance of pet, breeder, and wild rats is needed to identify the origin of the SEOV strain in Europe and to develop measures to prevent transmission to the human population.

1. Introduction

Hantaviruses are (re-)emerging zoonotic pathogens that are transmitted by small mammals, especially rats, voles, and mice [1]. Their genome consists of three RNA segments of negative polarity: the small (S) and medium (M) segments encode the structural proteins of the virion, the nucleocapsid (N) protein, and the glycoprotein precursor, respectively. The large (L) segment encodes the RNA-directed RNA polymerase [2]. In humans, hantavirus infection can cause hemorrhagic fever with renal syndrome (HFRS) and hantavirus cardiopulmonary syndrome. Seoul orthohantavirus (SEOV) causes a mild to moderate form of HFRS with initial symptoms such as fever, chills, headache, vomiting, and rash. The case fatality rate was estimated at about 1% [3,4]. Hantavirus isolation or viral RNA recovery from human samples are uncommon because the virus is often cleared in the organism before patient samples are collected for diagnostic approaches [5].
SEOV was first discovered in a Norway (or brown) rat (Rattus norvegicus) captured in Seoul, South Korea, in 1980 [6]. The virus was later detected in other parts of Asia, although with a lower prevalence than the Hantaan orthohantavirus (HTNV) [7,8]. Further investigations in Asia detected SEOV not only in Norway rats but also in black (or roof) rats (Rattus rattus) and other Rattus species [8]. Since then, SEOV has been detected on almost all continents following the global distribution of these synanthropic rodents [9,10,11,12].
In Europe, several studies evaluated the occurrence of SEOV in pet, breeder, and wild rats and detected human disease cases caused by SEOV transmission by pet rats [13,14]. SEOV, which was detected in pet and breeder rats from Great Britain [15] and the Netherlands [16] as well as in a pet rat in Sweden that had been imported from England [17], was associated with human disease cases [13,14]. Studies on wild rats from Belgium indicated the presence of SEOV-reactive antibodies [18]. Furthermore, studies in the Eastern province of Gelderland in the Netherlands and in pig farms in Great Britain in 2013 showed SEOV RNA in wild rats [19,20]. In Germany, a single human SEOV infection was reported in 2018 in a traveler returning from Southeast Asia [7], while other studies using human serum panels from Germany failed to detect SEOV-reactive antibodies using neutralization assays [21,22].
In October 2019, a female patient with a high fever who was admitted to the intensive care unit of a hospital in Lower Saxony (Germany) in critical condition was diagnosed with SEOV infection via conventional reverse transcription-polymerase chain reaction (RT-PCR). Her pet rat was presumed to be the source of infection [23]. An analysis of partial SEOV sequences from the patient and her pet rat indicated a large similarity to pet-and breeder-rat-derived sequences from the Netherlands, France, Great Britain, and the USA, but they differed from sequences found in wild rats from the Netherlands, Great Britain, and other countries. At that time, the occurrence of SEOV in breeder and wild rats in Germany was not investigated.
The two main objectives of this study were: (i) to identify more precisely the origin of the human SEOV infection in Germany based on epidemiological trace-back and trace-forward investigations; and (ii) to evaluate the potential origin of the SEOV strain in the patient by investigating pet, breeder, and wild rats from both Germany and the Netherlands.

2. Materials and Methods

2.1. Epidemiological Investigation of the Origin of the Human SEOV Infection in Germany

Once a pet rat was identified as the likely source of the human SEOV infection in 2019 [23], the origin of the pet rats was followed by tracing the distribution channels of the pet rats back and forth (Supplementary Figure S1). Local public health and veterinary authorities as well as members of the research team methodologically conducted interviews and home visits with the patient (pet-rat owner 1), the vendor, and another pet-rat owner (pet-rat owner 2) who purchased the rats from the same vendor.

2.2. Collection of Rats

In addition to the patient’s first pet rat (KS19/1354) [23], three more pet rats from pet-rat owner 1 and three from pet-rat owner 2 were collected, euthanized, and sent to the National Reference Laboratory for Hantaviruses at the Friedrich-Loeffler-Institut (FLI), Greifswald-Insel Riems, Germany. In addition, we included materials from earlier studies comprising samples of rats from breeding facilities and wild populations across Germany (2007 to 2019 [24,25]), from wild populations from rural and urban areas in the Netherlands (2013–2021), as well as from ongoing studies of commensal rodent pathogens in Germany (2020; see Figure 1 and Supplementary Table S1). The rats from breeding facilities in Germany comprised laboratory Norway rats, wild-trapped Norway and black rats, as well as feeder rats from a zoo. Rats captured in 2013–2014 in the Netherlands were live-trapped with methods described previously [26]; other rats were captured using snap-traps. All carcasses were stored frozen at −20 °C until dissection.

2.3. Dissection of Rats

Following standard protocols, the frozen rat carcasses were thawed at 4 °C and dissected. A standardized panel of organs (lung, heart, liver, spleen, blood/transudate, and kidney) was taken from rats in Germany for use in the screening for SEOV and stored at −20 °C until further investigation. In the Netherlands, lung samples were collected and stored in RNAlater® (Thermo Fisher Scientific, Landsmeer, The Netherlands) for 3 to 5 days before storage at −80 °C until they were further processed.

2.4. Rat Control Sera

SEOV-positive control sera (L and R) were generated via four-time immunization of two Norway rats in four-week intervals with 10 µg of yeast-expressed, purified, and His-tagged SEOV N protein [27]. For immunization, the antigen was mixed 1:1 with 100 µL of Adjuvans Gerbu MM (Gerbu Biotechnik GmbH, Heidelberg, Germany) and injected subcutaneously. Serum and thoracic lavage from a non-infected Norway rat from the breeding facility at FLI were used as a negative control.

2.5. Screening Strategy for the Different Rat Panels

As this study combined results from rats belonging to four different panels (I–IV) that were collected over a long period of time (15 years), slightly different procedures were used over the years (see Supplementary Table S1 and below). The pet rats from Germany (panel I) were investigated using a SEOV-immunoglobulin G (IgG) enzyme-linked immunosorbent assay (ELISA) (A), a hantavirus IgG ELISA (B), a line-blot assay, and conventional RT-PCR assay. The other rats from Germany (panels II and III) were investigated using an IgG ELISA (B), a line-blot assay, and RT-PCR. Rats collected in the Netherlands (panel IV) were all tested using an SEOV-IgG ELISA (C) and/or real-time RT-PCR [28]. All seropositive rat samples from the Netherlands (collected during 2016–2021) were also investigated via a virus neutralization test that used SEOV, Puumala orthohantavirus (PUUV), Dobrava-Belgrade orthohantavirus (DOBV), and Tula orthohantavirus (TULV) as previously described [29].

2.6. ELISA

2.6.1. SEOV-IgG ELISA (A)

For the samples collected from German pet rats (panel I), a SEOV-IgG ELISA was performed using thoracic lavage (diluted 1:10) and whole blood (diluted 1:200) of the rats. A yeast-expressed, His-tagged, and affinity-purified SEOV N-protein served as the antigen [27]. A microtiter plate (C96 PolySorp™, NUNC, VWR, Darmstadt, Germany) was coated with 100 µg/well of the antigen in carbonate buffer and incubated overnight at 4 °C. Thereafter, the plate was blocked with 1% bovine serum albumin (BSA; Serva, Heidelberg, Germany) in phosphate-buffered saline (PBS)/0.05% Tween20 (Merck, Darmstadt, Germany) for 1 h at room temperature. The plate was then inoculated with the samples (diluted with 0.5% BSA in PBS/0.05% Tween20) and incubated for 1 h at 37 °C, which was followed by a wash step with washing buffer (PBS/0.1% Tween20), then the addition of horseradish peroxidase (HRP)-labeled anti-rat IgG conjugate (1:12,500 100 µL) and a further 1 h incubation at 37 °C. The plate was washed again and then incubated with 100 µL of 3,3′,5,5′-tetramethylbenzidine (TMB) (BioRad, Munich, Germany) for 10 min, and the reaction was then stopped by adding 100 µL of 1 M sulfuric acid. Finally, the plate was measured in an ELISA plate reader Infinite® M200 Pro (Tecan, Männedorf, Switzerland) at 450 nm. The upper and lower cut-off values were defined as previously described [30]. Samples with an optical density (OD) value below the lower cut-off value were defined as negative and those with OD values above the upper cut-off value as positive. OD values lying between the upper and lower cut-off values were considered to be inconclusive. Positive and inconclusive samples were retested according to a previously published decision tree [30].

2.6.2. Hantavirus-IgG ELISA (B)

For the hantavirus IgG ELISA (NovaTec Immundiagnostica GmbH, Dietzenbach, Germany) that targeted multiple hantaviruses, medium-binding 96-well ELISA plates (Greiner, Nürtingen, Germany) were coated with a mixture of four antigens (DOBV, PUUV, HTNV, and SEOV) using 100 ng/well of each individual antigen in carbonate-buffer (NovaTec Immundiagnostica GmbH) and left overnight at 4 °C. The plates were washed once with 250 µL wash-after-coat buffer (NovaTec Immundiagnostica GmbH) and blocked with stabilization and blocking buffer (NovaTec Immundiagnostica GmbH) for 1 h at room temperature. After drying for 30 min at 37 °C, the plates were stored at 4 °C until further use. Thoracic lavage or blood samples of rat panels I-III were diluted 1:100 in sample dilution buffer (NovaTec Immundiagnostica GmbH). Aliquots of diluted samples (100 µL/well) were incubated for 1 h at 37 °C followed by 3 washes with washing buffer (0.2 M phosphate, pH = 7.2; NovaTec Immundiagnostica GmbH) and 30 min of incubation at room temperature with 100 µL/well of a titrated dilution of a HRP-labeled protein A/G conjugate (NovaTec Immundiagnostica GmbH; final dilution 1:900,000 of a 5 mg/mL stock solution). After three more washes with washing buffer, the plates were incubated for 15 min at room temperature with TMB substrate solution (NovaTec Immundiagnostica GmbH). The reaction was stopped by adding 100 µL/well of stop solution (0.2 M sulfuric acid; NovaTec Immundiagnostica GmbH) and evaluated by measuring at 450 nm in a plate reader (Anthos Labtec Instruments, Wals-Siezenheim, Austria).

2.6.3. SEOV-IgG ELISA (C)

Sera taken from rats from the Netherlands (panel IV) were tested for the presence of anti-SEOV IgG antibodies using a commercial human DOBV/HTNV IgG ELISA (Progen Biotechnik GmbH, Heidelberg, Germany), that was adapted to detect rat IgG antibodies by using rabbit-anti-rat HRP-labeled IgG conjugate (Sigma–Aldrich Chemie B.V. Zwijndrecht, The Netherlands) at a 1:5000 dilution. The optimal OD cut-off value was determined through the use of a binary mixture model as previously described [26].

2.7. Line-Blot Analysis

For the line-blot format, the recombinant antigens of DOBV, PUUV, HTNV, SEOV, Andes orthohantavirus, and Sin Nombre orthohantavirus were printed with a dispenser (FrontLine HR microliter contact; BioDot, Irvine, CA, USA) on a nitrocellulose membrane (GE Healthcare, Chicago, IL, USA) together with a control line for sample loading and conjugate function. After drying, the membranes were cut into 3 mm strips and stored at 4 °C. Prior to use, strips were equilibrated in 1 mL of sample dilution buffer (10 mM phosphate buffer, pH 7.2; NovaTec Immundiagnostica GmbH). Rat samples (panels I–III) were added to the membranes, and then the membranes were incubated with gentle shaking for 1 h at room temperature. After 3 washes with 1 mL of washing buffer (0.2 M phosphate, pH = 7.2; NovaTec Immundiagnostica GmbH) for 5 min each, the membranes were incubated with gentle shaking for 30 min with 1 mL of HRP-labeled protein A/G conjugate (NovaTec Immundiagnostica GmbH) at room temperature. The strips were then washed three times with 1 mL of washing buffer for 5 min. Signals were developed by incubating the membranes in 1 mL of TMB substrate solution (NovaTec Immundiagnostica GmbH) for 15 min with gentle shaking at room temperature. The reaction was stopped by the addition of at least 1 mL of distilled water. After drying the membranes for at least 30 min at room temperature, the strips were visually evaluated.

2.8. SEOV RNA Detection and Sequence Determination

For SEOV RNA detection and tissue-distribution analysis in rats from Germany, nucleic acids were extracted manually from homogenized heart, lung, liver, kidney, and spleen tissues using QIAzol Lysis Reagent (Qiagen, Hilden, Germany) or a NucleoMag® VET kit (Macherey-Nagel, Düren, Germany) according to the manufacturer’s instructions on a KingFisher™ Flex Purification System (Thermo Fisher Scientific, Waltham, MA, USA). The pet rats of the patient/pet-rat owner 1 were initially investigated using a hantavirus L segment-specific conventional one-step RT-PCR using primers Han-LF1: ATG TAY GTB AGT GCW GAT GC and Han-LR1: AAC CAD TCW GTY CCR TCA TC [31] and thereafter with conventional S segment-RT-PCR using the following primers: SEOV_415_alt_fw: TAC ATG TTA ACA ACA AGA GG and SEOV_1130_alt_rev: TCC AGT TGT ATT CCC ATT GAT TG for the S segment. These novel primers were designed based on an alignment of representative SEOV S segment sequences from Europe that included the pet-rat-derived sequence (KS19/1354). The three RT-PCR positive rats from pet-rat owner 1 were also investigated for the M segment using primers SEOV_M_47_fw: CTC TGG ACA GAT AAT GCT CAT G and SEOV_M_1068_rev: CCC AGC TAT TTT CAT ACT CAT AAT C. The rats from pet-rat owner 2 were investigated using the conventional S segment RT-PCR alone. The one-step RT-PCRs were performed using the SuperScript III RT-PCR Kit (Qiagen). Breeder and wild rats from Germany (panels II and III) were investigated using conventional S segment RT-PCR (above). RT-PCR products were resolved via electrophoresis in a 1.5% agarose gel with subsequent Safe View (Applied Biological Materials, Richmond, BC, Canada) or ethidium bromide (Sigma-Aldrich, St. Louis, MO, USA) staining. A BigDye™ Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems, Waltham, MA, USA) was used for Sanger sequencing. Consensus sequences were created using BioEdit v7.2.5 [32].
For the rats from the Netherlands (2013/2014, subset of panel IV), the lung tissue was manually disrupted in liquid nitrogen and RNA was extracted using a RNeasy isolation minikit (Qiagen) and including a Qiashredder step for homogenization and DNase I treatment. From 2016 onward, lung tissue samples were homogenized two times (40 sec at maximum speed) using MagNA Lyser Green Bead tube (Roche Diagnostics, GmbH, Mannheim, Germany) with 500 µL of lysis buffer (MagNa Pure 96 Total Nucleic Acid Isolation kit) on a FastPrep FP120 homogenizer (Thermo Fisher Scientific). Total nucleic acid was extracted using a MagNa Pure 96 Total Nucleic Acid Isolation kit on a MagNA pure 96 (Roche Diagnostics). Quality control of the RNA isolation and inhibition control was performed using β-actin RT-PCR. Samples were subsequently tested via real-time RT-PCR that targeted the S segment [28]. From 2013 to 2018, TaqMan Master (Roche Diagnostics) was used, while in subsequent years TaqMan Fast Virus 1-Step Master Mix was used (Thermo Fisher Scientific). The black rats that were captured in 2013–2014 were also tested using a nested conventional RT-PCR assay of the L segment as described previously [31].

2.9. Enrichment of Target Viral Sequences and High-Throughput Sequencing (HTS)

For the HTS libraries, 1 µg of total spleen and liver RNA from one rat from each pet-rat owner was processed using a NEBNext Ultra II RNA library preparation kit (New England Biolabs, Ipswich, MA, USA) according to the manufacturer’s instructions with eight minutes of RNA fragmentation and seven PCR cycles to amplify the resulting cDNA libraries. Amplified cDNA libraries were enriched for hantavirus sequences using a custom-made myBaits tiling array (Arbor Biosiences, Ann Arbor, MI, USA) designed to capture sequences from rodent-associated orthohantaviruses according to the manufacturer’s instructions. Briefly, 1 µg of HTS library was hybridized for 20 hours to biotinylated RNA baits, captured with streptavidin-coated paramagnetic beads, washed extensively to remove non-target sequences, and eluted from the beads through heat denaturation. Enriched sequences were amplified with 10 cycles of PCR using Q5 high-fidelity DNA polymerase (New England Biolabs) and primers annealed to the P5 and P7 library sequences. Enriched and amplified libraries were run on an Illumina MiSeq sequencer using v.3 chemistry for 2 × 300 bp paired-end reads.
Raw sequencing data were processed using Trimmomatic v.0.39 [33] and mapped against the genome of SEOV strain 201701093/SEOV/Illinois_US/Rat (GenBank accession: MK360784 (S segment), MK360797 (M segment), and MK360803 (L segment)) using the Burrows-Wheeler aligner v0.7.17 [34]. Mapping statistics were generated using Samtools v.1.10 [35], and alignments were visualized using IGV v2.9.4 for Linux [36]. For detection of single-nucleotide polymorphisms (SNPs), Freebayes (a Bayesian genetic variant detector) was used [37]. All SNPs with a minimum mapping quality of 5, a minimum count of 3, and a minimum fraction of 0.1 were considered. Consensus sequences for each sample were obtained using BCFtools [35].

2.10. Sequence Comparison and Phylogenetic Analyses

Sequences were identified using Basic Local Alignment Search Tool (BLAST, National Center for Biotechnology Information, NCBI) [38] and aligned using BioEdit v7.2.5 [32]. The best-fitting substitution model was determined using jModeltest 2.1.10 [39]. Phylogenetic tree reconstruction was performed using the CIPRES Science Gateway v3.3 with the implemented programs FastTreeMP on XSEDE v2.1.10 for maximum likelihood tree reconstruction and MrBayes on XSEDE v3.2.7a for Bayesian tree inference [40,41,42]. Subsequently, consensus trees were constructed based on Bayesian trees. For this purpose, bootstrap values were transferred from the maximum likelihood tree to the consensus tree if branches were congruent.

3. Results

3.1. Identification of the Origin of the Human SEOV Infection

The epidemiological investigation identified pet rats as a possible source of the laboratory-confirmed human SEOV infection in a patient (pet-rat owner 1) in Lower Saxony (Germany) [23]. A trace-back investigation revealed that the patient had purchased the infected index pet rat (KS19/1354) and two additional juvenile rats from a vendor in North Rhine-Westphalia (Germany) via an internet marketplace approximately three months before the onset of symptoms. The investigation also revealed that the patient already had owned an older female pet rat for some time before the onset of the disease that was kept separate from the new male rats (Supplementary Figure S1). After the first rat [23] and two additional rats (see below) tested positive, the patient agreed to euthanize and test the older fourth rat.
A trace-forward investigation at the vendor’s residence revealed that three other rats had been sold to another private customer (pet-rat owner 2) in Lower Saxony (Germany; Supplementary Figure S1). The responsible local health department contacted and informed pet-rat owner 2 about the potential risk of keeping pet rats that were potentially infected with SEOV. Pet-rat owner 2, who reported not having any signs of illness, agreed to sacrifice the rats for further investigations at the FLI, Germany, in December 2019.
The vendor himself had previously bought eight rats at a reptile market in North Rhine-Westphalia, Germany. There, the rats had been sold outside of the official market as feeder rats. No further trace-back was possible on the origin of the rats sold at the reptile market.

3.2. Analysis of Rats from Patient/Pet-Rat Owner 1 and Pet-Rat Owner 2

The analysis of the partial sequence of one of the patient’s pet rats (KS19/1354) revealed the potential origin of the patient’s infection due to a nucleotide sequence identity of 99.8% for the S segment and a 100% identity for the L segment (data not shown [23]). The patient’s other two juvenile rats (KS19/1352 and KS19/1353) also tested positive in the conventional S and M segment RT-PCRs, but only one (KS19/1352) was positive in the conventional L segment RT-PCR. The patient’s fourth, older rat (KS19/1414), which was obtained from a different source, tested negative in the S segment RT-PCR assay (Table 1).
The parallel analyses of five different tissues of the RT-PCR-positive pet rats revealed two to five RT-PCR-positive tissues (Supplementary Table S2). The SEOV-IgG ELISA (A) and hantavirus IgG ELISA (B) detected SEOV-reactive antibodies in all three rats of pet-rat owner 1 that originated from a vendor in North Rhine-Westphalia (Germany), but no antibodies were detected in the older rat (KS19/1414) obtained from a different source (Table 1). These results were confirmed by a line-blot assay (Supplementary Figure S2).
An analysis of the three rats from pet-rat owner 2 revealed that two of them tested positive for SEOV-reactive antibodies and SEOV RNA (Table 1 and Supplementary Figure S2).
A pairwise comparison of the partial S and L segment sequences revealed almost identical sequences for the SEOV-RNA-positive rats and the patient-derived sequences as well as a very high similarity to sequences from pet rats in the USA, Great Britain, France, and the Netherlands (Supplementary Table S3). The phylogenetic analysis of these partial L and S segment sequences confirmed the close phylogenetic relationship of all sequences from the German pet rats and the German patient with those from pet rats and patients in France, the Netherlands, Great Britain, and the USA (Figure 2a,b).

3.3. Complete Genome Identification and Phylogenetic Analysis

The complete coding sequences (CDSs) of two SEOV strains were determined via target enrichment-based HTS. The first strain originated from a pet rat that belonged to the patient/pet-rat owner 1 (KS19/1354; a partial sequence was already published [23]), and the second strain originated from a rat that belonged to pet-rat owner 2 (KS19/2195). The CDSs and amino acid (aa) sequences of both strains were almost identical for all three segments and the corresponding encoded proteins (Supplementary Table S4). The CDSs of the S, M, and L segments were 1290 nucleotides (nts), 3402 nts, and 6456 nts long (including the stop codon) and encoded proteins of 429 aa, 1133 aa, and 2151 aa residues, respectively.
A pairwise comparison of the S segment nt sequences of the rat of the patient/pet-rat owner 1 and the rat of pet-rat owner 2 confirmed their very high sequence similarity (99.7–100%) with sequences from breeder rats from the USA, Great Britain, France, and the Netherlands and from an SEOV patient in Great Britain (Supplementary Table S4). The S, M, and L segment sequences of both pet-rat-derived SEOV strains from Germany formed in the phylogenetic trees clusters with sequences from France, the Netherlands, Great Britain, and the USA (Figure 3a–c). These sequences were well separated in the phylogenetic trees from sequences from Belgium, China, Korea, and additional sequences from France and the USA. Similar results were obtained for partial M segment sequences of both pet rats deduced from the HTS dataset (Supplementary Figure S3).

3.4. Large-Scale Screening of Non-Pet Rats from Germany and The Netherlands

Six German scientific breeding facilities for Norway rats and black rats (panel II) were screened for SEOV RNA and SEOV-reactive antibodies. None of the investigated 404 Norway rats and 86 black rats showed any indication of an SEOV infection (Table 2).
In addition, 562 wild Norway rats and 1 wild black rat from urban and rural sites, livestock farms, and zoological gardens (pest rats) in Germany (panel III) were screened for SEOV RNA and SEOV-reactive antibodies. SEOV RNA was not detected in any of the investigated rats (Table 3). In the hantavirus IgG ELISA (B), a single rat from an agricultural farm in North Rhine-Westphalia tested positive, but only for the PUUV antigen in the line-blot assay (Table 3).
Of the 594 investigated rats from the Netherlands (panel IV), 6 tested positive in the IgG ELISA (C) (2013–2014: 1 R. rattus and 1 R. norvegicus; 2018: 4 R. norvegicus), all of which were negative in the real-time RT-PCR (Table 4). The four seropositive rats from 2018 were additionally tested via a virus neutralization assay using SEOV, DOBV, TULV, and PUUV that also turned out to be negative for SEOV and the other orthohantaviruses tested.

4. Discussion

This investigation of more than 1500 rats allowed for a broad overview of SEOV in pet, wild, and breeder rats in Germany and the Netherlands that was initiated after the first autochthonous human SEOV infection was discovered in Germany in 2019 [23]. In line with the detection of SEOV RNA in a patient’s rat [23], SEOV RNA was detected in two other pet rats of the same patient. The partial S and L segment SEOV sequences from the patient were almost identical to the novel rat-derived sequences detected here. In addition, all three pet rats tested positive for SEOV-reactive antibodies. The negative results in both SEOV RT-PCR and serology for a single older (female) pet rat in the household of the patient/pet-rat owner 1 suggested that the patient’s illness was caused by the recently purchased juvenile male rats. It might also indicate that the positive pet rats became infected at the breeding facility, whereas the patient’s older animal was refractory to infection due to the separate housing of male and female rats at the patient´s residence and possibly due to the lower infection rate in females [48].
Tracing back and forward activities resulted in the detection of additional SEOV-infected rats from pet-rat owner 2. The nucleotide sequences of strains from both the patient/pet-rat owner 1 and pet-rat owner 2 rats were almost identical, which suggested a common source of the virus strain. A phylogenetic analysis of the sequences revealed (nearly) identical SEOV sequences in breeder rats from France, the Netherlands and USA. This level of similarity was in contrast to previous findings concerning other human pathogenic hantaviruses; e.g., PUUV, for which a clear geographic clustering of sequences was observed that indicate virus evolution in isolated reservoirs of bank vole populations [49].
The high sequence similarity of the European pet and breeder rat SEOV sequences (even to sequences from breeder rats in the USA) suggested that there might be intensive private exchange of pet rats between neighboring European countries and even between Europe and the USA that is bypassing official trade channels. The trade of pet rats partly overlaps with the trade of feeder rats, as the history of the rats from pet-rat owners 1 and 2 showed, and was also reported previously [50]. However, the current study found no evidence of SEOV in breeder rats in Germany, which might be explained by the wild-rat origin of the vast majority of rats in these breedings, which are mainly dedicated to the phenotypic testing of rodenticide resistance. The finding in this study that no SEOV was found in wild rats from Germany or the Netherlands together with previously reported SEOV strains from wild rats in Europe that were divergent from strains from captive rats [19,20] suggested that the transmission of SEOV is maintained within the captive rat population itself.
None of the wild and breeder rats from Germany tested positive for SEOV RNA or anti-SEOV IgG antibodies except for a single rat from North Rhine-Westphalia that tested positive in the hantavirus IgG-ELISA (C) and for PUUV in the line-blot assay, which suggested a PUUV spillover infection when taking into account that North Rhine-Westphalia is an endemic region for PUUV [49]. Similarly, six rats from the Netherlands tested positive in the hantavirus IgG ELISA (B), but none of them tested positive in the real-time RT-PCR or in the virus neutralization test (used on 4/6 rats). The discrepancy between the positive IgG ELISA results and the lacking reactivity in the neutralization assay might have been caused by the lower sensitivity of the latter, in particular when a hantavirus was not used in the assay that caused the infection, or might have been due to a false-positive result in the ELISA [28]. The failed detection of viral RNA in seropositive rats might have been caused by virus-clearance-mediated absence of viral RNA. A false-negative result in the real-time RT-PCR cannot be excluded for SEOV strains circulating in the Netherlands because few nucleotide exchanges in the viral RNA might cause problems in primer or probe binding.
The SEOV strain previously detected in wild rats in the Netherlands [20] differed from the strains found in the pet rats from Germany, the Netherlands, and France but showed more similarity to the wild-rat-associated Humber strain in Great Britain. This finding might suggest a common origin of SEOV strains in wild rats; however, it remains speculative at this point whether the strain(s) in wild rats have their origin in pet or breeder rats. On the other hand, a potential wildlife origin of pet-rat-associated strain(s) in Germany, France, the Netherlands, Great Britain, and the USA still remains to be determined. Such One Health investigations should focus on monitoring of rats as well as the molecular characterization of the source of human HFRS cases worldwide. This rat monitoring could profit from the workflow developed here but might apply highly sensitive RT-PCR assays targeting different SEOV strains.

5. Conclusions

The results of this study clearly indicated that pet rats might pose a risk for SEOV transmission in Europe. This risk of zoonotic transmission is even higher because the contact between pet rats and humans is closer than between wild or feeder/breeder rats and humans. Therefore, locating the source of the SEOV infection and monitoring pet rats for infection with SEOV is important to prevent the spread of the virus and zoonotic transmission to humans. In addition, patients with suspected (SEOV) hantavirus infections should be routinely questioned about their previous contact with rats (specifically pet rats).

Supplementary Materials

The following Supplementary Materials are available online at https://www.mdpi.com/article/10.3390/v15020467/s1, Table S1: Screening strategies of the different rat panels; Table S2: Tissue distribution of SEOV RNA determined via conventional RT-PCR in pet rats from pet-rat owner 1/patient and pet-rat owner 2; Table S3: Pairwise partial nucleotide (nt) and amino acid (aa) sequence similarities of the two novel SEOV strains to the SEOV reference strains; Table S4: Pairwise nucleotide (nt) and amino acid (aa) sequence similarities of the complete coding sequences and the entire proteins of the two novel SEOV strains to the SEOV reference strains; Figure S1: Results of the trace-back and trace-forward investigations of the origin of the SEOV infection in the pet rats; Figure S2: Results of the line-blot assay of the rat sera from the patient (pet-rat owner 1) and pet-rat owner 2; Figure S3: Phylogenetic tree of partial M segment sequences of the pet rats of the patient and the other owner deduced from the HTS dataset.

Author Contributions

Conceptualization, R.G.U., A.L., J.D., J.H. and M.M.; methodology, M.M., A.L., G.H., J.T. and D.K.; validation, S.D., J.T. and M.M.; formal analysis, J.T., M.N. and M.P.d.C.; investigation, E.H., D.K., C.M., A.d.V., C.K., M.O., P.T., F.H., D.H., M.B. and M.N.; resources, J.T., D.K., R.P.-B., D.R., A.L. and G.H.; data curation, S.D., A.L., G.H. and M.M.; writing—original draft preparation, E.H., S.D., J.T., D.K., C.M., M.M., M.P.d.C., A.L., R.G.U. and A.d.V.; writing—review and editing, E.H., S.D., J.T., C.M., G.H., J.D., C.K., M.M., M.P.d.C., D.R. and R.G.U.; visualization, E.H., S.D. and M.P.d.C.; supervision, J.T., D.H.K., R.P.-B., M.M., A.L., J.H., G.H., J.D. and R.G.U.; project administration, M.M. and R.G.U.; funding acquisition, D.K., J.H. and R.G.U. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially supported by DZIF (TTU 01.808_00; “Emerging Infections”; awarded to R.G.U.). Additional funding was provided by the Federal Ministry of Education and Research through the Research Network Zoonotic Infectious Diseases (RoBoPub consortium; FKZ 01KI1721A awarded to R.G.U.; FKZ 01KI1721D awarded to J.D.) and the Federal Ministry of Education and Research through the National Research Platform for Zoonosis (FKZ 01KI1730 awarded to D.K.).

Institutional Review Board Statement

All animals were handled according to the applicable institutional, national, and international guidelines for the care and use of animals. Trappings of wild derived rats were conducted in accordance with the regular pest-control treatments in Germany according to German federal law (§ 18, Gesetz zur Verhütung und Bekämpfung von Infektionskrankheiten beim Menschen). The immunization of two Norway rats was performed in accordance with the general immunization program of the Friedrich-Loeffler-Institut (Landesamt für Landwirtschaft, Lebensmittelsicherheit und Fischerei, Mecklenburg-Vorpommern; permit: 28/17). Rats in the Netherlands were captured under permits from the Animal Ethics Committee (project numbers 201200208, 3 October 2012; 201400028, 26 February 2014, and 201500089, 10 September 2015) and from the Dutch Central Animal Experiments Committee (project number AVD3260020172104, 3 July 2017).

Informed Consent Statement

Not applicable.

Data Availability Statement

The relevant data are available upon reasonable request.

Acknowledgments

We thank Dörte Kaufmann for assistance with the rat dissections, Nicole Reimer for providing Figure 1, Bärbel Hammerschmidt and Sven Sander (all FLI) for the immunization of two rats, and Axel Blockesch (University of Applied Sciences, Frankfurt am Main) for support of the study. In addition, we acknowledge Jan Buijs and Huub Verlinden for rat collection in the Netherlands, Tryntsje Cuperus for the rat dissections, and Joke van der Giessen for project support. We also thank all collaborators who collected rats in Germany and provided them to us (as previously described in [23,24]); namely Erik Schmolz, Agnes Kalle, Anke Geduhn, Dania Richter, Franz-Rainer Matuschka, Jens van den Brandt, Jens Jacob, Uta Westerhüs, Christoph Langner, Lieselotte Minke, Johannes Lang, Nicole Kleemann, Wolfgang Fiedler, Henrike Gregersen, Sebastian Günther, Hermann Ansorge, Ulrich Riepe, Olaf Behlert, Eike Lange, René Driechciarz, Anita Plenge-Bönig, Thomas Scheibe, Jona Freise, Mechthild Budde, Anne Balkema Buschmann, Mario Heising, Richard Kruczewski, Marion Saathoff, Tessa Lohe, Lothar Stitz, Anja Globig, Odile Hecker, Tobias Eisenberg, Gabriele Ismer, Ute Ziegler, and Christin Trapp.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Krüger, D.H.; Figueiredo, L.T.M.; Song, J.-W.; Klempa, B. Hantaviruses--globally emerging pathogens. J. Clin. Virol. 2015, 64, 128–136. [Google Scholar] [CrossRef] [PubMed]
  2. Plyusnin, A.; Vapalahti, O.; Vaheri, A. Hantaviruses: Genome structure, expression and evolution. J. Gen. Virol. 1996, 77 Pt 11, 2677–2687. [Google Scholar] [CrossRef] [PubMed]
  3. Lee, H.W. Hemorrhagic Fever with Renal Syndrome in Korea. Clin. Infect. Dis. 1989, 11, S846–S876. [Google Scholar] [CrossRef]
  4. Clement, J.; LeDuc, J.W.; McElhinney, L.M.; Reynes, J.M.; Van Ranst, M.; Calisher, C.H. Clinical Characteristics of Ratborne Seoul Hantavirus Disease. Emerg. Infect. Dis. 2019, 25, 387–388. [Google Scholar] [CrossRef] [PubMed]
  5. Krüger, D.H.; Ulrich, R.G.; Hofmann, J. Hantaviruses as zoonotic pathogens in Germany. Dtsch. Arztebl. Int. 2013, 110, 461–467. [Google Scholar] [CrossRef]
  6. Lee, H.W.; Baek, L.J.; Johnson, K.M. Isolation of Hantaan virus, the etiologic agent of Korean hemorrhagic fever, from wild urban rats. J. Infect. Dis. 1982, 146, 638–644. [Google Scholar] [CrossRef]
  7. Hofmann, J.; Weiss, S.; Kuhns, M.; Zinke, A.; Heinsberger, H.; Krüger, D.H. Importation of Human Seoul Virus Infection to Germany from Indonesia. Emerg. Infect. Dis. 2018, 24, 1099–1102. [Google Scholar] [CrossRef]
  8. Zhang, Y.-Z.; Zou, Y.; Fu, Z.F.; Plyusnin, A. Hantavirus infections in humans and animals, China. Emerg. Infect. Dis. 2010, 16, 1195–1203. [Google Scholar] [CrossRef]
  9. Lin, X.-D.; Guo, W.-P.; Wang, W.; Zou, Y.; Hao, Z.-Y.; Zhou, D.-J.; Dong, X.; Qu, Y.-G.; Li, M.-H.; Tian, H.-F.; et al. Migration of Norway rats resulted in the worldwide distribution of Seoul hantavirus today. J. Virol. 2012, 86, 972–981. [Google Scholar] [CrossRef]
  10. Clement, J.; LeDuc, J.W.; Lloyd, G.; Reynes, J.-M.; McElhinney, L.; van Ranst, M.; Lee, H.-W. Wild Rats, Laboratory Rats, Pet Rats: Global Seoul Hantavirus Disease Revisited. Viruses 2019, 11, 652. [Google Scholar] [CrossRef]
  11. Castel, G.; Kant, R.; Badou, S.; Etougbétché, J.; Dossou, H.-J.; Gauthier, P.; Houéménou, G.; Smura, T.; Sironen, T.; Dobigny, G. Genetic Characterization of Seoul Virus in the Seaport of Cotonou, Benin. Emerg. Infect. Dis. 2021, 27, 2704–2706. [Google Scholar] [CrossRef]
  12. LeDuc, J.W.; Smith, G.A.; Childs, J.E.; Pinheiro, F.P.; Maiztegui, J.I.; Niklasson, B.; Antoniades, A.; Robinson, D.M.; Khin, M.; Shortridge, K.F.; et al. Global survey of antibody to Hantaan-related viruses among peridomestic rodents*. Bull. World Health Organ. 1986, 64, 139–144. [Google Scholar] [PubMed]
  13. Shepherd, J.G.; Blunsum, A.E.; Carmichael, S.; Smollett, K.; Maxwell-Scott, H.; Farmer, E.C.W.; Osborne, J.; MacLean, A.; Ashraf, S.; Shah, R.; et al. Seoul Virus Associated with Pet Rats, Scotland, UK, 2019. Emerg. Infect. Dis. 2021, 27, 2677–2680. [Google Scholar] [CrossRef]
  14. Reynes, J.-M.; Carli, D.; Bour, J.-B.; Boudjeltia, S.; Dewilde, A.; Gerbier, G.; Nussbaumer, T.; Jacomo, V.; Rapt, M.-P.; Rollin, P.E.; et al. Seoul Virus Infection in Humans, France, 2014–2016. Emerg. Infect. Dis. 2017, 23, 973–977. [Google Scholar] [CrossRef]
  15. McElhinney, L.; Fooks, A.R.; Featherstone, C.; Smith, R.; Morgan, D. Hantavirus (Seoul virus) in pet rats: A zoonotic viral threat. Vet. Rec. 2016, 178, 171–172. [Google Scholar] [CrossRef]
  16. Swanink, C.; Reimerink, J.; Gisolf, J.; de Vries, A.; Claassen, M.; Martens, L.; Waegemaekers, T.; Rozendaal, H.; Valkenburgh, S.; Hoornweg, T.; et al. Autochthonous Human Case of Seoul Virus Infection, the Netherlands. Emerg. Infect. Dis. 2018, 24, 2158–2163. [Google Scholar] [CrossRef]
  17. Lundkvist, Å.; Verner-Carlsson, J.; Plyusnina, A.; Forslund, L.; Feinstein, R.; Plyusnin, A. Pet rat harbouring Seoul hantavirus in Sweden, June 2013. Eurosurveillance 2013, 18, 20521. [Google Scholar] [CrossRef]
  18. Heyman, P.; Baert, K.; Plyusnina, A.; Cochez, C.; Lundkvist, A.; van Esbroeck, M.; Goossens, E.; Vandenvelde, C.; Plyusnin, A.; Stuyck, J. Serological and genetic evidence for the presence of Seoul hantavirus in Rattus norvegicus in Flanders, Belgium. Scand. J. Infect. Dis. 2009, 41, 51–56. [Google Scholar] [CrossRef]
  19. Murphy, E.G.; Williams, N.J.; Bennett, M.; Jennings, D.; Chantrey, J.; McElhinney, L.M. Detection of Seoul virus in wild brown rats (Rattus norvegicus) from pig farms in Northern England. Vet. Rec. 2019, 184, 525. [Google Scholar] [CrossRef] [PubMed]
  20. Verner-Carlsson, J.; Lõhmus, M.; Sundström, K.; Strand, T.M.; Verkerk, M.; Reusken, C.; Yoshimatsu, K.; Arikawa, J.; van de Goot, F.; Lundkvist, Å. First evidence of Seoul hantavirus in the wild rat population in the Netherlands. Infect. Ecol. Epidemiol. 2015, 5, 27215. [Google Scholar] [CrossRef] [PubMed]
  21. Meisel, H.; Wolbert, A.; Razanskiene, A.; Marg, A.; Kazaks, A.; Sasnauskas, K.; Pauli, G.; Ulrich, R.; Krüger, D.H. Development of novel immunoglobulin G (IgG), IgA, and IgM enzyme immunoassays based on recombinant Puumala and Dobrava hantavirus nucleocapsid proteins. Clin. Vaccine Immunol. 2006, 13, 1349–1357. [Google Scholar] [CrossRef] [PubMed]
  22. Ulrich, R.; Meisel, H.; Schütt, M.; Schmidt, J.; Kunz, A.; Klempa, B.; Niedrig, M.; Pauli, G.; Krüger, D.H.; Koch, J. Verbreitung von Hantavirusinfektionen in Deutschland. Bundesgesundheitsblatt Gesundh. Gesundh. 2004, 47, 661–670. [Google Scholar] [CrossRef] [PubMed]
  23. Hofmann, J.; Heuser, E.; Weiss, S.; Tenner, B.; Schoppmeyer, K.; Esser, J.; Klier, C.; Drewes, S.; Ulrich, R.G.; Krüger, D.H. Autochthonous Ratborne Seoul Virus Infection in Woman with Acute Kidney Injury. Emerg. Infect. Dis. 2020, 26, 3096–3099. [Google Scholar] [CrossRef] [PubMed]
  24. Fawzy, A.; Giel, A.-S.; Fenske, L.; Bach, A.; Herden, C.; Engel, K.; Heuser, E.; Boelhauve, M.; Ulrich, R.G.; Vogel, K.; et al. Development and validation of a triplex real-time qPCR for sensitive detection and quantification of major rat bite fever pathogen Streptobacillus moniliformis. J. Microbiol. Methods 2022, 199, 106525. [Google Scholar] [CrossRef] [PubMed]
  25. Raafat, D.; Mrochen, D.M.; Al’Sholui, F.; Heuser, E.; Ryll, R.; Pritchett-Corning, K.R.; Jacob, J.; Walther, B.; Matuschka, F.-R.; Richter, D.; et al. Molecular Epidemiology of Methicillin-Susceptible and Methicillin-Resistant Staphylococcus aureus in Wild, Captive and Laboratory Rats: Effect of Habitat on the Nasal S. aureus Population. Toxins 2020, 12, 80. [Google Scholar] [CrossRef]
  26. Maas, M.; de Vries, A.; Reusken, C.; Buijs, J.; Goris, M.; Hartskeerl, R.; Ahmed, A.; van Tulden, P.; Swart, A.; Pijnacker, R.; et al. Prevalence of Leptospira spp. and Seoul hantavirus in brown rats (Rattus norvegicus) in four regions in the Netherlands, 2011–2015. Infect. Ecol. Epidemiol. 2018, 8, 1490135. [Google Scholar] [CrossRef] [PubMed]
  27. Schmidt, J.; Jandrig, B.; Klempa, B.; Yoshimatsu, K.; Arikawa, J.; Meisel, H.; Niedrig, M.; Pitra, C.; Krüger, D.H.; Ulrich, R. Nucleocapsid protein of cell culture-adapted Seoul virus strain 80–39: Analysis of its encoding sequence, expression in yeast and immuno-reactivity. Virus Genes 2005, 30, 37–48. [Google Scholar] [CrossRef]
  28. Kramski, M.; Meisel, H.; Klempa, B.; Krüger, D.H.; Pauli, G.; Nitsche, A. Detection and typing of human pathogenic hantaviruses by real-time reverse transcription-PCR and pyrosequencing. Clin. Chem. 2007, 53, 1899–1905. [Google Scholar] [CrossRef]
  29. Hoornweg, T.E.; Zutt, I.; de Vries, A.; Maas, M.; Hoogerwerf, M.N.; Avšič-Županc, T.; Korva, M.; Reimerink, J.H.J.; Reusken, C.B.E.M. Development of a Comparative European Orthohantavirus Microneutralization Assay With Multi- Species Validation and Evaluation in a Human Diagnostic Cohort. Front. Cell. Infect. Microbiol. 2020, 10, 580478. [Google Scholar] [CrossRef]
  30. Mertens, M.; Wölfel, R.; Ullrich, K.; Yoshimatsu, K.; Blumhardt, J.; Römer, I.; Esser, J.; Schmidt-Chanasit, J.; Groschup, M.H.; Dobler, G.; et al. Seroepidemiological study in a Puumala virus outbreak area in South-East Germany. Med. Microbiol. Immunol. 2009, 198, 83–91. [Google Scholar] [CrossRef]
  31. Klempa, B.; Fichet-Calvet, E.; Lecompte, E.; Auste, B.; Aniskin, V.; Meisel, H.; Denys, C.; Koivogui, L.; ter Meulen, J.; Krüger, D.H. Hantavirus in African wood mouse, Guinea. Emerg. Infect. Dis. 2006, 12, 838–840. [Google Scholar] [CrossRef] [PubMed]
  32. Hall, T. BioEdit: An important software for molecular biology. GERF Bull. Biosci. 2011, 2, 60–61. [Google Scholar]
  33. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [PubMed]
  34. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef]
  35. Danecek, P.; Bonfield, J.K.; Liddle, J.; Marshall, J.; Ohan, V.; Pollard, M.O.; Whitwham, A.; Keane, T.; McCarthy, S.A.; Davies, R.M.; et al. Twelve years of SAMtools and BCFtools. Gigascience 2021, 10, giab008. [Google Scholar] [CrossRef]
  36. Robinson, J.T.; Thorvaldsdóttir, H.; Wenger, A.M.; Zehir, A.; Mesirov, J.P. Variant Review with the Integrative Genomics Viewer. Cancer Res. 2017, 77, e31–e34. [Google Scholar] [CrossRef]
  37. Garrison, E.; Marth, G. Haplotype-Based Variant Detection from Short-Read Sequencing. 2012. Available online: http://arxiv.org/pdf/1207.3907v2 (accessed on 17 August 2022).
  38. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  39. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. jModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef]
  40. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. In Proceedings of the 2010 Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; IEEE: New York, NY, USA, 2010; pp. 1–8, ISBN 978-1-4244-9751-5. [Google Scholar]
  41. Price, M.N.; Dehal, P.S.; Arkin, A.P. FastTree 2--approximately maximum-likelihood trees for large alignments. PLoS ONE 2010, 5, e9490. [Google Scholar] [CrossRef]
  42. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef]
  43. Ling, J.; Verner-Carlsson, J.; Eriksson, P.; Plyusnina, A.; Löhmus, M.; Järhult, J.D.; van de Goot, F.; Plyusnin, A.; Lundkvist, Å.; Sironen, T. Genetic analyses of Seoul hantavirus genome recovered from rats (Rattus norvegicus) in the Netherlands unveils diverse routes of spread into Europe. J. Med. Virol. 2019, 91, 724–730. [Google Scholar] [CrossRef] [PubMed]
  44. Lloyd, G.; Bowen, E.; Jones, N.; Pendry, A. Hfrs Outbreak Associated with Laboratory Rats In UK. The Lancet 1984, 323, 1175–1176. [Google Scholar] [CrossRef]
  45. Lloyd, G.; Jones, N. Infection of laboratory workers with hantavirus acquired from immunocytomas propagated in laboratory rats. J. Infect. 1986, 12, 117–125. [Google Scholar] [CrossRef]
  46. Shi, X.; McCaughey, C.; Elliott, R.M. Genetic characterisation of a Hantavirus isolated from a laboratory-acquired infection. J. Med. Virol. 2003, 71, 105–109. [Google Scholar] [CrossRef]
  47. Knust, B.; Brown, S.; de St Maurice, A.; Whitmer, S.; Koske, S.E.; Ervin, E.; Patel, K.; Graziano, J.; Morales-Betoulle, M.E.; House, J.; et al. Seoul Virus Infection and Spread in United States Home-Based Ratteries: Rat and Human Testing Results From a Multistate Outbreak Investigation. J. Infect. Dis. 2020, 222, 1311–1319. [Google Scholar] [CrossRef] [PubMed]
  48. Klein, S.L.; Marson, A.L.; Scott, A.L.; Ketner, G.; Glass, G.E. Neonatal sex steroids affect responses to Seoul virus infection in male but not female Norway rats. Brain Behav. Immun. 2002, 16, 736–746. [Google Scholar] [CrossRef]
  49. Drewes, S.; Ali, H.S.; Saxenhofer, M.; Rosenfeld, U.M.; Binder, F.; Cuypers, F.; Schlegel, M.; Röhrs, S.; Heckel, G.; Ulrich, R.G. Host-Associated Absence of Human Puumala Virus Infections in Northern and Eastern Germany. Emerg. Infect. Dis. 2017, 23, 83–86. [Google Scholar] [CrossRef] [PubMed]
  50. Cuperus, T.; de Vries, A.; Hoornweg, T.E.; Fonville, M.; Jaarsma, R.I.; Opsteegh, M.; Maas, M. Seoul Virus in Pet and Feeder Rats in The Netherlands. Viruses 2021, 13, 443. [Google Scholar] [CrossRef]
Figure 1. Location of Norway and black rat collection sites in Europe (a), Germany (wild-trapped rats and rats from breeding colonies) (b), and in the Netherlands (wild-trapped rats) (c). Map of the Netherlands (c): Amsterdam (blue); Rotterdam (yellow); Eindhoven (orange); Appeltern, Roermond, and Maastricht (green); and various municipalities of North Brabant and Limburg (purple).
Figure 1. Location of Norway and black rat collection sites in Europe (a), Germany (wild-trapped rats and rats from breeding colonies) (b), and in the Netherlands (wild-trapped rats) (c). Map of the Netherlands (c): Amsterdam (blue); Rotterdam (yellow); Eindhoven (orange); Appeltern, Roermond, and Maastricht (green); and various municipalities of North Brabant and Limburg (purple).
Viruses 15 00467 g001
Figure 2. Phylogenetic trees of partial Seoul virus (SEOV) sequences from patient and rats. (a) The partial L segment consensus tree was based on the HKY substitution model with gamma distribution for Bayesian analyses of 8 × 106 generations and the TPM1uf substitution model with gamma distribution and invariant sites for maximum likelihood analyses with an alignment length of 291 nucleotides and 1000 bootstrap replicates. (b) The partial S segment consensus tree was reconstructed using the TPM2 substitution model with gamma distribution for both the Bayesian (1 × 107 generations) and maximum likelihood analyses (1000 bootstrap replicates) using sequences of 515 nucleotides in length. The consensus trees were based on the Bayesian method. Bootstrap values from the maximum likelihood analyses were only transferred to the consensus tree if the branches were consistent. Posterior probabilities are given before and bootstrap values behind or below slashes. Relevant clades are highlighted by a light grey background. Sequences derived from the same outbreak are marked by a black bar, while novel sequences from the pet rats are indicated in darker grey. Country abbreviations: BEL, Belgium; BEN, Benin; CAM, Cambodia; CHN, China; FRA, France; GBR, Great Britain; GER, Germany; IDN, Indonesia; KOR, Korea; LAO, Laos; NED, the Netherlands; RUS, Russia; SGP, Singapore; SWE, Sweden; USA, United States of America; VNM, Vietnam. Other abbreviations: HU, human; PRO-1/P, pet-rat owner 1/patient; PRO-2, pet-rat owner 2; Rn, Rattus norvegicus; BR, breeder rat; CC, cell culture isolate; LR, laboratory rat; WR, wild rat. The categorization of the rats was according to [10,19,20,43,44,45,46,47].
Figure 2. Phylogenetic trees of partial Seoul virus (SEOV) sequences from patient and rats. (a) The partial L segment consensus tree was based on the HKY substitution model with gamma distribution for Bayesian analyses of 8 × 106 generations and the TPM1uf substitution model with gamma distribution and invariant sites for maximum likelihood analyses with an alignment length of 291 nucleotides and 1000 bootstrap replicates. (b) The partial S segment consensus tree was reconstructed using the TPM2 substitution model with gamma distribution for both the Bayesian (1 × 107 generations) and maximum likelihood analyses (1000 bootstrap replicates) using sequences of 515 nucleotides in length. The consensus trees were based on the Bayesian method. Bootstrap values from the maximum likelihood analyses were only transferred to the consensus tree if the branches were consistent. Posterior probabilities are given before and bootstrap values behind or below slashes. Relevant clades are highlighted by a light grey background. Sequences derived from the same outbreak are marked by a black bar, while novel sequences from the pet rats are indicated in darker grey. Country abbreviations: BEL, Belgium; BEN, Benin; CAM, Cambodia; CHN, China; FRA, France; GBR, Great Britain; GER, Germany; IDN, Indonesia; KOR, Korea; LAO, Laos; NED, the Netherlands; RUS, Russia; SGP, Singapore; SWE, Sweden; USA, United States of America; VNM, Vietnam. Other abbreviations: HU, human; PRO-1/P, pet-rat owner 1/patient; PRO-2, pet-rat owner 2; Rn, Rattus norvegicus; BR, breeder rat; CC, cell culture isolate; LR, laboratory rat; WR, wild rat. The categorization of the rats was according to [10,19,20,43,44,45,46,47].
Viruses 15 00467 g002
Figure 3. Phylogenetic trees of complete coding Seoul virus (SEOV) sequences of the S, M, and L segments. (a) The consensus tree of the complete coding S segment was derived from Bayesian analyses with 8 × 106 generations using the TIM1 substitution model with gamma distribution and invariant sites. The maximum likelihood reconstructions were performed using the TrN substitution model with gamma distribution and invariant sites on a sequence length of 1290 nucleotides and 1000 bootstrap replicates. (b) The complete coding M segment consensus tree was reconstructed using Bayesian algorithms with 8 × 106 generations and the TIM1 substitution model including gamma distribution and invariant sites. Maximum likelihood analyses used the GTR substitution model including gamma distribution and invariant sites based on a sequence length of 3402 nucleotides and 1000 bootstrap replicates. (c) The L segment complete coding consensus tree was reconstructed using the GTR substitution model with gamma distribution and invariant sites for Bayesian and maximum likelihood algorithms with 1 × 107 generations and 1000 bootstrap replicates, respectively, with an alignment length of 6456 nucleotides. Bootstrap values were only transferred from the maximum likelihood trees to establish consensus trees if the branches were identical. Posterior probabilities are given before and bootstrap values behind or below slashes. Relevant clades are highlighted in a light grey background. Sequences derived from the same outbreak are marked by a black bar, while novel sequences from the pet rats are indicated in a darker grey background. Country abbreviations: BEL, Belgium; BEN, Benin; CAM, Cambodia; CHN, China; FRA, France; GBR, Great Britain; GER, Germany; JPN, Japan; KOR, Korea; NED, the Netherlands; SGP, Singapore; USA, United States of America; VNM, Vietnam. Other abbreviations: HU, human; PRO-1/P, pet-rat owner 1/patient; PRO-2, pet-rat owner 2; Rn, Rattus norvegicus; BR, breeder rat; CC, cell culture isolate; LR, laboratory rat; WR, wild rat. The categorization of the rats was according to [10,19,20,43,44,45,46,47].
Figure 3. Phylogenetic trees of complete coding Seoul virus (SEOV) sequences of the S, M, and L segments. (a) The consensus tree of the complete coding S segment was derived from Bayesian analyses with 8 × 106 generations using the TIM1 substitution model with gamma distribution and invariant sites. The maximum likelihood reconstructions were performed using the TrN substitution model with gamma distribution and invariant sites on a sequence length of 1290 nucleotides and 1000 bootstrap replicates. (b) The complete coding M segment consensus tree was reconstructed using Bayesian algorithms with 8 × 106 generations and the TIM1 substitution model including gamma distribution and invariant sites. Maximum likelihood analyses used the GTR substitution model including gamma distribution and invariant sites based on a sequence length of 3402 nucleotides and 1000 bootstrap replicates. (c) The L segment complete coding consensus tree was reconstructed using the GTR substitution model with gamma distribution and invariant sites for Bayesian and maximum likelihood algorithms with 1 × 107 generations and 1000 bootstrap replicates, respectively, with an alignment length of 6456 nucleotides. Bootstrap values were only transferred from the maximum likelihood trees to establish consensus trees if the branches were identical. Posterior probabilities are given before and bootstrap values behind or below slashes. Relevant clades are highlighted in a light grey background. Sequences derived from the same outbreak are marked by a black bar, while novel sequences from the pet rats are indicated in a darker grey background. Country abbreviations: BEL, Belgium; BEN, Benin; CAM, Cambodia; CHN, China; FRA, France; GBR, Great Britain; GER, Germany; JPN, Japan; KOR, Korea; NED, the Netherlands; SGP, Singapore; USA, United States of America; VNM, Vietnam. Other abbreviations: HU, human; PRO-1/P, pet-rat owner 1/patient; PRO-2, pet-rat owner 2; Rn, Rattus norvegicus; BR, breeder rat; CC, cell culture isolate; LR, laboratory rat; WR, wild rat. The categorization of the rats was according to [10,19,20,43,44,45,46,47].
Viruses 15 00467 g003
Table 1. Molecular and serological analysis of rats from patient/pet-rat owner 1 and from pet-rat owner 2.
Table 1. Molecular and serological analysis of rats from patient/pet-rat owner 1 and from pet-rat owner 2.
OriginIDSexConventional RT-PCR *Serological InvestigationsHTSAccession Number
SEOV
S
SEOV
M
SEOV
L
IgG ELISA
SEOV (A)
Hantavirus
IgG ELISA (B)
Line-Blot **
Pet-rat owner 1KS19/1352m++++++n.d.OP870576 (S, par)
KS19/1353m++-+++n.d.n.a.
KS19/1354 ***m++++++spleen, liverOP870570 (S, CDS);
OP870572 (M, CDS); OP870574 (L, CDS)
KS19/1414f-n.d.----n.d.n.a.
Pet-rat owner 2KS19/2195m+n.d.n.d.+++spleen, liverOP870571 (S, CDS);
OP870573 (M, CDS); OP870575 (L, CDS)
KS19/2196m+n.d.n.d.+++n.d.n.a.
KS19/2197m-n.d.n.d.---n.d.n.a.
* Results of the conventional RT-PCR assays based on samples from lung, heart, liver, spleen, and kidney are given in Supplementary Table S2. ** Detailed line-blot reactivities of the samples are shown in Supplementary Figure S2. *** Initial detection of SEOV RNA from the first rat that caused the infection as previously described. Partial L and S segment sequences (accession numbers: MT123530 and MT123532, respectively) were already published [23]. Abbreviations: -, negative; +, positive; CDS, complete coding sequence; f, female; HTS, high-throughput sequencing; m, male; n.d., not done; par, partial sequence; n.a., not applicable.
Table 2. Results of serological and conventional RT-PCR analyses of different breeding colonies with wild rats or laboratory rats in Germany.
Table 2. Results of serological and conventional RT-PCR analyses of different breeding colonies with wild rats or laboratory rats in Germany.
Rat
Colony
NumberSpeciesHantavirus
IgG ELISA (B) **
Line-Blot
Assay
SEOV-S
RT-PCR
A39R. norvegicus0/390/10/39
33R. rattus0/330/10/33
C24R. norvegicus0/240/30/24
53R. rattus0/530/60/53
D *22R. norvegicus0/22n.d.0/22
E40R. norvegicus0/40n.d.0/40
F277R. norvegicus0/276n.d.0/277
G2R. norvegicus0/2n.d.0/2
Total490 0/4890/110/490
* Laboratory rats; ** used N proteins of DOBV, PUUV, HTNV, and SEOV; n.d., not done.
Table 3. Results of serological and conventional RT-PCR analyses of wild rats from different geographic regions in Germany.
Table 3. Results of serological and conventional RT-PCR analyses of wild rats from different geographic regions in Germany.
Region *NumberSpeciesHantavirus
IgG ELISA (B) **
Line-Blot
Assay
SEOV-S
RT-PCR
North45R. norvegicus0/330/40/45
West233R. norvegicus1 ***/2201 ***/160/228
1R. rattus0/1n.d.0/1
South121R. norvegicus0/1060/50/119
East + Central163R. norvegicus0/1500/90/152
Total563 1/5100/340/545
* North: Schleswig-Holstein and Hamburg; West: Lower Saxony, North Rhine-Westphalia, and Hesse; South: Baden-Wuerttemberg; East + Central: Mecklenburg-Western Pomerania, Brandenburg, Berlin, Saxony-Anhalt, and Saxony. ** Used N proteins of DOBV, PUUV, HTNV, and SEOV. *** Reactive in hantavirus IgG ELISA (B), but only reactive for PUUV in line-blot assay.
Table 4. Results of serological and real-time RT-PCR analyses of wildlife rats from different geographic regions in the Netherlands (Panel IV).
Table 4. Results of serological and real-time RT-PCR analyses of wildlife rats from different geographic regions in the Netherlands (Panel IV).
Study Location (s)Year of
Sample
Collection
Location TypeNumberSEOV-IgG ELISA (C)Real-Time RT-PCR
Amsterdam2013–2014Urban572/57 *0/2
Municipalities in the provinces of North-Brabant and Limburg2016Agricultural470/47n.d.
Appeltern, Roermond, and Maastricht2018Recreational water804/80 **0/80
Amsterdam2020Urban1360/1360/136
Eindhoven2020Urban50/50/5
Rotterdam2021Urban269n.d.0/269
Total 5946/3250/492
* One black rat and one Norway rat were seropositive and subsequently analyzed using RT-PCR. ** All four seropositive rats tested negative in the virus neutralization assay using SEOV, DOBV, PUUV, and TULV. Note: n.d., not done.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Heuser, E.; Drewes, S.; Trimpert, J.; Kunec, D.; Mehl, C.; de Cock, M.P.; de Vries, A.; Klier, C.; Oskamp, M.; Tenhaken, P.; et al. Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection. Viruses 2023, 15, 467. https://doi.org/10.3390/v15020467

AMA Style

Heuser E, Drewes S, Trimpert J, Kunec D, Mehl C, de Cock MP, de Vries A, Klier C, Oskamp M, Tenhaken P, et al. Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection. Viruses. 2023; 15(2):467. https://doi.org/10.3390/v15020467

Chicago/Turabian Style

Heuser, Elisa, Stephan Drewes, Jakob Trimpert, Dusan Kunec, Calvin Mehl, Marieke P. de Cock, Ankje de Vries, Christiane Klier, Martin Oskamp, Peter Tenhaken, and et al. 2023. "Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection" Viruses 15, no. 2: 467. https://doi.org/10.3390/v15020467

APA Style

Heuser, E., Drewes, S., Trimpert, J., Kunec, D., Mehl, C., de Cock, M. P., de Vries, A., Klier, C., Oskamp, M., Tenhaken, P., Hashemi, F., Heinz, D., Nascimento, M., Boelhauve, M., Petraityte-Burneikiene, R., Raafat, D., Maas, M., Krüger, D. H., Latz, A., ... Ulrich, R. G. (2023). Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection. Viruses, 15(2), 467. https://doi.org/10.3390/v15020467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop