The genetic origins of species boundaries at subtropical and temperate ecoregions in the North American racers (Coluber constrictor)

Phylogeographically structured lineages are a common outcome of range-wide population genetic studies. In the southeastern United States, disconnection between populations found at the intersection of the southeastern coastal plains of peninsular Florida and the southeastern plains of the adjacent continent is readily apparent among many plants and animals. However, the timing and maintenance of species boundaries between these distinctly different subtropical and temperate regions remains unknown for all organisms studied there. Using genome-scale data, we examine the timing of origins, gene flow, and the movement of genes under selection in unique ecoregions within the North American racers (Coluber constrictor). Isolation-migration models along with tests of genome-wide selection, locus-environment associations, and spatial and genomic clines demonstrate that two unrecognized species are present and are in contact at the boundary of these two ecoregions. We show that selection at several loci associated with unique environments have maintained species boundaries despite constant levels of gene flow between these lineages over thousands of generations. This research provides a new avenue of research to examine speciation processes in poorly studied biodiversity hotspots.