Genetic Characterization of Hantaviruses Transmitted by the Korean Field Mouse (Apodemus peninsulae), Far East Russia

In an epizootiologic survey of 122 rodents captured in Vladivostok, Russia, antibodies positive for hantavirus were found in Apodemus peninsulae (4/70), A. agrarius (1/39), and Clethrionomys rufocanus (1/8). The hantavirus sequences identified in two seropositive A. peninsulae and two patients with hemorrhagic fever with renal syndrome (HFRS) from the Primorye region of Far East Russia were designated as Solovey and Primorye, respectively. The nucleotide sequences of the Solovey, Primorye, and Amur (obtained through GenBank) sequences were closely related (>92% identity). Solovey and Primorye sequences shared 84% nucleotide identity with the prototype Hantaan 76-118. Phylogenetic analysis also indicated a close relationship between Solovey, Primorye, Amur, and other viruses identified in Russia, China, and Korea. Our findings suggest that the Korean field mouse (A. peninsulae) is the reservoir for a hantavirus that causes HFRS over a vast area of east Asia, including Far East Russia.

(ANDV), and other related viruses cause HPS in the New World and are carried by the deer mouse (Peromyscus maniculatus), white-footed mouse (P. leucopus), cotton rat (Sigmodon hispidus), marsh rice rat (Oryzomys palustris), and Oligoryzomys longicaudatus, respectively (2,6). Although the known genotypes and serotypes have increased in number with advances in the knowledge of epidemiology and epizootiology of hantavirus infection (2), some still-unidentified hantaviruses carried by specific rodent hosts may exist. HFRS is generally known to be endemic to Far East Russia. However, the genetics of hantaviruses that are pathogenic for humans are not well defined. Reed voles (Microtus fortis) in Far East Russia were found to harbor two novel hantaviruses, Khabarovsk virus (KHAB) and Vladivostok virus (7,8). Another hantavirus, Topografov virus (TOPV), was isolated from brown lemmings (Lemmus sibiricus). The correlation between these three viruses and their pathogenicity for humans are not yet known (9).
A recent study reported two novel hantaviruses, designated as Amur (AMR) and Far East (FE), that were identified from HFRS patients in Far East Russia (10). The natural reservoir of AMR genotype seems to be A. peninsulae, according to a recent study on nucleotide sequence comparisons by Yashina et al. (11).
In 1999, we carried out an epizootiologic survey in a suburb of Vladivostok, Russia, to determine the characteristics of hantaviruses circulating in Far East Russia and to examine the possibility that A. peninsulae is a carrier of pathogenic hantaviruses. We detected antibodies to hantaviruses in A. penin-RESEARCH sulae, and the viral genome characteristics were extremely similar to the newly identified genotype, AMR (10). Using phylogenetic analysis to characterize the sequences of viruses identified from HFRS patients and A. peninsulae, we were able to corroborate the assumption of Yashina et al. (11). We also found that A. peninsulae-related viruses are pathogenic for humans and are distributed over a large area of east Asia that includes Far East Russia.

Materials and Methods
We collected sera and organs from wild rodents captured during 1999. We also collected sera and autopsy materials from HFRS patients in two rural villages in the Primorye region of Russia, located 400 km and 600 km from Vladivostok.
Rodent sera were screened for antibodies to HTNV and PUUV or both by indirect immunofluorescent antibody assay (IFA). Vero E6 cells infected with the Hantaan 76-118 strain of HTNV or the Sotkamo strain of PUUV were used as antigen slides. Diluted sera (1:16 and 1:64) were spotted onto the antigen slides and incubated at 37°C for 1 h. After three washes with phosphate-buffered saline (PBS), protein G-conjugated fluorescein isothiocyanate (FITC) (Zymed Laboratories, Inc., San Francisco, CA) was spotted onto the slides. After incubation at 37°C for 1 h, the slides were washed and observed by fluorescence microscopy. Scattered, granular fluorescence in the cytoplasm of infected Vero E6 cells was considered a positive reaction. Antibodies in HFRS patient sera were detected by the same protocol, except for the substitution of FITC-conjugated antihuman immunoglobulin (Ig) G (ICN Pharmaceuticals, Inc., Aurora, OH).
Total RNA was extracted from lung tissues of seropositive A. peninsulae with Isogen (Nippon Gene Co., Ltd., Osaka, Japan), which is based on the acid guanidium-phenol-chloroform technique, according to manufacturer's instructions. Similarly, total RNA was extracted from lung, liver, kidney, spleen, and brain tissues of HFRS patients. Reverse transcription (RT) was carried out at 42°C for 30 min by using Superscript II and random primer (Gibco-BRL, Rockville, MD). Full-length S segments were amplified with Platinum Taq (Gibco-BRL) and HTNV-full S primer for 30 polymerase chain reaction (PCR) cycles of denaturation at 94°C for 30 s, annealing at 55°C for 30 s, and extension at 68°C for 2 min. Amplification of M segments was identical to that of S segments, except for the use of M genome-specific primers ( Table  1). Part of the M segment (232 nucleotides) and the entire S segment (except for the 5´ and 3´ ends) were sequenced with primers specific for HTNV or SEOV or both. Amplification of the partial M segment was achieved only with nested PCR. The PCR-amplified products were separated by using a Rapid Gel Extraction kit (Gibco-BRL) according to the manufacturer's instructions. Purified DNA fragments were cloned into the PCR 2.1 vector provided in the TA cloning kit (Invitrogen Corporation, Carlsbad, CA). The ligated products were transformed into Top 10 competent cells (Invitrogen Corporation) and purified with a Miniprep kit (QIAGEN GmbH, Hilden, Germany). DNA sequencing was performed with the ABI-PRISM Dye Terminator Sequencing kit (Applied Biosystems, Foster City, CA) and an ABI 373-A genetic analyzer.
We used the ClustalX program package (version 1.81; available from: URL: ftp://ftp-igbmc.u-strasbg.fr/pub/Clust-alX/) to generate the phylogenetic trees by using the neighborjoining method with 1000 bootstrap replicates. Hantavirus sequences used in the comparisons were obtained from Gen-Bank. The S and M genome sequences used in this study are listed in Table 2.
Formalin-fixed lung, liver, kidney, and brain tissues from an HFRS patient who died of acute renal failure were observed under light microscopy and subjected to immunohistochemical analysis with monoclonal antibodies against Hantaan virus.

Results
We carried out the epizootiologic survey on 122 rodents captured in a suburb of Vladivostok; results of serologic screening of rodent sera by IFA are shown in Table 3. Identified rodent species included (70) A. peninsulae, (39) A. agrarius, (8) C. rufocanus, (3) M. fortis, and (2) Tamias sibiricus. Screening by IFA showed that one A. agrarius (2.5%), four A. peninsulae (5.7%), and one C. rufocanus (12.5%) had antibodies to HTNV or PUUV or both. HTNV-antibody titers ranged from 1:32 to 1:512. All the seropositive rodents, except for C. rufocanus, lacked antibody against PUUV (Table 4). Lung tis- We obtained the clinical histories of two fatal cases of HFRS in the Primorye region. The patients, who lived in villages 400 km and 600 km from Vladivostok, died 8-13 days after the onset of illness; gastrointestinal bleeding and acute renal failure were the causes of death. Serologic screening showed that both patients were positive for hantaviral antibodies. Antibody titers to HTNV and SEOV were apparently higher than to PUUV. We used lung, liver, kidney, spleen, and brain tissues of these HFRS patients for RT-PCR analysis; the lung and kidney tissues of patient no. 1 and the spleen tissue of patient no. 2 were positive for hantaviral M segment.
To examine the histopathologic changes in HFRS patients, we used light microscopy to examine sections of formalinfixed lung, liver, kidney, spleen, and brain tissues from patient no. 2, who had died of acute renal failure (Figure 1). The kidney was the only tissue that showed the recognizable histopathologic changes. Salient changes included interstitial edema with mild infiltration of mononuclear cells (Figure 1, small arrow) and degeneration of renal tubules (Figure 1, large arrow) in the cortex (Figure 1, A). Although proteinaceous casts and exudates were observed in the lumina of renal tubules (Figure 1, arrowhead), there were no apparent glomerular changes. In addition, a prominent well-defined necrotic lesion (Figure 1, asterisk) was noted in the medulla (Figure 1, B). Viral antigens were not detected in these specimens by using monoclonal antibodies to HTNV. The entire S segments of the viruses from two seropositive A. peninsulae were amplified and sequenced. We designated these segments as Solovey/AP61/1999 and Solovey/AP63/ 1999 based on the name of the village closest to the survey point, the rodent species from which the sample was taken, and the year in which the epizootiologic survey was done. We compared the coding regions of these sequences with those of other hantaviruses ( Table 5). The S segments of the two Solovey sequences had 99.0% and 98.8% identities in nucleotide and amino acid sequences, respectively. Solovey sequences and Hantaan viruses had 78.2%-84.5% nucleotide sequence identity and 86.7%-93.3% amino acid sequence identity, regardless of their source or geographical origin. Lower nucleotide sequence identities were seen than in Solovey sequences and other viruses: DOBV (73.6%), SEOV (73.9%), and SNV (63.9%).
To explore the genetic diversity of hantaviruses identified in A. peninsulae in more detail, we sequenced the partial M segment of the G2 region (232 nt). We also sequenced the partial M segments of genetic lineages identified in the two HFRS patients from the Primorye region, designated as Primorye/H1/ 2000 and Primorye/H2/2000. The M segment of Solovey and Primorye sequences were compared with those of other hantaviruses (Table 6). Nucleotide sequence identities among these sequences were between 92.2% and 98.2%; amino acid sequence identities were almost identical (98.7%-100%). We also compared the M segment sequences of Solovey and Primorye with those of AMR genetic lineage, recently identified in HFRS patients and A. peninsulae in Far East Russia (10,11    RESEARCH genetic lineage is widely distributed throughout east Asia. The M segments of Solovey, Primorye, and AMR sequences formed a common phylogenetic lineage with high bootstrap support values, regardless of viral origin (Figure 2,  A). Furthermore, H8205 shared a common lineage with Solovey and Primorye sequences. Another phylogenetic analysis, based on a different region of the M segment, showed that Chinese virus isolates (H8205, H3, H5, and B78) formed a distinct lineage within the Hantaan clade (Figure 2, B). The phylogenetic tree constructed for the S sequences ( Figure 3) showed that Solovey sequences formed a single cluster, together with Maaji1 (a Korean isolate) and B78, in a common lineage with high bootstrap support values within the Hantaan clade.
To identify signature amino acids for each virus type, we compared the deduced partial amino acid sequences of their G2 regions using ClustalX multiple-sequence alignment (Figure 4). The presence of leucine or isoleucine at amino acid position (aa) 903 was unique to HTNV except for AMR lineage. The signature amino acids for SEOV were leucine at aa 918 and valine, isoleucine, and serine at aa 955-957. The sig-  nature amino acids for AMR lineage were methionine at aa 932 and aspartic acid at aa 967.

Discussion
Each hantavirus serotype or genotype is generally associated with a specific rodent host, and various rodent species act as reservoir animals and sources of human infection. Since contact between rodents and humans occurs frequently during agricultural and forestry activities, most infections have been reported in rural areas. However, an urban epidemic of HFRS caused by SEOV has also been reported (26). A large number of rodent species may serve as reservoir animals for pathogenic hantaviruses. For example, few researchers suspected that P. maniculatus could transmit highly virulent hantavirus to humans until SNV was identified (27,28). Later studies showed that the other viral agents of HPS such as NYV, BCCV, BAYV, and ANDV, were carried by P. leucopus, S. hispidus (29), O. palustris (30), and O. longicaudatus (31), respectively. We emphasize the importance of discovering the characteristics of hantaviruses found in endemic areas and identifying the primary hosts.
Although Far East Russia has long been considered an HFRS-endemic area, few reports describe the hantaviral sequences in this region, and information on reservoir animals carrying pathogenic hantaviruses is limited. Our studies therefore focused on determining the genetic characteristics of hantaviruses circulating in this geographic area. We identified A. peninsulae as the natural reservoir rodent for a hantavirus pathogenic for humans in Far East Russia. We also identified hantavirus sequences designated as Solovey and Primorye in A. peninsulae and HFRS patients, respectively; genetic analysis showed that these sequences were very closely related to each other. This information and the pathological findings from the HFRS case in which Primorye sequence was identified strongly suggest that the virus of Solovey sequence is the causative agent of HFRS. The nucleotide sequence and phylogenetic analysis also showed that Solovey and Primorye sequences were most closely related to AMR and H8205 sequences from patients in Russia and China, but were clearly distinguishable from the prototype of Hantaan virus. Genetic and phylogenetic analysis indicated that Solovey and Primorye sequences were closely related to AMR, Maaji1, H8205, and B78 sequences, viruses derived from distant areas. While Solovey sequences were identified in a suburb of Vladivostok and PRI sequences in two villages 400 km and 600 km from Vladivostok, the H8205 and B78 viruses were derived in

RESEARCH
China, and Maaji1 was isolated in Korea. A. peninsulae is distributed in the same region where A. agrarius is prevalent in Korea (PW Lee, pers. comm.). Recently, AMR sequences were found in both HFRS patients and A. peninsulae (11). We suggest that some of the viruses circulating in the area of this study cause severe HFRS and are carried by the same host species, A. peninsulae. Comparison of the deduced hantaviral amino acid sequences showed that aspartic acid and methionine represented signature amino acids for AMR genetic lineage, regardless of the region in which the virus was identified or its origin (Figure 4). These signature amino acids may be used to distinguish AMR genetic lineage from other hantaviruses. We conclude from our results that A. peninsulae carries a hantavirus that is pathogenic for humans. Since A. peninsulae is widely distributed in Far East Russia, China, Korea, and Japan, this hantavirus and associated cases of HFRS may also be widely distributed.
In the kidney tissue of one HFRS patient (no. 2) from Primorye region, we detected pathologic changes typical of severe HFRS caused by hantavirus infection (32)(33)(34)(35). We also detected and sequenced the partial M segment in the spleen of the same patient. However, we could not detect the viral antigen in the kidney samples, possibly because of low levels of the virus in the kidneys of this patient. Nested PCR allowed the amplification of viral M segments from the spleen, but not from kidney, of this patient.
Through epizootiologic, clinical, pathologic, and sequencing studies, we identified a hantavirus carried by A. peninsulae as one of the causative agents of HFRS. We think that this information may be helpful in preventing human infections in East Asia. Controversy persists over whether A. peninsulae carries a distinct virus type or a subtype of HTNV. A similar question arises with Dobrava/Slovenia and Dobrava/Saaremaa, which are carried by A. flavicollis and A. agrarius, respectively. The S segment identities between Dobrava/Slovenia and Dobrava/Saaremaa (both obtained from GenBank for comparison purposes) were 87.8% (nucleotide) and 92.7% (amino acid). Similarly, the nucleotide and amino acid sequence identities of the S segments of Solovey sequences and HTN 76-118 were 82.7% and 92.2%, respectively. We suggest that Solovey sequences belong to a sublineage within the HTNV clade.