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

Crowdsourced tick observation data from across 60 years reveals major increases and northwards shifts in tick contact areas in Finland

There is mounting evidence of increases in tick (Acari: Ixodidae) contacts in Finland during the past few decades, highlighted by increases in the incidence of Lyme borreliosis and tick-borne encephalitis (TBE). While nationwide field studies to map distributions of ticks are not feasible, crowdsourcing provides a comprehensive method with which to assess large-scale changes in tick contact areas. Here, we assess changes in tick contact areas in Finland between 1958 and 2021 using three different nationwide crowdsourced data sets. The data revealed vast increases in tick contact areas, with ticks estimated to be contacted locally approximately 400 km further north in western and approximately 100 km further north in eastern Finland in 2021 than 1958. Tick contact rates appeared to be highest along the coastline and on the shores of large lakes, possibly indicating higher tick abundance therein. In general, tick observations per inhabitant increased from 2015 to 2021. Tick contact areas have expanded in Finland over the past 60 years. It appears that taiga ticks (Ixodes persulcatus) are behind most of the northwards shifts in tick contact areas, with Ixodes ricinus contributing mostly to new contact areas in the south. While ticks are now present in most of Finland, there are still areas where tick abundance is low and/or establishment not possible, mainly in northern Finland.

Outdoor recreation, tick borne encephalitis incidence and seasonality in Finland, Norway and Sweden during the COVID-19 pandemic (2020/2021)

During the pandemic outdoor activities were encouraged to mitigate transmission risk while providing safe spaces for social interactions. Human behaviour, which may favour or disfavour, contact rates between questing ticks and humans, is a key factor impacting tick-borne encephalitis (TBE) incidence. We analyzed annual and weekly TBE cases in Finland, Norway and Sweden from 2010 to 2021 to assess trend, seasonality, and discuss changes in human tick exposure imposed by COVID-19. We compared the pre-pandemic incidence (2010-2019) with the pandemic incidence (2020-2021) by fitting a generalized linear model (GLM) to incidence data. Pre-pandemic incidence was 1.0, 0.29 and 2.8 for Finland, Norway and Sweden, respectively, compared to incidence of 2.2, 1.0 and 3.9 during the pandemic years. However, there was an increasing trend for all countries across the whole study period. Therefore, we predicted the number of cases in 2020/2021 based on a model fitted to the annual cases in 2010-2019. The incidences during the pandemic were 1.3 times higher for Finland, 1.7 times higher for Norway and no difference for Sweden. When social restrictions were enforced to curb the spread of SARS-CoV-2 there were profound changes in outdoor recreational behavior. Future consideration of public health interventions that promote outdoor activities may increase exposure to vector-borne diseases.

Ecological and environmental factors affecting the risk of tick-borne encephalitis in Europe, 2017 to 2021

BackgroundTick-borne encephalitis (TBE) is a disease which can lead to severe neurological symptoms, caused by the TBE virus (TBEV). The natural transmission cycle occurs in foci and involves ticks as vectors and several key hosts that act as reservoirs and amplifiers of the infection spread. Recently, the incidence of TBE in Europe has been rising in both endemic and new regions.AimIn this study we want to provide comprehensive understanding of the main ecological and environmental factors that affect TBE spread across Europe.MethodsWe searched available literature on covariates linked with the circulation of TBEV in Europe. We then assessed the best predictors for TBE incidence in 11 European countries by means of statistical regression, using data on human infections provided by the European Surveillance System (TESSy), averaged between 2017 and 2021.ResultsWe retrieved data from 62 full-text articles and identified 31 different covariates associated with TBE occurrence. Finally, we selected eight variables from the best model, including factors linked to vegetation cover, climate, and the presence of tick hosts.DiscussionThe existing literature is heterogeneous, both in study design and covariate types. Here, we summarised and statistically validated the covariates affecting the variability of TBEV across Europe. The analysis of the factors enhancing disease emergence is a fundamental step towards the identification of potential hotspots of viral circulation. Hence, our results can support modelling efforts to estimate the risk of TBEV infections and help decision-makers implement surveillance and prevention campaigns.

Potential drivers of human tick-borne encephalitis in the Örebro region of Sweden, 2010-2021

Incidence of tick-borne encephalitis (TBE) has increased during the last years in Scandinavia, but the underlying mechanism is not understood. TBE human case data reported between 2010 and 2021 were aggregated into postal codes within Örebro County, south-central Sweden, along with tick abundance and environmental data to analyse spatial patterns and identify drivers of TBE. We identified a substantial and continuing increase of TBE incidence in Örebro County during the study period. Spatial cluster analyses showed significant hotspots (higher number of cases than expected) in the southern and northern parts of Örebro County, whereas a cold spot (lower number of cases than expected) was found in the central part comprising Örebro municipality. Generalised linear models showed that the risk of acquiring TBE increased by 12.5% and 72.3% for every percent increase in relative humidity and proportion of wetland forest, respectively, whereas the risk decreased by 52.8% for every degree Celsius increase in annual temperature range. However, models had relatively low goodness of fit (R² < 0.27). Results suggest that TBE in Örebro County is spatially clustered, however variables used in this study, i.e., climatic variables, forest cover, water, tick abundance, sheep as indicator species, alone do not explain this pattern.

Leptospirosis modelling using hydrometeorological indices and random forest machine learning

Leptospirosis is a zoonosis that has been linked to hydrometeorological variability. Hydrometeorological averages and extremes have been used before as drivers in the statistical prediction of disease. However, their importance and predictive capacity are still little known. In this study, the use of a random forest classifier was explored to analyze the relative importance of hydrometeorological indices in developing the leptospirosis model and to evaluate the performance of models based on the type of indices used, using case data from three districts in Kelantan, Malaysia, that experience annual monsoonal rainfall and flooding. First, hydrometeorological data including rainfall, streamflow, water level, relative humidity, and temperature were transformed into 164 weekly average and extreme indices in accordance with the Expert Team on Climate Change Detection and Indices (ETCCDI). Then, weekly case occurrences were classified into binary classes "high" and "low" based on an average threshold. Seventeen models based on "average," "extreme," and "mixed" indices were trained by optimizing the feature subsets based on the model computed mean decrease Gini (MDG) scores. The variable importance was assessed through cross-correlation analysis and the MDG score. The average and extreme models showed similar prediction accuracy ranges (61.5-76.1% and 72.3-77.0%) while the mixed models showed an improvement (71.7-82.6% prediction accuracy). An extreme model was the most sensitive while an average model was the most specific. The time lag associated with the driving indices agreed with the seasonality of the monsoon. The rainfall variable (extreme) was the most important in classifying the leptospirosis occurrence while streamflow was the least important despite showing higher correlations with leptospirosis.

The range and habitat suitability of François' langur (Trachypithecus francoisi) in Mayanghe Nature Reserve, China

The François' langur (Trachypithecus francoisi) is an endangered primate living in limestone forests in Vietnam and China. From October 2017 to August 2018, the habitat preferences and the range of the Francois' langur were surveyed in the Mayanghe Nature Reserve, Guizhou, People's Republic of China. To estimate the range and predict suitable habitat of François' langur, a Gaussian normal kernel density estimation and species distribution models (BIOMOD2) were used along with data on environmental variables and records of the langur's occurrence. The total range of François' langur in the reserve is 68.76 km², accounting for 22.1% of the total area of the reserve. The elevation of the main utilisation area is 500-800 m, accounting for 48.53% of the total area of the reserve. The maximum slope utilised is 20°-30°, 30.19 km² and accounting for 30.56% of the total area. The habitat used is largely distributed along valleys, preferred broad leaf forest, lower elevation, and close to rivers. Broad leaf forest is the main habitat type utilised, totalling 25.57 km² and accounting for 37.19% of the total area. Our models predicted that the suitable habitat in the reserve is 62.46 km², accounting for 20.08% of the total reserve area, with 32.93-km² suitable habitat occurring in the core zone, 22.44 km² in the buffer zone, and 7.02 km² in the experimental zone. Our results indicate that only limited suitable habitat (51%) adverse reserve zoning exists in the Mayanghe Nature Reserve of François' langur.

Decoding the Geography of Natural TBEV Microfoci in Germany: A Geostatistical Approach Based on Land-Use Patterns and Climatological Conditions

Background: Tickborne-encephalitis (TBE) is a potentially life-threating neurological disease that is mainly transmitted by ticks. The goal of the present study is to analyze the potential uniform environmental patterns of the identified TBEV microfoci in Germany. The results are used to calculate probabilities for the present distribution of TBEV microfoci in Germany based on a geostatistical model. Methods: We aim to consider the specification of environmental characteristics of locations of TBEV microfoci detected in Germany using open access epidemiological, geographical and climatological data sources. We use a two-step geostatistical approach, where in a first step, the characteristics of a broad set of environmental variables between the 56 TBEV microfoci and a control or comparator set of 3575 sampling points covering Germany are compared using Fisher's Exact Test. In the second step, we select the most important variables, which are then used in a MaxEnt distribution model to calculate a high resolution (400 × 400 m) probability map for the presence of TBEV covering the entire area of Germany. Results: The findings from the MaxEnt prediction model indicate that multi annual actual evapotranspiration (27.0%) and multi annual hot days (22.5%) have the highest contribution to our model. These two variables are followed by four additional variables with a lower, but still important, explanatory influence: Land cover classes (19.6%), multi annual minimum air temperature (14.9%), multi annual sunshine duration (9.0%), and distance to coniferous and mixed forest border (7.0%). Conclusions: Our findings are based on defined TBEV microfoci with known histories of infection and the repeated confirmation of the virus in the last years, resulting in an in-depth high-resolution model/map of TBEV microfoci in Germany. Multi annual actual evapotranspiration (27%) and multi annual hot days (22.5%) have the most explanatory power in our model. The results may be used to tailor specific regional preventive measures and investigations.

Climate Changes Exacerbate the Spread of Ixodes ricinus and the Occurrence of Lyme Borreliosis and Tick-Borne Encephalitis in Europe-How Climate Models Are Used as a Risk Assessment Approach for Tick-Borne Diseases

Climate change has influenced the transmission of a wide range of vector-borne diseases in Europe, which is a pressing public health challenge for the coming decades. Numerous theories have been developed in order to explain how tick-borne diseases are associated with climate change. These theories include higher proliferation rates, extended transmission season, changes in ecological balances, and climate-related migration of vectors, reservoir hosts, or human populations. Changes of the epidemiological pattern have potentially catastrophic consequences, resulting in increasing prevalence of tick-borne diseases. Thus, investigation of the relationship between climate change and tick-borne diseases is critical. In this regard, climate models that predict the ticks’ geographical distribution changes can be used as a predicting tool. The aim of this review is to provide the current evidence regarding the contribution of the climatic changes to Lyme borreliosis (LB) disease and tick-borne encephalitis (TBE) and to present how computational models will advance our understanding of the relationship between climate change and tick-borne diseases in Europe.

Habitat suitability map of <em>Ixodes ricinus</em> tick in France using multi-criteria analysis

The tick Ixodes ricinus is widely distributed across Europe and is responsible for the transmission of several pathogens to humans and animals. In this study, we used a knowledge-based method to map variations in habitat suitability for I. ricinus ticks throughout continental France and Corsica. The multi-criteria decision analysis (MCDA) integrated four major biotic and abiotic factors known to influence tick populations: climate, land cover, altitude and the density of wild ungulates. For each factor, habitat suitability index (HSI) values were attributed to different locations based on knowledge regarding its impact on tick populations. For the MCDA, two methods of factor combination were tested, additive and multiplicative, both which were evaluated at the spatial scales of departments and local municipalities. The resulting habitat suitability maps (resolution=100x100 m) revealed that conditions are suitable for I. ricinus over most of France and Corsica. Particularly suitable habitats were located in central, north-eastern and south-western France, while less-suitable habitats were found in the Mediterranean and mountainous regions. To validate the approach, the HSI scores were compared to field data of I. ricinus nymph abundance. Regardless of scale, the correlation between abundance indicator and HSI score was stronger for the additive than for the multiplicative approach. Overall, this study demonstrates the value of MCDA for estimating habitat suitability maps for I. ricinus abundance, which could be especially useful in highlighting areas of the tick's distribution where preventive measures should be prioritised.

Changes in the Epidemiology of Zoonotic Infections in Children: A Nationwide Register Study in Finland

OBJECTIVES

Zoonotic infections are difficult to recognize in children. The age distributions and seasonal occurrences of these infections vary substantially, even among those transmitted by the same vectors, and their epidemiology may change over time. The aim was to report the incidences and trends of Borrelia burgdorferi, Puumala virus, Francisella tularensis and tick-borne encephalitis (TBE) virus infections in the pediatric population (age 0-19) of Finland.

METHODS

A nationwide survey based on the National Infectious Disease Register was conducted from 1996 to 2019 and all laboratory-confirmed cases were included. Age-stratified incidences per 100,000 person-years were calculated.

RESULTS

Cumulative incidences were B. burgdorferi 11.2, TBE 0.4, Puumala virus 6.4 and F. tularensis 2.5 per 100,000 person-years. An increasing trend in the incidences of B. burgdorferi and TBE was observed. Borrelia expanded geographically northward and inland. Tularemia follows a 2-4-year epidemic cycle and rates are similar across age groups. Puumala incidences are highest in the older children.

DISCUSSION

Borrelia infections increased most rapidly in children 5-9 years of age and overall expanded geographically in Finland. Tularemia epidemic cycles were shorter than those previously reported. These results will help clinicians to identify these infections in different geographic areas and age groups in Finland.

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.

Associating Land Cover Changes with Patterns of Incidences of Climate-Sensitive Infections: An Example on Tick-Borne Diseases in the Nordic Area

Some of the climate-sensitive infections (CSIs) affecting humans are zoonotic vector-borne diseases, such as Lyme borreliosis (BOR) and tick-borne encephalitis (TBE), mostly linked to various species of ticks as vectors. Due to climate change, the geographical distribution of tick species, their hosts, and the prevalence of pathogens are likely to change. A recent increase in human incidences of these CSIs in the Nordic regions might indicate an expansion of the range of ticks and hosts, with vegetation changes acting as potential predictors linked to habitat suitability. In this paper, we study districts in Fennoscandia and Russia where incidences of BOR and TBE have steadily increased over the 1995-2015 period (defined as 'Well Increasing districts'). This selection is taken as a proxy for increasing the prevalence of tick-borne pathogens due to increased habitat suitability for ticks and hosts, thus simplifying the multiple factors that explain incidence variations. This approach allows vegetation types and strengths of correlation specific to the WI districts to be differentiated and compared with associations found over all districts. Land cover types and their changes found to be associated with increasing human disease incidence are described, indicating zones with potential future higher risk of these diseases. Combining vegetation cover and climate variables in regression models shows the interplay of biotic and abiotic factors linked to CSI incidences and identifies some differences between BOR and TBE. Regression model projections up until 2070 under different climate scenarios depict possible CSI progressions within the studied area and are consistent with the observed changes over the past 20 years.

Public health preparedness and response synergies between institutional authorities and the community: a qualitative case study of emerging tick-borne diseases in Spain and the Netherlands

BACKGROUND

Communities affected by infectious disease outbreaks are increasingly recognised as partners with a significant role to play during public health emergencies. This paper reports on a qualitative case study of the interactions between affected communities and public health institutions prior to, during, and after two emerging tick-borne disease events in 2016: Crimean-Congo Haemorrhagic Fever in Spain, and Tick-Borne Encephalitis in the Netherlands. The aim of the paper is to identify pre-existing and emergent synergies between communities and authorities, and to highlight areas where synergies could be facilitated and enhanced in future outbreaks.

METHODS

Documentary material provided background for a set of semi-structured interviews with experts working in both health and relevant non-health official institutions (13 and 21 individuals respectively in Spain and the Netherlands), and focus group discussions with representatives of affected communities (15 and 10 individuals respectively). Data from all sources were combined and analysed thematically, initially independently for each country and then for both countries together.

RESULTS

Strong synergies were identified in tick surveillance activities in both countries, and the value of pre-existing networks of interest groups for preparedness and response activities was recognised. However, authorities also noted that there were hard-to-reach and potentially vulnerable groups, such as hikers, foreign tourists, and volunteers working in green areas. While the general population received preventive information about the two events, risk communication or other community engagement efforts were not seen as necessary specifically for these sub-groups. Post-event evaluations of community engagement activities during the two events were limited, so lessons learned were not well documented.

CONCLUSIONS

A set of good practices emerged from this study, that could be applied in these and other settings. They included the potential value of conducting stakeholder analyses of community actors with a stake in tick-borne or other zoonotic diseases; of utilising pre-existing stakeholder networks for information dissemination; and of monitoring community perceptions of any public health incident, including through social media. Efforts in the two countries to build on the community engagement activities that are already in place could contribute to better preparedness planning and more efficient and timely responses in future outbreaks.

Projecting the potential distribution of ticks in China under climate and land use change

Ticks are known as vectors of several pathogens causing various human and animal diseases including Lyme borreliosis, tick-borne encephalitis, and Crimean-Congo hemorrhagic fever. While China is known to have more than 100 tick species well distributed over the country, our knowledge on the likely distribution of ticks in the future remains very limited, which hinders the prevention and control of the risk of tick-borne diseases. In this study, we selected four representative tick species which have different regional distribution foci in mainland China. i.e., Dermacentor marginatus, Dermacentor silvarum, Haemaphysalis longicornis and Ixodes granulatus. We used the MaxEnt model to identify the key environmental factors of tick occurrence and map their potential distributions in 2050 under four combined climate and socioeconomic scenarios (i.e., SSP1-RCP2.6, SSP2-RCP4.5, SSP3-RCP7.0 and SSP5-RCP8.5). We found that the extent of the urban fabric, cropland and forest, temperature annual range and precipitation of the driest month were the main determinants of the potential distributions of the four tick species. Under the combined scenarios, with climate warming, the potential distributions of ticks shifted to further north in China. Due to a decrease in the extent of forest, the distribution probability of ticks declined in central and southern China. In contrast with previous findings on an estimated amplification of tick distribution probability under the extreme emission scenario (RCP8.5), our studies projected an overall reduction in the distribution probability under RCP8.5, owing to an expected effect of land use. Our results could provide new data to help identify the emerging risk areas, with amplifying suitability for tick occurrence, for the prevention and control of tick-borne zoonoses in mainland China. Future directions are suggested towards improved quantity and quality of the tick occurrence database, comprehensiveness of factors and integration of different modelling approaches, and capability to model pathogen spillover at the human-tick interface.

Game Animal Density, Climate, and Tick-Borne Encephalitis in Finland, 2007-2017

Tick-borne encephalitis (TBE) is an endemic infection of public health importance in Finland. We investigated the effect of ecologic factors on 2007-2017 TBE trends. We obtained domestic TBE case data from the National Infectious Diseases Register, weather data from the US National Oceanic and Atmospheric Administration, and data from the Natural Resources Institute in Finland on mammals killed by hunters yearly in game management areas. We performed a mixed-effects time-series analysis with time lags on weather and animal parameters, adding a random effect to game management areas. During 2007-2017, a total of 395/460 (86%) domestic TBE cases were reported with known place of exposure and date of sampling. Overall, TBE incidence increased yearly by 15%. After adjusting for the density of other animals and minimum temperatures, we found thatTBE incidence was positively associated with white-tailed deer density. Variation in host animal density should be considered when assessing TBE risks and designing interventions.

Tick-Borne Encephalitis Virus and Its European Distribution in Ticks and Endothermic Mammals

Tick-borne encephalitis (TBE) is the most common viral tick-borne disease in Europe causing thousands of human infections every year. Available risk maps in Europe are solely based on human incidences, but often underestimate areas with TBE virus circulation as shown by several autochthonous cases detected outside known risk areas. A dataset of more than 1300 georeferenced TBE virus detections in ticks and mammals except for humans was compiled and used to estimate the probability of TBE virus presence in Europe. For this, a random forests model was implemented using temperature- and precipitation-dependent bioclimatic variables of the WorldClim dataset, altitude, as well as land cover of the ESA GlobCover dataset. The highest probabilities of TBE virus presence were identified in Central Europe, in the south of the Nordic countries, and in the Baltic countries. The model performance was evaluated by an out-of-bag error (OOB) of 0.174 and a high area under the curve value (AUC) of 0.905. The TBE virus presence maps may subsequently be used to estimate the risk of TBE virus infections in humans and can support decision-makers to identify TBE risk areas and to encourage people to take appropriate actions against tick bites and TBE virus infections.

Mapping the Potential Distribution of Major Tick Species in China

Ticks are known as the vectors of various zoonotic diseases such as Lyme borreliosis and tick-borne encephalitis. Though their occurrences are increasingly reported in some parts of China, our understanding of the pattern and determinants of ticks' potential distribution over the country remain limited. In this study, we took advantage of the recently compiled spatial dataset of distribution and diversity of ticks in China, analyzed the environmental determinants of ten frequently reported tick species and mapped the spatial distribution of these species over the country using the MaxEnt model. We found that presence of urban fabric, cropland, and forest in a place are key determents of tick occurrence, suggesting ticks were likely inhabited close to where people live. Besides, precipitation in the driest month was found to have a relatively high contribution in mapping tick distribution. The model projected that theses ticks could be widely distributed in the Northwest, Central North, Northeast, and South China. Our results added new evidence on the potential distribution of a variety of major tick species in China and pinpointed areas with a high potential risk of tick bites and tick-borne diseases for raising public health awareness and prevention responses.