Batesian mimicry, in which harmless species (mimics) deter predators by deceitfully imitating the warning signals of noxious species (models), generates striking cases of phenotypic convergence that are classic examples of evolution by natural selection. However, mimicry of venomous coral snakes has remained controversial because of unresolved conflict between the predictions of mimicry theory and empirical patterns in the distribution and abundance of snakes. Here we integrate distributional, phenotypic and phylogenetic data across all New World snake species to demonstrate that shifts to mimetic coloration in nonvenomous snakes are highly correlated with coral snakes in both space and time, providing overwhelming support for Batesian mimicry. We also find that bidirectional transitions between mimetic and cryptic coloration are unexpectedly frequent over both long- and short-time scales, challenging traditional views of mimicry as a stable evolutionary ‘end point' and suggesting that insect and snake mimicry may have different evolutionary dynamics.

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community‐wide divergence should be explored before assuming the role of biogeographic barriers.

Color polymorphism in natural populations presents an ideal opportunity to study the evolutionary drivers of phenotypic diversity. Systems with striking spatial, temporal, and qualitative variation in color can be leveraged to study the mechanisms promoting the distribution of different types of variation in nature. We used the highly polymorphic ground snake (Sonora semiannulata), a putative coral snake mimic with both cryptic and conspicuous morphs, to compare patterns of neutral genetic variation and variation over space and time in color polymorphism to investigate the mechanistic drivers of phenotypic variation across scales. We found that strong selection promotes color polymorphism across spatial and temporal scales, with morph frequencies differing markedly between juvenile and adult age classes within a single population, oscillating over time within multiple populations, and varying drastically over the landscape despite minimal population genetic structure. However, we found no evidence that conspicuousness of morphs was related to which color pattern was favored by selection or to any geographic factors, including sympatry with coral snakes. We suggest that complex patterns of phenotypic variation in polymorphic systems may be a fundamental outcome of the conspicuousness of morphs and that explicit tests of temporal and geographic variation are critical to the interpretation of conspicuousness and mimicry.