Skip to main content

DNA virome of ticks in the Northeast and Hubei provinces of China reveals diverse single-stranded circular DNA viruses



Ticks are medically important vectors capable of transmitting a variety of pathogens to and between host species. Although the spectrum of tick-borne RNA viruses has been frequently investigated, the diversity of tick-borne DNA viruses remains largely unknown.


A total of 1571 ticks were collected from forests and infested animals, and the diversity of the viruses they harbored was profiled using a DNA-specific virome method. The viromic data were phylogenetically analyzed and validated by PCR assays.


Although diverse and abundant prokaryotic viruses were identified in the collected ticks, only eukaryotic DNA viruses with single-stranded circular genomes covering the anelloviruses and circular replication-associated (Rep) protein-encoding single-stranded (CRESS) DNA viruses were recovered from ticks. Anelloviruses were detected only in two tick pools, but CRESS DNA viruses were prevalent across these ticks except in one pool of Dermacentor spp. ticks. Phylogenetic analyses revealed that these tick-borne CRESS DNA viruses were related to viruses recovered from animal feces, tissues and even environmental samples, suggesting that their presence may be largely explained by environmental factors rather than by tick species and host blood meals.


Based on the results, tick-borne eukaryotic DNA viruses appear to be much less common than eukaryotic RNA viruses. Investigations involving a wider collection area and more diverse tick species are required to further support this speculation.

Graphical Abstract


Ticks are among the most important arthropod vectors capable of transmitting a wide range of microorganisms in humans, wildlife and domestic animals, including viruses, bacteria, protozoa, fungi and nematodes [1]. Due to their obligate propensity for feeding on the blood of different host species, ticks maintain the complex circulation of tick-borne pathogens in vertebrate hosts, seriously hindering the control and prevention of tick-borne diseases (TBDs) [2]. The changes in climate that are currently ongoing have significantly impacted the life-cycle of ticks, resulting in lower mortality during overwintering and a longer active period [3]. In addition, because the interfaces between tick habitats and human settlements are frequently disturbed by human activity [4], TBDs, especially those caused by a viral agent, are increasing; for example a number of new tick-borne viruses (TBVs) have been identified as causative agents of human febrile illnesses in northeast (NE) China over the past few years [5,6,7,8]. Therefore, continuous investigation and surveillance of genetically diverse TBVs is of significant importance for the prevention and control of emerging TBDs.

TBVs were first identified more than a century ago, but most studies have been confined to RNA viruses that are classified into two orders, eight families and at least 11 genera. It is known that lethal TBDs are predominantly caused by TBVs that are RNA viruses (RNA TBVs). For example, tick-borne encephalitis (TBE) and Omsk hemorrhagic fever are caused by flaviviruses [9, 10], while Crimean-Congo hemorrhagic fever virus (CCHFV) and severe fever with thrombocytopenia syndrome virus (SFTSV) are responsible for the deadly hemorrhagic fever diseases [11, 12]. As a result, RNA TBVs have received more attention, with tick viromic profiling focusing primarily on RNA viruses [13, 14]. Conversely, TBVs that are DNA viruses (DNA TBVs) have been relatively neglected, even though two well-known DNA viruses, namely African swine fever virus (ASFV) and lumpy skin disease virus (LSDV), are associated with ticks and cause serious diseases with significant economic consequences in the livestock industry. The former is transmitted by soft ticks from the genus Ornithodoros [15], while the latter belongs to the genus Capripoxvirus of the family Poxviridae [16] and is transmitted by Rhipicephalus appendiculatus and Amblyomma hebraeum ticks [17]. In addition, a few studies have identified several distinct small DNA viruses, such as circoviruses from Haemaphysalis longicornis and Haemaphysalis crenulatus ticks [18], genomoviruses and anelloviruses from Dermacentor variabilis and Ixodes scapularis ticks [19] and parvoviruses from Ixodes ricinus and Rhipicephalus sanguineus ticks [20, 21]. Also, murine gammaherpesvirus 68, a natural pathogen of murid rodents, was also identified in Dermacentor reticulatus and I. ricinus ticks, which has expanded the spectrum of DNA TBVs to the family Herpesviridae [22, 23].

The sporadic studies carried out on DNA TBVs are insufficient to capture a complete view of the diversity of DNA TBVs. To further understand the DNA viral spectrum harbored by ticks, we used a DNA-specific virome method to investigate the DNA virome harbored by free-living and engorged ticks collected from the NE province and Hubei province in China, respectively. The method used here employs the multiple displacement amplification (MDA) technique to magnify viral DNA and can capture DNA viruses of all genome types [24,25,26]. The two regions of this study are known to be natural habitats of abundant ticks and hotspots of many TBDs, such as TBE and SFTS [27,28,29]. The results revealed diverse tick-borne small circular DNA viruses, providing insight into the diversity of DNA TBVs.


Tick sampling and species identification

The most favorable climatic conditions for tick proliferation are during the spring and summer. Consequently, we collected ticks by dragging and flagging vegetation in forests adjacent to human settlements in Heilongjiang, Jilin and Liaoning provinces between May and June 2021. We also directly plucked ticks from infested goats in Hubei Province in May 2016 (Fig. 1a, b). Ticks were placed immediately upon collection in labeled vials and transported on dry ice to our laboratory. Adult ticks were initially identified to the species level using tick taxonomic keys by a trained expert under a stereo microscope (Nikon model H550S; Nikon Corp., Tokyo, Japan) based various morphological characters, such as the scutum, the basis capitula and palp, the anal groove and the shape and size of the spurs on the coxae [30,31,32,33]. Species identification was further confirmed by sequencing the mitochondrial (mt) cytochrome c oxidase subunit I gene (COI) [34], which was also used to identify the species of larval ticks.

Fig. 1
figure 1

Sample information and virome overview. a Geographical distribution of tick sampling sites. HLJ, Heilongjiang province; JL, Jilin province; LN, Liaoning province; HB, Hubei province; TH, Tahe; XK, Xunke; JY, Jiayin; LB, Luobei; DN, Dongning; HN, Huinan; LJ, Linjiang; HL, Helong; YJ, Yanji; DH, Dunhua; TM, Tumen; HC, Hunchun; ZH, Zhuanghe; ZA, Zhenan; HR, Huanren; XS, Xinshi; LS, Liusheng; YD, Yindian. b Numbers of ticks sampled at each site. c Overview of tick DNA virome. Left panel: Heat map of viral reads from 20 libraries. Right panel: amino acid identity range of these virus-like contigs with reference sequences in Genbank. Each box plot illustrates the estimated median (center line), and upper and lower quartiles (box limits) of the similarity. The identifier of the library was composed of four letters and a two-digit number, whereby the first two letters are abbreviations for tick species, the next two letters are indicative of the sampling sites and the number refers to the pool’s number. Euk, Eukaryotic; Prok, prokaryotic; Un, unclassified. d K-mer distance-based clustering of these pools. Pools with the same tick species are identified using the same color

Sample pretreatment and high-throughput sequencing

All ticks were grouped into pools according to sampling locations and species, with each pool containing 13–150 adult ticks or 800 larvae (Fig. 1c). Each pool’s name provides information on tick species (the first 2 letters), location (the second 2 letters) and pool number (the last 2 numbers). Ticks in each pool were first washed 3 times with sterile phosphate-buffered saline solution (PBS; pH = 7.4) and then homogenized using a ground-glass tissue grinder in 1 ml of sterile magnesium salt buffer (50 mM Tris, 10 mM MgSO4, 0.1 M NaCl, pH 7.5). Homogenates were centrifuged at 8000 g for 10 min at 4 °C to remove debris. The supernatants were filtered through a Millex filter (pore size: 0.45 μm) (MilliporeSigma, Burlington, MA, USA) and digested with 1 μl DNase I (TaKaRa, Dalian, China) and 1 μl RNase A (TaKaRa) at 37 °C for 60 min to eliminate host genomic DNA and other free nucleic acids. The total DNA was then extracted immediately using a DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The DNA-specific MDA method was employed to amplify the DNA using an Illustra GenomiPhi V2 DNA Amplification Kit (GE Healthcare, Chicago, IL, USA). A 1-μg sample of products was subjected to DNA library preparation using a NEB Next® Ultra™ DNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA), followed by paired-end (150-bp) sequencing on a NovaSeq sequencer (Illumina, Inc., San Diego, CA). At least 6 Gb data of each library were generated. Cross and foreign contamination are common in metagenomic sequencing [35]; six pools of bat rectums and double-distilled water were also processed synchronously to determine potential contamination events because ticks and bats belong to different classes and are thought to harbor distinct viromes [35].

Annotation and comparison of virome

The fastp v0.19.7 tool was used for the initial quality control of the raw data. Host sequences were removed from clean data by mapping against the whole genomic assemblies of H. longicorni (accession: GCA_013339765.2), Ixodes persulcatus (accession: GCA_013358835.1) and Dermacentor silvarum (accession: GCA_013339745.1) using the Bowtie2 v2.4.1 tool. Bacterial, archaeal and fungal reads in the remaining data were further subtracted using a fast taxonomic classification with Kraken2 (v2.0.9-beta) and a custom RefSeq-based database. The unclassified reads were de novo assembled separately using the metaSPAdes v3.14.9 metagenome assembler, and circular genomes were recovered using Recycler. Prior to annotation, tick viromic contigs were compared to those of bats using BLASTn. If an alignment was achieved with length > 500 bp and identity > 99%, the tick contig was considered to be a contaminant. This comparison procedure detected hundreds of tick contigs closely related to bat contigs. These sequences were finally validated to be eukaryotic sequences or genomic fragments of DNA viruses with low confidence after querying Genbank and, therefore, were deleted from the assembly. However, because the double-distilled water did not result in any nucleic acid products, no further action was taken. The remaining contigs were retained for BLASTn (v2.10.0) and diamond BLASTx (v0.9.35) searches against our refined eukaryotic viral reference database (EVRD)—nucleotides/amino acids (nt/aa) [35]. To investigate prokaryotic viruses, we also searched these contigs against the bacteriophage aa sequences of UniProt database. The read abundance was counted by mapping unclassified reads back to virus-like contigs using bowtie2 and SAMtools v1.10. To examine the viromic difference between libraries, we compared the unclassified data sets using MASH v2.2. A sketch file of each data set was created with sketch size of 10,000 and K-mer size of 21. These sketches were mutually compared to generate the MASH distances, based on which Ward’s minimum variance clustering was performed in the R environment v4.1.3 ® Foundation for Statistical Computing, Vienna, Austria).

PCR validation

The completeness of circular genome representatives was validated by PCR using partially overlapped back-to-back primer pairs. The DNA of the remaining supernatants of a particular pool was extracted as described above, from which the circular genome was recovered by PCR with Phanta® Max Super-Fidelity DNA Polymerase (Vazyme, Wuhan, China) and primer pairs with the following thermocycler program: pre-denaturation at 98 °C for 1 min; 35 cycles of denaturation at 98 °C for 5 s, annealing at 59 °C for 10 s and extension at 72 °C for 10 s; and a final extension at 72 °C for 5 min. The products were examined by electrophoresis in a 1% agarose gel with the expected products then directly sequenced on an ABI 3730xl DNA analyzer (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA). The sequence was further aligned against the initial circular genomes using MegAlign to check its accuracy.

Phylogenetic analyses

The putative open-reading frames (ORFs) of circular virus genomes were identified using NCBI ORF Finder ( References of the replication-associated protein genes (Rep) were withdrawn from Genbank after BLASTx search and aligned together with the Rep aa sequences identified here using the MAFFT L-INS-I algorithm. The alignments were used to infer maximum likelihood phylogenetic trees using IQ-TREE with the LG + G substitution model selected using ModelFinder. The resulting phylogenetic trees were mid-point rooted, and branches with < 70% bootstrap support were collapsed using the TreeGraph 2 graphical editor for phylogenetic trees.


Overview of tick-borne DNA virome

A total of 1571 ticks were collected from 18 locations in Heilongjiang, Jilin, Liaoning and Hubei provinces (Fig. 1a, b). Ticks from the first three provinces, all located in NE China, were all free-living and comprised adults (n = 591) and larvae (n = 800), while the ticks collected from Hubei Province were adults plucked from infested goats (n = 180) (Fig. 1a, b). These ticks were identified as Haemaphysalis concinna (n = 86), Haemaphysalis japonica (n = 864, including 800 larvae), H. longicornis (n = 480), I. persulcatus (n = 108) and Dermacentor spp. (n = 33) (Fig. 1b). They were grouped into 20 pools based on species and sampling location and subjected to DNA virome profiling, which revealed 875 virus-like contigs with 20 complete circular genomes. All 875 contigs were annotated to small circular DNA viruses, i.e. Anelloviridae, Circoviridae, Genomoviridae and unclassified viruses with single-stranded circular DNA (sscDNA) (Fig. 1c). The annotation using UniProt also uncovered diverse bacteriophages of different genomic types, including those with double-stranded linear and single-stranded circular genomes (Fig. 1c). These showed a wide range of aa identities to the reference sequences, with the majority distantly related to currently known viruses (Fig. 1c). Mapping of the reads back to contigs showed that virus-like reads composed 0 (DSHN01) to 54.5% (HCYJ01) of their total reads (Fig. 1c). Circoviridae and unclassified sscDNA viruses were widely distributed in these libraries with the highest abundances, while anelloviruses were only detected in two libraries with 0.9% and 0.1% relative abundances in HJDN01 and IPTM01, respectively (Fig. 1c). Interestingly, the larva pool LAHC01 also contained 1.6% reads related to Circoviridae, Genomoviridae and unclassified ssDNA viruses (Fig. 1C). Multiple tick species were involved in this study, providing an ideal model to inspect the DNA viromic composition feature of tick species. We abstracted the unassigned data into K-mer sketches using MASH and clustered these libraries based on the K-mer distances. Although one cluster was primarily composed of pools of H. longicornis ticks regardless of their sampling site or diet (Fig. 1d), the clustering of the remaining tick pools was not significantly correlated with tick species (Fig. 1d), implying that host species have little influence on the sequence composition feature of the tick DNA virome.

Characterization of circular Rep protein-encoding DNA molecules

Among the obtained 20 circular genomes, seven were 1105–1551 nt in length and only encoded a Rep protein (Fig. 2a). PCR amplification using partially overlapped back-to-back primer pairs and sequencing confirmed their circular structures (Fig. 2b). We termed these genomes tick-associated circular DNA molecules (TaCMs). The putative Rep proteins of these TaCMs were 101–309 aa long and showed a wide range of similarities to known sequences. Phylogenetic analyses showed that these Rep sequences were clustered into five clades within two major branches. Although circular DNA molecules in Genbank were among the genetic neighbors of some TaCMs, the closest relatives were all circular Rep-encoding single-stranded (CRESS) DNA viruses from a variety of sources, such as plants, bird feces and animal tissues (Fig. 2c). Most TaCM Rep were very divergent from their references, with 29.0–78.5% aa identities, whereas TaCM HJMX01-1 was highly similar to CRESS DNA viruses recovered from fecal metagenomes of Australian birds with up to 100% aa identities.

Fig. 2
figure 2

Genomic characterization of TaCMs and TaCVs. a Genomic structures of three TaCM and three TaCV representatives. b Back-to-back PCR validations of TaCMs LAHC01-1 and HJHN01-1. c Phylogenetic analysis of 7 TaCMs (highlighted in red) with their genetic neighbors retrieved from Genbank based on Rep proteins. d Mid-point rooted maximum likelihood phylogenetic tree of 13 TaCV Rep proteins (highlighted in red). Rep, Replication-associated protein; TaCVs, tick-associated circular DNA viruses; TACMs, tick-associated circular DNA molecules

Genetic diversity of circular DNA viruses

The remaining 13 circular genomes were 1669–2535 nt long and had the typical structure of circovirus, i.e. encoding Rep and catabolite activator proteins (Cap) (Fig. 2a). The two genes in most genomes were bi-directionally oriented, either detaching from each other or overlapping at their C-terminals (Fig. 2a). However, there were two genomes with uni-directionally oriented Rep and Cap (Fig. 2a), which were further verified by PCR amplification. These tick-associated circular DNA viruses (TaCVs) were distantly related to each other and showed low identities to known viruses (Fig. 2d). Phylogenetic analyses based on Rep proteins showed that these TaCVs were classified into the family Genomoviridae and unclassified viruses (Fig. 2d). Within the clade of Genomoviridae, the seven TaCVs genetically neighbored to genomoviruses (≥ 66.4% aa identities) originating from a wide range of samples, such as red panda feces (RpfaGmV), air-borne particulate matters (BfGmV 4 and PaGmV 14) and even citrus (CTGmV) (Fig. 2D). The highest similarity was observed for TaCV HLYD01-2/3, approximately 97% aa identical to CTGmVs. It should be noted that TaCV HLHC01-1 was 85.4% identical to TaGmV 2, a genomovirus recovered from an American Dermacentor variabilis tick. These TaCVs within the clades of unclassified viruses showed highly genetic diversities with each other and formed polyphyletic clades with CRESS DNA viruses from different hosts with 35.5–77.5% aa identities.


Tick-borne diseases pose a significant threat to public health and are mostly caused by infection with RNA viruses. Hence, RNA TBVs are often prioritized in investigations using either pathogen-specific or viromic methods, with the result that current knowledge of tick-borne RNA viruses is adequate. By contrast, studies on the diversity of DNA TBVs have been largely neglected. The tick-borne DNA virome in the present study revealed genetically diverse sscDNA viruses. The method used here is known to over-amplify small circular DNA genomes [25], but it also has an excellent capability to capture linear DNA genomes [26]. The presence of linear DNA bacteriophages highly implied that these ticks harbored rare linear eukaryotic DNA viruses. Viruses detected in ticks can be divided into tick-borne viruses that replicate in tick organs and viruses that are simply mechanically carried by ticks via feeding. The ticks collected from NE China were unfed and collected directly from the natural mixed forests. The detection of these viruses also in the tick larvae, suggests that the virome profiled here should be explained by tick-associated or environmental factors rather than host blood meals. The genetic diversity of a virus can be generally be partially shaped by its host species. For example, SFTSV is mainly maintained by H. longicornis [36], while CCHFV is associated with Hyalomma spp. [37]. However, the sscDNA viruses detected in the present study did not show any specificity to tick species, with some even having a close genetic relationship with viruses recovered from plants and animal feces. Hence, it is reasonable to infer that these tick-borne CRESS DNA viruses are associated primarily with environmental factors, rather than with viruses that actively replicate in ticks.

CRESS DNA viruses are among the most intriguing of viruses due to their ubiquity and diverse host range [38]. Their homologous Rep proteins facilitate the rolling-circle replication of genomes [39]. The number of known CRESS DNA viruses has greatly increased in recent years, largely due to the wide application of viral metagenomics. The latest International Committee on Taxonomy of Viruses (ICTV) taxonomy release groups CRESS DNA viruses into 11 families within the order Cressdnavircota, including the animal-associated Circoviridae, the plant-associated Geminiviridae, Metaxyviridae and Nanoviridae, the fungal-associated Genomoviridae, the diatom-associated Bacilladnaviridae, the protist-associated Naryaviridae, Nenyavirdae and Vilyaviridae, as well as two families with unknown host range: Smacoviridae and Redondoviridae [40, 41]. Interestingly, the Rep of CRESS DNA viruses is also associated with the sole protein encoded by a large number of circular, single-stranded DNA molecules. These molecules are 1.0–1.4 kb in length and have been assigned to the family Alphasatellitedae [42]. Of note is that all TaCMs identified in the present study were encompassed by CRESS DNA viruses, implying that TaCMs might originate from either defective genomes or small genome components of CRESS DNA viruses with multipartite genomes [43]. Most sequences recovered here were highly divergent from known references, with some capable of being classified as new species. These results suggest that tick-associated CRESS DNA viruses are much more diverse and abundant than expected.


This study provides important data to understand the spectrum of DNA viruses harbored by ticks in the NE provinces and Hubei Province of China. Although the ticks analyzed here were collected from a wide range of locations, covering five species with more than 1000 individuals, we did not find any linear DNA eukaryotic viruses, highly implying that the latter are only rarely harbored by ticks in these regions. However, considering that some DNA TBVs, such as ASFV and LSDV, pose severe threats to animal health, their surveillance should not be neglected but rather emphasized, similar to surveillance for RNA TBVs.

Availability of data and materials

All sequence reads generated in this study are available in the NCBI SRA database under BioProject accession PRJNA843368. The assembled complete genomes have been deposited in GenBank under accession numbers ON668986 to ON669005.



Amino acid


Circular Rep-encoding single-stranded


Eukaryotic Viral Reference Database


International Committee on Taxonomy of Viruses


Open reading frame


Single-stranded circular DNA


Replication-associated protein


Tick-associated circular DNA molecules


Tick-associated circular DNA viruses


Tick-borne disease


Tick-borne virus


  1. Sharifah N, Heo CC, Ehlers J, Houssaini J, Tappe D. Ticks and tick-borne pathogens in animals and humans in the island nations of Southeast Asia: a review. Acta Trop. 2020;209:105527.

    Article  PubMed  Google Scholar 

  2. Boulanger N, Boyer P, Talagrand-Reboul E, Hansmann Y. Ticks and tick-borne diseases. Med Mal Infect. 2019;49:87–97.

    Article  CAS  PubMed  Google Scholar 

  3. Sanchez-Vicente S, Tagliafierro T, Coleman JL, Benach JL, Tokarz R. Polymicrobial nature of tick-borne diseases. mBio. 2019;10:e02055-119.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Chomel B. Lyme disease. Rev Sci Tech. 2015;34:569–76.

    Article  CAS  PubMed  Google Scholar 

  5. 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:434–9.

    Article  CAS  PubMed  Google Scholar 

  6. Liu X, Zhang X, Wang Z, Dong Z, Xie S, Jiang M, et al. A tentative tamdy orthonairovirus related to febrile illness in Northwestern China. Clin Infect Dis. 2020;70:2155–60.

    Article  CAS  PubMed  Google Scholar 

  7. 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:1200–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. 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:2116–25.

    Article  CAS  PubMed  Google Scholar 

  9. Gritsun TS, Lashkevich VA, Gould EA. Tick-borne encephalitis. Antiviral Res. 2003;57:129–46.

    Article  CAS  PubMed  Google Scholar 

  10. Růžek D, Yakimenko VV, Karan LS, Tkachev SE. Omsk haemorrhagic fever.Lancet. 2010;376:2104–13.

    Article  PubMed  Google Scholar 

  11. Moming A, Yue X, Shen S, Chang C, Wang C, Luo T, et al. Prevalence and phylogenetic analysis of Crimean-Congo hemorrhagic fever virus in ticks from different ecosystems in Xinjiang. China Virol Sin. 2018;33:67–73.

    Article  CAS  PubMed  Google Scholar 

  12. 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:112–25.

    Article  CAS  PubMed  Google Scholar 

  13. 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:1975–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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:11480–92.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Burrage TG. African swine fever virus infection in Ornithodoros ticks. Virus Res. 2013;173:131–9.

    Article  CAS  PubMed  Google Scholar 

  16. Molini U, Aikukutu G, Khaiseb S, Haindongo NN, Lilungwe AC, Cattoli G, et al. Molecular characterization of lumpy skin disease virus in Namibia, 2017. Arch Virol. 2018;163:2525–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Lubinga JC, Tuppurainen ESM, Stoltsz WH, Ebersohn K, Coetzer JAW, Venter EH. Detection of lumpy skin disease virus in saliva of ticks fed on lumpy skin disease virus-infected cattle. Exp Appl Acarol. 2013;61:129–38.

    Article  CAS  PubMed  Google Scholar 

  18. Wang B, Sun L-D, Liu H-H, Wang Z-D, Zhao Y-K, Wang W, et al. Molecular detection of novel circoviruses in ticks in northeastern China. Ticks Tick Borne Dis. 2018;9:836–9.

    Article  PubMed  Google Scholar 

  19. Waits K, Edwards MJ, Cobb IN, Fontenele RS, Varsani A. Identification of an anellovirus and genomoviruses in ixodid ticks. Virus Genes. 2017;54:155–9.

    Article  PubMed  Google Scholar 

  20. Linden A, Gilliaux G, Paternostre J, Benzarti E, Rivas JF, Desmecht D, et al. A novel parvovirus, Roe deer copiparvovirus, identified in Ixodes ricinus ticks. Virus Genes. 2019;55:425–8.

    Article  CAS  PubMed  Google Scholar 

  21. Ravi A, Ereqat S, Al-Jawabreh A, Abdeen Z, Abu Shamma O, Hall H, et al. Metagenomic profiling of ticks: identification of novel rickettsial genomes and detection of tick-borne canine parvovirus. PLoS Negl Trop Dis. 2019;13:e0006805.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Kudelova M, Belvoncikova P, Vrbova M, Kovalova A, Stibraniova I, Kocakova P, et al. Detection of murine herpesvirus 68 (MHV-68) in Dermacentor reticulatus ticks. Microb Ecol. 2015;70:785–94.

    Article  PubMed  Google Scholar 

  23. Kudelova M, Janosova M, Belvoncikova P. First detection of murine herpesvirus 68 in adult Ixodes ricinus ticks. Folia Microbiol. 2018;63:511–5.

    Article  CAS  Google Scholar 

  24. Angly FE, Felts B, Breitbart M, Salamon P, Edwards RA, Carlson C, et al. The marine viromes of four oceanic regions. PLoS Biol. 2006;4:e368.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Parras-Molto M, Rodriguez-Galet A, Suarez-Rodriguez P, Lopez-Bueno A. Evaluation of bias induced by viral enrichment and random amplification protocols in metagenomic surveys of saliva DNA viruses. Microbiome. 2018;6:119.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Sun Y, Qu Y, Yan X, Yan G, Chen J, Wang G, et al. Comprehensive evaluation of RNA and DNA viromic methods based on species richness and abundance analyses using marmot rectal samples. mSystems. 2022;7:e0043022.

    Article  PubMed  Google Scholar 

  27. Zhang Y, Si B-Y, Liu B-H, Chang G-H, Yang Y-H, Huo Q-B, et al. Complete genomic characterization of two tick-borne encephalitis viruses isolated from China. Virus Res. 2012;167:310–3.

    Article  CAS  PubMed  Google Scholar 

  28. Zhang Y, Shen S, Shi J, Su Z, Li M, Zhang W, et al. Isolation, characterization, and phylogenic analysis of three new severe fever with thrombocytopenia syndrome bunyavirus strains derived from Hubei Province. China Virol Sin. 2017;32:89–96.

    Article  PubMed  Google Scholar 

  29. Gao X, Nasci R, Liang G. The neglected arboviral infections in mainland China. PLoS Negl Trop Dis. 2010;4:e624.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Deng G. Economic insects of China: class Arachnida, subclass Acari. 1st ed. Beijing: Science Press; 1978.

    Google Scholar 

  31. Yamaguti N, Tipton VJ, Keegan HL, et al. Ticks of Japan, Korea, and the Ryukyu Islands. Brigh Young Univ Sci Bull Biol. 1971;15:1–226.

    Google Scholar 

  32. Chen Z, Yang X, Bu F, Yang X, Yang X, Liu J. Ticks (acari: ixodoidea: argasidae, ixodidae) of China. Exp Appl Acarol. 2010;51:393–404.

    Article  PubMed  Google Scholar 

  33. Guglielmone AA, Robbins RG, Apanaskevich DA, Petney TN, EstradaPeña A, Horak IG. The hard ticks of the world: (Acari: Ixodida: Ixodidae). 1st ed. Dordrecht: Springer; 2014.

    Book  Google Scholar 

  34. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech. 1994;3:294–9.

    CAS  Google Scholar 

  35. Chen J, Sun Y, Yan X, Ren Z, Wang G, Liu Y, et al. Elimination of foreign sequences in eukaryotic viral reference genomes improves the accuracy of virome analysis. mSystems. 2022;7:e0090722.

    Article  PubMed  Google Scholar 

  36. Liu Q, He B, Huang S-Y, Wei F, Zhu X-Q. Severe fever with thrombocytopenia syndrome, an emerging tick-borne zoonosis. Lancet Infect Dis. 2014;14:763–72.

    Article  CAS  PubMed  Google Scholar 

  37. 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:141–7.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Zhao L, Rosario K, Breitbart M, Duffy S. Eukaryotic circular rep-encoding single-stranded DNA (CRESS DNA) viruses: ubiquitous viruses with small genomes and a diverse host range. Adv Virus Res. 2019;103:71–133.

    Article  CAS  PubMed  Google Scholar 

  39. Rosario K, Duffy S, Breitbart M. A field guide to eukaryotic circular single-stranded DNA viruses: insights gained from metagenomics. Arch Virol. 2012;157:1851–71.

    Article  CAS  PubMed  Google Scholar 

  40. Kinsella CM, Bart A, Deijs M, Broekhuizen P, Kaczorowska J, Jebbink MF, et al. Entamoeba and Giardia parasites implicated as hosts of CRESS viruses. Nat Commun. 2020;11:4620.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Hulo C, de Castro E, Masson P, Bougueleret L, Bairoch A, Xenarios I, et al. ViralZone: a knowledge resource to understand virus diversity. Nucleic Acids Res. 2011;39:D576–82.

    Article  CAS  PubMed  Google Scholar 

  42. Briddon RW, Martin DP, Roumagnac P, Navas-Castillo J, Fiallo-Olive E, Moriones E, et al. Alphasatellitidae: a new family with two subfamilies for the classification of geminivirus- and nanovirus-associated alphasatellites. Arch Virol. 2018;163:2587–600.

    Article  CAS  PubMed  Google Scholar 

  43. Kraberger S, Arguello-Astorga GR, Greenfield LG, Galilee C, Law D, Martin DP, et al. Characterisation of a diverse range of circular replication-associated protein encoding DNA viruses recovered from a sewage treatment oxidation pond. Infect Genet Evol. 2015;31:73–86.

    Article  CAS  PubMed  Google Scholar 

Download references


Not applicable.


This study was supported by the National Natural Science Foundation of China (32022083). The funder had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Author information

Authors and Affiliations



YL, LG and GW performed the experiments, completed the data analysis and drafted the manuscript. YL and GW performed morphological identifications. LG, FG, ZT, DX and LS participated in sample collection. LY, GZ and CT reviewed the manuscript. BH contributed to the conception and design of this work, supervised the research, analyzed the data and revised the manuscript. All authors have read and approved the final manuscript.

Corresponding author

Correspondence to Biao He.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

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

Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Guo, L., Wang, G. et al. DNA virome of ticks in the Northeast and Hubei provinces of China reveals diverse single-stranded circular DNA viruses. Parasites Vectors 16, 61 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Ticks
  • DNA virome
  • CRESS DNA viruses
  • Genetic diversity