Geographic isolation, genetic divergence, and ecological non-exchangeability define ESUs in a threatened sky-island rattlesnake

We utilized nine microsatellite loci and 865 bases from two mtDNA genes to estimate demographic parameters and visualize historic/contemporary connectivity among populations of a sky-island rattlesnake (New Mexico ridge-nosed rattlesnake, Crotalus willardi obscurus). This taxon is listed as threatened under the United States Endangered Species Act (ESA) and is distributed patchily within three borderland mountain ranges [Animas (ANM), Peloncillo (PEL), Sierra San Luis (SSL)] of southeastern Arizona, southwestern New Mexico and north-central México. Molecular data support a hypothesis of northward range expansion from México, with subsequent isolation on sky-islands through vicariant desertification that transformed intervening wooded valleys and low elevation passes into inhospitable habitat. Historic and recent movements have been within rather than among mountains. Although populations are linked via ancestral polymorphism, they do not represent a single mtDNA gene pool. All are genetically bottlenecked, with the Peloncillo reflecting deepest/sharpest declines and fewest captures per unit effort. Most recent population declines occurred ∼4 kybp (thousands of years before present). Thus, population reductions are historic and environmental rather than contemporary and anthropogenic. Our data demonstrate that PEL is ecologically non-exchangable with other sky-island populations, and thus comprises one ‘evolutionary significant unit’ (ESU), while SSL and ANM comprise ‘management units’ (MUs) within a second ESU. All three meet legal criteria for recognition as ‘distinct population segments’ (DPS) under the ESA.