Abstract:We present a set of tests to examine processes of speciation that will help delimit taxa given deep phylogeographic history and interaction within hybrid zones. Outcomes from distinct phylogeographic pairs of snakes in North America all have the signature of divergence before or near the mid‐Pleistocene often with low migration. For these pairs we demonstrate the presence of stable hybrid zones of varying size with a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between very distinct ecoregions. The outcomes of speciation across organismal dimensions (e.g., ecological, genetic, phenotypic) are often assessed using phylogeographic methods. At one extreme, reproductively isolated lineages represent easily delimitable species differing in many or all dimensions, and at the other, geographically distinct genetic segments introgress across broad environmental gradients with limited phenotypic disparity. In the ambiguous gray zone of speciation, where lineages are genetically delimitable but still interacting ecologically, it is expected that these lineages represent species in the context of ontology and the evolutionary species concept when they are maintained over time with geographically well‐defined hybrid zones, particularly at the intersection of distinct environments. As a result, genetic structure is correlated with environmental differences and not space alone, and a subset of genes fail to introgress across these zones as underlying genomic differences accumulate. We present a set of tests that synthesize species delimitation with the speciation process. We can thereby assess historical demographics and diversification processes while understanding how lineages are maintained through space and time by exploring spatial and genome clines, genotype‐environment interactions, and genome scans for selected loci. Employing these tests in eight lineage‐pairs of snakes in North America, we show that six pairs represent 12 "good" species and that two pairs represent local adaptation and regional population structure. The distinct species pairs all have the signature of divergence before or near the mid‐Pleistocene, often with low migration, stable hybrid zones of varying size, and a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between distinct ecoregions. Locally adapted populations are younger, exhibit higher migration, and less ecological differentiation. Our results demonstrate that interacting lineages can be delimited using phylogeographic and population genetic methods that properly integrate spatial, temporal, and environmental data.

Weak Relationship Between Divergence and Selection in a North American Toad Hybrid Zone

Selection leads to remarkable variability in the outcomes of hybridization across replicate hybrid zones

The maintenance of reproductive isolation in a mosaic hybrid zone between the fire-bellied toads Bombina bombina and B. variegata

Reproductive isolation in a three‐way contact zone

Species Persistence with Hybridization in Toad-Headed Lizards Driven by Divergent Selection and Low Recombination

Testing assertions of widespread introgressive hybridization in a clade of neotropical toads with low mate selectivity (Rhinella granulosa species group)

Multiple forms of selection shape reproductive isolation in a primate hybrid zone

Allopatric speciation in cattails: Genomic landscapes of divergence across spp. suggest balancing selection, introgressions, and the absence of adaptive divergence

Understanding species limits through the formation of phylogeographic lineages

Genomic regions of current low hybridisation mark long-term barriers to gene flow in scarce swallowtail butterflies

Finding Hybrid Incompatibilities Using Genome Sequences from Hybrid Populations

Hybrid zone dynamics under weak Haldane's rule

Differential Introgression among Loci Across a Hybrid Zone of the Intermediate Horseshoe Bat (rhinolophus Affinis).

Maintenance of a narrow hybrid zone between native and introduced red foxes (Vulpes vulpes) despite conspecificity and high dispersal capabilities

The genetic basis of developmental abnormalities in interpopulation hybrids of the moss Ceratodon purpureus

Coevolutionary Interactions between Sexual and Habitat Isolation during Reinforcement

Genetic basis of reproductive isolation in Torrey pine (Pinus torreyana Parry): insights from hybridization and adaptation

Geographic Variation in the Structure of Oak Hybrid Zones Provides Insights into the Dynamics of Speciation

High sexual display trait diversity without measured genetic divergence in a montane hybrid zone involving young species ( Habronattus americanus subgroup, Araneae: Salticidae)

Strong signatures of selection on genes underlying core reinforcement mechanisms in speciating desert tortoises

Untangling the contribution of adaptive versus non-adaptive processes in the evolution of reproductive isolation between Coenonympha butterflies.