Growth and life history of freshwater chelydrid turtles (Testudines: Cryptodira): A bone histological approach

The current study examines the growth pattern and lifestyle habits of the freshwater snapping turtles Chelydra and Macrochelys based on limb bone histology. Femora, humeri, and tibiae of 25 individuals selected from a range of ontogenetic stages were assessed to determine inter-element and intraskeletal histological variation. Osteohistological assessment of multiple elements is consistent with overall moderate growth rates as revealed by the dominance of parallel-fibered bone. However, the growth was cyclical as shown by deposition of multiple lines of arrested growths in the compacta. It appears that the bone tissue of C. serpentina is more variable through ontogeny with intermittent higher growth rates. M. temminckii appears to grow more slowly than C. serpentina possessing compact and thick cortices in accordance with their larger size. Overall, vascularization decreases through ontogeny with humeri and femora being well-vascularized in both species. Contrarily, epipodials are poorly vascularized, though simple longitudinal and radial canals are present, suggesting differences in growth patterns when compared with associated diaphyseal sections. The tibiae were found to be the least remodeled of the limb bones and therefore better suited for skeletochronology for snapping turtles. Intra-elementally, femora and humeri preserved higher cortical vascularity ventrally, suggestive of faster relative growth. We hypothesize that the differential growth pattern in limb bones of snapping turtles may relate to differential functional constraints, where forelimbs are operational in swimming while the hindlimbs provide stability.

The Phylogenomic and Biogeographic History of the Gartersnakes, Watersnakes, and Allies (Natricidae: Thamnophiini)

North American Thamnophiini (gartersnakes, watersnakes, brownsnakes, and swampsnakes) are an ecologically and phenotypically diverse temperate clade of snakes representing 61 species across 10 genera. In this study, we estimate phylogenetic trees using ∼3,700 ultraconserved elements (UCEs) for 76 specimens representing 75% of all Thamnophiini species. We infer phylogenies using multispecies coalescent methods and time calibrate them using the fossil record. We also conducted ancestral area estimation to identify how major biogeographic boundaries in North America affect broadscale diversification in the group. While most nodes exhibited strong statistical support, analysis of concordant data across gene trees reveals substantial heterogeneity. Ancestral area estimation demonstrated that the genus Thamnophis was the only taxon in this subfamily to cross the Western Continental Divide, even as other taxa dispersed southward toward the tropics. Additionally, levels of gene tree discordance are overall higher in transition zones between bioregions, including the Rocky Mountains. Therefore, the Western Continental Divide may be a significant transition zone structuring the diversification of Thamnophiini during the Neogene and Pleistocene. Here we show that despite high levels of discordance across gene trees, we were able to infer a highly resolved and well-supported phylogeny for Thamnophiini, which allows us to understand broadscale patterns of diversity and biogeography.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

Phylogeography and demographic expansion in the widely distributed horned passalus beetle, Odontotaenius disjunctus (coleoptera: Passalidae)

Dynamic climatic oscillations during the Pleistocene had profound effects on the distribution of species across North America. Although the role of historical climate change on speciation remains controversial, the impact on genetic variation within species has been well documented. Analyses of mtDNA sequences from the cytochrome oxidase I gene (911 bp) for 115 individuals of Odontotaenius disjunctus was combined with ecological niche modelling (ENM) to infer the demographic and population differentiation scenarios under present and past conditions. We inferred three lineages that diverged during the Pleistocene and replace each other geographically across the eastern United States. One of these lineages traverses previously identified genetic barriers for terrestrial animals including the Mississippi and Apalachicola Rivers and the Appalachian Mountains. We observed overlapping ranges between two haplotype groups as well as a region of secondary contact associated with ecological transition zone in northern Florida. The two continental lineages depict a genetic signature of a recent population increase associated with expanding niche envelope, whereas the clade restricted to peninsular Florida shows stable populations in a shrinking niche envelope. Given the lack of ecological separation, overlapping distribution of haplogroups and the presence of secondary contact zones, the taxonomic status of these lineages must await robust testing using multilocus DNA data to assess species boundaries.

Resolving spatial complexities of hybridization in the context of the gray zone of speciation in North American ratsnakes (Pantherophis obsoletus complex)

Inferring the history of divergence between species in a framework that permits the presence of gene flow has been crucial for characterizing the "gray zone" of speciation, which is the period of time where lineages have diverged but have not yet achieved strict reproductive isolation. However, estimates of both divergence times and rates of gene flow often ignore spatial information, for example when considering the location and width of hybrid zones with respect to changes in the environment between lineages. Using population genomic data from the North American ratsnake complex (Pantherophis obsoletus), we connected phylogeographic estimates of lineage structure, migration, historical demography, and timing of divergence with hybrid zone dynamics. We examined the spatial context of diversification by linking migration and timing of divergence to the location and widths of hybrid zones. Artificial neural network approaches were applied to understand how landscape features and past climate have influenced population genetic structure among these lineages. We found that rates of migration between lineages were associated with the overall width of hybrid zones. Timing of divergence was not related to migration rate or hybrid zone width across species pairs but may be related to the number of alleles weakly introgressing through hybrid zones. This research underscores how incomplete reproductive isolation can be better understood by considering differential allelic introgression and the effects of historical and contemporary landscape features on the formation of lineages as well as overall genomic estimates of migration rates through time.

Addressing the Linnean shortfall in a cryptic species complex

Biodiversity is being lost at an alarming rate, but the rate of this loss is likely to be underestimated as a result of a deficit in taxonomic knowledge (i.e. the Linnean shortfall). This knowledge gap is more extensive for morphologically indistinct taxa. The advancement of molecular techniques and delimitation methods has facilitated the identification of such cryptic species, but a majority of these taxa remain undescribed. To investigate the effects of taxonomic uncertainty on understanding of biodiversity, we applied the general lineage concept of species to an amphipod species complex, the Gammaruslacustris lineage that occupies springs of the northern Chihuahuan Desert, which is emerging in contemporary times. We investigated species boundaries using a validation-based approach and examined genetic structure of the lineage using a suite of microsatellite markers to identify independently evolving metapopulations. Our results show that each spring contains a genetically distinct population that is geographically isolated from other springs, suggesting evolutionary independence and status as separate species. Additionally, we observed subtle interspecific morphological variation among the putative species. We used multiple lines of evidence to formally describe four new species (Gammarus langi sp. nov., G. percalacustris sp. nov., G. colei sp. nov. and G. malpaisensis sp. nov.) endemic to the northern Chihuahuan Desert. Cryptic speciation is likely to be high in other aquatic taxa within these ecosystems, and across arid landscapes throughout North America and elsewhere, suggesting that the magnitude of the Linnean shortfall is currently underestimated in desert springs worldwide.

Connections and containers: Using genetic data to understand how watershed evolution and human activities influence cutthroat trout biogeography

Species with large geographic distributions often exhibit complex patterns of diversity that can be further complicated by human activities. Cutthroat trout (Oncorhynchus clarkii) are one of the most widely distributed freshwater fish species in western North America exhibiting substantial phenotypic and genetic variability; however, fish stocking practices have translocated populations outside of their native range and may have obscured intraspecific boundaries. This study focuses on cutthroat trout populations representing three distinct evolutionary clades that are found intermixed within a contact zone between the Bonneville and upper Snake River watersheds in the western United States. We used mitochondrial and microsatellite genetic data, as well as historical stocking records, to evaluate whether populations of cutthroat trout in the contact zone are native or are introduced. We found significant genetic differentiation and fine-scale genetic population structure that was organized primarily by watershed boundaries. While we detected increased genetic diversity in some areas in close proximity to the greatest number of stocking events, the highly organized population structure both within and between areas of the contact zone indicates that the populations are native to the watersheds. Intermixing of distinct evolutionary lineages of cutthroat trout appears to be the result of historical connections between paleodrainages. Our analyses provide a context for understanding how genetic data can be used to assess the status of populations as native or introduced.

Delineating modern variation from extinct morphology in the fossil record using shells of the Eastern Box Turtle (Terrapene carolina)

Characterization of morphological variation in the shells of extant Eastern Box Turtles, Terrapene carolina, provides a baseline for comparison to fossil populations. It also provides an example of the difficulties inherent to recognizing intraspecific diversity in the fossil record. The degree to which variation in fossils of T. carolina can be accommodated by extant variation in the species has been disagreed upon for over eighty years. Using morphometric analyses of the carapace, I address the relationship between modern and fossil T. carolina in terms of sexual dimorphism, geographic and subspecific variation, and allometric variation. Modern T. carolina display weak male-biased sexual size dimorphism. Sexual shape dimorphism cannot be reliably detected in the fossil record. Rather than a four-part subspecific division, patterns of geographic variation are more consistent with clinal variation between various regions in the species distribution. Allometric patterns are qualitatively similar to those documented in other emydid turtles and explain a significant amount of shape variation. When allometric patterns are accounted for, Holocene specimens are not significantly different from modern specimens. In contrast, several geologically older specimens have significantly different carapace shape with no modern analogue. Those large, fossilized specimens represent extinct variation occupying novel portions of morphospace. This study highlights the need for additional documentation of modern osteological variation that can be used to test hypotheses of intraspecific evolution in the fossil record.

Geographical features are the predominant driver of molecular diversification in widely distributed North American whipsnakes

Allopatric divergence following the formation of geographical features has been implicated as a major driver of evolutionary diversification. Widespread species complexes provide opportunities to examine allopatric divergence across varying degrees of isolation in both time and space. In North America, several geographical features may play such a role in diversification, including the Mississippi River, Pecos River, Rocky Mountains, Cochise Filter Barrier, Gulf of California and Isthmus of Tehuantepec. We used thousands of nuclear single nucleotide polymorphisms (SNPs) and mitochondrial DNA from several species of whipsnakes (genera Masticophis and Coluber) distributed across North and Central America to investigate the role that these geographical features have played on lineage divergence. We hypothesize that these features restrict gene flow and separate whipsnakes into diagnosable genomic clusters. We performed genomic clustering and phylogenetic reconstructions at the species and population levels using Bayesian and likelihood analyses and quantified migration levels across geographical features to assess the degree of genetic isolation due to allopatry. Our analyses suggest that (i) major genetic divisions are often consistent with isolation by geographical features, (ii) migration rates between clusters are asymmetrical across major geographical features, and (iii) areas that receive proportionally more migrants possess higher levels of genetic diversity. Collectively, our findings suggest that multiple features of the North American landscape contributed to allopatric divergence in this widely distributed snake group.

Variation of partial transferrin sequences and phylogenetic relationships among hares (Lepus capensis, Lagomorpha) from Tunisia

North African hares are currently included in cape hares, Lepus capensis sensu lato, a taxon that may be considered a superspecies or a complex of closely related species. The existing molecular data, however, are not unequivocal, with mtDNA control region sequences suggesting a separate species status and nuclear loci (allozymes, microsatellites) revealing conspecificity of L. capensis and L. europaeus. Here, we study sequence variation in the intron 6 (468 bp) of the transferrin nuclear gene, of 105 hares with different coat colour from different regions in Tunisia with respect to genetic diversity and differentiation, as well as their phylogenetic status. Forty-six haplotypes (alleles) were revealed and compared phylogenetically to all available TF haplotypes of various Lepus species retrieved from GenBank. Maximum Likelihood, neighbor joining and median joining network analyses concordantly grouped all currently obtained haplotypes together with haplotypes belonging to six different Chinese hare species and the African scrub hare L. saxatilis. Moreover, two Tunisian haploypes were shared with L. capensis, L timidus, L. sinensis, L. yarkandensis, and L. hainanus from China. These results indicated the evolutionary complexity of the genus Lepus with the mixing of nuclear gene haplotypes resulting from introgressive hybridization or/and shared ancestral polymorphism. We report the presence of shared ancestral polymorphism between North African and Chinese hares. This has not been detected earlier in the mtDNA sequences of the same individuals. Genetic diversity of the TF sequences from the Tunisian populations was relatively high compared to other hare populations. However, genetic differentiation and gene flow analyses (AMOVA, F_ST, Nm) indicated little divergence with the absence of geographically meaningful phylogroups and lack of clustering with coat colour types. These results confirm the presence of a single hare species in Tunisia, but a sound inference on its phylogenetic position would require additional nuclear markers and numerous geographically meaningful samples from Africa and Eurasia.

Influence of drainage divides versus arid corridors on genetic structure and demography of a widespread freshwater turtle, Emydura macquarii krefftii, from Australia

The influence of Pleistocene climatic cycles on Southern Hemisphere biotas is not yet well understood. Australia's eastern coastal margin provides an ideal setting for examining the relative influence of landscape development, sea level fluctuation, and cyclic climatic aridity on the evolution of freshwater biodiversity. We examined the impact of climatic oscillations and physical biogeographic barriers on the evolutionary history of the wide-ranging Krefft's river turtle (Emydura macquarii krefftii), using range-wide sampling (649 individuals representing 18 locations across 11 drainages) and analysis of mitochondrial sequences (∼1.3-kb control region and ND4) and nuclear microsatellites (12 polymorphic loci). A range of phylogeographic (haplotype networks, molecular dating), demographic (neutrality tests, mismatch distributions), and population genetic analyses (pairwise F ST, analysis of molecular variance, Bayesian clustering analysis) were implemented to differentiate between competing demographic (local persistence vs. range expansion) and biogeographic (arid corridor vs. drainage divide) scenarios. Genetic data reveal population genetic structure in Krefft's river turtles primarily reflects isolation across drainage divides. Striking north-south regional divergence (2.2% ND4 p-distance; c. 4.73 Ma, 95% higher posterior density (HPD) 2.08-8.16 Ma) was consistent with long-term isolation across a major drainage divide, not an adjacent arid corridor. Ancient divergence among regional lineages implies persistence of northern Krefft's populations despite the recurrent phases of severe local aridity, but with very low contemporary genetic diversity. Stable demography and high levels of genetic diversity are inferred for southern populations, where aridity was less extreme. Range-wide genetic structure in Krefft's river turtles reflects contemporary and historical drainage architecture, although regional differences in the extent of Plio-Pleistocene climatic aridity may be reflected in current levels of genetic diversity.

Venomics of New World pit vipers: genus-wide comparisons of venom proteomes across Agkistrodon

We report a genus-wide comparison of venom proteome variation across New World pit vipers in the genus Agkistrodon. Despite the wide variety of habitats occupied by this genus and that all its taxa feed on diverse species of vertebrates and invertebrate prey, the venom proteomes of copperheads, cottonmouths, and cantils are remarkably similar, both in the type and relative abundance of their different toxin families. The venoms from all the eleven species and subspecies sampled showed relatively similar proteolytic and PLA2 activities. In contrast, quantitative differences were observed in hemorrhagic and myotoxic activities in mice. The highest myotoxic activity was observed with the venoms of A. b. bilineatus, followed by A. p. piscivorus, whereas the venoms of A. c. contortrix and A. p. leucostoma induced the lowest myotoxic activity. The venoms of Agkistrodon bilineatus subspecies showed the highest hemorrhagic activity and A. c. contortrix the lowest. Compositional and toxicological analyses agree with clinical observations of envenomations by Agkistrodon in the USA and Central America. A comparative analysis of Agkistrodon shows that venom divergence tracks phylogeny of this genus to a greater extent than in Sistrurus rattlesnakes, suggesting that the distinct natural histories of Agkistrodon and Sistrurus clades may have played a key role in molding the patterns of evolution of their venom protein genes.

BIOLOGICAL SIGNIFICANCE

A deep understanding of the structural and functional profiles of venoms and of the principles governing the evolution of venomous systems is a goal of venomics. Isolated proteomics analyses have been conducted on venoms from many species of vipers and pit vipers. However, making sense of these large inventories of data requires the integration of this information across multiple species to identify evolutionary and ecological trends. Our genus-wide venomics study provides a comprehensive overview of the toxic arsenal across Agkistrodon and a ground for understanding the natural histories of, and clinical observations of envenomations by, species of this genus.

Evaluating the significance of paleophylogeographic species distribution models in reconstructing quaternary range-shifts of nearctic chelonians

The climatic cycles of the Quaternary, during which global mean annual temperatures have regularly changed by 5-10°C, provide a special opportunity for studying the rate, magnitude, and effects of geographic responses to changing climates. During the Quaternary, high- and mid-latitude species were extirpated from regions that were covered by ice or otherwise became unsuitable, persisting in refugial retreats where the environment was compatible with their tolerances. In this study we combine modern geographic range data, phylogeny, Pleistocene paleoclimatic models, and isotopic records of changes in global mean annual temperature, to produce a temporally continuous model of geographic changes in potential habitat for 59 species of North American turtles over the past 320 Ka (three full glacial-interglacial cycles). These paleophylogeographic models indicate the areas where past climates were compatible with the modern ranges of the species and serve as hypotheses for how their geographic ranges would have changed in response to Quaternary climate cycles. We test these hypotheses against physiological, genetic, taxonomic and fossil evidence, and we then use them to measure the effects of Quaternary climate cycles on species distributions. Patterns of range expansion, contraction, and fragmentation in the models are strongly congruent with (i) phylogeographic differentiation; (ii) morphological variation; (iii) physiological tolerances; and (iv) intraspecific genetic variability. Modern species with significant interspecific differentiation have geographic ranges that strongly fluctuated and repeatedly fragmented throughout the Quaternary. Modern species with low genetic diversity have geographic distributions that were highly variable and at times exceedingly small in the past. Our results reveal the potential for paleophylogeographic models to (i) reconstruct past geographic range modifications, (ii) identify geographic processes that result in genetic bottlenecks; and (iii) predict threats due to anthropogenic climate change in the future.

Phylogeography of North American herbaceous Smilax (Smilacaceae): Combined AFLP and cpDNA data support a northern refugium in the Driftless Area

PREMISE OF THE STUDY

The genetic structure of North American herbaceous Smilax, a clade of closely related understory herbs that has a wide and disjunct distribution, was investigated to test the hypothesis of whether a northern refugium in the upper Midwest may have existed for these plants during the last glaciation.

METHODS

We analyzed 33 populations sampled from California and throughout the eastern United States using AFLP and chloroplast DNA sequence data.

KEY RESULTS

Individuals of S. jamesii from northern California formed a clade sister to the eastern North American species, and they likely diverged from each other during the Pleistocene. Among the eastern species, two lineages were found on opposite sides of the Appalachian Mountains. The populations near the "Driftless Area" contained most of the chlorotypes found in the Midwest, including a unique one endemic to this area, and the AFLP data indicated that one population from this area has the highest frequency-down-weighted-marker value.

CONCLUSIONS

This study, and others that have focused on mammals, amphibians, and woody plants, provides further evidence for the debate over whether the upper Midwest's "Driftless Area" may have hosted a biologically diverse northern glacial refugium. Herbaceous species of Smilax from eastern North American exhibit an Appalachian discontinuity. They appear to have persisted in multiple refugia both east and west of the Appalachians, with several populations persisting in situ through the Last Glacial Maximum (LGM). We suggest that they experienced localized expansion after the LGM, possibly followed by subsequent contact between the Midwest and the East Coastal lineages.

Biogeography and body size shuffling of aquatic salamander communities on a shifting refuge

Freshwater habitats of coastal plains are refugia for many divergent vertebrate lineages, yet these environments are highly vulnerable to sea-level fluctuations, which suggest that resident communities have endured dynamic histories. Using the fossil record and a multi-locus nuclear phylogeny, we examine divergence times, biogeography, body size evolution and patterns of community assembly of aquatic salamanders from North American coastal plains since the Late Cretaceous. At least five salamander families occurred on the extensive Western Interior Coastal Plain (WICP), which existed from the Late Cretaceous through the Eocene. Four of these families subsequently colonized the emergent Southeastern Coastal Plain (SECP) by the Early Oligocene to Late Miocene. Three families ultimately survived and underwent extensive body size evolution in situ on the SECP. This included at least two major size reversals in recent taxa that are convergent with confamilial WICP ancestors. Dynamics of the coastal plain, major lineage extinctions and frequent extreme changes in body size have resulted in significant shuffling of the size structure of aquatic salamander communities on this shifting refuge since the Cretaceous.

Multi-locus phylogeographic and population genetic analysis of Anolis carolinensis: historical demography of a genomic model species

The green anole (Anolis carolinensis) has been widely used as an animal model in physiology and neurobiology but has recently emerged as an important genomic model. The recent sequencing of its genome has shed new light on the evolution of vertebrate genomes and on the process that govern species diversification. Surprisingly, the patterns of genetic diversity within natural populations of this widespread and abundant North American lizard remain relatively unknown. In the present study, we use 10 novel nuclear DNA sequence loci (N = 62 to 152) and one mitochondrial locus (N = 226) to delimit green anole populations and infer their historical demography. We uncovered four evolutionarily distinct and geographically restricted lineages of green anoles using phylogenetics, bayesian clustering, and genetic distance methods. Molecular dating indicates that these lineages last shared a common ancestor ∼2 million years ago. Summary statistics and analysis of the frequency distributions of DNA polymorphisms strongly suggest range-wide expansions in population size. Using Bayesian Skyline Plots, we inferred the timing of population size expansions, which differ across lineages, and found evidence for a relatively recent and rapid westward expansion of green anoles across the Gulf Coastal Plain during the mid-Pleistocene. One surprising result is that the distribution of genetic diversity is not consistent with a latitudinal shift caused by climatic oscillations as is observed for many co-distributed taxa. This suggests that the most recent Pleistocene glacial cycles had a limited impact on the geographic distribution of the green anole at the northern limits of its range.

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.

The genetic crisis of the Mexican Bolson Tortoise (Gopherus flavomarginatus: Testudinidae)

Gopherusflavomarginatus(Testudinidae) is endemic to a series of discontinuous, isolated basins collectively known as the Bolson de Mapimí in the Chihuahuan Desert. Its numbers declined after catastrophic levels of exploitation during the mid-20th century. However currently, the Bolson Tortoise appears to be on a path to recovery owing to intensive, sustained conservation efforts. We sequenced an 842-bp-long fragment of the D-loop from 76 individuals distributed throughout the species’ range. The results revealed only two haplotypes. An AMOVA showed that 95% of the variance occurred among populations, whereas the remaining 5% was explained by genetic differences within populations. Tectonic processes together with ecological transformation during the Pleistocene and Holocene may be responsible for the reduction in this species’ genetic variation. A bottleneck during which a significant percentage of the haplotype diversity was lost would result in genetic homogeneity. Although there is demographic growth, the lack of genetic diversity is indicative of the potential crisis that the Bolson Tortoise is facing, and awareness must be brought to this situation.