Current and Potential Spatial Distribution of Six Endangered Pine Species of Mexico : Towards a Conservation Strategy

Mexico is home to the highest species diversity of pines: 46 species out of 113 reported around the world. Within the great diversity of pines in Mexico, Pinus culminicola Andresen et Beaman, P. jaliscana Perez de la Rosa, P. maximartinenzii Rzed., P. nelsonii Shaw, P. pinceana Gordon, and P. rzedowskii Madrigal et M. Caball. are six catalogued as threatened or endangered due to their restricted distribution and low population density. Therefore, they are of special interest for forest conservation purposes. In this paper, we aim to provide up-to-date information on the spatial distribution of these six pine species according to different historical registers coming from different herbaria distributed around the country by using spatial modeling. Therefore, we recovered historical observations of the natural distribution of each species and modelled suitable areas of distribution according to environmental requirements. Finally, we evaluated the distributions by contrasting changes of vegetation in the period 1991–2016. The results highlight areas of distribution for each pine species in the northeast, west, and central parts of Mexico. The results of this study are intended to be the basis of in situ and ex situ conservation strategies for the endangered Mexican pines.


Introduction
Pines are the most significant members of the Pinaceae family; there are 113 species widely distributed around the world, most of them in the Northern Hemisphere [1].Mexico has the highest species diversity of pines (46 species), with contrasting geographical and intraspecific genetic patterns, because of adaptive responses to climate change in the past [2].The country hosts 29 endemic species, distributed in different ecological environments [1].In Mexico, pines are usually found in three types of communities, pure pine forests, pine-oak forests, and oak-pine forests, that occur at different elevations, climates, and exposure conditions [3], as representatives of the high ecological diversity of the Mexican landscape.Mexican temperate forests are mainly located in high altitudes over mountainous zones originating from volcanos [4]; however, pine communities can also be found in different vegetation types such as grasslands and scrublands, representatives of arid and semi-arid zones of Mexico [5].
Many coniferous forests have suffered a decrease in area, disappeared, or changed composition, allowing for substitution by angiosperms [5].The causes of the reduction of forests include the constant increase in the frequency of fires, due to both natural causes and human-made phenomena Forests 2018, 9, 767 2 of 21 and excessive logging, whether for the direct harvesting of timber and by-products or the expansion of agricultural and livestock activity.Also, recent trends in global warming pose a threat to many forest species [6], and pine species in particular are more vulnerable because of their specific environmental requirements and the limited distribution of their populations.Habitat fragmentation has straightforward consequences on plant populations as it significantly increases the risk of local extinction.Moreover, fragmentation leads to a reduction in effective population size and an increase in the degree of isolation between populations [7,8].The immediate effect of fragmentation on remaining populations tends to increase the degree of inbreeding as crosses between related individuals as well as self-fertilization increase would compromise the viability of the species [7].For instance, there is a need to develop research related to the status of endangered and endemic Mexican pine species populations.
Within the diversity of pines in Mexico, Pinus culminicola Andresen et Beaman, P. jaliscana Perez de la Rosa, P. maximartinenzii Rzed., P. nelsonii Shaw, P. pinceana Gordon, and P. rzedowskii Madrigal et M. Caball.are six Mexican pines catalogued as threatened or endangered species [9] due to their restricted distribution and low population density.These pines occupy diverse habitats, present a variety of ecological roles, and cover a wide range of environments differing in altitude, precipitation, temperature, and soil conditions [10][11][12][13].Spatial modeling and GIS techniques provide the tools to create current and future scenarios of possible suitable areas of distribution according to available databases and information about these species.
Spatial modeling and niche ecological models combined with other analytical tools (e.g., GIS) allow the evaluation of non-inventoried sites as well as the modeling of past or future scenarios in the distribution of species [14].Recently, numerous modeling methods and tools have been developed [15,16], although they are mainly used in ecology and biogeography [17,18].The use of these alternatives accomplishes a double function: first, they provide knowledge about the potential distribution of the species to know the richness and diversity of non-evaluated areas.Secondly, they use said predictions in the choice of sites of particular interest as biological conservation zones [19,20].
In this paper, using spatial modeling, we aim to provide up-to-date data on the spatial distribution of six endangered and endemic pine species of Mexico according to different historical registers coming from different herbaria distributed around the country.The objectives of this study were: (a) to recover historical observations of the natural distribution of each species, (b) model suitable areas of distribution according to environmental requirements and (c) evaluate the model predictions to highlight "hotspots" or areas of possible concern for the conservation of each pine species.All pine species selected are described as "in danger of extinction" by the Ministry of Natural Resources of Mexico.The geographic locations of each pine species were used to explain the distribution due to different conditions through a logistic regression approach and to highlight the most critical environmental factors that drive such distribution.

Materials and Methods
This study was divided into two steps.First, we conducted a literature review in scientific databases and performed visits to different herbaria located across the country to collect and retrieve historical information about the location and environmental characteristics of each pine species.Second, the potential distribution was modeled through a multiple logistic regression [21], using the information recovered from the literature and herbarium databases.Logistic regression modeling is a presence/absence-based method used here to model the potential distribution of six pine species.We estimated the probability of occurrence of several predictor variables coming from geographic layers and environmental and climatic databases.Later, the current and potential distribution of each pine species were compared against the official land cover cartography to validate the estimations and discuss whether land cover changes are having a negative impact on the current and potential distribution of the endemic pines of Mexico.

Study Area
The study area comprised the northern and central part of Mexico and included the States of Durango, Coahuila, Nuevo Leon, Tamaulipas, San Luis Potosi, Jalisco, Queretaro, Hidalgo and Michoacán.The study area mainly focused on temperate coniferous forests but also included other vegetation types such as grasslands, scrublands and secondary vegetation, according to each pine species distribution reported (Figure 1).We focused on the reports of distribution of six endangered and endemic pine species of Mexico: Pinus culminicola, P. jaliscana, P. maximartinezi, P. nelsonii, P. pinceana, and P. rzedowskii.These species were selected because, among all the Mexican pines, they are the only ones that are catalogued as endemic and endangered according to Mexican legislation (Table 1).Also, they represent a variety of bioclimatic conditions and vegetation types that can serve as a basis for further research on modeling species distribution.The study area comprised the northern and central part of Mexico and included the States of Durango, Coahuila, Nuevo Leon, Tamaulipas, San Luis Potosi, Jalisco, Queretaro, Hidalgo and Michoacán.The study area mainly focused on temperate coniferous forests but also included other vegetation types such as grasslands, scrublands and secondary vegetation, according to each pine species distribution reported (Figure 1).We focused on the reports of distribution of six endangered and endemic pine species of Mexico: Pinus culminicola, P. jaliscana, P. maximartinezi, P. nelsonii, P. pinceana, and P. rzedowskii.These species were selected because, among all the Mexican pines, they are the only ones that are catalogued as endemic and endangered according to Mexican legislation (Table 1).Also, they represent a variety of bioclimatic conditions and vegetation types that can serve as a basis for further research on modeling species distribution.

Pine Species Information
The spatial distribution and endangered status of each pine species are described in Table 1.All of the species below are listed as endemic to Mexico [9].Observations of presence of each species were collected from historical reports associated to a geographic reference.Meanwhile, observations of absence were randomly selected from areas where the vegetation type and conditions matched the distributions described in literature, but no reports of the species were found.
The historical information on pine samples across Mexico was retrieved by visiting different herbaria from different academic and research institutions (Table 2).The goal was to retrieve geographic locations for each pine species to construct a database that could serve as a basis for further modeling.Furthermore, an electronic search was also conducted to explore databases from different institutional websites.The websites consulted belong to the world network on biodiversity (REMIB), the virtual herbarium of the Biodiversity Conservation Commission (CONABIO; HVC), the Gymnosperm Database, the herbarium network of the northeast of Mexico (Herbanwmex), the Flora of Baja California, Global Biodiversity Information Facility (GBIF), Southeast Regional Network of Expertise and Collections (SERNEC), and the consortium of Midwest Herbaria.Additionally, we explored and analyzed the database of the national forest inventory provided by the National Forestry Commission (CONAFOR) of Mexico.
All the information was filtered and analyzed to avoid inconsistencies.We focused on natural populations.Thus, all the data coming from locations belonging to arboretums or plantations were discarded to avoid bias and misinterpretation during the modeling process.

Geographic Layers
The National Institute of Statistics and Geography of Mexico (INEGI) produces and publishes land cover and vegetation type maps at the national level at a scale of 1:250,000 using a 25-hectare minimum mapping unit [22].The first land cover map was published in 1991 and the latest one in 2016.Therefore, we used the first and sixth version of the land cover map to provide a spatial and temporal context for the analysis.A national Digital Elevation Model (DEM) raster dataset with approximately 30-m pixel size was used to create slope and aspect layers for further analysis.All digital geographic data described before are freely available (www.inegi.gob.mx) and were provided by INEGI.
Climatic data layers were retrieved from the Digital Climatic Atlas of Mexico (DCAM).DCAM was developed by the Informatics Unit for Atmospheric and Environmental Sciences (abbreviated as UNIATMOS in Spanish) of the National Autonomous University of Mexico (abbreviated as UNAM in Spanish) available at http://uniatmos.atmosfera.unam.mx/ACDM/.The bioclimatic parameters used in this study (Table 3) for the continental surface of Mexico were calculated per station from the 1902-2011 daily climatological database of the National Meteorological Service according to the ANUCLIM methodology [23].Additional bioclimatic variables were collected from WorldClim [24], an open climate dataset for ecological modeling and GIS (Table 3).WorldClim is a dataset of global Forests 2018, 9, 767 5 of 21 climate layers (gridded climate data) with a spatial resolution of about 1 km 2 , available at www. worldclim.org.Table 3. Bioclimatic variables used in this study [24].

Species Distribution Modeling
All pine observations were divided into two groups: 70% of the data were used to develop the models and the remaining 30% were retained to validate the models independently.
We used correlation matrices between the predictors to indicate if there were correlations/ associations between predictors.To reduce multicollinearity issues, only bioclimatic regions that had a low correlation between each other were used in the modeling.All data analyses carried out in this study were performed using the R software [25] packages "CAR" [26], "PerformanceAnalytics" [27], and "dplyr" [28].

Logistic Regression
To model the spatial distribution of the species, logistic regression was selected because it assumes the probability distribution as the response variable, and hence for the error terms from the fitted model, is adequately described by the random component chosen [21].In this study, we considered the response variable Y as the Bernoulli random variable with parameter E(Y) = π, taking only the values of 1 (when the species is present) and 0 (when the species is not present) with probabilities π and 1 − π, respectively.Then, the logistic regression model is Y = E(Y) + ε, and E(Y) = π is expressed as [21]: where, if we consider p − 1 predictor variables, vectors β and X are expressed as: Forests 2018, 9, 767 6 of 21 and Then, we have the linear predictor as: Numerical search procedures were used to estimate the values of β 0 , β 1 , . . ., β p−1 that maximize the log-likelihood.These maximum likelihood estimates are denoted by b 0 , b 1 , . . ., b p−1 .Therefore, the fitted logistic response function can be expressed as follows: where: The estimated odds ratio for predictor variable X k , denoted by π, is exp(bk).For continuous predictors this is the change in the ratio of probability of the species of pine occurring at any given location; when there is a unit change at X, we assume that all other predictor variables are held constant.

Validation of the Models
As measures of predictive models, the sensitivity, specificity and concordance indices were calculated.These indices were calculated on the validation sample that was not used to train the models.Finally, the overall accuracy was calculated through a confusion matrix.Also, the beta coefficients, standard errors, z values, and p values were retrieved for model diagnosis.Additionally, the receiver operating characteristic curve (ROC) was traced as a measure of fit.

Identifying Areas of Possible Conservation Concern
Once the current and potential distribution maps were prepared, we used the land use and land cover maps that the Mexican government provides (Series I, published in 1973, and VI, published in 2016) to evaluate changes in the land cover where the endangered Mexican pines are distributed.The main idea was to highlight possible areas of concern by contrasting areas where vegetation is drastically changing against the distribution of pine trees.First, we homogenized land cover classes between epochs.Then, both layers were transformed to raster files and we performed a bi-temporal analysis using the land change modeler of the TerrSet ® Version 18.11 (Clark Labs, Clark University, Worcester, MA, USA).The difference between layers was measured directly from values of the pixel image.The expression of image differencing is as follows: where I 1 and I 2 are layers from time t 1 and t 2 and (x, y) are coordinates and I d is the difference image.Pixels with no change are distributed around the mean, while pixels with change are distributed in the tails of the distribution curve.Since change can occur in both directions, it is therefore up to the analyst to decide which image to subtract from which [29].All spatial analyses using vector files were carried out using the software Arc GIS ® Version 10.6 (ESRI, Redlands, CA, USA).

Pine Species Distribution
Pinus culminicola, P. rzedowskii and P. maximartinezii presented minor presence observations, while P. jaliscana, P. pinceana, and P. nelsonii presented the most substantial presence observations.This is due to the number of localities where these species live.The observations of pine species had a description of the location, geographic coordinates, validation of the observation by a taxonomist and morphological description.
From the National Forest Inventory provided by the National Forestry Commission of Mexico, only observations of Pinus pinceana were reported.A total of 722 registers of the presence of endangered Mexican pines were used to define current distribution.Once the information was prepared, the natural distribution of Mexican pines was mapped out (Figure 2).

Pine Species Distribution
Pinus culminicola, P. rzedowskii and P. maximartinezii presented minor presence observations, while P. jaliscana, P. pinceana, and P. nelsonii presented the most substantial presence observations.This is due to the number of localities where these species live.The observations of pine species had a description of the location, geographic coordinates, validation of the observation by a taxonomist and morphological description.
From the National Forest Inventory provided by the National Forestry Commission of Mexico, only observations of Pinus pinceana were reported.A total of 722 registers of the presence of endangered Mexican pines were used to define current distribution.Once the information was prepared, the natural distribution of Mexican pines was mapped out (Figure 2).

Spatial Distribution of Pines Species According to Logistic Regression
By using geographic information, we developed models that describe in spatial context the probability of the presence of each pine species evaluated.The next section describes the results for each species.We set a threshold of 75% (p > 0.75) as a way of representing the likelihood of presence for each pine.

Pinus culminicola
According to the results, the dwarf pine (P.culminicola) (Figure 3) is found in the east mountains of Saltillo city, and southeast of Monterrey City, as well as small areas in the southeast of Nuevo Leon state.The modeling suggests the mean precipitation of the coldest quarters (Bio19) as the most important environmental variable (Table 4).The area estimated as suitable for P. culminicola is about 421 km 2 .

Spatial Distribution of Pines Species According to Logistic Regression
By using geographic information, we developed models that describe in spatial context the probability of the presence of each pine species evaluated.The next section describes the results for each species.We set a threshold of 75% (p > 0.75) as a way of representing the likelihood of presence for each pine.

Pinus culminicola
According to the results, the dwarf pine (P.culminicola) (Figure 3) is found in the east mountains of Saltillo city, and southeast of Monterrey City, as well as small areas in the southeast of Nuevo Leon state.The modeling suggests the mean precipitation of the coldest quarters (Bio19) as the most important environmental variable (Table 4).The area estimated as suitable for P. culminicola is about 421 km 2 .

Pinus Jaliscana
Distribution of P. jaliscana was mainly driven by precipitation variables (Bio18, Bio12, Bio16) as can be seen in Table 5.The results showed that the area suitable for the development of the species was 5502 km 2 and was well distributed in Jalisco state (Figure 4).

Pinus jaliscana
Distribution of P. jaliscana was mainly driven by precipitation variables (Bio18, Bio12, Bio16) as can be seen in Table 5.The results showed that the area suitable for the development of the species was 5502 km 2 and was well distributed in Jalisco state (Figure 4).

Pinus maximartinenzii
In the case of P. maximartinenzii, the environmental variables that were significant to the model were the mean precipitation in the warmer (Bio18) and colder (Bio19) quarters of the year (Table 6).The distribution is the western part of the state of Durango and the southwest of Zacatecas, which are dominated by pine forests (Figure 5).The potential suitable area is estimated at 18,355 km 2 .

Pinus maximartinenzii
In the case of P. maximartinenzii, the environmental variables that were significant to the model were the mean precipitation in the warmer (Bio18) and colder (Bio19) quarters of the year (Table 6).The distribution is the western part of the state of Durango and the southwest of Zacatecas, which are dominated by pine forests (Figure 5).The potential suitable area is estimated at 18,355 km 2 .

Pinus nelsonii
This species (P.nelsonii) is well adapted to the topographic and weather conditions of the northeastern Mexican mountains, between the states of Nuevo Leon, San Luis Potosi, and Tamaulipas,  This species (P.nelsonii) is well adapted to the topographic and weather conditions of the northeastern Mexican mountains, between the states of Nuevo Leon, San Luis Potosi, and Tamaulipas, where the main distribution occurs (Figure 6).The modeling procedure selected variables for both precipitation (Bio16, Bio17, Bio19) and topographic (aspect and slope) variables as they had significant contributions to modeling spatial distribution (Table 7).The estimated area was about 814 km 2 .

Pinus pinceana
This species has a wide distribution in the Sierra Madre Oriental of Mexico (30,745 km 2 ), mainly in the states of Coahuila, Nuevo León, San Luis Potosí, Zacatecas, and part of Durango (Figure 7), as well as scattered areas of Hidalgo and Querétaro.Precipitation (Bio17), topography (aspect), and temperature (Bio6, Bio5, and Bio11) variables significantly determined the distribution of this species (Table 7).

Pinus pinceana
This species has a wide distribution in the Sierra Madre Oriental of Mexico (30,745 km 2 ), mainly in the states of Coahuila, Nuevo León, San Luis Potosí, Zacatecas, and part of Durango (Figure 7), as well as scattered areas of Hidalgo and Querétaro.Precipitation (Bio17), topography (aspect), and temperature (Bio6, Bio5, and Bio11) variables significantly determined the distribution of this species (Table 8).
This species has a wide distribution in the Sierra Madre Oriental of Mexico (30,745 km 2 ), mainly in the states of Coahuila, Nuevo León, San Luis Potosí, Zacatecas, and part of Durango (Figure 7), as well as scattered areas of Hidalgo and Querétaro.Precipitation (Bio17), topography (aspect), and temperature (Bio6, Bio5, and Bio11) variables significantly determined the distribution of this species (Table 7).

Pinus rzedowskii
This species had a restricted distribution on the Pacific side of the Sierra Madre Occidental (2936 km 2 ), in the state of Michoacán, and some areas to the north of Cerro Tancítaro (Figure 8).Topography (aspect), precipitation (Bio19), and temperature variables (Bio2) determined the distribution of this species (Table 9).

Pinus rzedowskii
This species had a restricted distribution on the Pacific side of the Sierra Madre Occidental (2,936 km 2 ), in the state of Michoacán, and some areas to the north of Cerro Tancítaro (Figure 8).Topography (aspect), precipitation (Bio19), and temperature variables (Bio2) determined the distribution of this species (Table 8).

Validation of the Models
The beta coefficients, standard errors, z values, and p values are presented in Table 9.Additionally, the receiver operating characteristic curve (ROC) was used as a measure of fit (Figure 9).

Validation of the Models
The beta coefficients, standard errors, z values, and p values are presented in Table 10.Additionally, the receiver operating characteristic curve (ROC) was used as a measure of fit (Figure 9).

Hotspots or Areas of Possible Concern for Conservation
The areas of biggest concern due to a change of vegetation were coniferous forest distributed in the state of Durango, Jalisco, and Michoacán within the mountain range known as "The Sierra Madre Occidental" (Figure 10).This land cover change would affect potential areas of distribution of P. pincena, P. maximartinenzi, P. jaliscana, and P. rzedowskii.Scrubland from the northeast of Mexico is also suffering land cover changes and possible impacts on P. culminicola and P. nelsonii populations and suitable areas of distribution.

Hotspots or Areas of Possible Concern for Conservation
The areas of biggest concern due to a change of vegetation were coniferous forest distributed in the state of Durango, Jalisco, and Michoacán within the mountain range known as "The Sierra Madre Occidental" (Figure 10).This land cover change would affect potential areas of distribution of P. pincena, P. maximartinenzi, P. jaliscana, and P. rzedowskii.Scrubland from the northeast of Mexico is also suffering land cover changes and possible impacts on P. culminicola and P. nelsonii populations and suitable areas of distribution.During the period evaluated, coniferous forests and grasslands had the most significant loss regarding the area changing to a different vegetation type (Table 11, Figure 11).As was expected, agriculture, secondary vegetation, and urban areas increased due to the loss of forest area, grassland, and scrubland.Still, scrubland covers 42% of the total area evaluated (Table 11, Figure 11).During the period evaluated, coniferous forests and grasslands had the most significant loss regarding the area changing to a different vegetation type (Table 11, Figure 11).As was expected, agriculture, secondary vegetation, and urban areas increased due to the loss of forest area, grassland, and scrubland.Still, scrubland covers 42% of the total area evaluated (Table 11, Figure 11).

Discussion
The modeling of the ecological niche of the species is a useful tool for biological conservation [30] because it allows us to determine the situation and define conservation strategies, especially for threatened species.To advance the conservation of six endangered and endemic pine species from Mexico, in this paper we defined their natural distribution using historical observations and modeled the suitable areas of distribution based on their environmental requirements.We identified locations where the species are distributed naturally across the country but also the suitable areas that can serve as habitat for those species.Both distributions (natural and potential) let us detect differences, with strong implications for forest conservation.The estimated potential distribution areas varied according to the records of the species; P. rzedowskii and P. pinceana had the smallest and largest distribution areas, respectively.In a similar study [31], the smaller and largest distribution areas were for P. culminicola and P. pinceana, respectively; the differences may be due to the number of records, the method used, the number of interactions, and the different response variables.

Implications and Conservation Strategies
The modeled distribution of P. culminicola coincided with previous records at Cerro Potosí, Sierra La Viga, and Sierra La Marta in the states of Coahuila and Nuevo León [32,33].The estimated potential distribution areas reaffirm the restricted distribution of this species, which is strongly determined by precipitation in the coldest quarter (Bio19); however, it has been found that the minimum temperature of the coldest month, precipitation in the wettest quarter, and total precipitation also determine the distribution of this species [33].Even with the protection of Mexican and international laws, the populations of this species are affected by land use change, road construction, cattle grazing, and forest wildfires [10], causing the population of this species to undergo drastic reductions.For example, in the Cerro Potosí from 1960 to 1970 the population decreased 66% [34].According to some reports [33,35,36], P. culminicola covers between 10 and 20 km 2 ; however, the potential area determined in this study is higher (421 km 2 ), which suggests that there are large areas with adequate conditions for ex situ conservation of this species.Therefore, the specific in situ conservation activities for P. culminicola should be aimed at reducing the effect of fire and livestock, since no seedling of this species survived the grazing and trampling of cattle after four years [10].
Pinus jaliscana is distributed in nine populations in the west of the state of Jalisco [37,38].Agricultural activities, livestock, and wildfires put natural regeneration at risk in some populations [38].It is estimated that this species covers more than 4185 ha [37], which is lower than what was estimated in this study (5502 km 2 ).The above suggests a vast area with potential for ex situ conservation of this species, including areas of the north of the state of Jalisco where there are no reports of P. jaliscana.On the other hand, the seed efficiency rate is 4% [38], which indicates reproductive problems associated with self-fertilization due to the isolation and small size of the populations.Therefore, in situ conservation activities for this species should be oriented towards increasing the genetic variability of the populations to avoid inbreeding problems.
Pinus maximartinezii has only two natural populations, one in the southern state of Zacatecas and another in southern Durango [39].The population located in the state of Durango has an area of 110 ha, while the population of southern Zacatecas can have up to 925 ha [39]; however, the potential distribution area determined in this study is considerably larger (18,355 km 2 ) and extends north of its distribution, on the Sierra Madre Oriental, in the state of Durango and Zacatecas.Our results suggest the possible existence of other populations that have not yet been located, but also large areas with conditions suitable to establish plantations for ex situ conservation strategies.The seeds of this species are larger and more nutritious than other pine species [40], so they are used as food; however, even with the use of seeds, the population growth rate of the southern Zacatecas population is 1117, which suggests that the use of seeds does not put this population at risk [41].Still, the population of southern Zacatecas has a low percentage of polymorphic loci (30.3%) and only two alleles per polymorphic locus, suggesting a bottleneck effect less than 100 years ago [42].The in situ conservation activities for this species should be oriented towards favoring the survival of young and adult individuals during the first reproductive events, protecting the habitat of the species [41], and increasing genetic diversity through the exchange of genetic material of both populations.
Pinus nelsonii is endemic and restricted in distribution to the Sierra Madre Oriental, in the states of Tamaulipas, San Luis Potosí, Nuevo Leon, and Coahuila [43,44].Due to the decline in its populations, it is categorized as an endangered species by NOM-059-SEMARNAT-2010 (SEMARNAT 2010) and IUCN [45].Although, among the pines of the Sierra Madre Oriental, P. nelsonii is the second most important species [12], there is little information on the number and size of its populations.It is estimated that it covers only 84 km 2 [45], smaller than the potential distribution area (814 km 2 ) defined in this study, but partially coinciding with the current distribution; if we know the zones we can carry out ex situ conservation programs such as protection plantations.The natural populations of this species are severely affected by pests and forest fires.For example, in a population of Tamaulipas, 98% of the plants were affected by the pinyon pitch nodule moth (Retinia arizonensis) [44], and 30% of the area was affected by droughts, which reduced their rate of population growth (from 1.083 to 0.990) in a single year [12].On the other hand, the variation and genetic differentiation of several populations of this species are low to moderate [43].Therefore, the in situ conservation activities of P. nelsonii should be oriented to avoid forest fires to favor natural regeneration, monitoring and control of sanitation, as well as increasing the genetic diversity of the populations.
Among the species studied in this work, P. rzedowskii has a more restricted distribution; it is only distributed in the state of Michoacán [5].Only one population with 12 localities is recognized, which in total hold approximately 6500 individuals [46].Although the genetic structure is marked, and the genetic diversity is moderate to high [46,47], because the populations are tiny, the risk of extinction due to significant forest fires is high [48].Therefore, in situ conservation activities for this species should be oriented towards forest fire prevention measures, as well as avoiding activities related to agriculture and grazing.Of the six species studied, P. rzedowskii has a big potential area of modeled distribution (2936 km 2 ), which can be explained by its restricted distribution; this area is enough to establish plantations for ex situ conservation.
Among the species studied, P. pinceana has the widest distribution; its distribution extends 750 km along the Sierra Madre Oriental [49] in the states of Hidalgo, Querétaro, San Luis Potosí, Zacatecas, and Coahuila [50,51].Despite its wide distribution, its populations are fragmented; there is low natural regeneration and potential threats due to the collection of firewood and seeds, but there are no signs of significant recent or current decline, which is why IUCN considers it to be of less concern [49].Also, the differentiation and genetic diversity of the populations are high [11].Within the natural distribution of P. pinceana, according to the composition and coefficient of floristic similarity, they differ among three sets of communities [50], which coincides partially with genetic isolation between populations from the north, center, and south [11].The exact area covered by this species is unknown, but it is estimated at 2000 km 2 [49], which is lower than the potential area modeled in the present study (30,745 km 2 ).On the other hand, the northeast of Mexico presents conditions of high aptitude for this species, with 10,268 km 2 suitable for the distribution of the species [52].These results suggest a large potential area for the establishment of conservation plantations of this species as an ex situ conservation strategy.
For target species, there are seed dispersal limitations due to low tree density, wildfires, grazing, small animals, and extreme climatic conditions: e.g., P. culminicola seedlings were not detected in Cerro El Potosí for four years due to them being consumed by rodents [10], and the survival of the species was strongly affected by low temperatures, strong winds, scarce rainfall, low relative humidity, little soil depth, and loss of edaphic material [53].Seedlings of P. maximartinezi have growth problems in natural stands because of soil degradation due to grazing, fires, agricultural activities, and desiccation [41].P. nelsonii populations have fewer than 10,000 mature trees and face two main problems: increased cattle and incidence of destructive fires [54], both of which reduce seed germination and seedling growth.
In Mexico, the primary strategy of in situ conservation is the establishment of Protected Natural Areas (ANPs).In this sense, the studied species can be conserved in the existing ANPs, and, in the case of P. pinceana, P. culminicola, and P. nelsonii, also in proposed ANPs [55].In the present study, although the percentage of the potential distribution area corresponding to ANPs was not determined, between 1.4% (P.rzedowskii) and 17.0% (P.maximartinezii) of this surface could encompass ANPs [31].However, the network of ANPs in Mexico does not adequately protect the Pinus species [31], so the in situ conservation activities discussed for each species are necessary, as is the primary strategy of in situ conservation.Special attention should be given to strategies related to the use of seeds for human consumption; the establishment of Management Units for the Conservation of Wildlife (UMAs) could be proposed as an option because the owners can make sustainable use of the habitat, and at the same time work towards the conservation of the species [56].Although in situ conservation is essential to renew genetic diversity and face future environmental changes, ex situ conservation is operationally convenient for short-term results [57].It is suggested that the main ex situ conservation strategy should be the establishment of conservation plantations in the areas determined in this study.
The establishment of conservation plantations exclusively in the areas determined for each species should enhance the success and adaptation of the plantations.Also, ex situ conservation strategies should include the collection, storage, and conservation of germplasm in botanical gardens, germplasm banks, and plant tissue culture laboratories in Mexico [58].Due to the genetic problems that several of these species present, the germplasm collection for ex situ conservation must be carried out in all the localities or populations of each species.Although Mexico has made significant progress in regulating the use of plant species [56], legislation must be revised and adapted to favor in situ and ex situ conservation activities with the species studied in this work, as well as others that are at risk.

Final Considerations
In this study, the primary potential distribution of the target species was located within the temperate coniferous forest from the central part of Mexico (part of the "Sierra Madre Occidental" and "Sierra Madre Oriental" systems), but also within other vegetation communities such as mixed forests and scrubland.Our results are in agreement with some previous reports on the distribution of endangered Mexican pines [5,31,33,39,44].It is important to highlight that the spatial distribution and number of observations collected and retrieved from the diverse databases were an essential factor in the performance of all models and the ability to predict the potential spatial distribution of Mexican pines accurately.
Although it was not part of this study, it would be worth performing an analysis by region including other crucial variables such as distance to roads or population density, as they are considered factors in the vulnerability of Mexican forests [59].It is also important to highlight that land cover change is affecting potential suitable areas of distribution of the species.Forest net loss was about 25%.It is imperative to consider land use change and land cover dynamics to promote conservation strategies for each species.The Mexican territory suffers fragmentation, disturbances, and some degradation processes [59][60][61], causing changes in the structure, function, composition, productivity, and extent of forests that may contribute to the diminishing of forest species [62,63].
In this paper we used spatial modeling and niche ecological models combined with other analytical tools (e.g., GIS) to evaluate non-inventoried sites as well as model the distribution of species [14].The selection of methods used accomplished a double function: first, to provide knowledge about the potential distribution of the species and use said predictions in the choice of sites of particular interest as biological conservation zones [19,20].Although species distribution models (e.g., MaxEnt algorithm) can fit with limited presence-only data [64], logistic regression was considered a better option as it allows better control of modeling compared to MaxEnt [16,65].Still some issues need to be considered before selecting the appropriate modeling method: the clarification of the niche concept, sampled data, parametrization strategies, model selection and predictor contribution, and evaluation strategies [66].
Our approach, as we stated before, took advantage of freely available data and, although there was quality control during the modeling, we may have carried over some errors from the data collection, identification of the species, and registration of the geographic data.Still, our results showed good agreement with the performance of all models according to the validation process.As far as we know, there is no maximum number of observations suitable for spatial modeling; however, species distributions models (SDMs) provide the most useful performance using a minimum of 14 to 25 observations for narrow-ranged and widespread species, respectively [67].In this study, we had a variety in the number of observations from the field, and assumed that the spatial distribution and number for each pine species were important to the performance of each of the models developed in this study.Our results corroborated the importance and usefulness of diagnosis data for occurrence prediction and their contribution to a better understanding of a species' reaction to climatic variables [68].
Finally, we would like to highlight, as future research endeavors, the inclusion of the effects of disturbances, landscape fragmentation and land use dynamics, and climate change on species distribution modeling to provide more comprehensive and accurate results for the implementation of conservation strategies.

Conclusions
This study presents valuable information about the distribution of six Mexican pine species catalogued as endemic and endangered.The approach followed here takes advantage of the freely available data and resources that academic institutions can provide in monitoring critical species such as P. culminicola, P. jaliscana, P. maximartinenzi, P. nelsonii, P. pinceana, and P. rzedoswki.We were able to model suitable potential areas of distribution of each species evaluated.The results of this study are intended to be the basis of in situ and ex situ conservation strategies.

Figure 1 .
Figure 1.Vegetation types from the study area.

Figure 1 .
Figure 1.Vegetation types from the study area.

Figure 2 .
Figure 2. The natural distribution of the six Mexican pines evaluated in this study.

Figure 2 .
Figure 2. The natural distribution of the six Mexican pines evaluated in this study.

Forests 22 Figure 10 .
Figure 10.Areas of land cover change in the period 1991-2016.Colors indicate changes from one vegetation class to another during the period evaluated.

Figure 10 .
Figure 10.Areas of land cover change in the period 1991-2016.Colors indicate changes from one vegetation class to another during the period evaluated.

Table 1 .
Characteristics of the pine species selected in this study.

Table 1 .
Characteristics of the pine species selected in this study.

Table 2 .
Herbaria consulted in this study.

Table 9 .
Indices used for validation of the models.

Table 10 .
Indices used for validation of the models.

Table 11 .
Land cover changes in the study area.