Geographic distribution and modeling of ticks in the Republic of Korea and the application of tick models towards understanding the distribution of associated pathogenic agents

Tick Diversity and Pathogen Transmission in Daejeon, Korea: Implications from Companion Animals and Walking Trails

With the ongoing global warming-induced climate change, there has been a surge in vector-borne diseases, particularly tick-borne diseases (TBDs). As the population of companion animals grows, there is growing concern from a One Health perspective about the potential for these animals to spread TBDs. In this study, ticks were collected from companion animals and the surrounding environment in Daejeon Metropolitan City, Korea, using flagging and dragging, and CO2 trap methods. These ticks were then subjected to conventional (nested) PCR for severe fever with thrombocytopenia syndrome virus (SFTSV), Anaplasma spp., Ehrlichia spp., and Borrelia spp. We identified a total of 29,176 ticks, consisting of three genera and four species: H. longicornis, H. flava, I. nipponensis, and A. testudinarium. Notably, H. longicornis was the predominant species. The presence of A. testudinarium suggested that the species traditionally found in southern regions are migrating northward, likely as a result of climate change. Our PCR results confirmed the presence of all four pathogens in both companion animals and the surrounding environment, underscoring the potential for the indirect transmission of tick-borne pathogens to humans through companion animals. These findings emphasize the importance of the ongoing surveillance of companion animals in the management and control of TBDs.

Transcriptome analysis of Haemaphysalis flava female using Illumina HiSeq 4000 sequencing: de novo assembly, functional annotation and discovery of SSR markers

BACKGROUND

Ticks are ectoparasites capable of directly damaging their hosts and transmitting vector-borne diseases. The ixodid tick Haemaphysalis flava has a broad distribution that extends from East to South Asia. This tick is a reservoir of severe fever with thrombocytopenia syndrome virus (SFTSV) that causes severe hemorrhagic disease, with cases reported from China, Japan and South Korea. Recently, the distribution of H. flava in South Korea was found to overlap with the occurrence of SFTSV.

METHODS

This study was undertaken to discover the molecular resources of H. flava female ticks using the Illumina HiSeq 4000 system, the Trinity de novo sequence assembler and annotation against public databases. The locally curated Protostome database (PANM-DB) was used to screen the putative adaptation-related transcripts classified to gene families, such as angiotensin-converting enzyme, aquaporin, adenylate cyclase, AMP-activated protein kinase, glutamate receptors, heat shock proteins, molecular chaperones, insulin receptor, mitogen-activated protein kinase and solute carrier family proteins. Also, the repeats and simple sequence repeats (SSRs) were screened from the unigenes using RepeatMasker (v4.0.6) and MISA (v1.0) software tools, followed by the designing of SSRs flanking primers using BatchPrimer 3 (v1.0) software.

RESULTS

The transcriptome produced a total of 69,822 unigenes, of which 46,175 annotated to the homologous proteins in the PANM-DB. The unigenes were also mapped to the EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) specializations. Promiscuous presence of protein kinase, zinc finger (C2H2-type), reverse transcriptase, and RNA recognition motif domains was observed in the unigenes. A total of 3480 SSRs were screened, of which 1907 and 1274 were found as tri- and dinucleotide repeats, respectively. A list of primer sequences flanking the SSR motifs was detailed for validation of polymorphism in H. flava and the related tick species.

CONCLUSIONS

The reference transcriptome information on H. flava female ticks will be useful for an enriched understanding of tick biology, its competency to act as a vector and the study of species diversity related to disease transmission.

Geographical distribution of Ixodes persulcatus and associated pathogens: Analysis of integrated data from a China field survey and global published data

The increasing incidence and range expansion of tick-borne diseases have caused global threats to human and animal health under the background of climate and socioeconomic changes. As an efficient vector in transmission of tick-borne diseases, a growing burden caused by Ixodes persulcatus and associated pathogens could not be underestimated. This study summarized the distribution, hosts, and pathogens of I. persulcatus, and predicted the suitable habitats of this tick species worldwide. An integrated database involving a field survey, reference book, literature review, and related website was constructed. Location records of I. persulcatus and associated pathogens were incorporated into distribution maps using ArcGIS software. Positive rates for I. persulcatus-associated agents were estimated by meta-analysis. The global distribution of the tick species was predicted using Maxent model. I. persulcatus was distributed in 14 countries across the Eurasian continent, involving Russia, China, Japan, and several Baltic Sea states, which ranged between 21°N to 66°N. The tick species fed on 46 species of hosts, and 51 tick-borne agents could be harbored by I. persulcatus. The predictive model showed that I. persulcatus could be predominantly distributed in northern Europe, western Russia, and northern China. Our study fully clarified the potential public health risks posed by I. persulcatus and I. persulcatus-borne pathogens. Surveillance and control measures of tick-borne diseases should be enhanced to promote the health of humans, animals, and ecosystems.

Reconstruction of mitochondrial genomes from raw sequencing data provides insights on the phylogeny of Ixodes ticks and cautions for species misidentification

High-throughput sequencing (HTS) technology has profoundly been involved in sequencing whole genomes of several organisms in a fast and cost-effective manner. Although HTS provides an alternative biomonitoring method to the time-consuming and taxonomy-expertise dependent morphological approach, still we cannot rule out the possibility of the impediment and misidentification biases. In this article we aim to retrieve whole mitochondrial genome (mitogenome) sequences from publicly available raw sequencing data for phylogenetic comparison of Ixodes persulcatus. For this comparison, we sequenced whole mitogenomes of four I. persulcatus ticks from Japan and constructed mitogenomes from raw sequencing data of 74 I. persulcatus ticks from China. Bayesian phylogenetic trees were inferred by the concatenated fifteen mitochondrial genes. We further tested our results by the phylogenetic analysis of cytochrome c oxidase subunit 1 (cox1) gene and internal transcribed spacer 2 (ITS2) sequences. Our findings showed that 70 constructed mitogenomes from China were clustered with the sequenced four mitogenomes of I. persulcatus from Japan. We also revealed that mitogenome sequences retrieved from two data sets CRR142297 and CRR142298 were clustered with Ixodes nipponensis. Moreover, other two mitogenome sequences from CRR142310 and CRR142311 formed a clade with Ixodes pavlovskyi. The phylogenetic analysis of cox1 gene and ITS2 sequences confirmed the identification errors of these four samples. The overall phylogenetics in our study concluded that accurate morphological identification is necessary before implementing HTS to avoid any misidentification biases.