Molecular evidence of Tula virus in Microtus obscurus in the region of Yili, Xinjiang, China
BMC Infectious Diseases volume 19, Article number: 527 (2019)
Hantaviruses are important zoonotic pathogens, and they pose a profound risk to public health. So far, there has been no evidence showing that Tula virus (TULV), one species of hantavirus, is endemic in China. In this study, we captured rodents and found that the Tula virus had infected voles in Yili region, Xinjiang, China.
Rodents were captured by flooding their burrows in mountain pasture areas in Narati, Xinyuan County, Xinjiang, China. Hantavirus L gene fragments were amplified by nest RT-PCR using genus-specific primers. Positive samples were further identified by sequencing of RT-PCR products of S gene fragment for species identification. To identify the species of captured small mammals, the rodents’ cytochrome b (Cytb) was amplified by PCR and sequenced. Phylogenetic analysis was used to show the clustering and evolution relationship of the viral nucleic acids.
Here, 31 out of 198 voles captured (16%) were infected with TULV. Host sequencing analysis showed these voles were Microtus obscurus (M. obscurs). Alignment and phylogenetic analysis of the exon region (1191 bp) of the hantavirus S gene confirmed that all of the detected amplicons were TULV, which was similar to one strain of TULV identified in Kazakhstan.
This is the first identification of Tula virus in China, and we found that M. obscurus acts as a natural reservoir for carrying the virus. Although the infection rate in the local human population remains unknown, the high prevalence of TULV in the small mammals in the region constitutes a risk that this putative pathogen may spread to the local population.
Tula virus (TULV) is a Puumala-like virus, belongs to the genus hantavirus of family Bunyaviridae, and is primarily carried by Microtus arvalis . TULV was initially discovered in M. arvarlis in Tula, Central Russia in the early 1990s . Since then, M. arvarlis (obscurus ) has been found to act as a reservoir host of TULV in many European countries [4, 5], West Siberia , and Eastern Kazakhstan . A few reports have shown human cases presenting a fever syndrome with severe headache  or hemorrhagic fever with renal syndrome and pulmonary involvement [9, 10]. However, most individuals with serological and/or pathogenic positive test were asymptomatic in European countries . So far, no study has shown the epidemiological situation of TULV in China, even though voles are common rodents in some pasture areas in northern Xinjiang, China.
Sample collection and preparation
The investigation was performed in September 2015 in Narati mountain pastures in Xinyuan County, Yili region, Xinjiang. Voles were captured by flooding their burrows with a water pump. The voles that were captured alive were euthanized with barbiturate (100 mg/ kg). The voles were initially identified morphologically, then autopsied and approximately 1 g of lung tissue from each specimen was isolated and cut into 0.5-cm-thick slices using a pair of aseptic operation scissors and then placed to a sterile tube containing 5 mL of RNA stabilization reagent (Ambion RNAlater®, Life Technologies, Carlsbad, CA, US). Each tube was stored at 5 °C overnight, then transferred to − 80 °C for long-term storage. The capture of rodents in fields and protocols for using animals were approved by the Ethics Committee of the First Affiliated Hospital of Xinjiang Medical University, Urumqi, China (approval IACUC-2015).
RNA extraction, cDNA synthesis, and PCR
Total RNA was extracted from the lung tissues of voles using TRIzol reagent (Thermo Fisher) according to the manufacturer’s instructions. The reverse transcription primer P14 (Table 1) was used for the first strand cDNA synthesis with the GoScript Reverse Transcription System (Promega, Beijing, China) [12, 16]. In brief, A PCR tube containing 2.0 μL of extracted RNA (500 ng/μL) and 1.0 μL of primer P14 (10 pmol/l) was heated in a 70 °C water bath for 5 min, then chilled in ice water for 5 min. After brief centrifugation, 2.0 μL of GoScript 5 × Reverse Transcriptase buffer, 0.5 μL of GoScript Reverse Transcriptase (200 U/μL), 0.5 μL of Recombinant RNasin Ribonuclease inhibitor (40 U/μL), 1 μL of PCR Nucleotide Mix (10 mM), 2.0 μL of MgCl2 (25 mM), and 1.0 μL of RNase-free deionized water were added to the tube. After incubation at 25 °C for 5 min, the tube was placed in a 42 °C water bath for 60 min. Reverse transcriptase was inactivated by heating to 70 °C for 15 min.
To identify hantavirus, the L gene segment (438 bp) was amplified by nest-PCR with viral-genus-specific primers, including outer PCR primers HAN-L-F1 and HAN-L-R1 and inner PCR primers HAN-L-F2 and HAN-L-R2 (Table 1) [13, 17]. All of the PCR products were analyzed using 2% agarose gel. To determine the species of the virus, S gene fragments were amplified by nest-PCR with cDNAs from L gene positive samples using reported primers in Table 1 .
DNA extraction and molecular identification for species of hosts
Genomic DNA was extracted from approximately 30 mg of lungs of small mammals infected with hantavirus using a Tiangen TIANamp Genomic DNA Kit (Tiangen Biotech, Beijing, China) according to the manufacturer’s instructions. Complete cytochrome b (Cytb) gene as a phylogenetic marker was amplified by PCR and the PCR products were subjected to sequencing .
The target nucleotide sequences were compared to sequences that were available in public databases using BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The other hantavirus sequences for HTNV, DOBV, SEOV, SNV, ANDV, BAYV, PUUV, PHV, ISLAV, and LX307 (which should not be identified as TULV ) were set as outgroups. The phylogenetic tree was constructed using MEGA neighbor joining method (MEGA, version 7.0; http://www.megasoftware.net/)  and bootstrap values were 70%, calculated from 1000 replicates.
Vole species identification
A total of 198 voles were captured in a Narati pasture (43°23′ N, 83°96′ E, altitude: above 1800 m) in Xinyuan County, Xinjiang, China. These small mammals came from 20 burrow-groups near the streams and each group contained 15–30 occupied burrows.
All of the voles were M. obscurus, initially classified by morphological identification. Two voles from hantavirus positive-groups that showed slight differences in morphology were further identified by PCR amplification and sequencing of Cytb gene, which confirmed that the voles were M. obscurus with identical nucleotide sequences (GenBank accession No: KX058268-KX058269).
Hantavirus detection and evolution analysis
Hantavirus L gene fragments were successfully amplified and sequenced from 31 out of 198 voles, such that 16% of voles were infected with hantavirus. We further amplified and sequenced the S gene of hantavirus from the two L-gene-positive samples to determine the species of the viruses. The sequence analysis showed the nucleotide sequences of two isolates are identical, called Tula xinjiang4 (GenBank accession No: KX270414). Combining with the sequences published in the GenBank, the phylogenetic analysis showed 5 genetic lineages (I–V) likely based on the geographic distribution in Eurasia (Fig. 1). Tula xinjiang4 is very similar to Karata322 (isolated from eastern Kazakhstan, altitude: 1200 m)  with which it has been found to share 92.1% nucleotide identity and 98.2% amino acid identity, and Omsk23 (isolated from southwest region of Siberia in Russia), with which it has been found to share 87.5% nucleotide identity and 97.5% amino acid identity. The referenced sequences were isolated from European common voles, M. arvalis (obscurus), M. gregalis, and Pitymys subterraneus, respectively. (Additional file 1: Table S1 and Table S2).
Few studies have explored the epidemiology of hantaviruses in Xinjiang. Only Seoul virus (SEOV, a hantavirus), found in Rattus norvegicus in Urumqi, Xinjiang, has been assessed in this way [20, 21]. The present study showed for the first time that M. obscurus [3, 22] carries TULV in China. The distribution of hantavirus is likely related to the distribution of rodents [1, 23]. In Xinjiang, there are 69 species of rodents, which belong to 10 families and 34 genera, accounting for 40% of rodent species in China [23, 24]. Among the small mammals, there are 19 species in subfamily Microtinae and 7 species in the genera Microtus . In this way, it is likely that other voles or small mammal species may also serve as hosts for TULV or hantavirus in Xinjiang. This needs further investigation.
Hantavirus infection is a zoonosis and it remains a severe public health problem. However, very few patients infected with TULV in Europe have shown symptoms [8,9,10] with the most serologically positive cases being asymptomatic . It may be relevant whether the populations in these areas where M. obscurus are likely to be infected with TULV.
The two locations, Yili in China and eastern Kazakhstan, where the TULV-infected voles were found, are geographically close and have similar landscapes; both are in the Tianshan Mountains. Viral sequence analysis showed TULV Xinjiang4 to be very similar to Karata322 and the overall lineages appeared closely related to those isolates from central Russia and Siberia, indicating these viruses share a common ancestor (bootstrap support value 98%, on the NJ-tree).
Hantavirus contains three genes, negative-stranded large (L), middle (M), and small (S) genomic RNA segments encode a polymerase, M proteins composed of two surface glycoproteins (Gc and Gn), and the nucleocapsid protein (N protein), respectively. The complete S segment sequence of TULV strain is 1830 nt in size and the sequence includes the following structures: a 5′-untranslated region (UTR, nt 1–42), an open reading frame for the 430 amino acids (AA-long N protein, nt 43–1335) and a 3′-UTR (nt 1336–1830) . In our study, the predicted N-protein sequence revealed that the N protein sequence region (amino acid, AA 244–273) of TULV contains a highly variable region. There are 3 AA (AA 256–258) deletions in Tula Xinjiang4 and other TULV strains, 4 AAs deleted in hantavirus (HTNV) and SEOV (AA 257–260) compared with Prospect Hill virus (PHV) and Puumala virus (PUUV) (Fig. 2). It may be useful to determine whether the deletion of these sequence impacts the function of N protein.
TULV was identified in Narati pasture, Yili region, Xinjiang, western China by molecular detection. The high prevalence of the virus in M. obscurus indicates Tula virus is naturally transmitted in this area. It may be useful to determine whether this virus infects other hosts, such as humans.
Availability of data and materials
All of the data generated or analysed during this study are included in this published article. The sequences generated have been submitted to GenBank under accession numbers KX058268 to KX058269 for Cytb gene sequences and KX270414 for the partial segment of S genomic exon sequence.
Molecular Evolutionary Genetic Analysis
Open Reading Frame
Reverse Transcription-Polymerase Chain Reaction
Jonsson CB, Figueiredo LT, Vapalahti O. A global perspective on hantavirus ecology, epidemiology, and disease. Clin Microbiol Rev. 2010;23(2):412–41.
Plyusnin A, Vapalahti O, Lankinen H, Lehvaslaiho H, Apekina N, Myasnikov Y, et al. Tula virus: a newly detected hantavirus carried by European common voles. J Virol. 1994;68(12):7833–9.
Mahmoudi A, Darvish J, Aliabadian M, Moghaddam FY, Krystufek B. New insight into the cradle of the grey voles (subgenus Microtus) inferred from mitochondrial cytochrome b sequences. Mammalia. 2017;81(6):583–93.
Schlegel M, Kindler E, Essbauer SS, Wolf R, Thiel J, Groschup MH, et al. Tula virus infections in the Eurasian water vole in Central Europe. Vector Borne Zoonotic Dis. 2012;12(6):503–13.
Reusken C, de Vries A, Adema J, Vos W, van der Giessen J, Bekker D, et al. First genetic detection of Tula hantavirus in wild rodents in the Netherlands. J Inf Secur. 2008;57(6):500–3.
Iakimenko VV, Dekonenko AE, Malkova MG, Kuzmin IV, Tantsev AK, Dzagurova TK, et al. The spread of hantaviruses in western Siberia. Med Parazitol (Mosk). 2000;(3):21–8.
Plyusnina A, Laakkonen J, Niemimaa J, Henttonen H, Plyusnin A. New genetic lineage of Tula hantavirus in Microtus arvalis obscurus in eastern Kazakhstan. Open Virol J. 2008;2:32–6.
Reynes JM, Carli D, Boukezia N, Debruyne M, Herti S. Tula hantavirus infection in a hospitalised patient, France, June 2015. Euro Surveill. 2015. https://doi.org/10.2807/1560-7917.ES.2015.20.50.30095.
Klempa B, Meisel H, Rath S, Bartel J, Ulrich R, Kruger DH. Occurrence of renal and pulmonary syndrome in a region of Northeast Germany where Tula hantavirus circulates. J Clin Microbiol. 2003;41(10):4894–7.
Zelena H, Mrazek J, Kuhn T. Tula hantavirus infection in immunocompromised host, Czech Republic. Emerg Infect Dis. 2013;19(11):1873–5.
Mertens M, Hofmann J, Petraityte-Burneikiene R, Ziller M, Sasnauskas K, Friedrich R, et al. Seroprevalence study in forestry workers of a non-endemic region in eastern Germany reveals infections by Tula and Dobrava-Belgrade hantaviruses. Med Microbiol Immunol. 2011;200(4):263–8.
Schmaljohn CS, Jennings GB, Hay J, Dalrymple JM. Coding strategy of the s genome segment of hantaan virus. Virology. 1986;155(2):633–43.
Klempa B, Fichet-Calvet E, Lecompte E, Auste B, Aniskin V, Meisel H, et al. Hantavirus in African wood mouse, Guinea. Emerg Infect Dis. 2006;12(5):838–40.
Tkachenko EA, Witkowski PT, Radosa L, Dzagurova TK, Okulova NM, Yunicheva YV, et al. Adler hantavirus, a new genetic variant of Tula virus identified in Major's pine voles (Microtus majori) sampled in southern European Russia. Infect Genet Evol. 2015;29:156–63.
Schlegel M, Ali HS, Stieger N, Groschup MH, Wolf R, Ulrich RG. Molecular identification of small mammal species using novel cytochrome b gene-derived degenerated primers. Biochem Genet. 2012;50(5–6):440–7.
Wang H, Yoshimatsu K, Ebihara H, Ogino M, Araki K, Kariwa H, et al. Genetic diversity of hantaviruses isolated in China and characterization of novel hantaviruses isolated from Niviventer confucianus and Rattus rattus. Virology. 2000;278(2):332–45.
Cao ZW, Zuo SQ, Ma CF. Sequence analysis of partail L and S segments of hantavirus in Xi’an, China. Acta Parasitol Et Med Entomol. 2010;17(2):98–102 (Chinese).
Wang JL, Zhang HL, Zhou JH, Yang WH, Zhang YZ, Mi ZQ, et al. First detection of Tula-like hantavirus from Eothenomys milelus in China. CHN J Zoonoses. 2010;26(5):408–12 (Chinese).
Sohpal VK, Dey A, Singh A. MEGA biocentric software for sequence and phylogenetic analysis: a review. Int J Bioinforma Res Appl. 2010;6(3):230–40.
Guo G, Sheng J, Wu X, Wang Y, Guo L, Zhang X, et al. Seoul virus in the Brown rat ( Rattus norvegicus ) from Urumqi, Xinjiang, northwest of China. J Wildl Dis. 2016;52(3):705–8.
Zhang D, Mahmut JT. Spreading again along with the railway of Rattus norvegicus and its regulatiag methods in Xinjiang. Endemic Dis Bull. 1997;12(2):106–8 (Chinese).
Tougard C, Montuire S, Volobouev V, Markova E, Contet J, Aniskin V, et al. Exploring phylogeography and species limits in the Altai vole (Rodentia: Cricetidae). Biol J Linn Soc. 2013;108(2):434–52.
Zhang YZ, Zou Y, Fu ZF, Plyusnin A. Hantavirus infections in humans and animals, China. Emerg Infect Dis. 2010;16(8):1195–203.
Wang S, Yang G. Xinjiang Rodents. 1st ed. Urumqi: Xinjiang People's Publishing House; 1983 (Chinese).
This work was supported by the Natural Science Foundation of China (31860703) and co-funded by the Project of State Key Laboratory of Pathogenesis, Prevention, and Treatment of High Incidence Diseases in Central Asia (SKL-HIDCA-2017-5). The grantors were not involved in the study design, analysis or interpretation, or decision to submit this manuscript.
Ethics approval and consent to participate
The entire rodent trapping campaign and all of the animal experiments were reviewed and approved by Ethics Committee of the First Affiliated Hospital of Xinjiang Medical University. In addition, permission to capture the voles and small mammals was obtained from each farmer before entering their farm. All of the experimental procedures were carried out in accordance with laws for the prevention of cruelty to animals and with approved guidelines. None of the rodent species investigated in the present study are protected in China and none of the species captured are included in the China Species Red List.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Identities and differences among nucleotide sequences of TULV virus strains from Eurasia. Table S2. Identities and difference among amino acid sequences deduced for TULV virus strains from Eurasia. (DOC 1366 kb)
About this article
Cite this article
Guo, G., Guo, B., Wu, X. et al. Molecular evidence of Tula virus in Microtus obscurus in the region of Yili, Xinjiang, China. BMC Infect Dis 19, 527 (2019). https://doi.org/10.1186/s12879-019-4133-x