Skip to main content
  • Research article
  • Open access
  • Published:

Molecular detection and genetic characterization of Wenzhou virus in rodents in Guangzhou, China



Wenzhou virus (WENV), a newly discovered mammarenavirus in rodents, is associated with fever and respiratory symptoms in humans. This study was aimed to detect and characterize the emerging virus in rodents in Guangzhou, China.


A total of 100 small mammals, including 70 Rattus norvegicus, 22 Suncus murinus, 4 Bandicota indica, 3 Rattus flavipectus, and 1 Rattus losea, were captured in Guangzhou, and their brain tissues were collected and pooled for metagenomic analysis, which generated several contigs targeting the genome of WENV. Two R. norvegicus (2.9%) were further confirmed to be infected with WENV by RT-PCR. The complete genome (RnGZ37-2018 and RnGZ40-2018) shared 85.1–88.9% nt and 83.2–96.3% aa sequence identities to the Cambodian strains that have been shown to be associated with human disease. Phylogenetic analysis showed that all identified WENV could be grouped into four different lineages, and the two Guangzhou strains formed an independent clade. We also analyzed the potential recombinant events occurring in WENV strains.


Our study showed a high genetic diversity of WENV strains in China, emphasizing the relevance of surveillance of this emerging mammarenavirus in both natural reservoirs and humans.


Mammarenaviruses, belonging to the Mammarenavirus genus in the Arenaviridae family, include several members responsible for severe hemorrhagic fever diseases. Of them, Lassa virus (LASV) may cause an approximate 100,000-300,000 case and 5000 deaths annually in Western Africa (Lassa Fever CDC, 2019,, while Junín virus (JUNV), Guanarito virus (GTOV), Machupo virus (MACV), and Sabía virus (SABV) can cause fatal hemorrhagic fever diseases in Latin America [1]. In Asia, the known mammarenaviruses are lymphacytic choriomeningitis virus (LCMV), Wenzhou virus (WENV), Loei River virus (LORV), Ryukyo virus (RVKV), Xingyi and Lijiang viruses. LCMV infection can lead to various clinical manifestations, ranging from asymptomatic to influenza-like illness, aseptic meningitis, or meningoencephalitis. LCMV is also associated with congenital microencephaly in pregnant women [2, 3]. WENV, first discovered in 2014 in Wenzhou, Zhejiang Province, China [4], has been shown to be related to fever and respiratory symptoms in humans [5]. LORV was reported in 2016 from Thailand, and RVKV, Xingyi and Lijiang viruses were all found in China during 2018–2019, whose relevance to public health remains unclear [5,6,7].

Mammarenaviruses are negative-sense single-stranded RNA viruses that include the small (S) and large (L) genomic segments. The S segment encodes an envelope glycoprotein precursor (GPC) and a nucleoprotein (NP), and the L segment encodes a RNA-dependent RNA polymerase (RdRp) and a zinc binding matrix protein (ZP) [8]. Rodents are considered as the primary reservoir for mammarenaviruses, with several exceptions. For example, Tacaribe virus (TCRV) has been isolated from fruit-eating bats (Artibeus jamaicensis) and ticks, LCMV has been isolated from ticks, and WENV has also been detected in shrews [4, 9,10,11]. Based on genetic and geographic relationships, mammarenaviruses are divided into two monophyletic groups, including the New World (NW) and the Old World (OW). The NW arenaviruses, including JUNV, GTOV, MACV, SABV, WWAV, and TCRV, are geographically confined to Americas. While the OW arenaviruses are generally distributed in Africa (LASV) and Asia (WENV, LORV, RVKV, Xingyi and Lijiang virus) [8]. The high genetic diversity has been shown in the NW and OW arenaviruses [12, 13], and WENVs from different geographic area have formed independent clades [7, 14].

Diverse zoonotic pathogens (i.e., hantaviruses, hepatitis viruses, leptospira spp., bartonella spp.) have been reported in rodents in Guangdong Province [15, 16]. However, little is known about WENV circulation status in Guangzhou. Herein, we attempt to detect and characterize of WENV, which may contribute to understanding the evolution of WENV.


Detection of WENV in rodents in Guangzhou

In 2018, a total of 100 rodent samples were collected from Conghua and Haizhu Districts in Guangzhou, Guangdong Province, China (Fig. 1). The rodents included 3 genera and 5 species, including 70 Rattus norvegicus, 22 Suncus murinus, 4 Bandicota indica, 3 Rattus flavipectus and 1 Rattus losea (Table 1). Brain tissues of each rodent were collected for viral metagenomic analysis, which produced 2,722 contigs, with an average length of 211 nt. There were 12 (0.44%) contigs annotated to the L and S segment of WENV, with 88.5–92.9% and 90.8–95.8% nt identities, respectively (Table S1). Semi-nested RT-PCR targeting the 290 nt fragment of S segment was performed to screen WENV infection in all collected samples, showing that two R. norvegicus (2/70, 2.9%) in Conghua District were positive for WENV. One of the two R. norvegicus was male, adult; the other was female, juvenile. No WENV infection was detected in B. indica, R. flavipectus, R. losea, or S. murinus.

Fig. 1
figure 1

The sampling sites in Guangzhou, Guangdong Province, China. This map was drawn by us for this study, and plotted by combination of Surfer software version 4 (Golden software, USA) and Adobe illustrator CS6 (Adobe Systems, USA)

Table 1 Detection of WENV in rodents captured in Guangzhou, China

Genetic analysis of WENV

The complete S and L segments were recovered from the two R. norvegicus samples (RnGZ37-2018 and RnGZ40-2018), which comprised the S (3,350 nt) and L (7,118 nt) segments, each containing two open reading frames (ORFs) in an ambisense organization with an intergenic region (IGR). L segment encoded a ZP of 276 nt (91 aa) and a RdRp of 6672 nt (2,223 aa), connected by an IGR of 122 nt. S segment encoded a NP of 1704 nt (567 aa) and a GPC of 1479 nt (492 aa), connected by an IGR of 62 nt (Table S4). Interestingly, the 122 nt IGR of the L segment from WENV in this study was longer than other WENV isolates (~ 119 nt) [7].

Sequence comparison revealed that RnGZ37-2018 and RnGZ40-2018 were closely related with each other, with 93.5–98.6% nt and 95.1–97.6% aa identities, and had 83.5–93.3% nt and 84.2–91.5% identity for the S and L segments in comparison with other WENV strains, respectively (Table 2). The nt and aa identity for RdRp, ZP, NP and GPC were also shown in Table 2. RnGZ37-2018 and RnGZ40-2018 were most closely related to the 9–24 strain, with 90.9–93.9% nt and 93.3–98.0% aa identity. Notably, both RnGZ37-2018 and RnGZ40-2018 were also closely related to the Cambodian strains, which were identified in flu-like patients, with 85.1–88.9% nt and 83.2–96.3% aa identity [5].

Table 2 Genome and ORFs comparisons between RnGZ37-2018, RnGZ40-2018 and other representative WENV strains

Phylogenetic analysis of WENV

Phylogenetic analysis based on the full-length of the S and L segments showed that WENVs formed a monophyletic clade distinct from other OW arenaviruses [4, 14]. RnGZ37-2018 and RnGZ40-2018 were clustered together and formed an independent clade within WENV (Fig. 2). According to the evolutionary relationships, both S and L segment of all WENV strains (except WENV G107) could be grouped into 4 different lineages. Based on the S segment tree, lineage I included WENV strains from Zhejiang and Hainan; the novel WENVs identified in Guangzhou were clustered in lineage II; lineage III contained the strains identified in patients in Cambodia; and lineage IV was comprised the virus isolated in Yunnan and Xinjiang (Fig. 2A). Phylogenetic topology of L segment was the same as the S segment (Fig. 2B). Genetic diversity and protein differences within and between the viral lineages were shown in Table 3. Genetic distance of within or between lineages on the L segment was like the S segment, with the mean p-distance ranging from 0.07 to 0.75. The aa divergence of RdRp, ZP, GPC, and NP were 0.33, 0.27, 0.25, and 0.09, respectively, suggesting that NP is highly conserved in comparison with other proteins.

Fig. 2
figure 2

The maximum likelihood trees based on the full-length genomic sequences of WENV. The phylogenetic trees of the L (A) and S (B) segments of WENV were constructed using ML methods based on GTR model with gamma distribution and invariant sites, and tested with 1000 bootstrap replicates. Sequence were identified by the GenBank accession number and the strain name, followed by their origin. The newly identified WENV variants in Guangzhou were indicated with black circles

Table 3 Genetic and protein diversity of WENV strains

Recombination analysis of WENV

Both RDP4 and Simplot methods were used to analyze potential recombination events of WENVs. No obvious recombination events were found in RnGZ37-2018 and RnGZ40-2018 (Figure S1 and Table S5). Interestingly, the potential recombinant events in WENV isolates were found using RDP4 (p < 0.05 and the RDP recombination consensus score > 0.6), showing that the potential recombinant events occurred in the S segment of RnYM3-2016 between RtYM51-2015 and RnYL4-2016, and the Rn242 between DK and WZ140510, also as in the L segment of WZ140510 between WZ140512 and an unknown isolate (Table S5).


WENV was first identified in rodents and shrews by nested RT-PCR using the universal primers of arenaviruses in Wenzhou, Zhejiang Province in China in 2014, where R. flavipectus (15.4%), R. losea (11.8%), R. norvegicus (17.1%), R. rattus (75.0%), Niviventer niviventer (1.9%), and S. murinus (4.4%) were tested positive [4]. Since then, the virus has been found in R. norvegicus in Shandong, Hainan and Xinjiang Provinces [7, 17, 18], also found in R. losea in Hunan Province [6]. In Shenzhen, Guangdong Province, WENV was found in R. norvegicus and S. murinus with an infection rate of 6.7% and 0.5%, respectively [15]. In Yunnan Province, WENV has been found in R. nitidus, R. norvegicus, R. tanezumi, and Tupaia belangeri, with an infection rate of 8.3%, 17.4%, 1.5%, and 10%, respectively [14]. WENV has also been detected in Southeast Asia, including in 14.6% R. exulans and 21.1% R. norvegicus from one Cambodian province, whilst the closely related LORV was found in 25% B. indica, 28.6% Bandicota savilei and 1.9% Niviventer fulvescens in the Thai province of Loei [5]. In this study, two WENV isolates RnGZ37-2018 and RnGZ40-2018 were identified in the brain tissues of R. norvegicus (2.9%) but none was found in other rodent species in Guangzhou (Table 1), probably owing to the limited numbers and the brain tissue used. S. murinus was the carrier of WENV with a rate varying from 0.5–4.4%, and our results were consistent with previous studies showing a rate of 0.5% in Shenzhen [4, 15]. WENV has been shown to replicate in brain tissue in an experimental rat model, with higher viral loads in liver, lung and spleen [5, 7], which may explain the relatively lower rate in the brain tissues [4, 14, 15].

R. norvegicus has been considered as the host of WENV, which is also infected with LCMV [19, 20]. LCMV has been divided into 3 and 4 lineages based on L and S segment, respectively, and shows no clear correlation between the genetic diversity and geographic locations [10, 12]. To date, WENV has been reported only in Asia, and clusters together according to their geographic distribution [7, 14]. Remarkably, though Shenzhen and Guangzhou are over 100 kilometers away, the WENV forms two different sister clades (Figure S2), showing a high genetic diversity of WENV, which is consistent with previous study [21]. The 9–24 strain was most closely related to the Guangzhou strains but the geographic origin of this strain is unknown. Thus, the geographical relationship between the 9–24 strain and the Guangzhou strains was not discussed in this study. Interestingly, the IGR of L segment of WENV identified in Guangzhou was longer than most known WENV (Table S4), and substitution or deletion of IGR in L segment in LASV and MACV can reduce viral titer compared to wide type [22, 23]. However, the effect of the longer IGR on viral replication needs further study. Recombination and reassortment may contribute to forming high genetic diversity in mammarenaviruses [24, 25]. The first evidence of recombination in S segments within WENV was sampled in Yunnan Province with RnYM3-2016 isolates [14]. The potential recombinant event in S segment with Rn242 (Zhejiang Province) and L segment with WZ140510 (Zhejiang Province) were also observed (Table S5). These data have shown a complex evolutionary history of WENVs.

WENV infection may cause diffuse pneumonia in the experimental rat models [5, 7], and also cause a small cranial cavity in R. exulans, suggesting microcephaly as observed in LCMV [26, 27]. Moreover, WENV has been found seropositive (17.4%) in dengue-like/influenza-like patients and healthy individuals in Southeastern Asia, and 4.6% in healthy individuals in China [5, 28]. It is likely to show that the impact of WENV infections on public health may be underestimated, as the symptoms are mild and similar to those of respiratory diseases.


New genotypes of WENV in Rattus norvegicus from Guangzhou of China were identified. The results showed that WENVs from Guangzhou clustered together, and they were closely related to the Cambodian strains associated with human disease, and WENV from Guangzhou also showed highly genetic diversity in China. These findings highlight the importance of surveillance of WENV in both natural reservoirs and humans.


Specimen collection

In 2018, wild rodents were trapped using cages in residential areas in two districts (Conghua and Haizhu) of Guangzhou (Fig. 1). Trap cages were placed next to garbage cans and human settlements. The rodents were euthanized by the inhalation of carbon dioxide by trained personnel with appropriate technique, equipment, and agents [29]. The trapped rodent transferred to euthanasia chamber connected with a CO2 tank which was equipped with a pressure regulator and flowmeter. Exposed to 30% CO2 for at least 3 minutes or longer at a flow rate of 3.4 L/min until rodents became unconsciousness, even stopping breathing. Leaving rodents exposed to 70% CO2 for an additional 2 minutes at a flow rate of 7.9 L/min and turning off the CO2. Death was ascertained by absence of movement, cardiac and respiratory arrest, fixed and dilated pupils. Death was further confirmed by decapitation (> 200 g) or cervical dislocation (≤ 200 g). The blood and organs were collected by necropsy, then stored at -80℃. Rodents species were classified by morphology and subsequently confirmed by sequencing the mitochondria cytochrome b gene [30]. The species, sex, weight, and age category (juvenile or adult) were recorded in Table 1.

Viral metagenomic analysis

The rodent brain tissues (~ 50 mg) were homogenized with Qiagen tissue lyser II (Qiagen), and used for viral metagenomic analysis as previously described [31, 32]. Briefly, the viral nucleic acids were extracted and reverse transcribed with anchored random primers. The double-strand cDNA was amplified using sequence-independent, single-primer amplification, then PCR products were purified. The equal amount of the purified PCR products from 100 samples were pooled together and sent to the Beijing Genome Institute (BGI, Shenzhen, China) for high-throughput sequencing. All generated sequences were subjected to BLASTx and BLASTn search against the viral reference database for further analysis.

WENV screening and complete genome recovery

Based on the WENV-like contigs from the metagenomic analysis, we designed the semi-nested RT-PCR primers that targeted a 290 nt fragment of the S segment (F1748: 5'-TTCTTCTTTTCAACAACCAC-3'; R1-2055: 5'-GAGCCAACAGACGCCAAG-3'; R2-2037: 5'-CAGTTCAAGCAAGACTCT-3'). Total RNA was extracted with MiniBest viral RNA/DNA extraction kit (TaKaRa), and reverse transcription was conducted with the 1st cDNA synthesis kit (TaKaRa). cDNA was amplified by using the ExTaq (TaKaRa) with the following PCR parameters: 94℃ for 5 min, followed by 35 cycles of 94℃ for 35 s, 51℃ for 35 s, 72℃ for 25 s, and a final elongation step at 72℃ for 10 min. The PCR products were analyzed on gel electrophoresis, and positive samples were further sequenced.

To get complete genomic sequences, overlapping primers (Table S2) that covered the terminal ends were designed according to the WENV-like contigs and reference sequences (GenBank accession numbers MF414208 and MF414207). The amplicons were cloned to the pMD18-T vector (TaKaRa), and at least two clones for each PCR products were sequenced to obtain a consensus sequence. The resulting sequences were assembled into the complete genome using SeqMan in DNASTAR.

Genomic and phylogenetic analysis

The complete genome sequences of WENV representative strains were retrieved from GenBank and listed in Table S3. Pairwise sequence alignment was performed using DNASTAR [33]. Phylogenetic analysis was conducted with the maximum likelihood (ML) method based on the general time reversible model with gamma distribution and invariant sites, and tested with 1000 bootstrap replicates in MEGA X program, also the genetic distance was tested by p-distance with transitions and transversions method [34]. The recombination of S and L segments was analyzed by Recombination Detection Program version 4 (RDP4) and Simplot [35, 36].

Availability of data and materials

The data generated or analyzed in this study are included in the manuscript and the supplementary information files. The complete genome sequences obtained in this study have been submitted to GenBank with accession numbers: MW174777, MW174778, MW174779, and MW174780.



Wenzhou virus


New World arenaviruses


Old World arenaviruses


Glycoprotein precursor




RNA-dependent RNA polymerase


Zinc binding matrix protein


Open reading frames


Amino acid




Intergenic region


Lassa virus


Lymphocytic choriomeningitis virus


Junín virus


Guanarito virus


Machupo virus


Sabía virus


Loei River virus


Ryukyo virus


Maximum likelihood method


Recombination Detection Program version 4

B. indica :

Bandicota indica; R. norvegicus

Rattus norvegicus :

R. rattus:Rattus rattus

R. flavipectus :

Rattus flavipectus

R.losea :

Rattus losea

R. exulans :

Rattus exulans

R.nitidus :

Rattus nitidus

R. tanezumi :

Rattus tanezumi

S. murinus :

Suncus murinus


  1. Hallam SJ, Koma T, Maruyama J, Paessler S. Review of Mammarenavirus Biology and Replication. Front Microbiol. 2018;9:1751.

    Article  Google Scholar 

  2. Jamieson DJ, Kourtis AP, Bell M, Rasmussen SA. Lymphocytic choriomeningitis virus: an emerging obstetric pathogen? Am J Obstet Gynecol. 2006;194(6):1532–6.

    Article  Google Scholar 

  3. Bonthius DJ. Lymphocytic choriomeningitis virus: an underrecognized cause of neurologic disease in the fetus, child, and adult. Semin Pediatr Neurol. 2012;19(3):89–95.

    Article  Google Scholar 

  4. Li K, Lin XD, Wang W, Shi M, Guo WP, Zhang XH, et al. Isolation and characterization of a novel arenavirus harbored by Rodents and Shrews in Zhejiang province, China. Virology. 2015;476:37–42.

    Article  CAS  Google Scholar 

  5. Blasdell KR, Duong V, Eloit M, Chretien F, Ly S, Hul V, et al. Evidence of human infection by a new mammarenavirus endemic to Southeastern Asia. Elife. 2016;5:e13135.

    Article  Google Scholar 

  6. Wu Z, Lu L, Du J, Yang L, Ren X, Liu B, et al. Comparative analysis of rodent and small mammal viromes to better understand the wildlife origin of emerging infectious diseases. Microbiome. 2018;6(1):178.

    Article  Google Scholar 

  7. Tan Z, Yu H, Xu L, Zhao Z, Zhang P, Qu Y, et al. Virome profiling of rodents in Xinjiang Uygur Autonomous Region, China: Isolation and characterization of a new strain of Wenzhou virus. Virology. 2019;529:122–34.

    Article  CAS  Google Scholar 

  8. Radoshitzky SR, Bao Y, Buchmeier MJ, Charrel RN, Clawson AN, Clegg CS, et al. Past, present, and future of arenavirus taxonomy. Arch Virol. 2015;160(7):1851–74.

    Article  CAS  Google Scholar 

  9. Downs WG, Anderson CR, Spence L, Aitken TH, Greenhall AH. Tacaribe virus, a new agent isolated from Artibeus bats and mosquitoes in Trinidad, West Indies. Am J Trop Med Hyg. 1963;12:640–6.

    Article  CAS  Google Scholar 

  10. Zhang L, Li S, Huang SJ, Wang ZD, Wei F, Feng XM, et al. Isolation and genomic characterization of lymphocytic choriomeningitis virus in ticks from northeastern China. Transbound Emerg Dis. 2018;65(6):1733–9.

    Article  CAS  Google Scholar 

  11. Sayler KA, Barbet AF, Chamberlain C, Clapp WL, Alleman R, Loeb JC, et al. Isolation of Tacaribe virus, a Caribbean arenavirus, from host-seeking Amblyomma americanum ticks in Florida. PLoS One. 2014;9(12):e115769.

    Article  Google Scholar 

  12. Albarino CG, Palacios G, Khristova ML, Erickson BR, Carroll SA, Comer JA, et al. High diversity and ancient common ancestry of lymphocytic choriomeningitis virus. Emerg Infect Dis. 2010;16(7):1093–100.

    Article  Google Scholar 

  13. Archer AM, Rico-Hesse R. High genetic divergence and recombination in Arenaviruses from the Americas. Virology. 2002;304(2):274–81.

    Article  CAS  Google Scholar 

  14. Wang J, Yang X, Liu H, Wang L, Zhou J, Han X, et al. Prevalence of Wenzhou virus in small mammals in Yunnan Province, China. PLoS Negl Trop Dis. 2019;13(2):e0007049.

    Article  CAS  Google Scholar 

  15. Wang B, Cai CL, Li B, Zhang W, Zhu Y, Chen WH, et al. Detection and characterization of three zoonotic viruses in wild rodents and shrews from Shenzhen city, China. Virol Sin. 2017;32(4):290–7.

    Article  CAS  Google Scholar 

  16. Su Q, Chen Y, Wang B, Huang C, Han S, Yuan G, et al. Epidemiology and genetic diversity of zoonotic pathogens in urban rats (Rattus spp.) from a subtropical city, Guangzhou, southern China. Zoonoses Public Health. 2020;67(5):534–45.

  17. Zhang Y, Rao L, Lv G, Zhang L, Fu R, Wang S, et al. Complete Genome Sequence of a Mammarenavirus Harbored by Rodents on Hainan Island, China. Genome Announc. 2018;6(10):e00129–18.

  18. Liu MM, Li LL, Wang XF, Duan ZJ. Complete Genome Sequence of a Novel Variant of Wenzhou Mammarenavirus. Genome Announc. 2017;5(47):e01303–17.

    Article  Google Scholar 

  19. Moolla N, Weyer J. The Arenaviridae. In: Ennaji MM, editor. Emerging and Reemerging Viral Pathogens. London: Academic Press; 2020. p. 69–100.

    Chapter  Google Scholar 

  20. Patterson BD. A New 'Age of Discovery' for Mammals Mammal Species of the World: A Taxonomic and Geographic Reference, 3rd Edition. Don E. Wilson and DeeAnn M. Reeder, editors. J Mamm Evol. 2007;14(1):67–9.

  21. Li K, Lin XD, Li MH, Wang MR, Sun XY, Zhang YZ. Genomic analysis of Wenzhou virus in rodents from Zhejiang province. Zhonghua Liu Xing Bing Xue Za Zhi. 2017;38(3):384–7.

    PubMed  Google Scholar 

  22. Golden JW, Beitzel B, Ladner JT, Mucker EM, Kwilas SA, Palacios G, et al. An attenuated Machupo virus with a disrupted L-segment intergenic region protects guinea pigs against lethal Guanarito virus infection. Sci Rep. 2017;7(1):4679.

    Article  Google Scholar 

  23. Cai Y, Iwasaki M, Motooka D, Liu DX, Yu S, Cooper K, et al. A Lassa Virus Live-Attenuated Vaccine Candidate Based on Rearrangement of the Intergenic Region. mBio. 2020;11(2):e00186–20.

    Article  Google Scholar 

  24. Pontremoli C, Forni D, Sironi M. Arenavirus genomics: novel insights into viral diversity, origin, and evolution. Curr Opin Virol. 2019;34:18–28.

    Article  CAS  Google Scholar 

  25. Emonet SF, de la Torre JC, Domingo E, Sevilla N. Arenavirus genetic diversity and its biological implications. Infect Genet Evol. 2009;9(4):417–29.

    Article  CAS  Google Scholar 

  26. Bonthius DJ, Perlman S. Congenital viral infections of the brain: lessons learned from lymphocytic choriomeningitis virus in the neonatal rat. PLoS Pathog. 2007;3(11):e149.

    Article  Google Scholar 

  27. Artois J, Blasdell K, Duong V, Buchy P, Hul V, Morand S, et al. Effects of mammarenavirus infection (Wenzhou virus) on the morphology of Rattus exulans. Infect Genet Evol. 2018;63:404–9.

    Article  CAS  Google Scholar 

  28. Guo L, Liu S, Song J, Han L, Zhang H, Wu C, et al. Seroprevalence of Wenzhou virus in China. Biosaf Health. 2020;2(3):152–6.

    Article  Google Scholar 

  29. AVMA. AVMA Guidelines for the Euthanasia of Animals: 2020 Edition. 2020.

  30. 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.

    Article  CAS  Google Scholar 

  31. He B, Li Z, Yang F, Zheng J, Feng Y, Guo H, et al. Virome profiling of bats from Myanmar by metagenomic analysis of tissue samples reveals more novel Mammalian viruses. PLoS One. 2013;8(4):e61950.

    Article  CAS  Google Scholar 

  32. Wang ZD, Wang B, Wei F, Han SZ, Zhang L, Yang ZT, et al. A New Segmented Virus Associated with Human Febrile Illness in China. N Engl J Med. 2019;380(22):2116–25.

    Article  CAS  Google Scholar 

  33. Higgins DG, Thompson JD, Gibson TJ. Using CLUSTAL for multiple sequence alignments. Methods Enzymol. 1996;266:383–402.

    Article  CAS  Google Scholar 

  34. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  Google Scholar 

  35. Norberg P, Tyler S, Severini A, Whitley R, Liljeqvist JA, Bergstrom T. A genome-wide comparative evolutionary analysis of herpes simplex virus type 1 and varicella zoster virus. PLoS One. 2011;6(7):e22527.

    Article  CAS  Google Scholar 

  36. Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1(1):vev003.

    Article  Google Scholar 

Download references


We are grateful to Guangdong Medical Waste Treatment Center Co., Ltd. for treatment with the autoclaved rodent carcass.


This work was supported by the Science and technology Project of Guangzhou, Guangdong Province of China (202103000008), the Pearl River Talent Plan in Guangdong Province of China (2019CX01N111) and the Key Laboratory for Preventive Research of Emerging Animal Diseases in Foshan University (KLPREAD201801-13). The funding bodies had no role in the design, analysis and writing the manuscript.

Author information

Authors and Affiliations



QL and NW conceived and designed the experiments. SC collected the samples. NW, LY, GL performed the experiment. NW, XZ, JS, and JM analyzed the data and wrote the manuscript. QL revised the manuscript. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Shouyi Chen or Quan Liu.

Ethics declarations

Ethics approval and consent to participate

The sampling sites in Guangzhou were approved by local community, and during specimen collection, all animal work was complied with the guidelines relating to the recommendations from the Guide for the Care and Use of Laboratory Animals of the Ministry of Health, China. The study protocol was approved by the Research Ethics Committee of the College of Life Science and Engineering, Foshan University, also approved by the Animal Ethics and Welfare Committee of Guangzhou Center for Disease Control and Prevention, and adhered to the guidelines for the Rules for the Implementation of Laboratory Animal Medicine (1998) from the Ministry of Health, and Law on Animal Epidemic Prevention (2015) from the Ministry of Agriculture, China. All surgical procedures were performed in efforts to minimize suffering. Endangered or protected species were not involved in this study.

Consent for publication

Not applicable.

Competing interests

The authors have declared no competing financial or non-financial interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Supplementary Table 1. The contigs annotated to WENV. Supplementary Table 2. The primers used for amplifying the complete genome of WENV. Supplementary Table 3. The reference sequences used in this study. Supplementary Table 4. The genome organization of Wenzhou virus strains identified in Guangzhou. Supplementary Table 5. WENV recombination events detected using the RDP package. Supplementary Fig. 1. The recombination analysis of RnGZ37-2018 and RnGZ40-2018 by similarity plots. The recombination events were analyzed based on the full-length of L/S from RnGZ37-2018, RnGZ40-2018 and other representative WENV strains. RnGZ37-2018 and RnGZ40-2018 were as query sequence. The recombination analysis: (A) S segment of RnGZ37-2018; (B) L segment of RnGZ37-2018; (C) S segment of RnGZ40-2018; (D) L segment of RnGZ40-2018. Supplementary Fig. 2. Phylogenetic analysis of the partial RdRp gene (589 bp) from WENV in Guangzhou and Shenzhen, China. The phylogenetic tree was performed using ML methods based on GTR model with gamma distribution and invariant sites, and evaluated with 1000 bootstrap replicates. Sequences were identified by the GenBank accession number and the strain name, followed by their origin. The partial sequences of WENV isolated from Shenzhen and Guangzhou were indicated.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wang, N., Yang, L., Li, G. et al. Molecular detection and genetic characterization of Wenzhou virus in rodents in Guangzhou, China. BMC Vet Res 17, 301 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: