Next Article in Journal
Prediction of the Hind-Leg Muscles Weight of Yearling Dairy-Beef Steers Using Carcass Weight, Wither Height and Ultrasound Carcass Measurements
Previous Article in Journal
Echocardiographic Findings in Canine Model of Chagas Disease Immunized with DNA Trypanosoma cruzi Genes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 10
Issue 4
10.3390/ani10040649
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ecological Niche Models of Four Hard Tick Genera (Ixodidae) in Mexico

by
Emilio Clarke-Crespo
1,2,
Claudia N. Moreno-Arzate
3 and
Carlos A. López-González
2,*
1
Tecnológico de Monterrey, Escuela de Ingeniería y Ciencias, Centro de Bioingeniería, Queretaro 76130, Mexico
2
Facultad de Ciencias Naturales, Universidad Autónoma de Querétaro, Queretaro 76230, Mexico
3
Instituto de Ecología, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Animals 2020, 10(4), 649; https://doi.org/10.3390/ani10040649
Submission received: 11 February 2020 / Revised: 7 April 2020 / Accepted: 7 April 2020 / Published: 9 April 2020
(This article belongs to the Section Ecology and Conservation)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Vector-borne diseases currently represent a significant threat to public health, mainly due to the changes that humans are producing in ecosystems and climates. Analyzing the environmental conditions that allow the establishment and survival of ticks could help determine possible sites for the appearance of infectious outbreaks. In this study, nine ecological niche models were generated from different algorithms to determine the current potential distribution of four tick genera in Mexico. Temperature and moisture have been considered as the main factors limiting tick distribution. However, the analysis of the ecological niche models determined that the four genera exhibited different distribution patterns, which may be associated with their physiological and ecological differences. This type of analysis can improve our understanding of the dynamics of ticks and, therefore, can be very useful in monitoring programs of the diseases they transmit.

Abstract

Ticks are vectors of a large number of pathogens of medical and veterinary importance, and in recent years, they have participated in the rise of multiple infectious outbreaks around the world. Studies have proposed that temperature and precipitation are the main variables that limit the geographical distribution of ticks. The analysis of environmental constraints with ecological niche modeling (ENM) techniques can improve our ability to identify suitable areas for emergence events. Algorithms used in this study showed different distributional patterns for each tick genera; the environmental suitability for Amblyomma includes warm and humid localities below 1000 m above the sea level, while Ixodes is mainly associated with ecosystems with high vegetation cover. Dermacentor and Rhipicephalus genus presented wider distribution patterns; the first includes species that are well adapted to resist desiccation, whereas the latter includes generalist species that are mostly associated with domestic hosts in Mexico. Ecological niche models have proven to be useful in estimating the geographic distribution of many taxa of ticks. Despite our limited knowledge of tick’s diversity, ENM can improve our understanding of the dynamics of vector-borne diseases and can assist public health decision-making processes.

1. Introduction

Vector-borne diseases are complex systems that require the interactions between arthropod vectors hosts and pathogens, constrained by a set of environmental variables [1]. Particularly, vector reproduction, dispersal, and survival are strongly influenced by environmental conditions [2]. Ecological niche models (ENMs) have become powerful tools in the study of vector ecology. These are mainly used to determine their current and/or future potential distribution, to identify risk areas prone to the emergence of infectious outbreaks, or to distinguish environmental variables that regulate the dynamics of infectious diseases at a landscape level [3]. When applied to the estimation of potential geographic distribution, ecological niche models identify nonrandom patterns that usually emerge from the superposition analysis of a collection of presence-only or presence/absence data with a set of environmental variables [4,5] using different machine learning techniques [6,7].
An important arthropod disease carrier is ticks, which, in addition to being hematophagous ectoparasites of wild and domesticated animals and humans, are also known to be key vectors of a great variety of pathogens, such as protozoa, rickettsia, spirochaetes, and viruses [8,9,10]. Human-induced environmental changes appear to be important drivers that enhance tick distribution and survival. These vectors have triggered important human outbreaks such as tick-borne encephalitis (TBE) in Europe, Kyasanur forest disease (KFD) in India, Crimean-Congo hemorrhagic fever (CCHF) in Turkey and Russia, Q fever in the Netherlands, and Rocky Mountain spotted fever (RMSF) in the southern United States and in northern Mexico [11,12,13]. Overall, these events have supported the recognition of tick and tick-borne diseases as important emerging threats to humans and animals [14,15,16].
Hard ticks’ survival off hosts relies heavily on their ability to resist desiccation either by physiological or behavioral adaptations or by the environmental characteristics of the ecosystem in which they are found. For example, Illoldi-Rangel et al. [17] determined that the ecological suitability determined by climatic and environmental variables of ticks related to the transmission of Lyme disease of the genus Ixodes and Amblyomma cajennense can be used to assess the risk of vector-borne diseases even in poorly studied sites such as Mexico. Feria-Arrollo et al. [18] describe how the distribution of Ixodes scapularis responds positively to climate change, expanding its geographical distribution and therefore increasing the risk of contact between this important vector with humans and livestock.
Many vector-borne pathogens such as Lyme disease, malaria, Chagas disease, or Leishmaniosis simultaneously use several arthropod species to reach their ultimate hosts [19,20,21,22]. The assessment of ecological niche models at a genus level can be a useful tool to identify areas of importance for the dynamics of vector-borne diseases, for the identification of areas of higher risk of contact, or for the delimitation of sanitary control fences to prevent the emergence of a potential infectious outbreak. Strategies intended to contain outbreak episodes are usually expensive and difficult to implement [23], and modeling techniques can be very useful in the design of more efficient sentinel strategies.
In this study, we estimated the potential distribution of four genera of ticks of medical and veterinary importance in Mexico using nine modeling algorithms in order to identify areas of greater suitability for their presence.

2. Materials and Methods

2.1. Tick Occurrence Data

For this study, an exhaustive search was conducted of the presence data of the genera Amblyomma, Ixodes, Dermacentor, and Rhipicephalus in Mexico. This search was complicated by the controversial systematics and the difficult taxonomic diagnosis of the group as well as the general lack of study of the tick species of Mexico, with many species having a single record for the entire country. To reduce problems associated with misidentification, we cross-checked all records with the publications of the national collection of mites of the National Autonomous University of Mexico (CNAC-UNAM) [24,25,26]. It is important to note that many tick records in Mexico do not possess appropriate georeferencing, ambiguously indicating their location at the municipal or state level, data useless for the generation of species-specific ecological niche models. The last criterion that we used to refine the database for each genus of tick was its occurrence on a native host. This assumption eliminated many records that are associated with the movement and concentration of livestock in central feedlots collection centers, which could strongly over-fit the predicted distribution in the resulting models.
Each of the four data sets was calibrated and validated according to the training ratio proposed by Huberty [27], which consists of using 70% of the randomly selected records as the training set and the remaining 30% as the evaluation set of each model.

2.2. Pseudoabsence Data

Since real absences are difficult to obtain and various ENM need absences, we calculated random pseudoabsence records which is considered a robust method for the algorithms used [28]. We utilized Ecospat Package in R [29] with the following criteria: the proportion of one presence to 10 pseudoabsences to avoid extremely unbalanced samples of presence-absence [30], excluding points of presence, with a minimum distance of 1 km between them and the presence points, and within the calibration area [31].

2.3. Environmental Variables

We calibrated the ecological niche models of the four tick genera in Mexico using two sets of bioclimatic variables. The bioclimatic variables were downloaded from WorldClim (www.worldclim.org) at ~1 km2 spatial resolution [32]. Both climatic sets were chosen based on tick biology [33], and of these, highly correlated variables were removed from the analysis following a Pearson correlation test [34]. In addition to the climatic information, we also included two categorical variables that have been proposed as important environmental factors related to tick survival: the type of soil (www.inegi.org.mx/temas/edafologia/) and land use and vegetation (www.inegi.org.mx/temas/usosuelo/). We resampled both of these categorical layers to match the spatial resolution of the WoldClim layers.
The mean diurnal range, the annual precipitation, the precipitation of driest month, the precipitation seasonality, the type of soil, and the type of land use and vegetation were used to generate the ecological niche models for the four tick genera. To complement the aforementioned additional variables for the Amblyomma genus included the minimum temperature of the coldest month and the precipitation of the warmest quarter. For the three remaining genera, we added the maximum temperature of the warmest month, the mean temperature of the wettest quarter, and the precipitation of the wettest quarter (Table 1). The extent of the model calibration area had a strong effect on ENM results [35]. We selected the same polygon for all genera considering Mexico’s political borders. Although the four genera have a larger distribution than our surveyed area, we considered that expanding the models outside the area of interest would result in the loss of the fine-scale aspects of the analysis for Mexico.

2.4. Ecological Niche Modeling

The ecological niche models for the four tick genera were estimated with nine different algorithms using seven models included in SDM package in R software v. 3.5 [36], maximum entropy (MAXENT) with ENMEval package [37], and generalized additive model (GAM) with the GAM package [38]. GAM was calibrated out of the sdm package because the number of species-specific tick records considered was too small to predict successfully their distribution. The modeling algorithms developed were Bioclimatic envelope (BIOCLIM) [20], maximum entropy (MAXENT) [39,40], generalized linear models (GLM) [41,42], multivariate adaptive regression splines (MARS) [43], classification and regression trees (CART) [44], mixture discriminant analysis (MDA) [45], random forest (RF) [46], boosted regression trees (BRT) [47,48], and generalized additive model (GAM) [38].
BIOCLIM is a simple algorithm that describes the n-dimensional niche of a species in terms of a quadrangular climatic envelope defined by the range of occurrence values at each variable [49,50]. MAXENT is an algorithm that uses the maximum entropy principle and a Bayesian procedure to produce a probability surface where entropy is maximized to reflect the environmental suitability of the geographic area for the species [51]. GLM is an extension of linear models without force data into unnatural scales and related the mean of the response variable and their linear combination of the explanatory variables. GAM is a semi-parametric extension of GLM and the relationship between the mean of the response variable and a “smoothed” function of the explanatory variable(s) [52]. BRT is an ensemble method and combines regression trees and boosting [48]. RF is also a combination of tree predictors and is a classifier. MARS combines regression trees and spline fitting and describes nonlinear relationships between species and environmental variables [53].
We executed Bioclim in the DISMO package [54] and GLM, GAM, BRT, RF, and MARS in SDM package [36]. We ran MAXENT in the ENEval package using default parameters to select feature types and regulation multiplier based on changes in Akaike´s Information Criterion (AIC) [37]. All models were developed in environment R v. 3.5 [55].
The accuracy of each predicted model was evaluated with the area under the curve index (AUC) [56]. Despite the enormous discussion about its relevance as a discriminatory method in the species distribution modeling community, it remains as one of the most used criteria for their evaluation [57,58] because this metric has the ability to discriminate presence from absence (or background) and provides a single measure of performance. The AUC criterion ranges from 0.5 to 1. Values that vary between 0.9 and 1 imply that models have remarkable predictive accuracy, while values that vary between 0.8 and 0.9 have good accuracy and between 0.7 and 0.8 suggest that the models have regular performance [59]. AUC values below 0.7 indicate models of poor and/or failed predictive capacity and should be unconsidered. For all algorithms, we transformed each resulting model into a binary (presence-absence) map by selecting the threshold value where 10% of the occurrence records were left out. We decided to use this threshold criterion to avoid overprediction due to potentially erroneous occurrences [60].
Finally, we added the nine models generated for each genus of ticks in QGIS v. 3.6 to generate consensus maps that show the most suitable sites for these genera of ticks in Mexico.

3. Results

We found 85 total records for Amblyomma spp., 44 for Dermacentor spp., 71 for Ixodes spp., and 72 for Rhipicephalus spp. (Figure 1) (Table S1). We selected only the ticks recorded from native hosts because our main concern was the natural distribution of the four genera of ticks in Mexico. This criterion reduced the final number of occurrences because of numerous reports on livestock of uncertain geographic origin. The accuracy for ENM estimated that the models for the tick genera varied from very good (AUC values ≥ 0.8) to regular (AUC values between 0.7 and 0.8) for the nine algorithms used (Table 2). We removed Bioclim from the final model because of its low accuracy for all genera.
The models of distribution estimated for the genera Amblyomma and Rhipicephalus had the highest AUC values (Table 2). Overall, the algorithms used to estimate the potential distribution of Amblyomma spp. suggested that the areas with greater suitability in Mexico corresponded with the neotropical region of the country (Figure 2). The most important environmental variables for the genus Amblyomma were, in order of importance, land use and vegetation, minimum temperature of the coldest month, and precipitation of the warmest quarter (Table S1), which together accounted for about 85% of the registered presences. In contrast, the distribution models for Rhipicephalus spp. suggested that the areas of greatest suitability are found in the northeastern portion of the country, in the Yucatan peninsula, and in some areas around the sea of Cortez (Figure 2). The environmental variables that were strongly associated with the presence of the records of Rhipicephalus spp. were land use and vegetation, and annual precipitation (Table S2).
The models generated for the Ixodes and Dermacentor genera had a lower predictive precision. The areas of greatest suitability for the Ixodes genus in Mexico are mainly associated with forested areas in the central and southern mountain ranges (Neovolcanic axis and the Sierra Madre Oriental and Sur) and the coastal areas of the Yucatan peninsula (Figure 2). The environmental variables that were associated with the presence of almost 90% of the records are land use, annual rainfall, and average daytime range (Table S2). Finally, the areas suitable for the genus Dermacentor in Mexico are more widely distributed than the other three genera (Figure 2). However, the areas of greatest suitability for this genus were mainly located in the Nearctic portion of the country.

4. Discussion

The ecological niche models indicated that the four tick genera have different distribution patterns in Mexico (Figure 2). The genus Amblyomma in Mexico is represented by 26 species, seven having a Neotropical distribution (A. nodosum, A. oblongoguttatum, A. longirostre, A. pacar, A. pecarium, A. rotundatum, and A. sabanerae) and three with a Nearctic distribution (A. americanum, A. coelebs, and A. (=Robertsicus) elaphense) [61]. Eleven have presence in both biogeographic realms (A. auricularium, A. cajennense, A. calcaratum, A. dissimile, A. imitator, A. inornatum, A. maculatum, A. ovale, A. parvum, A. triste, and A. scutatum). Five species (A. humerale, A. multipuntum, A. tigrinum, A. tuberculatum, and A. varium) were not considered in this study because no reliable information was found about their specific location [25].
Amblyomma cajennense has the most extensive distribution in the country [25], and according to the literature, this species is also one of the most widely distributed in the Americas, ranging from the southern United States to northern Argentina [62]. Ticks from the genus Amblyomma are aggressive and are generalist species that participate as vectors of pathogens of medical and veterinary importance. Their presence in the northern hemisphere is associated with warm and humid localities that range in elevation up to 1000 m above sea level and with high ecosystem productivity (NDVI values < 0.56) [63,64,65]. The genus Amblyomma has been recorded in 30 of the 32 states of Mexico [25]. The states with the highest number of records are Chiapas, Tamaulipas, and Veracruz, followed by Tabasco, Yucatan, Sinaloa, Nayarit, Colima, and Oaxaca [25,62]. The ecological characteristics described above were consistent with our model, including the prediction being favorable, and the areas associated with the states bordering the Gulf of Mexico and the Pacific Ocean within the Neotropical biogeographical region were predicted as favorable (Figure 2).
The geographical distribution of Amblyomma tends to expand in response to climate change, shifts in land-use, the movement of humans and domestic animals, and the introduction of alien species [66]. Importantly, individuals that have been collected in localities that range between the 1000 and 1500 m above sea levels or in semi-dessert habitats are not documented as reproductive populations [62]. These group of species have a low resistance to desiccation [67] and poor tolerance to temperature variations [62,68,69,70].
The genus Ixodes is the most diverse of the family Ixodidae [71]. In Mexico, there are 26 species of the 243 recognized worldwide. However, most of these species are poorly represented in Mexico, represented by only one record within the country, and many of the immature stages of these species remain unknown [24]. Ixodes ricinus is recognized as a complex of several species of ticks that together have an almost cosmopolitan distribution (I. ricinus, I. scapularis, I. pacificus, I. affinis, I. pavloski, I. persulcatus, I. nipponensis, I. gibbosus, I. jellisoni, I. pararicinus, and I. nuttallianus) [72]. Subsequent phylogenetic studies determined that I. muris, I. minor, and I. granulatus should also be included in the I. ricinus complex [73,74]. The most widely distributed species from the complex in Europe are I. ricinus and I. persulcatus, while in North America, they are I. scapularis and I. pacificus. These species occupy similar niches in their respective distribution areas [75].
In North America, many of the immature stages of Ixodes species can parasitize a large number of hosts that include reptiles, birds, and mammals, whereas Ixodes adults tend to be restricted to large mammals (e.g., cervids and carnivores) [76,77,78]. Most Ixodes ticks perform vertical movements on the vegetation to reach their host. However, within the genus, it has also been observed that some species search for their hosts in relatively open environments. Regardless of their behavior, our results are consistent with the associations found previously of these ticks to ecosystems with a relatively high percentage of vegetation covering the ground since they have shown to be sensitive to extreme temperatures and dry conditions (Figure 2). For example, it was determined that unfed nymphs reached a 50% mortality at an exposure of −11.6 °C during eight consecutive hours under laboratory conditions [79]. Vegetation was found to generate microenvironmental conditions that prevent the temperature from falling below 0 °C, even during the intensely cold winter periods of the forested environments where these ticks are regularly found. Furthermore, it has been shown that Ixodes members are particularly sensitive to high temperatures (~30 °C) and water loss [80], with higher temperatures associated with increased mortality, reduced oviposition success, and host-seeking activity.
Globally, 36 species of Dermacentor are recognized [81], 11 of which are documented in Mexico. In general, Dermacentor ticks are found in all biogeographic regions [82], but they are better represented in the Holarctic region. However, D. albipictus, D. dissimilis, D. halli, D. nitens, and D. variabilis occur in both the Nearctic and Neotropical realms, and D. imitans is a Neotropical species. Native artiodactyl and perissodactyl mammals are common hosts of this species. However, D. albipictus and D. variabilis, which are the most common species in Mexico, are also associated with native carnivora and lagomorpha that inhabit many states of the country [26,83].
As with other tick genera, the distribution of Dermacentor is the result of a complex interaction between climate variables, and host and habitat availability [83,84]. Nonetheless, climate constraints seem not to be a limiting factor for Dermacentor spp. Dispersal because they have an efficient water balance, which enables them to colonize new environments [80,85]. This resistance to desiccation combined with the association of adult Dermacentor ticks with wide distributed hosts [26] could explain the board distribution in Mexico (Figure 2).
Rhipicephalus is present in all biogeographic realms. However, there are no endemic species of this genus in the Nearctic and the Neotropical regions [82], which suggests that this genus was probably introduced to the Americas with livestock or other animals from Eurasia and Africa. The genus Rhipicephalus contains what are probably the most generalist species, parasitizing amphibians, reptiles, birds, and mammals. Besides high diversity of hosts, it has been suggested that Rhipicephalus genus in the Americas also has species that have adapted well to both tropical and temperate ecosystems [86]. Within Mexico, there is a greater proportion of species adapted to tropical ecosystems, with the more suitable areas for this genus in costal ecosystems (Figure 2).
Overall, ticks have physiological requirements which require certain environmental characteristics that are optimal for their development and survival [87]. The localities where these characteristics converge can potentially be part of their geographical distribution [88], and due to their ability to maintain and transmit pathogens of veterinary and medical importance, these localities can be considered as important risk areas [89]. Although potential distribution models are extensively used, it is important to consider that different algorithm can estimate different outcomes [90]. However, the generation of an assembly of the predicted models can help identify areas of agreed suitability for ticks and for the establishment of surveillance and control programs.
Management recommendations should consider that the main limitation of ecological niche models when used to estimate the potential distribution of tick species is the taxonomic inconsistencies widely discussed by different authors [91,92,93]. The geographic distribution of a species is determined by a set of complex ecological, geological, and evolutionary processes of each taxon, so the correct geographic location of a record to a species is crucial for the generation of successful models [18,88,94,95]. We strongly believe that the estimation of genus-based ecological niche models can be useful in approaching the identification of suitable areas for these vectors of medical and veterinary importance.
The knowledge about tick diversity and the ecological features that determine their distribution in Mexico is far from completion. We are currently at a decisive moment for humanity, in which we must not only be able to deal the consequences of the changes we have exerted directly on biodiversity, land use, and climate but also must understand how these modifications influence complex processes such as vector dynamics and the diseases they transmit. The technical improvements to ecological niche modeling and the generation of quality databases may assist the decision-making processes during this period of uncertainty to prevent the emergence of infectious outbreaks that threaten human wellbeing and animal health.

Supplementary Materials

The following are available online at https://www.mdpi.com/2076-2615/10/4/649/s1, Table S1: Records used for the generation of ecological niche models for the four genera of ticks. Table S2: Climatic and environmental variables analyzed to make ecological niche models for the four tick genera.

Author Contributions

Conceptualization, C.A.L.-G. and E.C.-C.; methodology, C.N.M.-A. and E.C.-C.; software, C.M. and E.C.-C.; validation, E.C.-C. and C.N.M.-A.; formal analysis, E.C.-C., C.N.M.-A., and C.A.L.-G.; investigation, E.C.-C.; resources, E.C.-C.; data curation, E.C.-C. and C.N.M.-A.; writing—original draft preparation, E.C.-C.; writing—review and editing, E.C.-C., C.N.M.-A., and C.A.L.-G.; visualization, E.C.-C., C.N.M.-A., and C.A.L.-G.; supervision, C.A.L.-G.; project administration, C.A.L.-G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Programa para el Desarrollo Profesional Docente, para el Tipo Superior in the project 511-6/17-11904 of the Cuerpo Académico Ecología y Diversidad Faunisitica of the Universiad Autónoma de Querétaro.

Acknowledgments

We would like the social service students of the Zoology Laboratory of the Autonomous University of Queretaro for their assistance in data collection. Many thanks to Robert W. Jones for reviewing the manuscript format in English and to the anonymous reviewers of this research.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Tabachnick, W.J. Challenges in predicting climate and environmental effects on vector-borne disease episystems in a changing world. J. Exp. Biol. 2010, 213, 946–954. [Google Scholar] [CrossRef] [Green Version]
  2. Ramos, N.; Rodrigues, S.; Piovezan, U.; Pablo, M.; Szabó, J. Microhabitat determines uneven distribution of Amblyomma parvum but not of Amblyomma sculptum ticks within forest patches in the Brazilian Pantanal. Exp. Appl. Acarol. 2019, 79, 405–410. [Google Scholar] [CrossRef]
  3. Peterson, A.T. Ecologic niche modeling and spatial patterns of disease transmission. Emerg. Infect. Dis. 2006, 12, 1822–1826. [Google Scholar] [CrossRef]
  4. Soberón, J.; Peterson, A.T. Interpretation of models of fundamental ecological niches and species distributional areas. Biodivers. Inf. 2005, 2, 1–10. [Google Scholar] [CrossRef] [Green Version]
  5. Soberón, J.M. Niche and area of distribution modeling: A population ecology perspective. Ecogr. Cop. 2010, 33, 159–167. [Google Scholar] [CrossRef]
  6. Rangel, T.F.; Loyola, R.D. Labeling Ecological Niche Models. Nat. Conserv. 2012, 10, 119–126. [Google Scholar] [CrossRef]
  7. Olden, J.D.; Lawler, J.J.; Poff, N.L. Machine Learning Methods without Tears: A primer for Ecologists. Q. Rev. Biol. 2008, 83, 171–193. [Google Scholar] [CrossRef] [Green Version]
  8. Parola, P.; Raoult, D. Ticks and tickborne bacterial diseases in humans: An emerging infectious threat. Clin. Infect. Dis. 2001, 32, 897–928. [Google Scholar] [CrossRef] [PubMed]
  9. Smith, K.F.; Guégan, J.-F. Changing Geographic Distributions of Human Pathogens. Annu. Rev. Ecol. Evol. Syst. 2010, 41, 231–250. [Google Scholar] [CrossRef] [Green Version]
  10. Jongejan, F.; Uilenberg, G. The global importance of ticks. Parasitology 2004, 129, S3–S14. [Google Scholar] [CrossRef] [PubMed]
  11. Amsden, J.R.; Warmack, S.; Gubbins, P.O. Tick-borne bacterial, rickettsial, spirochetal, and protozoal infectiouss diseases in the United States: A comprehenssive review. Pharmacotherapy 2005, 25, 191–210. [Google Scholar] [CrossRef]
  12. Walker, A.R. Ticks and associated diseases: A retrospective review. Med. Vet. Entomol. 2014, 28, 1–5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sprong, H.; Tijsse-Klasen, E.; Langelaar, M.; De Bruin, A.; Fonville, M.; Gassner, F.; Takken, W.; Van Wieren, S.; Nijhof, A.; Jongejan, F.; et al. Prevalence of Coxiella Burnetii in ticks after a large outbreak of Q fever. Zoonoses Public Health 2011, 59, 69–75. [Google Scholar] [CrossRef]
  14. Chomel, B. Tick-borne infections in dogs-an emerging infectious threat. Vet. Parasitol. 2011, 179, 294–301. [Google Scholar] [CrossRef] [PubMed]
  15. Dantas-Torres, F. Climate change, biodiversity, ticks and tick-borne diseases: The butterfly effect. Int. J. Parasitol. 2015, 4, 452–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Esteve-Gassent, M.D.; Castro-Arellano, I.; Feria-Arroyo, T.P.; Patino, R.; Li, A.Y.; Medina, R.F.; de León, A.A.P.; Rodríguez-Vivas, R.I. Translating ecology, physiology, biochemistry, and population genetics research to meet the challenge of tick and tick-borne diseases in North America. Arch. Insect Biochem. Physiol. 2016, 92, 38–64. [Google Scholar] [CrossRef] [Green Version]
  17. Illoldi-rangel, P.; Rivaldi, C.; Sissel, B.; Fryxell, R.T.; 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, 2012, 959101. [Google Scholar] [CrossRef] [Green Version]
  18. Feria-Arroyo, T.P.; Castro-Arellano, I.; Gordillo-Perez, G.; Cavazos, A.L.; Vargas-Sandoval, M.; Grover, A.; Torres, J.; Medina, R.F.; de León, A.A.P.; Esteve-Gassent, M.D. Implications of climate change on the distribution of the tick vector Ixodes scapularis and risk for Lyme disease in the Texas-Mexico transboundary region. Parasit. Vectors 2014, 7, 199. [Google Scholar] [CrossRef] [Green Version]
  19. González, C.; Wang, O.; Strutz, S.E.; González-Salzar, 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. [Google Scholar] [CrossRef] [Green Version]
  20. Gurgel-Gonçalves, R.; Galvão, C.; Costa, J.; Peterson, A.T. Geographic distribution of chagas disease vectors in Brazil based on ecological niche modeling. J. Trop. Med. 2012, 2012, 705326. [Google Scholar] [CrossRef] [Green Version]
  21. Gorla, D.; Noireau, F. Geographic distribution of Triatominae vectors in America. In American Trypanaosomiasis Chagas Disease: One Hundred Years of Research; Telleria, J., Tibayrenc, M., Eds.; Elsevier: Amsterdam, The Netherlands, 2010; pp. 209–231. [Google Scholar]
  22. Moffett, A.; Strutz, S.; Guda, N.; Gonález, C.; Ferro, M.C.; Sánchez-Cordero, V.; Sarkar, S. A global public database of disease vector and reservoir distributions. PLoS Negl. Trop. Dis. 2009, 3, 1–3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Stevens, K.B.; Gilbert, M.; Pfeiffer, D.U. Modeling habitat suitability for occurrence of highly pathogenic avian influenza virus H5N1 in domestic poultry in Asia: A spatial multicriteria decision analysis approach. Spat. Spatio-Tempor. Epidemiol. 2013, 4, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Guzmán-Cornejo, C.; Robbins, R.G. The genus Ixodes (Acari: Ixodidae) in Mexico: Adult identification keys, diagnoses, hosts, and distribution. Rev. Mex. Biodivers. 2010, 81, 289–298. [Google Scholar]
  25. Guzmán-Cornejo, C.; Robbins, R.G.; Guglielmone, A.A.; Montiel, G.; Pérez, T.M. The Amblyomma (Acari: Ixodida: Ixodidae) of Mexico: Identification keys, distribution and hosts. Zootaxa 2011, 2998, 16–38. [Google Scholar] [CrossRef]
  26. Guzmán-Cornejo, C.; Robbins, R.G.; Guglielmone, A.A.; Montiel-Parra, G.; Rivas, G.; Pérez, T.M. The Dermacentor (Acari, Ixodida, Ixodidae) of Mexico: Hosts, geographical distribution and new records. Zookeys 2016, 22, 1–22. [Google Scholar] [CrossRef]
  27. Huberty, C.J. Applied Discriminant Analysis; Wiley-Interscience: New York, NY, USA; University of Michigan: Ann Arbor, MI, USA, 1994; 496p. [Google Scholar]
  28. Barbet-Massin, M.; Jiguet, F.; Albert, C.H.; Thuiller, W. Selecting pseudo-absences for species distribution models: How, where and how many? Methods Ecol. Evol. 2012, 3, 327–338. [Google Scholar] [CrossRef]
  29. Di Cola, V.; Broennimann, O.; Petitpierre, B.; Breiner, F.T.; Amen, M.D.; Randin, C.; Engler, R.; Pottier, J.; Pio, D.; Dubuis, A.; et al. ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecogr. Cop. 2016, 40, 774–787. [Google Scholar] [CrossRef]
  30. Jiménez-Valverde, A.; Lobo, J.M.; Hortal, J. The effect of prevalence and its interaction with sample size on the reliability of species distribution models. Community Ecol. 2008, 10, 196–205. [Google Scholar] [CrossRef]
  31. Stokland, J.N.; Halvorsen, R.; Støa, B. Species distribution modelling—Effect of design and sample size of pseudo-absence observations. Ecol. Model. 2011, 222, 1800–1809. [Google Scholar] [CrossRef]
  32. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Clim. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  33. Randolph, S.E. Ecology of non-nidicolous ticks. In Biology of Ticks, 2nd ed.; Sonenshine, D.E., Roe, R.M., Eds.; Oxford University Press: Oxford, UK, 2014; pp. 3–37. [Google Scholar]
  34. Warren, D.L.; Seifert, S.N. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model. Ecol. Appl. 2011, 21, 335–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Barve, N.; Barve, V.; Jiménez-Valverde, A.; Lira-Noriega, A.; Maher, S.P.; Peterson, A.T.; Soberón, J.; Villalobos, F. The crucial role of the accessible area in ecological niche modeling and species distribution modeling. Ecol. Model. 2011, 222, 1810–1819. [Google Scholar] [CrossRef]
  36. Naimi, B.; Araújo, M.B. Sdm: A reproducible and extensible R platform for species distribution modelling. Ecogr. Cop. 2016, 39, 368–375. [Google Scholar] [CrossRef] [Green Version]
  37. Muscarella, R.; Galante, P.J.; Soley-guardia, M.; Boria, R.A.; Kass, J.M.; Anderson, R.P. ENMeval: An R package for conducting spatially independent evaluations and estimating optimal model complexity for MAXENT ecological niche models. Methods Ecol. Evol. 2014, 5, 1198–1205. [Google Scholar] [CrossRef]
  38. Hastie, T. Generalized Additive Models [Internet]. R Package. 2019. Available online: http://cran.r-project.org/web/packages/gam/gam.pdf (accessed on 8 January 2020).
  39. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  40. Phillips, S.J.; Ddudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Echography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  41. McCullagh, P. Generalized linear models. Eur. J. Oper. Res. 1984, 16, 285–292. [Google Scholar] [CrossRef]
  42. McCullagh, P.; Nelder, J.A. Generalized Linear Models; Chapman and Hall: Boca Raton, FL, USA, 1989; 532p. [Google Scholar]
  43. Friedman, H.H. Multivariate adaptive regression splines. Ann. Stat. 1991, 19, 1–141. [Google Scholar] [CrossRef]
  44. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Classification and Regression Trees; Chapman and Hall: Boca Raton, FL, USA, 1984; 386p. [Google Scholar]
  45. Hastie, T.; Tibshirani, R.; Buja, A. Flexible discriminant analysis by optimal scoring. J. Am. Stat. Assoc. 1994, 89, 1255–1270. [Google Scholar] [CrossRef]
  46. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  47. Friedman, J.H. Greedy function approximation: A gradient boosting machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  48. Elith, J.; Leathwick, J.R.; Hastie, T. A working guide to boosted regression trees. J. Anim. Ecol. 2008, 77, 802–813. [Google Scholar] [CrossRef] [PubMed]
  49. Nix, H.A. A biogeographic analysis of Australian Elapid snakes. In Australian Flora and Fauna Series Number 7: Atlas of Elapid Snakes of Australia; Lonngmore, R., Ed.; Australian Goverment Publishing Service: Canberra, Australia, 1986; pp. 4–5. [Google Scholar]
  50. Booth, T.H.; Nix, H.A.; Busby, J.R.; Hutchinson, M.F. BIOCLIM: The firs species distribution modelling packages, its early applications and relevance to most current MAXENT studies. Divers. Distrib. 2014, 20, 1–9. [Google Scholar] [CrossRef]
  51. Elith, J.; Graham, C.H.; Anderson, R.P.; Duudík, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.; Leathwick, J.R.; Lehmann, A.; et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecogr. Cop. 2006, 2, 129–151. [Google Scholar] [CrossRef] [Green Version]
  52. Guisan, A.; Edwards, T.C.; Hastie, T. Generalized linear and generalized additive models in studies of species distributions: Setting the scene. Ecol. Model. 2002, 157, 89–100. [Google Scholar] [CrossRef] [Green Version]
  53. Leathwick, J.R.; Elith, J.; Hastie, T. Comparative performance of generalized additive models and multivariate adaptive regression splines for statistical modelling of species distributions. Ecol. Model. 2006, 199, 188–196. [Google Scholar] [CrossRef]
  54. Hijmans, R.J.; Phillips, S.; Leathwick, J.R.; Elith, J. Package ‘Dismo‘. Available online: https://cran.r-project.org/web/packages/dismo/index.html (accessed on 8 January 2020).
  55. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; Available online: http://www.r-project.org/ (accessed on 8 January 2020).
  56. Swets, J.A. Measuring the accuracy of diagnostic systems. Science 1988, 240, 1285–1293. [Google Scholar] [CrossRef] [Green Version]
  57. Fiellding, A.H.; Bell, J.F. A review of methods for the assessment of prediction errors in conservation presence/absence models. Environ. Conserv. 1997, 24, 38–49. [Google Scholar] [CrossRef]
  58. Jiménez-Valverde, A. Insights into the area under the receiver operating characteristic curve (AUC) as a discrimination measure in species distribution modelling. Glob. Biogeogr. 2012, 21, 498–507. [Google Scholar] [CrossRef]
  59. Guisan, A.; Graham, C.H.; Elith, J.; Huettmann, F. NCEAS Species Distribution Modelling Group. Sensitivity of predictive species distribution models to change in grain size. Divers. Distrib. 2007, 13, 332–340. [Google Scholar] [CrossRef]
  60. Radosavljevic, A.; Anderson, R.P. Making better Maxent models of species distributions: Complexity, overfitting and evaluation. J. Biogeogr. 2014, 41, 629–643. [Google Scholar] [CrossRef]
  61. Barker, S.C.; Burger, T.D. Two new genera of hard ticks, Robertsicus n. gen. and Archaeocroton n. gen.; and the solution to the mystery of Hoogstraal’s and Kaufman’s “primitive” tick from the Carpathian Mountains. Zootaxa 2018, 4500, 543–552. [Google Scholar] [CrossRef] [PubMed]
  62. Estrada-Peña, A.; Guglielmone, A.A.; Mangold, A.J. The distribution and ecological “preferences” of the tick Amblyomma cajennense (Acari: Ixodidae), an ectoparasite of humans and other mammals in the Americas. Ann. Trop. Med. Parasitol. 2004, 98, 283–292. [Google Scholar] [CrossRef]
  63. Bishopp, F.C.; Hixson, H. Biology and economic importance of the Gulf Coast tick. J. Econ. Entomol. 1936, 29, 1068–1076. [Google Scholar] [CrossRef]
  64. Goddard, J.; Norment, B.R. Notes on the geographical distribution of the Gulf coast tick, Amblyomma maculatum (Koch) [Acari, Ixodidae]. Entomol. News 1983, 94, 103–104. [Google Scholar]
  65. Estrada-peña, A.; Tarragona, E.L.; Vesco, U.; Meneghi, D.; De Mastropaolo, M.; Mangold, A.J.; Guglielmone, A.A.; Nava, S. Divergent environmental preferences and areas of sympatry of tick species in the Amblyomma cajennense complex (Ixodidae). Int. J. Parasitol. 2014, 44, 1081–1089. [Google Scholar]
  66. Pascoe, E.L.; Marcantonio, M.; Caminade, C.; Foley, J.E. Modeling Potential Habitat for Amblyomma Tick Species in California. Insects 2019, 10, 201. [Google Scholar] [CrossRef] [Green Version]
  67. Strey, O.F.; Teel, P.D.; Longnecker, M.T.; Needham, G.R. Survival and Water-Balance Characteristics of Unfed Adult Amblyomma cajennense (Acari: Ixodidae). J. Med. Entomol. 1996, 33, 63–73. [Google Scholar] [CrossRef]
  68. Needham, G.R.; Teel, P.D. Off-host physiological ecology of ixodid ticks. Annu. Rev. Entomol. 1991, 36, 659–681. [Google Scholar] [CrossRef]
  69. Ketchum, H.R.; Teel, P.D.; Coates, C.J.; Strey, O.F.; Longnecker, M.T. Genetic variation in 12S and 16 S mitochondrial r DNA genes of four geographically isolated populations of Gulf Coast ticks (Acari: Ixodidae). J. Med. Entomol. 2009, 46, 482–489. [Google Scholar] [CrossRef]
  70. Sanders, D.M.; Schuster, A.L.; McCardle, P.W.; Strey, O.F.; Blankenship, T.L.; Teel, P.D. Ixodid ticks associated with feral swine in Texas. J. Vector Ecol. 2013, 38, 361–373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Guglielmone, A.A.; Robbins, R.G.; Apanaskevich, D.A.; Petney, T.N.; Estrada-Peña, A.; Hoorak, I.G.; Shao, R.; Barker, S.C. The Argasidae, Ixodidae and Nuttalliellidae (Acari: Ixodida) of the world: A list of valid species names. Zootaxa 2010, 2528, 1–28. [Google Scholar] [CrossRef] [Green Version]
  72. Keirans, J.E.; Needham, G.R.; Oliver, J.H. The Ixodes (Ixodes) ricinus complex worldwide: Diagnosis of the species in the complex, hosts and distribution. In Acarology IX Symposia; University of Wisconsin-Madison: Madison, WI, USA, 1999; pp. 341–347. [Google Scholar]
  73. Xu, G.; Fang, Q.Q.; Keirans, J.E.; Durden, L.A. Molecular phylogenetic analyses indicate that the Ixodes ricinus complex is a paraphyletic group. J. Parasitol. 2003, 89, 452–457. [Google Scholar] [CrossRef]
  74. Beati, L.; Klompen, H. Phylogeography of ticks (Acari: Ixodida). Annu. Rev. Entomol. 2019, 64, 379–397. [Google Scholar] [CrossRef] [PubMed]
  75. Gray, J.S.; Kahl, O.; Lane, R.S.; Levin, M.L.; Tsao, J.I. Diapause in ticks on the medically important Ixodes ricinus species complex. Ticks Tick-Borne Dis. 2016, 7, 992–1003. [Google Scholar] [CrossRef] [Green Version]
  76. LoGiudice, K.; Ostfeld, R.S.; Schmidt, K.A.; Keesing, F. The ecology of infectious disease: Effects of host diversity and community composition on Lyme disease risk. Proc. Natl. Acad. Sci. USA 2003, 100, 567–571. [Google Scholar] [CrossRef] [Green Version]
  77. Eisen, R.J.; Eisen, L.; Lane, R.-S. Habitat related variation in infestation of lizards and rodents with Ixodes ticks in dense woodlands in Mendocino County, California. Exp. Appl. Acarol. 2004, 33, 215–233. [Google Scholar] [CrossRef]
  78. Salkeld, D.J.; Lane, R.S. Community ecology and disease risk: Lizards, squirrels, and the Lyme disease spirochete in California, USA. Ecology 2010, 91, 293–298. [Google Scholar] [CrossRef]
  79. Vandyk, J.K.; Bartholomew, D.M.; Rowley, W.A.; Platt, K.B. Survival of Ixodes scapularis (Acari: Ixodidae) exposed to cold. J. Med. Entomol. 1996, 33, 6–10. [Google Scholar] [CrossRef]
  80. Yoder, J.A.; Chris, B.S.; Croxall, T.J.; Laura, K.S.; Tank, J.L. Moisture requirements for activity/survival of the gulf coast tick Amblyomma maculatum Koch (Acari: Ixodidae), based on a water balance study of all life cycle stages. Int. J. Acarol. 2009, 34, 285–292. [Google Scholar] [CrossRef]
  81. Wang, F.; Wang, D.; Guo, G.; Hu, Y.; Wei, J.; Liu, J. Species delimitation of the Dermacentor ticks based on phylogenetic clustering and niche modeling. Peer J. 2019, 7, e6911. [Google Scholar] [CrossRef] [Green Version]
  82. Guglielmone, A.A.; Robbins, R.G.; Apanaskevich, D.A.; Petney, T.N.; Estrada-Peña, A.; Horak, I.G. The Hard Ticks of the World; Springer: New York, NY, USA, 2013; 738p. [Google Scholar]
  83. Minigan, J.N.; Hager, H.A.; Peregrine, A.S.; Newman, J.A. Current and potential future distribution of the American dog tick (Dermacentor variabilis, Say) in North America. Tick Tick-Borne Dis. 2018, 9, 354–362. [Google Scholar] [CrossRef] [PubMed]
  84. Süss, J.; Klaus, C.; Gerstengarbe, F.; Werner, P.C. What makes ticks tick? Climate change, ticks, and tick-borne diseases. J. Travel Med. 2008, 15, 39–45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Holmes, C.J.; Dobrotka, C.J.; Farrow, D.W.; Rosendale, A.J.; Benoit, J.B.; Pekins, P.J.; Yoder, J.A. Low and high thermal tolerance characteristics for unfed larvae of the winter tick Dermacentor albipictus (Acari: Ixodidae) with special reference to moose. Ticks Tick-Borne Dis. 2018, 9, 25–30. [Google Scholar] [CrossRef] [PubMed]
  86. Laburna, M.B.; Gerardi, M.; Krawczak, F.S.; Moraes-Filho, J. Comparative biology of the tropical and temperate species of Rhipicephalus sanguineus sensu lato (Acari: Ixodidae) under different laboratory conditions. Ticks Tick-Borne Dis. 2017, 8, 146–156. [Google Scholar] [CrossRef] [PubMed]
  87. Kessler, W.H.; Ganser, C.; Glass, G.E. Modelling the distribution of medically important tick species in Florida. Insects 2019, 10, 190. [Google Scholar] [CrossRef] [Green Version]
  88. Guisan, A.; Zimmermann, N. Predictive habitat distribution models in ecology. Ecol. Model. 2000, 135, 147–186. [Google Scholar] [CrossRef]
  89. Bouchard, C.; Dibernardo, A.; Koffi, J.; Wood, H.; Leighton, P.A.; Lindssay, L.R. Increased risk of tick-borne diseases with climate and environmental changes. Can. Commun. Dis. Rep. 2019, 45, 81–89. [Google Scholar] [CrossRef]
  90. Marmion, M.; Parviainen, M.; Luoto, M.; Heikkinen, R.K.; Thuiller, W. Evaluation of consensus methods in predictive species distribution modelling. Divers. Distrib. 2009, 15, 59–69. [Google Scholar] [CrossRef]
  91. Camicas, J.-L.; Hervy, J.-P.; Adam, F.; Morel, P.-C. Les tiques du monde (Acarida, Ixodida). Nomenclature, stades décrits, hotes, répartion; Orstom: Paris, France, 1998; 233p. (In French) [Google Scholar]
  92. Horak, I.G.; Horak, J.L.; Keiranss, J.E. The Argasidae, Ixodidae, and Nuttalliellidae (Acari: Ixodida): A world list of valid tick names. Exp. Appl. Acarol. 2002, 28, 27–54. [Google Scholar] [CrossRef]
  93. Barker, S.C.; Murrell, A. Systematics and Evolution of ticks a list of valid genus and species names. In Ticks: Biology Disease and Control; Bowman, A.S., Nuttall, P., Eds.; Cambridge University Press: Cambridge, UK, 2008; pp. 1–39. [Google Scholar]
  94. Austin, M. Species distribution models and ecological theory: A critical assessment and some possible new approaches. Ecol. Model. 2007, 200, 1–19. [Google Scholar] [CrossRef]
  95. Mori, E.; Menchetti, M.; Zozzoli, R.; Milaanesi, P. The importance of taxonomy in species distribution models at a global scale: The case of an overlooked alien squirrel facing taxonomic revision. J. Zool. 2018, 307, 43–52. [Google Scholar] [CrossRef]
Animals 10 00649 g001 550
Figure 1. Records used to estimate the ecological niche models for the four tick genera: Black circles indicate the presence of Rhipicephalus spp., red triangles indicate the presence of Ixodes spp., and green and blue boxes show the presence of Dermacentor spp. and Amblyomma spp., respectively.
Figure 1. Records used to estimate the ecological niche models for the four tick genera: Black circles indicate the presence of Rhipicephalus spp., red triangles indicate the presence of Ixodes spp., and green and blue boxes show the presence of Dermacentor spp. and Amblyomma spp., respectively.
Animals 10 00649 g001
Animals 10 00649 g002 550
Figure 2. Assembly that brings together the estimated potential distribution of the four tick genera considering the nine algorithms: The darkest areas indicate agreement of the prediction of habitat suitability for each tick genus. The values indicate the number of methods that coincided in determining a certain site as suitable for each genus of tick.
Figure 2. Assembly that brings together the estimated potential distribution of the four tick genera considering the nine algorithms: The darkest areas indicate agreement of the prediction of habitat suitability for each tick genus. The values indicate the number of methods that coincided in determining a certain site as suitable for each genus of tick.
Animals 10 00649 g002
Table
Table 1. Climatic and environmental variables used for the generation of ecological niche models of the four tick genera.
Table 1. Climatic and environmental variables used for the generation of ecological niche models of the four tick genera.
GeneraAmblyommaDermacentor, Ixodes and Rhipicephalus
Annual mean temperature •
Mean diurnal range ••
Temperature Seasonality •
Max temperature of warmest month •
Min temperature of coldest month•
Mean temperature of wettest quarter •
Annual precipitation••
Precipitation of driest month••
Precipitation seasonality ••
Precipitation of wettest quarter •
Precipitation of the warmest quarter•
Type of soil••
Type of land use and vegetation••
Table
Table 2. Area under the curve index (AUC) values for each of the models generated by tick genus.
Table 2. Area under the curve index (AUC) values for each of the models generated by tick genus.
AlgorithmAmblyomma spp.Dermacentor spp.Ixodes spp.Rhipicephalus spp.
BIOCLIM0.7060.6640.6690.772
BRT0.9050.8720.8920.921
CART0.9130.8560.8830.965
MDA0.8830.7970.7890.869
GAM0.930.8710.9620.98
GLM0.8880.8040.7920.878
MARS0.920.9410.9470.962
MAXENT0.9010.8400.9180.931
RF0.9990.9960.9940.999
BRT: Boosted Regression Trees, CART: Classification And Regression Tree, MDA: Mixture Discriminant Analysis, GAM: Generalized Additive Models, GLM: Generalized Linear Models; MARS: Multivariate adaptive regression spline, Maxent: Maximum entropy, RF: Random Forest

Share and Cite

MDPI and ACS Style

Clarke-Crespo, E.; Moreno-Arzate, C.N.; López-González, C.A. Ecological Niche Models of Four Hard Tick Genera (Ixodidae) in Mexico. Animals 2020, 10, 649. https://doi.org/10.3390/ani10040649

AMA Style

Clarke-Crespo E, Moreno-Arzate CN, López-González CA. Ecological Niche Models of Four Hard Tick Genera (Ixodidae) in Mexico. Animals. 2020; 10(4):649. https://doi.org/10.3390/ani10040649

Chicago/Turabian Style

Clarke-Crespo, Emilio, Claudia N. Moreno-Arzate, and Carlos A. López-González. 2020. "Ecological Niche Models of Four Hard Tick Genera (Ixodidae) in Mexico" Animals 10, no. 4: 649. https://doi.org/10.3390/ani10040649

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