Integrative species delimitation and biogeography of the Rhinella margaritifera species group (Amphibia, Anura, Bufonidae) suggest an intense diversification throughout Amazonia during the last 10 million years

The accumulation of studies delimiting species in Amazonia has not only shed light on the patterns of its outstanding species richness but also allowed a better understanding of the processes of diversification within this immense region. Nevertheless, vast knowledge gaps remain even for prominent anuran species complexes, such as the Rhinella margaritifera species group. This clade of toads comprises 23 valid species-level taxa, mainly distributed in Amazonia but also in South America’s Dry Diagonal and Atlantic and trans-Andean rainforests. Species boundaries and taxonomy in this group are notoriously complex, with studies suggesting the existence of several unnamed species. Available phylogenetic information suggests an Andean-western Amazonian origin of the group with subsequent diversification within Amazonian lowlands during the last 10 Myr and secondary dispersals into other Neotropical regions. To further test this biogeographic scenario and improve knowledge on species diversity, we used an unprecedentedly large mtDNA sampling (>800 16S sequences) across the clade’s distribution and comprising all but one described species. We delimited 54 Molecular Operational Taxonomic Units, which we tested further based on patterns of variation of a nuclear locus and acoustic and morphological data. This approach confirmed the existence of at least 25 candidate species, 19 of which correspond to currently recognized taxa whereas 30 remained ‘unconfirmed’. Our results clarify the taxonomic status of some species but also suggest multiple introgression events that blur some mtDNA-based species boundaries. Lastly, to provide a temporal framework for the clade’s diversification, we generated a time-calibrated phylogenetic tree based on a mitogenomic matrix, which confirmed a Miocene (∼9 Ma) western Amazonian origin and six major clades in the group, each having initially diversified in different regions within Amazonia. Most of these clades have later dispersed throughout Amazonia during the establishment of the modern Amazonian hydrographic system, i.e., in the last 6 Myr.