Perspectives on modelling the distribution of ticks for large areas: so far so good?

Stochastic modelling to predict a priori the potential outcomes of different control schemes against the cattle tick Rhipicephalus microplus

Stochastic models are valuable tools to describe and analyze tick population dynamic in a given area, and to evaluate different control schemes. The objective of this study was to provide a tool to predict a priori how a control scheme could affect the abundance of Rhipicephalus microplus in an area highly favourable for its development, through the building of a stochastic model of the population dynamic of this tick. The dynamic was stochastically modelled using field data of the parasitic and non-parasitic phases of R. microplus. The host susceptibility to tick infestations was also considered by including representative values of natural resistance of three different breeds (Bos indicus, B. taurus and B. indicus x B. taurus). Two different control schemes (strategic and threshold methods) using three annual applications of synthetic acaricides in different moments were evaluated. Furthermore, we have evaluated the impact of different moments of pasture spelling as a tool for tick control on the dynamic of the non-parasitic phases of R. microplus. The results showed that the model appropriately fits to field data and can simulate the annual variability in female ticks count in animals from different ecological characteristics. The model has also captured the variability that host susceptibility has on tick abundance. This fact was expressed in the different levels of tick abundance obtained for the three breeds of bovines. According to the model, strategic control applications of chemical acaricides have more efficacy than the threshold control method. The model indicates that a pasture spelling started in early or mid-summer generates a higher reduction in the number of available eggs in pastures. Analysis of the results allows to conclude that the model developed in this study was useful to predict a priori the performance of different control methods based on the strategic application of chemical acaricides or by applying pasture spelling.

Mapping the potential distribution of the principal vector of Crimean-Congo haemorrhagic fever virus Hyalomma marginatum in the Old World

Crimean-Congo haemorrhagic fever (CCHF) is the most widely distributed tick-borne viral disease in humans and is caused by the Crimean-Congo haemorrhagic fever virus (CCHFV). The virus has a broader distribution, expanding from western China and South Asia to the Middle East, southeast Europe, and Africa. The historical known distribution of the CCHFV vector Hyalomma marginatum in Europe includes most of the Mediterranean and the Balkan countries, Ukraine, and southern Russia. Further expansion of its potential distribution may have occurred in and out of the Mediterranean region. This study updated the distributional map of the principal vector of CCHFV, H. marginatum, in the Old World using an ecological niche modeling approach based on occurrence records from the Global Biodiversity Information Facility (GBIF) and a set of covariates. The model predicted higher suitability of H. marginatum occurrences in diverse regions of Africa and Asia. Furthermore, the model estimated the environmental suitability of H. marginatum across Europe. On a continental scale, the model anticipated a widespread potential distribution encompassing the southern, western, central, and eastern parts of Europe, reaching as far north as the southern regions of Scandinavian countries. The distribution of H. marginatum also covered countries across Central Europe where the species is not autochthonous. All models were statistically robust and performed better than random expectations (p < 0.001). Based on the model results, climatic conditions could hamper the successful overwintering of H. marginatum and their survival as adults in many regions of the Old World. Regular updates of the models are still required to continually assess the areas at risk using up-to-date occurrence and climatic data in present-day and future conditions.

Biogeography of the theileriosis vector, Rhipicephalus appendiculatus under current and future climate scenarios of Zimbabwe

Climate directly influences the epidemiology of vector-borne diseases at various spatial and temporal scales. Following the recent increased incidences of theileriosis in Zimbabwe, a disease mainly transmitted by Rhipicephalus appendiculatus, we determined lethal temperatures for the species and current and possible future distribution using the machine learning algorithm 'Maxent'. Rhipicephalus appendiculatus larvae had an upper lethal temperature (ULT₅₀) of about 44 ± 0.5 °C and this was marginally higher for nymphs and adults at 46 ± 0.5 °C. Environmental temperatures recorded in selected zonal tick microhabitats were below the determined lethal limits, indicating the ability of the tick to survive these regions. The resultant model under current climatic conditions showed areas with high suitability indices to the eastern, northeastern and southeastern parts of the country, mainly in Masvingo, Manicaland and Mashonaland Central provinces. Future predictions as determined by 2050 climatic conditions indicate a reduction in suitable habitats with the tick receding to presently cooler high elevation areas such as the eastern Highlands of Zimbabwe and a few isolated pockets in the interior of the country. Lowveld areas show low suitability under current climatic conditions and are expected to remain unsuitable in future. Overall, the study shows that R. appendiculatus distribution is constrained by climatic factors and helps identify areas of where occurrence of the species and the disease it transmits is highly likely. This will assist in optimizing disease surveillance and vector management strategies targeted at the species.

An Agenda for Research of Uncovered Epidemiological Patterns of Tick-Borne Pathogens Affecting Human Health

The panorama of ticks and tick-borne pathogens (TBP) is complex due to the many interactions among vertebrates, vectors, and habitats, occurring at different scales. At a broad spatial range, climate and host availability regulate most tick processes, including questing activity, development, and survival. At a local scale, interactions are obscured by a high indeterminacy, making it arduous to record in field surveys. A solid modelling framework could translate the local/regional empirical findings into larger scales, shedding light on the processes governing the circulation of TBP. In this opinion paper, we advocate for a re-formulation of some paradigms in the research of these outstanding cycles of transmission. We propose revisiting concepts that faced criticisms or lacked solid support, together with the development of a conceptual scheme exploring the circulation of TBP under a range of conditions. We encourage (i) an adequate interpretation of the niche concept of both ticks and vertebrate/reservoir hosts interpreting the (a)biotic components that shape the tick's niche, (ii) an assessment of the role played by the communities of wild vertebrates on the circulation of pathogens, and (iii) the development of new approaches, based on state-of-the-art epidemiological concepts, to integrate findings and modelling efforts on TBP over large regions.

Widespread exposure to Crimean-Congo haemorrhagic fever in Uganda might be driven by transmission from Rhipicephalus ticks: Evidence from cross-sectional and modelling studies

BACKGROUND

Crimean-Congo haemorrhagic fever (CCHF) is a widespread tick-borne viral infection, present across Africa and Eurasia, which might pose a cryptic public health problem in Uganda. We aimed to understand the magnitude and distribution of CCHF risk in humans, livestock and ticks across Uganda by synthesising epidemiological (cross-sectional) and ecological (modelling) studies.

METHODS

We conducted a cross-sectional study at three urban abattoirs receiving cattle from across Uganda. We sampled humans (n = 478), livestock (n = 419) and ticks (n = 1065) and used commercially-available kits to detect human and livestock CCHF virus (CCHFV) antibodies and antigen in tick pools. We developed boosted regression tree models to evaluate the correlates and geographical distribution of expected tick and wildlife hosts, and of human CCHF exposures, drawing on continent-wide data.

FINDINGS

The cross-sectional study found CCHFV IgG/IgM seroprevalence in humans of 10·3% (7·8-13·3), with antibody detection positively associated with reported history of tick bite (age-adjusted odds ratio = 2·09 (1·09-3·98)). Cattle had a seroprevalence of 69·7% (65·1-73·4). Only one Hyalomma tick (CCHFV-negative) was found. However, CCHFV antigen was detected in Rhipicephalus (5·9% of 304 pools) and Amblyomma (2·9% of 34 pools) species. Modelling predicted high human CCHF risk across much of Uganda, low environmental suitability for Hyalomma, and high suitability for Rhipicephalus and Amblyomma.

INTERPRETATION

Our epidemiological and ecological studies provide complementary evidence that CCHF exposure risk is widespread across Uganda. We challenge the idea that Hyalomma ticks are consistently the principal reservoir and vector for CCHFV, and postulate that Rhipicephalus might be important for CCHFV transmission in Uganda, due to high frequency of infected ticks and predicted environmental suitability.

FUNDING

UCL Global Challenges Research Fund (GCRF) and Pan-African Network on Emerging and Re-Emerging Infections (PANDORA-ID-NET) funded by the European and Developing Countries Clinical Trials Partnership (EDCTP) under the EU Horizon 2020 Framework Programme for Research and Innovation.

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.

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.

Hard tick species (Acari: Ixodidae) and infestation in two livestock agroecosystems from Antioquia, Colombia

Tick infestation affects about 80% of livestock globally while transmitting various pathogens causing high economic losses. This study aimed to determine the degree of tick infestation in two regions, North and Middle Magdalena in Antioquia, Colombia, to identify the ixodid tick species found and the associated risk factors. A cross-sectional study was carried out in 48 farms distributed in six municipalities of Antioquia. Two paddocks and eight bovines per farm were evaluated to estimate tick infestation (adults, nymphs, and larvae). Tick species were identified through a morphological and molecular analysis based on partial sequences of data obtained from DNA molecular markers, two mitochondrial (16S rRNA and COI), and one genomic DNA gene (18S rRNA). A multivariate Poisson regression model was applied to estimate the associated risk factors with ticks in cattle. Rhipicephalus microplus, Amblyomma patinoi and Dermacentor nitens were present in the livestock agroecosystems in the Middle Magdalena region; the highest incidence of tick infestation in cows and paddocks was reported in the municipality of Puerto Triunfo. The livestock agroecosystems in Middle Magdalena were characterized by a higher presence of adult R. microplus in cattle. Larval infestation of R. microplus, followed by D. nitens, was also found in paddocks. The multivariate analysis showed that the origin of cattle was the main risk factor associated with the presence of ticks (i.e., when cattle came from outside the farm). Cattle movement between farms in Middle Magdalena can contribute to the spread of ticks in this region.

Use of meteorological data in biosecurity

Pests, pathogens and diseases cause some of the most widespread and damaging impacts worldwide — threatening lives and leading to severe disruption to economic, environmental and social systems. The overarching goal of biosecurity is to protect the health and security of plants and animals (including humans) and the wider environment from these threats. As nearly all living organisms and biological systems are sensitive to weather and climate, meteorological, ‘met’, data are used extensively in biosecurity. Typical applications include, (i) bioclimatic modelling to understand and predict organism distributions and responses, (ii) risk assessment to estimate the probability of events and horizon scan for future potential risks, and (iii) early warning systems to support outbreak management. Given the vast array of available met data types and sources, selecting which data is most effective for each of these applications can be challenging. Here we provide an overview of the different types of met data available and highlight their use in a wide range of biosecurity studies and applications. We argue that there are many synergies between meteorology and biosecurity, and these provide opportunities for more widespread integration and collaboration across the disciplines. To help communicate typical uses of meteorological data in biosecurity to a wide audience we have designed the ‘Meteorology for biosecurity’ infographic.

Models for Studying the Distribution of Ticks and Tick-Borne Diseases in Animals: A Systematic Review and a Meta-Analysis with a Focus on Africa

Ticks and tick-borne diseases (TTBD) are constraints to the development of livestock and induce potential human health problems. The worldwide distribution of ticks is not homogenous. Some places are ecologically suitable for ticks but they are not introduced in these areas yet. The absence or low density of hosts is a factor affecting the dissemination of the parasite. To understand the process of introduction and spread of TTBD in different areas, and forecast their presence, scientists developed different models (e.g., predictive models and explicative models). This study aimed to identify models developed by researchers to analyze the TTBD distribution and to assess the performance of these various models with a meta-analysis. A literature search was implemented with PRISMA protocol in two online databases (Scopus and PubMed). The selected articles were classified according to country, type of models and the objective of the modeling. Sensitivity, specificity and accuracy available data of these models were used to evaluate their performance using a meta-analysis. One hundred studies were identified in which seven tick genera were modeled, with Ixodes the most frequently modeled. Additionally, 13 genera of tick-borne pathogens were also modeled, with Borrelia the most frequently modeled. Twenty-three different models were identified and the most frequently used are the generalized linear model representing 26.67% and the maximum entropy model representing 24.17%. A focus on TTBD modeling in Africa showed that, respectively, genus Rhipicephalus and Theileria parva were the most modeled. A meta-analysis on the quality of 20 models revealed that maximum entropy, linear discriminant analysis, and the ecological niche factor analysis models had, respectively, the highest sensitivity, specificity, and area under the curve effect size among all the selected models. Modeling TTBD is highly relevant for predicting their distribution and preventing their adverse effect on animal and human health and the economy. Related results of such analyses are useful to build prevention and/or control programs by veterinary and public health authorities.

Effects of white-tailed deer habitat use preferences on southern cattle fever tick eradication: simulating impact on "pasture vacation" strategies

BACKGROUND

Rhipicephalus (Boophilus) annulatus and Rhipicephalus (Boophilus) microplus (southern cattle fever tick; SCFT), collectively known as cattle-fever ticks (CFTs), are vectors of protozoal parasites (Babesia bigemina and Babesia bovis) that cause bovine babesiosis (also known as cattle fever). One traditional strategy for CFT eradication involves the implementation of a "pasture vacation," which involves removing cattle (Bos taurus) from an infested pasture for an extended period of time. However, vacated pastures are often inhabited by wildlife hosts, such as white-tailed deer (WTD; Odocoileus virginianus), which can serve as alternate hosts for questing CFTs. We hypothesized that the distribution of host-seeking larvae among habitat types post-pasture vacation would reflect habitat use patterns of WTD, and in turn, affect the subsequent rate of pasture infestation by CFT.

METHODS

We adapted a spatially explicit, individual-based model to simulate interactions among SCFT, cattle, and WTD as a tool to investigate the potential effects of WTD habitat use preferences on the efficacy of a pasture vacation. We parameterized the model to represent conditions typical of rangelands in south Texas, USA, simulated a 1-year pasture vacation under different assumptions regarding WTD habitat use preferences, and summarized effects on efficacy through (1) time post-vacation to reach 100% of pre-vacation densities of host-seeking larvae, and (2) the ecological conditions that resulted in the lowest host-seeking larval densities following pasture vacation.

RESULTS

Larval densities at the landscape scale varied seasonally in a similar manner over the entire simulation period, regardless of WTD habitat use preferences. Following the removal of cattle, larval densities declined sharply to < 100 larvae/ha. Following the return of cattle, larval densities increased to > 60% of pre-vacation densities ≈ 21 weeks post-vacation, and reached pre-vacation levels in less than a year. Trends in larval densities in different habitat types paralleled those at the landscape scale over the entire simulation period, but differed quantitatively from one another during the pasture vacation. Relative larval densities (highest to lowest) shifted from (1) wood/shrub, (2) grass, (3) mixed-brush during the pre-vacation period to (1) mixed-brush, (2) wood/shrub, (3) grass or (1) wood/shrub, (2) mixed-brush, (3) grass during the post-vacation period, depending on WTD habitat use preferences.

CONCLUSIONS

By monitoring WTD-driven shifts in distributions of SCFT host-seeking larvae among habitat types during simulated pasture vacation experiments, we were able to identify potential SCFT refugia from which recrudescence of infestations could originate. Such information could inform timely applications of acaricides to specific refugia habitats immediately prior to the termination of pasture vacations.

Modelling the potential spread of Hyalomma marginatum ticks in Europe by migratory birds

This study modelled the probability of introduction of Hyalomma marginatum into Europe by predicting the potential migratory routes of 28 bird species and the probability to carry immatures of the tick. Flyways were modelled as a spatio-temporal feature, at weekly intervals, using satellite-derived data of temperature and vegetal phenology, together with cost surfaces derived from speed and direction of the wind (years 2002-2018). The expected period of activity of tick immatures defined the probability of ticks being carried by birds along the modelled flyways. The probability of moulting of the engorged nymphs was modelled as a linear relationship of the daily sum of temperatures after tick introduction by birds. Positive probabilities of tick introduction extend the known northern range of permanent populations to central and western France, and large portions of central Europe. The flight of birds into an area and thence the risk of introduction of H. marginatum is very heterogeneous, with sites receiving "waves" of different bird species at diverse times of the year. Therefore, there is not a clear period of time for introduction, as it depends on the modelled behaviour of the bird species. The probability of introduction into Baltic and Nordic countries is small. We hypothesise that conditions of a warmer climate might support permanent populations of H. marginatum if a high number of immatures is introduced. Active surveys in risky territories, where the tick is not yet established, are advisable for rapid intervention.

MaxEnt Modeling of Dermacentor marginatus (Acari: Ixodidae) Distribution in Xinjiang, China

Dermacentor marginatus Sulkzer is a common tick species found in the Xinjiang Uygur Autonomous Region (XUAR) of China, and is a vector for a variety of pathogens. To determine the potential distribution of this tick species in Xinjiang, a metadata containing 84 D. marginatus presence records combined with four localities from field collection were used for MaxEnt modeling to predict potential distribution of this tick species. Identification of tick samples showed 756 of 988 (76%) were D. marginatus. MaxEnt modeling results indicated that the potential distribution of this tick species was mainly confined to northern XUAR. Highly suitable areas included west side of Altay mountain, west rim of Junggar basin, and Yili River valley in the study area. The model showed an AUC value of 0.838 ± 0.063 (SD), based on 10-fold cross-validation. Although tick presence records used for modeling were limited, this is the first regional tick distribution model for D. marginatus in Xinjiang. The model will be helpful in assessing the risk of tick-borne diseases to human and animals in the region.

Knowledge and perception on ticks and tick-borne diseases among veterinary medicine students from the North African countries of Algeria, Egypt, and Tunisia

Ticks are important vectors of both animal and human pathogens. The epidemiology of tick-borne diseases (TBDs) has dramatically changed in several regions in the world. As parasitology is a continuously growing field, assessing the knowledge of veterinary medicine students provides useful indicators and information on the level of intervention required to adapt parasitological courses to meet the demands in a changing world.

This study aimed to assess, in three North African veterinary education establishments, the basic parasitology knowledge of veterinary medicine students. Such a study is essential to build up core competencies regarding ticks and TBDs, and to suggest suitable adjustment measures to parasitology courses.

The present study was based on a self-administered and anonymous questionnaire on ticks and TBDs basic knowledge and perception. The survey was completed by 558 veterinary medicine students in Algeria, Egypt, and Tunisia in 2018. The students were divided into two groups: the “before” group – students who had not yet completed the parasitology course, and the “after” group – students who had already completed it.

In all studied countries, the “after” students' group gave significantly more correct answers (83.16%) than the “before” students' group (16.84%) (p < 0.001). Similarly, the percentage of “no answer” was higher in the “before” students' group (51.02%) compared to the “after” students' group (48.98%) (p < 0.001). The most frequent false answers provided by the “after” students' group regarded the number of tick species present in their own countries (5.14% of correct answers), and the most common tick species in their countries (18.11% of correct answers).

Almost 58.38% (216/370) of the “after” students' group knew that ticks transmit zoonotic pathogens to humans; among them, only 63 (29.17%) gave the correct names of the zoonotic diseases in their country.

Among the three countries, more than 80% of the “after” students' group thought that climate has an influences on ticks. According to this group, the most frequent factor that has influences on ticks' abundance is heat (53.02%).

As North African countries share several similitudes, we suggest creating a network of parasitological teachers where common teaching sources and resources could be developed for both teachers and students in the region. This network would allow the exchange of teaching approaches and materials to introduce harmonization into veterinary parasitological courses across North African countries. This is particularly important when considering the increasing incidence of ticks and TBDs in the region.

A Retrospective Assessment of Temperature Trends in Northern Europe Reveals a Deep Impact on the Life Cycle of Ixodes ricinus (Acari: Ixodidae)

This study modelled the changes in the development processes of the health-threatening tick Ixodes ricinus in Northern Europe as driven by the trends of temperature (1950–2018). We used the ECA&D dataset to calculate the annual accumulated temperature to obtain the development rates of the oviposition, incubation, larva–nymph, and nymph–adult molts. Annual values were used to ascertain the trend in development rates of each stage. The ecological classification of Northern Europe (LANMAP2) was used to summarize results. The temperature in 1950–2018 clearly increased in the target territory. The development rates of every tested life cycle process were faster along the time series. Faster oviposition and incubation rates resulted in central Sweden, Baltic countries, and parts of Finland. Faster molting rates were observed in the same territories and in large areas of Western Norway. The trend of temperature in the period 1950–2018 shows a consistent inflection point around 1990, demonstrating that the increased annual accumulated temperature has a deeper impact on the life cycle of I. ricinus since approximately 1990. Faster development rates could be part of the processes driving the reported spread of the tick in the target area and should be considered as a serious threat to human health.

Mapping the environmental suitability of etiological agent and tick vectors of Crimean-Congo hemorrhagic fever

Crimean-Congo hemorrhagic fever (CCHF) is one of the most important public health threats in many regions across Africa, Europe, and Asia. This study used ecological niche modeling analyses to map the environmental suitability of both CCHF virus (CCHFV), and its tick vectors (Amblyomma variegatum, Dermacentor marginatus, Hyalomma marginatum, Hyalomma rufipes, Hyalomma truncatum, Rhipicephalus appendiculatus, and Rhipicephalus evertsi evertsi) in the Old World countries. The CCHFV was anticipated to occur with high environmental suitability across southern and central Europe, northwestern Africa, central Asia, and western Mediterranean region. Ecological niche models of tick vectors anticipated diverse patterns based on the tick species in question; D. marginatus and H. marginatum showed high environmental suitability in southern and central Europe, and North Africa. The remaining vector species were anticipated to occur in Africa. All models were statistically robust and performed better than random (P < 0.001). Finally, we tested the niche similarities between CCHFV and diverse tick vectors and could not reject the null hypotheses of niche similarity in all vector-virus combinations (P > 0.05) except the combinations of CCHFV with A. variegatum, R. evertsi evertsi and R. appendiculatus (P < 0.05).

An updated meta-analysis of the distribution and prevalence of Borrelia burgdorferi s.l. in ticks in Europe

Background

The bacteria of the group Borrelia burgdorferi s.l. are the etiological agents of Lyme borreliosis in humans, transmitted by bites of ticks. Improvement of control measures requires a solid framework of the environmental traits driving its prevalence in ticks.

Methods

We updated a previous meta-analysis of the reported prevalence of Borrelia burgdorferi s.l. in questing nymphs of Ixodes ricinus with a literature search from January 2010–June 2017. This resulted in 195 new papers providing the prevalence of Bb for 926 geo-referenced records. Previously obtained data (878 records, years 2000–2010) were appended for modelling. The complete dataset contains data from 82,004 questing nymphs, resulting in 558 records of B. afzelii, 404 of B. burgdorferi s.s. (only 80 after the year 2010), 552 of B. garinii, 78 of B. lusitaniae, 61 of B. spielmanii, and 373 of B. valaisiana. We associated the records with explicit coordinates to environmental conditions and to a categorical definition of European landscapes (LANMAP2) looking for a precise definition of the environmental niche of the most reported species of the pathogen, using models based on different classification methods.

Results

The most commonly reported species are B. afzelii, B. garinii and B. valaisiana largely overlapping across Europe. Prevalence in ticks is associated with portions of the environmental niche. Highest prevalence occurs in areas of 280°–290° (Kelvin) of mean annual temperature experiencing a small amplitude, steady spring slope, together with high mean values and a moderate spring rise of vegetation vigor. Low prevalence occurs in sites with low and a noteworthy annual amplitude of temperature and the Normalized Difference Vegetation Index (colder areas with abrupt annual changes of vegetation). Models based on support vector machines provided a correct classification rate of the habitat and prevalence of 89.5%. These results confirm the association of prevalence of the three most commonly reported species of B. burgdorferi s.l. in Europe to parts of the environmental niche and provide a statistically tractable framework for analyzing trends under scenarios of climate change.

Electronic supplementary material

The online version of this article (10.1186/s12942-018-0163-7) contains supplementary material, which is available to authorized users.

Effect of Climate and Land Use on the Spatio-Temporal Variability of Tick-Borne Bacteria in Europe

The incidence of tick-borne diseases caused by Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum and Rickettsia spp. has been rising in Europe in recent decades. Early pre-assessment of acarological hazard still represents a complex challenge. The aim of this study was to model Ixodes ricinus questing nymph density and its infection rate with B. burgdorferi s.l., A. phagocytophilum and Rickettsia spp. in five European countries (Italy, Germany, Czech Republic, Slovakia, Hungary) in various land cover types differing in use and anthropisation (agricultural, urban and natural) with climatic and environmental factors (Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Land Surface Temperature (LST) and precipitation). We show that the relative abundance of questing nymphs was significantly associated with climatic conditions, such as higher values of NDVI recorded in the sampling period, while no differences were observed among land use categories. However, the density of infected nymphs (DIN) also depended on the pathogen considered and land use. These results contribute to a better understanding of the variation in acarological hazard for Ixodes ricinus transmitted pathogens in Central Europe and provide the basis for more focused ecological studies aimed at assessing the effect of land use in different sites on tick–host pathogens interaction.

Climate change influences on the potential geographic distribution of the disease vector tick Ixodes ricinus

Background

Ixodes ricinus is a species of hard tick that transmits several important diseases in Europe and North Africa, including Lyme borreliosis and tick-borne encephalitis. Climate change is affecting the geographic distributions and abundances of arthropod vectors, which in turn influence the geographic distribution and epidemiology of associated vector-borne diseases. To date, few studies have investigated effects of climate change on the spatial distribution of I. ricinus at continental extents. Here, we assessed the potential distribution of I. ricinus under current and future climate conditions to understand how climate change will influence the geographic distribution of this important tick vector in coming decades.

Method

We used ecological niche modeling to estimate the geographic distribution of I. ricinus with respect to current climate, and then assessed its future potential distribution under different climate change scenarios. This approach integrates occurrence records of I. ricinus with six relevant environmental variables over a continental extent that includes Europe, North Africa, and the Middle East. Future projections were based on climate data from 17 general circulation models (GCMs) under 2 representative concentration pathway emissions scenarios (RCPs), for the years 2050 and 2070.

Result

The present and future potential distributions of I. ricinus showed broad overlap across most of western and central Europe, and in more narrow zones in eastern and northern Europe, and North Africa. Potential expansions were observed in northern and eastern Europe. These results indicate that I. ricinus populations could emerge in areas in which they are currently lacking, posing increased risks to human health in those areas. However, the future of I. ricinus ticks in some important regions such the Mediterranean was unclear owing to high uncertainty in model predictions.

Host Distribution Does Not Limit the Range of the Tick Ixodes ricinus but Impacts the Circulation of Transmitted Pathogens

Ticks, pathogens, and vertebrates interact in a background of environmental features that regulate the densities of ticks and vertebrates, affecting their contact rates and thence the circulation of the pathogens. Regional scale studies are invaluable sources of information about the regulation of these interactions, but a large-scale analysis of the interaction of communities of ticks, hosts, and the environment has been never modeled. This study builds on network analysis, satellite-derived climate and vegetation, and environmental modeling, quantifying the interactions between the tick Ixodes ricinus and the transmitted bacteria of the complex Borrelia burgdorferi s.l. in the Western Palaearctic. We derived the rates of contact of the tick with 162 species of vertebrates recorded as hosts, and the relative importance of each vertebrate in the circulation of the pathogen. We compiled more than 11 millions of pairs of coordinates of the vertebrates, deriving distribution models of each species and the relative faunal composition in the target territory. The results of the modeling of the distribution of the tick and its hosts, weighted by their importance in the circulation of Borrelia captured the spatial patterns of interactions that allow the circulation of the pathogen. Results indicate that both I. ricinus and B. burgdorferi s.l. are supported in the Western Palaearctic by complex communities of vertebrates, which have large distribution ranges. This high functional redundancy results in the pervasiveness of B. burgdorferi s.l., which depends on the gradient of contributions of the large community of vertebrates, instead of relying on a few dominant vertebrates, which was the prevailing paradigm. Most prominent reservoirs of the pathogen are distributed in specific regions of the environmental niche. However, literally dozens of potential reservoirs can colonize many other environmental regions, marginally but efficiently contributing to the circulation of the pathogen. These results consistently point to the need of evaluating the beta-diversity of the community of vertebrates acting as reservoirs of the pathogen to better know the interactions with the vector. They also demonstrate why the pathogen is so resilient to perturbations in the composition of the reservoirs.