Diversity patterns of Mexican land and freshwater snails: a spatiotemporal approach

Patrones de diversidad de los caracoles terrestres y dulceacuícolas mexicanos: una aproximación espacio-temporal

Resumen

Taxonomic distribution patterns arise from interactions among historical, ecological, and stochastic processes (Halffter & Morrone, 2017; Haydon et al., 1994; Wiens & Donoghue, 2004). The analysis of these events within a temporal framework may explain how geological and climatic changes promote diversification patterns (Bryson et al., 2012; Martínez-Méndez et al., 2019). The integration of both historical and ecological components is required to explain species richness patterns and community composition at local and regional scales (Holloway, 2003; Wiens & Donoghue, 2004).

In historical biogeography, the shared patterns among organisms with different dispersal capabilities allow to evaluate the relationships between different geographic units through identifying generalized tracks (GTs), which arise from the geographic congruence of different individual tracks (ITs). GTs suggest the existence of ancestral biotic components that were affected by vicariant events (Morrone, 2001). Furthermore, panbiogeographic nodes (PNs) that arise as a product of the convergence of individual or generalized tracks represent regions where different biogeographic histories overlap, indicating areas of high species richness (i.e., hotspots) that could lead to the establishment of conservation policies (Morrone, 2001). Analyses of spatial diversity patterns at different scales (e.g., alpha diversity, or species richness within a community; gamma diversity, or species richness within a landscape; and beta diversity, or variation in species composition among communities within a landscape) are fundamental for understanding community patterns and processes (Hawkins et al., 2003). Such analyses also link local community structure and composition to regional and continental assemblages (Cornell & Lawton, 1992).

Studies in Mexico have evaluated diversity patterns using many different criteria (Espinosa-Organista et al., 2000; Morrone et al., 2017; Olson et al., 2001). The existence of general distribution patterns can be determined from geological history, climatic changes, GTs, and PNs (Aguilar-Aguilar et al., 2008; Contreras-Medina et al., 2007; Corona et al., 2009; Escalante et al., 2018; Huidobro et al., 2006). On the other hand, by using grid-cells in analyses of spatial diversity, researchers have recognized that most species occupy small areas; finally, other studies highlighted the role of the Mexican Transition Zone and the occurrence of diversity hotspots (Aguilar-Aguilar et al., 2008; Koleff et al., 2008; Ochoa-Ochoa et al., 2012).

Globally, continental gastropods are one of the most diverse animal groups; their low vagility and dispersal potential, combined with their high endemicity, make them an ideal model for biogeographic studies (Hylander et al., 2005). However, beyond a few regional inventories, no studies have addressed the diversity patterns of continental gastropods in Mexico (Correa-Sandoval et al., 2007; Naranjo-García, 2003; Naranjo-García & Fahy, 2010; Van Devender et al., 2012). Diversity patterns in continental gastropods are different compared to other terrestrial taxa such as vertebrates, besides, the patterns differ between land and freshwater snails, in the case of land snails, most of the diversity is concentrated in islands and limestone mountain ranges (Solem, 1984); on the other hand, freshwater snails richness is concentrated in springs and groundwaters, large rivers and their tributaries, ancient oligotrophic lakes and monsoonal wetlands (Strong et al., 2008). In this study, a comprehensive biogeographic approach was conducted to identify GTs and PNs and evaluate spatial diversity patterns in Mexican continental gastropods as a whole; additionally, we evaluate land and freshwater snails regardless to evaluate if there are similarities in their distribution patterns despite their different environments. 

Materials and methods

Records of continental gastropods in Mexico were gathered from a checklist by Thompson (2011), an updated checklist of freshwater snails by Czaja et al. (2020), and the GBIF database (http://www.gbif.org/). In addition, a search in the available literature was carried out and new species described from 2008 to the present were incorporated. Collection records with missing coordinates were georeferenced using maps and gazettes. Introduced species recognized by Thompson (2011) and new introduced species reported by Naranjo-García and Castillo-Rodríguez (2017), as well as ambiguous and duplicated geographic records, were not considered. The taxonomic status of each species was verified using Molluscabase (http://www.molluscabase.org/) and original descriptions. The below mentioned analyses were performed considering 3 sets of data, the complete continental gastropods fauna and land snails and freshwater snails separately. 

The majority of gastropod species in Mexico have restricted distributions, and most of them are only known from their type locality. Nevertheless, a panbiogeographic analysis involve the comparative study of several unrelated taxa and the integration of these with the tectonic history (Heads, 2005). For this reason, panbiogeographic track analysis was performed at the genus level as has been used in other studies where the number of records diminish considerably if used at species level (Palma-Ramírez et al., 2014). Individual tracks were obtained using minimum spanning trees constructed from the geographic distances among collection records and the Prim algorithm in the ‘fossil’ package (Vavrek, 2011) in R (R Development Core Team, 2017); this analysis was based on the methods proposed by Escalante et al. (2017). Generalized tracks (GTs) were obtained using a parsimony analysis of endemicity with progressive character elimination (PAE-PCE; Echeverry & Morrone, 2010). The 14 biogeographic provinces defined by Morrone et al. (2017) were used as geographic units. ITs were codified as multistate non additive characters (absent = 0, present = 1, and crossing through the province = 2). Phylogenetic analyses were performed in TNT version 1.1 (Goloboff et al., 2008) using a traditional search based on 1,000 replicates of random sequence addition and branch exchange via tree bisection and reconnection (TBR). Ten trees were retained in each run, and a final strict consensus tree was generated from the most parsimonious trees. GTs in the tree were identified as nodes supported by at least 2 geographic synapomorphies (ITs) shared by a biogeographic province. This procedure was iterated until no more GTs were recovered. Finally, GTs and panbiogeographic nodes (PNs) were plotted using the extension Trazos2004 (Rojas-Parra, 2007) distinguishing between individual nodes (IN) that represent the intersection of at least 2 individual tracks and generalized nodes (GN) that represent the intersection of 2 generalized tracks (Morrone, 2015).

The geographic area of Mexico was divided into 184 grid-cells of 1° latitude by 1° longitude to control for the effect of area size on diversity values (Koleff et al., 2008). One-degree cells were used instead of the 30-second grid-cells used in other studies to avoid a large number of empty cells where no gastropod records are available. A presence-absence matrix (PAM) was constructed from the 184 grid-cells and the species recorded in this study. Alpha diversity and average proportional range size were generated (Koleff et al., 2008) and edited in QGIS version 3.14.0 (Team, 2016), meanwhile, the Shannon index was calculated in DIVA-GIS version 7.5.0 (Hijmans et al., 2001). Cells with proportion of species higher than 10% were recognized as diversity hotspots for the continental gastropods, meanwhile, for the land and freshwater snails datasets the threshold value was fixed at 7% due to lower number of species recorded. After it, the cells were intersected with the terrestrial ecoregions of the world and the freshwater ecoregions of the world to have a more precise location of the hotspots (Abell et al., 2008; Olson et al., 2001). 

Finally, to compare the spatial distribution of continental gastropods, as well as land and freshwater snails among the Nearctic and Neotropical regions and the Mexican Transition Zone (Morrone et al., 2017), we recorded the number of total and endemic species in each region. 

Results

A total of 4,379 records for 1,075 continental gastropod species were recovered, of them, 890 were terrestrial species and 185 were freshwater species. The species were from 187 genera in 63 families (Supplementary material Table S1, available at: https://github.com/HOmarMejiaG/diversity-patterns). 

The panbiogeographic analysis for the continental gastropods included 159 ITs from 187 recorded genera; 28 of the genera were restricted to a single biogeographic province. Three PAE-PCE runs were performed until no more GTs were recovered. The first PAE-PCE run produced 1 tree with a length (L) of 786, a consistency index (CI) of 0.40, and a retention index (RI) of 0.47. The second run produced 1 tree with a L of 713, a CI of 0.41, and a RI of 0.48. The last run produced 1 tree with a L of 905, a CI of 0.30, and a RI of 0.18. From these 3 cladograms, 5 GTs (GT1-5) were generated (Table 1, Fig. 1). In the Nearctic region, GT1 included the Sonora biogeographic province and was supported by the occurrence of the genera Eremarionta and Helmintoglypta (Table 1, Fig. 1A). GT3 was located between the California and Baja California biogeographic provinces and was supported by the occurrence of 5 genera (Table 1, Fig. 1B). GT2 and GT4 were present throughout the Nearctic and Neotropical regions, as well as the Mexican Transition Zone. GT2 crossed 10 biogeographic provinces and was supported by 7 genera (Table 1, Fig. 1A); GT4 was supported by 4 genera and was present in 12 of the 14 biogeographic provinces, absent only in the Balsas Basin and Sierra Madre del Sur provinces (Table 1, Fig. 1B). GT5 was supported by 8 genera, being the only track that is supported by the distribution of freshwater gastropods (Aroapyrgus and Balconorbis) and was located near the Atlantic slope; it included the Tamaulipas, Veracruzan, Sierra Madre Oriental, and Transmexican Volcanic Belt provinces (Table 1, Fig. 1C). A total of 86 individual nodes (IN) were identified: 13 in the Nearctic zone, 36 in the Mexican Transition Zone, and 37 in the Neotropical zone, on the other hand, a single generalized node (GN) was recovered by the intersection of the Generalized Tracks 1 and 4 in northwestern Mexico (Fig. 1D). The analysis for the land snails was performed with 149 genera and the results were similar to that obtained with the complete dataset (data not shown), finally, for the analysis of freshwater taxa a total of 38 genera were used, a single tree was recovered for the PAE-PCE analysis with a length (L) of 137, a consistency index (CI) of 0.44, and a retention index (RI) of 0.47. Two generalized tracks were recovered, the first extends from southern Chiapas towards Tabasco and Veracruz and is supported by the distribution of 3 genera, meanwhile, the second track extends from northern Mexico through the Atlantic coast to the Yucatán Peninsula and is supported by the distribution of 12 genera, only 2 individual nodes were recovered (Supplementary material Table S2, Fig. S1, available at: https://github.com/HOmarMejiaG/diversity-patterns).

 

Table 1

Number of generalized tracks (GT) recovered in this study of the 3 iterations performed on PAE-PCE analysis. The IT represent the name of the genera that support each one of the 3 the generalized tracks

Iteration	GT	IT
1	1	Eremarionta, Helminthoglypta
	2	Allopeas, Drymaeus, Euglandina, Guppya, Helicina, Orthalicus, Pupisoma
2	3	Binneya, Bulimulus, Haplotrema, Herpeteros, Nearctula
	4	Echinix, Hawaiia, Polygyra, Thysanophora
3	5	Aroapyrgus, Balconorbis, Chondropoma, Epirobia, Gastrocopta, Microconus, Pallifera, Pycnogyra

Of the 189 grid-cells, 13 (7%) of them had no species records. Most of the other grid-cells were characterized by low diversity values (<6% of total species). However, 3 grid-cells contained over 10% of continental gastropod species and could be considered hotspots. The first hotspot contained 108 species (10% of the total) and was located in the Veracruz moist forest and the Oaxacan montane forest ecoregions. The second hotspot is located in Veracruz moist forest and Veracruz montane forest ecoregions containing 130 species (12% of the total). The third hotspot was located in the Oaxaca montane forest, Sierra Madre de Oaxaca oak-pine forest and Petén Veracruz moist forest ecoregions and include 229 species (21%), a fifth of Mexico’s continental gastropod diversity (Fig. 2A). For the land snails, 4 grid cells could be considered as hotspots, the number of species ranged from 7.1% to 13.4%, the location of the grid cells match the previously described hotspots, but an additional one is located in Tehuacán Valley matorral and Sierra Madre de Oaxaca pine-oak forest ecoregions (Supplementary material Fig. S2, available at: https://github.com/HOmarMejiaG/diversity-patterns). In the freshwater snails, 3 hotspots were recognized, the first one contained 18 species (9.7 % of the total) and was located in Pánuco tropical and subtropical coastal rivers freshwater ecoregion, the second hotspot contained 15 species (8.1%) and was located in the Río Conchos xeric freshwater and endorheic basins freshwater ecoregion, finally, the third hotspot comprised 13 species (7.1%) and is located in the Grijalva-Usumacinta tropical and subtropical coastal rivers freshwater ecoregion (Supplementary material Fig. S3, available at: https://github.com/HOmarMejiaG/diversity-patterns) 

The Shannon diversity values tended to be more geographically homogeneous, but this analysis indicates grid-cells with high diversity values (3.6-4.4) in Hidalgo, Puebla, Queretaro, Nuevo León, Coahuila, Oaxaca, and Veracruz for the complete dataset and the land snails. On the other hand, for the freshwater snails the highest diversity cells (2.1-2.6) are located in Coahuila, San Luis Potosí and Veracruz (Supplementary material Fig. S4, available at: https://github.com/HOmarMejiaG/diversity-patterns)

Of the 1,075 continental gastropod species studied, 814 (75.7%) are endemic to Mexico, but only 15 species (1.4%) (all terrestrial) had widespread distributions, occupying 10 or more grid-cells: Drymaeus multilineatus (10 grid-cells), Praticolella berlandieriana (10), Bulimulus unicolor (10), Allopeas beckianum (11), Choanopoma largillierti (11), Rabdotus schiedeanus (11), Gastrocopta pellucida (11), Drymaeus emeus (12), Helicina chrysocheila (15), Neocyclotus dysoni (15), Orthalicus princeps (17), Thysanophora horni (17), Succinea luteola (18), Hawaiia minuscula (19) and Rabdotus alternatus (21). Meanwhile, 433 (39.8%) of the species occupied less than 10 grid-cells, and the majority of species (586, or 58.7%) are restricted to only 1 grid-cell (Fig. 2B; Supplementary material Table S1, available at: https://github.com/HOmarMejiaG/diversity-patterns)

The Nearctic region comprises 5 biogeographic provinces with an estimated area of 959,205 km² and includes 280 species of continental gastropods of which 204 are endemics; the Neotropical region comprises 4 biogeographic provinces with an estimated area of 433,558 km² and include 452 species of continental gastropods of which 326 are endemics, finally, the Mexican Transition Zone comprises 5 biogeographic provinces with an estimated area of 591,901 km² and include 343 species of continental gastropods of which 284 are endemics. For the land snails, the Nearctic region included 250 species (156 endemic), the Neotropical region 460 species (281 endemic) and the MTZ 185 species (175 endemic); meanwhile, for the freshwater snails, the Nearctic region included 50 species (45 endemic), the Neotropical region 90 (44 endemic) and the MTZ 40 species (13 endemic).

Discussion

There are several critiques for PAE as a panbiogeographic method for recovering historical biogeographic patterns, such as the method does not take into account the phylogenetic relationships (Humphries, 1989), the use of a posteriori optimality criterion (Szumik, 2002) and the use of artificial delimited areas (Peterson, 2008), among others. Nevertheless, most of the critiques arise from the use of PAE as a method for recovery of historic relationships among areas and not as a method for the recognition of areas of endemism (Nihei, 2006). Areas of endemism are strong evidence of the historical structure of biotic components and represent a basic step in evolutionary biogeography, where the PAE PCE method appears as the best approach for the identification of generalized tracks (Morrone, 2014).

The biogeographic patterns recovered in this study were similar to those previously established in other biological groups with the same methodological approach of PAE PCE (Contreras-Medina et al., 2007). Beyond agreement with previous research, there are also other pieces of evidence that support the patterns recovered in this study. First, the GTs of northwestern Mexico (GT1 and GT3; Fig. 1, Table 1) were supported by 7 genera of Nearctic affinity restricted to the arid ecosystems of the Sonora, California, and Baja California provinces (Naranjo-García & Fahy, 2010) and are similar to GTs recovered for other taxa (Escalante et al., 2018; Rojas-Soto et al., 2003). The heterogeneous landscape, as well as geological changes that occurred in the Baja California Peninsula 6–15 MYA can explain these patterns (Escalante et al., 2018; Naranjo-García, 1988; Riddle et al., 2000).

Second, the geographic distribution of GT2 and GT4 (Fig. 1) are congruent with the typical Neotropical pattern proposed by Halffter (1987), the lowland Mesoamerican track proposed by Morrone and Márquez (2001), and the patterns recovered in other groups (Escalante et al., 2018; Gonzalez-Ávila et al., 2017; Luna-Vega et al., 2001; Márquez & Morrone, 2003; Rosas-Valdez & Pérez-Ponce de León, 2008). The Neotropical affinity groups may extend beyond the coastal regions and low altitudes, as evidenced by the presence of Neotropical affinity species above 1,400 m such as the genera Allopeas and Drymaeus in GT2 that could reach up to 2,059 and 3,698 m, respectively, and the genera Echinix and Polygyra in GT4 that occur up to 1,485 and 2643 m, respectively. This distribution pattern in highland regions has also been found in other groups, such as volant mammals, and is characteristic of a transitional fauna that may have arisen either by the dominance of endemic or widely distributed species or by a balance between the number of Neartic and Neotropical species in grid cells (Ortega & Arita, 1998). 

Furthermore, the support of GT2 by the genera Drymaeus, Euglandina, and Orthalicus demonstrate the role of Sierra Madre del Sur, the Sierra Madre Occidental, and the Transmexican Volcanic Belt as important biogeographic barriers to dispersal for several terrestrial gastropod families (Mejía & Zúñiga, 2007; Naranjo-García & Fahy, 2010). 

The panbiogeographic tracks recovered for the freshwater snails were similar to the tracks recovered for the land snails; although an incomplete species sampling over the whole country could lead to misleading historical patterns (Brooks & Van Veller, 2003; Nihei, 2006), our results support that a panbiogeographic analysis is useful as a preliminary approach to state primary biogeographic homology hypothesis in groups where phylogenetic analysis are scarce.

The species richness in tropical regions is explained by both ecological and historical factors. The number of species recovered in this paper in the Neotropical region (5 times greater than the Nearctic region and 2.5 times higher than the MTZ) were also recovered in other biological groups such as different phyla of freshwater organisms (Quiroz-Martínez et al., 2014) and mammals (Gómez-Ortiz et al., 2017), supporting that the Neotropical region comprises the higher diversity in Mexican fauna. Additionally, the results of this study highlight the complexity of the MTZ (Corona et al., 2009; Escalante et al., 2018; Halffter & Morrone, 2017) that has been considered either as a high richness species area (Aguilar-Méndez et al., 2021) or as an area of high phylogenetic diversity, (Gómez-Ortiz et al., 2017) that are reflected in the high number of total (343) and endemic (284) species of continental gastropods. Thus, the Mexican Transvolcanic Belt provide a set of factors that are related with high speciation rates in land snails such as a topographic mosaic, rain shadow effects, heat budget differences, vegetation cover and effective refugia from floods that serve as shelter areas (Solem, 1984). The diversity hotspots of freshwater snails recovered in this study were similar to the previously reported by Czaja et al. (2020), for Cuatro Cienégas-Sabinas system and the Río Verde-Guayalejo, with the exception of the Río Conchos-Guzmán system hotspot that has not been recovered in this work, although we were able to recovered a grid in the Grijalva that has been postulated as a potential hotspot (Czaja et al., 2020). The main difference between both approaches is that Czaja et al. (2020) used the complete freshwater ecoregions meanwhile in this work we used grid cells that reduce the size area effect allowing a more precise comparison. 

In particular, the nodes gathered in this study were geographically similar to the ones proposed by Morrone and Márquez (2008), who defined 15 nodes in the boundaries of some Mexican biogeographic provinces; they are also similar to 3 of the 5 freshwater snails endemicity hotspots suggested by Czaja et al. (2020) and to the highest species diversity areas for terrestrial gastropods recovered in previous regional studies (Correa-Sandoval et al., 2007; Naranjo-García & Fahy, 2010).

In this study, a single GN was recovered in Northwestern Mexico and only 2 INs were recovered for the Mexican Pacific slope. In contrast, 2 regions along the Gulf of Mexico basin had very high node concentrations. The intersection of the Tamaulipas, Sierra Madre Oriental, Veracruzan, and Transmexican Volcanic Belt biogeographic provinces had up to 60 INs; this coincides with the proposal by Correa-Sandoval (2003), who suggested the Sierra Madre Oriental malacological province for its distinctive fauna. The nodes identified in the central area of the Veracruzan province are geographically similar to the endemicity hotspots for freshwater gastropods in the Río Verde, Río Guayalejo, and Río Papaloapan regions, as previously defined by Czaja et al. (2020). 

The results of spatial richness patterns obtained in this study are similar to those recovered in other taxa (Aguilar-Aguilar et al., 2008; Koleff et al., 2008; Ochoa-Ochoa et al., 2012). However, only 3 grid-cells contained more than 10% of the species diversity, which is in contrast to previous research on vertebrates, where there were a large number of grid-cells with at least these diversity values (Koleff et al., 2008). This could be the product of environmental factors, which have a considerable effect on species richness patterns for taxa with low dispersal rates because they have geographically restricted distributions and high levels of endemism (Clements et al., 2008). Finally, the INs and highest species diversity values detected in this study coincide with conservation priority areas previously defined either by species richness or by the occurrence of biotic elements of different origins (Corona et al., 2009). 

Acknowledgments

This work was supported by SIP project number 2021044. We appreciate the comments of Enrique Ortiz, two anonymous reviewers and the comments of the Editor Fernando Álvarez that allowed us to improve the manuscript. 

References

Abell, R. A., Thieme M., & Revenga, C. (2008). Freshwater ecoregions of the world: a new map of ecoregions of the world: a new map of biogeographic units for freshwater biodiversity conservation. Bioscience, 58, 403 –414. https://doi.org/10.1641/b580507 

Aguilar-Aguilar, R., Salgado-Maldonado, G., Contreras-Medina, R., & Martínez-Aquino, A. (2008). Richness and endemism of helminth parasites of freshwater fishes in Mexico. Biological Journal of the Linnean Society, 94, 435 –444. https://doi.org/10.1111/j.1095-8312.2008.00994.x 

Aguilar-Méndez, M. J., Rosas-Mejía, M., Vásquez-Bolaños, M., González-Hernández, G. A., & Janda, M. (2021). New distributional records for ants and the evaluation of ant species richness and endemism patterns in Mexico. Biodiversity Data Journal, 9, e60630. https://doi.org/10.3897/bdj.9.e60630 

Brooks, D. R., & Van Veller, M. G. (2003). Critique of parsimony analysis of endemicity as a method of historical biogeography. Journal of Biogeography, 30, 819 –825. https://doi.org/10.1046/j.1365-2699.2003.00848.x 

Bryson Jr, R. W., García-Vázquez, U. O., & Riddle, B. R. (2012). Diversification in the Mexican horned lizard Phrynosoma orbiculare across a dynamic landscape. Molecular Phylogenetics and Evolution, 62, 87–96. https://doi.org/10.1016/j.ympev.2011.09.007 

Clements, R., Ng, P. K., Lu, X. X., Ambu, S., Schilthuizen, M., & Bradshaw, C. J. (2008). Using biogeographical patterns of endemic land snails to improve conservation planning for limestone karsts. Biological Conservation, 141, 2751–2764. https://doi.org/10.1016/j.biocon.2008.08.011 

Contreras-Medina, R., Luna-Vega, I., & Morrone, J. J. (2007). Application of parsimony analysis of endemicity to Mexican gymnosperm distributions: grid-cells, biogeographical provinces and track analysis. Biological Journal of the Linnean Society, 92, 405–417. https://doi.org/10.1111/j.1095-8312.2007.00844.x 

Cornell, H. V., & Lawton, J. H. (1992). Species interactions, local and regional processes, and limits to the richness of ecological communities: a theoretical perspective. Journal of Animal Ecology, 61, 1–12. https://doi.org/10.2307/5503 

Correa-Sandoval, A. (2003). Gastrópodos terrestres del noreste de México. Revista de Biología Tropical, 51, 507–522.

Correa-Sandoval, A., Strenth, N. E., & Salazar Rodríguez, M. D. C. (2007). Zoogeografía de los gasterópodos terrestres del sur de Nuevo León, México. Acta Zoológica Mexicana, 23, 143–162. https://doi.org/10.21829/azm.2017.3311015 

Corona, A. M., Toledo, V. H., & Morrone, J. J. (2009). Track analysis of the Mexican species of Buprestidae (Coleoptera): testing the complex nature of the Mexican Transition Zone. Journal of Biogeography, 36, 1730–1738. https://www.jstor.org/stable/40305939 

Czaja, A., Meza-Sánchez, I. G., Estrada-Rodríguez, J. L., Romero-Méndez, U., Sáenz-Mata, J., Ávila-Rodríguez, V. et al. (2020). The freshwater snails (Mollusca: Gastropoda) of Mexico: updated checklist, endemicity hotspots, threats and conservation status. Revista Mexicana de Biodiversidad, 91, 1–22. http://dx.doi.org/10.22201/ib.20078706e.2020.91.2909 

Echeverry, A., & Morrone, J. J. (2010). Parsimony analysis of endemicity as a panbiogeographical tool: an analysis of Caribbean plant taxa. Biological Journal of the Linnean Society, 101, 961-976. https://doi.org/10.1111/j.1095-8312.2010.01535.x 

Escalante, T., Noguera-Urbano, E. A., Pimentel, B., & Aguado-Bautista, O. (2017). Methodological issues in modern track analysis. Evolutionary Biology, 44, 284–293. https://doi.org/10.1007/s11692-016-9401-8 

Escalante, T., Noguera-Urbano, E. A., & Corona, W. (2018). Track analysis of the Nearctic region: Identifying complex areas with mammals. Journal of Zoological Systematics and Evolutionary Research, 56, 466–477. http://orcid.org/0000-0002-9127-8969 

Espinosa-Organista, D., Morrone, J. J., Aguilar, C., & Llorente, J. (2000). Regionalización biogeográfica de México: Provincias bióticas. In J. Llorente Bousquets, E. González, & N. Papavero (Eds.), Biodiversidad, taxonomía y biogeografía de artrópodos de México: hacia una síntesis de su conocimiento, Vol II. (pp. 61–94). Facultad de Ciencias, UNAM, Mexico City.

GBIF.org (2020). Accessed 08 Mar, 2020 from https://doi.org/10.15468/dl.example-doNotCite

Goloboff, P. A., Farris, J. S., & Nixon, K. C. (2008). TNT, a free program for phylogenetic analysis. Cladistics, 24, 774–786. https://doi.org/10.1111/j.1096-0031.2008.00217.x 

Gómez-Ortiz, Y., Domínguez-Vega, H., & Moreno, C. E. (2017). Spatial variation of mammal richness, functional and phylogenetic diversity in the Mexican Transition Zone. Community Ecology, 18, 121–127. https://doi.org/10.1556/168.2017.18.2.1 

Gonzalez-Ávila, A., Contreras-Medina, R., Espinosa, D., & Luna-Vega, I. (2017). Track analysis of the order Gomphales (Fungi: Basidiomycota) in Mexico. Phytotaxa, 316, 22–38. https://doi.org/10.11646/phytotaxa.316.1.2 

Halffter, G. (1987). Biogeography of the montane entomofauna of Mexico and Central America. Annual Review of Entomology, 32, 95–114.

Halffter, G., & Morrone, J. J. (2017). An analytical review of Halffter’s Mexican transition zone, and its relevance for evolutionary biogeography, ecology and biogeographical regionalization. Zootaxa, 4226, 1–46. https://doi.org/10.11646/zootaxa.4226.1.1 

Haydon, D. T., Crother, B. I., & Pianka, E. R. (1994). New directions in biogeography? Trends in Ecology and Evolution, 9, 403–405.

Hawkins, B. A., Field, R., Cornell, H. V., Currie, D. J., Guégan, J. F., Kaufman, D. M. et al. (2003). Energy, water, and broad-scale geographic patterns of species richness. Ecology, 84, 3105-3117. https://doi.org/10.1890/03-8006 

Heads, M. (2005). Towards a panbiogeography of the seas. Journal of Biogeography, 84, 675–723.

Hijmans, R. J., Guarino, L., Cruz, M., & Rojas, E. (2001). Computer tools for spatial analysis of plant genetic resources data: 1. DIVA-GIS. Plant Genetic Resources Newsletter, 127, 15–19.

Holloway, J. D. (2003). Biological images of geological history: through a glass darkly or brightly face to face? Journal of Biogeography, 30, 165–179. https://doi.org/10.1046/j.1365-2699.2003.00826.x 

Huidobro, L., Morrone, J. J., Villalobos, J. L., & Álvarez, F. (2006). Distributional patterns of freshwater taxa (fishes, crustaceans and plants) from the Mexican Transition Zone. Journal of Biogeography, 33, 731–741. https://doi.org/10.1111/j.1365-2699.2005.01400.x 

Humphries, C. J. (1989). Any advance on assumption 2? Journal of Biogeography, 16, 101–102.

Hylander, K., Nilsson, C., Gunnar-Jonsson, B., & Göthner, T. (2005). Differences in habitat quality explain nestedness in a land snail meta-community. Oikos, 108, 351–361. https://doi.org/10.1111/j.0030-1299.2005.13400.x 

Koleff, P., Soberón, J., Arita, H. T., Dávila, P., Flores-Villela, O., Golubov, J. et al. (2008). Patrones de diversidad espacial en grupos selectos de especies. In J. Soberón, G. Halffter, & J. Llorente-Bousquets (comp.) Capital natural de México. Vol. 1. Conocimiento actual de la biodiversidad (pp. 323–346). México D.F.: Comisión Nacional para el Conocimiento y Uso de la Biodiversidad. 

Luna-Vega, I., Morrone, J. J., Ayala, O. A., & Espinosa Organista, D. (2001). Biogeographical affinities among Neotropical cloud forests. Plant Systematics and Evolution, 228, 229–239. https://doi.org/10.1007/s006060170031 

Martínez-Méndez, N., Mejía, O., Ortega, J., & Méndez-de la Cruz, F. (2019). Climatic niche evolution in the viviparous Sceloporus torquatus group (Squamata: Phrynosomatidae). PeerJ, 6, e6192. https://doi.org/10.7717/peerj.6192 

Márquez, J., & Morrone, J. J. (2003). Análisis panbiogeográfico de las especies de Heterolinus y Homalolinus (Coleoptera: Staphylinidae: Xantholinini). Acta Zoológica Mexicana (N.S.), 90, 15–25. https://doi.org/10.21829/azm.2003.902549 

Mejía, O., & Zúñiga, G. (2007). Phylogeny of the three brown-banded land snail genus Humboldtiana (Pulmonata: Humboldtianidae). Molecular Phylogenetics and Evolution, 45, 587–595. https://doi.org/10.1016/j.ympev.2007.07.013 

MolluscaBase (2020). Accessed 08 Jun, 2020 from http://www.molluscabase.org

Morrone, J. J. (2001). Homology, biogeography and areas of endemism. Diversity and Distributions, 7, 297–300. https://doi.org/10.1046/j.1366-9516.2001.00116.x 

Morrone, J. J. (2014). Parsimony analysis of endemicity (PAE) revisited. Journal of Biogeography, 41, 842–854. https://doi.org/10.1111/jbi.12251 

Morrone, J. J. (2015). Track analysis beyond panbiogeography. Journal of Biogeography, 42, 413–425. https://doi.org/10.1111/jbi.12467 

Morrone, J. J., Escalante, T., & Rodríguez-Tapia, G. (2017). Mexican biogeographic provinces: map and shapefiles. Zootaxa, 4277, 277–279. https://doi.org/10.11646/zootaxa.
4277.2.8 

Morrone, J. J., & Márquez, J. (2001). Halffter’s Mexican Transition Zone, beetle generalized tracks, and geographical homology. Journal of Biogeography, 28, 635–650. https://doi.org/10.1046/j.1365-2699.2001.00571.x 

Morrone, J. J., & Márquez, J. (2008). Biodiversity of Mexican terrestrial arthropods (Arachnida and Hexapoda): a biogeographical puzzle. Acta Zoológica Mexicana, 24, 15-41. https://doi.org/10.21829/azm.2008.241613

Naranjo-García, E. (1988). Systematics and biogeography of the Helminthoglyptidae of Sonora (Ph. D. Thesis). Faculty of the Department of Ecology and Evolutionary Biology, The Arizona University. E.U.A.

Naranjo-García, E. (2003). Moluscos continentales de México: dulceacuícolas. Revista de Biología Tropical, 51, 495–505. 

Naranjo-García, E., & Castillo-Rodríguez, Z. G. (2017). First inventory of the introduced and invasive mollusks in Mexico. Nautilus, 131, 107–126.

Naranjo-García, E., & Fahy, N. E. (2010). The lesser families of Mexican terrestrial molluscs. American Malacological Bulletin, 28, 59–80. https://doi.org/10.4003/006.028.0222 

Nihei, S. S. (2006). Misconceptions about parsimony analysis of endemicity. Journal of Biogeography, 33, 2099–2106. https://doi.org/10.1111/j.1365-2699.2006.01619.x 

Ochoa-Ochoa, L. M., Rodríguez, P., Mora, F., Flores-Villela, O., & Whittaker, R. J. (2012). Climate change and amphibian diversity patterns in Mexico. Biological Conservation, 150, 94–102. https://doi.org/10.1016/j.biocon.2012.03.010 

Olson, D. M., Dinerstein, E., Wikramanayake, E. D., Burgess, N. D., Powell, G. V., Underwood, E. C. et al. (2001). Terrestrial ecoregions of the world: a new map of life on earth a new global map of terrestrial ecoregions provides an innovative tool for conserving biodiversity. BioScience, 51, 933–938. https://doi.org/10.1641/0006-3568(2001)051[0933:TEOTWA]2.0.CO;2 

Ortega, J., & Arita, H. T. (1998). Neotropical-Nearctic limits in Middle America as determined by distributions of bats. Journal of Mammalogy, 79, 772–783. https://doi.org/10.2307/1383088 

Palma-Ramírez, A., Goyenechea, I., & Castillo-Cerón, J. M. (2014). Panbiogeography of the Santa María Amajac area, Hidalgo, Mexico. Revista Mexicana de Biodiversidad, 85, 1228–1234. https://doi.org/ 10.7550/rmb.44392

Peterson, A. T. (2008). Parsimony analysis of endemicity (PAE) and studies of Mexican Biogeography. Revista Mexicana de Biodiversidad, 79, 541–542. https://doi.org/10.22201/ib.20078706e.2008.002.544

Quiroz-Martínez, B., Álvarez, F., Espinosa, H., & Salgado-Maldonado, G. (2014). Concordant biogeographic patterns among multiple taxonomic groups in the Mexican freshwater biota. Plos One, e105510. https://doi.org/10.1371/journal.pone.0105510 

R Development Core Team. (2017). R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.

Riddle, B. R., Hafner, D. J., Alexander, L. F., & Jaeger, J. R. (2000). Cryptic vicariance in the historical assembly of a Baja California Peninsular Desert biota. Proceedings of the National Academy of Sciences, 97, 14438–14443. https://doi.org/10.1073/pnas.250413397 

Rosas-Valdez, R., & Pérez-Ponce de León, G. (2008). Composición taxonómica de los helmintos parásitos de ictalúridos y heptaptéridos (Osteichthyes: Siluriformes) de México, con una hipótesis de homología biogeográfica primaria. Revista Mexicana de Biodiversidad, 79, 473–499. https://doi.org/10.22201/ib.20078706e.2008.002.552

Rojas-Soto, O. R., Alcántara-Ayala, O., & Navarro, A. G. (2003). Regionalization of the avifauna of the Baja California Peninsula, Mexico: a parsimony analysis of endemicity and distributional modeling approach. Journal of Biogeography, 30, 449–461. https://doi.org/10.1046/j.1365-2699.2003.00836.x 

Rojas-Parra, C. A. (2007). Una herramienta automatizada para realizar análisis panbiogeográficos. Biogeografía, 1, 31–33.

Solem, A., (1984). A world model of land snail diversity and abundance. In A. Solem, & A. C. Van Bruggen (Eds.), World-wide snails: biogeographical studies on non-marine Mollusca (pp. 6–22). Leiden: E.J. Brill/W. Backhuys.

Strong, E. E., Gargominy, O., Ponder, W. F., & Bouchet, P. (2008). Global diversity of gastropods (Gastropoda; Mollusca) in freshwater. Hydrobiologia, 595, 149–166. https://doi.org/10.1007/s10750-007-9012-6 

Szumik, C. A., Cuezzo, F., Goloboff, P. A., & Chalup, P. E. (2002). An optimality criterion to determine areas of endemism. Systematic Biology, 51, 806–816.

Team, Q. D. (2016). QGIS geographic information system.
Open Source Geospatial Foundation Project. http://qgis.osgeo.org/

Thompson, F. G. (2011). An annotated checklist and bibliography of the land and freshwater snails of Mexico and Central America. Bulletin of the Florida Museum of Natural History, 50, 1–299. 

Van Devender, A. S., Van Devender, R. W., Rivera-García, A., Jiménez-Maldonado, R. E., & Van Devender, M. N. (2012). Report on the terrestrial mollusks of the Sierra de la Madera (Oposura), Sonora, Mexico— The Caracoleros. American Malacological Bulletin, 30, 315-322. https://doi.org/10.4003/006.030.0211 

Vavrek, M. J. (2011). Fossil: palaeoecological and palaeogeographical analysis tools. Palaeontologia Electronica, 14, 1–16.

Wiens, J. J., & Donoghue, M. J. (2004). Historical biogeography, ecology and species richness. Trends in Ecology & Evolution, 19, 639–644. https://doi.org/10.1016/j.tree.2004.09.011