Integrative Species Delimitation in a Speciation Continuum: Phylogenomics, Cryptic Diversity, Diversification and Historical Biogeography of Sinocyclocheilus Cavefish

The transition from structured populations to distinct species often unfolds along a speciation continuum. However, empirically dissecting this continuum poses challenges for species delimitation, particularly in rapid radiations marked by recent divergence, incomplete lineage sorting, and gene flow. The species-rich Sinocyclocheilus cavefish radiation of Southwest China, an emerging evolutionary multi-species model system, shows a notable mismatch between its high morphological diversity and the limited divergence observed in commonly used mitochondrial DNA markers. This pattern suggests that true species richness may be underestimated. Yet, comprehensive genome-wide approaches to resolve species boundaries in this group are still lacking. To address this, we combine phylogenomics, coalescent-based species delimitation, population genetics, and historical biogeography, using both genome-wide RAD-seq and multi-locus Sanger data (including nuclear and mitochondrial DNA). Our phylogenomic analyses resolve major clades and reveal substantial cryptic diversity, well beyond what was detected by earlier markers. Species delimitation based on RAD-seq genomic data identifies several cryptic evolutionary lineages. We recover a dynamic divergence history shaped by isolation and episodic connectivity. Biogeographic reconstruction supports a mid-Miocene origin following a major vicariance event, with subsequent founder-event dispersals into subterranean habitats. This long history of fragmentation is further complicated by reticulation and ancestral polymorphism, and is reflected in present-day patterns of restricted gene flow across river valleys. These results highlight the utility of integrative genomic approaches in resolving species boundaries and uncovering the evolutionary processes that underlie high diversity in large and complex radiations.