Next Article in Journal
Transplantation of Fibroblast Sheets with Blood Mononuclear Cell Culture Exerts Cardioprotective Effects by Enhancing Anti-Inflammation and Vasculogenic Potential in Rat Experimental Autoimmune Myocarditis Model
Next Article in Special Issue
Predicting Current Potential Distribution and the Range Dynamics of Pomacea canaliculata in China under Global Climate Change
Previous Article in Journal
Phenotypic and Molecular Characteristics of the MDR Efflux Pump Gene-Carrying Stenotrophomonas maltophilia Strains Isolated in Warsaw, Poland
Previous Article in Special Issue
Projecting the Potential Distribution of Glossina morsitans (Diptera: Glossinidae) under Climate Change Using the MaxEnt Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Projecting the Potential Distribution Areas of Ixodes scapularis (Acari: Ixodidae) Driven by Climate Change

1
State Key Laboratory of Infectious Disease Prevention and Control, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 102206, China
2
School of Public Health and Health Management, Shandong First Medical University and Shandong Academy of Medical Sciences, Jinan 250117, China
3
Shandong University Climate Change and Health Center, School of Public Health, Shandong University, Jinan 250012, China
*
Author to whom correspondence should be addressed.
Biology 2022, 11(1), 107; https://doi.org/10.3390/biology11010107
Submission received: 29 November 2021 / Revised: 6 January 2022 / Accepted: 7 January 2022 / Published: 10 January 2022
(This article belongs to the Special Issue The Impact of Climate Change on Ecosystems and Human Health)

Abstract

:

Simple Summary

Ixodes scapularis is a vector of tick-borne diseases. Climate change is one of the main factors affecting the distribution of I. scapularis. We used the maximum entropy model to project the potential distribution and future trends of I. scapularis. The potential suitable area of I. scapularis is dynamically changing in the context of climate change and precipitation in May makes the greatest contribution to such expansion.

Abstract

Ixodes scapularis is a vector of tick-borne diseases. Climate change is frequently invoked as an important cause of geographic expansions of tick-borne diseases. Environmental variables such as temperature and precipitation have an important impact on the geographical distribution of disease vectors. We used the maximum entropy model to project the potential geographic distribution and future trends of I. scapularis. The main climatic variables affecting the distribution of potential suitable areas were screened by the jackknife method. Arc Map 10.5 was used to visualize the projection results to better present the distribution of potential suitable areas. Under climate change scenarios, the potential suitable area of I. scapularis is dynamically changing. The largest suitable area of I. scapularis is under SSP3-7.0 from 2081 to 2100, while the smallest is under SSP5-8.5 from 2081 to 2100, even smaller than the current suitable area. Precipitation in May and September are the main contributing factors affecting the potential suitable areas of I. scapularis. With the opportunity to spread to more potential suitable areas, it is critical to strengthen surveillance to prevent the possible invasion of I. scapularis.

1. Introduction

Ixodes scapularis (I. scapularis) is a vector of many tick-borne diseases [1]. In the United States, for example, I. scapularis is a vector of seven human pathogens [2]. Climate change is believed to be a factor in the increased incidence of tick-borne diseases [3,4]. Since Lyme disease was first described in the 1970s in the United States, its geographical range has gradually expanded from northeastern to western and southern, and even to Eurasia and other regions [5,6]. Evidence is mounting that climate change affects the spread of Lyme disease [7]. As the main vector of Lyme disease, I. scapularis is also affected by climate change. Studies have found that the geographic distribution of I. scapularis and the pathogens it carrys has been expanding over the past 20 years [2]. The impact of climate change on tick vectors and diseases is expected to be profound [8].
Currently, climate change is one of the major international concerns. The Intergovernmental Panel on Climate Change (IPCC) published Shared Socioeconomic Pathways (SSPs) [9], in 2010, which were developed on the basis of representative concentration pathways (RCPs) to quantitatively describe the relationship between climate change and socioeconomic development paths and reflect future climate change adaptation and mitigation challenges for society. Many researchers have used SSPs to simulate future climate change to project suitable areas of different species [10,11,12,13,14]. It is critical to understand the impact of climate change on vector distribution, which can provide a theoretical basis for the prevention and control of vector-borne diseases.
Species distribution models are widely used in various research fields. The most common species distribution models include GARP, EnFA, BIOCLIM, DOMAIN, MaxEnt, etc. [15]. The MaxEnt model is the ecological model with the best prediction effect at present. Compared with other models, its prediction is more consistent with the real world [16]. The MaxEnt model can accurately predict the potential distribution of species through the data of the species’ distribution points and environmental variables, providing an important basis for scientific research. Even a few distribution points can be accurately predicted as compared with other models [17]. The MaxEnt model can self-check the simulated prediction by using the receiver operating characteristic curve (ROC) of the jackknife test. The jackknife method has also been selected to measure the contribution of each environmental variable and to judge the different effects of different environments on species distribution [18,19].
This study used the MaxEnt model, based on the actual distribution sites of I. scapularis and environmental variables, to project the potential distribution of I. scapularis internationally and to find the main contributing environmental variables under future climate scenarios.

2. Materials and Methods

2.1. Software and Geographic Data

Maxent software (version 3.4.1, https://biodiversityinformatics.amnh.org/, accessed on 23 May 2021) was used to project the potential distribution of I. scapularis. ArcGIS software (version 10.5, https://www.arcgis.com/, accessed on 23 May 2021) was used to visualize the projection results. R software (version 4.0.0, http://www.r-project.org/, accessed on 31 May 2021) and DIVA-GIS software (version 7.5.0, https://www.diva-gis.org/, accessed on 31 May 2021) were used to adjust model parameters. The world map (scale 1:10,000,000) was downloaded from the Natural Earth (http://www.naturalearthdata.com, accessed on 13 June 2021). The scale of the map in this article represents mid-latitude distances.

2.2. Distribution Data Collection and Processing of I. scapularis

The distribution data of I. scapularis were obtained from the Global Biodiversity Information Facility (GBIF) (https://www.gbif.org/. accessed on 13 June 2021). Via Google Earth (https://earth.google.com, accessed on 15 June 2021), the accuracy of coordinate points were checked and points with apparent errors in geographical coordinates were removed. To eliminate the possible spatial sampling bias caused by too dense distribution points in a certain region, Proximity in ArcMap was used to create a 2.5 km buffer with points as factors to screen the distribution data. After the above processing, we obtained 1140 points of I. scapularis. All sample point data were saved in a CSV file, which contained three columns, i.e., species name, longtitude and latitude coordinate values. It can be seen from the figure that I. scapularis is mainly distributed in the eastern part of the United States (Figure 1).

2.3. Environmental Data

The following were downloaded: 19 bioclimatic variables, 36 historical climate data (spatial resolution 2.5 min) of monthly minimum temperature (Tmin1–12), monthly maximum temperature (Tmax1–12), monthly total precipitation (Prec1–12), and elevation from Worldclim v2.1 released in January 2020 (http://worldclim.org/version2, accessed on 16 June 2021) (Table S1). The shared socioeconomic paths (SSPs) of BCC-CSM2-MR, i.e., SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5 (spatial resolution 2.5 min), were simulated as future climate conditions.

2.4. Selection of Model Parameters

Recent studies have shown that the MaxEnt model using default parameters may not be optimal, especially when there are few species distribution points [20]. We need to adjust the parameters to prevent the model from being too complex and/or overfitting [21,22]. Maxent optimizes the model by modifying the two parameters of feature classes (FC) and regularization multiplier (RM). The FC consist of lineal (L), quadratic (Q), product (P), threshold (T), and hinge (H) [23]. The RM parameters are set to eight levels: 0.5, 1, 1.5, 2, 2.5, 3, 3.5, and 4. To avoid overfitting/complexity, we used the R software package “ENMeval” to obtain the best model through Akaike information criterion (AIC). When AICc is the lowest (i.e., AICc = 0), the model is considered to be the most appropriate [24]. The FC and RM of the MaxEnt model in this study were LQHPT and 4.

2.5. Project the Potential Suitable Area

The occurrence points and environment variables were imported into the MaxEnt model; 25% distribution points were randomly selected as testing data, and the remaining 75% as training data. The maximum iteration mode was “bootstrap”, and the maximum number of repetitions was 10,000. The number of repeated training data was set to 20 times. The threshold rule was used to select the “10 percentile training presence”. There may be multicollinearity between environmental variables [25]. In order to avoid overfitting the MaxEnt model, we used the Pearson coefficient to test the potential correlation of environmental variables [26]. When Pearson’s |r| > 0.8, it was considered to be highly correlated. According to the contribution of environmental variables, we selected them from high to low, and eliminated the highly related variables. The receiver operating characteristic curve (ROC) was used to test the projection accuracy [18,23]. When the area under curve (AUC) was less than 0.7, the model accuracy was considered to be low; when it was between 0.7 and 0.9, the model accuracy was considered to be medium; and when it was greater than 0.9, the model accuracy is considered to be high [27,28]. The output format was ASCII grid layer, and the range of fitness index was 0–1. The running results of MaxEnt were loaded into ArcMap 10.5 for reclassification and visual expression. The distribution area of I. scapularis was divided into unsuitable, high suitable, middle suitable, and low suitable areas by the Jenks natural breaks classification method in ArcGIS. The criteria are: 0–0.07, unsuitable area; 0.08–0.29, low suitable area; 0.30–0.56, middle suitable area; 0.57–1.00, high suitable area.

3. Results

3.1. The Main Contributing Environmental Variables

The geographical distribution of I. scapularis was most affected by five climatic factors. They were precipitation in May (Prec05), precipitation in September (Prec09), precipitation of the driest month (Bio14), temperature seasonality (Bio4) and the mean diurnal range (Bio2) (Table 1). Prec05 contains the most useful information, while Bio4 contains the most information not present in the other variables (Figure 2).

3.2. The Suitable Areas of I. scapularis under Near Current Climatic Condition

Under current climatic condition, the suitable areas of the I. scapularis in the world are mainly in parts of America, Europe, and Asia. In North America, it is mainly distributed in the eastern United States, the southeast and southwest of Canada, the northeast of Mexico, and the Bahamas. In South America, it is mainly distributed in Uruguay, southern Brazil, northeastern Argentina, and southeastern Paraguay. In Europe, it is mainly distributed in Switzerland, the southern part of France, Italy, Austria, Germany, etc. In Asia, it is mainly distributed in the southeastern part of China, South Korea, and Japan. It is also distributed in New Zealand and southeastern Australia in Oceania (Figure 3).

3.3. The Suitable Areas of I. scapularis under Future Climate Change Scenarios

Our projected results show a dynamic change in the suitable area of I. scapularis around the world from 2021 to 2100 under four future climate change scenarios (Figure 4, Figure 5, Figure 6 and Figure 7). Compared with the current suitable area, the potential suitable area will increase under the scenarios of SSP1-2.6, SSP2-4.5, and SSP3-7.0 (1.17–10.41%), while under SSP5-8.5, the potential suitable areas in the periods 2021–2040 and 2041–2060 will increase (4.92% and 5.66%), but in the periods 2061–2080 and 2081–2100 will decrease (−1.79% and −3.38%) (Table 2). The potential suitable area of I. scapularis will be the largest in the period 2081–2100 SSP3-7.0. Compared with the current suitable area, it will expand northward in Canada and Greenland. In addition, western Russia, southwest and northeast China, northeast Argentina, and southeast Chile whill also have significant expansion (Figure 8). In the period 2081–2100 SSP5-8.5, the area of the potential suitable area will be the smallest. Compared with the current suitable area, southern Europe, South America, Brazil, Uruguay, and other regions, as well as central Asia, will have a significant contraction (Figure 9).
Under various climate change scenarios, the minimum high suitable area was projected in the period 2021–2040 SSP1-2.6, reducing to an area of 2.50 × 106 km2. Compared with the current high suitable area, it decreased by 11.73% (Figure 10a). Under the SSP3-7.0, the projected high suitable area is the largest from 2081 to 2100, reaching 3.31 × 106 km2, an increase of 16.79% over the current climate conditions (Figure 10c). Middle suitable areas show a downward trend under SSP2-4.5 (Figure 10b). Under the SSP5-8.5 scenario, there will be a significant contraction of the low suitability area (Figure 10d).

3.4. The Accuracy of the MaxEnt Model

The ROC was used to test the projection accuracy. The AUC of the MaxEnt model was 0.950 (Figure 11), and the standard deviation was 0.001, which showed that the accuracy of the model was high.

4. Discussion

In the two-year life cycle of I. scapulalis, more than 95% of the time is spent outside the host, either looking for the host or diapause [29]. Because of their high surface area to volume ratio, I. scapulalis is susceptible to desiccation when questing for hosts [30]. Therefore, I. scapulalis needs a suitable habitat to survive and reproduce. We found that the most important environmental variables affecting the potential distribution of I. scapulalis were precipitation in May (Prec05), precipitation in September (Prec09), precipitation in the driest month (Bio14), temperature seasonality (Bio4), and the mean diurnal range (Bio2). Micah et al. [31] found that the main environmental variables affecting the distribution of I. scapulalis in the USA were the max temperature of the warmest month (Bio5) and precipitation of the warmest quarter (Bio18). I. scapularis expand its range into northwestern Minnesota, as well as central and northern Michigan in the USA. Johnson et al. [32] predicted the habitat suitability of I. scapularis in Minnesota, in which land cover, temperature annual range (Bio7), and precipitation of the wettest quarter (Bio16) contributed the most. Their model predicted that the habitat of I. scapularis extended up into the Minnesota River Valley. Although the main environmental variables derived from this study are slightly different from previous studies, the projected distribution of suitable areas in the Americas is consistent. One possible reason is that the study only included bioclimatic variables, not variables such as monthly total precipitation. The second reason may be that it refers to WorldClim1.4, which is the newer version 2.1 used in this study. Micah et al. [31] and Johnson et al. [32] mainly focused on the Americas, while this study projected the distribution of the global suitable areas of I. scapularis. Other studies have shown that Ontario, Canada, will become a high-risk area for I. scapulalis in the 2050s [4,33]. This is consistent with our projections.
Dry years may result in reduced nymphal activity and mortality as a result of desiccation [30]. In hot summer, enough precipitation is needed to offset the death of I. scapulalis from heat [32]. However, heavy precipitation associated with tropical cyclones is projected to be higher at 2 °C as compared with 1.5 °C of global warming [34], that is, global warming will bring more extreme weather such as droughts and floods [35]. If precipitation and temperature exceed the tolerance of I. scapulalis, then, these may be the reasons for the projected potential contraction of distribution areas under the SSP5-8.5 scenario from 2061 to 2100. When studying the ecological factors related to I. scapulalis, it is critical to study its microhabitat characteristics [36]. Some data show that soil properties, such as litter depth, can affect the survival and density of I. scapulalis. Proper litter and canopy cover are helpful to stabilize the microclimate of ticks [37,38].
The distribution of I. scapularis is indeed affected by many factors, including vegetation and host. Such variables including soil properties, vegetation types, and the main host distribution are not available in this study, which is a limitation.
In this study, the MaxEnt was used to project the potential distribution of I. scapulalis. The accuracy of the model was greater than 0.9 and the fitting was good. The MaxEnt model is popular and has been used by many researchers to predict the distribution of species. In the future, we need more data for simulation scenario analysis to make the MaxEnt model more accurate and reliable. In addition, MaxEnt needs accurate longitude and latitude information, but these can hardly be obtained from the previous data or literature. Subsequent research should consider taking the region as the research object [39] and using the enhanced regression tree model for projection.
With global climate change, the potential distribution of I.scapulalis is also changing. The projection results show that some areas such as Greenland and southern Chile have suitable areas, that is, the potential suitable areas of I.scapulalis is expanding to higher latitudes. Climate change is indeed affecting the spread of vector-borne infectious diseases, so we should actively respond to climate change and reduce the harm brought by climate change to human beings.

5. Conclusions

We used climate conditions and occurrence points to project the potential distribution of I.scapulalis under MaxEnt. It is found that the range of suitable areas generally showed an upward trend under the future climate change scenarios, and the largest area is in the periods 2081–2100 SSP3-7.0, but the suitable area of SSP5-8.5 shows a downward trend in the periods 2061–2080 and 2081–2100. Precipitation and temperature are the main factors affecting the distribution of suitable areas. With global climate change, the potential suitable areas of I.scapulalis is likely expanding to higher latitudes. Climate change is indeed affecting the distribution of ticks’ suitable areas.. It is critical to actively respond to climate change.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology11010107/s1, Table S1: Climate variables used in projecting the potential geographic distribution of I. scapularis, Table S2: The distribution points of I. scapulalis for model operation.

Author Contributions

Conceptualization, L.Z.; methodology, L.Z. and D.M.; software, R.Z. and D.M.; formal analysis, L.Z. and D.M.; data curation, L.Z. and D.M.; writing—original draft preparation, L.Z. and C.L.; writing—review and editing, C.L., R.Z., and Q.L.; supervision, Q.L.; funding acquisition, J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (grant number 2020YFC1200101) and the National Natural Science Foundation of China (grant number 32090023).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Environmental variables can be downloaded from Worldclim v2.1 released in January 2020 (http://worldclim.org/version2, accessed on 16 June 2021). Occurrence records can be downloaded from GBIF (https://www.gbif.org, accessed on 13 June 2021).

Acknowledgments

We acknowledge Juan Shi and Tai Gao, working at Beijing Forestry University, for their assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ogden, N.H.; Maarouf, A.; Barker, I.K.; Bigras-Poulin, M.; Lindsay, L.R.; Morshed, M.G.; O’Callaghan, C.J.; Ramay, F.; Waltner-Toews, D.; Charron, D.F. Climate change and the potential for range expansion of the Lyme disease vector Ixodes scapularis in Canada. Int. J. Parasitol. 2006, 36, 63–70. [Google Scholar] [CrossRef]
  2. Eisen, R.J.; Eisen, L. The Blacklegged Tick, Ixodes scapularis: An Increasing Public Health Concern. Trends Parasitol. 2018, 34, 295–309. [Google Scholar] [CrossRef]
  3. Sanchez-Vicente, S.; Tagliafierro, T.; Coleman, J.L.; Benach, J.L.; Tokarz, R. Polymicrobial Nature of Tick-Borne Diseases. mBio 2019, 10, e02055-19. [Google Scholar] [CrossRef] [Green Version]
  4. Ogden, N.H.; Radojevic, M.; Wu, X.; Duvvuri, V.R.; Leighton, P.A.; Wu, J. Estimated effects of projected climate change on the basic reproductive number of the Lyme disease vector Ixodes scapularis. Environ. Health Perspect. 2014, 122, 631–638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bacon, R.M.; Kugeler, K.J.; Mead, P.S. Surveillance for Lyme disease—United States, 1992–2006. MMWR Surveill. Summ. 2008, 57, 1–9. [Google Scholar] [PubMed]
  6. Semenza, J.C.; Suk, J.E. Vector-borne diseases and climate change: A European perspective. FEMS Microbiol. Lett. 2018, 365, fnx244. [Google Scholar] [CrossRef]
  7. Couper, L.I.; MacDonald, A.J.; Mordecai, E.A. Impact of prior and projected climate change on US Lyme disease incidence. Glob. Chang. Biol. 2021, 27, 738–754. [Google Scholar] [CrossRef] [PubMed]
  8. Levy, S. Ticking time bomb? Climate change and Ixodes scapularis. Environ. Health Perspect. 2014, 122, A168. [Google Scholar] [CrossRef] [Green Version]
  9. O’Neill, B.; Kriegler, E.; Ebi, K.L.; Kemp-Benedict, E.; Riahi, K.; Rothman, D.S.; Ruijven, B.V.; Vuuren, D.V.; Birkmann, J.; Kok, K. The roads ahead: Narratives for shared socioeconomic pathways describing world futures in the 21st century. Glob. Environ. Change 2017, 42, 169–180. [Google Scholar] [CrossRef] [Green Version]
  10. Duan, R.Y.; Kong, X.Q.; Huang, M.Y.; Varela, S.; Ji, X. The potential effects of climate change on amphibian distribution, range fragmentation and turnover in China. PeerJ 2016, 4, e2185. [Google Scholar] [CrossRef] [Green Version]
  11. Alkishe, A.A.; Peterson, A.T.; Samy, A.M. Climate change influences on the potential geographic distribution of the disease vector tick Ixodes ricinus. PLoS ONE 2017, 12, e0189092. [Google Scholar] [CrossRef] [PubMed]
  12. Ma, D.; Lun, X.; Li, C.; Zhou, R.; Zhao, Z.; Wang, J.; Zhang, Q.; Liu, Q. Predicting the Potential Global Distribution of (Acari: Ixodidae) under Near Current and Future Climatic Conditions, Using the Maximum Entropy Model. Biology 2021, 10, 1057. [Google Scholar] [CrossRef] [PubMed]
  13. Zhou, R.; Gao, Y.; Chang, N.; Gao, T.; Ma, D.; Li, C.; Liu, Q. Projecting the Potential Distribution of (Diptera: Glossinidae) under Climate Change Using the MaxEnt Model. Biology 2021, 10, 1150. [Google Scholar] [CrossRef] [PubMed]
  14. Gao, T.; Shi, J. The Potential Global Distribution of Sirex juvencus (Hymenoptera: Siricidae) under Near Current and Future Climatic Conditions as Predicted by the Maximum Entropy Model. Insects 2021, 12, 222. [Google Scholar] [CrossRef] [PubMed]
  15. Cauwer, V.D.; Muys, B.; Revermann, R.; Trabucco, A. Potential, realised, future distribution and environmental suitability for Pterocarpus angolensis DC in southern Africa. For. Ecol. Manag. 2014, 315, 211–226. [Google Scholar] [CrossRef]
  16. Gastón, A.; Garcia-Vinas, J.I.; Bravo-Fernandez, A.J.; López-Leiva, C.; Oliet, J.A.; Roig, S.; Serrada, R. Species distribution models applied to plant species selection in forest restoration: Are model predictions comparable to expert opinion? New Forest 2014, 45, 641–653. [Google Scholar] [CrossRef]
  17. Phillips, S.J.; Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  18. Wang, Y.; Xie, B.; Wan, F.; Xiao, Q.; Dai, L. Application of ROC curve analysis in evaluating the performance of alien species’ potential distribution models. Biodivers. Sci. 2007, 15, 365–372. [Google Scholar]
  19. Xiang, L.; Fucheng, Z.; Shuai, H.; Yuqing, W.; Lin, M. Research Progress Analysison the Comprehensive Application of MaxEnt Model. J. Green Sci. Technol. 2020, 18, 14–17. [Google Scholar]
  20. Soucy, J.R.; Slatculescu, A.M.; Nyiraneza, C.; Ogden, N.H.; Leighton, P.A.; Kerr, J.T.; Kulkarni, M.A. High-Resolution Ecological Niche Modeling of Ixodes scapularis Ticks Based on Passive Surveillance Data at the Northern Frontier of Lyme Disease Emergence in North America. Vector Borne Zoonotic Dis. 2018, 18, 235–242. [Google Scholar] [CrossRef] [Green Version]
  21. Elith, J.; Kearney, M.; Phillips, S. The art of modelling range-shifting species. Methods Ecol. Evol. 2010, 1, 330–342. [Google Scholar] [CrossRef]
  22. Shcheglovitova, M.; Anderson, R.P. Estimating optimal complexity for ecological niche models: A jackknife approach for species with small sample sizes. Ecol. Model. 2013, 269, 9–17. [Google Scholar] [CrossRef]
  23. Morales, N.S.; Fernández, I.C.; Baca-González, V. MaxEnt’s parameter configuration and small samples: Are we paying attention to recommendations? A systematic review. PeerJ 2017, 5, e3093. [Google Scholar] [CrossRef]
  24. Dan, L.W.; Seifert, S.N. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecol. Appl. 2011, 21, 335–342. [Google Scholar]
  25. Raghavan, R.K.; Peterson, A.T.; Cobos, M.E.; Ganta, R.; Foley, D. Current and Future Distribution of the Lone Star Tick, Amblyomma americanum (L.) (Acari: Ixodidae) in North America. PLoS ONE 2019, 14, e0209082. [Google Scholar] [CrossRef]
  26. Graham, M.H. Confronting multicollinearity in ecological multiple regression. Ecology 2003, 84, 2809–2815. [Google Scholar] [CrossRef] [Green Version]
  27. Franklin, J. Ecological understanding of species distributions. In Mapping Species Distributions (Spatial Inference and Prediction); Cambridge University Press: Cambridge, UK, 2010; pp. 34–52. [Google Scholar]
  28. Franklin, J. Mapping Species Distributions (Spatial Inference and Prediction) || Ecological Understanding of Species Distributions; Cambridge University Press: Cambridge, UK, 2010; pp. 34–52. [Google Scholar]
  29. Ostfeld, R.S.; Brunner, J.L. Climate change and Ixodes tick-borne diseases of humans. Philos. Trans. R Soc. Lond. B Biol. Sci. 2015, 370, 20140051. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Berger, K.A.; Ginsberg, H.S.; Gonzalez, L.; Mather, T.N. Relative humidity and activity patterns of Ixodes scapularis (Acari: Ixodidae). J. Med. Entomol. 2014, 51, 769–776. [Google Scholar] [CrossRef]
  31. Hahn, M.B.; Jarnevich, C.S.; Monaghan, A.J.; Eisen, R.J. Modeling the Geographic Distribution of Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae) in the Contiguous United States. J. Med. Entomol. 2016, 53, 1176–1191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Johnson, T.L.; Bjork, J.K.H.; Neitzel, D.F.; Dorr, F.M.; Schiffman, E.K.; Eisen, R.J. Habitat Suitability Model for the Distribution of Ixodes scapularis (Acari: Ixodidae) in Minnesota. J. Med. Entomol. 2016, 53, 598–606. [Google Scholar] [CrossRef] [Green Version]
  33. Mcpherson, M.Y.; Garcíagarcía, A.; Cuestavalero, F.J.; Beltrami, H.; Hansenketchum, P.; Macdougall, D.; Ogden, N.H. Expansion of the Lyme Disease Vector Ixodes scapularis in Canada inferred from CMIP5 Climate Projections. In Proceedings of the EGU General Assembly Conference, Vienna, Austria, 23–28 April 2017. [Google Scholar]
  34. Hoegh-Guldberg, O.; Jacob, D.; Taylor, M.; Guillén Bolaños, T.; Bindi, M.; Brown, S.; Camilloni, I.A.; Diedhiou, A.; Djalante, R.; Ebi, K.; et al. The human imperative of stabilizing global climate change at 1.5 °C. Science 2019, 365, eaaw6974. [Google Scholar] [CrossRef] [Green Version]
  35. Feng, X.; Liu, C.; Xie, F.; Lu, J.; Chen, B. Precipitation Characteristic Changes due to Global Warming in a High-Resolution (16-km) ECMWF Simulation. Q. J. R. Meteorol. Soc. 2018, 145, 303–317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Estrada-PeñA, A.; José, D. The ecology of ticks and epidemiology of tick-borne viral diseases. Antiviral Res. 2014, 108, 104–128. [Google Scholar] [CrossRef] [PubMed]
  37. Clow, K.M.; Ogden, N.H.; Lindsay, L.R.; Michel, P.; Pearl, D.L.; Jardine, C.M. The influence of abiotic and biotic factors on the invasion of Ixodes scapularis in Ontario, Canada. Ticks Tick Borne Dis. 2017, 8, 554–563. [Google Scholar] [CrossRef] [PubMed]
  38. Lubelczyk, C.B.; Elias, S.P.; Rand, P.W.; Holman, M.S.; Lacombe, E.H.; Smith, R.P. Habitat Associations of Ixodes scapularis (Acari: Ixodidae) in Maine. Environ. Entomol. 2004, 33, 900–906. [Google Scholar] [CrossRef]
  39. Guoping, Z. Spatial Distribution of Ticks and Their Hazard Prediction in China. Ph.D. Thesis, Academy of Military Sciences, Beijing, China, 2018. [Google Scholar]
Figure 1. Current occurrence of Ixodes scapularis.
Figure 1. Current occurrence of Ixodes scapularis.
Biology 11 00107 g001
Figure 2. The results of the jackknife test of variable importance.
Figure 2. The results of the jackknife test of variable importance.
Biology 11 00107 g002
Figure 3. The suitable areas of I. scapularis under near current climatic condition.
Figure 3. The suitable areas of I. scapularis under near current climatic condition.
Biology 11 00107 g003
Figure 4. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 1–2.6 (sustainability and radiative forcing of 2.6 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Figure 4. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 1–2.6 (sustainability and radiative forcing of 2.6 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Biology 11 00107 g004
Figure 5. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 2–4.5 (middle of the road and forcing of 4.5 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Figure 5. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 2–4.5 (middle of the road and forcing of 4.5 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Biology 11 00107 g005
Figure 6. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 3–7.0 (regional rivalry and forcing of 7.0 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Figure 6. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 3–7.0 (regional rivalry and forcing of 7.0 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Biology 11 00107 g006
Figure 7. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 5–8.5 (fossil-fuel development and forcing of 8.5 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Figure 7. The suitable areas of I. scapularis around the world under shared socioeconomic pathway 5–8.5 (fossil-fuel development and forcing of 8.5 W/m2 to 2100) during different periods of the 21st century: (a) 2021–2040; (b) 2041–2060; (c) 2061–2080; (d) 2081–2100.
Biology 11 00107 g007
Figure 8. Comparison of the current suitable area with the period 2081–2100 under SSP3-7.0.
Figure 8. Comparison of the current suitable area with the period 2081–2100 under SSP3-7.0.
Biology 11 00107 g008
Figure 9. Comparison of the current suitable area with the period 2081–2100 under SSP5-8.5.
Figure 9. Comparison of the current suitable area with the period 2081–2100 under SSP5-8.5.
Biology 11 00107 g009
Figure 10. The change rate of future suitable areas as compared with current suitable areas and the four shared socio-economic pathways: (a) SSP1-2.6; (b) SSP2-4.5; (c) SSP3-7.0; (d) SSP5-8.5.
Figure 10. The change rate of future suitable areas as compared with current suitable areas and the four shared socio-economic pathways: (a) SSP1-2.6; (b) SSP2-4.5; (c) SSP3-7.0; (d) SSP5-8.5.
Biology 11 00107 g010
Figure 11. The AUC projected by the MaxEnt model of I. scapularis.
Figure 11. The AUC projected by the MaxEnt model of I. scapularis.
Biology 11 00107 g011
Table 1. The high contributing environmental variables of the potential suitable area of Ixodes scapularis.
Table 1. The high contributing environmental variables of the potential suitable area of Ixodes scapularis.
VariablesDescriptionUnitContribution (%)
Prec05Precipitation in Maymm39.6
Prec09Precipitation in Septembermm22.0
Bio14Precipitation of the driest monthmm19.6
Bio4Temperature seasonality (standard deviation × 100)\17.0
Bio2Mean diurnal range (mean of monthly (max temp—min temp))°C1.80
Table 2. Current and future suitable areas of I. scapularis under different climatic conditions (×106 km2).
Table 2. Current and future suitable areas of I. scapularis under different climatic conditions (×106 km2).
Climate ScenarioPeriodLow
Suitable Area
Middle
Suitable Area
High
Suitable Area
Total AreaArea ChangeArea Change Ratio (%)
current1970–20006.632.642.8312.100.00
SSP1-2.62021–20406.802.942.5012.240.141.17
2041–20607.592.642.6212.860.756.21
2061–20806.912.532.8212.260.151.25
2081–21007.312.873.0513.221.129.26
SSP2-4.52021–20407.252.782.6412.670.574.67
2041–20606.972.752.9812.700.604.97
2061–20806.942.552.9412.430.332.69
2081–21007.062.562.9012.530.423.48
SSP3-7.02021–20406.912.712.6512.270.171.39
2041–20607.023.062.9613.040.947.77
2061–20807.042.772.8912.700.604.93
2081–21007.412.643.3113.361.2610.41
SSP5-8.52021–20407.102.742.8612.700.604.92
2041–20607.122.762.9112.790.695.66
2061–20806.322.872.7011.89−0.22−1.79
2081–21006.172.433.0911.69−0.41−3.38
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, L.; Ma, D.; Li, C.; Zhou, R.; Wang, J.; Liu, Q. Projecting the Potential Distribution Areas of Ixodes scapularis (Acari: Ixodidae) Driven by Climate Change. Biology 2022, 11, 107. https://doi.org/10.3390/biology11010107

AMA Style

Zhang L, Ma D, Li C, Zhou R, Wang J, Liu Q. Projecting the Potential Distribution Areas of Ixodes scapularis (Acari: Ixodidae) Driven by Climate Change. Biology. 2022; 11(1):107. https://doi.org/10.3390/biology11010107

Chicago/Turabian Style

Zhang, Lu, Delong Ma, Chao Li, Ruobing Zhou, Jun Wang, and Qiyong Liu. 2022. "Projecting the Potential Distribution Areas of Ixodes scapularis (Acari: Ixodidae) Driven by Climate Change" Biology 11, no. 1: 107. https://doi.org/10.3390/biology11010107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop