Species distribution models predict genetic isolation of <i>Hetaerina vulnerata</i> Hagen in Selys, 1853 (Odonata, Calopterygidae)

Understanding how past and current environmental conditions shape the demographic and genetic distributions of organisms facilitates our predictions of how future environmental patterns may affect populations. The Canyon Rubyspot damselfly (Odonata: Zygoptera: <i>Hetaerina vulnerata</i>) is an insect with a range distribution from Colombia to the arid southwestern United States, where it inhabits shaded mountain streams in the arid southwestern United States. Past spatial fragmentation of habitat and limited dispersal capacity of <i>H. vulnerata</i> may cause population isolation and genetic differentiation, and projected climate change may exacerbate isolation by further restricting the species' distribution. We constructed species distribution models (SDMs) based on occurrences of <i>H. vulnerata</i> and environmental variables characterizing the species' niche. We inferred seven current potential population clusters isolated by unsuitable habitat. Paleoclimate models indicated habitat contiguity in past conditions; projected models indicated some habitat fragmentation in future scenarios. Seventy-eight <i>H. vulnerata</i> individuals from six of the current clusters were sequenced via ddRADseq and processed with Stacks. Principal components and phylogeographic analyses resolved three subpopulations; <i>Structure</i> resolved four subpopulations. <i>F</i> _ST values were low (<0.05) for nearby populations and >0.15 for populations separated by expanses of unsuitable habitat. Isolation by distance was an existing but weak factor in determining genomic structure; isolation by environment and the intervening landscape explained a significant proportion of genetic distance. <i>Hetaerina vulnerata</i> populations were shown to be isolated by a lack of tree canopy coverage, an important habitat predictor for oviposition and territoriality. Thus, <i>H. vulnerata</i> populations are likely separated and are genetically isolated. Integrating SDMs with landscape genetics allowed us to identify populations separated by distance and unsuitable habitat, explaining population genetic patterns and probable fates for populations under future climate scenarios.