Extensive Admixture Among Karst‐Obligate Salamanders Reveals Evidence of Recent Divergence and Gene Exchange Through Aquifers

Karst ecosystems often contain extraordinary biodiversity, but the complex underground aquifers of karst regions present challenges for assessing and conserving stygobiont diversity and investigating their evolutionary history. We examined the karst-obligate salamanders of the <i>Eurycea neotenes</i> species complex in the Edwards Plateau region of central Texas using population genomics data to address questions about population connectivity and the potential for gene exchange within the underlying aquifer system. The <i>E. neotenes</i> species complex has historically been divided into three nominal species, but their status, and spatial extent of species ranges, have remained uncertain. We discovered evidence of extensive admixture among species within the complex and with adjacent lineages. We observed relatively low levels of differentiation among all sampling localities which supports the hypothesis of recent divergence. Nominal taxonomy, aquifer region and geography each accounted for a modest amount of the overall population genomic variation; however, these predictors were largely confounded and difficult to disentangle. Importantly, current taxonomy of the three nominal species does not reflect the admixture apparent in clustering analyses. Inference of migration events revealed a complex pattern of gene exchange, suggesting that <i>Eurycea</i> salamanders have a dynamic history of dispersal through the aquifer system. These results highlight the need for greater understanding of how stygobiont populations are connected via dispersal and gene exchange through karst aquifers. These results also highlight the applicability of population genomics data as a powerful lever for investigating connectivity among populations in systems where direct detection of dispersal paths is difficult, as in underground, aquatic systems.