Neogene diversification and taxonomic stability in the snake tribe Lampropeltini (Serpentes: Colubridae)

Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae)

Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.

Comparative analysis of chloroplast genome and new insights into phylogenetic relationships of Ajuga and common adulterants

Introduction

The potential contamination of herbal medicinal products poses a significant concern for consumer health. Given the limited availability of genetic information concerning Ajuga species, it becomes imperative to incorporate supplementary molecular markers to enhance and ensure accurate species identification.

Methods

In this study, the chloroplast (cp) genomes of seven species of the genus Ajuag were sequenced, de novo assembled and characterized.

Results

exhibiting lengths ranging from 150,342 bp to 150,472 bp, encompassing 86 - 88 protein-coding genes (PCGs), 35 - 37 transfer RNA, and eight ribosomal RNA. The repetitive sequences, codon uses, and cp genomes of seven species were highly conserved, and PCGs were the reliable molecular markers for investigating the phylogenetic relationship within the Ajuga genus. Moreover, four mutation hotspot regions (accD-psaI, atpH-atpI, ndhC-trnV(UAC), and ndhF-rpl23) were identified within cp genomes of Ajuga, which could help distinguish A. bracteosa and its contaminants. Based on cp genomes and PCGs, the phylogenetic tree preliminary confirmed the position of Ajuga within the Lamiaceae family. It strongly supported a sister relationship between Subsect. Genevense and Subsect. Biflorae, suggesting the merger of Subsect. Biflorae and Subsect. Genevenses into one group rather than maintaining separate categorizations. Additionally, molecular clock analysis estimated the divergence time of Ajuga to be around 7.78 million years ago.

Discussion

The species authentication, phylogeny, and evolution analyses of the Ajuga species may benefit from the above findings.

Considering admixture when producing draft genomes: an example in North American ratsnakes (Pantherophis alleghaniensis/Pantherophis obsoletus)

The number of reference genomes of snakes lags behind several other vertebrate groups (e.g. birds and mammals). However, in the last two years, a concerted effort by researchers from around the world has produced new genomes of snakes representing members from several new families. Here, we present a high-quality, annotated genome of the central ratsnake (Pantherophis alleghaniensis), a member of the most diverse snake lineage, Colubroidea. Pantherophis alleghaniensis is found in the central part of the Nearctic, east of the Mississippi River. This genome was sequenced using 10X Chromium synthetic long reads and polished using Illumina short reads. The final genome assembly had an N50 of 21.82 Mb and an L50 of 22 scaffolds with a maximum scaffold length of 82.078 Mb. The genome is composed of 49.24% repeat elements dominated by long interspersed elements. We annotated this genome using transcriptome assemblies from 14 tissue types and recovered 28,368 predicted proteins. Finally, we estimated admixture proportions between two species of ratsnakes and discovered that this specimen is an admixed individual containing genomes from the western (Pantherophis obsoletus) and central ratsnakes (P. alleghaniensis). We discuss the importance of considering interspecific admixture in downstream approaches for inferring demography and phylogeny.

A new, rare, small-ranged, and endangered mountain snake of the genus Elaphe from the Southern Levant

The genus Elaphe Fitzinger, 1833 includes 17 species of charismatic, large-sized, non-venomous, Eurasian snakes. In the Western Palearctic, the genus is represented by three species from the Elaphe quatuorlineata group ranging from the Apennine peninsula to Central Asia. The southernmost population of this group is distributed in the mountains of the Southern Levant, with more than 400 km gap to other Elaphe populations. This population has been known to science for only 50 years and is virtually unstudied due to its extreme rarity. We studied these snakes' morphological and genetic variation from the three countries where they are known to occur, i.e., Israel (Hermon, the Israeli-controlled Golan Heights), Lebanon, and Syria. We used nine mitochondrial and nuclear genes, complete mitogenome sequences, and a comprehensive morphological examination including published data, our own field observations, and museum specimens, to study its relationship to other species in the group. The three currently recognized species of the group (E. quatuorlineata, E. sauromates, E. urartica), and the Levant population, form four deeply divergent, strongly supported clades. Three of these clades correspond to the abovementioned species while the Southern Levant clade, which is genetically and morphologically distinct from all named congeners, is described here as a new species, Elaphe druzei sp. nov. The basal divergence of this group is estimated to be the Late Miocene with subsequent radiation from 5.1 to 3.9 Mya. The revealed biogeography of the E. quatuorlineata group supports the importance of the Levant as a major center of endemism and diversity of biota in Eurasia. The new species is large-sized and is one of the rarest snakes in the Western Palearctic. Because of its small mountain distribution range, in an area affected by land use and climate change, the new Elaphe urgently needs strict protection. Despite political issues, we hope this will be based on the cooperation of all countries where the new species occurs.

Phylogenetic relationships and biogeographic range evolution in cat-eyed snakes, Boiga (Serpentes: Colubridae)

The genus Boiga includes 35, primarily arboreal snake species distributed from the Middle East to Australia and many islands in the western Pacific, with particularly high species diversity in South-East Asia. Despite including the iconic mangrove snakes (Boiga dendrophila complex) and the brown tree snake (Boiga irregularis; infamous for avian extinctions on small islands of the Pacific), species-level phylogenetic relationships and the biogeographic history of this ecologically and morphologically distinct clade are poorly understood. In this study, we sequenced mitochondrial and nuclear DNA for 24 Boiga species and used these data to estimate a robust phylogenetic inference, in order to (1) test the hypothesis that Boiga is monophyletic, (2) evaluate the validity of current species-level taxonomy and (3) examine whether geographic range evolution in Boiga is consistent with expectations concerning dispersal and colonization of vertebrates between continents and islands. Our results support the prevailing view that most dispersal events are downstream – from continents to oceanic islands – but we also identify a role for upstream dispersal from oceanic islands to continents. Additionally, the novel phylogeny of Boiga presented here is informative for updating species-level taxonomy within the genus.

Out of the Hengduan Mountains: Molecular phylogeny and historical biogeography of the Asian water snake genus Trimerodytes (Squamata: Colubridae)

The Asian water snake genus Trimerodytes (formerly Sinonatrix) is endemic to East and Southeast Asia. Although several species have been included in various phylogenetic studies previously, the evolution and relationships among members of this genus as a whole remain unexplored. In this study, we report the sequencing two protein-coding mitochondrial gene fragments (MTCYB and ND2) and three nuclear genes (c-mos, NT3, and Rag1), reconstruct interspecific phylogeny, and explore biogeography for the genus Trimerodytes. Both Bayesian inference and maximum likelihood analyses consistently recover the monophyly of Trimerodytes with strong support, with T. yapingi the sister-group to the remaining species. The divergence date and ancestral area estimation suggest that Trimerodytes likely originated in Hengduan Mountains (eastern Tibetan Plateau) in western China at 23.93 Ma (95% HPD: 17.09-31.30), and intraspecific divergence began at about 4.23 Ma (95% HPD: 2.74-6.10). Analyses support the validity of T. yunnanensis.

Phylogeny of caucasian rock lizards (Darevskia) and other true lizards based on mitogenome analysis: Optimisation of the algorithms and gene selection

We generated a phylogeny for Caucasian rock lizards (Darevskia), and included six other families of true lizards (Lacertini), based on complete mitochondrial genome analysis. Next-generation sequencing (NGS) of genomic DNA was used to obtain 16 new mitogenomes of Darevskia. These, along with 35 sequences downloaded from GenBank: genera Darevskia, Zootoca, Podarcis, Phoenicolacerta, Takydromus, Lacerta, and Eremias-were used in the analysis. All four analytical methods (Bayesian Inference, BI; Maximum Likelihood, ML; Maximum Parsimony, MP; and Neighbor-Joining, NJ) showed almost congruent intra-generic topologies for Darevskia and other lizard genera. However, ML and NJ methods on one side, and BI and MP methods on the other harvested conflicting phylogenies. The ML/NJ topology supports earlier published separation of Darevskia into three mitochondrial clades (Murphy, Fu, Macculloch, Darevsky, and Kupinova, 2000), but BI and MP topologies support that the basal branching occurred between D. parvula from the western Lesser Caucasus and the rest of Darevskia. All topologies altered the phylogenetic position of some individual species, including D. daghestanica, D. derjugini, and D. chlorogaster. Reanalysis after excluding four saturated genes from the data set, and excluding genus Eremias gives fully convergent topologies. The most basal branching for true lizards was between Far Eastern Takydromus and the Western Eurasian genera (BI). Comparing phylogenetic performance of individual genes relative to whole mitogenome data, concatenated 16S RNA (the least saturated gene in our analyses) and Cytochrome b genes generate a robust phylogeny that is fully congruent with that based on the complete mitogenome.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

Hidden diversity within the depauperate genera of the snake tribe Lampropeltini (Serpentes, Colubridae)

Accurate representation of lineage diversity through complete taxon sampling is crucial to understanding the evolution of biodiversity, particularly when using molecular phylogenetics to estimate evolutionary relationships. In this interest, taxonomic diversity is often used as a proxy for lineage diversity even though the two concepts are not synonymous. We explore this within the snake tribe Lampropeltini which includes some of the most conspicuous and heavily studied snakes in North America. Both the taxonomy and hypothesized relationships within this tribe have been in flux. The number of species has increased from 23 to 51 over the last thirty years, predominately within three of the nine genera (Lampropeltis, Pantherophis, Pituophis). The remaining six depauperate genera (Arizona, Bogertophis, Cemophora, Pseudelaphe, Rhinocheilus, and Senticolis) have been poorly represented in phylogenetic studies. To estimate evolutionary relationships and determine if the dichotomy in depauperate and speciose genera within Lampropeltini is a function of taxon sampling or truly represents the lineage diversity, we estimated the phylogeny of this group using nuclear and mitochondrial loci in a concatenated and coalescent framework with the largest sampling of the six depauperate genera to date. In addition, we estimated the divergence dates among the genera to assess whether the instability of Lampropeltini phylogenetic relationships is due to an adaptive radiation. While some nodes still remain unresolved, the generic-level relationships we recovered agree with those of a recent next-generation study that used a much larger set of loci for fewer individuals. We also tested two putative species, Arizona pacata and Pseudelaphe phaescens, for the first time phylogenetically and find evidence that they are distinct lineages. Overall, we find that the taxonomic and genetic diversity are not correlated in Lampropeltini and that representing putative diversity in phylogenies will lead to a better estimate of evolutionary histories, especially in groups with complex radiations.

Independent Demographic Responses to Climate Change among Temperate and Tropical Milksnakes (Colubridae: Genus Lampropeltis)

The effects of Late Quaternary climate change have been examined for many temperate New World taxa, but the impact of Pleistocene glacial cycles on Neotropical taxa is less well understood, specifically with respect to changes in population demography. Here, we examine historical demographic trends for six species of milksnake with representatives in both the temperate and tropical Americas to determine if species share responses to climate change as a taxon or by area (i.e., temperate versus tropical environments). Using a multilocus dataset, we test for the demographic signature of population expansion and decline using non-genealogical summary statistics, as well as coalescent-based methods. In addition, we determine whether range sizes are correlated with effective population sizes for milksnakes. Results indicate that there are no identifiable trends with respect to demographic response based on location, and that species responded to changing climates independently, with tropical taxa showing greater instability. There is also no correlation between range size and effective population size, with the largest population size belonging to the species with the smallest geographic distribution. Our study highlights the importance of not generalizing the demographic histories of taxa by region and further illustrates that the New World tropics may not have been a stable refuge during the Pleistocene.

A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes

BACKGROUND

The extant squamates (>9400 known species of lizards and snakes) are one of the most diverse and conspicuous radiations of terrestrial vertebrates, but no studies have attempted to reconstruct a phylogeny for the group with large-scale taxon sampling. Such an estimate is invaluable for comparative evolutionary studies, and to address their classification. Here, we present the first large-scale phylogenetic estimate for Squamata.

RESULTS

The estimated phylogeny contains 4161 species, representing all currently recognized families and subfamilies. The analysis is based on up to 12896 base pairs of sequence data per species (average = 2497 bp) from 12 genes, including seven nuclear loci (BDNF, c-mos, NT3, PDC, R35, RAG-1, and RAG-2), and five mitochondrial genes (12S, 16S, cytochrome b, ND2, and ND4). The tree provides important confirmation for recent estimates of higher-level squamate phylogeny based on molecular data (but with more limited taxon sampling), estimates that are very different from previous morphology-based hypotheses. The tree also includes many relationships that differ from previous molecular estimates and many that differ from traditional taxonomy.

CONCLUSIONS

We present a new large-scale phylogeny of squamate reptiles that should be a valuable resource for future comparative studies. We also present a revised classification of squamates at the family and subfamily level to bring the taxonomy more in line with the new phylogenetic hypothesis. This classification includes new, resurrected, and modified subfamilies within gymnophthalmid and scincid lizards, and boid, colubrid, and lamprophiid snakes.

Genus-level phylogeny of snakes reveals the origins of species richness in Sri Lanka

Snake diversity in the island of Sri Lanka is extremely high, hosting at least 89 inland (i.e., non-marine) snake species, of which at least 49 are endemic. This includes the endemic genera Aspidura, Balanophis, Cercaspis, Haplocercus, and Pseudotyphlops, which are of uncertain phylogenetic affinity. We present phylogenetic evidence from nuclear and mitochondrial loci showing the relationships of 40 snake species from Sri Lanka (22 endemics) to the remaining global snake fauna. To determine the phylogenetic placement of these species, we create a molecular dataset containing 10 genes for all global snake genera, while also sampling all available species for genera with endemic species occurring in Sri Lanka. Our sampling comprises five mitochondrial genes (12S, 16S, cyt-b, ND2, and ND4) and five nuclear genes (BDNF, c-mos, NT3 RAG-1, and RAG-2), for a total of up to 9582bp per taxon. We find that the five endemic genera represent portions of four independent colonizations of Sri Lanka, with Cercaspis nested within Colubrinae, Balanophis in Natricinae, Pseudotyphlops in Uropeltidae, and that Aspidura+Haplocercus represents a distinct, ancient lineage within Natricinae. We synonymize two endemic genera that render other genera paraphyletic (Haplocercus with Aspidura, and Cercaspis with Lycodon), and discover that further endemic radiations may be present on the island, including a new taxon from the blindsnake family Typhlopidae, suggesting a large endemic radiation. Despite its small size relative to other islands such as New Guinea, Borneo, and Madagascar, Sri Lanka has one of the most phylogenetically diverse island snake faunas in the world, and more research is needed to characterize the island's biodiversity, with numerous undescribed species in multiple lineages.

Understanding the formation of ancient intertropical disjunct distributions using Asian and Neotropical hinged-teeth snakes (Sibynophis and Scaphiodontophis: Serpentes: Colubridae)

Numerous taxa show ancient intertropical disjunct distributions. Many can be explained by well-known processes of historical vicariance, such as the breakup of Gondwanaland. Others, such as Asian-Neotropical divergences are not as well understood. To clarify the phylogenetic position and understand biogeographic and temporal origins of the geographically disjunct and morphologically unique genera of hinged-teeth snakes, Scaphiodontophis (n=1) and Sibynophis (n=9; Colubridae), we inferred a time-calibrated phylogeny with additional 107 taxa representing the superfamily Colubroidea using four genes (c-mos, cyt-b, ND2, RAG-1; 3085 bp). We used this tree to estimate ancestral areas for the group. The results show that Scaphiodontophis is sister to Sibynophis, both originated in the late Eocene/Oligocene in Asia and likely dispersed through Beringia to the New World, but unlike other snake groups left no extant species in temperate North America. Current recognition of Scaphiodontophiinae renders Colubrinae paraphyletic, and we resurrect the previously named subfamily Sibynophiinae to encompass both genera and use the tribes Sibynophiini (Sibynophis) and Scaphiodontophiini (Scaphiodontophis) to highlight the geographically distinct areas occupied by these taxa. These results suggest that intercontinental dispersal with extinction in intermediate areas can explain puzzling patterns of ancient intertropical disjunct distributions.

Trait-dependent diversification and the impact of palaeontological data on evolutionary hypothesis testing in New World ratsnakes (tribe Lampropeltini)

For studies investigating trait evolution, there are at least two important questions. First, have traits under consideration influenced cladogenesis and extinction in the group? Second, how do fossil data alter inferences about trait evolution or diversification-rate dynamics? However, relatively few studies have assessed these questions. Here, we use recently developed methods to test for trait-dependent diversification in the New World colubrid snake tribe Lampropeltini. We also integrate data from fossil taxa into phylogenetic estimation of evolutionary parameters using a simple Monte Carlo randomization test. These analyses suggest that ecological conditions in temperate regions are tied to higher rates of cladogenesis, but that body size is not related to diversification in the group. We also find that the inclusion of fossil taxa alters absolute estimates of size and the rate of size evolution, but not the overall pattern of ecomorphological diversification, as well as estimates of evolutionary rates, particularly extinction.

Refugial isolation and divergence in the Narrowheaded Gartersnake species complex (Thamnophis rufipunctatus) as revealed by multilocus DNA sequence data

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.

Mitochondrial DNA sequence variation and phylogeography of the scarlet kingsnake (Lampropeltis elapsoides)

BACKGROUND AND AIMS. With the goal of assessing population structure and geographic distribution of haplotype lineages among Lampropeltis elapsoides, we sequenced the ND4 mitochondrial DNA locus from 96 specimens of this snake across its area of distribution. MATERIALS AND METHODS. We relied heavily on formalin-fixed museum specimens to accomplish this analysis. RESULTS. The sequence alignment consisted of 491 bp of the selected gene, with 28% missing data. A simulation used to assess the effect of missing data on population genetic and phylogenetic resolution indicated increased character conflict, but with minimal loss of phylogenetic structure. CONCLUSION. This limited dataset suggests that L. elapsoides constitutes a largely unstructured population, with both widespread haplotypes and large number of private haplotypes, a moderate level of nucleotide diversity, and a low, but significant, degree of north-south population differentiation. Haplotype structure and frequency, nucleotide frequency, and values for Tajima's D and Fu's F(S) indicate a recent range or population expansion following a historic bottleneck.

The phylogeny of advanced snakes (Colubroidea), with discovery of a new subfamily and comparison of support methods for likelihood trees

The superfamily Colubroidea (> 2500 species) includes the majority of snake species and is one of the most conspicuous and well-known radiations of terrestrial vertebrates. However, many aspects of the phylogeny of the group remain contentious, and dozens of genera have yet to be included in molecular phylogenetic analyses. We present a new, large-scale, likelihood-based phylogeny for the colubroids, including 761 species sampled for up to five genes: cytochrome b (93% of 761 species sampled), ND4 (69%), ND2 (28%), c-mos (54%), and RAG-1 (13%), totaling up to 5814bp per species. We also compare likelihood bootstrapping and a recently proposed ultra-fast measure of branch support (Shimodaira-Hasegawa-like [SHL] approximate likelihood ratio), and find that the SHL test shows strong support for several clades that were weakly-supported by bootstrapping in this or previous analyses (e.g., Dipsadinae, Lamprophiidae). We find that SHL values are positively related to branch lengths, but show stronger support for shorter branches than bootstrapping. Despite extensive missing data for many taxa (mean=67% per species), neither bootstrap nor SHL support values for terminal species are related to their incompleteness, and that most highly incomplete taxa are placed in the expected families from previous taxonomy, typically with very strong support. The phylogeny indicates that the Neotropical colubrine genus Scaphiodontophis represents an unexpectedly ancient lineage within Colubridae. We present a revised higher-level classification of Colubroidea, which includes a new subfamily for Scaphiodontophis (Scaphiodontophiinae). Our study provides the most comprehensive phylogeny of Colubroidea to date, and suggests that SHL values may provide a useful complement to bootstrapping for estimating support on likelihood-based trees.

Arboreal habitat structure affects route choice by rat snakes

In arboreal habitats gaps between branches and branch structure profoundly affect the ability of animals to move; hence, an ability to perceive such attributes could facilitate choosing routes that enhance the speed and ease of locomotion. Although many snakes are arboreal, no previous study has determined whether they can perceive structural variation of branches that is mechanically relevant to their locomotion. We tested whether the gap distance, location, and attributes of two destination perches on the far side of a crossable gap affected the route travelled by North American rat snakes (Pantherophis), which are proficient climbers. Snakes usually chose routes with shorter gaps. Within a horizontal plane, the snakes usually went straight rather than crossing an equal distance gap with a 90° turn, which was consistent with our finding that crossing a straight gap was easier. However, decreasing the distance of the gap with a 90° turn eliminated the preference for going straight. Additional factors, such as the width of the landing surface and the complexity of branching of the destination perches, resulted in non-random route choice. Thus, many of the observed biases in route choice suggested abilities to perceive structural variation and select routes that are mechanically beneficial.

High-model abundance may permit the gradual evolution of Batesian mimicry: an experimental test

In Batesian mimicry, a harmless species (the 'mimic') resembles a dangerous species (the 'model') and is thus protected from predators. It is often assumed that the mimetic phenotype evolves from a cryptic phenotype, but it is unclear how a population can transition through intermediate phenotypes; such intermediates may receive neither the benefits of crypsis nor mimicry. Here, we ask if selection against intermediates weakens with increasing model abundance. We also ask if mimicry has evolved from cryptic phenotypes in a mimetic clade. We first present an ancestral character-state reconstruction showing that mimicry of a coral snake (Micrurus fulvius) by the scarlet kingsnake (Lampropeltis elapsoides) evolved from a cryptic phenotype. We then evaluate predation rates on intermediate phenotypes relative to cryptic and mimetic phenotypes under conditions of both high- and low-model abundances. Our results indicate that where coral snakes are rare, intermediate phenotypes are attacked more often than cryptic and mimetic phenotypes, indicating the presence of an adaptive valley. However, where coral snakes are abundant, intermediate phenotypes are not attacked more frequently, resulting in an adaptive landscape without a valley. Thus, high-model abundance may facilitate the evolution of Batesian mimicry.