To split or not to split? Multilocus phylogeny and molecular species delimitation of southeast Asian toads (family: Bufonidae)

BACKGROUND: Recent studies have demonstrated that Bayesian species delimitation based on the multispecies coalescent model can produce inaccurate results by misinterpreting population splits as species divergences. An approach based on the genealogical divergence index (gdi) was shown to be a viable alternative, especially for delimiting allopatric populations where gene flow is low. We implemented these analyses to assess species boundaries in Southeast Asian toads, a group that is understudied and characterized by numerous unresolved species complexes. RESULTS: Multilocus phylogenetic analyses showed that deep evolutionary relationships including the genera Sigalegalephrynus, Ghatophryne, Parapelophryne, Leptophryne, Pseudobufo, Rentapia, and Phrynoides remain unresolved. Comparison of genetic divergences revealed that intraspecific divergences among allopatric populations of Pelophyrne signata (Borneo vs. Peninsular Malaysia), Ingerophrynus parvus (Peninsular Malaysia vs. Myanmar), and Leptophryne borbonica (Peninsular Malaysia, Java, Borneo, and Sumatra) are consistent with interspecific divergences of other Southeast Asian bufonid taxa. Conversely, interspecific divergences between Pelophryne guentheri/P. api, Ansonia latiffi/A. leptopus, and I. gollum/I. divergens were low (< 3%) and consistent with intraspecific divergences of other closely related taxa. The BPP analysis produced variable results depending on prior settings and priors estimated from empirical data produced the best results that were also congruent with the gdi analysis. CONCLUSIONS: This study showed that the evolutionary history of Southeast Asian toads is difficult to resolve and numerous relationships remain ambiguous. Although some results from the species delimitation analyses were inconclusive, they were nevertheless efficacious at identifying potential new species and taxonomic incompatibilities for future in-depth investigation. We also demonstrated the sensitivity of BPP to different priors and that careful selection priors based on empirical data can greatly improve the analysis. Finally, the gdi can be a robust tool to complement other species delimitation methods. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1186/s12862-019-1422-3) contains supplementary material, which is available to authorized users.