Testing hypotheses of diversification in Panamanian frogs and freshwater fishes using hierarchical approximate Bayesian computation with model averaging

Abstract

Most Neotropical frog and freshwater fish species sampled to date show phylogeographic breaks along the Pacific coast of the Isthmus of Panama, with lineages in Costa Rica and western Panama isolated from central Panama. We examine temporal patterns of diversification of taxa across this ‘western Panama isthmus’ (WPI) break to test hypotheses about the origin of species geographical distributions and genetic structuring in this region. We tested for synchronous diversification of four codistributed frog taxon-pairs and three fish taxon-pairs sharing the WPI break using hierarchical approximate Bayesian computation with model averaging based on mitochondrial DNA sequences. We also estimated lineage divergence times using full-Bayesian models. Several of our results supported synchronous divergences within the frog and freshwater fish assemblages; however, Bayes factor support was equivocal for or against synchronous or asynchronous diversification. Nevertheless, we infer that frog populations were likely isolated by one or multiple Pliocene–Pleistocene events more recently than predicted by previous models, while fish genetic diversity was structured by Pleistocene events. By integrating our results with external information from geology and elevational sea level modeling, we discuss the implications of our findings for understanding the biogeographical scenario of the diversification of Panamanian frogs and fishes. Consistent with the ‘Bermingham/Martin model’ (Mol. Ecol. 1998, 7: 499-517), we conclude that the regional fish assemblage was fractured by processes shaping isthmian landscapes during the Pleistocene glaciations, including drainage basin isolation during lowered sea levels.