Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus , we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

No abstract

Biodiversity elements with narrow niches and restricted distributions (i.e., ‘short range endemics,’ SREs) are particularly vulnerable to climate change. The New Mexico Ridge-nosed Rattlesnake (Crotalus willardi obscurus, CWO), an SRE listed under the U.S. Endangered Species Act within three sky islands of southwestern North America, is constrained at low elevation by drought and at high elevation by wildfire. We combined long-term recapture and molecular data with demographic and niche modeling to gauge its climate-driven status, distribution, and projected longevity. The largest population (Animas) is numerically constricted (N = 151), with few breeding adults (Nb = 24) and an elevated inbreeding coefficient (ΔF = 0.77; 100 years). Mean home range (0.07km2) is significantly smaller compared to other North American rattlesnakes, and movements are within, not among sky islands. Demographic values, when gauged against those displayed by other endangered/Red-Listed reptiles [e.g., Loggerhead Sea Turtle ( Caretta caretta )], are either comparable or markedly lower. Survival rate differs significantly between genders (female

Glacial-interglacial cycles of the Pleistocene are hypothesized as one of the foremost contributors to biological diversification. This is especially true for cold-adapted montane species, where range shifts have had a pronounced effect on population-level divergence. Gartersnakes of the Thamnophis rufipunctatus species complex are restricted to cold headwater streams in the highlands of the Sierra Madre Occidental and southwestern USA. We used coalescent and multilocus phylogenetic approaches to test whether genetic diversification of this montane-restricted species complex is consistent with two prevailing models of range fluctuation for species affected by Pleistocene climate changes. Our concatenated nuDNA and multilocus species analyses recovered evidence for the persistence of multiple lineages that are restricted geographically, despite a mtDNA signature consistent with either more recent connectivity (and introgression) or recent expansion (and incomplete lineage sorting). Divergence times estimated using a relaxed molecular clock and fossil calibrations fall within the Late Pleistocene, and zero gene flow scenarios among current geographically isolated lineages could not be rejected. These results suggest that increased climate shifts in the Late Pleistocene have driven diversification and current range retraction patterns and that the differences between markers reflect the stochasticity of gene lineages (i.e. ancestral polymorphism) rather than gene flow and introgression. These results have important implications for the conservation of T. rufipunctatus (sensu novo), which is restricted to two drainage systems in the southwestern US and has undergone a recent and dramatic decline.

No abstract