Skip to main content

Effects of global changes on the climatic niche of the tick Ixodes ricinus inferred by species distribution modelling

Abstract

Background

Global climate change can seriously impact on the epidemiological dynamics of vector-borne diseases. In this study we investigated how future climatic changes could affect the climatic niche of Ixodes ricinus (Acari, Ixodida), among the most important vectors of pathogens of medical and veterinary concern in Europe.

Methods

Species Distribution Modelling (SDM) was used to reconstruct the climatic niche of I. ricinus, and to project it into the future conditions for 2050 and 2080, under two scenarios: a continuous human demographic growth and a severe increase of gas emissions (scenario A2), and a scenario that proposes lower human demographic growth than A2, and a more sustainable gas emissions (scenario B2). Models were reconstructed using the algorithm of “maximum entropy”, as implemented in the software Maxent 3.3.3e; 4,544 occurrence points and 15 bioclimatic variables were used.

Results

In both scenarios an increase of climatic niche of about two times greater than the current area was predicted as well as a higher climatic suitability under the scenario B2 than A2. Such an increase occurred both in a latitudinal and longitudinal way, including northern Eurasian regions (e.g. Sweden and Russia), that were previously unsuitable for the species.

Conclusions

Our models are congruent with the predictions of range expansion already observed in I. ricinus at a regional scale and provide a qualitative and quantitative assessment of the future climatically suitable areas for I. ricinus at a continental scale. Although the use of SDM at a higher resolution should be integrated by a more refined analysis of further abiotic and biotic data, the results presented here suggest that under future climatic scenarios most of the current distribution area of I. ricinus could remain suitable and significantly increase at a continental geographic scale. Therefore disease outbreaks of pathogens transmitted by this tick species could emerge in previous non-endemic geographic areas. Further studies will implement and refine present data toward a better understanding of the risk represented by I. ricinus to human health.

Background

Global climate change can seriously impact on the epidemiological dynamics of vector-borne diseases [15]. However, because of the influence of several biotic and abiotic factors on hosts, arthropods and pathogens they vector, future spatial and temporal distribution of vector-borne diseases is still difficult to predict [4, 6]. Under the above circumstances, predicting how climatic changes will affect geographic distribution of the vector is an obligate step. Indeed, vectors play a key role in infectious disease areas where they may transmit pathogens to a variety of animal hosts, often representing the bridge between zoonotic reservoirs and humans [3, 7, 8]. Likewise, since pathogens may disperse through arthropods into previously non-endemic areas, climate-induced changes in vector distribution ultimately affect the epidemiology of vector-borne diseases [3, 9]. In addition, a distributional shift of vectors may also lead to spatial overlap of different vector species [10], thus changing the impact of pathogen transmission, interfering with their epidemic cycles [4] or as a consequence of vector interspecific competition [1113].

In the last few years, Species Distribution Modelling (SDM) has greatly contributed to understanding the effect of climatic changes on vector distribution, combining known occurrence points of a species with a set of climatic variables [14, 15]. Indeed, such an approach allows developing the potential geographic distribution under current and future climatic conditions [1619]. SDM has been used in several studies on arthropod vectors [2025], allowing to infer future range expansion, such as in the case of sand fly species Lutzomyia anthophora and L. diabolica in North America [22], or a potential distributional shift and/or a reduction of geographic range, as suggested for the malaria vectors in Africa [21].

In this study we investigated how future climatic changes could affect the climatic niche of the tick Ixodes ricinus (Acari, Ixodidae). This species is regarded as the most important vector occurring in Eurasian regions since it is a multi-competent vector of bacteria (i.e. Borrelia burgdorferi s.l., Babesia spp., Anaplasma and Erlichia spp.) and viruses (i.e. Flavivirus spp.), to humans and animals [2630]. The above tick species is regarded as a major vector responsible for many zoonotic diseases, such as Lyme disease, babesiosis, anaplasmosis, erlichiosis and Tick-Borne Encephalitis (TBE) [31]. Ixodes ricinus has been recorded in Europe, Russia, up to the Caspian Sea on the east, and North Africa [32]. Along the western boundary of Russia and the neighboring countries, its range overlaps with the range of I. persulcatus tick, distributed in Eastern Europe and across Asia [30, 33, 34]. Recently, a latitudinal and altitudinal shift has been reported in I. ricinus distribution in European regions, and temperature rise was suggested to be among the factors responsible for this phenomenon [3540]. In spite of the epidemiological implications that vector distribution changes might have on vectored pathogens, no studies have investigated the impact of future climatic changes on the geographic distribution of this tick species at a continental geographic scale.

By using SDM we aimed to reconstruct the climatic niche of I. ricinus and projected it into the future conditions for 2050 and 2080, under two possible scenarios: i) a continuous human demographic growth and a severe increase of gas emissions (A2 scenario), ii) a lower human demographic growth than A2, and a more sustainable gas emissions (scenario B2) [41].

Methods

To reconstruct the current and future climatic niche of I. ricinus, the algorithm of “maximum entropy” has been used, as implemented in the software Maxent 3.3.3e. This technique, using presence-only points in conjunction with environmental variables, estimates the potential distribution of the species finding the probability distribution that approximates the distribution of maximum entropy [42]. This approach is largely used to reconstruct SDM because its performances are highly competitive with the other modeling methods. Indeed it showed the better calibration when compared to 16 other algorithms, including several traditional tools using presence-absence data, such as general linear models (GLM) and general additive models (GAM) [18, 4346]. To model the current and future climatic niche of I. ricinus, we used 4,544 occurrence points (Figure 1) obtained from the digital dataset available at https://sites.google.com/site/palticks/home/download[47]. The dataset includes occurrence points obtained by a systematic search of the published, peer-reviewed literature since, approximately, 1970 to 2010; records of ticks available in some curated collections were also included. It covers most of the I. ricinus range and does not include records considered as potentially incorrect, based on the known distribution of I. ricinus[47]. Environmental data were downloaded from WorldClim database (http://www.worldclim.org) with a resolution of 2.5 arc-minutes (~ 5 km). Nineteen bioclimatic variables derived from monthly temperature and rainfall values were downloaded, that represent annual trends, seasonality and extreme or limiting environmental factors. The runs have been made using 75% of the occurrence points to construct the model (training data) and the remaining 25% to test it. The default parameters of Maxent have been used with the exception of the Regularization parameter β. This parameter acts as a multiplier for the default values and it regulates the smoothness and regularity of the model [44, 48]. Recently, Warren & Seifert [48] have highlighted the importance to test different values of this parameter to improve the Maxent’s performance. We developed ten models using different values of β (1, 3, 5, 7, 9, 11, 13, 15, 17 and 19) and then we chose the model that outperformed all the others using the sample sized corrected Akaike Information Criterion (AICc) score, as implemented in the software ENMTools. The model built with the value of β equal to 1 gained the lowest score for the AICc, thus outperforming all the remaining nine models tested, therefore it has been used to reconstruct the final models. To evaluate the accuracy of the developed models we used the area under the curve (AUC) of the receiver operating characteristic (ROC) as suggested by Peterson et al.[49]. Explorative analyses have been done using all 19 variables, then the environmental variables outside the range present in the training data have been excluded. The final models have been reconstructed using 15 variables (Additional file 1). Ten replicates were run using the cross-validation form of replication. This approach randomly split the data into equal-size groups (“folds”) and creates models leaving out each fold in turn and using them for evaluation [50]. The model developed for the present-day conditions has been projected onto the future climatic conditions. Two scenarios proposed by the Special Report of Emission Scenarios (SRES) of the Inter-governmental Panel on Climate Change (IPCC) [51], have been considered, namely: scenario A2 (that proposes a continuous human demographic growth and a severe increase of gas emissions) and scenario B2 (that proposes lower human demographic growth than A2, and a more sustainable gas emissions) [17, 41]. Potential climatic niche under these scenarios has been predicted for 2050 and 2080 using the Canadian Center for Climate Modeling and Analysis CCCAM-CGCM3.1-T47 model [52].

Figure 1
figure 1

Occurrence points. Map shows the 4,544 records of I. ricinus used to develop the Species Distribution Modelling.

All SDM predictions were visualized in Quantum-GIS 1.8.0 (http://download.qgis.org). To estimate the future increase of the niche area respect to current model, presence/absence maps using the minimum training presence threshold were constructed [53]. Then future increase of the niche area respect to current model was calculated using the software ImageJ (http://rsb.info.nih.gov/ij/).

Results

The averaged climatic niche of the species for current conditions is shown in Figure 2a. The averaged value of AUC for this model was 0.860 (± 0.004), indicating an optimal performance of the models. Among the 15 variables used BIO6 (Min Temperature in the coldest period) and BIO17 (Precipitation of the driest quarter) give the highest percent contribution to construct the model (Additional file 1), according to the biology of I. ricinus, that lives in biotopes that offer moderate temperatures and high relative humidity [54]. The climatically suitable area predicted by our model under current conditions encompasses the known geographic distribution of I. ricinus and that previously inferred for this species using both climatic features and vegetation index [55], supporting the validity of our reconstruction.

Figure 2
figure 2

Geographic distribution predicted by species distribution modelling for Ixodes ricinus . (a) Current conditions; (b-c) 2050, scenario A2 and B2, respectively; (d-e) 2080, scenario A2 and B2, respectively.

In Figure 2b-e the models predicted under future climatic conditions are shown. In both scenarios an increase of a climatically suitable area of about two times greater than the current area was predicted (Additional file 2). Such an increase occurred both in a latitudinal and longitudinal way, including northern Eurasian regions (e.g. Sweden and Russia), that were previously unsuitable for the species. Notably, climatically suitable conditions appeared also in the easternmost Asian region, such as Central China, South Korea, Japan and Kamchatka Peninsula. While the two scenarios showed no difference about the extension of the predicted suitable area, they instead showed some difference in its climatic suitability. Indeed, under the scenario B2 most of the predicted area showed a higher climatic suitability than A2 in both 2050 and 2080. Furthermore, under the scenario B2 an increase of climatic suitability for I. ricinus has been predicted between 2050 and 2080 in several areas (i.e. Ireland and Eastern Europe) (Figure 2c,e).

Discussion

The increase in temperature, which has been predicted to occur in the future years [41], could drastically affect the ecology and geographic distribution of I. ricinus in Eurasian regions. Climatic changes have been shown to affect seasonal activity and feeding behavior of I. ricinus, at different life stages [3, 32, 40]. Furthermore, because of the I. ricnus vulnerability to drought [54], global warming may affect negatively the species in southernmost part of its range of distribution. On the other hand, higher temperature could lead to milder winter and extended spring and autumn seasons than those actually characterizing northern regions, making these areas climatically suitable for the species. Evidence that this latter phenomenon is in progress in Europe has been already recorded for I. ricinus, in which an expansion of the northern distribution limits in Sweden and Norway has been reported since the 1980s [35, 3840, 56]. Likewise, a shift of the limit to higher altitudes northward in Czech Republic and Switzerland respectively, has been shown [36, 37]. Our results by using SDM are congruent with the predictions of range expansion and provide a qualitative and quantitative assessment of the future areas climatically suitable for I. ricinus at a continental scale.

While interpreting inferred expansion of climatic niche in the context of possible changes in geographic distribution of I. ricinus[1518] the following two issues should be considered. First, the climatic features are only one of the ecological features that affect the geographic distribution of tick species [57, 58]. However, I. ricinus, as all other ixodid ticks, spend most of their life off the host, so that climate is an essential determinant of their occurrence [59, 60]. At a smaller geographic scale, other abiotic factors should be considered to assess the effective occurrence of I. ricinus, such as landscape physical features, or landscape use [15] and biotic factors (competition, hosts abundance etc.) [15, 61, 62]. The integration of these factors could show discontinuity in areas that our model predicted as large continuous areas of climate suitability. With respect to biotic factors, the hosts are key components for the ecological niche of ticks [7, 63, 64]. Ixodes ricinus is able to exploit a large variety of terrestrial vertebrates [7, 6365], so that host occurrence should not be a limiting factor to their persistence under climatic change scenarios [66]. Its wide ecological plasticity with respect to host choice, for example, was a key factor in allowing I. ricinus to survive during the last glacial phases without significant range reduction across the European continent [67].

Secondly, it should be considered that in order to have an effective range expansion, the new climatically suitable areas should be reached by the species. Host movements largely determine tick dispersal during infestation [65, 6872]. Among the several hosts exploited by I. ricinus, some of them, such as birds and cervids, are characterized by high dispersal ability. Host mediated dispersal of ticks also across long distances and more distant geographic areas have been demonstrated in several studies both for I. ricinus[7274] and for other tick species [40, 67, 7577]. Following the above considerations, although the use of SDM models at a higher resolution should be integrated by a more refined analysis of further abiotic and biotic data, the general significance of results presented here suggest that under future climatic scenarios, most of the current distribution area of I. ricinus could remain suitable and significantly increase at a continental geographic scale.

Our results suggest also that the more favourable climatic conditions for I. ricinus will occur under the scenario B2 (i.e., lower human demographic growth, and a more sustainable gas emissions than A2). Interestingly, reduced demographic growth and gas emissions are objectives of several international policies [78]. In this context it is further advisable to investigate the factors that could affect the epidemiological dynamics in which this tick is involved. For example, future inferred expansion eastward of the climatic niche of I. ricinus could increase the overlapping areas with the Eurasian tick I. persulcatus[79], increasing or decreasing pathogen transmission due to vector interspecific competition [1113], or to interferences among pathogen epidemic cycles [4]. This may be the case for Borrelia spp. in which a relationship has been established between bacterial complexes and tick species. Indeed, while the Eurasian type of B. garinii and B. afzelii are carried by both I. ricinus and I. persulcatus in North Asian regions, B. burgdorferi s.s. seems to be vectored exclusively by I. ricinus, and the Asian type of B. garinii by I. persulcatus[80]. The expansion of the areas of co-presence of the two tick species could enhance the diffusion of the types for which both of them are competent.

Conclusions

In the epidemiological systems involving I. ricinus as a vector, climatic changes have been shown to have several effects. First, the seasonal activity of ticks could undergo some changes [40, 81]. Indeed, previous studies conducted on I. ricinus showed that this species could prolong its questing season, usually spanning from March to November, until January as a response of milder winters due to the temperature increase [3, 8184]. As a consequence, more abundant populations of ticks could survive the winter, thus a higher probability of tick bites, and in turn of disease transmission, could be expected [82]. Second, different life stages could become active and search for a host simultaneously [32, 40]. Larvae and nymphs could parasitize the same host individual at the same time, favoring the trans-stadial transmission of pathogens by co-feeding and enhancing the efficiency of disease transmission [40, 85]. In addition, our results suggest a potential significant increase of climatic niche of I. ricinus in the future years under both scenarios. Therefore disease outbreaks of pathogens transmitted by this tick species could emerge in previous non-endemic geographic areas. Further studies will implement and refine present data toward a better understanding of the risk represented by I. ricinus for human health.

References

  1. Githeko AK, Lindsay SW, Confalonieri UE, Patz JA: Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ. 2000, 78: 1136-1147.

    PubMed Central  CAS  PubMed  Google Scholar 

  2. Patz JA, Olson SH: Climate change and health: global to local influences on disease risk. Ann Trop Med Parasit. 2006, 100: 535-549. 10.1179/136485906X97426.

    Article  CAS  PubMed  Google Scholar 

  3. Gage KL, Burkot TR, Eisen RJ, Hayes EB: Climate and vector borne disease. Am J Prev Med. 2008, 35: 436-450. 10.1016/j.amepre.2008.08.030.

    Article  PubMed  Google Scholar 

  4. Tabachnick WJ: Challenges in predicting climate and environmental effects on vector-borne disease episystems in a changing world. J Exp Biol. 2009, 213: 946-954.

    Article  Google Scholar 

  5. Mirski T, Bartoszcze M, Bielawska-Dròzd A: Impact of climate change on infectious diseases. Pol J Environ Stud. 2012, 3: 525-532.

    Google Scholar 

  6. Mills JN, Gage KL, Khan AS: Potential influence of climate change on vector-borne and zoonotic diseases: a review and proposed research plan. Environ Health Perspect. 2010, 118: 1507-1514. 10.1289/ehp.0901389.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Sonenshine DE: The biology of tick vectors of humans disease. Tick-borne diseases of humans. Edited by: Goodman JL, Dennis DT, Sonenshine DE. 2003, Washington: ASM Press, 12-31.

    Google Scholar 

  8. Eisen RJ, Eisen L: Spatial modeling of human risk of exposure to vector-borne pathogens based on epidemiological versus arthropod vector data. J Med Entomol. 2008, 45: 181-192. 10.1603/0022-2585(2008)45[181:SMOHRO]2.0.CO;2.

    Article  PubMed  Google Scholar 

  9. Sutherst RW: Global change and human vulnerability to vector-borne diseases. Clin Microbiol Rev. 2004, 17: 136-173. 10.1128/CMR.17.1.136-173.2004.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Brooks DR, Hoberg EP: How will global climate change affect parasite–host assemblages?. Trends Parasitol. 2007, 23: 571-574. 10.1016/j.pt.2007.08.016.

    Article  PubMed  Google Scholar 

  11. Reitz SR, Trumble JT: Competitive displacement among insects and arachnids. Annu Rev Entomol. 2002, 47: 435-465. 10.1146/annurev.ento.47.091201.145227.

    Article  CAS  PubMed  Google Scholar 

  12. Little E, Barrera R, Seto KC, Diuk-Wasser M: Co-occurrence patterns of the dengue vector aedes aegypti and aedes mediovitattus, a dengue competent mosquito in Puerto Rico. Ecohealth. 2011, 8: 365-375. 10.1007/s10393-011-0708-8.

    Article  PubMed  Google Scholar 

  13. Kweka EJ, Zhou G, Beilhe LB, Dixit A, Afrane Y, Gilbreath TM, Munga S, Nyindo M, Githeko AK, Yan G: Effects of co-habitation between anopheles gambiae s.s. And Culex quinquefasciatus aquatic stages on life history traits. Parasit Vectors. 2012, 5: 33-10.1186/1756-3305-5-33.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Guisan A, Thuiller W: Predicting species distribution: offering more than simple habitat models. Ecol Lett. 2005, 8: 993-1009. 10.1111/j.1461-0248.2005.00792.x.

    Article  Google Scholar 

  15. Zimmermann NE, Edwards TC, Graham CH, Pearman PB, Svenning JC: New trends in species distribution modelling. Ecography. 2010, 33: 985-989. 10.1111/j.1600-0587.2010.06953.x.

    Article  Google Scholar 

  16. Peterson AT: Ecologic niche modeling and spatial patterns of disease transmission. Emerg Infecti Dis. 2006, 12: 1822-1826. 10.3201/eid1212.060373.

    Article  Google Scholar 

  17. Quante M: The changing climate: past, present, future. Relict species: Phylogeography and Conservation biology. Edited by: Habel JC, Assmann T. 2010, Berlin: Springer, 9-56.

    Chapter  Google Scholar 

  18. Liu X, Guo Z, Ke Z, Wang S, Li Y: Increasing potential risk of a global aquatic invader in Europe in contrast to other continents under future climate change. PLoS One. 2011, 6: e18429-10.1371/journal.pone.0018429.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Caminade C, Medlock JM, Ducheyne E, McIntyre KM, Leach S, Baylis M, Morse AP: Suitability of European climate for the Asian tiger mosquito Aedes albopictus: recent trends and future scenarios. J R Soc Interface. 2013, 9: 2708-2717.

    Article  Google Scholar 

  20. Brownstein JS, Holford TR, Fish D: Effect of climate change on Lyme disease risk in North America. Eco Health. 2005, 2: 38-46.

    PubMed Central  PubMed  Google Scholar 

  21. Peterson AT: Shifting suitability for malaria vectors across Africa with warming climates. BMC Infect Dis. 2009, 9: 59-10.1186/1471-2334-9-59.

    Article  PubMed Central  PubMed  Google Scholar 

  22. Gonzàlez C, Wang O, Strutz SE, Gonzàlez-Salazar C, Sànchez-Cordero V, Sarkar S: Climate change and risk of leishmaniasis in North America: predictions from ecological niche models of vector and reservoir species. PLoS Negl Trop Dis. 2010, 4: e585-10.1371/journal.pntd.0000585.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Larson SR, De Groote JP, Bartholomay LC, Sugumaran R: Ecological niche modeling of potential West Nile virus vector mosquito species in Iowa. J Insect Sci. 2010, 10: 110-

    Article  PubMed Central  PubMed  Google Scholar 

  24. Rose H, Wall R: Modelling the impact of climate change on spatial patterns of disease risk: sheep blowfly strike by Lucilia sericata in Great Britain. Int J Parasitol. 2011, 41: 739-746. 10.1016/j.ijpara.2011.01.012.

    Article  PubMed  Google Scholar 

  25. Illoldi-Rangel P, Rivaldi CL, Sissel B, Trout Fryxell R, Gordillo-Pérez G, Rodríguez-Moreno A, Williamson P, Montiel-Parra G, Sánchez-Cordero V, Sarkar S: Species distribution models and ecological suitability analysis for potential tick vectors of lyme disease in mexico. J Trop Med. 2012, doi:10.1155/2012/959101

    Google Scholar 

  26. Parola P, Raoult D: Tick-borne bacterial diseases emerging in Europe. Clin Microbiol Infect. 2001, 7: 80-83. 10.1046/j.1469-0691.2001.00200.x.

    Article  CAS  PubMed  Google Scholar 

  27. Stojek NM, Dutkiewicz J: Studies on the occurrence of gram negative bacteria in ticks: ixodes ricinus as a potential vector of pasturella. Ann Agric Environ Med. 2004, 11: 319-22.

    PubMed  Google Scholar 

  28. Süss J, Klaus C, Diller R, Schrader C, Wohanka N, Abel U: TBE incidence versus virus prevalence and increased prevalence of the TBE virus in ixodes ricinus removed from humans. Int J Med Microbiol. 2006, 296: 63-68.

    Article  PubMed  Google Scholar 

  29. Eremeeva ME, Oliveira A, Moriarity J, Robinson JB, Tokarevich NK, Antyukova LP, Pyanyh VA, Emeljanova ON, Ignatjeva VN, Buzinov R, Pyankova V, Dasch GA: Detection and identification of bacterial agents in Ixodes persulcatus Schulze ticks from the north western region of Russia. Vector Borne Zoonotic Dis. 2007, 7: 426-436. 10.1089/vbz.2007.0112.

    Article  PubMed  Google Scholar 

  30. Jääskeläinen AE, Tonteri E, Sironen T, Pakarinen L, Vaheri A, Vapalahti O: European subtype tick-borne encephalitis virus in ixodes persulcatus ticks. Emerg Infect Dis. 2011, 17: 323-325. 10.3201/eid1702.101487.

    Article  PubMed Central  PubMed  Google Scholar 

  31. Dantas-Torres F, Chomel BB, Otranto D: Ticks and tick-borne diseases: a one health perspective. Trends Parasitol. 2012, 28: 437-446. 10.1016/j.pt.2012.07.003.

    Article  PubMed  Google Scholar 

  32. Gray JS: The development and seasonal activity of the tick ixodes ricinus, a vector of Lyme borreliosis. Rev Med Vet Entomol. 1991, 79: 323-339.

    Google Scholar 

  33. Korenberg EI, Gorban’ LY, Kovalevskii YV, Frizen VI, Karavanov AS: Risk for human tick-borne encephalitis, borrelioses, and double infection in the pre-ural region of Russia. Emerg Infect Dis. 2001, 7: 459-462.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Bugmyrin SV, Bespyatova LA, Korotkov YS, Burenkova LA, Belova OA, Romanova LI, Kozlovskaya LI, Karganova GG, Ieshko EP: Distribution of ixodes ricinus and I. Persulcatus ticks in southern Karelia (Russia). Ticks Tick Borne Dis. 2013, 4: 57-62. 10.1016/j.ttbdis.2012.07.004.

    Article  PubMed  Google Scholar 

  35. Lindgren E, Tälleklint L, Polfeldt T: Impact of climatic change on the northern distribution limit and population density of the disease-transmitting European tick ixodes ricinus. Environ Health Perspect. 2000, 108: 119-123. 10.1289/ehp.00108119.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  36. Materna J, Daniel M, Danielová V: Altitudinal distribution limit of the tick Ixodes ricinus shifted considerably towards higher altitudes in central Europe: results of three years monitoring in the Krkonose Mts. (Czech Republic). Cent Eur J Publ Health. 2005, 13: 24-28.

    CAS  Google Scholar 

  37. Cadenas FM, Rais O, Jouda F, Douet V, Humair PF, Moret J, Gern L: Phenology of Ixodes ricinus and infection with Borrelia burgdorferi sensu lato along a north- and south-facing altitudinal gradient on Chaumont Mountain, Switzerland. J Med Entomol. 2007, 44: 683-693. 10.1603/0022-2585(2007)44[683:POIRAI]2.0.CO;2.

    Article  Google Scholar 

  38. Lindgren E, Jaenson TGT: Lyme borreliosis in Europe: influences of climate and climate change, epidemiology, ecology and adaptation measures. 2006, Copenhagen, Denmark: WHO Regional Office for Europe

    Google Scholar 

  39. Jore S, Viljugrein H, Hofshagen M, Brun-Hansen H, Kristoffersen AB, Nygard K, Brun E, Ottesen P, Saevik BK, Ytrehus B: Multi-source analysis reveals latitudinal and altitudinal shifts in range of Ixodes ricinus at its northern distribution limit. Parasit Vectors. 2011, 4: 84-10.1186/1756-3305-4-84.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Jaenson TGT, Jaenson DGE, Eisen L, Petersson E, Lindgren E: Changes in the geographical distribution and abundance of the tick Ixodes ricinus during the past 30 years in Sweden. Parasit Vectors. 2012, 5: 8-10.1186/1756-3305-5-8.

    Article  PubMed Central  PubMed  Google Scholar 

  41. IPCC: Climate change 2007: the physical science basis. Contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change. Edited by: Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL. 2007, United Kingdom and New York: Cambridge University Press, Cambridge

    Google Scholar 

  42. Phillips SJ, Anderson RP, Schapire RE: Maximum entropy modeling of species geographic distributions. Ecol Model. 2006, 190: 231-259. 10.1016/j.ecolmodel.2005.03.026.

    Article  Google Scholar 

  43. Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, Hijmans RJ, Huettmann F, John R: Novel methods improve prediction of species’ distributions form occurrence data. Ecography. 2006, 29: 129-151. 10.1111/j.2006.0906-7590.04596.x.

    Article  Google Scholar 

  44. Elith J, Phillips SJ, Hastie T, Dudìk M, En Chee Y, Yates CJ: A statistical explanation of MaxEnt for ecologists. Diversity Distrib. 2011, 17: 43-57. 10.1111/j.1472-4642.2010.00725.x.

    Article  Google Scholar 

  45. Wisz MS, Hijmans RJ, Peterson AT, Graham CH, Guisan A, NCEAS Predicting Species Distributions Working Group: Effects of sample size on the performance of species distribution models. Divers Distrib. 2008, 14: 763-773. 10.1111/j.1472-4642.2008.00482.x.

    Article  Google Scholar 

  46. Oppel S, Meirinho A, Ramírez I, Gardner B, O’Connell AF, Miller PI, Louzao M: Comparison of five modelling techniques to predict the spatial distribution and abundance of seabirds. Biol Conserv. 2012, 156: 94-104.

    Article  Google Scholar 

  47. Estrada-Peña A, Farkas R, Jaenson TG, Koenen F, Madder M, Pascucci I, Salman M, Tarres-Call J, Jongejan F: Association of environmental traits with the geographic ranges of ticks (acari: ixodidae) of medical and veterinary importance in the western palearctic: a digital data set. Exp Appl Acarol. 2013, 59: 351-366. 10.1007/s10493-012-9600-7.

    Article  PubMed Central  PubMed  Google Scholar 

  48. Warren DL, Seifert SN: Ecological niche modeling in maxent: the importance of model complexity and the performance of model selection criteria. Ecol Appl. 2011, 21: 335-342. 10.1890/10-1171.1.

    Article  PubMed  Google Scholar 

  49. Peterson AT, Papeş M, Eaton M: Transferability and model evaluation in ecological niche modeling: a comparison of GARP and maxe nt. Ecography. 2007, 30: 550-560.

    Article  Google Scholar 

  50. Phillips SJ, Dudik M: Modelling of species distributions with maxent: new extensions and a comprehensive evaluation. Ecography. 2008, 31: 161-175. 10.1111/j.0906-7590.2008.5203.x.

    Article  Google Scholar 

  51. Nakicenovic N, Swart R: Emissions scenarios. 2000, United Kingdom, Cambridge: University Press., Cambridge

    Google Scholar 

  52. Scinocca JF, Mcfarlane NA, Lazare M, Li J, Plummer D: Technical note: the CCCma third generation AGCM and its extension into the middle atmosphere. Atmos Chem Phys. 2008, 8: 7055-7074. 10.5194/acp-8-7055-2008.

    Article  CAS  Google Scholar 

  53. Pearson RG: Species’ distribution modeling for conservation educators and practitioners. Synthesis. 2007, New York Am Mus Natl Hist,http://ncep.amnh.org,

    Google Scholar 

  54. Gern L, Jaulhac B: Life cycle of Borrelia burgdorferi sensu lato and trasmission to humans. Lyme borreliosis. Current problem in dermatology. Edited by: Lipsker D. 2009, Switzerland: Karger (Basel), 19-30.

    Google Scholar 

  55. Estrada-Peña A, Ayllón N, de la Fuente J: Impact of climate trends on tick-borne pathogen transmission. Front Physiol. 2012, 3: 1-12.

    Article  Google Scholar 

  56. Andreassen A, Jore S, Cuber P, Dudman S, Tengs T, Isaksen K, Hygen HO, Viljugrein H, Ånestad G, Ottesen P, Vainio K: Prevalence of tick borne encephalitis virus in tick nymphs in relation to climatic factors on the southern coast of Norway. Parasit Vectors. 2012, 5: 177-10.1186/1756-3305-5-177.

    Article  PubMed Central  PubMed  Google Scholar 

  57. Soberon JM, Peterson AT: Interpretation of models of fundamental ecological niches and species’ distributional areas. Biodiv Inform. 2005, 2: 1-10.

    Article  Google Scholar 

  58. Medlock JM, Hansford KM, Bormane A, Derdakova M, Estrada-Peña A, George JC, Golovljova I, Jaenson TG, Jensen JK, Jensen PM: Driving forces for changes in geographical distribution of Ixodes ricinus ticks in Europe. Parasit Vectors. 2013, 6: 1-10.1186/1756-3305-6-1.

    Article  PubMed Central  PubMed  Google Scholar 

  59. Filippova NA: Forms of sympatry and possible ways of microevolution of closely related species of the group Ixodes ricinus-persulcatus (Ixodidae). Acta Zool Lit. 2002, 12: 1392-1657.

    Article  Google Scholar 

  60. Swanson SJ, Neitzel D, Reed KD, Belongia EA: Coinfections acquired from Ixodes Ticks. Clin Microbiol Rev. 2006, 19: 708-727. 10.1128/CMR.00011-06.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Araújo MB, Luoto M: The importance of biotic interactions for modelling species distributions under climate change. Glob Ecol Biogeogr. 2007, 16: 743-753. 10.1111/j.1466-8238.2007.00359.x.

    Article  Google Scholar 

  62. Fischer D, Moeller P, Thomas SM, Naucke TJ, Beierkuhnlein C: Combining climatic projections and dispersal ability: a method for estimating the responses of sandfly vector species to climate change. PLoS Negl Trop Dis. 2011, 5: e1407-10.1371/journal.pntd.0001407.

    Article  PubMed Central  PubMed  Google Scholar 

  63. Sonenshine DE: Biology of ticks, Volume 1. 1991, New York: Oxford University Press

    Google Scholar 

  64. Sonenshine DE: Biology of ticks, Volume 2. 1993, New York: Oxford University Press

    Google Scholar 

  65. Milne A: The ecology of the sheep tick, Ixodes ricinus–host relationships of the tick part 1: review of previous work in Britain. Parasitology. 1949, 39: 167-172. 10.1017/S0031182000083724.

    Article  CAS  PubMed  Google Scholar 

  66. Poulin R, Keeney DB: Host specificity under molecular and experimental scrutiny. Trends Parasitol. 2008, 24: 24-28. 10.1016/j.pt.2007.10.002.

    Article  CAS  PubMed  Google Scholar 

  67. Porretta D, Mastrantonio V, Mona S, Epis S, Montagna M, Sassera D, Bandi C, Urbanelli S: The integration of multiple independent data reveals an unusual response to Pleistocene climatic changes in the hard tick Ixodes ricinus. Mol Eco. 2013, 22 (6): 1666-82. 10.1111/mec.12203.

    Article  Google Scholar 

  68. Korch GW: Geographic dissemination of tick-borne zoonoses. 1994, New York: Oxford University Press

    Google Scholar 

  69. Estrada-Peña A: The relationships between habitat topology, critical scales of connectivity and tick abundance Ixodes ricinus in a heterogeneous landscape in northern Spain. Ecography. 2003, 26: 661-671. 10.1034/j.1600-0587.2003.03530.x.

    Article  Google Scholar 

  70. Estrada-Peña A: Climate, niche, ticks, and models: what they are and how we should interpret them. Parasitol Res. 2008, 103 (Suppl. 1): S87-S95.

    Article  PubMed  Google Scholar 

  71. Santos-Silva MM, Sousa R, Santos AS, Melo P, Encarnaҫão V, Bacellar F: Ticks parasitizing wild birds in Portugal: detection of Rickettsia aeschlimannii, R. helveticaand, R. massiliae. Exp Appl Acarol. 2006, 39: 331-338. 10.1007/s10493-006-9008-3.

    Article  PubMed  Google Scholar 

  72. Hasle G, Bjune G, Edvardsen E, Jakobsen C, Linnehol B, Røer JE, Mehl R, Røed KH, Pedersen J, Leinaas HP: Transport of ticks by migratory passerine birds to Norway. J Parasitol. 2009, 95: 1342-1351. 10.1645/GE-2146.1.

    Article  PubMed  Google Scholar 

  73. De Meeûs T, Beati L, Delaye C, Aeschlimann A, Renaud F: Sex-biased genetic structure in the vector of Lyme disease, Ixodes ricinus. Evolution. 2002, 56: 1802-1807.

    Article  PubMed  Google Scholar 

  74. Paulauskas A, Radzijevskaja J, Rosef O, Turcinaviciene J, Ambrasiene D, Makareviciute D: Genetic variation of ticks (Ixodes ricinus L.) in the Lithuanian and Norwegian populations. Exp Appl Acarol. 2006, 40: 259-270.

    Article  CAS  PubMed  Google Scholar 

  75. Hoogstraal H, Kaiser MN, Traylor MA, Guindy E, Gaber S: Ticks (Ixodidae) on birds migrating from Europe and Asia to Africa, 1959–61. Bull Wld Hlth Org. 1963, 28: 235-262.

    CAS  Google Scholar 

  76. Hoogstraal H, Traylor MA, Gaber S, Malakatis G, Guindy E, Helmy I: Ticks (ixodidae) on migrating birds in Egypt, spring and fall 1962. Bull Wld Hlth Org. 1964, 30: 355-367.

    CAS  Google Scholar 

  77. Ogden NH, Bigras-Poulin M, O’Callaghan CJ, Barker IK, Lindsay LR, Maarouf A, Smoyer-Tomic KE, Waltner-Toews D, Charron D: A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. Int J Parasitol. 2005, 35: 375-389. 10.1016/j.ijpara.2004.12.013.

    Article  CAS  PubMed  Google Scholar 

  78. Bostrom A, O’Connor R, Böhm G, Hanss D, Bodi O, Ekström F, Halder P, Jeschke S, Mack B, Qu M, Rosentrater L, Sandve A, Sælensminde I: Causal thinking and support for climate change policies: international survey findings. Glob Environ Chang. 2012, 22: 210-222. 10.1016/j.gloenvcha.2011.09.012.

    Article  Google Scholar 

  79. Mastrantonio V, Amendolia S, Porretta D, Gaiarsa S, Epis S, Genchi C, Bandi C, Otranto D, Urbanelli S: Annual meeting of the Italian society of ecology, 16th-19th September 2013. Expansion of the climatic niche in the Eurasian tick Ixodes ricinus and I. persulcatus under future climatic changes. Ancona, Italy

  80. Masuzawa T: Terrestrial distribution of the Lyme borreliosis agent borrelia burgdorferi sensu lato in east Asia. Jpn J Infect Dis. 2004, 57: 229-235.

    PubMed  Google Scholar 

  81. Gray JS, Dautel H, Estrada-Peña A, Kahl O, Lindgren E, Gray JS, Dautel H, Estrada-Peña A, Kahl O, Lindgren E: Effects of climate change on ticks and tick-borne diseases in Europe. Interdiscip Perspect Infect Dis. 2009, 2009: 593232-

    PubMed Central  CAS  PubMed  Google Scholar 

  82. Gray JS: Ixodes ricinus seasonal activity: implications of global warming indicated by revisiting tick and weather data. Int J Med Microbiol. 2008, 298: 19-24.

    Article  Google Scholar 

  83. Dautel H, Dippel C, Kämmer D, Werkhausen A, Kah O: Winter activity of ixodes ricinus in a Berlin forest. Int J Med Microbiol. 2008, 298: 50-54.

    Article  Google Scholar 

  84. Dantas-Torres F, Otranto D: Seasonal dynamics of ixodes ricinus on ground level and higher vegetation in a preserved wooded area in southern Europe. Vet Parasitol. 2013, 192: 253-258. 10.1016/j.vetpar.2012.09.034.

    Article  PubMed  Google Scholar 

  85. Gern L, Rais O: Efficient transmission of borrelia burgdorferi between cofeeding ixodes ricinus ticks (acari: ixodidae). J Med Entomol. 1996, 33: 189-92.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgment

We thank Alessandra Spanò for technical assistance; the editors and four anonymous reviewers for their comments and suggestions. The work was funded through grants from the Ministero dell'Istruzione, dell'Università e della Ricerca (MIUR).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sandra Urbanelli.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

DP, VM, SA SU, conceived and designed the study. DP, VM, SA drafted the manuscript. VM, SA analyzed the data. SG, SE, CG, CB, DO critically revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

13071_2013_1046_MOESM1_ESM.doc

Additional file 1: List of bioclimatic variables used to develop the Species Distribution Models and variable percent contribution to construct the models.(DOC 38 KB)

13071_2013_1046_MOESM2_ESM.doc

Additional file 2: Increase of geographic distribution of Ixodes ricinus , predicted by species distribution modelling. The values under Current, 2050 and 2080 (A2 and B2 scenarios) columns represent the areas in square km. The values in the Ratio column represent the increase of area as ratio between future and current conditions. (DOC 30 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Porretta, D., Mastrantonio, V., Amendolia, S. et al. Effects of global changes on the climatic niche of the tick Ixodes ricinus inferred by species distribution modelling. Parasites Vectors 6, 271 (2013). https://doi.org/10.1186/1756-3305-6-271

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1756-3305-6-271

Keywords