Research Article
Print
Research Article
Factors related to species richness, endemism, and conservation status of the herpetofauna (Amphibia and Reptilia) of Mexican states
expand article infoGeoffrey R. Smith, Julio A. Lemos-Espinal§
‡ Denison University, Granville, United States of America
§ FES Iztacala UNAM, Tlalnepantla, Mexico
Open Access

Abstract

Mexico is a megadiverse country with high endemicity in its herpetofauna. We examine how species richness, proportion of state and country endemic species, and proportion of species in a category of conservation concern using listings in the International Union for Conservation of Nature (IUCN) Red List and the Secretaría del Medio Ambiente y Recursos Naturales (SEMARNAT) in 27 of 32 Mexican states are related to environmental and human demographic and socioeconomic variables. Amphibian and reptile species richness were positively related to latitude range and number of physiographic regions and negatively related to latitude. The proportion of state endemic amphibian species in a state was negatively related to latitude whereas no variables influenced the proportion in reptiles. The proportion of country endemics in a state was positively related to human population density and the number of physiographic regions and negatively related to per capita gross domestic product (GDP) and latitude range for amphibians; it was positively related to human population density and elevation range and negatively related to latitude range for reptiles. The proportion of amphibian species in an IUCN category of concern in a state was positively related to human population density and negatively related to latitude; for reptiles, it was negatively related to human population density. The proportion of SEMARNAT-listed species in a state was positively related to human population density for both amphibians and reptiles and negatively related to latitude range for amphibians. Our analyses found that larger macroecological patterns (e.g., latitudinal species gradient, heterogeneity-richness relationships) and human population density play important roles in determining the richness and conservation status of Mexican amphibians and reptiles.

Keywords

amphibians, environmental factors, human demographic factors, IUCN status, reptiles, SEMARNAT listing, socioeconomic factors

Introduction

Much of Mexico consists of a transition zone (i.e., the Mexican Transition Zone) between the Nearctic and Neotropical zones (Villaseñor et al. 2020). As a consequence, Mexico has one of the highest levels of biodiversity of any country in the world (Morrone 2019), including amphibian and reptile richness (Wilson and Johnson 2010; Chen and Peng 2017). Mexico also has a high level of endemicity in its herpetofauna (Wilson and Johnson 2010; see also Murali et al. 2021).

Unfortunately, Mexico is not immune to global environmental decline. Indeed, Mexico is an area with high extinction debt and risk for amphibians and reptiles (Chen and Peng 2017). In addition, Mexico is one of the countries where the decline in the conservation status of amphibians is greatest (Rodrigues et al. 2014) and has a high proportion (ca 80%) of species of amphibians showing population declines (Becker and Loyola 2008). The biodiversity of Mexico is subject to a variety of pressures, such as climate change, land use change (including agriculture and livestock, mining, deforestation, and urbanization), invasive species, disease, exploitation, and pollution (Rodrigues et al. 2014; Lazcano et al. 2019; Ramírez-Bautista et al. 2020; González-Sánchez et al. 2021; Masés-García et al. 2021), which are the consequence of human demographic changes (e.g., population growth), socioeconomics, and governmental actions (Challenger et al. 2009; Sarukhán et al. 2015).

Human pressure on biodiversity has increased in the Neotropics (Geldmann et al. 2014). More specifically, mean annual temperatures in Mexico have increased around 0.2 °C from 1970–2000, with greater increases in northern Mexico than in southern Mexico (see also Pavia et al. 2009; Cuervo-Robayo et al. 2020). This climate change has driven changes in the vegetation and distribution of habitats, especially in the mountains of Mexico (Téllez-Valdés et al. 2006; Gómez-Mendoza and Arriaga 2007; Jiménez-García et al. 2021), and continued climate change is likely to result in the loss of suitable habitat in the future (Chacón-Prieto et al. 2021). Parts of Mexico are also undergoing rapid land use change through burning, human settlement, and conversion to agriculture, with extensive loss of forest, including in protected areas (Lorenzo et al. 2017; Hu et al. 2021). Habitat loss has impacted several species of terrestrial vertebrates in Mexico, especially endemic species in the Transvolcanic Mexican Belt, Mexican High Plateau, and the Humid Coastal Plains and Hills of the Gulf of Mexico, and the effects appear to be cumulative (i.e., not just recent habitat loss) (Mayani-Parás et al. 2021). The loss of habitat, and in particular forest habitats, has negatively affected the amphibians and reptiles of Mexico (Lara-Tufiño et al. 2019; Mayani-Parás et al. 2019), and a high proportion of endangered amphibians in Mexico are found in areas that have experienced transformation to agriculture or urbanization (Londoño-Murcia and Sánchez-Cordero 2011) and these trends are likely to continue or increase in the future (Mendoza-Ponce et al. 2020).

Here we examine the distribution of amphibian and reptile species richness among 27 of 32 Mexican states. In particular, we examine relationships between species richness, proportion of state and country endemic species (i.e., the proportion of the species in a state that are state or country endemics), and proportion of a state’s amphibian or reptile species in a category of conservation concern using the Interational Union for Conservation of Nature (IUCN) Red List (i.e., vulnerable, threatened, endangered, critically endangered, near extinction), and the proportion of a state’s amphibian or reptile species listed in Secretaría del Medio Ambiente y Recursos Naturales (SEMARNAT) (2019) with environmental variables (state area, proportion of land protected, latitude, latitude range, elevation range, and number of physiographic regions) and human demographic and socioeconomic variables [human population, human population density, and per capita gross domestic product (GDP)].

Methods

We collected species lists for amphibians and reptiles of Mexican states from the available literature and updated these species lists using additional literature through November 2020 (see Suppl. material 1: Table S1 for sources used for base species lists and updates). We generally followed Frost (2020) and AmphibiaWeb (2020) for amphibian taxonomy and Uetz and Hošek (2019) for reptile species. We were able to compile updated species lists for 27 of the 32 Mexican federal entities (i.e., states), with the five remaining states lacking published updated species checklists (Suppl. material 2: Table S2). We include Mexico City (formerly known as Mexico, Distrito Federal) that comprises the urban area of Mexico City proper to the south and mountains and valleys with fragmented forests and grasslands to the north. For each species in our list, we obtained their global conservation status from the IUCN Red List version 2021-3 (https//:www.iucnredlist.org) and their Mexican conservation status from SEMARNAT (2019) (Suppl. material 2: Table S2). From these lists we gathered the following information for each state: species richness, proportion of state and country endemics, proportion of species in an IUCN category of concern (critically endangered, endangered, threatened, near threatened, and vulnerable), and the proportion of SEMARNAT-listed species under the categories of threatened and endangered for amphibians and reptiles separately (Table 1). For each state we collected data on human demographic and socioeconomic variables (human population, human population density, and per capita GDP) and geographic and climatic variables [state area, proportion of land protected, latitude (midpoint of state), latitude range (difference between minimum and maximum latitude), elevation range (difference between minimum and maximum elevations), and number of physiographic regions] (Table 2).

Table 1.

Amphibians and reptile species richness, proportion of species in a state that are state and country endemics, proportion of species that are in an IUCN category of concern, and the proportion of species that are SEMARNAT listed for Mexican states.

State Amphibians Reptiles
Species Richness Prop. State Endemic Prop. Country Endemic Prop. IUCN Prop. SEMARNAT Species Richness Prop. State Endemic Prop. Country Endemic Prop. IUCN Prop. SEMARNAT
Aguascalientes 19 0 0.68 0.16 0.16 65 0 0.54 0.05 0.15
Baja California 17 0 0 0.29 0.06 103 0.23 0.32 0.21 0.23
Baja California Sur 3 0 0 0 0 84 0.44 0.57 0.13 0.32
Campeche 24 0 0.04 0 0 103 0 0.12 0.13 0.16
Chiapas 108 0.10 0.21 0.53 0.04 219 0.06 0.18 0.11 0.16
Chihuahua 37 0.03 0.35 0.11 0.03 140 0.01 0.33 0.06 0.15
Coahuila 24 0 0.25 0.17 0.04 117 0.06 0.30 0.15 0.22
Colima 39 0.03 0.68 0.13 0.08 117 0.04 0.66 0.11 0.18
Durango 36 0 0.6 0.11 0.03 119 0.02 0.45 0.05 0.18
Guerrero 78 0.32 0.76 0.45 0.11 181 0.15 0.68 0.08 0.16
Hidalgo 52 0.08 0.65 0.46 0.15 126 0.01 0.49 0.10 0.18
Jalisco 55 0.06 0.72 0.22 0.07 173 0.01 0.66 0.10 0.17
Mexico 49 0.08 0.77 0.42 0.21 99 0.01 0.71 0.08 0.17
Mexico City 18 0.11 0.83 0.5 0.39 45 0 0.76 0.04 0.22
Michoacán 58 0.09 0.74 0.26 0.09 161 0 0.7 0.09 0.18
Morelos 38 0 0.70 0.27 0.14 97 0 0.68 0.04 0.17
Nayarit 37 0 0.6 0.08 0.03 117 0.01 0.60 0.07 0.14
Nuevo León 25 0 0.32 0.2 0.04 119 0.02 0.34 0.09 0.20
Oaxaca 152 0.39 0.69 0.60 0.10 299 0.11 0.55 0.11 0.18
Puebla 92 0.03 0.7 0.49 0.16 176 0.02 0.61 0.07 0.14
Querétaro 34 0.60 0.56 0.32 0.12 104 0.02 0.49 0.07 0.18
Quintano Roo 23 0 0.09 0.04 0 108 0.01 0.12 0.12 0.16
San Luis Potosí 42 0 0.41 0.34 0.12 138 0 0.41 0.09 0.17
Sinaloa 39 0 0.55 0.13 0 119 0.01 0.51 0.10 0.14
Sonora 36 0 0.39 0.11 0.03 159 0.10 0.37 0.13 0.18
Tamaulipas 44 0.11 0.41 0.27 0.14 137 0.04 0.34 0.11 0.21
Yucatán 17 0 0.12 0.06 0 85 0 0.12 0.12 0.18
Table 2.

Human demographic and socioeconomic variables and environmental variables for Mexican states.

State State Area (km2)1 Human Population (2018)1 Human Population Density (N/km2)1 Per capita GDP (US$)2 Proportion Protected Territory3 GPS Coordinates1 (°) Elevation Range (m)1 Number of Physiographic Regions4
Aguascalientes 5618 1.337,792 238.1 9975 26.0 22.1243, 1.0042 1666 3
Baja California 71.450 3.633,772 50.9 9449 19.0 30.3593, 4.7186 3100 2
Baja California Sur 73.909 832.827 11.3 11.060 42.0 25.4360, 5.1280 2080 1
Campeche 57.507 948.459 16.5 51.460 39.6 19.6167, 0.7667 390 2
Chiapas 73.311 5.445,233 74.3 3592 18.0 17.2588, 5.4530 4080 3
Chihuahua 247.460 3.816,865 15.4 8833 8.1 28.6843, 6.4175 3050 2
Coahuila 151.595 3.063,662 20.2 12.838 19.0 27.2114, 5.3372 3380 3
Colima 5627 759.686 135 9177 6.6 19.0983, 0.8283 3820 2
Durango 123.317 1.815,966 14.5 7888 22.1 24.5950, 4.5000 3240 4
Guerrero 63.596 3.625,040 5.7 4586 0.15 17.6018, 2.5719 3550 2
Hidalgo 20.813 2.980,532 143.2 6508 6.9 20.4982, 1.8008 3251 3
Jalisco 78.588 8.197,483 104.3 9239 11.2 20.8380, 3.8244 4339 4
Mexico 22.500 17.604,619 782.4 6199 43.8 19.3264, 1.9189 5268 2
Mexico City 1495 8.788,141 5878.4 21.079 14.1 19.3206, 0.5444 1702 1
Michoacán 58.599 4.687,211 80 5522 5.9 19.1547, 2.4794 4100 2
Morelos 4879 1.987,596 407.4 6961 26.8 18.7319, 0.7994 4580 2
Nayarit 27.857 1.290,519 46.3 6220 30.8 21.8439, 2.4811 2760 4
Nuevo León 64.156 5.300,619 82.6 16.228 8.9 25.4810, 4.6364 3660 3
Oaxaca 93.757 4.084,674 43.6 4446 7.1 17.1635, 3.0125 3720 5
Puebla 34.306 6.371,381 185.7 5890 19.5 19.3500, 2.9667 5530 4
Querétaro 11.699 2.091,823 178.8 12.502 33.6 20.8425, 1.655 2600 3
Quintano Roo 50.212 1.709,479 34 11.381 32.6 19.7000, 3.7667 230 1
San Luis Potosí 61.137 2.824,976 46.2 8118 6.6 22.8258, 3.3311 3160 3
Sinaloa 58.328 3.059,322 52.5 8108 7.6 24.7547, 4.5750 2520 2
Sonora 179.355 3.050,473 17 11.543 10.3 29.3954, 6.1969 2620 4
Tamaulipas 80.249 3.661,162 45.5 9347 13.7 24.9430, 5.4722 3280 3
Yucatán 39.524 2.199,618 55.7 8.011 25.6 20.5667, 2.0667 210 1

We used generalized linear models (Normal distribution, Identity link) for amphibians and reptiles separately to examine the relationships of the human demographic and socioeconomic variables and the geographic and climatic variables and species richness, proportion of country endemics, proportion of state endemics, proportion of species in an IUCN category of concern, and proportion SEMARNAT-listed species. We used JMP Pro 15.1 (SAS Institute, Cary, NC). for statistical analyses.

Results

Amphibian species richness was positively related to latitude range and number of physiographic regions and negatively related to latitude (Table 3; Fig. 1A-C). The proportion of state endemics was negatively related to latitude (Table 3; Fig. 1D). The proportion of country endemics was positively related to human population density and the number of physiographic region and negatively related to per capita GDP and latitude range (Table 3; Fig. 1E-H). The proportion of amphibian species in an IUCN category of concern was positively related to human population density and negatively related to latitude (Table 3; Fig. 2A, B). The proportion of SEMARNAT-listed amphibian species was positively related to human population density and negatively related to latitude range (Table 3; Fig. 2C, D).

Table 3.

Results of generalized linear models examining the relationship between human demographic and socioeconomic variables and environmental variables and species richness, proportion of species in a state that are state and country endemics, proportion of species that are in an IUCN category of concern, and the proportion of species that are SEMARNAT listed for amphibians in 27 Mexican states. Values on first line are coefficients, values on second line are P-values. Bolded entries are significant at α = 0.05.

Variable Intercept State area Human population Human population density per capita GDP Proportion Protected Latitude Latitude Range Elevation range Physiographic regions Overall P
Species Richness 153.03 < 0.0001 0.00011 0.16 < 0.00001 0.63 -0.0010 0.725 -0.0002 0.18 -44.61 0.074 -7.37 < 0.0001 5.91 0.043 0.0025 0.39 10.62 0.0006 < 0.0001
Prop State endemic 0.447 0.0003 < 0.000001 0.17 <0.0000001 0.39 <0.000001 0.65 -<0.000001 0.49 -0.247 0.057 -0.018 0.0009 0.004 0.81 -<0.00001 0.52 0.013 0.36 0.0108
Prop Country endemic 0.568 0.0046 0.000001 0.10 <<0.000001 0.36 0.000062 0.024 -0.0000010 0.004 -0.41 0.066 0.00055 0.95 -0.116 < 0.0001 0.000047 0.078 0.069 0.008 <0.0001
Prop IUCN 0.469 0.0065 -<0.000001 0.35 <<0.000001 0.26 0.000050 0.0336 -0.0000016 0.55 -0.33 0.086 -0.018 0.017 0.022 0.30 0.000029 0.20 0.035 0.10 0.0002
Prop SEMARNAT 0.0045 0.93 -0.0000001 0.56 <<0.0000001 0.14 0.000050 < 0.0001 -0.0000012 0.19 -0.042 0.19 0.0034 0.18 -0.018 0.016 0.000012 0.12 0.0096 0.19 <0.0001
Figure 1. 

The relationships between amphibian species richness and a state’s latitude A latitude range B and number of physiographic regions C between the proportion of a state’s amphibian species that are state endemics and a state’s latitude D and between the proportion of a state’s amphibian species that are country endemics and the state’s human population density E per capita GDP F latitude range G and the number of physiographic regions for Mexico H.

Figure 2. 

The relationships between the proportion of a state’s amphibian species that are in an IUCN category of concern and a state’s human population density A and latitude B and between the proportion of a state’s amphibian species that are listed in SEMARNAT (2019) and a state’s human population density C and latitude range for Mexico D.

Reptile species richness was positively related to latitude range and the number of physiographic regions and negatively related to latitude and the proportion of land protected (Table 4; Fig. 3A-D). The proportion of state endemics of reptiles was not related to any variables (Table 4). The proportion of country endemic reptile species was positively related to human population density and elevation range and negatively related to latitude range (Table 4; Fig. 3E-G). The proportion of a state’s reptile species in an IUCN category of concern was negatively related to human population density (Table 4; Fig. 3H). The proportion of a state’s reptile species that are SEMARNAT listed was positively related to human population density (Table 4; Fig. 3I).

Table 4.

Results of generalized linear models examining the relationship between human demographic and socioeconomic variables and environmental variables and species richness, proportion of species in a state that are state and country endemics, proportion of species that are in an IUCN category of concern, and the proportion of species that are SEMARNAT listed for reptiles in 27 Mexican states. Values on first line are coefficients, values on second line are P-values. Bolded entries are significant at α = 0.05.

Variable Intercept State area Human population Human population density per capita GDP Proportion protected Latitude Latitude Range Elevation range Physiographic regions Overall P
Species Richness 279.65 < 0.0001 0.00023 0.06 0.0000005 0.73 -0.0084 0.061 0.00026 0.62 -109.7 0.0046 -10.74 < 0.0001 13.35 0.0031 0.0013 0.76 16.57 0.0003 <0.0001
Prop State endemic -0.070 0.58 -0.0000002 0.61 <0.0000001 0.074 0.000015 0.40 -<0.00001 0.79 0.272 0.072 0.0009 0.87 0.031 0.08 0.000025 0.17 -0.024 0.15 0.22
Prop Country endemic 0.245 0.18 0.00000035 0.61 <0.0000001 0.39 0.000068 0.011 -0.0000065 0.039 0.082 0.69 0.0094 0.26 -0.067 0.009 0.000091 0.0011 -0.0016 0.95 0.0002
Prop IUCN 0.051 0.28 -0.0000003 0.095 <0.0000001 0.34 -0.000015 0.035 0.0000014 0.10 -0.064 0.24 0.0024 0.28 0.0122 0.063 -0.000005 0.42 -0.0060 0.34 0.076
Prop SEMARNAT 0.079 0.087 -0.0000002 0.23 -<<0.00001 0.084 0.000015 0.028 0.00000025 0.74 0.088 0.10 0.0031 0.14 0.010 0.09 0.000011 0.10 -0.010 0.083 0.054
Figure 3. 

The relationships between reptile species richness and a state’s proportion of protected land A latitude B latitude range C and number of physiographic regions D between the proportion of a state’s reptile species that are country endemics and the state’s human population density E latitude range F and elevation range G between the proportion of a state’s reptile species that are in an IUCN category of concern and a state’s human population density H and between the proportion of a state’s reptile species that are listed in SEMARNAT (2019) and a state’s human population density for Mexico I.

Discussion

Our analyses found that species richness, endemism, and conservation status of amphibians and reptiles in Mexican states are related to both human demographic and socioeconomic variables and environmental variables. Below we discuss our observations on the factors related to species richness, endemism, and conservation status of amphibians and reptiles in Mexico.

Species richness and endemicity

For amphibians, species richness was positively related to latitude range and the number of physiographic regions in a state and negatively related to latitude (i.e., species richness decreased with latitude), whereas reptile species richness was positively related to latitudinal range and the number of physiographic regions and negatively related to latitude and the proportion of land protected. The proportion of a state’s amphibian species that are state endemics was negatively related to latitude whereas none of the state variables we examined were related to the proportion of a state’s reptile species that are state endemics. For the proportion of a state’s species of amphibians that are country endemics there was a positive relationship with human population density and the number of physiographic regions and a negative relationship with per capita GDP and latitude range. For reptiles, this proportion was positively related to human population density and elevation range and negatively related to latitude range.

The positive relationships between species richness and latitude range and the number of physiographic regions for both amphibians and reptiles likely result from the increased variety of habitats and climates in a state leading to a greater number of niches, which can result in increased number of species occurring in a state. Our results for the herpetofauna of Mexico are similar to studies on Sceloporus lizards (Rivera et al. 2021) and on mammals finding that species richness at a variety of geographic scales increases with habitat or environmental heterogeneity (Amori et al. 2013, 2019; Udy et al. 2021).

The decrease in species richness of amphibians and reptiles with increasing latitude is consistent with the latitudinal species gradient (Pianka 1966; Willig et al. 2003; Hillebrand 2004; Pontarp et al. 2019). Species richness of amphibians in Mesoamerica is highest to the south (Wilson and Johnson 2010). Similarly, mammalian species richness in Mexico decreased with increasing latitude (Ceballos et al. 1998). These results are also similar to those found for amphibians and reptiles in a variety of regions and countries, such as Europe (Assunção-Albuquerque et al. 2012), United States (Schall and Pianka 1978), Australia (Schall and Pianka 1978), and North America (Rivera et al. 2021; Whiting and Fox 2021). The latitudinal gradient of species richness in both amphibians and reptiles in Mexico might be related to latitudinal gradients in climate related factors. For example, the species richness of reptiles and amphibians is often related to solar radiation, temperature, precipitation, annual potential and actual evapotranspiration (e.g., Schall and Pianka 1978; Rodríguez et al. 2005; Powney et al. 2010; Kafesh et al. 2020), all of which vary with latitude.

In addition, for reptiles, species richness decreased with increasing human population density and the proportion of the state’s area protected. This relationship suggests there might be a negative impact of human population on species richness, perhaps due to the impact on amphibian and reptile populations. For example, in Europe, climate and human related factors explained 41–42% of variation in species richness of amphibians and reptiles (Assunção-Albuquerque et al. 2012). The species richness of reptiles in Mediterranean France is correlated with climate, elevation, and land use (Barnagaud et al. 2021). Reptile species richness is negatively affected by livestock production and urbanization (Cordier et al. 2021). Thus, the species richness of reptiles appears to be particularly susceptible to human pressures.

The patterns of state endemism that we observed show limited effects of the state variables we examined, with only a negative relationship between latitude and the proportion of a state’s amphibian species being state endemics. In part this general lack of relationships may reflect the artificial nature of state boundaries (see Caveats below) such that states often share such physiographic regions or habitats and so likely share species, even over a small area along borders. For country-level endemism, we found positive relationships with human population density and measures related to habitat heterogeneity (e.g., number of physiographic regions or elevation range) for both amphibians and reptiles. Such relationships likely reflect the high level of endemism found in the Transvolcanic Mexican Belt (Flores-Villela et al. 2010) which is found in an area of high human population density (e.g., Mexico City and its environs). The negative relationship between country-wide endemism and latitude range likely also reflects the smaller states and federal entities found in central and southern Mexico where high levels of endemism are found.

Conservation status

Both human demographic and environmental variables affected the proportion of both amphibians and reptiles of conservation concern in a state. The proportion of amphibian species in an IUCN category of concern decreased with latitude and increased with human population density. For SEMARNAT, the proportion of listed species increased with human population density and decreased with latitudinal range. For reptiles, the proportion of a state’s species in an IUCN category of concern decreased with human population density, but the proportion that is SEMARNAT listed increased with human population density. The importance of human population density in determining conservation status is consistent with the impacts that anthropogenic effects on the environment have on amphibians and reptiles. For example, a high proportion of endangered amphibians in Mexico being found in areas that have experienced transformation to agriculture or urbanization (Londoño-Murcia and Sánchez-Cordero 2011). In addition, amphibian species diversity in central Mexico is reduced with the loss of canopy (Lara-Tufiño et al. 2019). Mayani-Parás et al. (2019) demonstrated that many species of amphibians and reptiles in Mexico have suffered major reductions in their distributions, and that this is particularly the result of the combined effects of mining and habitat loss. In addition, reptiles, and, to a lesser extent, amphibians, are subject to illegal trade and collection in Mexico (Masés-García et al. 2021). These patterns in Mexico are also consistent with patterns of conservation status and human pressures at the global and more regional scales for amphibians and reptiles. Globally, amphibian species richness is susceptible to deforestation, timber harvesting and production, and urbanization (Cordier et al. 2021). Reptile species richness is negatively affected by livestock production and urbanization (Cordier et al. 2021). In addition, reptile abundance globally is negatively affected by anthropogenic habitat changes (Doherty et al. 2020). Declines in European amphibians and reptiles were associated with the number of alien species and loss of habitat due to climate change (Falaschi et al. 2019). In China, the distribution of biodiversity loss is driven primarily by climate and anthropogenic sources (Lu et al. 2020). Similarly, in Australia extinction risk in reptiles is related to anthropogenic pressures and proximity to human populations as well as climatic variables (Senior et al. 2021).

Caveats

We recognize that our analysis is a snapshot in time of both taxonomic knowledge and conservation status. Species lists are dynamic and changing as new species are described, new localities are found, and populations are extirpated. In addition, conservation status for species, whether at the global (i.e., IUCN) or national (i.e., SEMARNAT) scale, are frequently reassessed. Thus, we realize that our analysis represents our knowledge at the time we generated our species lists and conducted our analyses. However, obtaining a snapshot now will provide a baseline that can be monitored and evaluated as our understanding of taxonomy, species distributions, and conservation status change over time.

We also recognize that by using global and national conservation status we are not taking into account variation in conservation status of species populations in each state such that species may be doing well in some states but poorly in other states. This is a challenge, but unfortunately conservation status at the state level is known for even fewer species than for the national and global measures. In addition, many regulations are focused on, or use, national or global level assessments rather than state level assessments.

Finally, we further recognize that the use of states in our analyses does not necessarily reflect ecological or biogeographical reality (i.e., they are not natural regions). However, given the nature of governmental processes, conservation efforts are usually a function of state or political boundaries and thus we argue that understanding patterns at the state level is pragmatic.

Acknowledgements

We thank L. Wilson, J. Sigala, and an anonymous reviewer for helpful comments on the manuscript. Support for this study was provided by Dirección General de Asuntos del Personal Académico – Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica (DGAPA-PAPIIT) through Project IN202021.

References

  • Amori G, Chiozza F, Patterson BD, Rondinini C, Schipper J, Luiselli L (2013) Species richness and distribution of Neotropical rodents, with conservation implications. Mammalia 77(1): 1–19. https://doi.org/10.1515/mammalia-2012-0050
  • Amori G, Boitani L, Milana G, Maiorano L, Luiselli L (2019) Endemism and diversity in European montane mammals: Macro-ecological patterns. Biological Journal of the Linnean Society. Linnean Society of London 128(1): 225–237. https://doi.org/10.1093/biolinnean/blz081
  • AmphibiaWeb (2020) AmphibiaWeb. University of California, Berkeley, CA, USA. https://amphibiaweb.org [Accessed on 2020-11-30]
  • Assunção-Albuquerque MJT, Rey Benayas JM, Rodríguez MÅ, Albuquerque FS (2012) Geographic patterns of vertebrate diversity and identification of relevant areas for conservation in Europe. Animal Biodiversity and Conservation 35(1): 1–11. https://doi.org/10.32800/abc.2012.35.0001
  • Barnagaud J-Y, Geniez P, Cheylan M, Crocket PA (2021) Climate overrides the effects of land use on the functional composition and diversity of Mediterranean reptile assemblages. Diversity & Distributions 27(1): 50–64. https://doi.org/10.1111/ddi.13176
  • Becker CG, Loyola RD (2008) Extinction risk assessments at the population and species level: Implications for amphibian conservation. Biodiversity and Conservation 17(9): 2297–2304. https://doi.org/10.1007/s10531-007-9298-8
  • Chacón-Prieto F, Rodríguez-Soto C, Cuervo-Robayo AP, Carbajal Monroy JC, Alagador D (2021) Protected areas in Central Mexico – are they fit in promoting species persistence under climate and land use changes? Biological Conservation 260: е109186. https://doi.org/10.1016/j.biocon.2021.109186
  • Challenger A, Dirzo R, López Acosta JC, Mendoza E, Lira-Noriega A, Cruz I, Flores Martínez A, González Espinosa M (2009) Factores de cambio y estado de la biodiversidad. In: Capital Natural de México, vol. II: Estado de Conservación y Tendencias de Cambio. CONABIO, Mexico City, 37–73.
  • Chen Y, Peng S (2017) Evidence and mapping of extinction debts for global forest-dwelling reptiles, amphibians and mammals. Scientific Reports 7(1): е44305. https://doi.org/10.1038/srep44305
  • Cordier JM, Aguilar R, Lescano JN, Leynaud GC, Bonino A, Miloch D, Loyola R, Nori J (2021) A global assessment of amphibian and reptile response to land-use changes. Biological Conservation 253: е108863. https://doi.org/10.1016/j.biocon.2020.108863
  • Cuervo-Robayo AP, Ureta C, Goméz-Albores MA, Meneses-Mosquero AK, Téllez-Valdés O, Martínez-Meyer E (2020) One hundred years of climate change in Mexico. PLoS ONE 15(7): e0209808. https://doi.org/10.1371/journal.pone.0209808
  • Doherty TS, Balouch S, Bell K, Burns TJ, Feldman A, Fist C, Garvey TF, Jessop TS, Meiri S, Driscoll DA (2020) Reptile responses to anthropogenic habitat modification: A global meta-analysis. Global Ecology and Biogeography 29(7): 1265–1279. https://doi.org/10.1111/geb.13091
  • Falaschi M, Manenti R, Thuiller W, Ficetola GF (2019) Continental-scale determinants of population trends in European amphibians and reptiles. Global Change Biology 25(10): 3504–3515. https://doi.org/10.1111/gcb.14739
  • Flores-Villela O, Canseco-Márquez L, Ochoa-Ochoa LM (2010) Geographic distribution and conservation of the Mexican Central Highlands herpetofauna. In: Wilson LD, Townsend JH, Johnson JD (Eds) Conservation of Mesoamerican Amphibians and Reptiles. Eagle Mountain Publishing, Utah, 303–321.
  • Frost DR (2020) Amphibian Species of the World: an Online Reference. Version 6.1 (10 January 2020). American Museum of Natural History, New York. https://doi.org/10.5531/db.vz.0001
  • Geldmann J, Joppa LN, Burgess ND (2014) Mapping change in human pressure globally on land and within protected areas. Conservation Biology 28(6): 1604–1616. https://doi.org/10.1111/cobi.12332
  • González-Sánchez VH, Johnson JD, González-Solís D, Fucsko LA, Wilson LD (2021) A review of the introduced herpetofauna of Mexico and Central America, with comments on the effects of invasive species and biosecurity methodology. ZooKeys 1022: 79–154. https://doi.org/10.3897/zookeys.1022.51422
  • Hu X, Huang B, Verones F, Cavalett O, Cherubin F (2021) Overview of recent land-cover changes in biodiversity hotspots. Frontiers in Ecology and the Environment 19(2): 91–97. https://doi.org/10.1002/fee.2276
  • INEGI - Instituto Nacional de Estadística y Geografía (México) (2018) Anuario Estadísticos y Geográfico por Entidad Federative. Instituto Nacional de Estadistica Geografía, México.
  • Jiménez-García D, Li X, Lira-Noriega A, Townsend Peterson A (2021) Upward shifts in elevational limits of forest and grassland for Mexican volcanoes over three decades. Biotropica 53(3): 798–807. https://doi.org/10.1111/btp.12942
  • Kafesh A, Ashrafi S, Yousefi M, Rastegar-Pouyani E, Rajabizadeh M, Ahmadzadeh F, Grünig M, Pellissier L (2020) Reptile species richness associated to ecological and historical variables in Iran. Scientific Reports 10(1): e18167. https://doi.org/10.1038/s41598-020-74867-3
  • Lara-Tufiño JD, Badillo-Saldaña LM, Hernández-Austria R, Ramírez-Bautista A (2019) Effects of traditional agroecosystems and grazing areas on amphibian diversity in a region of central Mexico. PeerJ 7: e6390. https://doi.org/10.7717/peerj.6390
  • Lazcano D, Nevárez-de los Reyes M, García-Padilla E, Johnson JD, Mata-Silva V, DeSantis DL, Wilson LD (2019) The herpetofauna of Coahuila, Mexico: Composition, distribution, and conservation status. Amphibian & Reptile Conservation 13: 31–94.
  • Londoño-Murcia MC, Sánchez-Cordero V (2011) Distribución y conservación de especies amenazadas en Mesoamérica, Chocó y Andes tropicales. Revista Mexicana de Biodiversidad 82(3): 926–950. https://doi.org/10.22201/ib.20078706e.2011.3.1214
  • Lorenzo C, Sántiz EC, Bolaños-Citalán J, Navarrete-Gutiérrez D (2017) Influences of long-term changes in land cover on mammal populations: An example from Mexico. Oryx 53(2): 321–328. https://doi.org/10.1017/S0030605317000485
  • Lu Y, Yang Y, Sun B, Yuan J, Yu M, Stenseth NC, Bullock JM, Obersteiner M (2020) Spatial variation in biodiversity loss across China under multiple environmental stressors. Science Advances 6(47): eebd0952. https://doi.org/10.1126/sciadv.abd0952
  • Masés-García CA, Briones-Salas M, Sosa-Escalante JE (2021) Assessment of wildlife crime in a high-biodiversity region of Mexico. Journal for Nature Conservation 59: e125932. https://doi.org/10.1016/j.jnc.2020.125932
  • Mayani-Parás F, Botello F, Castañeda S, Sánchez-Cordero V (2019) Impact of habitat loss and mining on the distribution of endemic species of amphibians and reptiles in Mexico. Diversity (Basel) 11(11): e210. https://doi.org/10.3390/d11110210
  • Mayani-Parás F, Botello F, Castañeda S, Munguía-Carrana M, Sánchez-Cordero V (2021) Cumulative habitat loss increases conservation threats on endemic species of terrestrial vertebrates in Mexico. Biological Conservation 253: e108864. https://doi.org/10.1016/j.biocon.2020.108864
  • Murali G, Gumbs R, Meiri S, Roll U (2021) Global determinants and conservation of evolutionary and geographic rarity in land vertebrates. Science Advances 7(42): eabe5582. https://doi.org/10.1126/sciadv.abe5582
  • Pavia EG, Graef F, Reyes J (2009) Annual and seasonal surface air temperature trends in Mexico. International Journal of Climatology 29(9): 1324–1329. https://doi.org/10.1002/joc.1787
  • Pianka E (1966) Latitudinal gradients in species diversity: A review of concepts. American Naturalist 100(910): 33–46. https://doi.org/10.1086/282398
  • Pontarp M, Bannefeld L, Sarmento Cabral J, Etienne RS, Fritz SA, Gillespie R, Graham CH, Hagen O, Hartig F, Huang S, Jansson R, Maliet O, Münkemüller T, Pellissier L, Rangel TF, Storch D, Wiegand T, Hurlbert AH (2019) The latitudinal diversity gradient: Novel understanding through mechanistic eco-evolutionary models. Trends in Ecology & Evolution 34(3): 211–223. https://doi.org/10.1016/j.tree.2018.11.009
  • Powney GD, Grenyer R, Orme CDL, Owens LPF, Meiri S (2010) Hot, dry and different: Australian lizard richness is unlike that of mammals, amphibians, and birds. Global Ecology and Biogeography 19(3): 386–396. https://doi.org/10.1111/j.1466-8238.2009.00521.x
  • Ramírez-Bautista A, Hernández-Salinas U, Cruz-Elizalde R, Berriozabal-Islas C, Moreno-Lara I, DeSantis DL, Johnson JD, García-Padilla E, Mata-Silva V, Wilson LD (2020) The herpetofauna of Hidalgo, Mexico: Composition, distribution, and conservation status. Amphibian & Reptile Conservation 14: 63–118.
  • Rivera JA, Rich HN, Lawing AM, Rosenberg MS, Martins EP (2021) Occurrence data uncover patterns of allopatric divergence and interspecies interactions in the evolutionary history of Sceloporus lizards. Ecology and Evolution 11(6): 2796–2813. https://doi.org/10.1002/ece3.7237
  • Rodrigues ASL, Brooks TM, Butchart SHM, Chanson J, Cox N, Hoffmann M, Stuart SN (2014) Spatially explicit trends in the global conservation status of vertebrates. PLoS ONE 9(11): e113934. https://doi.org/10.1371/journal.pone.0113934
  • Rodríguez MA, Belmontes JA, Hawkins BA (2005) Energy, water and large-scale patterns of reptile and amphibian species richness in Europe. Acta Oecologica 28(1): 65–70. https://doi.org/10.1016/j.actao.2005.02.006
  • Sarukhán J, Urquiza-Haas T, Koleff P, Carabias J, Dirzo R, Ezcurra E, Cerdeira-Estrada S, Soberón J (2015) Strategic actions to value, conserve, and restore the natural capital of megadiversity countries: The case of Mexico. Bioscience 65(2): 164–173. https://doi.org/10.1093/biosci/biu195
  • SEMARNAT (2019) Modificación al anexo normative III, lista de especies en riesgo de la Norma Oficial Mexicana NOM-059-SEMARNAT-2010. Protección Ambiental-Especies nativas de México de flora y fauna silvestres- Categorías de riesgo y especificaciones para su inclusion, exclusion o cambio- Lista de especies en riesgo, publicada el 30 de diciembre de 2010 (14 noviembre 2019). Secrearía de Medio Ambiente y Recursos Naturales, México. https://www.dof.gob.mx/nota_detalle.Php?codigo=5578808&fecha=14/11/2019.
  • Senior AF, Böhm M, Johnstone CP, McGee MD, Meiri S, Chapple DG, Tingley R, Byrne M (2021) Correlates of extinction risk in Australian squamate reptiles. Journal of Biogeography 48(9): 2144–2152. https://doi.org/10.1111/jbi.14140
  • Téllez-Valdés O, Dávila-Aranda P, Lira-Saade R (2006) The effects of climate change on the long-term conservation of Fagus grandifolia var. mexicana, an important species of the Cloud Forest in eastern Mexico. Biodiversity and Conservation 15(4): 1095–1107. https://doi.org/10.1007/s10531-004-1868-4
  • Udy K, Fritsch M, Meyer KM, Grass I, Hanß S, Hartig F, Kneib T, Kreft H, Kukunda CB, Pe’er G, Reininghaus H, Tietjen B, Tscharntke T, van Waveren C-S, Wiegand K (2021) Environmental heterogeneity predicts global species richness patterns better than area. Global Ecology and Biogeography 30(4): 842–851. https://doi.org/10.1111/geb.13261
  • Villaseñor JL, Ortiz E, Delgadillo-Moya C, Juárez D (2020) The breadth of the Mexican Transition Zone as defined by its flowering plant genetic flora. PLoS ONE 15(6): e0235267. https://doi.org/10.1371/journal.pone.0235267
  • Whiting ET, Fox DL (2021) Latitudinal and environmental patterns of species richness in lizards and snakes across continental North America. Journal of Biogeography 48(2): 291–304. https://doi.org/10.1111/jbi.13996
  • Wilson LD, Johnson JD (2010) Distributional patterns of the herpetofauna of Mesoamerica, a biodiversity hotspot. In: Wilson LD, Townsend JH, Johnson JD (Eds) Conservation of Mesoamerican Amphibians and Reptiles. Eagle Mountain Publishing, Utah, 31–235.

Supplementary materials

Supplementary material 1 

Table S1

Geoffrey R. Smith, Julio A. Lemos-Espinal

Data type: Summary of sources for data used.

Explanation note: This file contains a list of the checklists from which we obtained our species lists for each Mecixan state. It also provides information on sources used to update those checklists.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (39.12 kb)
Supplementary material 2 

Table S2

Geoffrey R. Smith, Julio A. Lemos-Espinal

Data type: Occurrences.

Explanation note: This file contains the data matrix used to generate the data used in our analyses. Basically it is a matrix of presence and absence of herpetofaunal species in each of the considered Mexican states. This matrix was created using the sources provided in Table S1.

This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.
Download file (199.87 kb)
login to comment