Conclusions
Our results substantiate species-level genomic divergences betweenS. atlantis , S. hesperis , and S. nausicaa , and analyses using ENMs and landscape resistance surfaces have enabled us to attribute the maintenance of genomic integrity in each of these lineages to a likely mechanism - ecologically-based divergent selection. In addition to our results, we suggest that there is great opportunity to build on this approach in subsequent systematic work. In particular, the incorporation of ecological modelling in genetic species delimitations may be used to reciprocally inform Bayesian priors for phylogeographic and multispecies coalescent-based modelling. Such an integration of genetic and ecological analyses should be particularly useful for species delimitations that have conservation impacts or high visibility in citizen science, and could be readily integrated into a conservation framework that links genomic divergence to ecological distinctiveness, informing prioritization of habitat and biodiversity management. This multidisciplinary method of species delimitation is thus positioned to contribute to stable taxonomies, delimit meaningful biodiversity units for conservation, and characterize extrinsic factors that influence patterns of lineage diversification in a broad range of taxa.