Mapping the distribution of tick-borne encephalitis in mainland China

Behavioral and cognitive factors influencing tick-borne disease risk in northeast China: Implications for prevention and control strategies

The growth in ecotourism and nature-based recreational activities in China has resulted in an increased frequency of visits to green spaces, thereby elevating exposure to ticks and the subsequent risk of tick-borne diseases. This study comprehensively investigate individual behavioral and cognitive factors associated with the risk of contracting tick-borne diseases to facilitate the development of effective prevention and control strategies, supporting public health initiatives in high-prevalence regions. We conducted an extensive questionnaire survey among 3000 residents from three northeastern provinces in China (Heilongjiang, Jilin, and Liaoning), where tick-borne diseases exhibit relatively high prevalence. The survey focused on gathering information regarding participants' tick bite history, perception of tick-borne disease risks, and outdoor activity patterns. Using structural equations analysis, we explored the pathways and strengths of the associations between these factors. Our findings revealed an average self-reported tick bite rate of 14% among the participants. Notably, tick-borne encephalitis exhibited the highest self-reported prevalence of infection (4%) among tick-borne diseases, while both Lyme disease and Severe fever with thrombocytopenia syndrome had a prevalence of 2%. The average rate of tick bites among respondents' pets was 14%, with bites predominantly located on the ears, back, and abdomen. The strongest correlation was observed between tick bite rate and subsequent infections, emphasizing its role as the primary contributing factors to infectious status. Moreover, our results indicated that the causal structure of tick-borne disease infections varied across different cities, underscoring the significance of considering the ecological environment and regional knowledge on ticks. This study provides valuable insights into the current landscape of tick-borne disease infections in northeast China and identifies potential behavioral and cognitive factors, an aspect that has not been previously investigated. Our findings enable predictions on the future impact of knowledge dissemination efforts and improved urban facilities on mitigating tick bites and reducing tick-borne disease infections.

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.

Risk factors for tick attachment in companion animals in Great Britain: a spatiotemporal analysis covering 2014-2021

BACKGROUND

Ticks are an important driver of veterinary health care, causing irritation and sometimes infection to their hosts. We explored epidemiological and geo-referenced data from > 7 million electronic health records (EHRs) from cats and dogs collected by the Small Animal Veterinary Surveillance Network (SAVSNET) in Great Britain (GB) between 2014 and 2021 to assess the factors affecting tick attachment in an individual and at a spatiotemporal level.

METHODS

EHRs in which ticks were mentioned were identified by text mining; domain experts confirmed those with ticks on the animal. Tick presence/absence records were overlaid with a spatiotemporal series of climate, environment, anthropogenic and host distribution factors to produce a spatiotemporal regression matrix. An ensemble machine learning spatiotemporal model was used to fine-tune hyperparameters for Random Forest, Gradient-boosted Trees and Generalized Linear Model regression algorithms, which were then used to produce a final ensemble meta-learner to predict the probability of tick attachment across GB at a monthly interval and averaged long-term through 2014-2021 at a spatial resolution of 1 km. Individual host factors associated with tick attachment were also assessed by conditional logistic regression on a matched case-control dataset.

RESULTS

In total, 11,741 consultations were identified in which a tick was recorded. The frequency of tick records was low (0.16% EHRs), suggesting an underestimation of risk. That said, increased odds for tick attachment in cats and dogs were associated with younger adult ages, longer coat length, crossbreeds and unclassified breeds. In cats, males and entire animals had significantly increased odds of recorded tick attachment. The key variables controlling the spatiotemporal risk for tick attachment were climatic (precipitation and temperature) and vegetation type (Enhanced Vegetation Index). Suitable areas for tick attachment were predicted across GB, especially in forests and grassland areas, mainly during summer, particularly in June.

CONCLUSIONS

Our results can inform targeted health messages to owners and veterinary practitioners, identifying those animals, seasons and areas of higher risk for tick attachment and allowing for more tailored prophylaxis to reduce tick burden, inappropriate parasiticide treatment and potentially TBDs in companion animals and humans. Sentinel networks like SAVSNET represent a novel complementary data source to improve our understanding of tick attachment risk for companion animals and as a proxy of risk to humans.

Metatranscriptomics Reveals the RNA Virome of Ixodes Persulcatus in the China-North Korea Border, 2017

In recent years, numerous viruses have been identified from ticks, and some have been linked to clinical cases of emerging tick-borne diseases. Chinese northeast frontier is tick infested. However, there is a notable lack of systematic monitoring efforts to assess the viral composition in the area, leaving the ecological landscape of viruses carried by ticks not clear enough. Between April and June 2017, 7101 ticks were collected to perform virus surveillance on the China-North Korea border, specifically in Tonghua, Baishan, and Yanbian. A total of 2127 Ixodes persulcatus were identified. Further investigation revealed the diversity of tick-borne viruses by transcriptome sequencing of Ixodes persulcatus. All ticks tested negative for tick-borne encephalitis virus. Transcriptome sequencing expanded 121 genomic sequence data of 12 different virus species from Ixodes persulcatus. Notably, a new segmented flavivirus, named Baishan Forest Tick Virus, were identified, closely related to Alongshan virus and Harz mountain virus. Therefore, this new virus may pose a potential threat to humans. Furthermore, the study revealed the existence of seven emerging tick-borne viruses dating back to 2017. These previously identified viruses included Mudanjiang phlebovirus, Onega tick phlebovirus, Sara tick phlebovirus, Yichun mivirus, and three unnamed viruses (one belonging to the Peribunyaviridae family and the other two belonging to the Phenuiviridae family). The existence of these emerging tick-borne viruses in tick samples collected in 2017 suggests that their history may extend further than previously recognized. This study provides invaluable insights into the virome of Ixodes persulcatus in the China-North Korea border region, enhancing our ongoing efforts to manage the risks associated with tick-borne viruses.

Co-infection of tick-borne bacterial pathogens in ticks in Inner Mongolia, China

Tick-borne infectious diseases pose a serious health threat in certain regions of the world. Emerging infectious diseases caused by novel tick-borne pathogens have been reported that are causing particular concern. Several tick-borne diseases often coexist in the same foci, and a single vector tick can transmit two or more pathogens at the same time, which greatly increases the probability of co-infection in host animals and humans and can lead to an epidemic of tick-borne disease. The lack of epidemiological data and information on the specific clinical symptoms related to co-infection with tick-borne pathogens means that it is not currently possible to accurately and rapidly distinguish between a single pathogen infection and co-infection with multiple pathogens, which can have serious consequences. Inner Mongolia in the north of China is endemic for tick-borne infectious diseases, especially in the eastern forest region. Previous studies have found that more than 10% of co-infections were in host-seeking ticks. However, the lack of data on the specific types of co-infection with pathogens makes clinical treatment difficult. In our study, we present data on the co-infection types and the differences in co-infection among different ecological regions through genetic analysis of tick samples collected throughout Inner Mongolia. Our findings may aid clinicians in the diagnosis of concomitant tick-borne infectious diseases.

Mapping the risk distribution of Borrelia burgdorferi sensu lato in China from 1986 to 2020: a geospatial modelling analysis

Lyme borreliosis, recognized as one of the most important tick-borne diseases worldwide, has been increasing in incidence and spatial extent. Currently, there are few geographic studies about the distribution of Lyme borreliosis risk across China. Here we established a nationwide database that involved Borrelia burgdorferi sensu lato (B. burgdorferi) detected in humans, vectors, and animals in China. The eco-environmental factors that shaped the spatial pattern of B. burgdorferi were identified by using a two-stage boosted regression tree model and the model-predicted risks were mapped. During 1986−2020, a total of 2,584 human confirmed cases were reported in 25 provinces. Borrelia burgdorferi was detected from 35 tick species with the highest positive rates in Ixodes granulatus, Hyalomma asiaticum, Ixodes persulcatus, and Haemaphysalis concinna ranging 20.1%−24.0%. Thirteen factors including woodland, NDVI, rainfed cropland, and livestock density were determined as important drivers for the probability of B. burgdorferi occurrence based on the stage 1 model. The stage 2 model identified ten factors including temperature seasonality, NDVI, and grasslands that were the main determinants used to distinguish areas at high or low-medium risk of B. burgdorferi, interpreted as potential occurrence areas within the area projected by the stage 1 model. The projected high-risk areas were not only concentrated in high latitude areas, but also were distributed in middle and low latitude areas. These high-resolution evidence-based risk maps of B. burgdorferi was first created in China and can help as a guide to future surveillance and control and help inform disease burden and infection risk estimates.

Analysis of the spatial-temporal components driving transmission of the severe fever with thrombocytopenia syndrome in Shandong Province, China, 2016-2018

Existing models about the spatial-temporal distribution of the severe fever with thrombocytopenia syndrome (SFTS) entirely concentrate on aggregation, which provides limited knowledge to develop effective measures to control the epidemic of SFTS. This study aimed to identify the main spatial-temporal components and heterogeneity in different regions in Shandong Province, China. We applied the spatial-temporal multicomponent model to detect the spatial-temporal component values. A total of 2814 cases were reported from 2016 to 2018 in Shandong Province. The prevalence rate was 0.627 per 100,000, with an overall case fatality rate of 8.99%. SFTS cases were mostly clustered in central and eastern regions of Shandong Province. The total effect values of the autoregressive component, the spatiotemporal component and the endemic component were 0.586, 0.244 and 0.084, respectively, which demonstrated that the autoregressive component was the main factor driving the incidence of SFTS, followed by the spatiotemporal component. Gross domestic product per capita and weekly mean atmospheric pressure contributed to the incidence of SFTS with inverse effects. Obvious heterogeneity across regions for the autoregressive component and the spatiotemporal component was identified. In conclusion, the autoregressive and spatiotemporal components play a key role in driving the transmission of SFTS in Shandong Province. Based on the main component values, targeted measures should be formulated to control SFTS epidemics in different regions.

Mapping the viruses belonging to the order Bunyavirales in China

BACKGROUND

Viral pathogens belonging to the order Bunyavirales pose a continuous background threat to global health, but the fact remains that they are usually neglected and their distribution is still ambiguously known. We aim to map the geographical distribution of Bunyavirales viruses and assess the environmental suitability and transmission risk of major Bunyavirales viruses in China.

METHODS

We assembled data on all Bunyavirales viruses detected in humans, animals and vectors from multiple sources, to update distribution maps of them across China. In addition, we predicted environmental suitability at the 10 km × 10 km pixel level by applying boosted regression tree models for two important Bunyavirales viruses, including Crimean-Congo hemorrhagic fever virus (CCHFV) and Rift Valley fever virus (RVFV). Based on model-projected risks and air travel volume, the imported risk of RVFV was also estimated from its endemic areas to the cities in China.

RESULTS

Here we mapped all 89 species of Bunyavirales viruses in China from January 1951 to June 2021. Nineteen viruses were shown to infect humans, including ten species first reported as human infections. A total of 447,848 cases infected with Bunyavirales viruses were reported, and hantaviruses, Dabie bandavirus and Crimean-Congo hemorrhagic fever virus (CCHFV) had the severest disease burden. Model-predicted maps showed that Xinjiang and southwestern Yunnan had the highest environmental suitability for CCHFV occurrence, mainly related to Hyalomma asiaticum presence, while southern China had the highest environmental suitability for Rift Valley fever virus (RVFV) transmission all year round, mainly driven by livestock density, mean precipitation in the previous month. We further identified three cities including Guangzhou, Beijing and Shanghai, with the highest imported risk of RVFV potentially from Egypt, South Africa, Saudi Arabia and Kenya.

CONCLUSIONS

A variety of Bunyavirales viruses are widely distributed in China, and the two major neglected Bunyavirales viruses including CCHFV and RVFV, both have the potential for outbreaks in local areas of China. Our study can help to promote the understanding of risk distribution and disease burden of Bunyavirales viruses in China, and the risk maps of CCHFV and RVFV occurrence are crucial to the targeted surveillance and control, especially in seasons and locations at high risk.

A Scoping Review of Species Distribution Modeling Methods for Tick Vectors

Background

Globally, tick-borne disease is a pervasive and worsening problem that impacts human and domestic animal health, livelihoods, and numerous economies. Species distribution models are useful tools to help address these issues, but many different modeling approaches and environmental data sources exist.

Objective

We conducted a scoping review that examined all available research employing species distribution models to predict occurrence and map tick species to understand the diversity of model strategies, environmental predictors, tick data sources, frequency of climate projects of tick ranges, and types of model validation methods.

Design

Following the PRISMA-ScR checklist, we searched scientific databases for eligible articles, their references, and explored related publications through a graphical tool (www.connectedpapers.com). Two independent reviewers performed article selection and characterization using a priori criteria.

Results

We describe data collected from 107 peer-reviewed articles that met our inclusion criteria. The literature reflects that tick species distributions have been modeled predominantly in North America and Europe and have mostly modeled the habitat suitability for Ixodes ricinus (n = 23; 21.5%). A wide range of bioclimatic databases and other environmental correlates were utilized among models, but the WorldClim database and its bioclimatic variables 1–19 appeared in 60 (56%) papers. The most frequently chosen modeling approach was MaxEnt, which also appeared in 60 (56%) of papers. Despite the importance of ensemble modeling to reduce bias, only 23 papers (21.5%) employed more than one algorithm, and just six (5.6%) used an ensemble approach that incorporated at least five different modeling methods for comparison. Area under the curve/receiver operating characteristic was the most frequently reported model validation method, utilized in nearly all (98.9%) included studies. Only 21% of papers used future climate scenarios to predict tick range expansion or contraction. Regardless of the representative concentration pathway, six of seven genera were expected to both expand and retract depending on location, while Ornithodoros was predicted to only expand beyond its current range.

Conclusion

Species distribution modeling techniques are useful and widely employed tools for predicting tick habitat suitability and range movement. However, the vast array of methods, data sources, and validation strategies within the SDM literature support the need for standardized protocols for species distribution and ecological niche modeling for tick vectors.

Geo-Spatial Characteristics of 567 Places of Tick-Borne Encephalitis Infection in Southern Germany, 2018-2020

Tick-borne encephalitis (TBE) is a growing public health problem with increasing incidence and expanding risk areas. Improved prevention requires better understanding of the spatial distribution and ecological determinants of TBE transmission. However, a TBE risk map at sub-district level is still missing for Germany. We investigated the distribution and geo-spatial characteristics of 567 self-reported places of probable TBE infection (POI) from 359 cases notified in 2018-2020 in the study area of Bavaria and Baden-Wuerttemberg, compared to 41 confirmed TBE foci and 1701 random comparator places. We built an ecological niche model to interpolate TBE risk to the entire study area. POI were distributed heterogeneously at sub-district level, as predicted probabilities varied markedly across regions (range 0-93%). POI were spatially associated with abiotic, biotic, and anthropogenic geo-spatial characteristics, including summer precipitation, population density, and annual frost days. The model performed with 69% sensitivity and 63% specificity at an optimised probability threshold (0.28) and an area under the curve of 0.73. We observed high predictive probabilities in small-scale areas, consistent with the known circulation of the TBE virus in spatially restricted microfoci. Supported by further field work, our findings may help identify new TBE foci. Our fine-grained risk map could supplement targeted prevention in risk areas.

Near-term forecasting of companion animal tick paralysis incidence: An iterative ensemble model

Tick paralysis resulting from bites from Ixodes holocyclus and I. cornuatus is one of the leading causes of emergency veterinary admissions for companion animals in Australia, often resulting in death if left untreated. Availability of timely information on periods of increased risk can help modulate behaviors that reduce exposures to ticks and improve awareness of owners for the need of lifesaving preventative ectoparasite treatment. Improved awareness of clinicians and pet owners about temporal changes in tick paralysis risk can be assisted by ecological forecasting frameworks that integrate environmental information into statistical time series models. Using an 11-year time series of tick paralysis cases from veterinary clinics in one of Australia's hotspots for the paralysis tick Ixodes holocyclus, we asked whether an ensemble model could accurately forecast clinical caseloads over near-term horizons. We fit a series of statistical time series (ARIMA, GARCH) and generative models (Prophet, Generalised Additive Model) using environmental variables as predictors, and then combined forecasts into a weighted ensemble to minimise prediction interval error. Our results indicate that variables related to temperature anomalies, levels of vegetation moisture and the Southern Oscillation Index can be useful for predicting tick paralysis admissions. Our model forecasted tick paralysis cases with exceptional accuracy while preserving epidemiological interpretability, outperforming a field-leading benchmark Exponential Smoothing model by reducing both point and prediction interval errors. Using online particle filtering to assimilate new observations and adjust forecast distributions when new data became available, our model adapted to changing temporal conditions and provided further reduced forecast errors. We expect our model pipeline to act as a platform for developing early warning systems that can notify clinicians and pet owners about heightened risks of environmentally driven veterinary conditions.

Spatial and temporal distribution patterns of tick-borne diseases (Tick-borne Encephalitis and Lyme Borreliosis) in Germany

BACKGROUND

In the face of ongoing climate warming, vector-borne diseases are expected to increase in Europe, including tick-borne diseases (TBD). The most abundant tick-borne diseases in Germany are Tick-Borne Encephalitis (TBE) and Lyme Borreliosis (LB), with Ixodes ricinus as the main vector.

METHODS

In this study, we display and compare the spatial and temporal patterns of reported cases of human TBE and LB in relation to some associated factors. The comparison may help with the interpretation of observed spatial and temporal patterns.

RESULTS

The spatial patterns of reported TBE cases show a clear and consistent pattern over the years, with many cases in the south and only few and isolated cases in the north of Germany. The identification of spatial patterns of LB disease cases is more difficult due to the different reporting practices in the individual federal states. Temporal patterns strongly fluctuate between years, and are relatively synchronized between both diseases, suggesting common driving factors. Based on our results we found no evidence that weather conditions affect the prevalence of both diseases. Both diseases show a gender bias with LB bing more commonly diagnosed in females, contrary to TBE being more commonly diagnosed in males.

CONCLUSION

For a further investigation of of the underlying driving factors and their interrelations, longer time series as well as standardised reporting and surveillance system would be required.

Wild animal and zoonotic disease risk management and regulation in China: Examining gaps and One Health opportunities in scope, mandates, and monitoring systems

Emerging diseases of zoonotic origin such as COVID-19 are a continuing public health threat in China that lead to a significant socioeconomic burden. This study reviewed the current laws and regulations, government reports and policy documents, and existing literature on zoonotic disease preparedness and prevention across the forestry, agriculture, and public health authorities in China, to articulate the current landscape of potential risks, existing mandates, and gaps. A total of 55 known zoonotic diseases (59 pathogens) are routinely monitored under a multi-sectoral system among humans and domestic and wild animals in China. These diseases have been detected in wild mammals, birds, reptiles, amphibians, and fish or other aquatic animals, the majority of which are transmitted between humans and animals via direct or indirect contact and vectors. However, this current monitoring system covers a limited scope of disease threats and animal host species, warranting expanded review for sources of disease and pathogen with zoonotic potential. In addition, the governance of wild animal protection and utilization and limited knowledge about wild animal trade value chains present challenges for zoonotic disease risk assessment and monitoring, and affect the completeness of mandates and enforcement. A coordinated and collaborative mechanism among different departments is required for the effective monitoring and management of disease emergence and transmission risks in the animal value chains. Moreover, pathogen surveillance among wild animal hosts and human populations outside of the routine monitoring system will fill the data gaps and improve our understanding of future emerging zoonotic threats to achieve disease prevention. The findings and recommendations will advance One Health collaboration across government and non-government stakeholders to optimize monitoring and surveillance, risk management, and emergency responses to known and novel zoonotic threats, and support COVID-19 recovery efforts.

Molecular detection and phylogenetic analysis of tick-borne encephalitis virus in ticks in northeastern China

Tick-borne encephalitis virus (TBEV) is an important causative agent that causes neurological infections in humans and animals. In recent years, only few epidemiological surveys on TBEV have been conducted in China. The purpose of this study was to determine the prevalence and subtype of TBEV in ticks in northeastern (NE) China. A total of 3799 questing ticks were collected in NE China between April 2015 and June 2016. Ticks were pooled and tested for TBEV RNA using semi-nested reverse transcriptase-polymerase chain reaction. Positive pools were used to isolate the virus and amplify complete sequences, followed by sequence identity and phylogenetic analysis. TBEV RNA was detected in Ixodes persulcatus ticks at a total prevalence of 2.9% (6/143; 95% confidence interval: 1.2%-5.9%). Three TBEV strains were isolated (JL-T75, HLB-T74, and DXAL-T83) and showed 93.9%-99.1% nucleotide identities and 97.1%-99.5% amino acid identities in Far Eastern (FE) TBEV subtypes, and 82.9%-87.6% nucleotide identities and 92.9%-96.4% amino acid identities in other subtypes. For polyprotein, the JL-T75, HLB-T74, and DXAL-T83 strains showed 29, 50, and 55 amino acid residues, respectively, different from those in the TBEV vaccine (Senzhang) strain in China. Phylogenetic analysis revealed that these viruses were clustered in the FE-TBEV branch but formed distinct clades depending on the natural foci. The results of this study suggest that the FE-TBEV subtype is still endemic in I. persulcatus ticks in NE China, and the viruses in different natural foci in NE China are more likely to have genetic differences.

Recovery of a Far-Eastern Strain of Tick-Borne Encephalitis Virus with a Full-Length Infectious cDNA Clone

Tick-borne encephalitis virus (TBEV) is a pathogenic virus known to cause central nervous system (CNS) diseases in humans, and has become an increasing public health threat nowadays. The rates of TBEV infection in the endemic countries are increasing. However, there is no effective antiviral against the disease. This underscores the urgent need for tools to study the emergence and pathogenesis of TBEV and to accelerate the development of vaccines and antivirals. In this study, we reported an infectious cDNA clone of TBEV that was isolated in China (the WH2012 strain). A beta-globin intron was inserted in the coding region of nonstructural protein 1 (NS1) gene to improve the stability of viral genome in bacteria. In mammalian cells, the inserted intron was excised and spliced precisely, which did not lead to the generation of inserted mutants. High titers of infectious progeny viruses were generated after the transfection of the infectious clone. The cDNA-derived TBEV replicated efficiently, and caused typical cytopathic effect (CPE) and plaques in BHK-21 cells. In addition, the CPE and growth curve of cDNA-derived virus were similar to that of its parental isolate in cells. Together, we have constructed the first infectious TBEV cDNA clone in China, and the clone can be used to investigate the genetic determinants of TBEV virulence and disease pathogenesis, and to develop countermeasures against the virus.

Mapping ticks and tick-borne pathogens in China

Understanding ecological niches of major tick species and prevalent tick-borne pathogens is crucial for efficient surveillance and control of tick-borne diseases. Here we provide an up-to-date review on the spatial distributions of ticks and tick-borne pathogens in China. We map at the county level 124 tick species, 103 tick-borne agents, and human cases infected with 29 species (subspecies) of tick-borne pathogens that were reported in China during 1950-2018. Haemaphysalis longicornis is found to harbor the highest variety of tick-borne agents, followed by Ixodes persulcatus, Dermacentor nutalli and Rhipicephalus microplus. Using a machine learning algorithm, we assess ecoclimatic and socioenvironmental drivers for the distributions of 19 predominant vector ticks and two tick-borne pathogens associated with the highest disease burden. The model-predicted suitable habitats for the 19 tick species are 14‒476% larger in size than the geographic areas where these species were detected, indicating severe under-detection. Tick species harboring pathogens of imminent threats to public health should be prioritized for more active field surveillance.

Vectors of disease at the northern distribution limit of the genus Dermacentor in Eurasia: D. reticulatus and D. silvarum

The two ixodid tick species Dermacentor reticulatus (Fabricius) and Dermacentor silvarum Olenev occur at the northern distribution limit of the genus Dermacentor in Eurasia, within the belt of [Formula: see text] latitude. Whilst the distribution area of D. reticulatus extends from the Atlantic coast of Portugal to Western Siberia, that of D. silvarum extends from Western Siberia to the Pacific coast. In Western Siberia, the distribution areas of the two Dermacentor species overlap. Although the two tick species are important vectors of disease, detailed information concerning the entire distribution area, climate adaptation, and proven vector competence is still missing. A dataset was compiled, resulting in 2188 georeferenced D. reticulatus and 522 D. silvarum locations. Up-to-date maps depicting the geographical distribution and climate adaptation of the two Dermacentor species are presented. To investigate the climate adaptation of the two tick species, the georeferenced locations were superimposed on a high-resolution map of the Köppen-Geiger climate classification. The frequency distribution of D. reticulatus under different climates shows two major peaks related to the following climates: warm temperate with precipitation all year round (57%) and boreal with precipitation all year round (40%). The frequency distribution of D. silvarum shows also two major peaks related to boreal climates with precipitation all year round (30%) and boreal winter dry climates (60%). Dermacentor silvarum seems to be rather flexible concerning summer temperatures, which can range from cool to hot. In climates with cool summers D. reticulatus does not occur, it prefers warm and to a lesser extent hot summers. Lists are given in this paper for cases of proven vector competence for various agents of both Dermacentor species. For the first time, the entire distribution areas of D. reticulatus and D. silvarum were mapped using georeferenced data. Their climate adaptations were quantified by Köppen profiles.

Modelling habitat suitability for occurrence of human tick-borne encephalitis (TBE) cases in Finland

The numbers of reported human tick-borne encephalitis (TBE) cases in Europe have increased in several endemic regions (including Finland) in recent decades, indicative of an increasing threat to public health. As such, it is important to identify the regions at risk and the most influential factors associated with TBE distributions, particularly in understudied regions. This study aimed to identify the risk areas of TBE transmission in two different datasets based on human TBE disease cases from 2007 to 2011 (n = 86) and 2012-2017 (n  = 244). We also examined which factors best explain the presence of human TBE cases. We used ensemble modelling to determine the relationship of TBE occurrence with environmental, ecological, and anthropogenic factors in Finland. Geospatial data including these variables were acquired from several open data sources and satellite and aerial imagery and, were processed in GIS software. Biomod2, an ensemble platform designed for species distribution modelling, was used to generate ensemble models in R. The proportion of built-up areas, field, forest, and snow-covered land in November, people working in the primary sector, human population density, mean precipitation in April and July, and densities of European hares, white-tailed deer, and raccoon dogs best estimated distribution of human TBE disease cases in the two datasets. Random forest and generalized boosted regression models performed with a very good to excellent predictive power (ROC = 0.89-0.96) in both time periods. Based on the predictive maps, high-risk areas for TBE transmission were located in the coastal regions in Southern and Western Finland (including the Åland Islands), several municipalities in Central and Eastern Finland, and coastal municipalities in Southern Lapland. To explore potential changes in TBE distributions in future climate, we used bioclimatic factors with current and future climate forecast data to reveal possible future hotspot areas. Based on the future forecasts, a slightly wider geographical extent of TBE risk was introduced in the Åland Islands and Southern, Western and Northern Finland, even though the risk itself was not increased. Our results are the first steps towards TBE-risk area mapping in current and future climate in Finland.

Changes in Protein Phosphorylation during Salivary Gland Degeneration in Haemaphysalis longicornis

The ticks feed large amount of blood from their hosts and transmit pathogens to the victims. The salivary gland plays an important role in the blood feeding. When the female ticks are near engorgement, the salivary gland gradually loses its functions and begins to rapidly degenerate. In this study, data-independent acquisition quantitative proteomics was used to study changes in the phosphorylation modification of proteins during salivary gland degeneration in Haemaphysalis longicornis. In this quantitative study, 400 phosphorylated proteins and 850 phosphorylation modification sites were identified. Trough RNA interference experiments, we found that among the proteins with changes in phosphorylation, apoptosis-promoting Hippo protein played a role in salivary gland degeneration.

Effect of climatic factors on the seasonal fluctuation of human brucellosis in Yulin, northern China

Background

Brucellosis is a serious public health problem primarily affecting livestock workers. The strong seasonality of the disease indicates that climatic factors may play important roles in the transmission of the disease. However, the associations between climatic variability and human brucellosis are still poorly understood.

Methods

Data for a 14-year series of human brucellosis cases and seven climatic factors were collected in Yulin City from 2005 to 2018, one of the most endemic areas in northern China. Using cross-correlation analysis, the Granger causality test, and a distributed lag non-linear model (DLNM), we assessed the quantitative relationships and exposure-lag-response effects between monthly climatic factors and human brucellosis.

Results

A total of 7103 cases of human brucellosis were reported from 2005 to 2018 in Yulin City with a distinct peak between April and July each year. Seasonal fluctuations in the transmission of human brucellosis were significantly affected by temperature, sunshine duration, and evaporation. The effects of climatic factors were non-linear over the 6-month period, and higher values of these factors usually increased disease incidence. The maximum separate relative risk (RR) was 1.36 (95% confidence interval [CI], 1.03–1.81) at a temperature of 17.4 °C, 1.12 (95% CI, 1.03–1.22) with 311 h of sunshine, and 1.18 (95% CI, 0.94–1.48) with 314 mm of evaporation. In addition, the effects of these three climatic factors were cumulative, with the highest RRs of 2.27 (95% CI, 1.09–4.57), 1.54 (95% CI, 1.10–2.18), and 1.27 (95% CI, 0.73–2.14), respectively.

Conclusions

In Yulin, northern China, variations in climatic factors, especially temperature, sunshine duration, and evaporation, contributed significantly to seasonal fluctuations of human brucellosis within 6 months. The key determinants of brucellosis transmission and the identified complex associations are useful references for developing strategies to reduce the disease burden.

Epidemiology of tick-borne encephalitis in China, 2007- 2018

Tick-borne encephalitis (TBE) is endemic to Europe and some Asian countries and is prevalent in northeast China. We analyzed the epidemiology of TBE in China from 2007 to 2018. A total of 3,364 TBE cases were reported in mainland China from 2007 to 2018, for an annual incidence of 0.09 to 0.44/100,000. Among the TBE cases, 89.92% were reported in forest areas (41.94% in DaXingAnLing, 8.70% in XiaoXingAnLing, and 39.21% in ChangBaiShan) in northeast China. The TBE cases were primarily male with a proportion of 67.15% (2,259/3,364 cases) and in 40–49-year age group with a proportion of 31.89% (1,073/3,364 cases). The epidemiology of TBE differed slightly among the three forest regions. Domestic workers and forestry workers accounted for the most of the TBE cases in DaXingAnLing, and domestic workers and farmers in XiaoXingAnLing and ChangBaiShan, respectively. The TBE cases mainly occurred from April to August with a peak in June. The TBE laboratory confirmed rate in DaXingAnLing (84.14%, 1,189/1,413 cases) was highest, compared with XiaoXingAnLing and ChangBaiShan (13.99% and 11.37%, respectively). Moreover, the hospital with the highest laboratory confirmed rate (88.01%, 1,336/1,518 cases) was Inner Mongolia Forestry General Hospital of DaXingAnling region. Systematic enhanced TBE surveillance and a vaccination program are needed to improve the laboratory confirmed rate and reduce the incidence of TBE in northeast China.

Cross protection induced by combined Subolesin-based DNA and protein immunizations against adult Haemaphysalis longicornis

Vaccination against ticks is an environmentally friendly alternative control method compared to chemical acaricide applications. Subolesin is a conserved protein in ticks, which can provide protection against some tick species. In this study, we evaluated the capacity of cocktail vaccination with Subolesin and ribosomal acidic protein 0 (P0) peptide against adults of Haemaphysalis longicornis. Priming with DNA vaccine expressing subolesin, followed by boosters of a single antigen (rRhSub) or a chimeric polypeptide (rRhSub/P0), provided cross protection. This treatment resulted in significant mortality, reduced blood ingestion and reduced reproduction in H. longicornis adults. Vaccination efficacies of 79.3% and 86.6% are reported in groups supplemented with rRhSub and rRhSub/P0, respectively. Conserved antigens, such as subolesin, formulated as DNA vaccine and enhanced with chimeric polypeptides, could be used as an anti-tick vaccine application, especially for control of infestation involving several tick species.

A dataset of distribution and diversity of ticks in China

While tick-borne zoonoses, such as Lyme disease and tick-borne encephalitis, present an increasing global concern, knowledge of their vectors' distribution remains limited, especially for China. In this paper, we present the first comprehensive dataset of known tick species and their distributions in China, derived from peer-reviewed literature published between 1960 and 2017. We searched for journal articles, conference papers and degree thesis published in both English and Chinese, extracted geographic information associated with tick occurrence, and applied quality-control procedures to remove duplicates and ensure accuracy. The dataset contains 5731 records of geo-referenced occurrences for 123 tick species distributed over 1141 locations distinguished at four levels of scale i.e., provincial, prefectural, county, and township and finer. The most frequently reported tick species include Haemaphysalis longicornis, Dermacentor silvarum, Ixodes persulcatus, Haemaphysalis conicinna, Rhipicephalus microplus, and Rhipicephalus sanguineus sensu lato. The geographical dataset provides an improved map of where ticks inhabit China and can be used for a variety of spatial analyses of ticks and the risk of zoonoses they transmit.

Epidemiological Characteristics and Spatial Analysis of Tick-Borne Encephalitis in Jilin Province, China

Tick-borne encephalitis (TBE) is a viral infectious disease and has become a reemerging public health threat in recent years in northeastern China. However, no studies has characterized the epidemiologic features and explored the spatial dynamics and environmental factors of TBE cases in Jilin Province. In this study, we have described the epidemiological features of 846 reported human TBE cases from 2006 to 2016 in Jilin Province. There was an obvious single peak pattern of TBE cases from May to July in Jilin Province. More than 60% of TBE cases occurred in farmers, and the people in 50- to 59-year-old group had the high incidence of the disease. The results of Getis-Ord Gi* statistics demonstrated that the human TBE cases were more clustered in the northeastern border including Dunhua and Yanji cities and Antu and Wangqing counties, and southern areas including Huinan, Jingyu, Jiangyuan, and Liuhe counties in Jilin Province. We demonstrated that the temporal dynamics of TBE in Jilin was significantly associated with the dynamics of meteorological factors especially after 2009. The results from the auto-logistic regression analysis showed that the percentage coverage of forest, temperature, and autoregressive term were significantly associated with the occurrence of human TBE cases in Jilin Province. Our findings will provide a scientific evidence for the targeted prevention and control programs.

Ixodid ticks and tick-borne encephalitis virus prevalence in the South Asian part of Russia (Republic of Tuva)

The most significant processes of arbovirus evolution can be expected to occur in the territories where ticks of different species cohabitate and at the boundaries of virus occurrence, where the probability of the appearance of new virus variants is high due to the possible shift in the main vectors and/or vertebrate hosts. One of the most interesting regions in this regard is the Republic of Tuva. Since most of its territory is covered by mountain ranges and intermountain basins, we were able to study the distribution of vectors and viruses in geographically isolated areas at different altitudes and in various landscapes. From 2008 to 2017, we conducted six expeditions to Tuva and collected 3,077 adult ticks and 24 nymphs. The distribution of tick species was confined to specific landscapes, as follows: Dermacentor nuttalli occurred in steppes, D. silvarum inhabited forest-steppe areas, and Ixodes persulcatus inhabited mixed forests. All three species of ticks were collected on plains and mountain slopes. The range of D. silvarum was shown to be lower than 1300 m above sea level (a.s.l.). Only D. nuttalli and I. persulcatus were collected at higher altitudes. According to our observations, single nymphs of D. nuttalli appear on animals one month before larvae appear. This finding confirms the hypothesis that the immature forms of D. nuttalli are able to overwinter under favourable conditions. We isolated 9 strains and 3 isolates of tick-borne encephalitis virus (TBEV) from I. persulcatus, one strain from D. nuttalli and one strain from D. silvarum. The TBEV strain from D. nuttalli was isolated from the territory inhabited only by Dermacentor ticks. All isolated strains belong to the Siberian subtype of TBEV. TBEV was detected in ticks from all the investigated altitudes. There were no statistically significant differences in the virus prevalence between the Dermacentor and Ixodes ticks. The results of our work provide additional support for the hypothesis of the existence of TBEV foci in areas with an absolute dominance of D. nuttalli.

Virome analysis of tick-borne viruses in Heilongjiang Province, China

Ticks are implicated in the transmission of various human and livestock pathogens worldwide. This study aimed to understand the geographical distribution of tick species, along with tick-associated viruses, in Heilongjiang Province, northeast China. Molecular methods were used to classify tick species, with next-generation sequencing and polymerase chain reaction-based analyses used to assess the viromes of ticks from four representative sampling locations in the Greater Khingan Mountains. Five species of ixodid ticks were identified, including Ixodes persulcatus, Dermacentor nuttalli, Dermacentor silvarum, Haemaphysalis longicornis, and Haemaphysalis concinna. From the 1102 ticks, 3,568,561 high-quality reads were obtained by next-generation sequencing. Following trimming, 302,540 reads were obtained, of which 6577 (2.16%) reads were annotated to viruses. Phylogenetic analysis revealed that the viral sequences shared a close relationship with Orthonairovirus, Phlebovirus, deer tick Mononegavirales-like virus, and Jingmen tick virus sequences, but the significance of these newly-identified tick-borne viruses to human and animal health requires further investigation. The results of this study provide a basis not only for further studies on the relationship between ticks and tick-borne viruses, but also for preventing future tick-borne epidemic outbreaks by means of vector control.

A new subtype of eastern tick-borne encephalitis virus discovered in Qinghai-Tibet Plateau, China

Tick-borne encephalitis virus (TBEV) has been classified into three subtypes, namely the European (Eu-TBEV), Far Eastern (FE-TBEV), and Siberian (Sib-TBEV). In this study, we discovered a new subtype of TBEV in wild rodent Marmota himalayana in Qinghai-Tibet Plateau in China, proposed as subtype Himalayan (Him-TBEV). Two complete genomes of TBEV were obtained from respiratory samples of 200 marmots. The phylogenetic analysis using the E protein and polyprotein demonstrated that the two strains of Him-TBEV formed an independent branch, separated from Eu-TBEV, Sib-TBEV, and FE-TBEV. The nomenclature of Him-TBEV as a new subtype was also supported by comparative analysis using nucleotide and amino acid sequences of E protein and polyprotein. For E protein, The Him-TBEV showed 82.6–84.6% nucleotide identities and 92.7–95.0% amino acid identities with other three subtypes. For polyprotein, the Him-TBEV showed 83.5–85.2% nucleotide identities and 92.6–94.2% amino acids identities with other three subtypes. Furthermore, of 69 amino acid substitutions profiles detected in complete polyprotein of 112 strains of TBEV, Him-TBEV subtype displayed unique amino acids in the 36 positions. Notably, for the subtype-specific amino acid position 206 of E protein, Him-TBEV shared the Val with Eu-TBEV, but differed from FE-TBEV and Sib-TBEV. The evolutionary analysis with BEAST suggested that Him-TBEV diverged from other subtypes of eastern TBEV group about 2469 years ago. It should be mentioned that Qinghai-Tibet Plateau in China is the plague endemic region where Marmota himalayana is the primary host. The public health significance of discovery of Him-TBEV in Marmota himalayana must be carefully evaluated.

Geographical distribution, climate adaptation and vector competence of the Eurasian hard tick Haemaphysalis concinna

The ixodid tick Haemaphysalis concinna Koch, 1844 is a proven vector of tick-borne encephalitis (TBE) virus and Francisella tularensis, the causative agent of tularaemia. In the present study, up-to-date maps depicting the geographical distribution and climate adaptation of H. concinna are presented. A dataset was compiled, resulting in 656 georeferenced locations in Eurasia. The distribution of H. concinna ranges from the Spanish Atlantic coast to Kamchatka, Russia, within the belt of 28-64° N latitude. H. concinna is the second most abundant tick species after Ixodes ricinus collected from birds, and third most abundant tick species flagged from vegetation in Central Europe. To investigate the climate adaptation of H. concinna, the georeferenced locations were superimposed on a high-resolution map of the Köppen-Geiger climate classification. A frequency distribution of the H. concinna occurrence under different climates shows three peaks related to the following climates: warm temperate with precipitation all year round, boreal with precipitation all year round and boreal, winter dry. Almost 87.3 % of all H. concinna locations collected are related to these climates. Thus, H. concinna prefers climates with a warm and moist summer. The remaining tick locations were characterized as cold steppes (6.2%), cold deserts (0.8%), Mediterranean climates (2.7%) or warm temperate climates with dry winter (2.9%). In those latter climates H. concinna occurs only sporadically, provided the microclimate is favourable. Beyond proven vector competence pathogen findings in questing H. concinna are compiled from the literature.

Tick-Borne Viruses

Ticks are important vectors for the transmission of pathogens including viruses. The viruses carried by ticks also known as tick-borne viruses (TBVs), contain a large group of viruses with diverse genetic properties and are concluded in two orders, nine families, and at least 12 genera. Some members of the TBVs are notorious agents causing severe diseases with high mortality rates in humans and livestock, while some others may pose risks to public health that are still unclear to us. Herein, we review the current knowledge of TBVs with emphases on the history of virus isolation and identification, tick vectors, and potential pathogenicity to humans and animals, including assigned species as well as the recently discovered and unassigned species. All these will promote our understanding of the diversity of TBVs, and will facilitate the further investigation of TBVs in association with both ticks and vertebrate hosts.