Extensive diversity of RNA viruses in ticks revealed by metagenomics in northeastern China

Exploring the potential of tick transcriptomes for virus screening: A data reuse approach for tick-borne virus surveillance

BACKGROUND

We set out to investigate the utility of publicly available tick transcriptomic data to identify and characterize known and recently described tick-borne viruses, using de novo assembly and subsequent protein database alignment and taxonomical binning.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 127 virus contigs were recovered from 35 transcriptomes, originating from cell lines (40%), colony-reared ticks (25.7%) or field-collected ticks (34.2%). Generated virus contigs encompass DNA (n = 2) and RNA (n = 13) virus families, with 3 and 28 taxonomically distinct isolates, respectively. Known human and animal pathogens comprise 32.8% of the contigs, where Beiji nairovirus (BJNV) was the most prevalent tick-borne pathogenic virus, identified in 22.8% of the transcriptomes. Other pathogens included Nuomin virus (NUMV) (2.8%), African swine fever virus (ASFV) (5.7%), African horse sickness virus 3 (AHSV-3) (2.8%) and Alongshan virus (ALSV) (2.8%).

CONCLUSIONS

Previously generated transcriptome data can be leveraged for detecting tick-borne viruses, as exemplified by new descriptions of ALSV and BJNV in new geographic locations and other viruses previously detailed in screening reports. Monitoring pathogens using publicly available data might facilitate biosurveillance by directing efforts to regions of preliminary spillover and identifying targets for screening. Metadata availability is crucial for further assessments of detections.

Development of a loop-mediated isothermal amplification (LAMP) assay for the rapid detection of Beiji nairovirus

Beiji nairovirus (BJNV) is a recently discovered tick-borne RNA virus associated with human febrile illness. This study aimed to develop a loop-mediated isothermal amplification (LAMP) assay for the precise detection of BJNV, with a specific focus on assessing its effectiveness with clinical samples. The optimal molecular target was identified as the BJNV small (S) segment gene, and the ideal reaction conditions were established at 65 °C for 50 min. A neutral red stain concentration of 300 μM was determined to be optimal for visualizing the LAMP reaction. The LAMP method demonstrated an impressive lower limit of detection at 10 copies/μL, highlighting a sensitivity level 10,000 times higher than traditional PCR methods. Moreover, it showed no cross-reactivity with three other viruses. The LAMP assay demonstrated results consistent with those of semi-nested PCR when applied to clinical samples, thereby validating its suitability for field testing. In conclusion, the LAMP assay developed in this study represents a significant advancement in the rapid detection of BJNV. Its high sensitivity, specificity, and ease of use make it a valuable tool for the establishment of effective prevention and control strategies against this emerging pathogen.

Detection and characterization of bovine hepacivirus in cattle and sheep from Hulunbuir, northeastern China

Background

Bovine hepacivirus (BovHepV) is a recently identified member of the expanding genus Hepacivirus within the family Flaviviridae. However, the genetic diversity, geographical distribution, and host range of the virus remains poorly understood.

Methods

In this study, serum samples from cattle and sheep were obtained in Hulunbuir and pooled to establish RNA libraries, which were then analyzed using transcriptome sequencing. BovHepV-positive samples were confirmed using semi-nested PCR with primers designed based on the obtained viral sequences. Comprehensive bioinformatics analyses were employed to assess sequence similarity, phylogenetic evolution, and recombination of the obtained viral strains.

Results

A total of 988 serum samples from sheep (520) and cattle (468) were collected from 12 administrative districts in Hulunbuir from June to August, 2022. Semi-nested PCR revealed 6 BovHepV-positive districts with prevalence ranging from 2.0% to 35.0% in cattle, and one BovHepV-positive district with prevalence of 2.5% in sheep. The nucleotide sequence identities between viral strains from sheep and cattle ranged from 91.3% to 93.8%, while the amino acid sequence identities were between 95.4% and 96.7%. Phylogenetic analyses classified the obtained BovHepV strains within subtype G, genotype 1. Recombination analysis revealed the intergenerational relationships among the viral strains obtained from cattle and sheep.

Conclusion

We identified genetic diversity in subtype G strains in cattle and detected a BovHepV strain in a sheep for the first time in northeastern China, confirming cross-species transmission and co-circulation between cattle and sheep, thus expanding the virus's host range.

Metatranscriptomic analysis reveals the diversity of RNA viruses in ticks in Inner Mongolia, China

BACKGROUND

Ticks are widely distributed throughout China and are the second most prevalent pathogen vectors in the world, following only mosquitoes. Tick bites can lead to Lyme disease, forest encephalitis, and other illnesses that may result in death under severe circumstances. Materials and methods: Ticks collected from March 2021 to May 2023 were pooled and used in metatranscriptomic analyses to gain insight into the diversity and distribution of tick-borne viruses in Inner Mongolia. Next-generation sequencing (NGS) outcomes were validated, and viral prevalence across distinct tick species was determined through the application of polymerase chain reaction (PCR) paired with Sanger sequencing.

RESULTS

A total of 20 RNA viruses belonging to at least 8 families, including Chuviridae, Flaviviridae, Solemoviridae, Nairoviridae, Partitiviridae, Phenuiviridae, Rhabdoviridae, and Totiviridae, and to unclassified families were identified by NGS. Five of the identified RNA viruses (Nuomin virus, Yezo virus, tick-borne encephalitis virus, Alongshan virus, and Beiji nairovirus) are considered human pathogens. A potential human pathogen, Mukawa virus, was also among the identified viruses. Ixodes persulcatus carried a significantly greater number of viral species than did Dermacentor nuttalli, Hyalomma marginatum, and Haemaphysalis concinna. The prevalence of coinfection with multiple viruses differed in I. persulcatus from Hinggan League and Hulun Buir, and Beiji nairovirus was the codominant virus species.

CONCLUSIONS

There is a remarkable diversity of RNA viruses harboured by ticks in Inner Mongolia, with variations observed in the distribution of these tick-borne viruses across different regions and tick hosts.

Identification and characterization of Jingmen tick virus from Rhipicephalus microplus in Hunan, China

Jingmen tick virus (JMTV) is a tick-borne pathogen known to affect human beings, characterized by a segmented genome structure that defies the conventional understanding of the Flaviviridae family. In the present study, we employed metagenomic analysis to screen for tick-borne viruses in Hunan Province, China, and identified five JMTV variants with complete genomes from Rhipicephalus microplus ticks sampled from cattle. These viral strains exhibited the highest sequence similarity to JMTV isolates previously reported in Hubei Province, China. However, evidence of genomic reassortment was detected, particularly with the S2 segment showing greater similarity to the strains from Japan. Phylogenetic analysis demonstrated that JMTV strains cluster predominantly based on their geographic origin. In agreement with the homology data, the S1, S3, and S4 segments of the strains identified in this study grouped with those from Hubei Province, while the S2 segment displayed a distinct topological structure. Moreover, JMTV displayed limited replication in mammal-derived cells, but thrived in tick-derived cell lines. In addition to the commonly used R. microplus-derived BME/CTVM23 cells, we found that JMTV also proliferated robustly in both Ixodes scapularis-derived ISE6 and Ixodes ricinus-derived IRE/CTVM19 cells, offering new avenues for in vitro production of the virus. In summary, this study expands the known geographic distribution and genetic diversity of JMTV, providing valuable insights into its epidemiology and potential for in vitro cultivation.

Development of a LAMP assay for the rapid visual detection of the emerging tick-borne Songling virus

BACKGROUND

Songling virus (SGLV) within the genus Orthonairovirus, family Nairoviridae, is an emerging tick-borne virus associated with human febrile illness. However, no rapid detection method for SGLV has been established.

METHODS

In this study, four primer sets targeting the nucleocapsid protein gene of SGLV were designed for use in the LAMP assay and evaluated to identify the optimal primer set. Recombinant plasmids were constructed and utilized for assessing the sensitivity of the assay. Tacheng tick virus 1 (TcTV-1)-, Beiji nairovirus (BJNV)-, Yezo virus (YEZV)-, severe fever with thrombocytopenia syndrome virus (SFTSV)-, and tick-borne encephalitis virus (TBEV)-positive tick samples were utilized to assess the specificity. Field-collected ticks were also evaluated as biological specimens to validate the assay.

RESULTS

A SGLV-specific LAMP assay was established with a detection limit of 1 × 10-2 copies/μl and could be visually confirmed by a color change from purple to blue in SGLV-positive samples. No cross-reactivity was observed in the detection of TcTV-1, BJNV, YEZV, SFTSV, and TBEV using the LAMP assay. In addition to the detection of the same seven high-copy numbers of SGLV as the SYBR Green quantitative RT-PCR assay within a reduced timeframe, the developed LAMP method also effectively identified an additional sample with a low copy number in the field-collected tick samples.

CONCLUSIONS

We successfully developed a sensitive, specific, and cost-effective visual method for the rapid detection of SGLV using the LAMP assay, which can be applied in pathogenesis and epidemiological surveillance studies of SGLV.

A wide distribution of Beiji nairoviruses and related viruses in Ixodes ticks in Japan

Beiji nairovirus (BJNV), in the family Nairoviridae, the order Bunyavirales, was recently reported as a causative agent of an emerging tick-borne zoonotic infection in China. This study investigated the prevalence of BJNV in ticks in Japan. Screening of over 2,000 ticks from multiple regions revealed a widespread distribution of BJNV and BJNV-related viruses in Japan, particularly in the northern island, and in other high altitude areas with exclusive occurrence of Ixodes ticks. Phylogenetic analysis identified three distinct groups of nairoviruses in ticks in Japan: BJNV, Yichun nairovirus (YCNV) and a newly identified Mikuni nairovirus (MKNV). BJNV and YCNV variants identified in ticks in Japan exhibited high nucleotide sequence identities to those in China and Russia with evidence of non-monophyletic evolution among BJNVs, suggesting multiple cross-border transmission events of BJNV between the Eurasian continent and Japan. Whole genome sequencing of BJNV and MKNV revealed a unique GA-rich region in the S segment, the significance of which remains to be determined. In conclusion, the present study has shown a wide distribution and diversity of BJNV-related nairoviruses in Ixodes ticks in Japan and has identified unique genomic structures. The findings demonstrate the significance of BJNV as well as related viruses in Japan and highlight the necessity of monitoring emerging nairovirus infections and their potential risks to public health.

A New Orthonairovirus Associated with Human Febrile Illness

BACKGROUND

In June 2019, a patient presented with persistent fever and multiple organ dysfunction after a tick bite at a wetland park in Inner Mongolia. Next-generation sequencing in this patient revealed an infection with a previously unknown orthonairovirus, which we designated Wetland virus (WELV).

METHODS

We conducted active hospital-based surveillance to determine the prevalence of WELV infection among febrile patients with a history of tick bites. Epidemiologic investigation was performed. The virus was isolated, and its infectivity and pathogenicity were investigated in animal models.

RESULTS

WELV is a member of the orthonairovirus genus in the Nairoviridae family and is most closely related to the tickborne Hazara orthonairovirus genogroup. Acute WELV infection was identified in 17 patients from Inner Mongolia, Heilongjiang, Jilin, and Liaoning, China, by means of reverse-transcriptase-polymerase-chain-reaction assay. These patients presented with nonspecific symptoms, including fever, dizziness, headache, malaise, myalgia, arthritis, and back pain and less frequently with petechiae and localized lymphadenopathy. One patient had neurologic symptoms. Common laboratory findings were leukopenia, thrombocytopenia, and elevated d-dimer and lactate dehydrogenase levels. Serologic assessment of convalescent-stage samples obtained from 8 patients showed WELV-specific antibody titers that were 4 times as high as those in acute-phase samples. WELV RNA was detected in five tick species and in sheep, horses, pigs, and Transbaikal zokors (Myospalax psilurus) sampled in northeastern China. The virus that was isolated from the index patient and ticks showed cytopathic effects in human umbilical-vein endothelial cells. Intraperitoneal injection of the virus resulted in lethal infections in BALB/c, C57BL/6, and Kunming mice. The Haemaphysalis concinna tick is a possible vector that can transovarially transmit WELV.

CONCLUSIONS

A newly discovered orthonairovirus was identified and shown to be associated with human febrile illnesses in northeastern China. (Funded by the National Natural Science Foundation of China and the Chinese Academy of Medical Sciences Innovation Fund for Medical Sciences.).

Metagenome reveals the midgut microbial community of Haemaphysalis qinghaiensis ticks collected from yaks and Tibetan sheep

BACKGROUND

Haemaphysalis qinghaiensis is a tick species distributed only in China. Due to its ability to transmit a variety of pathogens, including species of the genera Anaplasma, Rickettsia, Babesia, and Theileria, it seriously endangers livestock husbandry. However, the microbial community of the midgut of H. qinghaiensis females collected from yaks and Tibetan sheep has not yet been characterized using metagenomic sequencing technology.

METHODS

Haemaphysalis qinghaiensis were collected from the skins of yaks and Tibetan sheep in Gansu Province, China. Genomic DNA was extracted from the midguts and midgut contents of fully engorged H. qinghaiensis females collected from the two hosts. Metagenomic sequencing technology was used to analyze the microbial community of the two groups.

RESULTS

Fifty-seven phyla, 483 genera, and 755 species were identified in the two groups of samples. The ticks from the two hosts harbored common and unique microorganisms. At the phylum level, the dominant common phyla were Proteobacteria, Firmicutes, and Mucoromycota. At the genus level, the dominant common genera were Anaplasma, Ehrlichia, and Pseudomonas. At the species level, bacteria including Anaplasma phagocytophilum, Ehrlichia minasensis, and Pseudomonas aeruginosa along with eukaryotes such as Synchytrium endobioticum and Rhizophagus irregularis, and viruses such as the orf virus, Alphadintovirus mayetiola, and Parasteatoda house spider adintovirus were detected in both groups. In addition, the midgut of H. qinghaiensis collected from yaks had unique microbial taxa including two phyla, eight genera, and 23 species. Unique microorganisms in the midgut of H. qinghaiensis collected from Tibetan sheep included two phyla, 14 genera, and 32 species. Kyoto Encyclopedia of Genes and Genomes enrichment analysis revealed that the functional genes of the microbiome of H. qinghaiensis were annotated to six pathways, and the metabolic pathways included 11 metabolic processes, in which the genes involved in carbohydrate metabolism were the most abundant, followed by the genes involved in lipid metabolism.

CONCLUSIONS

These findings indicate that most of the microbial species in the collected H. qinghaiensis ticks were the same in both hosts, but there were also slight differences. The analytical data from this study have enhanced our understanding of the midgut microbial composition of H. qinghaiensis collected from different hosts. The database of H. qinghaiensis microbe constructed from this study will lay the foundation for predicting tick-borne diseases. Furthermore, a comprehensive understanding of tick microbiomes will be useful for understanding vector competency and interactions with ticks and midgut microorganisms.

Viral Diversity and Ecological Impact of DNA Viruses in Dominant Tick Species in China

Ticks are blood-feeding ectoparasites that also transmit various pathogens, posing severe risks to human and animal health. DNA viruses play a crucial role in the microbial ecology of ticks, but their distribution and ecological significance remain largely undetermined. Here, we assembled an extensive catalog encompassing 4320 viral operational taxonomic units (vOTUs) from six main dominant tick species in China, of which 94.8% have not been found in any other environment. To bridge the knowledge gap in tick DNA virus research and provide a crucial resource platform, we developed the Tick DNA Virus Database. This database includes the vOTUs that are known to cause diseases. Most of the predicted vOTUs are associated with dominant bacterial and archaeal phyla. We identified 105 virus-encoded putative auxiliary metabolic genes (AMGs) that are involved in host metabolism and environmental adaptation, potentially influencing ticks through both top-down and bottom-up mechanisms. The identification of microbial communities and antibiotic resistance in wild tick species suggests that wild ticks are reservoirs of antibiotic resistance and potential spreaders of antibiotic resistance. These findings reveal the potential role of tick viruses in ecosystems, highlighting the importance of monitoring tick microbiomes to address global public health challenges.

Alongshan Virus Infection in Rangifer tarandus Reindeer, Northeastern China

We investigated Alongshan virus infection in reindeer in northeastern China. We found that 4.8% of the animals were viral RNA-positive, 33.3% tested positive for IgG, and 19.1% displayed neutralizing antibodies. These findings suggest reindeer could serve as sentinel animal species for the epidemiologic surveillance of Alongshan virus infection.

Xinyang flavivirus, from Haemaphysalis flava ticks in Henan Province, China, defines a basal, likely tick-only Orthoflavivirus clade

Tick-borne orthoflaviviruses (TBFs) are classified into three conventional groups based on genetics and ecology: mammalian, seabird and probable-TBF group. Recently, a fourth basal group has been identified in Rhipicephalus ticks from Africa: Mpulungu flavivirus (MPFV) in Zambia and Ngoye virus (NGOV) in Senegal. Despite attempts, isolating these viruses in vertebrate and invertebrate cell lines or intracerebral injection of newborn mice with virus-containing homogenates has remained unsuccessful. In this study, we report the discovery of Xinyang flavivirus (XiFV) in Haemaphysalis flava ticks from Xìnyáng, Henan Province, China. Phylogenetic analysis shows that XiFV was most closely related to MPFV and NGOV, marking the first identification of this tick orthoflavivirus group in Asia. We developed a reverse transcriptase quantitative PCR assay to screen wild-collected ticks and egg clutches, with absolute infection rates of 20.75 % in adult females and 15.19 % in egg clutches, suggesting that XiFV could be potentially spread through transovarial transmission. To examine potential host range, dinucleotide composition analyses revealed that XiFV, MPFV and NGOV share a closer composition to classical insect-specific orthoflaviviruses than to vertebrate-infecting TBFs, suggesting that XiFV could be a tick-only orthoflavivirus. Additionally, both XiFV and MPFV lack a furin cleavage site in the prM protein, unlike other TBFs, suggesting these viruses might exist towards a biased immature particle state. To examine this, chimeric Binjari virus with XIFV-prME (bXiFV) was generated, purified and analysed by SDS-PAGE and negative-stain transmission electron microscopy, suggesting prototypical orthoflavivirus size (~50 nm) and bias towards uncleaved prM. In silico structural analyses of the 3'-untranslated regions show that XiFV forms up to five pseudo-knot-containing stem-loops and a prototypical orthoflavivirus dumbbell element, suggesting the potential for multiple exoribonuclease-resistant RNA structures.

Diversity of species and geographic distribution of tick-borne viruses in China

Introduction

Tick-borne pathogens especially viruses are continuously appearing worldwide, which have caused severe public health threats. Understanding the species, distribution and epidemiological trends of tick-borne viruses (TBVs) is essential for disease surveillance and control.

Methods

In this study, the data on TBVs and the distribution of ticks in China were collected from databases and literature. The geographic distribution of TBVs in China was mapped based on geographic locations of viruses where they were prevalent or they were detected in vector ticks. TBVs sequences were collected from The National Center for Biotechnology Information and used to structure the phylogenetic tree.

Results

Eighteen TBVs from eight genera of five families were prevalent in China. Five genera of ticks played an important role in the transmission of TBVs in China. According to phylogenetic analysis, some new viral genotypes, such as the Dabieshan tick virus (DTV) strain detected in Liaoning Province and the JMTV strain detected in Heilongjiang Province existed in China.

Discussion

TBVs were widely distributed but the specific ranges of viruses from different families still varied in China. Seven TBVs belonging to the genus Orthonairovirus of the family Nairoviridae such as Nairobi sheep disease virus (NSDV) clustered in the Xinjiang Uygur Autonomous Region (XUAR) and northeastern areas of China. All viruses of the family Phenuiviridae except Severe fever with thrombocytopenia syndrome virus (SFTSV) were novel viruses that appeared in the last few years, such as Guertu virus (GTV) and Tacheng tick virus 2 (TcTV-2). They were mainly distributed in the central plains of China. Jingmen tick virus (JMTV) was distributed in at least fourteen provinces and had been detected in more than ten species of tick such as Rhipicephalus microplus and Haemaphysalis longicornis, which had the widest distribution and the largest number of vector ticks among all TBVs. Parainfluenza virus 5 (PIV5) and Lymphatic choriomeningitis virus (LCMV) were two potential TBVs in Northeast China that could cause serious diseases in humans or animals. Ixodes persulcatus carried the highest number of TBVs, followed by Dermacentor nuttalli and H. longicornis. They could carry as many as ten TBVs. Three strains of Tick-borne encephalitis (TBEV) from Inner Mongolia Province clustered with ones from Russia, Japan and Heilongjiang Province, respectively. Several SFTSV strains from Zhejiang Province clustered with strains from Korea and Japan. Specific surveillance of dominant TBVs should be established in different areas in China.

Characterisation of putative novel tick viruses and zoonotic risk prediction

Tick-associated viruses remain a substantial zoonotic risk worldwide, so knowledge of the diversity of tick viruses has potential health consequences. Despite their importance, large amounts of sequences in public data sets from tick meta-genomic and -transcriptomic projects remain unannotated, sequence data that could contain undocumented viruses. Through data mining and bioinformatic analysis of more than 37,800 public meta-genomic and -transcriptomic data sets, we found 83 unannotated contigs exhibiting high identity with known tick viruses. These putative viral contigs were classified into three RNA viral families (Alphatetraviridae, Orthomyxoviridae and Chuviridae) and one DNA viral family (Asfarviridae). After manual checking of quality and dissimilarity towards other sequences in the data set, these 83 contigs were reduced to five contigs in the Alphatetraviridae from four putative viruses, four in the Orthomyxoviridae from two putative viruses and one in the Chuviridae which clustered with known tick-associated viruses, forming a separate clade within the viral families. We further attempted to assess which previously known tick viruses likely represent zoonotic risks and thus deserve further investigation. We ranked the human infection potential of 133 known tick-associated viruses using a genome composition-based machine learning model. We found five high-risk tick-associated viruses (Langat virus, Lonestar tick chuvirus 1, Grotenhout virus, Taggert virus and Johnston Atoll virus) that have not been known to infect human and two viral families (Nairoviridae and Phenuiviridae) that contain a large proportion of potential zoonotic tick-associated viruses. This adds to the knowledge of tick virus diversity and highlights the importance of surveillance of newly emerging tick-associated diseases.

Spatial risk of Haemaphysalis longicornis borne Dabieshan tick virus (DBTV) in China

The uncertainty and unknowability of emerging infectious diseases have caused many major public health and security incidents in recent years. As a new tick-borne disease, Dabieshan tick virus (DBTV) necessitate systematic epidemiological and spatial distribution analysis. In this study, tick samples from Liaoning Province were collected and used to evaluate distribution of DBTV in ticks. Outbreak points of DBTV and the records of the vector Haemaphysalis longicornis in China were collected and used to establish a prediction model using niche model combined with environmental factors. We found that H. longicornis and DBTV were widely distributed in Liaoning Province. The risk analysis results showed that the DBTV in the eastern provinces of China has a high risk, and the risk is greatly influenced by elevation, land cover, and meteorological factors. The risk geographical area predicted by the model is significantly larger than the detected positive areas, indicating that the etiological survey is seriously insufficient. This study provided molecular and important epidemiological evidence for etiological ecology of DBTV. The predicted high-risk areas indicated the insufficient monitoring and risk evaluation and the necessity of future monitoring and control work.