Species limits and diversification of the Dendropsophus rubicundulus subgroup (Anura, Hylidae) in Neotropical savannas

Understanding the processes that generate and maintain biodiversity at and below the species level is a central goal of evolutionary biology. Here we explore the spatial and temporal drivers of diversification of the treefrog subgroup Dendropsophus rubicundulus, a subgroup of the D. microcephalus species group, over periods of pronounced geological and climatic changes in the Neotropical savannas that they inhabit. This subgroup currently comprises 11 recognized species distributed across the Brazilian and Bolivian savannas, but the taxonomy has been in a state of flux, necessitating reexamination. Using newly generated single nucleotide polymorphism (SNP) data from restriction-site associated DNA sequencing (RADseq) and mitochondrial 16S sequence data for-150 specimens, we inferred phylogenetic relationships, tested species limits using a model-based approach, and estimated divergence times to gain insights into the geographic and climatic events that affected the diversification of this subgroup. Our results recognized at least nine species: D. anataliasiasi, D. araguaya, D. cerradensis, D. elianeae, D. jimi, D. rubicundulus, D. tritaeniatus, D. rozenmani, and D. sanborni. Although we did not collect SNP data for the latter two species, they are likely distinct based on mitochondrial data. In addition, we found genetic structure within the widespread species D. rubicundulus, which comprises three allopatric lineages connected by gene flow upon secondary contact. We also found evidence of population structure and perhaps undescribed diversity in D. elianeae, which warrants further study. The D. rubicundulus subgroup is estimated to have originated in the Late Miocene (-5.45 million years ago), with diversification continuing through the Pliocene and Early Pleistocene, followed by the most recent divergence of D. rubicundulus lineages in the Middle Pleistocene. The epeirogenic uplift followed by erosion and denudation of the central Brazilian plateau throughout the Pliocene and Pleistocene, in combination with the increasing frequency and amplitude of climatic fluctuations during the Pleistocene, was important for generating and structuring diversity at or below the species level in the D. rubicundulus subgroup.