Speciation with gene flow in a narrow endemic West Virginia cave salamander (Gyrinophilus subterraneus)

Due to their limited geographic distributions and specialized ecologies, cave species are often highly endemic and can be especially vulnerable to habitat degradation within and surrounding the cave systems they inhabit. We investigated the evolutionary history of the West Virginia Spring Salamander (Gyrinophilus subterraneus), estimated the population trend from historic and current survey data, and assessed the current potential for water quality threats to the cave habitat. Our genomic data (mtDNA sequence and ddRADseq-derived SNPs) reveal two, distinct evolutionary lineages within General Davis Cave corresponding to G. subterraneus and its widely distributed sister species, Gyrinophilus porphyriticus, that are also differentiable based on morphological traits. Genomic models of evolutionary history strongly support asymmetric and continuous gene flow between the two lineages, and hybrid classification analyses identify only parental and first generation cross (F1) progeny. Collectively, these results point to a rare case of sympatric speciation occurring within the cave, leading to strong support for continuing to recognize G. subterraneus as a distinct and unique species. Due to its specialized habitat requirements, the complete distribution of G. subterraneus is unresolved, but using survey data in its type locality (and currently the only known occupied site), we find that the population within General Davis Cave has possibly declined over the last 45 years. Finally, our measures of cave and surface stream water quality did not reveal evidence of water quality impairment and provide important baselines for future monitoring. In addition, our unexpected finding of a hybrid zone and partial reproductive isolation between G. subterraneus and G. porphyriticus warrants further attention to better understand the evolutionary and conservation implications of occasional hybridization between the species.