Metagenomic Analysis of Tick-Borne Viruses in Hulunbuir, Inner Mongolia, China: Epidemiological Risk of Potential Novel Pathogenic Viruses Relevant to Public Health

Abstract

Objective: Inner Mongolia, with its vast grasslands providing ideal habitats for ticks, has recently emerged as a hotspot of novel tick-borne virus (TBV) prevalence. However, the diversity of TBVs and their potential spillover risks remain to be elucidated.

Methods: In 2021, we collected 500 Haemaphysalis longicornis, Ixodes persulcatus, and Dermacentor silvarum ticks from Hulunbuir, Inner Mongolia, and conducted virome analysis and viral surveillance to evaluate TBV exposure risk in Hulunbuir.

Results: A total of 26 viral species belonging to nine families were identified, including four novel viruses (Yiliekede tick virus 1 [YLTV1, Rhabdoviridae], Meitian tick virus [MtTV,Phenuiviridae], Wunuer pestivirus, and Wunuer tick phlebovirus) and 11 new strains of human-pathogenic viruses (including Songling virus [SGLV], Beiji nairovirus [BJNV], and Nuomin virus [NOMV]), from three dominant tick species. Virome profiling revealed striking species-specific patterns:I. persulcatus exhibited the highest viral richness and abundance, dominated by Chuviridae (NOMV) and Phenuiviridae; D. silvarum uniquely carried Phenuiviridae (MtTV); and H. longicornis showed Nairoviridae predominance (SGLV, BJNV). Phylogenetic analyses demonstrated that MtTV clustered with human-pathogenic Tacheng tick virus 2, and NOMV-YL strains showed high similarity to clinical NOMV isolates from febrile patients. Regional surveillance demonstrated MtTV’s pan-regional circulation inD. silvarum across nine locations (MIR=4–8%), contrasting with localized Yiliekede foci for SGLV (MIR=0.6%), BJNV (0.4%), and NOMV (1.0%).

Conclusions: This study not only mapped the landscape of TBVs in a key pastoral region, but also highlighted the potential spillover risk of TBVs, including MtTV, YLTV1, SGLV, BJNV, and NOMV. Our findings underscore the urgent need for One Health strategies for monitoring emerging TBVs at the tick-human-livestock interface in northern China.

Keywords

Tick-borne virus / metagenomics / virome / epidemical investigation / ecological distribution / Hulunbuir

References

[1] Pielnaa P, Al-Saadawe M, Saro A, Dama MF, Zhou M, Huang Y, et al. Zika virus-spread, epidemiology, genome, transmission cycle, clinical manifestation, associated challenges, vaccine and antiviral drug development. Virology. 2020; 543: 34-42.

[2] Kularatne SA, Dalugama C. Dengue infection: global importance, immunopathology and management. Clin Med (Lond). 2022; 22(1): 9-13.

[3] Tuells J, Henao-Martínez AF, Franco-Paredes C. Yellow fever: a perennial threat. Arch Med Res. 2022; 53(7): 649-657.

[4] Kuehnert PA, Stefan CP, Badger CV, Ricks KM. Crimean-Congo hemorrhagic fever virus (CCHFV): a silent but widespread threat. Curr Trop Med Rep. 2021; 8(2): 141-147.

[5] Pulkkinen LIA, Butcher SJ, Anastasina M. Tick-borne encephalitis virus: a structural view. Viruses. 2018; 10(7): 350.

[6] Li J, Li S, Yang L, Cao P, Lu J. Severe fever with thrombocytopenia syndrome virus: a highly lethal bunyavirus. Crit Rev Microbiol. 2021; 47(1): 112-125.

[7] Krasteva S, Jara M, Frias-De-Diego A, Machado G. Nairobi Sheep Disease virus: a historical and epidemiological perspective. Front Vet Sci. 2020; 7: 419.

[8] Gaudreault NN, Madden DW, Wilson WC, Trujillo JD, Richt JA. African swine fever virus: an emerging DNA arbovirus. Front Vet Sci. 2020; 7: 215.

[9] Moming A, Bai Y, Wang J, Zhang Y, Tang S, Fan Z, et al. The known and unknown of global tick-borne viruses. Viruses. 2024; 16(12): 1807.

[10] 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-2125.

[11] Ma J, Lv XL, Zhang X, Han SZ, Wang ZD, Li L, et al. Identification of a new orthonairovirus associated with human febrile illness in China. Nat Med. 2021; 27(3): 434-439.

[12] Zhang XA, Ma YD, Zhang YF, Hu ZY, Zhang JT, Han S, et al. A new orthonairovirus associated with human febrile illness. N Engl J Med. 2024; 391(9): 821-831.

[13] Dong Z, Yang M, Wang Z, Zhao S, Xie S, Yang Y, et al. Human Tacheng tick virus 2 infection, China, 2019. Emerg Infect Dis. 2021; 27(2): 594-598.

[14] Wang YC, Wei Z, Lv X, Han S, Wang Z, Fan C, et al. A new nairo-like virus associated with human febrile illness in China. Emerg Microbes Infect. 2021; 10(1): 1200-1208.

[15] Wang ZD, Wang W, Wang NN, Qiu K, Zhang X, Tana G, et al. Prevalence of the emerging novel Alongshan virus infection in sheep and cattle in Inner Mongolia, northeastern China. Parasit Vectors. 2019; 12(1): 450.

[16] Kong Y, Yan C, Liu D, Jiang L, Zhang G, He B, et al. Phylogenetic analysis of Crimean-Congo hemorrhagic fever virus in inner Mongolia, China. Ticks Tick Borne Dis. 2022; 13(1): 101856.

[17] Li Y, Yan C, Liu D, He B, Tu C. Seroepidemiological Investigation of Crimean-Congo hemorrhagic fever virus in sheep and camels of Inner Mongolia of China. Vector Borne Zoonotic Dis. 2020; 20(6): 461-467.

[18] Qin H, Xin X, Tang Q, Feng X, Yin B. The Prevalence of tick-borne encephalitis virus in the ticks and humans of China from 2000 to 2023: a systematic review and meta-analysis. Vet Sci. 2025; 12: 146.

[19] Li CX, Shi M, Tian JH, Lin XD, Kang YJ, Chen LJ, et al. Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. Elife. 2015; 4: e05378.

[20] Tokarz R, Williams SH, Sameroff S, Sanchez Leon M, Jain K, Lipkin WI. Virome analysis of Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis ticks reveals novel highly divergent vertebrate and invertebrate viruses . J Virol. 2014; 88(19): 11480-11492.

[21] Cross ST, Kapuscinski ML, Perino J, Maertens BL, Weger-Lucarelli J, Ebel GD, et al. Co-infection patterns in individual ixodes scapularis ticks reveal associations between viral, eukaryotic and bacterial microorganisms. Viruses. 2018; 10(7): 388.

[22] Ma R, Li C, Tian H, Zhang Y, Feng X, Li J, et al. The current distribution of tick species in Inner Mongolia and inferring potential suitability areas for dominant tick species based on the MaxEnt model. Parasit Vectors. 2023; 16(1): 286.

[23] Zhang Y, Hu B, Agwanda B, Fang Y, Wang J, Kuria S, et al. Viromes and surveys of RNA viruses in camel-derived ticks revealing transmission patterns of novel tick-borne viral pathogens in Kenya. Emerg Microbes Infect. 2021; 10(1): 1975-1987.

[24] Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011; 29(7): 644-652.

[25] Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie2. Nat Methods. 2012; 9(4): 357-359.

[26] Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009; 10: 421.

[27] Ni XB, Cui XM, Liu JY, Ye RZ, Wu YQ, Jiang JF, et al. Metavirome of 31 tick species provides a compendium of 1,801 RNA virus genomes. Nat Microbiol. 2023; 8(1): 162-173.

[28] Wille M, Harvey E, Shi M, Gonzalez-Acuña D, Holmes EC, Hurt AC. Sustained RNA virome diversity in Antarctic penguins and their ticks. ISME J. 2020; 14(7): 1768-1782.

[29] Pettersson JH, Ellström P, Ling J, Nilsson I, Bergström S, González-Acuña D, et al. Circumpolar diversification of the Ixodes uriae tick virome . PLoS Pathog. 2020; 16(8): e1008759.

[30] Simmonds P, Becher P, Bukh J, Gould EA, Meyers G, Monath T, et al. ICTV virus taxonomy profile: Flaviviridae. J Gen Virol. 2017; 98(1): 2-3.

[31] Matsuno K, Kajihara M, Nakao R, Nao N, Mori-Kajihara A, Muramatsu M, et al. The unique phylogenetic position of a novel tick-borne phlebovirus ensures an ixodid origin of the genus Phlebovirus. mSphere. 2018; 3(3): e00239-18.

[32] Longdon B, Murray GG, Palmer WJ, Day JP, Parker DJ, Welch JJ, et al. The evolution, diversity, and host associations of rhabdoviruses. Virus Evol. 2015; 1(1): vev014.

[33] Pettersson JH, Shi M, Bohlin J, Eldholm V, Brynildsrud OB, Paulsen KM, et al. Characterizing the virome of Ixodes ricinus ticks from northern Europe . Sci Rep. 2017; 7(1): 10870.

[34] Černý J, Buyannemekh B, Needham T, Gankhuyag G, Oyuntsetseg D. Hard ticks and tick-borne pathogens in Mongolia - a review. Ticks Tick Borne Dis. 2019; 10(6): 101268.

[35] Kholodilov IS, Belova OA, Morozkin ES, Litov AG, Ivannikova AY, Makenov MT, et al. Geographical and tick-dependent distribution of flavi-like Alongshan and Yanggou tick viruses in Russia. Viruses. 2021; 13(3): 458.

[36] Ma C, Zhang R, Zhou H, Yu G, Yu L, Li J, et al. Prevalence and genetic diversity of Dabieshan tick virus in Shandong Province, China. J Infect. 2022; 85(1): 90-122.

[37] Ye RZ, Li YY, Xu DL, Wang BH, Wang XY, Zhang MZ, et al. Virome diversity shaped by genetic evolution and ecological landscape of Haemaphysalis longicornis. Microbiome. 2024; 12(1): 35.

[38] Di Paola N, Dheilly NM, Junglen S, Paraskevopoulou S, Postler TS, Shi M, et al. Jingchuvirales: a new taxonomical framework for a rapidly expanding order of unusual monjiviricete viruses broadly distributed among Arthropod subphyla. Appl Environ Microbiol. 2022; 88(6): e0195421.

[39] Temmam S, Chrétien D, Bigot T, Dufour E, Petres S, Desquesnes M, et al. Monitoring silent spillovers before emergence: a pilot study at the tick/human interface in Thailand. Front Microbiol. 2019; 10: 2315.

[40] Tufts DM, Sameroff S, Tagliafierro T, Jain K, Oleynik A, VanAcker MC, et al. A metagenomic examination of the pathobiome of the invasive tick species, Haemaphysalis longicornis, collected from a New York City borough, USA. Ticks Tick Borne Dis. 2020; 11(6): 101516.

[41] Harvey E, Rose K, Eden JS, Lo N, Abeyasuriya T, Shi M, et al. Extensive diversity of RNA viruses in Australian ticks. J Virol. 2019; 93(3): e01358-18.

[42] Li X, Ji H, Wang D, Che L, Zhang L, Li L, et al. Molecular detection and phylogenetic analysis of tick-borne encephalitis virus in ticks in northeastern China. J Med Virol. 2022; 94(2): 507-513.

[43] Sun MH, Ji YF, Li GH, Shao JW, Chen RX, Gong HY, et al. Highly adaptive Phenuiviridae with biomedical importance in multiple fields. J Med Virol. 2022; 94(6): 2388-2401.

[44] Liu Z, Li L, Xu W, Yuan Y, Liang X, Zhang L, et al. Extensive diversity of RNA viruses in ticks revealed by metagenomics in northeastern China. PLoS Negl Trop Dis. 2022; 16(12): e0011017.

[45] Wang YN, Jiang RR, Ding H, Zhang XL, Wang N, Zhang YF, et al. First detection of Mukawa Virus in Ixodes persulcatus and Haemaphysalis concinna in China . Front Microbiol. 2022; 13: 791563.

[46] Ergönül Ö. Crimean-Congo haemorrhagic fever. Lancet Infect Dis. 2006; 6(4): 203-214.

[47] Xia H, Hu C, Zhang D, Tang S, Zhang Z, Kou Z, et al. Metagenomic profile of the viral communities in Rhipicephalus spp. ticks from Yunnan, China. PLoS One. 2015; 10(3): e0121609.