Shifts in Climatic Niche Occupation in Astrophytum Coahuilense ( H . Möller ) Kayser and Its Potential Distribution in Mexico

It has frequently been reported that species with strong niche conservatism will not be able to adapt to new climatic conditions, so they must migrate or go extinct. We have evaluated the shifts in climatic niche occupation of the species Astrophytum coahuilense and its potential distribution in Mexico. We understand niche occupation as the geographic zones with available habitats and with the presence of the species. To assess shifts in climatic niche occupation, we used niche overlap analysis, while potential distribution modeling was performed based on the principle of maximum entropy. The results indicate that this species presents a limited amplitude in its climate niche. This restriction of the climatic niche of A. coahuilense limits its ability to colonize new geographical areas with different climatic environments. On the other hand, the potential distribution models obtained from the present study allow us to identify potential zones based on the climatic requirements of the species. This information is important to identify high priority areas for the conservation of A. coahuilense.


Introduction
Mexico is the most important center of diversity for genus and species of cactus.Most of the species occur in arid or semiarid regions of the country, particularly in the southeastern part of the Chihuahuan Desert [1].Cactaceae is one of the most threatened groups of plants, and is faced with conservation problems due to overexploitation and illegal trade [2].Astrophytum coahuilense is an endemic cactus from Mexico that occurs in southwestern Coahuila and adjacent areas of Durango, México [3].This species grows predominately on rock-strewn, calcareous terrain, together with conspicuous xerophytes, which include Larrea tridentata, Prosopis juliflora, Fouquiera splendens, Agave lechuguilla, Euphorbia antisyphilitica, and Hechtia sp.[4].Moreover, A. coahuilense has a relatively small geographic range and faces several threats to its survival [4].Its principal method of dispersion is hydrochloric, and its current range of distribution was established through lakes and rivers of the Pleistocene.Therefore, its current dispersal capacity is limited [5].It has been mentioned that the main factors responsible for the decline of its populations are the fragmentation and reduction of habitat quality [6].
It is generally expected that endemic species with restricted distribution ranges would be particularly threatened by alterations like climate [7,8].It has been suggested that climatic variations would severely affect all species of Astrophytum, as most of their populations may experience extreme climatic conditions in the next decades [7,9].Nevertheless, in response to these threats and in seeking to evade extinction, species may track suitable environmental conditions in space or time, or adapt to the modified environment through acclimation [10].
Assuming that species need to migrate or adapt to changing environmental conditions, establishment success, and subsequent distribution within a given geographic range, is determined by the species' association to abiotic variables such as climate.These associations can be interpreted through the niche concept [11].The ecological niche of a species can be defined as the range of environmental and biotic conditions within which its populations can persist [12].Two distinct formulations of the niche concept have been considered: 1) fundamental niche, which is genetically and physiologically determined, and 2) realized niche, which includes, additionally, constraints arising from interspecific competition [13].
Species distribution has been studied from a niche conservatism perspective (the tendency of species to retain ancestral ecological characteristics [14]) and niche shift perspective (any change in the position of either the fundamental or realized niche of a species [15]), or both [13].In this way, niche conservatism and niche shifts can have important implications for understanding the effects of climate change on species distribution [14,15].
Niche conservatism has been evaluated in terms of statistical difference, i.e., "niche similarity" (the geographical areas occupied by a species are more similar in climate than expected by chance) and "niche equivalence" (the geographic areas occupied by a species are climatically indistinguishable)" [16].Nevertheless, the similarity and equivalence tests only reveal statistically significant differences in species climate niches, but not the causes underlying such niche shifts [17,18].Hence, absolute values of "equivalence" and "similarity" have little biological significance unless niche shifts are disentangled from their main drivers of change, namely "niche unfilling" and "niche expansion" [19,20].
The threat of environmental disruption for biodiversity can be analyzed from the perspective of niche conservatism [21].Species with this characteristic will not have sufficiently broad climatic tolerance to acclimatize to the new environmental conditions of the habitat, and consequently they will have to emigrate or go extinct [22].Those species that do not exhibit niche conservatism could increase their range beyond their native climatic environment, or adapt to new environmental conditions [13].Considering the above, the question arises of whether A. coahuilense will have the capacity to occupy niches with different climatic conditions, and adapt to the new niches generated by global climate change and local climatic variations resulting from degradation of the habitat and human activities.Therefore, the aims of the present study are: i) to evaluate the shifts in climatic niche occupation in Astrophytum coahuilense, and ii) to model its potential distribution in Mexico.This allows us to infer the vulnerability of this species to climatic variations, and to identify potential sites for the conservation of this species.

Geographical Records
During 10 years of intermittent sampling, 32 records of A. coahuilense were identified inside and outside the distribution area proposed by Fitz-Maurice et al. [3,6].Each record corresponds to a population of this species (Figure 1).Within this area, three zones with the presence of this cactus were classified (considering the political division and geographical distances): zone A (15 records) included Sierra de Jimulco and Sierra de las Noas in Coahuila, zone B (9 records) included Sierra el Sarnoso and Sierra de Lerdo in Durango, and zone C (4 records) included Sierra de Parras and Sierra de Mayrán in Coahuila.Four records were located outside the distribution area proposed by Fitz-Maurice et al. [3,6], three to the southwest in the state of Nuevo León (Sierra de las Mitras) and one to the southeast in the state of Coahuila (Sierra de la Paila).In the zones identified, the minimum distance between the records of presence of A. coahuilense was at least one kilometer.All geographic records were visualized in the Google Earth program in order to verify their accuracy (Figure 1).allows us to infer the vulnerability of this species to climatic variations, and to identify potential sites for the conservation of this species.

Geographical Records
During 10 years of intermittent sampling, 32 records of A. coahuilense were identified inside and outside the distribution area proposed by Fitz-Maurice et al. [3,6].Each record corresponds to a population of this species (Figure 1).Within this area, three zones with the presence of this cactus were classified (considering the political division and geographical distances): zone A (15 records) included Sierra de Jimulco and Sierra de las Noas in Coahuila, zone B (9 records) included Sierra el Sarnoso and Sierra de Lerdo in Durango, and zone C (4 records) included Sierra de Parras and Sierra de Mayrán in Coahuila.Four records were located outside the distribution area proposed by Fitz-Maurice et al. [3,6], three to the southwest in the state of Nuevo León (Sierra de las Mitras) and one to the southeast in the state of Coahuila (Sierra de la Paila).In the zones identified, the minimum distance between the records of presence of A. coahuilense was at least one kilometer.All geographic records were visualized in the Google Earth program in order to verify their accuracy (Figure 1).In the study areas, the black triangles correspond to the records of zone A, the black circles refer to the records of zone B, the black boxes indicate the records of zone C, and the black stars refer to the points identified outside the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6].The Black dotted line shows the distribution proposed by Fitz-Maurice et al. [3,6], while solid black lines refer to the state limits.In the study areas, the black triangles correspond to the records of zone A, the black circles refer to the records of zone B, the black boxes indicate the records of zone C, and the black stars refer to the points identified outside the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6].The Black dotted line shows the distribution proposed by Fitz-Maurice et al. [3,6], while solid black lines refer to the state limits.

Climatic Variables
Climatic information was obtained from 19 current climatic layers and elevation data, available in WorldClim database version 1.4 [23].These layers contain climatic averages of weather conditions recorded from 1950-2000 with a spatial resolution of 30 arc-seconds (~1 km).For the selection of environmental variables, 10,000 background points were added to the polygon distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] (Figure 1).The information of the 19 environmental variables from the current climate was added to these points.With the generated information, a bivariate correlation analysis was conducted in order to reduce the multicollinearity between the input variables [24,25]; predictor variables that were highly correlated (|r| ≥ 0.7) were excluded.The addition of 10,000 random points was based on the criteria of not discriminating (non-repetitive) relevant information, but segregated geographical areas within the range of the A. coahuilense climatic information [26].With the aforementioned variables, the climatic profile of A. coahuilense was generated.Using a principal component analysis (PCA), the distribution of the species in the climatic environment was identified using the software R, (R version 3.1.3,library ecospat; [27]).

Niche Overlap
From geographical records of A. coahuilense distributed in zone A and zone B, and from the six climatic variables retained in the correlation analysis, the climatic niche differentiation was evaluated.Zone C and records outside of the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] were not considered for this analysis because they did not reach the minimum number of geographic registers recommended by Broennimann et al. [27].The principal component analysis (PCA) approach proposed by Broennimann et al. [19] was used to measure equivalence and similarity between the realized niche of zone A and the realized niche of zone B (with the variables retained in bivariate correlation analysis).This method compares the environmental conditions available for a species within a defined study extent (background) with its observed occurrences, and it calculates the available environmental space defined by the first two axes from the PCA.This method corrects for sampling bias using a smooth kernel density function [19].
The niche overlap between zone A and B was calculated using Schoener's D-metric [28], which varies from 0 (no overlap between niches) to 1 (complete overlap).The niche equivalency and similarity tests were built from the methodology described in Broennimann et al. [19].The niche equivalency test determines whether niches of two entities in two geographical ranges are equivalent (i.e., whether the niche overlap is constant when randomly reallocating the occurrences of both entities among the two ranges).All occurrences are pooled and randomly split into two datasets, maintaining the number of occurrences as in the original dataset.This process is repeated 100 times to ensure that the null hypothesis can be rejected with high confidence.
If the observed value of D falls within the density of 95% of the simulated values, the null hypothesis of niche equivalency cannot be rejected.Therefore, the niche similarity test differs from the equivalency test because the former examines whether the overlap between observed niches in two ranges is different from the overlap between the observed niche in one range and niches selected at random from the other range.The test of niche similarity is also based on 100 repetitions.If the observed overlap is greater than 95% of the simulated values, the entity occupies environments in both of its ranges that are more similar to each other than expected.The proportion of the niche in zone B that was also measured was stable (i.e., the niche B overlapping with the niche zone A), unfilled (i.e., a partial filling of the zone A niche in the zone B range), and expanding (i.e., the invaded niche non-overlapping with the native niche) compared to its native niche [17,20].The niche overlap analysis works with a minimum of five geographic records for each evaluated species.All analyses were computed with the ecospat [27] package in R 3.1.3[29].

Niche Modeling
A maximum entropy model (MaxEnt, version 3.3.3k)[30,31] was used to represent the potential distribution of A. coahuilense.The MaxEnt model was chosen because it uses presence-background data (i.e., randomly selected absences from areas that have been accessible to the species).It generally has a better performance than presence-absence models, but absence data may not be reliable.A species may go undetected or it may not have had sufficient time to disperse to new locations yet [32].Another advantage of MaxEnt relevant to this study is that it can handle small samples of a species' presence data [33].
MaxEnt predicts habitat suitability as a function of environmental variables and species occurrence data.This habitat suitability is represented by a scale ranging from 0 (low suitability) to 1 (high fitness) [32,[34][35][36].Proper calibration and evaluation are necessary to reduce the complexity of the model [37], considering the choice of: i) accessible area (background or M area), ii) the type of variables that MaxEnt constructs (features), and iii) the type of model output (raw, cumulative, or logistic), as these considerations affect the inferences to be made [24].Proper calibration and evaluation is also especially important for data sets suffering from sampling bias, and for studies that require transfer models through space or time [38,39].
In this study, the calibration and evaluation method for A. coahuilense modeling were carried out using the library "ENMeval" [40] in the statistical software R 3. 1.3 [29].This study considered the 28 geographic records located within the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] and the six climatic variables identified above, in addition to the variable of altitude.The calibrated model was evaluated by calculating the coefficient standardized Akaike information criterion (AICc).The AICc provides information on the relative quality of a model [37].Because the AICc is calculated using the data set, it is not affected by the method chosen for the data partition [40].The model with the lowest AICc was selected as the best fit for the species.
The information obtained from the calibrated model was projected to Mexico, considering the environmental variables described above and the soil layers of litosol and calcareous regosol, using the software MaxEnt [30].There were 100 repetitions performed [41] to obtain an ecological niche model, geographically represented as a map of habitat suitability under current climatic conditions for A. coahuilense.This model was converted to an absence/presence map using maximum test sensitivity plus specificity, and considering only areas of presence that coincided with the distribution of rosetophilous desert scrub [42] and limestone soil [43].This species requires the presence of this type of soil and vegetation (4).For the validation of the model, the value of AUC (Area under the ROC curve) was considered.Also, records located outside the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] were used to validate the model.
The importance of each bioclimatic variable in the observed distribution of A. coahuilense was evaluated according to the relative importance of each variable, which was obtained by adding the percentage of contribution (PC) and the importance of permutation (IP), evaluated by MaxEnt, and the result was divided by two (average contribution (PC + IP)/2) [44].

Climatic Variables
The following six climatic variables were retained from the bivariate correlation analysis: Annual Mean Temperature (bio1), Mean Diurnal Range (bio2), Max Temperature of Warmest Month (bio5), Precipitation of Driest Quarter (bio17), Precipitation of Warmest Quarter (bio18), and Precipitation of Coldest Quarter (bio19).The climatic ranges and distribution of A. coahuilense in the environmental context are shown in Table 1 and Figure 2.

Niche Overlap
The equivalence of the realized niche of zone A and zone B was significant (p < 0.05).The similarity test between the realized niche of zone A and zone B gave a value of D = 0.007 and was not significantly different (p = 0.148).Zone B, compared to zone A, gave a value of D = 0.007 and was significantly different (p = 0.009).In addition, the analysis shows that zone B, with respect to zone A, had a value of expansion of 0.921, a stability value of 0.078, and a niche unfilled value of 0.009 (Figure 3).

Niche Overlap
The equivalence of the realized niche of zone A and zone B was significant (p < 0.05).The similarity test between the realized niche of zone A and zone B gave a value of D = 0.007 and was not significantly different (p = 0.148).Zone B, compared to zone A, gave a value of D = 0.007 and was significantly different (p = 0.009).In addition, the analysis shows that zone B, with respect to zone A, had a value of expansion of 0.921, a stability value of 0.078, and a niche unfilled value of 0.009 (Figure 3).

Niche Modeling
The model had an AUC value greater than 0.95 (which indicated a model with good reliability), and predicted areas with habitat suitability at the geographic points used to validate the model.According to the analysis of relative importance, the elevation variable presented a high percentage of contribution and importance of permutation, which allowed it to obtain 86.9% in the relative importance for the distribution of A. coahuilense.On the other hand, the variable Precipitation of Coldest Quarter had lower scores than those of Precipitation of Warmest Quarter, with respect to the contribution percentage (19.9%and 24.9%, respectively).However, the values of permutation importance were higher in Precipitation of Coldest Quarter than in Precipitation of Warmest Quarter (14.5% and 0.5%, respectively), which allowed it to obtain a 26.95% of relative importance.The variable with the lowest relative importance was Max Temperature of Warmest Month (Table 2).

Niche Modeling
The model had an AUC value greater than 0.95 (which indicated a model with good reliability), and predicted areas with habitat suitability at the geographic points used to validate the model.According to the analysis of relative importance, the elevation variable presented a high percentage of contribution and importance of permutation, which allowed it to obtain 86.9% in the relative importance for the distribution of A. coahuilense.On the other hand, the variable Precipitation of Coldest Quarter had lower scores than those of Precipitation of Warmest Quarter, with respect to the contribution percentage (19.9%and 24.9%, respectively).However, the values of permutation importance were higher in Precipitation of Coldest Quarter than in Precipitation of Warmest Quarter (14.5% and 0.5%, respectively), which allowed it to obtain a 26.95% of relative importance.The variable with the lowest relative importance was Max Temperature of Warmest Month (Table 2).The model indicated a habitat suitability area of 3909.778km 2 , distributed in the northeast of the state of Durango, south of Coahuila, north of Zacatecas, and northwest of Nuevo León.Within the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] the model indicated a habitat suitability area of 914.03 km 2 (Figure 4).

Discussion
The successful establishment and subsequent distribution of a species correlates with climatic conditions more strongly than with other factors, at least at wide spatial resolutions [11,[45][46][47].However, one species can occupy different climatic environments through its geographical distribution.For example, different analyses have reported that plant species with a wide distribution occupy different niches [48,49].Likewise, research of plant species across seven regions (located along an ~1800 km latitudinal gradient, from northern France to central Sweden and Estonia via Belgium, western and eastern Germany, and southern Sweden) showed that these species occupy different niches along a latitudinal gradient [10].In this regard, the occupation of different niches by a species may have important implications for understanding their vulnerability to the effects of climate change [14].
We found that, within the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6], the equivalence analysis shows that zone A and zone B were statistically significant; therefore, there is no evidence to accept the null hypothesis of equivalence proposed by Warren, et al. [16], enabling us to interpret that the evaluated niches are ecologically distinct.Niche overlap between zone A and zone B was low.Climatic similarity tests indicate that zone A did not present statistically significant differences when compared with zone B, while zone B showed significant statistical values compared to zone A. In this regard, it has been mentioned that the climate niche will be similar only if the test is significant in both directions [16].Under these criteria, the results of the similarity test allow us to reject the null hypothesis of climatic similarity proposed by Broennimann et al. [19], suggesting that the realized niche of A. coahuilense in zone A and zone B is climatically distinct.In this regard, most of the studies that indicate the occupation of different niches by a species have been documented in widely distributed species [50].Our results provide evidence that this characteristic can also occur in species that have a restricted distribution.

Discussion
The successful establishment and subsequent distribution of a species correlates with climatic conditions more strongly than with other factors, at least at wide spatial resolutions [11,[45][46][47].However, one species can occupy different climatic environments through its geographical distribution.For example, different analyses have reported that plant species with a wide distribution occupy different niches [48,49].Likewise, research of plant species across seven regions (located along an ~1800 km latitudinal gradient, from northern France to central Sweden and Estonia via Belgium, western and eastern Germany, and southern Sweden) showed that these species occupy different niches along a latitudinal gradient [10].In this regard, the occupation of different niches by a species may have important implications for understanding their vulnerability to the effects of climate change [14].
We found that, within the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6], the equivalence analysis shows that zone A and zone B were statistically significant; therefore, there is no evidence to accept the null hypothesis of equivalence proposed by Warren, et al. [16], enabling us to interpret that the evaluated niches are ecologically distinct.Niche overlap between zone A and zone B was low.Climatic similarity tests indicate that zone A did not present statistically significant differences when compared with zone B, while zone B showed significant statistical values compared to zone A. In this regard, it has been mentioned that the climate niche will be similar only if the test is significant in both directions [16].Under these criteria, the results of the similarity test allow us to reject the null hypothesis of climatic similarity proposed by Broennimann et al. [19], suggesting that the realized niche of A. coahuilense in zone A and zone B is climatically distinct.In this regard, most of the studies that indicate the occupation of different niches by a species have been documented in widely distributed species [50].Our results provide evidence that this characteristic can also occur in species that have a restricted distribution.
Previous studies have recognized niche shifts in an extensive variety of taxa, but causes of these changes remain poorly understood [20]; however, the differences in niche occupation of a species can arise through niche expansion or niche unfilling [51].The niche expansion is assumed a true niche change, while the unfilled niche is partially filled by the native niche in the invaded range [19].Our results indicate that A. coahuilense, within the realized niche of zone B, only covers 7% of the environmental conditions occupied in the realized niche of zone A. Regarding niche expansion, our results indicate high values in zone B (92%).This suggests that the absence of climatic similarity obtained in this study can be explained, to a greater extent, by the occupation of new climatic environments in the realized niche in zone B. On the other hand, it is considered that the unfilled niche seems to be larger for species introduced more recently and into a smaller number of locations, compared to those with ancient colonization histories introduced at several points in the geographic timeline [51].
Based on the aforementioned, our results indicate that the colonization potential of A. coahuilense in this region might be extremely low, which may explain the highly disjunct distribution in Coahuila, Durango, and Nuevo León.This can be explained by the historic biogeography of A. coahuilense.It has already been documented that species is hydrochoric [5], and its principal method of dispersion was through extinct paleolithic lakes and rivers of the Pleistocene [52,53].Likewise, it is possible to indicate that the occupation of new climatic niches by A. coahuilense occurs over thousands or millions of years.According to Webb et al. [54] and Ackerly [55], this characteristic corresponds to the species that presents with conservatism of its climatic niche.
The threat of global climate change to biodiversity can be viewed from a niche conservatism perspective [56,57].If the climatic tolerance of a species is not wide enough to encompass the new conditions or to acclimatize to them, the species with strong climatic niche conservatism must either migrate or go extinct, whereas, more evolutionarily labile species can potentially adapt [56,57].In this way, the presence of niche conservatism in A. coahuilense, in addition to the low dispersion capacity, could mean that this species does not have the capacity to adapt to the accelerated climate change that is currently occurring, increasing the risk of extinction of this species.As mentioned for several endemic species of mammals, birds, amphibians, reptiles, invertebrates, and plants [56,57].In this regard, the vegetation structure plays a key role in the activity, feeding, and distribution of some organisms [58][59][60][61].In the case of A. coahuilense, nurse plants provide microclimate to this organism, influencing its distribution [62].Thus, it is possible to point out that the viability of the species depends, to a large extent, on the conservation of its habitat that allows it to cushion the adverse effects of climate change.
On the other hand, the model generated from records inside the distribution area proposed by Fitz-Maurice et al. [3,6] predicts the populations of A. coahuilense that were recorded in outside areas, according to said area proposed by Fitz-Maurice.This indicates that these records are within the climatic niche considered in this study, and adding them will hardly improve the model.However, they are a good option to be used as validation criteria in testing the efficiency of the model.Likewise, this indicates that a better sampling effort was needed by Fitz-Maurice et al. [3,6] to take into consideration the new records of A. coahuilense proposed in this study.
The potential distribution models obtained from the present study allow us to identify potential zones based on the climatic requirements, elevation, and soil type of the species.Our results indicate a potential distribution area of 3909.778km 2 for A. coahuilense in northern Mexico.The estimated extent of occupation areas proposed by Fitz-Maurice et al. [3,6] is approximately 6000 km 2 .However, our model indicates that the potential range within the distribution of A. coahuilense proposed by Fitz-Maurice et al. [3,6] is 914.03 km 2 , that is, over 5000 km 2 is being overestimated.These differences may be due to the fact that these authors considered the hole polygon as habitat ocupated for the species.The area proposed from our results only considers areas with climate, soil, and adequate elevation for the presence of the species within the polygon that delimits this distribution.
Also, the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6] does not take into account records outside the considered range, which leads to underestimating the distribution of the species, whereas, the potential distribution model obtained in this study predicts records of A. coahuilense that were not considered by Fitz-Maurice et al. [3,6].Based on the above, it is possible to point out that it is necessary to use the implementation of new methods, such as the modeling of species distribution, to allow a more efficient evaluation of the known range of the species.However, it is important to consider that to develop potential distribution models, it is necessary to take into account the largest number of subsets of the climatic niche that a species may have, which make up the niche of the species.In this way, we will avoid, as much as possible, modeling only a subset of its climate niche.
In conclusion, we can show that A. coahuilense occupies distinct niches inside the distribution area proposed by Fitz-Maurice, et al. [3,6].However, the presence of niche conservatism cannot be ruled out.On the other hand, the potential distribution models obtained from the present study allow us to identify potential zones based on the climatic requirements of the species.This information is important to identify high priority areas for the conservation of A. coahuilense, as well as to identify new records of their presence, which may directly contribute to the conservation and management plans of this species.

Figure 1 .
Figure 1.Geographical records of Astrophytum coahuilense.In the study areas, the black triangles correspond to the records of zone A, the black circles refer to the records of zone B, the black boxes indicate the records of zone C, and the black stars refer to the points identified outside the distribution area of A. coahuilense proposed by Fitz-Maurice et al.[3,6].The Black dotted line shows the distribution proposed by Fitz-Maurice et al.[3,6], while solid black lines refer to the state limits.

Figure 1 .
Figure 1.Geographical records of Astrophytum coahuilense.In the study areas, the black triangles correspond to the records of zone A, the black circles refer to the records of zone B, the black boxes indicate the records of zone C, and the black stars refer to the points identified outside the distribution area of A. coahuilense proposed by Fitz-Maurice et al.[3,6].The Black dotted line shows the distribution proposed by Fitz-Maurice et al.[3,6], while solid black lines refer to the state limits.

Figure 2 .
Figure 2. Distribution of A. coahuilense in the climate space of Northern Mexico.The gray dots represent the multivariate climatic space in Northern Mexico, the red dots indicate the distribution of A. coahuilense in the climate space.

Figure 2 .
Figure 2. Distribution of A. coahuilense in the climate space of Northern Mexico.The gray dots represent the multivariate climatic space in Northern Mexico, the red dots indicate the distribution of A. coahuilense in the climate space.

Figure 3 .
Figure 3. Niche overlap.Panel 1 shows the evaluation frequencies A vs. B and B vs. A. Panel 2 presents the dimensions of the climatic niche of A. coahuilense in zone A and zone B, represented in a Principle Component Analysis (PCA).Grey shading shows the density of the occurrences of the species per cell.The dotted line illustrates 50%, and the solid line 100%, of the available environment.

Figure 3 .
Figure 3. Niche overlap.Panel 1 shows the evaluation frequencies A vs. B and B vs. A. Panel 2 presents the dimensions of the climatic niche of A. coahuilense in zone A and zone B, represented in a Principle Component Analysis (PCA).Grey shading shows the density of the occurrences of the species per cell.The dotted line illustrates 50%, and the solid line 100%, of the available environment.

Figure 4 .
Figure 4. Potential distribution model of A. coahuilense.In red, the areas with habitat suitability for the species are shown.The black dotted line is the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6].The state limits are indicated in black continuous lines.

Figure 4 .
Figure 4. Potential distribution model of A. coahuilense.In red, the areas with habitat suitability for the species are shown.The black dotted line is the distribution area of A. coahuilense proposed by Fitz-Maurice et al. [3,6].The state limits are indicated in black continuous lines.

Table 1 .
Distribution of A. coahuilense in the climatic environment.

Table 1 .
Distribution of A. coahuilense in the climatic environment.

Table 2 .
Relative importance of climatic variables in the distribution of A. coahuilense.

Table 2 .
Relative importance of climatic variables in the distribution of A. coahuilense.
Note: The value of each variable is the average contribution (contribution percentage + importance of permutation)/2.