Ecogeography of Dioscorea remotiflora Kunth Reveals Its Potential in Regions Already Affected by Climate Change

1. Introduction

Dioscorea remotiflora is a monocotyledonous, perennial, and climbing plant, with heart-shaped leaves, dioecious flowers, and seeds in axillary clusters. It belongs to the Dioscoreaceae family [1] and is one of the native species of México from the genus Dioscorea [2,3]; currently, it is considered an endemic species [4]. Although D. remotiflora is not a cultivated plant, its roots, tubers, and rhizomes are collected since prehistoric times to be used as food [5], they contain 85% carbohydrates, 7.35% proteins, 3.76% lipids, and 3.68% minerals (K, Fe, Na, and Mg), therefore, it is used as a healthy snack or even as a gourmet dish [6,7]. The tubers contain secondary metabolites such as steroids, saponins, and diosgenin, though in low concentration [7], thus, they are not used for the extraction of diosgenin. Currently, D. remotiflora is considered an underutilized species; being a factor that contributes to this, the lack of research on this species; Up to now, the reports regarding the description of the ecology, climatic adaptation, and potential distribution of D. remotiflora are bare and insufficient; this hinders the design of strategies for its conservation and optimal use [8,9]. Therefore, the objective of this research was to use the occurrence data of D. remotiflora, in its natural distribution areas, to conduct an eco-geographical analysis in order to elucidate the contribution of several ecological descriptors in determine its current distribution, identify its adaptation patterns, and develop optimal Maxent models of potential geographic distribution, through the use of the Kuenm R-package, which automates the creation, calibration, and evaluation of ecological niche models [10].

Plant growth and distribution are the results of the species’ response to the environmental complex that prevails in the occurrence sites [11,12], where aspects of climate, soil, vegetation, and others concur. As a species moves from its center of origin to other geographic regions, it finds it necessary to adapt to different environmental conditions. Whether the adaptation process is successful, then the species will have colonized new territories, extending its distribution, and thus, its adaptation environmental scope, which usually triggers the increase of its tolerance to abiotic and biotic stress [13]; this aspect may not be manifested uniformly in all the eco-geographic populations of the species [14], and, moreover, the novel traits may be the result of the action of genetic or epigenetic adaptation processes [15,16] with the possibility of simultaneous and complementary action of both mechanisms [16].

D. remotiflora is distributed mainly in the central, southern, and western regions of México [17], generally in dry deciduous tropical forests [18,19], which indicates its adaptive capacity to diverse environments [20]. This is why characterizing the ecogeography of a plant species considering their all accessions and populations is a key aspect.

With the presence of climate change in agricultural areas, and the increasing environmental stress for crops, food production must face the challenge with productive plant species that are resistant to drought and heat (main types of abiotic stress in México; [21]. However, in various regions of México, the crop pattern is reduced due to a short growing season with a low availability of rainwater [22]. In these cases, it is necessary to explore new crop species, for which a characterization study is required such as the one addressed in the present investigation for D. remotiflora.

2. Results

2.1. Selection of environmental variables

The Spearman correlation analysis enabled to reduce from 42 to 20 environmental variables to be used in further analysis [23,24]. Furthermore, preliminary tests with MaxEnt revealed that of the 20 variables, seven are the most relevant in determining the distribution of D. remotiflora.

2.2. Current distribution, climatic adaptation and ecological descriptors

Figure 1 shows the current distribution of D. remotiflora in the agroclimatic regions of México, as shown, this species is found predominantly in the areas near the Mexican Pacific coast, mainly in the central and southern portions. A greater presence of D. remotiflora is observed in the agroclimatic regions: Humid-subhumid semi-warm (98 accessions), Dry-subhumid semi-warm (88 accessions), Humid-subhumid warm (67 accessions, and Dry-subhumid warm (66 accessions). However, D. remotiflora is present in 15 of the 29 agroclimatic regions of México (Table 3), which enables this species to develop in thermal zones from semi-cold to very warm and in hydric zones from semi-arid to humid.

Table 4 shows the FAO soil units and the textural classes in which D. remotiflora is distributed. As can be seen, most of the occurrence sites of this species are distributed in the following soil types: Lithosol (144 accessions), Calcaric Regosol (98 accessions), Eutric Regosol (99 accessions) and Haplic Faozem (45 accessions).

Figure 2. Current distribution of D. remotiflora in soil units of México.

Figure 2. Current distribution of D. remotiflora in soil units of México.

Table 5 shows the ecological descriptors of D. remotiflora, which contain the environmental adaptation ranges (RAA) of the species (minimum and maximum values of environmental variables that allow the presence of D. remotiflora, and the optimal environmental range (RAO), which makes possible the higher frequency of occurrence sites of this plant. According to the information in Table 5, D. remotiflora is distributed in environments with an annual moisture availability index (MAI) ranging from 0.27 to 2.32, with an optimal range of 0.40 to 0.99, which corresponds to an RAA and RAO of an annual precipitation interval of 444 to 2886 and 700 to 1299 mm, respectively. This plant prefers areas where precipitation in the wettest quarter ranges from 400 to 884 mm, although it tolerates conditions from 240 to 1204 mm. This species can tolerate a long season (November-April) with low precipitation, even with an accumulation of only 23 mm in those six months (Table 5). The optimum growing season goes from 120 to 150 days, although it grows in regions with a growing season as long as 120-190 days.

Altitudinally, D. remotiflora is present from 6 to 4295 masl, but most of the presence sites occur between 200 and 800 masl. This encourages this species to develop in areas with an average annual temperature between 14.7 and 28.5°C, with an optimum of 19 to 27°C, an extreme monthly average maximum temperature of 41.2°C (Maximum maximorum temperature) and an extreme monthly average minimum temperature of 1.7°C (Minimum minimorum temperature). Both in the seasonal periods May-October and November-April, the interval in which the greatest number of accessions occurs is from 19 to 26°C, very similar to the optimum average annual temperature. Regarding soil properties, D. remotiflora is present in coarse, medium and fine textured soils, but prefers its presence is more abundant in coarse-textured soils, the number of occurrence sites decreases markedly in medium and fine-textured soils, indicating a clear preference of the species for soils with excellent drainage.

2.3. Modeling distribution niches of D. remotiflora

Due to Maxent not being exempt from the effects of collinearity and the fact that these can intervene in the estimation of the factors, as well as inferring the uncertainties when the models are transferred spatially and temporally [25], and with the purpose of obtaining an accurate model with a reduced number of variables, eliminating the possibility of obtaining an over-fit or over-parameterized model[10]; the Kuenm R-package allowed, through ecological niche modeling, optimizing the number of environmental variables, leaving only 7, due to their greater contribution to the presence and distribution of D. remotiflora. The results of the Jackknife statistical test of factorial importance reported that the most determining environmental variables in the presence and distribution of D. remotiflora are: precipitation of the warmest quarter (42.4%), mean precipitation of the driest month (17.5%), minimum temperature of the coldest month (15%), November-April mean solar radiation (10%), annual mean relative humidity (8.5%), annual moisture availability index (5.7%), and May-October mean temperature (0.9%) (Table 6).

Based on the results obtained through the Jackknife test, three sets of cases were analyzed: "only with variable", "without variable" and "with all variables" [26], thus revealing the effects of environmental variables in the appropriate range of D. remotiflora (Figure 3). Among these environmental variables, it was found that solar radiation from November to April was the most relevant, with a regularized training gain greater than 0.8. In addition, other variables of significant importance were identified, such as the precipitation of the warmest quarter, the annual humidity index, relative humidity, the minimum temperature of the coldest month, the precipitation of the driest month, and the average temperature from May to October. In all these cases, the regularized training gains exceeded 0.6.

Kuenm R-package enabled to obtain 372 models (Figure 4), all of them significant, however, only two models positively met the Akaike criterion (AIC=0) and a maximum omission rate of 5%. The Kuenm R-package implemented with MaxEnt determined that 3.0 was the optimum regularization multiplier. The model finally selected to depict the distribution of D. remotiflora was judged excellent since the AUC of the ROC curve accounted for 0.935.

Figure 5 shows the current and potential distribution of D. remotiflora based on thresholding the environmental suitability with the “Balance training omission, predicted area, and threshold value” criterion. The map in Figure 4 shows four types of areas: areas not suitable (white color), areas with low environmental suitability (light green), areas with medium suitability (dark green), and areas highly suitable for D. remotiflora presence (dark grey), which matches the largest number of occurrence sites. The total potential distribution area for D. remotiflora accounted 428,747.68 km², distributed as follows: The areas with low suitability account for 150,547.104 Km2 and are located in Sonora, north Sinaloa, Guanajuato, Querétaro, México State, coasts of Guerrero, center Oaxaca, Chiapas, Campeche, San Luis Potosí, south Tamaulipas, and Veracruz. The areas with medium suitability cover 127,817,096 km2 and locate in southern Sonora, southern Chihuahua, northern Sinaloa, the coasts of Nayarit, Jalisco, Michoacán, Guerrero, Oaxaca, and Chiapas. The highly suitable areas resulted in 150,383.48 km2 (dark gray), and are located in Sinaloa, Nayarit, Jalisco, Colima, Michoacán, Morelos, little portions of Puebla and Tamaulipas, northern Guerrero, and the coasts of Oaxaca. As can be seen in Figure 4, there are potential areas in regions where the presence of D. remotiflora has not yet been reported, which indicates territories that could be the target of future exploration of new populations. This is the case of the Peninsula of Yucatán, center Veracruz, southeastern San Luis Potosí, southern Tamaulipas, Guanajuato, northeastern Jalisco, northern Sinaloa, northern Hidalgo, and western México state.

3. Discussion

3.1. Current distribution, climate adaptation and ecological descriptors

Most of the D. remotiflora occurrence sites are concentrated amongst the dry-subhumid and humid-subhumid agroclimatic regions with the warm and semiwarm variants, with an annual mean temperature between 19 and 27°C and an annual moisture availability index (MAI) between 0.4 and 0.99 (Table 4). These results agree with the ones reported by [2] who indicate that D. remotiflora is distributed in areas with semi-warm, warm and sub-humid climates, ranging from the northern portion to the central region of México. However, an interesting aspect is the climatic extremes in which D. remotiflora can adapt, meaning that even when most of the occurrence sites of this species are located between 700 and 1299 mm of annual rainfall, there are populations that subsist with 444 mm per year, and on the other extreme, other populations do it in sites with 2886 mm of precipitation per year. When combined with annual potential evapotranspiration data, these precipitation values translate into MAI values from 0.27 to 2.32, which according to arid zones scheme from UNEP (Middleton & Thomas, 1997; adapted by Ruiz et al., 2004); match with semiarid to very wet lands. Regarding this, [28] mentions that some species of Dioscorea adapt to dry periods and can survive under conditions of water deficit better than many other species and crops. This seems to be the case of D. remotiflora, according to the ecological indicators obtained. The presence of populations of D. remotiflora in extremely humid sites (2886 annual mm and MAI > 2) could be explained in sites with excellent soil drainage or to consider those data as an indicator of the ability of this species to colonize habitats typically unsuitable for its development [18].

Regarding the variable precipitation of the warmest quarter, which was the most significant for the presence of D. remotiflora in the MaxEnt modelling, Table 4 indicates a range of 240 to 1,204 mm, with an optimal range of 400-884 mm. Considering then that D. remotiflora is present in sites with a minimum annual rainfall of 444 mm, and that, from those millimeters at least 240 must occur in the warmest quarter [29](Avalo and García, 2021), we can conclude that this plant has the ability to adapt and develop in environments with an irregular distribution of precipitation, which may also represent a comparative advantage of D. remotiflora in relation to other species. This also leads to conclude that D. remotiflora is a species that adapts more to moderate to low humidity conditions, thus, the very humid orarid environments are not conducive to high productivity of this species. The results obtained also coincide with previous reports related to the adaptation of the genus Dioscorea to tropical and subtropical zones, with tolerance of some species to conditions of water deficit [28,30].

The environmental characterization of the presence sites of D. remotiflora shows that there is a greater number of accessions of this plant in Lithosol, calcaric Regosol and eutric Regosol soils, which do not offer the best conditions for the development of vegetation and crops [31]. According to [32], Leptosols (Litosols and Rendzinas) represent 28.3% of the Mexican territory and are characterized by very thin, stony and poorly developed soils, that can contain a large amount of calcareous material, which immobilizes mineral nutrients. They are common in mountainous areas and on shallow limestone plains. Their agricultural potential is limited by their shallow depth and high compaction, which makes them difficult to work on. On the other hand, Regosols are considered very young soils that develop on unconsolidated material, light in color and poor in organic matter. They are common in arid, semi-arid (including the dry tropics) and montaneous regions, they can be found associated to Leptosols and with rock or tepetate outcrops. From the above, it can be deduced that D. remotiflora has a great adaptive capacity to poor soils, with non-optimal agroclimatic conditions for the rest of the species, which may be an attribute that makes this species an alternative for cultivation in regions where climate change is deteriorating the environmental conditions of agricultural production systems [33].

Regarding soil texture, according to Table 4, D. remotiflora prefers coarse-textured soils, typical of Lithosol and Regosol soils, which do not store a large amount of water. Unlike other species of the Dioscorea genus, D. remotiflora is susceptible to tuber putrescence [18], therefore it requires well-drained soils. The presence of D. remotiflora detected in medium-textured soils is related to the occurrence of lower annual precipitation levels (735.6 mm on average) than those of sites with coarse-textured soils (1034 mm on average), which ensures that even in soils that store more moisture, it is possible for this species to adapt, as long as the volumes of precipitation are not high; this compensatory effect has been previously reported for diverse crop species [34].

According to table 4, D. remotiflora can tolerate an extreme monthly average minimum temperature of 1.7°C and, on the other hand, it can survive an extreme monthly average maximum temperature of 41.7°C, with annual thermal oscillations ranging from 10 to 20° C and with an optimal annual thermal range of 13 to 16°C, fact that shows the wide range of thermal conditions in which this plant can survive, and this includes temperature regimes that are classified as very extreme [35]. Other species of the genus Dioscorea have shown tolerance to extreme thermal environments, such as D. divaricata, which can tolerate temperatures as low as -18°C. The occurrence of these extreme minimum temperatures, however, causes delayed full maturity and reached until 3 to 4 years after [36].

On the other hand, extreme temperatures are generally considered a source of dormancy in postharvest tubers and seeds [37]. [28] mentions that the dormancy of the tubers of the Dioscorea species lasts 120 days, which limits their agricultural production; this indicates an aspect that should be worked on in the immediate future to make D. remotiflora a more promising agricultural species. [37] report that tuberization is induced by environmental cues such as short days, low temperatures, and higher soil moisture content.

Another important aspect in the development of this species is photoperiod since it intervenes in the formation and growth of leaves and tubers; there are differences among Dioscorea species in relation to their response to photoperiod. In the case of D. remotiflora, in the long-day season, which is from May to October, foliar growth is favored, and in short days, the growth and swelling of the tuber are stimulated, promoting the production and storage of starches [38]. In Table 4, it can be observed that the optimal range of photoperiod for exhibiting adequate leaf growth is from 12.6 to 12.9 h, and for good tuber development the optimal range is from 11.10 to 11.39 h.

3.2. Modeling of distribution niches of D. remotiflora

The Kuenm R-package implemented with MaxEnt, allowed to obtain an optimal niche model to appropriately depict the D. remotiflora distribution in México; based on the requirements of statistical significance, the optimal regularization multiplier, the feature classes, and the omission rate established, the analysis process with Kuenm produced two possible models (Figure 4). One of them was selected, which can be considered a good decision since it fulfills the requirements established [39], and the AUC of the ROC curve accounted for a value greater than 0.93 [40]. The AUC value is an important tool to assess model performance, the higher the AUC value (closer to 1) the better the model performance is [41,42]. The Jackknife test identified the most determining variables for the presence of D. remotiflora, the precipitation of the warmest quarter, precipitation of the driest month, minimum temperature of the coldest month, November-April solar radiation, annual mean relative humidity, annual moisture availability index, May-October mean temperature. However, the Jackknife test indicates that the most important variable is solar radiation from November to April (Figure 3). These results agree with what was reported by [43], who mentions that the production of diosgenin, corticosteroids, carbohydrates, and other compounds produced by these species are subject to solar radiation. On the other hand, these results agree with those reported by Viruel et al. (2014) who pointed out that the environmental variables of major influence for the presence of Dioscorea humilis are the precipitation of the wettest quarter, the precipitation of the warmest quarter, and the precipitation of the coldest quarter. On the other hand, [43] report that for 10 species of the Dioscorea genus, the variables that are most important and intervene in the production of secondary metabolites are annual precipitation and average annual radiation.

The capabilities observed in the present research by D. remotiflora, to develop in adverse environments, could be explained through the ability to adjust its metabolism and modify its morphological characteristics, such as the size and thickness of the leaf, allowing it to adapt to extreme conditions, where other species fail to thrive [44,45].

The current distribution of D. remotiflora obtained in the present research (Figure 4), agrees with that reported by [46], who mentions that this species is distributed in Chiapas, Chihuahua, Colima, Durango, Guanajuato, Jalisco, Oaxaca, Guerrero, Michoacán, Guerrero, Nayarit, Puebla, Tamaulipas, Tabasco, and Zacatecas; which represents areas along the foot of the Sierra Madre Occidental, the Sierra Madre del Sur and its confluence with the Transversal Neovolcanic Axis, where coniferous and oak forests are found, these two, are home to herbaceous and forest communities with presence of endemism [47]. However, the distribution, diversity, and structure of populations are strongly influenced by historical, geographical, and climatic events [48], which also trigger speciation processes [49,50]. These arguments explain the number of species of the genus Dioscorea and the differences between them in terms of environmental ranges. On the other hand, the areas of high environmental suitability are located on the Mexican Pacific coast; this location suggests that optimal conditions of agroclimatic variables exist in these areas for the development of this plant. In contrast, in the areas of medium and low environmental suitability, it is possible that the species is adapting to agroclimatic ranges that go from very arid to very humid, moving to the central part of México. These areas, therefore, have the potential to establish crops of the species. However, it is important to note that the classification of environmental suitability areas may be modified by the effects of climate change [51], thus, a complementary research regarding potential distribution áreas under climate change scenarios would be required.

3.3. D. remotiflora cultivation prospects

Currently, D. remotiflora is a plant that is exploited for human consumption, but basically through the collection of tubers in its natural habitat (Figure 5). However, it presents favorable characteristics to integrate into regional crop patterns, given its comparative advantages regarding other plant species, such as its adaptation and production in poor, shallow, and infertile soils, as well as its tolerance to drought, since it does not have a high water requirement for its development. Besides, some characteristics inherent to the species should be taken into account before considering it, as a cultivation option, such as the fact that sexual propagation is not considered a viable alternative up to now, due to the fact that the seeds it produces are attacked by pests, damaging 60 to 80% of them; thus, the health of the seeds should be ensured through pest control or opting for asexual propagation through the sprouting of tubers. However, the tubers have a prolonged dormancy of up to 120 days, which is a limitation to program their propagation. This leads to the need of exploring some of the available techniques to reduce the duration of dormancy and accelerate the germination process.

4. Materials and Methods

4.1. Occurrence data

1,030 geo-referenced accessions from herbaria, floristic inventories, scientific articles and databases were identified, later on, they were reviewed to eliminate repeated records, records with poor geographic coordinates, and those out of the study area [52]. These strategies were used to avoid considering accessions that correspond to introductions outside the natural areas of distribution [53]. Finally, 480 records were selected (Table 1).

4.2. Climatic data

Monthly, quarterly, seasonal, and annual rasters of precipitation, temperature, solar radiation and relative humidity were used to determine potential distribution areas of D. remotiflora. These climatic data were obtained from the Agroclimatic Information System for México-Central America (SIAMEXCA) [21]. The raster images have a resolution of 30” of arc and correspond to the period 1961-2010. From the SIAMEXCA rasters, other additional variables were generated adding a total of 42 environmental variables (Table 2).

4.3. Environmental characterization of the occurrence sites

Based on the geographic coordinates of the D. remotiflora occurrence sites, an environmental characterization of these sites was carried out using the 42 variables mentioned in Table 2. For this, data extraction procedures were done, out from the raster images, using the ArcMap software version 10.8 [54]. With data extracted from the 42 variables, an environmental data matrix (EDM) was built in Microsoft Excel.

4.4. Selection of environmental variables

Prior to the execution of the statistical analysis, the Shapiro-Wilk test was applied to verify the data normality, not finding normality (P<0.05) for any of the data series within the 42 environmental variables included in the EDM.

Diverse studies have shown that multicollinearity is a problem that can cause high correlations among independent variables, a fact that can lead to unreliable and unstable estimations of the regression coefficients [55]. For determining the presence of multicollinearity, a Spearman correlation coefficient r > 0.9 was established as a threshold value [39]. In this way, the correlated variables with an r <0.9 coefficient were selected, and among the variables with collinearity, the one considered the most relevant for the presence of the species was selected. Data from the EDM were used to perform the correlation analyses; this statistical analysis was carried out with programs developed in the R software version 4.05 [56] and with normalized data. The result of these analyses reported 20 useful variables, later Preliminary analyses were carried out in MaxEnt individually and in conjunction with the Kuenm R-package in order to use the Jackknife test to select the most relevant variables in the distribution of D. remotiflora and thus carry out the final modeling of the ecological niche.

4.5. Characterization of the adaptive capacity of D. remotiflora.

Based on the geographic location of the D. remotiflora occurrence sites, the agroclimatic regions, soil units, and vegetation types where the species is currently distributed were characterized. For this purpose, the map of agroclimatic regions of México and Central America was used [13], as well as the soil units map and the vegetation types map for México [57]. This made possible the elaboration of a list of the edapho-climatic conditions to which D. remotiflora currently adapts and its preference for certain hábitats [58.

In addition, the ecological descriptors for D. remotiflora were determined by using the EDM information, but only taking into account the 20 environmental variables that were selected after the correlation-colinearity analysis. The ecological descriptors were established in terms of environmental ranges for D. remotiflora adaptation, and optimal environmental ranges for D. remotiflora presence, which corresponded to the highest frequency of occurrence sites.

4.6. Ecological niche modeling

Ecological niche modeling was done by using Maxent model, which uses the principle of maximum entropy with species presence and environmental data to create a correlative model that relates the ecological requirements of a species with the regional environmental availabilities to predict the relative habitat suitability, also it allows deriving specific descriptors for enriching the ecological characterization of the territories [59,60,61]. We used the Kuenm R-package [10,56,62] to automate and optimize the Ecological Niche Modeling (ENM) process. Kuenm's Kuenm_ceval function creates preliminary models [63] with occurrence data and environmental predictors, and also evaluates the efficiency of these models through their statistical significance with the cal_eval function, in addition, the relative quality of the model is evaluated by using the Akaike Information Criterion (AIC), corrected for small sample sizes (AICc) [10]. Kuenm also enables to assess diverse regularization muliplier (RM) factors, combinations of feature classes (FC) and different groups of environmental predictors, as well as to establish the allowable omission rate (OR). For each parameter setting, two models are created: one based on the complete set of occurrences, and the other based only on the training data. For this research, models were tested using a sequential order of the FC (L, LQ, H, LQH, LQHP, LQHPT) and RM values of 0.1 to 5 with 0.1 increases, a maximum omission rate of 5%, and run 50 k-fold replicates of each configuration; 500 iterations were used [64]. The D. remotiflora occurrence sites and the ASCII files for the 20 variables selected after correlation-collinearity analysis were used as inputs in the ENM process with Kuenm R-package.

5. Conclusions

D. remotiflora is an endemic species from México, whose tubers are collected in their natural habitat for human consumption. D. remotiflora is mostly adapted to semi-warm to warm environments and to semi-arid to sub-humid climates with a long dry season. However, the whole range of environmental adaptation of D. remotiflora also includes extremely humid, temperate, and very warm habitats, which indicates a high adaptive capacity of this species. In addition, D. remotiflora has a greater presence in regions where very thin, stony, poorly developed, and low-fertility soils predominate. All these characteristics make D. remotiflora a species with comparative advantages to develop in edapho-climatic environments that are adverse for most crop species. However, the greater presence of D. remotiflora occurrence sites indicates optimal environmental conditions, under which the development and productivity of the species would be expected to be maximum. In this way, populations of D. remotiflora proliferate in soils of coarse texture, in low to midlands, with annual precipitation of 700 to 1300 mm, and annual mean temperature of 19 to 27°C.

D. remotiflora is currently distributed in the western portion of México, along regions bordering the Mexican Pacific, which hints at its center of geographic origin. Niche modeling identified precipitation of the warmest quarter, precipitation of the driest month, mean minimum temperature of the coldest month, November-April mean solar radiation, Annual mean relative humidity, Annual moisture availability index, and May-October mean temperature as the variables with the greatest contribution to explain the presence of D. remotiflora. Furthermore, the Kuenm R package enabled the selection of a niche model that optimized the depiction of the potential distribution areas for this species. Thus, potential areas with high environmental suitability for D. remotiflora were located in the Mexican States where this species is already present, and potential areas with low to mid environmental suitability were identified in regions with the current presence of D. remotiflora as well as regions where it has not yet been reported, as in the Yucatan Peninsula, Northeast region, and the Gulf of México.

The cultivation perspectives of D. remotiflora are favorable considering its capacity to adapt to harsh environments, and that its nutritional and medicinal properties are valuable, however, the prolonged dormancy of its tubers is one of the intrinsic aspects of this plant that limits its offer as a cultivation option.

The information generated in this research will make it possible to learn more about the distribution, adaptation, and ecology of this species, which can also contribute to the design of conservation strategies for D. remotiflora and to the planning and selection of potential areas for its cultivation.

Currently, plant genetic resources worldwide are facing the pressure of overexploitation and environmental change, resulting in habitat fragmentation and biodiversity threats. To ensure the permanence of these natural resources, it is essential to maintain or increase the resilience of these systems against these pressures (Ten Broeke et al., 2019; Gülçin, 2021). Therefore, an efficient conservation management model must be adopted to address these changes and adequately preserve this plant genetic resources.