Seasonal cycles of the TBE and Lyme borreliosis vector Ixodes ricinus modelled with time-lagged and interval-averaged predictors

Correlation between airborne pollen data and the risk of tick-borne encephalitis in northern Italy

Tick-borne encephalitis (TBE) is caused by a flavivirus that infects animals including humans. In Europe, the TBE virus circulates enzootically in natural foci among ticks and rodent hosts. The abundance of ticks depends on the abundance of rodent hosts, which in turn depends on the availability of food resources, such as tree seeds. Trees can exhibit large inter-annual fluctuations in seed production (masting), which influences the abundance of rodents the following year, and the abundance of nymphal ticks two years later. Thus, the biology of this system predicts a 2-year time lag between masting and the incidence of tick-borne diseases such as TBE. As airborne pollen abundance is related to masting, we investigated whether inter-annual variation in pollen load could be directly correlated with inter-annual variation in the incidence of TBE in human populations with a 2-year time lag. We focused our study on the province of Trento (northern Italy), where 206 TBE cases were notified between 1992 and 2020. We tested the relationship between TBE incidence and pollen load collected from 1989 to 2020 for 7 different tree species common in our study area. Through univariate analysis we found that the pollen quantities recorded two years prior for two tree species, hop-hornbeam (Ostrya carpinifolia) and downy oak (Quercus pubescens), were positively correlated with TBE emergence (R² = 0.2) while a multivariate model with both tree species better explained the variation in annual TBE incidence (R² = 0.34). To the best of our knowledge, this is the first attempt at quantifying the correlation between pollen quantities and the incidence of TBE in human populations. As pollen loads are collected by widespread aerobiological networks using standardized procedures, our study could be easily replicated to test their potential as early warning system for TBE and other tick-borne diseases.

Revealing the Tick Microbiome: Insights into Midgut and Salivary Gland Microbiota of Female Ixodes ricinus Ticks

The ectoparasite Ixodes ricinus is an important vector for many tick-borne diseases (TBD) in the northern hemisphere, such as Lyme borreliosis, rickettsiosis, human granulocytic anaplasmosis, or tick-borne encephalitis virus. As climate change will lead to rising temperatures in the next years, we expect an increase in tick activity, tick population, and thus in the spread of TBD. Consequently, it has never been more critical to understand relationships within the microbial communities in ticks that might contribute to the tick's fitness and the occurrence of TBD. Therefore, we analyzed the microbiota in different tick tissues such as midgut, salivary glands, and residual tick material, as well as the microbiota in complete Ixodes ricinus ticks using 16S rRNA gene amplicon sequencing. By using a newly developed DNA extraction protocol for tick tissue samples and a self-designed mock community, we were able to detect endosymbionts and pathogens that have been described in the literature previously. Further, this study displayed the usefulness of including a mock community during bioinformatic analysis to identify essential bacteria within the tick.

The Emergence and Dynamics of Tick-Borne Encephalitis Virus in a New Endemic Region in Southern Germany

Tick-borne encephalitis (TBE) is the most important viral tick-borne infection in Europe and Asia. It is emerging in new areas. The mechanisms of emergence are fairly unknown or speculative. In the Ravensburg district in southern Germany, TBE emerged, mainly over the last five years. Here, we analyzed the underlying epidemiology in humans. The resulting identified natural foci of the causal TBE virus (TBEV) were genetically characterized. We sampled 13 potential infection sites at these foci and detected TBEV in ticks (Ixodes ricinus) at eight sites. Phylogenetic analysis spurred the introduction of at least four distinct TBEV lineages of the European subtype into the Ravensburg district over the last few years. In two instances, a continuous spread of these virus strains over up to 10 km was observed.

Predicting habitat suitability for Ixodes ricinus and Ixodes persulcatus ticks in Finland

Background

Ticks are responsible for transmitting several notable pathogens worldwide. Finland lies in a zone where two human-biting tick species co-occur: Ixodesricinus and Ixodespersulcatus. Tick densities have increased in boreal regions worldwide during past decades, and tick-borne pathogens have been identified as one of the major threats to public health in the face of climate change.

Methods

We used species distribution modelling techniques to predict the distributions of I.ricinus and I.persulcatus, using aggregated historical data from 2014 to 2020 and new tick occurrence data from 2021. By aiming to fill the gaps in tick occurrence data, we created a new sampling strategy across Finland. We also screened for tick-borne encephalitis virus (TBEV) and Borrelia from the newly collected ticks. Climate, land use and vegetation data, and population densities of the tick hosts were used in various combinations on four data sets to estimate tick species’ distributions across mainland Finland with a 1-km resolution.

Results

In the 2021 survey, 89 new locations were sampled of which 25 new presences and 63 absences were found for I.ricinus and one new presence and 88 absences for I.persulcatus. A total of 502 ticks were collected and analysed; no ticks were positive for TBEV, while 56 (47%) of the 120 pools, including adult, nymph, and larva pools, were positive for Borrelia (minimum infection rate 11.2%, respectively). Our prediction results demonstrate that two combined predictor data sets based on ensemble mean models yielded the highest predictive accuracy for both I.ricinus (AUC = 0.91, 0.94) and I.persulcatus (AUC = 0.93, 0.96). The suitable habitats for I.ricinus were determined by higher relative humidity, air temperature, precipitation sum, and middle-infrared reflectance levels and higher densities of white-tailed deer, European hare, and red fox. For I.persulcatus, locations with greater precipitation and air temperature and higher white-tailed deer, roe deer, and mountain hare densities were associated with higher occurrence probabilities. Suitable habitats for I.ricinus ranged from southern Finland up to Central Ostrobothnia and North Karelia, excluding areas in Ostrobothnia and Pirkanmaa. For I.persulcatus, suitable areas were located along the western coast from Ostrobothnia to southern Lapland, in North Karelia, North Savo, Kainuu, and areas in Pirkanmaa and Päijät-Häme.

Conclusions

This is the first study conducted in Finland that estimates potential tick species distributions using environmental and host data. Our results can be utilized in vector control strategies, as supporting material in recommendations issued by public health authorities, and as predictor data for modelling the risk for tick-borne diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05410-8.

Meteorological and climatic variables predict the phenology of Ixodes ricinus nymph activity in France, accounting for habitat heterogeneity

Ixodes ricinus ticks (Acari: Ixodidae) are the most important vector for Lyme borreliosis in Europe. As climate change might affect their distributions and activities, this study aimed to determine the effects of environmental factors, i.e., meteorological, bioclimatic, and habitat characteristics on host-seeking (questing) activity of I. ricinus nymphs, an important stage in disease transmissions, across diverse climatic types in France over 8 years. Questing activity was observed using a repeated removal sampling with a cloth-dragging technique in 11 sampling sites from 7 tick observatories from 2014 to 2021 at approximately 1-month intervals, involving 631 sampling campaigns. Three phenological patterns were observed, potentially following a climatic gradient. The mixed-effects negative binomial regression revealed that observed nymph counts were driven by different interval-average meteorological variables, including 1-month moving average temperature, previous 3-to-6-month moving average temperature, and 6-month moving average minimum relative humidity. The interaction effects indicated that the phenology in colder climates peaked differently from that of warmer climates. Also, land cover characteristics that support the highest baseline abundance were moderate forest fragmentation with transition borders with agricultural areas. Finally, our model could potentially be used to predict seasonal human-tick exposure risks in France that could contribute to mitigating Lyme borreliosis risk.

Model-based extrapolation of ecological systems under future climate scenarios: The example of Ixodes ricinus ticks

Models can be applied to extrapolate consequences of climate change for complex ecological systems in the future. The acknowledged systems' behaviour at present is projected into the future considering climate projection data. Such an approach can be used to addresses the future activity and density of the castor bean tick Ixodes ricinus, the most widespread tick species in Europe. It is an important vector of pathogens causing Lyme borreliosis and tick-borne encephalitis. The population dynamics depend on several biotic and abiotic factors. Such complexity makes it difficult to predict the future dynamics and density of I. ricinus and associated health risk for humans. The objective of this study is to force ecological models with high-resolution climate projection data to extrapolate I. ricinus tick density and activity patterns into the future. We used climate projection data of temperature, precipitation, and relative humidity for the period 1971-2099 from 15 different climate models. Tick activity was investigated using a climate-driven cohort-based population model. We simulated the seasonal population dynamics using climate data between 1971 and 2099 and observed weather data since 1949 at a specific location in southern Germany. We evaluated derived quantities of local tick ecology, e.g. the maximum questing activity of the nymphal stage. Furthermore, we predicted spatial density changes by extrapolating a German-wide tick density model. We compared the tick density of the reference period (1971-2000) with the counter-factual densities under the near-term scenario (2012-2041), mid-term scenario (2050-2079) and long-term scenario (2070-2099). We found that the nymphal questing peak would shift towards early seasons of the year. Also, we found high spatial heterogeneity across Germany, with predicted hotspots of up to 2,000 nymphs per 100 m2 and coldspots with constant density. As our results suggest extreme changes in tick behaviour and density, we discuss why caution is needed when extrapolating climate data-driven models into the distant future when data on future climate drive the model projection.

Beech tree masting explains the inter-annual variation in the fall and spring peaks of Ixodes ricinus ticks with different time lags

BACKGROUND

The tick Ixodes ricinus is an important vector of tick-borne diseases including Lyme borreliosis. In continental Europe, the nymphal stage of I. ricinus often has a bimodal phenology with a large spring peak and a smaller fall peak. There is consensus about the origin of the spring nymphal peak, but there are two alternative hypotheses for the fall nymphal peak. In the direct development hypothesis, larvae quest as nymphs in the fall of the same year that they obtained their larval blood meal. In the developmental diapause hypothesis, larvae overwinter in the engorged state and quest as nymphs one year after they obtained their larval blood meal. These two hypotheses make different predictions about the time lags that separate the larval blood meal and the density of questing nymphs (DON) in the spring and fall.

METHODS

Inter-annual variation in seed production (masting) by deciduous trees is a time-lagged index for the density of vertebrate hosts (e.g., rodents) which provide blood meals for larval ticks. We used a long-term data set on the masting of the European beech tree and a 15-year study on the DON at 4 different elevation sites in western Switzerland to differentiate between the two alternative hypotheses for the origin of the fall nymphal peak.

RESULTS

Questing I. ricinus nymphs had a bimodal phenology at the three lower elevation sites, but a unimodal phenology at the top elevation site. At the lower elevation sites, the DON in the fall was strongly correlated with the DON in the spring of the following year. The inter-annual variation in the densities of I. ricinus nymphs in the fall and spring was best explained by a 1-year versus a 2-year time lag with the beech tree masting index. Fall nymphs had higher fat content than spring nymphs indicating that they were younger. All these observations are consistent with the direct development hypothesis for the fall peak of I. ricinus nymphs at our study site. Our study provides new insight into the complex bimodal phenology of this important disease vector.

CONCLUSIONS

Public health officials in Europe should be aware that following a strong mast year, the DON will increase 1 year later in the fall and 2 years later in the spring. Studies of I. ricinus populations with a bimodal phenology should consider that the spring and fall peak in the same calendar year represent different generations of ticks.

The complex interplay of climate, TBEV vector dynamics and TBEV infection rates in ticks-Monitoring a natural TBEV focus in Germany, 2009-2018

BACKGROUND

Tick-borne encephalitis (TBE) is the most important tick-borne viral disease in Eurasia and causes disease in humans and in a number of animals, among them dogs and horses. There is still no good correlation between tick numbers, weather conditions and human cases. There is the hypothesis that co-feeding due to simultaneous occurrence of larvae and nymphs may be a factor for the increased transmission of the virus in nature and for human disease. Based on long-term data from a natural TBEV focus, phylogenetic results and meteorological data we sought to challenge this hypothesis.

METHODS

Ticks from an identified TBE natural focus were sampled monthly from 04/2009 to 12/2018. Ticks were identified and pooled. Pools were tested by RT-qPCR. Positive pools were confirmed by virus isolation and/or sequencing of additional genes (E gene, NS2 gene). Temperature data such as the decadal (10-day) mean daily maximum air temperature (DMDMAT) were obtained from a nearby weather station and statistical correlations between tick occurrence and minimal infection rates (MIR) were calculated.

RESULTS

In the study period from 04/2009 to 12/2018 a total of 15,530 ticks (2,226 females, 2,268 males, 11,036 nymphs) were collected. The overall MIR in nymphs over the whole period was 77/15,530 (0.49%), ranging from 0.09% (2009) to 1.36% (2015). The overall MIR of female ticks was 0.76% (17/2,226 ticks), range 0.14% (2013) to 3.59% (2016). The overall MIR of males was 0.57% (13/2,268 ticks), range from 0.26% (2009) to 0.97% (2015). The number of nymphs was statistically associated with a later start of spring/vegetation period, indicated by the onset of forsythia flowering.

CONCLUSION

There was no particular correlation between DMDMAT dynamics in spring and/or autumn and the MIR of nymphs or adult ticks detected. However, there was a positive correlation between the number of nymphs and the number of reported human TBE cases in the following months, but not in the following year. The hypothesis of the importance of co-feeding of larvae and nymphs for the maintenance of transmission cycle of TBEV in nature is not supported by our findings.

Long-term monitoring of the seasonal density of questing ixodid ticks in Vienna (Austria): setup and first results

The first long-term monitoring to document both activity and density of questing ixodid ticks in Vienna, Austria, is introduced. It was started in 2017 and is planned to run over decades. Such long-term monitorings are needed to quantify possible effects of climate change or to develop tick density forecast models. The monthly questing tick density at three sites has been observed by using a standardized sampling method by dragging an area of [Formula: see text] at each occasion. Popular recreational areas were chosen as study sites. These are the Prater public park, the wooded Kahlenberg, and a wildlife garden in Klosterneuburg. First results show a 3-year time series of nymphs and adults of the Ixodes ricinus species complex and Haemaphysalis concinna for the period 2017-2019. Whereas questing nymphs of the I. ricinus species complex were collected from February to November, H. concinna nymphs were only dragged from May to October. The peak of nymphal activity of the I. ricinus species complex was in May, that of H. concinna in August. In addition, a brief overview is given about ticks and tick-borne pathogens occurring in urban and suburban areas of Vienna.

Tick-borne encephalitis (TBE) cases are not random: explaining trend, low- and high-frequency oscillations based on the Austrian TBE time series

BACKGROUND

Why human tick-borne encephalitis (TBE) cases differ from year to year, in some years more 100%, has not been clarified, yet. The cause of the increasing or decreasing trends is also controversial. Austria is the only country in Europe where a 40-year TBE time series and an official vaccine coverage time series are available to investigate these open questions.

METHODS

A series of generalized linear models (GLMs) has been developed to identify demographic and environmental factors associated with the trend and the oscillations of the TBE time series. Both the observed and the predicted TBE time series were subjected to spectral analysis. The resulting power spectra indicate which predictors are responsible for the trend, the high-frequency and the low-frequency oscillations, and with which explained variance they contribute to the TBE oscillations.

RESULTS

The increasing trend can be associated with the demography of the increasing human population. The responsible GLM explains 12% of the variance of the TBE time series. The low-frequency oscillations (10 years) are associated with the decadal changes of the large-scale climate in Central Europe. These are well described by the so-called Scandinavian index. This 10-year oscillation cycle is reinforced by the socio-economic predictor net migration. Considering the net migration and the Scandinavian index increases the explained variance of the GLM to 44%. The high-frequency oscillations (2-3 years) are associated with fluctuations of the natural TBE transmission cycle between small mammals and ticks, which are driven by beech fructification. Considering also fructification 2 years prior explains 64% of the variance of the TBE time series. Additionally, annual sunshine duration as predictor for the human outdoor activity increases the explained variance to 70%.

CONCLUSIONS

The GLMs presented here provide the basis for annual TBE forecasts, which were mainly determined by beech fructification. A total of 3 of the 5 years with full fructification, resulting in high TBE case numbers 2 years later, occurred after 2010. The effects of climate change are therefore not visible through a direct correlation of the TBE cases with rising temperatures, but indirectly via the increased frequency of mast seeding.

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.

Report of the 20th annual meeting of the International Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE): ISW-TBE: 20 years of commitment and still challenges ahead

The 20th Meeting of the International Scientific Working Group on Tick-Borne Encephalitis (ISW-TBE) - a group of neurologists, general practitioners, clinicians, travel physicians, virologists, pediatricians, ecologists, and epidemiologists - was held under the theme "ISW-TBE: 20 years of commitment & still challenges ahead". On the occasion of the 20^th anniversary of the ISW-TBE several key topics in TBE were extensively discussed, among them current epidemiological developments and investigations, expansion of risk areas, clinical aspects and cases, traveling and mobility, and latest news on TBE vaccination. The main goals that have been achieved by the ISW-TBE so far are, among others, an increased awareness in endemic and non-endemic countries, an increase of vaccination rates in various countries, getting TBE acknowledged and established as a travel-related risk, and building contact with the European Centre for Disease Prevention and Control.

Forecasting next season's Ixodes ricinus nymphal density: the example of southern Germany 2018

The castor bean tick, Ixodes ricinus (L.) (Ixodida: Ixodidae), is the principal vector of pathogens causing tick-borne encephalitis or Lyme borreliosis in Europe. It is therefore of general interest to make an estimate of the density of I. ricinus for the whole year at the beginning of the tick season. There are two necessary conditions for making a successful prediction: a long homogeneous time series of observed tick density and a clear biological relationship between environmental predictors and tick density. A 9-year time series covering the period 2009-2017 of nymphal I. ricinus flagged at monthly intervals in southern Germany has been used. With the hypothesis that I. ricinus density is triggered by the fructification of the European beech 2 years before, the mean annual temperature of the previous year, and the current mean winter temperature (December-February), a forecast of the annual nymphal tick density has been made. Therefore, a Poisson regression model was generated resulting in an explained variance of 93.4% and an error of [Formula: see text] ticks per [Formula: see text] (annual [Formula: see text] collected ticks/[Formula: see text]). An independent verification of the forecast for the year 2017 resulted in 187 predicted versus 180 observed nymphs per [Formula: see text]. For the year 2018 a relatively high number of 443 questing I. ricinus nymphs per [Formula: see text] is forecasted, i.e., a "good" tick year.