Systematic review of distribution models for Amblyomma ticks and Rickettsial group pathogens
============================================================================================

* Catherine A. Lippi
* Holly D. Gaff
* Alexis L. White
* Sadie J. Ryan

## Abstract

The rising prevalence of tick-borne diseases in humans in recent decades has called attention to the need for more information on geographic risk for public health planning. Species distribution models (SDMs) are an increasingly utilized method of constructing potential geographic ranges.There are many knowledge gaps in our understanding of risk of exposure to tick-borne pathogens, particularly for those in the rickettsial group. Here, we conducted a systematic review of the SDM literature for rickettsial pathogens and tick vectors in the genus *Amblyomma*. Of the 174 reviewed papers, only 24 studies used SDMs to estimate the potential extent of vector and/or pathogen ranges. The majority of studies (79%) estimated only tick distributions using vector presence as a proxy for pathogen exposure. Studies were conducted at different scales and across multiple continents. Few studies undertook original data collection, and SDMs were mostly built with presence-only datasets from public database or surveillance sources. While we identify agap in knowledge, this may simply reflect a lag in new data acquisition and a thorough understanding of the tick-pathogen ecology involved.

Keywords
*   *Amblyomma*
*   *Rickettsia*
*   species distribution model
*   PRISMA

## Introduction

Tick-borne diseases are a global threat to public health, posing risks to both humans and domesticated animals. In recent years there have been documented increases in tick-borne diseases both in the United States and around the world. Much of this burden can be attributed to Lyme disease in the United States, Europe, and northern Asia. However, in the past 20 years,identification of previously unrecognized pathogens has revealed a great diversity in tick-borne viruses and bacteria (Paddock et al., 2016). Increases in tick-borne pathogen transmission and case detection have garnered a great deal of attention, triggering greater funding, resources, and agency responses (CDC, 2018; Couzin-Frankel, 2019). Nevertheless, the expanding burden of tick-borne disease has also highlighted crucial gaps in knowledge, particularly with regards to geographic risk mapping, an area of great interest to public health agencies. This is particularly evident in the case of rickettsial pathogens, comprising the ehrlichiosis, anaplasmosis, and spotted fever rickettsioses, which compared to Lyme disease remain understudied. Rickettsial pathogens of medical importance may be encountered worldwide, and ticks from the *Amblyomma* genus are competent vectors for many of these pathogens (Levin et al., 2018).Although numbers of documented cases have been increasing in recent years, the true extent of geographic risk for rickettsial pathogens is challenging to delineate due to a lack of consistent,long-term, and widespread surveillance data, and regionally low case detection.

Species distributions models (SDMs), also commonly referred to as ecological niche models (ENMs), are becoming routinely used in vector-borne disease systems to model the potential geographic distribution of risk (e.g. Baak-Baak et al., 2017; Carvalho et al., 2015; Lippiet al., 2019; Peterson et al., 2002; Thomas and Beierkuhnlein, 2013). Broadly, this is accomplished by correlating locations where a species of interest is known to occur with the underlying environmental characteristics (e.g. climate, elevation, land cover). The resulting model can then be projected to unsampled areas on the landscape, providing a spatial prediction of areas that are ecologically suitable for species presence. In addition to predicting contemporary species distributions, SDMs are also employed to estimate the extent of potentially suitable habitat for invasive species, and potential shifts in geographic distributions due to climate change (Lippi et al., 2019). There are many methodological approaches to estimating species distributions, and some of the more commonly encountered approaches include Maximum Entropy (MaxEnt), Generalized Additive Models (GAM), Boosted Regression Trees(BRT), and Random Forests (RF) (Elith et al., 2008; Elith and Leathwick, 2009; Evans et al.,2011; Phillips et al., 2006). Although SDMs are a commonly used tool in estimating species ranges, the diversity in modeling approaches and applications makes it challenging to compare results across models.

For vector-borne disease SDMs, records of vector and/or pathogen presence (either the vector, the pathogen, the vector and pathogen, or even simply human case data) are often used as proxies for risk of exposure, and therefore transmission. Species distribution models have been used in a risk mapping capacity for many vector-borne disease systems, spanning a range of pathogens vectored by arthropods including mosquitoes, gnats, phlebotomine flies, fleas,triatomine bugs, and ticks (Crkvencic and Šlapeta, 2019). This framework is particularly useful in determining species limits for vectored transmission owing to the very close relationships between ectotherm life histories, pathogen replication, and environmental drivers such a temperature. Underlying distributions of reservoir hosts, another requisite component of zoonotic transmission cycles, are also determined by land cover and environmental conditions.

This work provides a comprehensive review of the published, peer-reviewed literature of studies that estimated species distribution, or ecological niche, of *Amblyomma* ticks, the rickettsial pathogens they vector, or their combined distributions. Following Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA) guidelines, we identified and compiled studies that used occurrence records and environmental predictors to estimate the geographic range of target organisms (Moher et al., 2009). Additionally, we provide a synthesis of current knowledge in the field, identifying the range of regions, spatial scales, and environmental determinants used to define risk in these systems. This work serves as a baseline for identifying knowledge gaps and guiding new studies of geographic risk mapping in understudied tick-borne disease systems.

## Materials and Methods

Literature searches were conducted following the guidelines in the PRISMA Statement, a checklist and flow diagram to ensure transparency and reproducibility in systematic reviews and meta-analyses (Liberati et al., 2009; Moher et al., 2009). Initial searches for peer-reviewed studies were conducted through September 2019. Five online databases were searched including Web of Science (Web of Science Core Collection, MEDLINE, BIOSIS Citation Index,Zoological Record) and Google Scholar. Searches were performed with combinations of key terms including “Amblyomma”, “Rickettsia*”, “niche model”, “ecological niche model”, and “species distribution model”. No restrictions were placed on geographic region of study or date of publication. Additional novel records for screening were identified via literature cited sections in records identified via database searches.

Duplicate records from the initial database searches were removed, and the remaining abstracts were screened for relevance. Records were excluded in this first screening based on publication type (i.e. literature reviews, opinion pieces, and synthesis papers were excluded), methodology (i.e. studies that did not examine geographic distributions or risk were excluded), and target organisms (i.e. vectors and diseases other than the *Amblyomma* spp. and rickettsial pathogens (*Ehrlichia* spp. and *Rickettsia* spp.) were excluded). *Anaplasma* spp. were excluded from our analysis as these are primarily transmitted by *Ixodes* ticks.

The remaining articles were then assessed in full for eligibility, where information on vectors, pathogens, modeling methods, geographic region, geographic scale of study, time period of study, data sources, and major findings were extracted from the text. Full-text articles in this final screening were flagged for possible exclusion based on focus the study (e.g. modeling tick distributions as an exercise to compare modeling logistics and methodologies), mismatch in target organisms (e.g. the distribution of a target pathogen in a vector of a non-target genus), or status as grey literature (e.g. unpublished theses).

## Results

### Reviewed Studies

The initial search of the peer-reviewed literature yielded 174 studies published between 1994 and 2019. After removal of duplicate search results, 109 unique studies remained. Studies were screened for topic relevance and abstract content, after which 33 studies remained. These publications were then assessed in full, resulting in 24 studies on species distribution models for *Amblyomma* ticks and/or rickettsial pathogens that were included in the literature synthesis (Table 1). The PRISMA flow diagram outlining our literature search and screening process is shown in Fig. 1.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/04/10/2020.04.07.20057083/T1)

Table 1. 
List of final publications on *Amblyomma* ticks and rickettsial group pathogens featured in literature review.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/10/2020.04.07.20057083/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/04/10/2020.04.07.20057083/F1)

Figure 1. PRISMA flow diagram outlining the literature search and screening process.

The geographic extent of suitability models varied with study focus and stated research goals, ranging from local and regional foci (e.g. county and state-level) to national and global species distributions. Studies were primarily conducted for the United States (46%), often limited to a single state or regional boundary. Other common geographic foci included Africa (29%) and Latin America (21%). The majority of studies (83%) reported the spatial scales of predicted distributions, which varied considerably across studies ranging from fine-scale (1 km) to coarse resolution (50 km) gridded models.

### Species Occurrence Data

The majority of reviewed studies (79%) modeled geographic distributions only for *Amblyomma* ticks, using vector presence as a proxy for pathogen transmission and disease risk. *Amblyomma americanum* was featured in 33% of papers, making it the most commonly studied vector, followed by *Amblyomma variegatum* (29%) and *Amblyomma hebraeum* (25%). Studies that estimated pathogen distributions, using health department data or wildlife blood samples to determine presence, accounted for 20% of the reviewed literature. Only 8% of studies used both tick occurrence and pathogen presence to model geographic risk of transmission.

A majority of studies (63%) obtained positive records of species occurrence from previously published literature. The time period of sample collections from previously published sources typically ranged from the 1950s through the early 2000s, though one study incorporated historical records dating back to the 1900s. In contrast, 29% of studies primarily obtained georeferenced data points from public databases or entomological collections. While using pre-existing databases of tick records yield higher occurrence frequencies that span greater periods of time, few details are typically provided regarding the nature of sample collections (e.g. active versus passive surveillance, transects versus convenience sampling, etc.). Only 21% of studies collected occurrence records solely through field sampling, with field collections spanning one to thirteen years. Few studies (13%) used true absence data collected via field sampling, instead opting to use modeling approaches that take advantage of presence-only datasets.

### Environmental Data

Environmental predictor datasets used to build SDMs were generally chosen in accordance with the specified geographic extent, scale, and goals of a given study. Many of the data products used for localized studies are only available for a given region or country (e.g. Daymet climate data for North America, USGS National Land Cover Database with coverage for the United States, etc.). Despite the wide range of environmental inputs across studies, the WorldClim dataset of long-term climate averages, and derived bioclimatic variables, were the most commonly utilized source of climatological data, featured as input data in 36% of reviewed papers (Hijmans et al., 2005). Six papers estimated potential shifts in vector ranges driven by future climate change and required modeled climate data at given time horizons, with half of these studies using WorldClim scenarios of future climatic conditions.

### Modeling Approaches and Output

The reviewed literature primarily consisted of studies that used presence-only, correlative modeling approaches. A variety of SDM methods were used to estimate tick distributions, including logistic regression (LR), geographically weighted regression (GWR), ecological niche factor analysis (ENFA), and generalized additive models (GAM) (Table 1). However, the MaxEnt algorithm was the most commonly used method, with 50% of studies using MaxEnt to estimate species distributions. Regression analyses were also frequently included in SDM studies (42%), even when more advanced statistical models were also used.

Environmental factors that were most influential in published SDMs were reported in 80% of reviewed studies. Generally, studies included some measure of temperature, precipitation, soil moisture, or land cover in models of tick distributions. The temporal scale of environmental predictors varied considerably between studies, ranging from daily temperature estimates and monthly ranges to annual and long-term climate averages. Covariates that contributed most prominently to distribution estimates were often reported in the reviewed literature, yet actual values and ranges for environmental predictors were seldom reported. Indicators of seasonality for precipitation (66% of studies) and temperature (50% of studies) were among the most consistent covariates included in final SDMs. The Normalized Difference Vegetation Index (NDVI), an index derived from remote sensing data to measure green vegetation cover, was used as an environmental predictor in 25% of studies. Despite the prevalence of NDVI in these studies, compared to climate variables there was little consensus regarding the reliability of this predictor in defining vector niches across studies. Only one study incorporated tick host density as an environmental predictor of habitat suitability.

## Discussion

Species distribution modeling has become a widely used tool for estimating ranges of organisms, or in the case of pathogens and their vectors, the geographic risk of disease exposure. However, the potential geographic distributions of rickettsial pathogens transmitted by *Amblyomma* spp. are still relatively understudied, compared to other vector-borne disease systems. In contrast with the 174 candidate publications identified for screening in this literature review, similar search terminology applied to other vector-borne disease systems yielded raw publication counts of 1,126 for *Aedes* spp. and dengue fever, 728 for *Anopheles* spp. and malaria, and 366 for *Ixodes* spp. and Lyme disease. This lag is likely in part to more recent recognition of the presence of these disease transmission systems for the spotted fever rickettsioses and ehrlichiosis tick-borne diseases (Childs and Paddock, 2003). Correlative models, such as SDMs, are useful in estimating potential species limits, particularly when data for mechanistic or process-driven models are lacking. This is the case for many tick species, where physiological mechanisms and limits are often poorly understood. For *Amblyomma* spp. only a few papers address the effects of temperature and humidity in a controlled setting (Guglielmone, 1992; Koch, 1983). At extreme temperatures where insect vectors may die, ticks will become quiescent until conditions are more favorable. Thresholds that are known to cause instantaneous tick mortality are prohibitive to life such as -22°C (Burks et al., 1996). Presence-only modeling approaches, which predominated in the reviewed literature, are also convenient when studying organisms that are difficult to extensively sample throughout their range. These methods are attractive in that they allow us to take advantage of pre-existing collection datasets, seemingly obviating the need for labor and resource intensive field sampling. However, with the exception of studies which explicitly conducted surveys for ticks, the full methods for originally obtaining presence points (e.g. the original collection strategies) are not always clearly defined. Biases introduced via sampling protocol (e.g. convenience sampling, or targeting a single life stage or behavior) may not adequately represent the true realized niche for species, dramatically influencing SDM predictions. In many tick models, data are collected through dragging or flagging which only samples questing ticks that have not found hosts. The number of successful host-seeking ticks is not known and for ticks that are not possible to collect on drags/flags the surveillance method may drastically underrepresent this life stage (Gaff et al., 2020). If the purpose of the model is to measure disease risk, questing data may be appropriate, but for tick control a greater understanding of the species life history would be needed. While these limitations may be logistically unavoidable, we recommend more detailed reporting of sampling methods and their associated limitations in future studies.

A further difficulty presents itself when estimating species distributions for ticks based on environmental conditions. Other arthropod vectors, such as mosquitoes, are extremely sensitive to fluctuations in climatic conditions, which in turn dictate the suitability of an area for survival and reproduction. In many instances, the physiological responses and environmental limits of insect vectors are well understood, via both laboratory experiments and empirical field studies (Mordecai et al., 2019; Paaijmans et al., 2013; Reuss et al., 2018). In contrast with insect vectors, ticks are resilient to many of the climatic factors that would limit other species. In other words, broad-scale patterns in temperature and precipitation are not necessarily primary drivers of tick presence on the landscape. This may contribute to the relatively low agreement in niche-defining environmental parameters across the reviewed studies. With the exception of indicators of seasonality, the major climatic and land cover predictors in the literature vary greatly with species and geographic extent. Given the close association between ticks and their vertebrate hosts, this may indicate that it is not the vector’s niche that is being modeled, but rather the niche of the host organisms that support tick populations. Furthermore, reaching consensus across SDMs is notoriously problematic, owing largely to the abundance of methodological approaches and lack of standardized reporting practices in presence data and final models (Carlson et al., 2018; Hao et al., 2019; Merow et al., 2013; Mordecai et al., 2019; Rund et al., 2019). We find similar issues when comparing published SDMs for *Amblyomma* ticks and rickettsial pathogens, where there is considerable diversity in methods and primary findings despite the small number of studies performed. While major environmental predictors are typically reported for SDMs, most studies do not report values or numerical ranges for suitability.

## Conclusions

Species distribution modeling of *Amblyomma* ticks and the rickettsial group pathogens they vector is underrepresented in the literature compared to other vector-borne disease systems. Even among a limited number of published studies, there is considerable variation in the methods and reported environmental influences for these models. This systematic literature review highlights a knowledge gap in our understanding of potential geographic risk for this transmission system. Given the recent public health interest in tick-borne diseases, the dearth of studies may result from lags in new data acquisition and limitations in our knowledge of the tick-pathogen ecology involved.

## Data Availability

All data are within the manuscript

## Funding

CAL, HDG, and SJR were funded by NIH 1R01AI136035-01. ALW and SJR were additionally funded by CDC grant 1U01CK000510-01: Southeastern Regional Center of Excellence in Vector-Borne Diseases: The Gateway Program. This publication was supported by the Cooperative Agreement Number above from the Centers for Disease Control and Prevention. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of the Centers for Disease Control and Prevention.

## Acknowledgements

## Abbreviations

SDM
:   species distribution model
ENM
:   ecological niche model
PRISMA
:   Preferred Reporting Items for Systematic Reviews and Meta-analyses
GAM
:   generalized additive model
BRT
:   boosted regression trees
RF
:   random forests
LR
:   logistic regression (LR)
GWR
:   geographically weighted regression
ENFA
:   ecological niche factor analysis

*   Received April 7, 2020.
*   Revision received April 7, 2020.
*   Accepted April 10, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## Literature Cited

1.  Acevedo-Gutierrez, L.Y., Paternina, L.E., Londoño, A.F., Parra-Henao, G., Rodas, J.D., 2018. Potential models of the geographic and climatic distribution of the Amblyomma cajennense complex (Acari: Ixodidae), a potential vector of Rickettsia rickettsii in Colombia. Biomédica 38.
    
    

2.  Baak-Baak, C.M., Moo-Llanes, D.A., Cigarroa–Toledo, N., Puerto, F.I., Machain-Williams, C., Reyes-Solis, G., Nakazawa, Y.J., Ulloa-Garcia, A., Garcia-Rejon, J.E., 2017. Ecological Niche Model for Predicting Distribution of Disease-Vector Mosquitoes in Yucatán State, México. J. Med. Entomol. 54, 854–861. [https://doi.org/10.1093/jme/tjw243](https://doi.org/10.1093/jme/tjw243)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jme/tjw243&link_type=DOI) 

3.  Behravesh, C.B., Paddock, C.D., Dahlgren, F.S., Eisen, R.J., Springer, Y.P., 2016. Expanding Range of Amblyomma americanum and Simultaneous Changes in the Epidemiology of Spotted Fever Group Rickettsiosis in the United States. Am. J. Trop. Med. Hyg. 94, 35–42. [https://doi.org/10.4269/ajtmh.15-0580](https://doi.org/10.4269/ajtmh.15-0580)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo3OiI5NC8xLzM1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMTAvMjAyMC4wNC4wNy4yMDA1NzA4My5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

4.  Bermúdez, S.E., Castro, A.M., Trejos, D., García, G.G., Gabster, A., Miranda, R.J., Zaldívar, Y., Paternina, L.E., 2016. Distribution of Spotted Fever Group Rickettsiae in Hard Ticks (Ixodida: Ixodidae) from Panamanian Urban and Rural Environments (2007–2013). EcoHealth 13, 274–284. [https://doi.org/10.1007/s10393-016-1118-8](https://doi.org/10.1007/s10393-016-1118-8)
    
    

5.  Burks, C.S., Stewart Jr, R.L., Needham, G.R., Lee Jr, R.E., 1996. The role of direct chilling injury and inoculative freezing in cold tolerance of Amblyomma americanum, Dermacentor variabilis and Ixodes scapularis. Physiol. Entomol. 21, 44–50.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-3032.1996.tb00833.x&link_type=DOI) 

6.  Carlson, C.J., Dougherty, E., Boots, M., Getz, W., Ryan, S.J., 2018. Consensus and conflict among ecological forecasts of Zika virus outbreaks in the United States. Sci. Rep. 8. [https://doi.org/10.1038/s41598-018-22989-0](https://doi.org/10.1038/s41598-018-22989-0)
    
    

7.  Carvalho, B.M., Rangel, E.F., Ready, P.D., Vale, M.M., 2015. Ecological Niche Modelling Predicts Southward Expansion of Lutzomyia (Nyssomyia) flaviscutellata (Diptera: Psychodidae: Phlebotominae), Vector of Leishmania (Leishmania) amazonensis in South America, under Climate Change. PLOS ONE 10, e0143282. [https://doi.org/10.1371/journal.pone.0143282](https://doi.org/10.1371/journal.pone.0143282)
    
    

8.  CDC, 2018. Illnesses on the rise from mosquito, tick, and flea bites.
    
    

9.  Childs, J.E., Paddock, C.D., 2003. The ascendancy of Amblyomma americanum as a vector of pathogens affecting humans in the United States. Annu. Rev. Entomol. 48, 307–337.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.ento.48.091801.112728&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12414740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180709200015&link_type=ISI) 

10. Couzin-Frankel, J., 2019. Long-underfunded Lyme disease research gets an injection of money—and ideas. Science. [https://doi.org/10.1126/science.aax7340](https://doi.org/10.1126/science.aax7340)
    
    

11. Crkvencic, N., Šlapeta, J., 2019. Climate change models predict southerly shift of the cat flea (Ctenocephalides felis) distribution in Australia. Parasit. Vectors 12. [https://doi.org/10.1186/s13071-019-3399-6](https://doi.org/10.1186/s13071-019-3399-6)
    
    

12. Cumming, G.S., 2002. COMPARING CLIMATE AND VEGETATION AS LIMITING FACTORS FOR SPECIES RANGES OF AFRICAN TICKS. Ecology 83, 255–268. [https://doi.org/10.1890/0012-9658(2002)083[0255:CCAVAL]2.0.CO;2](https://doi.org/10.1890/0012-9658(2002)083[0255:CCAVAL]2.0.CO;2)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1890/0012-9658(2002)083[0255:CCAVAL]2.0.CO;2&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173117700022&link_type=ISI) 

13. Cumming, G.S., 2000. Using habitat models to map diversity: pan-African species richness of ticks (Acari: Ixodida). J. Biogeogr. 27, 425–440. [https://doi.org/10.1046/j.1365-2699.2000.00419.x](https://doi.org/10.1046/j.1365-2699.2000.00419.x)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2699.2000.00419.x&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000089193500016&link_type=ISI) 

14. Cumming, G.S., Van Vuuren, D.P., 2006. Will climate change affect ectoparasite species ranges? Glob. Ecol. Biogeogr. 15, 486–497. [https://doi.org/10.1111/j.1466-822X.2006.00241.x](https://doi.org/10.1111/j.1466-822X.2006.00241.x)
    
    

15. Elith, J., Leathwick, J.R., 2009. Species Distribution Models: Ecological Explanation and Prediction Across Space and Time. Annu. Rev. Ecol. Evol. Syst. 40, 677–697. [https://doi.org/10.1146/annurev.ecolsys.110308.120159](https://doi.org/10.1146/annurev.ecolsys.110308.120159)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.ecolsys.110308.120159&link_type=DOI) 

16. Elith, J., Leathwick, J.R., Hastie, T., 2008. A working guide to boosted regression trees. J. Anim. Ecol. 77, 802–813. [https://doi.org/10.1111/j.1365-2656.2008.01390.x](https://doi.org/10.1111/j.1365-2656.2008.01390.x)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2656.2008.01390.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18397250&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256539800020&link_type=ISI) 

17. Estrada-Peña, A., 2003. Climate change decreases habitat suitability for some tick species (Acari: Ixodidae) in South Africa. Onderstepoort J. Vet. Res. 70, 79–93.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12967169&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

18. Estrada-Peña, A., Horak, I.G., Petney, T., 2008. Climate changes and suitability for the ticks Amblyomma hebraeum and Amblyomma variegatum (Ixodidae) in Zimbabwe (1974– 1999). Vet. Parasitol. 151, 256–267. [https://doi.org/10.1016/j.vetpar.2007.11.014](https://doi.org/10.1016/j.vetpar.2007.11.014)
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18083309&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

19. Estrada-Peña, A., Pegram, R.G., Barré, N., Venzal, J.M., 2007. Using invaded range data to model the climate suitability for Amblyomma variegatum (Acari: Ixodidae) in the New World. Exp. Appl. Acarol. 41, 203–214. [https://doi.org/10.1007/s10493-007-9050-9](https://doi.org/10.1007/s10493-007-9050-9)
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17347921&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

20. Estrada-Peña, A., Tarragona, E.L., Vesco, U., Meneghi, D. de, Mastropaolo, M., Mangold, A.J., Guglielmone, A.A., Nava, S., 2014. Divergent environmental preferences and areas of sympatry of tick species in the Amblyomma cajennense complex (Ixodidae). Int. J. Parasitol. 44, 1081–1089. [https://doi.org/10.1016/j.ijpara.2014.08.007](https://doi.org/10.1016/j.ijpara.2014.08.007)
    
    

21. 1.  Drew, C.A., 
    2.  Wiersma, Y.F., 
    3.  Huettmann, F.
    
    Evans, J.S., Murphy, M.A., Holden, Z.A., Cushman, S.A., 2011. Modeling Species Distribution and Change Using Random Forest, in: Drew, C.A., Wiersma, Y.F., Huettmann, F. (Eds.), Predictive Species and Habitat Modeling in Landscape Ecology. Springer New York, New York, NY, pp. 139–159. [https://doi.org/10.1007/978-1-4419-7390-0_8](https://doi.org/10.1007/978-1-4419-7390-0_8)
    
    

22. Gaff, H., Eisen, R.J., Eisen, L., Nadolny, R., Bjork, J., Monaghan, A.J., 2020. LYMESIM 2.0: An Updated Simulation of Blacklegged Tick (Acari: Ixodidae) Population Dynamics and Enzootic Transmission of Borrelia burgdorferi (Spirochaetales: Spirochaetaceae). J. Med. Entomol. [https://doi.org/10.1093/jme/tjz252](https://doi.org/10.1093/jme/tjz252)
    
    

23. Guglielmone, A.A., 1992. The effect of temperature and humidity on development and longevity of Amblyomma triguttatum triguttatum (Acarina: Ixodidae). Bull. Entomol. Res. 82, 203–208. [https://doi.org/10.1017/S0007485300051737](https://doi.org/10.1017/S0007485300051737)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0007485300051737&link_type=DOI) 

24. Hao, T., Elith, J., Guillera□Arroita, G., Lahoz□Monfort, J.J., 2019. A review of evidence about use and performance of species distribution modelling ensembles like BIOMOD. Divers. Distrib. 25, 839–852. [https://doi.org/10.1111/ddi.12892](https://doi.org/10.1111/ddi.12892)
    
    

25. Hijmans, R.J., Cameron, S.E., Parra, J.L., Jones, P.G., Jarvis, A., 2005. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. [https://doi.org/10.1002/joc.1276](https://doi.org/10.1002/joc.1276)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/joc.1276&link_type=DOI) 

26. Illoldi-Rangel, P., Rivaldi, C.-L., Sissel, B., Trout Fryxell, R., Gordillo-Pérez, G., Rodríguez-Moreno, A., Williamson, P., Montiel-Parra, G., Sánchez-Cordero, V., Sarkar, S., 2012. Species Distribution Models and Ecological Suitability Analysis for Potential Tick Vectors of Lyme Disease in Mexico. J. Trop. Med. 2012, 1–10. [https://doi.org/10.1155/2012/959101](https://doi.org/10.1155/2012/959101)
    
    

27. Kessler, William H, Blackburn, J.K., Sayler, K.A., Glass, G.E., 2019. Estimating the Geographic Distribution of Host-Seeking Adult Amblyomma americanum (Acari: Ixodidae) in Florida. J. Med. Entomol. 56, 55–64. [https://doi.org/10.1093/jme/tjy147](https://doi.org/10.1093/jme/tjy147)
    
    

28. Kessler, William H., Ganser, C., Glass, G.E., 2019. Modeling the Distribution of Medically Important Tick Species in Florida. Insects 10, 190. [https://doi.org/10.3390/insects10070190](https://doi.org/10.3390/insects10070190)
    
    

29. Koch, H.G., 1983. Lone star tick: oviposition, egg hatch, and molting under naturally fluctuating ambient temperatures and humidities in the field. Southwest. Entomol. 6, 240–244.
    
    

30. Levin, M.L., Schumacher, L.B.M., Snellgrove, A., 2018. Effects of Rickettsia amblyommatis Infection on the Vector Competence of Amblyomma americanum Ticks for Rickettsia rickettsii. Vector Borne Zoonotic Dis. Larchmt. N 18, 579–587. [https://doi.org/10.1089/vbz.2018.2284](https://doi.org/10.1089/vbz.2018.2284)
    
    

31. Liberati, A., Altman, D.G., Tetzlaff, J., Mulrow, C., Gøtzsche, P.C., Ioannidis, J.P.A., Clarke, M., Devereaux, P.J., Kleijnen, J., Moher, D., 2009. The PRISMA Statement for Reporting Systematic Reviews and Meta-Analyses of Studies That Evaluate Health Care Interventions: Explanation and Elaboration. PLoS Med. 6, e1000100. [https://doi.org/10.1371/journal.pmed.1000100](https://doi.org/10.1371/journal.pmed.1000100)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1000100&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19621070&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

32. Lippi, C.A., Stewart-Ibarra, A.M., Loor, M.E.F.B., Zambrano, J.E.D., Lopez, N.A.E., Blackburn, J.K., Ryan, S.J., 2019. Geographic shifts in Aedes aegypti habitat suitability in Ecuador using larval surveillance data and ecological niche modeling: Implications of climate change for public health vector control. PLoS Negl. Trop. Dis. 13, e0007322. [https://doi.org/10.1371/journal.pntd.0007322](https://doi.org/10.1371/journal.pntd.0007322)
    
    

33. Lynen, G., Zeman, P., Bakuname, C., Di Giulio, G., Mtui, P., Sanka, P., Jongejan, F., 2007. Cattle ticks of the genera Rhipicephalus and Amblyomma of economic importance in Tanzania: distribution assessed with GIS based on an extensive field survey. Exp. Appl. Acarol. 43, 303–319. [https://doi.org/10.1007/s10493-007-9123-9](https://doi.org/10.1007/s10493-007-9123-9)
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18044004&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

34. Manangan, J.S., Schweitzer, S.H., Nibbelink, N., Yabsley, M.J., Gibbs, S.E.J., Wimberly, M.C., 2007. Habitat Factors Influencing Distributions of Anaplasma phagocytophilum and Ehrlichia chaffeensis in the Mississippi Alluvial Valley. Vector-Borne Zoonotic Dis. 7, 563–574. [https://doi.org/10.1089/vbz.2007.0116](https://doi.org/10.1089/vbz.2007.0116)
    
    

35. Merow, C., Smith, M.J., Silander, J.A., 2013. A practical guide to MaxEnt for modeling species’ distributions: what it does, and why inputs and settings matter. Ecography 36, 1058–1069. [https://doi.org/10.1111/j.1600-0587.2013.07872.x](https://doi.org/10.1111/j.1600-0587.2013.07872.x)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0587.2013.07872.x&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325114500002&link_type=ISI) 

36. Moher, D., Liberati, A., Tetzlaff, J., Altman, D.G., The PRISMA Group, 2009. Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. PLoS Med. 6, e1000097. [https://doi.org/10.1371/journal.pmed.1000097](https://doi.org/10.1371/journal.pmed.1000097)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1000097&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19621072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

37. Mordecai, E.A., Caldwell, J.M., Grossman, M.K., Lippi, C.A., Johnson, L.R., Neira, M., Rohr, J.R., Ryan, S.J., Savage, V., Shocket, M.S., Sippy, R., Stewart Ibarra, A.M., Thomas, M.B., Villena, O., 2019. Thermal biology of mosquito□borne disease. Ecol. Lett. 22, 1690–1708. [https://doi.org/10.1111/ele.13335](https://doi.org/10.1111/ele.13335)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ele.13335&link_type=DOI) 

38. Norval, R.A.I., Perry, B.D., Meltzer, M.I., Kruska, R.L., Booth, T.H., 1994. Factors affecting the distributions of the ticks Amblyomma hebraeum and A. variegatum in Zimbabwe: implications of reduced acaricide usage. Exp. Appl. Acarol. 18, 383–407. [https://doi.org/10.1007/BF00051522](https://doi.org/10.1007/BF00051522)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00051522&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7628253&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 

39. Oliveira, S.V. de, Romero-Alvarez, D., Martins, T.F., Santos, J.P. dos, Labruna, M.B., Gazeta, G.S., Escobar, L.E., Gurgel-Gonçalves, R., 2017. Amblyomma ticks and future climate: Range contraction due to climate warming. Acta Trop. 176, 340–348. [https://doi.org/10.1016/j.actatropica.2017.07.033](https://doi.org/10.1016/j.actatropica.2017.07.033)
    
    

40. Paaijmans, K.P., Heinig, R.L., Seliga, R.A., Blanford, J.I., Blanford, S., Murdock, C.C., Thomas, M.B., 2013. Temperature variation makes ectotherms more sensitive to climate change. Glob. Change Biol. 19, 2373–2380. [https://doi.org/10.1111/gcb.12240](https://doi.org/10.1111/gcb.12240)
    
    

41. Paddock, C.D., Lane, R.S., Staples, J.E., Labruna, M.B., 2016. Changing paradigms for tick-borne diseases in the Americas, in: Mack, A. (Ed.), Global Health Impacts of Vector-Borne Diseases: Workshop Summary. National Academies of Sciences, Engineering, and Medicine, Washington, D.C. [https://doi.org/10.17226/21792](https://doi.org/10.17226/21792)
    
    

42. Pascoe, E.L., Marcantonio, M., Caminade, C., Foley, J.E., 2019. Modeling Potential Habitat for Amblyomma Tick Species in California. Insects 10, 201. [https://doi.org/10.3390/insects10070201](https://doi.org/10.3390/insects10070201)
    
    

43. Peterson, A.T., Sánchez-Cordero, V., Ben Beard, C., Ramsey, J.M., 2002. Ecologic Niche Modeling and Potential Reservoirs for Chagas Disease, Mexico. Emerg. Infect. Dis. 8, 662–667. [https://doi.org/10.3201/eid0807.010454](https://doi.org/10.3201/eid0807.010454)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid0807.010454&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12095431&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F10%2F2020.04.07.20057083.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000176394200004&link_type=ISI) 

44. Phillips, S.J., Anderson, R.P., Schapire, R.E., 2006. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190, 231–259. [https://doi.org/10.1016/j.ecolmodel.2005.03.026](https://doi.org/10.1016/j.ecolmodel.2005.03.026)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ecolmodel.2005.03.026&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233859600001&link_type=ISI) 

45. Raghavan, R.K., Goodin, D.G., Hanzlicek, G.A., Zolnerowich, G., Dryden, M.W., Anderson, G.A., Ganta, R.R., 2016a. Maximum Entropy-Based Ecological Niche Model and Bio-Climatic Determinants of Lone Star Tick (Amblyomma americanum) Niche. Vector-Borne Zoonotic Dis. 16, 205–211. [https://doi.org/10.1089/vbz.2015.1837](https://doi.org/10.1089/vbz.2015.1837)
    
    

46. Raghavan, R.K., Goodin, D.G., Neises, D., Anderson, G.A., Ganta, R.R., 2016b. Hierarchical Bayesian Spatio–Temporal Analysis of Climatic and Socio–Economic Determinants of Rocky Mountain Spotted Fever. PLOS ONE 11, e0150180. [https://doi.org/10.1371/journal.pone.0150180](https://doi.org/10.1371/journal.pone.0150180)
    
    

47. Raghavan, R.K., Peterson, A.T., Cobos, M.E., Ganta, R., Foley, D., 2019. Current and Future Distribution of the Lone Star Tick, Amblyomma americanum (L.) (Acari: Ixodidae) in North America. PLOS ONE 14, e0209082. [https://doi.org/10.1371/journal.pone.0209082](https://doi.org/10.1371/journal.pone.0209082)
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0209082&link_type=DOI) 

48. Reese, S.M., Yates, K.F., Mead, P.S., Nicholson, W.L., Blount, K., Graham, C.B., Brown, H.E., Helterbrand, W.S., Dietrich, G., Patrick, S.L., MacMillan, K., Eisen, R.J., 2011. An Acarologic Survey and Amblyomma americanum Distribution Map with Implications for Tularemia Risk in Missouri. Am. J. Trop. Med. Hyg. 84, 411–419. [https://doi.org/10.4269/ajtmh.2011.10-0593](https://doi.org/10.4269/ajtmh.2011.10-0593)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI4NC8zLzQxMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzEwLzIwMjAuMDQuMDcuMjAwNTcwODMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

49. Reuss, F., Wieser, A., Niamir, A., Bálint, M., Kuch, U., Pfenninger, M., Müller, R., 2018. Thermal experiments with the Asian bush mosquito (Aedes japonicus japonicus) (Diptera: Culicidae) and implications for its distribution in Germany. Parasit. Vectors 11. [https://doi.org/10.1186/s13071-018-2659-1](https://doi.org/10.1186/s13071-018-2659-1)
    
    

50. Rund, S.S.C., Braak, K., Cator, L., Copas, K., Emrich, S.J., Giraldo-Calderón, G.I., Johansson, M.A., Heydari, N., Hobern, D., Kelly, S.A., Lawson, D., Lord, C., MacCallum, R.M., Roche, D.G., Ryan, S.J., Schigel, D., Vandegrift, K., Watts, M., Zaspel, J.M., Pawar, S., 2019. MIReAD, a minimum information standard for reporting arthropod abundance data. Sci. Data 6. [https://doi.org/10.1038/s41597-019-0042-5](https://doi.org/10.1038/s41597-019-0042-5)
    
    

51. Springer, Y.P., Jarnevich, C.S., Monaghan, A.J., Eisen, R.J., Barnett, D.T., 2015. Modeling the Present and Future Geographic Distribution of the Lone Star Tick, Amblyomma americanum (Ixodida: Ixodidae), in the Continental United States. Am. J. Trop. Med.Hyg. 93, 875–890. [https://doi.org/10.4269/ajtmh.15-0330](https://doi.org/10.4269/ajtmh.15-0330)
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI5My80Lzg3NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzEwLzIwMjAuMDQuMDcuMjAwNTcwODMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

52. Thomas, S.M., Beierkuhnlein, C., 2013. Predicting ectotherm disease vector spread—benefits from multidisciplinary approaches and directions forward. Naturwissenschaften 100, 395–405. [https://doi.org/10.1007/s00114-013-1039-0](https://doi.org/10.1007/s00114-013-1039-0)
    
    

53. Wimberly, M.C., Yabsley, M.J., Baer, A.D., Dugan, V.G., Davidson, W.R., 2008. Spatial heterogeneity of climate and land-cover constraints on distributions of tick-borne pathogens. Glob. Ecol. Biogeogr. 17, 189–202. [https://doi.org/10.1111/j.1466-8238.2007.00353.x](https://doi.org/10.1111/j.1466-8238.2007.00353.x)