- Short report
- Open Access
Fugong virus, a novel hantavirus harbored by the small oriental vole (Eothenomys eleusis) in China
© Ge et al. 2016
- Received: 16 December 2015
- Accepted: 2 February 2016
- Published: 16 February 2016
The Erratum to this article has been published in Virology Journal 2016 13:75
Rodents are natural reservoirs of hantaviruses, which cause two disease types: hemorrhagic fever with renal syndrome in Eurasia and hantavirus pulmonary syndrome in North America. Hantaviruses related human cases have been observed throughout Asia, Europe, Africa, and North America. To date, 23 distinct species of hantaviruses, hosted by reservoir, have been identified. However, the diversity and number of hantaviruses are likely underestimated in China, and hantavirus species that cause disease in many regions, including Yunnan province, are unknown.
In August 2012, we collected tissue samples from 189 captured animals, including 15 species belonging to 10 genera, 5 families, and 4 orders in Fugong county, Yunnan province, China. Seven species were positive for hantavirus: Eothenomys eleusis (42/94), Apodemus peninsulae (3/25), Niviventer eha (3/27), Cryptotis montivaga (2/8), Anourosorex squamipes (1/1), Sorex araneus (1/1), and Mustela sibirica (1/2). We characterized one full-length genomic sequence of the virus (named fugong virus, FUGV) from a small oriental vole (Eothenomys eleusis). The full-length sequences of the small, medium, and large segments of FUGV were 1813, 3630, and 6531 nt, respectively. FUGV was most closely related to hantavirus LX309, a previously reported species detected in the red-backed vole in Luxi county, Yunnan province, China. However, the amino acid sequences of nucleocapsid (N), glycoprotein (G), and large protein (L) were highly divergent from those of Hantavirus LX309, with amino acid differences of 11.2, 15.3, and 12.7 %, respectively. In phylogenetic trees, FUGV clustered in the lineage corresponding to hantaviruses carried by rodents in the subfamily Arvicolinae.
High prevalence of hantavirus infection in small mammals was found in Fugong county, Yunnan province, China. A novel hantavirus species FUGV was identified from the small oriental vole. This virus is phylogenetic clustering with another hantavirus LX309, but shows highly genomic divergence.
- Eothenomys eleusis
- Hemorrhagic fever with renal syndrome
Hantaviruses (genus Hantavirus, family Bunyaviridae) cause two disease types: hemorrhagic fever with renal syndrome (HFRS) in Eurasia and hantavirus pulmonary syndrome (HPS) in North America . The hantavirus genome is composed of three segments, designated the large (L), medium (M), and small (S) segments, which encode an RNA-dependent RNA polymerase [2, 3], a glycoprotein precursor that is cotranslationally processed to produce two envelope glycoproteins (Gn and Gc), and a nucleocapsid (N) protein, respectively [3, 4]. HFRS is a typical rodent-borne disease and a severe public health issue in China . In addition to rodents, hantaviruses have recently been detected in shrews (order Soricomorpha, family Soricidae and Talpidae) [6–11] and bats (order Chiroptera, family Rhinolophidae, Hipposideridae, Vespertilionidae, and Nycteridae) [12–16], but it is not clear whether these hantaviruses lead to human illness. It is believed that only hantaviruses in rodent hosts (order Rodentia, family or subfamily Muridae, Arvicolinae, and Sigmodontinae) cause HFRS and HPS in Asia, Europe, Africa, and North America . To date, 23 distinct hantavirus species are recognized by the ICTV . In China, two main hantaviruses, Hantaan virus (HTNV) carried by the striped field mouse (Apodemus agrarius) and Seoul virus (SEOV) carried by the Norway rat (Rattus norvegicus), cause HFRS . In 2012, we reported a novel hantavirus harbored by the Yunnan red-backed vole (Eothenomys miletus) . In this study, we identified another novel hantavirus in the small oriental vole (Eothenomys eleusis, subfamily Arvicolinae), designated Fugong virus (FUGV) because it was recovered in Fugong county, Yunnan Province, China. The small oriental vole is distributed widely in Hunan, Hubei, Sichuan, Yunnan, and Guizhou, China .
In August 2012, 189 small mammals were captured in forests of the suburbs of Fugong county, Yunnan province. Animal species were identified based on morphology, and animal weight and sex were recorded. Species were further identified by DNA sequencing of the mitochondrial cytochrome b (CytB) gene following previously described methods . Animal lung tissues were collected for further analysis.
Direct immunofluorescence assay (IFA) detection of hantavirus
Hantavirus antigens present in the lung tissue samples were detected by direct immunofluorescence assay (IFA) as previously described . Briefly, animal lung tissue samples were fixed and cut into 4-μm sections using a Cryocut microtome (KD-1508 Frigocut; Zhejiang, China) at −25 °C and fixed in cold acetone. Sections were stained with rabbit anti-SEOV/Z37 or HTNV/Z10 antibodies labeled with fluorescein isothiocyanate (FITC) (Zhejiang CDC) . Scattered, granular fluorescence in the cytoplasm was considered positive staining.
DNA and RNA extraction and reverse transcription PCR
Total DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) from tissue samples according to the manufacturer’s protocol. Viral RNA was extracted from lung tissues using the QIAamp Viral RNA Mini Kit (Qiagen) following the manufacturer’s instructions, and cDNA was synthesized with M-MLV reverse transcriptase (Promega, Madison, WI, USA) in the presence of the primer P14 .
Tissues from positive samples were screened by reverse transcription PCR (RT-PCR) using degenerate primers based on the conserved domain of the L fragment of Hantavirus genomes . Standard precautions were taken to avoid PCR contamination, and no false-positives were observed in negative controls. PCR products were gel purified and sequenced with both forward and reverse primers using the 3100 Sequencer (ABI, Waltham, MA, USA).
One positive sample of E. eleusis was chosen for hantavirus genomic sequencing. The complete genomic sequence was amplified by PCR using primers designed from published hantavirus sequences or from sequences obtained in this study (available upon request). PCR products were also gel purified and sequenced with both forward and reverse primers using the 3100 Sequencer (ABI). The sequencing chromatograms were inspected carefully for overlapping multicolor peaks, which are an indicator of sequence heterogeneity in the amplicons. The PCR products of these samples were cloned into the pGEM-T Easy Vector (Promega) and at least 5 clones for each PCR fragment were sequenced to obtain a consensus sequence.
Preliminary sequence management and analysis were carried out using Geneious (Version 5.5.9, Biomatters Limited, Aukland, New Zealand) and sequence alignment and editing were performed using ClustalW (Version 2.0)  and BioEdit (Version 7.1.9) . The potential in vitro recombinant sequences (i.e., PCR artifacts) were screened and discarded using Recombination Detection Program v2.0 . The consensus sequences were compared with known hantavirus sequences available in GenBank. Phylogenetic trees were constructed using the maximum likelihood (ML) algorithm with bootstrap values determined by 1000 replicates in Geneious.
Detection of hantaviruses in several animals
Detection of hantavirus in small mammals in Fugong County, Yunnan Province, China in 2012
Virus Counts (Positive/Tested)
Rodentia, Cricetidae, Eothenomys
Rodentia, Muridae, Apodemus
Rodentia, Muridae, Niviventer
Rodentia, Muridae, Vernaya
Soricomorpha, Soricidae, Crocidura
Soricomorpha, Soricidae, Cryptotis
Soricomorpha, Soricidae, Anourosorex
Soricomorpha, Soricidae, Sorex
Lagomorpha, Ochotomidae, Ochotona
Carnivora, Mustelidae, Mustela
Hantavirus RNA was detected in 21 of 42 IFA-positive samples of E. eleusis. PCR amplification of the hantavirus sequence from samples obtained from 6 additional species failed. We sequenced all of the amplified fragments and found that the 420-nt (nucleotide) sequences were highly similar, with >99 % nt sequence identities. The hantavirus was named Fugong virus (FUGV).
Genetic characterization and analysis
Nucleotide and amino acid sequence similarity (%) between FUGV virus and representative hantaviruses (23 ICTV-nominated species and two closely related viruses LX309 and YN06-862)
nt and aa identities (%)
Virus species (strain)
Andes virus Chile-9717869
Bayou virus F0260003
Black Creek Canal virus
Cano Delgadito virus VHV-574
Dobrava-Belgrade virus Afl9/1999
El Moro Canyon virus RM-97
Hantaan virus 76–118
Isla Vista virus MC-SB-1
Khabarovsk virus Fuyuan-Mm-312
Laguna Negra virus 510B
Muleshoe virus SH-Tx-339
New York virus NY-2
Prospect Hill virus PH-1
Puumala virus Sotkamo
Rio Mamore virus HTN-007
Rio Segundo virus RMx-Costa-1
Saaremaa virus Saaremaa-160v
Sangassou virus SA14
Seoul virus 80–39
Sin Nombre virus NM H10
Thailand virus Anjozorobe
Thottapalayam virus VRC 66412
Topografov virus Ls136V
Tula virus Moravia/5302v/95
The full-length L-genomic segment was 6531 nt and encoded an RNA-dependent RNA polymerase of 2152 aa; it contained NCRs of 37 and 35 nt at the 5'- and 3'-ends, respectively. Previous studies have identified five conserved motifs (A, B, C, D, and E) in hantavirus RNA polymerases [31, 32]. These motifs were all detected in the FUGV L protein. In addition, other typical conserved motifs among hantavirus RNA polymerases, like premotif A and the acidic C-terminal domain , were found (data not shown).
A pairwise alignment and comparison with known hantavirus genomes (from previously detected species or strains) showed low nt sequence similarity in the S, M, and L segments, ranging from 41.5 to 63.1 %, 51.1 to 75.5 %, and 61.9 to 76.3 %, respectively (Table 2). Deduced aa similarities were higher (45.6–88.8 %, 41.5–84.7 %, and 61.5–87.3 % for the N, Gn/Gc, and L proteins, respectively). Additionally, all of the FUGV segments were most closely related to Hantavirus LX309 sequences, which is a putative novel hantavirus species found in the red-backed vole in Luxi County, Yunnan, China . However, the aa sequences were highly divergent between the N, G, and L proteins of FUGV and Hantavirus LX309, with differences reaching 11.2, 15.3, and 12.7 %, respectively.
To infer the evolutionary relationships for FUGV and its host, an ML tree was built using mitochondrial CytB of the Rodentia hosts of representative hantaviruses, and Soricomorpha hosts were used as an outgroup (Fig. 2c). The tree of the hosts was in agreement with that of previous studies; the Murinae subfamily and Cricetidae family formed two monophyletic groups, with Soricomorpha as the outgroup. The Cricetidae were further subdivided into the subfamilies Arvicolinae, Neotominae, and Sigmodontinae . To understand virus-host co-divergence, the L tree of viruses were compared with the CytB tree of hosts. We found that FUGV and LX309 diverged earlier than Prospect Hill virus, Tula virus, Puumala virus, and others, which was not congruent with their Myodes spp. host evolutionary statuses in the CytB tree. This phenomenon may be explained by differences in the evolutionary rates of hantaviruses in different hosts. In the tree, the same phenomenon can also can be observed for Sangassou virus and Hylomyscus simus (African wood mouse), and for Dobrava-Belgrade virus and Apodemus flavicollis (yellow-necked mouse) [33, 34].
Hantavirus antigens were detected by direct IFA in 7 out of 15 total species collected in Fugong county, Yunnan Province. We successfully amplified the nucleotide sequences of the virus obtained from E. eleusis. These results suggested a high cross-reactivity of SEOV and HTNV with other hantaviruses in serological level, but their genomic sequences were highly divergent. We thus suspected that there were other novel hantaviruses that cannot be detected with the degenerate primers used in this study. We also cannot exclude potential RT-PCR sensitivity issues.
The three segments of FUGV harbored by E. eleusis were most closely related to Hantavirus LX309 carried by E. miletus. However, the aa sequences of the N, G, and L proteins of FUGV were highly divergent with respect to Hantavirus LX309, with similarities reaching 11.2, 15.3, and 12.7 %, respectively. The new ICTV classification guidelines indicate that a novel Hantavirus species should show at least a 10 % difference in S segment similarity and a 12 % difference in M segment similarity based on complete amino acid sequences . Accordingly, the hantavirus harbored by E. eleusis represents a novel species. Additionally, no reassortment or recombination events were detected among FUGV and other hantaviruses. We propose that this novel species be named Fugong virus (FUGV) based on the location where it was first recovered.
In phylogenetic analyses of the three segments of the FUGV genome, it grouped with hantaviruses harbored by rodents in the subfamily Arvicolinae, which can cause HFRS in Eurasia. For example, the Puumala virus (PUUV) harbored by Arvicolinae (Myodes glareolus) causes most HFRS cases, with 35,424 reported cases in 2006 [36, 37]. Recently, more than 100 HFRS cases have been reported in Yunnan province (unpublished data). However, the types of hantaviruses that cause HFRS in Yunnan province are still unknown. It is not clear whether FUGV can infect humans and lead to illness.
Phylogenetic analyses of hantaviruses and their rodent hosts suggested that the viruses have a long history of co-evolution with their predominant rodent carriers [38, 39]. FUGV was most closely related to hantavirus LX309, which is carried by different rodent species within the same genus. Additionally, we detected differences in the evolutionary rates of hantaviruses in different hosts. The Eothenomys genus consists of eight rodent species, but only two have been investigated. Considering the high prevalence and diversity of hantaviruses in the genus, additional investigations should be performed in the future to extend the number of species considered.
This work was jointly funded by the National Natural Science Foundation of China (81260437; 81060132; 81290341), the Project of State Key Laboratory for Infectious Diseases Prevention and Control (2013SKLID302), the Scientific and Technological Basis Special Project from the Minister of Science and Technology of the People’s Republic of China (2013FY113500), and the China Mega-Project for Infectious Disease (2014ZX10004001-003) from the Minister of Science and Technology of the People’s Republic of China.
This study was reviewed and approved by the Ethics Committee of the Yunnan Institute of Endemic Disease Control and Prevention. All animals were treated in strict according to the guidelines for the Laboratory Animal Use and Care from the Chinese CDC and the Rules for the Implementation of Laboratory Animal Medicine (1998) from the Ministry of Health, China, under the protocols approved by the National Institute for Communicable Disease Control and Prevention. All surgery was performed under ether anesthesia, and all efforts were made to minimize suffering.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Schmaljohn C, Hjelle B. Hantaviruses: a global disease problem. Emerg Infect Dis. 1997;3:95–104.View ArticlePubMedPubMed CentralGoogle Scholar
- Kukkonen SK, Vaheri A, Plyusnin A. Completion of the Tula hantavirus genome sequence: properties of the L segment and heterogeneity found in the 3’ termini of S and L genome RNAs. J Gen Virol. 1998;79:2615–22.View ArticlePubMedGoogle Scholar
- Schmaljohn CS, Hasty SE, Harrison SA, Dalrymple JM. Characterization of Hantaan virions, the prototype virus of hemorrhagic fever with renal syndrome. J Infect Dis. 1983;148:1005–12.View ArticlePubMedGoogle Scholar
- Plyusnin A, Vapalahti O, Vaheri A. Hantaviruses: genome structure, expression and evolution. J Gen Virol. 1996;77:2677–87.View ArticlePubMedGoogle Scholar
- Huang X, Yin H, Yan L, Wang X, Wang S. Epidemiologic characteristics of haemorrhagic fever with renal syndrome in Mainland China from 2006 to 2010. Western Pac Surveill Response J. 2012;3:12–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Arai S, Ohdachi SD, Asakawa M, Kang HJ, Mocz G, et al. Molecular phylogeny of a newfound hantavirus in the Japanese shrew mole (Urotrichus talpoides). Proc Natl Acad Sci U S A. 2008;105:16296–301.View ArticlePubMedPubMed CentralGoogle Scholar
- Gu SH, Kang HJ, Baek LJ, Noh JY, Kim HC, et al. Genetic diversity of Imjin virus in the Ussuri white-toothed shrew (Crocidura lasiura) in the Republic of Korea, 2004–2010. Virol J. 2011;8:56.View ArticlePubMedPubMed CentralGoogle Scholar
- Klempa B, Fichet-Calvet E, Lecompte E, Auste B, Aniskin V, et al. Novel hantavirus sequences in Shrew, Guinea. Emerg Infect Dis. 2007;13:520–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Radosa L, Schlegel M, Gebauer P, Ansorge H, Heroldova M, et al. Detection of shrew-borne hantavirus in Eurasian pygmy shrew (Sorex minutus) in Central Europe. Infect Genet Evol. 2013;19:403–10.View ArticlePubMedGoogle Scholar
- Song JW, Kang HJ, Gu SH, Moon SS, Bennett SN, et al. Characterization of Imjin virus, a newly isolated hantavirus from the Ussuri white-toothed shrew (Crocidura lasiura). J Virol. 2009;83:6184–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Zuo SQ, Gong ZD, Fang LQ, Jiang JF, Zhang JS, et al. A new hantavirus from the stripe-backed shrew (Sorex cylindricauda) in the People’s Republic of China. Virus Res. 2014;184:82–6.View ArticlePubMedGoogle Scholar
- Weiss S, Witkowski PT, Auste B, Nowak K, Weber N, et al. Hantavirus in bat, Sierra Leone. Emerg Infect Dis. 2012;18:159–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Sumibcay L, Kadjo B, Gu SH, Kang HJ, Lim BK, et al. Divergent lineage of a novel hantavirus in the banana pipistrelle (Neoromicia nanus) in Cote d’Ivoire. Virol J. 2012;9:34.View ArticlePubMedPubMed CentralGoogle Scholar
- Guo WP, Lin XD, Wang W, Tian JH, Cong ML, et al. Phylogeny and origins of hantaviruses harbored by bats, insectivores, and rodents. PLoS Pathog. 2013;9:e1003159.View ArticlePubMedPubMed CentralGoogle Scholar
- Gu SH, Lim BK, Kadjo B, Arai S, Kim JA, et al. Molecular phylogeny of hantaviruses harbored by insectivorous bats in Cote d’Ivoire and Vietnam. Viruses. 2014;6:1897–910.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu L, Wu J, He B, Qin S, Xia L, et al. Novel hantavirus identified in black-bearded tomb bats, China. Infect Genet Evol. 2015;31:158–60.View ArticlePubMedGoogle Scholar
- Jonsson CB, Figueiredo LT, Vapalahti O. A global perspective on hantavirus ecology, epidemiology, and disease. Clin Microbiol Rev. 2010;23:412–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Plyusnin A, Beatty BJ, Elliott RM, Goldbach R, Kormelink R, et al. In: King AMQ, Lefkowitz EJ, Adams MJ, Carstens EB, editors. Ninth report of the International Committee on Taxonomy of Viruses. Bunyaviridae. San Diego, CA, USA: Elsevier Academic Press; 2012. p. 725–41.Google Scholar
- Zhang Y, Zhang H, Dong X, Yuan J, Zhang H, et al. Hantavirus outbreak associated with laboratory rats in Yunnan, China. Infect Genet Evol. 2010;10:638–44.View ArticlePubMedGoogle Scholar
- Zhang Y, Yuan J, Yang X, Zhou J, Yang W, et al. A novel hantavirus detected in Yunnan red-backed vole (Eothenomys miletus) in China. J Gen Virol. 2011;92:1454–7.View ArticlePubMedGoogle Scholar
- Luo J, Yang D, Suzuki H, Wang Y, Chen WJ, et al. Molecular phylogeny and biogeography of Oriental voles: genus Eothenomys (Muridae, Mammalia). Mol Phylogenet Evol. 2004;33:349–62.View ArticlePubMedGoogle Scholar
- Ge XY, Li JL, Yang XL, Chmura AA, Zhu G, et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature. 2013;503:535–8.View ArticlePubMedGoogle Scholar
- Wang H, Yoshimatsu K, Ebihara H, Ogino M, Araki K, et al. Genetic diversity of hantaviruses isolated in china and characterization of novel hantaviruses isolated from Niviventer confucianus and Rattus rattus. Virology. 2000;278:332–45.View ArticlePubMedGoogle Scholar
- Jiang JF, Zuo SQ, Zhang WY, Wu XM, Tang F, et al. Prevalence and genetic diversities of hantaviruses in rodents in Beijing, China. Am J Trop Med Hyg. 2008;78:98–105.PubMedGoogle Scholar
- Schmaljohn CS, Jennings GB, Hay J, Dalrymple JM. Coding strategy of the S genome segment of Hantaan virus. Virology. 1986;155:633–43.View ArticlePubMedGoogle Scholar
- Klempa B, Fichet-Calvet E, Lecompte E, Auste B, Aniskin V, et al. Hantavirus in African wood mouse, Guinea. Emerg Infect Dis. 2006;12:838–40.View ArticlePubMedPubMed CentralGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23:2947–8.View ArticlePubMedGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Martin DP, Williamson C, Posada D. RDP2: recombination detection and analysis from sequence alignments. Bioinformatics. 2005;21:260–2.View ArticlePubMedGoogle Scholar
- Lober C, Anheier B, Lindow S, Klenk HD, Feldmann H. The Hantaan virus glycoprotein precursor is cleaved at the conserved pentapeptide WAASA. Virology. 2001;289:224–9.View ArticlePubMedGoogle Scholar
- Kukkonen SK, Vaheri A, Plyusnin A. L protein, the RNA-dependent RNA polymerase of hantaviruses. Arch Virol. 2005;150:533–56.View ArticlePubMedGoogle Scholar
- Yadav PD, Vincent MJ, Nichol ST. Thottapalayam virus is genetically distant to the rodent-borne hantaviruses, consistent with its isolation from the Asian house shrew (Suncus murinus). Virol J. 2007;4:80.View ArticlePubMedPubMed CentralGoogle Scholar
- Klempa B, Witkowski PT, Popugaeva E, Auste B, Koivogui L, et al. Sangassou virus, the first hantavirus isolate from Africa, displays genetic and functional properties distinct from those of other murinae-associated hantaviruses. J Virol. 2012;86:3819–27.View ArticlePubMedPubMed CentralGoogle Scholar
- Nemirov K, Vapalahti O, Papa A, Plyusnina A, Lundkvist A, et al. Genetic characterization of new Dobrava hantavirus isolate from Greece. J Med Virol. 2003;69:408–16.View ArticlePubMedGoogle Scholar
- Maes P, Klempa B, Clement J, Matthijnssens J, Gajdusek DC, et al. A proposal for new criteria for the classification of hantaviruses, based on S and M segment protein sequences. Infect Genet Evol. 2009;9:813–20.View ArticlePubMedGoogle Scholar
- Heyman P, Vaheri A, ENIVD Members. Situation of hantavirus infections and haemorrhagic fever with renal syndrome in European countries as of December 2006. Euro Surveill. 2008;13. http://www.eurosurveillance.org/ViewArticle.aspx?ArticleId=18925
- Heyman P, Vaheri A, Lundkvist A, Avsic-Zupanc T. Hantavirus infections in Europe: from virus carriers to a major public-health problem. Expert Rev Anti Infect Ther. 2009;7:205–17.View ArticlePubMedGoogle Scholar
- Plyusnin A, Morzunov SP. Virus evolution and genetic diversity of hantaviruses and their rodent hosts. Curr Top Microbiol Immunol. 2001;256:47–75.PubMedGoogle Scholar
- Easterbrook JD, Klein SL. Immunological mechanisms mediating hantavirus persistence in rodent reservoirs. PLoS Pathog. 2008;4:e1000172.View ArticlePubMedPubMed CentralGoogle Scholar