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 (