Species delimitation with gene flow: a methodological comparison and population genomics approach to elucidate cryptic species boundaries in Malaysian Torrent Frogs

Accurately delimiting species boundaries is a non-trivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly-used species delimitation methods (SDMs) and a population genomics approach based on genome-wide single nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genome-wide SNPs to identify putative species boundaries by implementing non-coalescent and coalescent-based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2, ABBA-BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean Fst = 0.9) to differentiating populations involvin...