Rapid speciation in Harlequin Toads ( Anura : Bufonidae ) endemic to the Sierra Nevada de Santa Marta

The Neotropical genus Atelopus has experienced a drastic population decline in recent decades. Despite this, a knowledge gap remains regarding the conservation genetic status, phylogeography, and demographic history of most of its species, especially those endemic to areas with limited access like the Sierra Nevada de Santa Marta in Colombia. In this genomic study, we inferred phylogenetic relationships, demographic history, and gene flow among four endemic Atelopus morphospecies in the Sierra Nevada de Santa Marta-SNSM. Additionally, we compared the effective population size (N _e ) estimates with the available population census data. NextRAD was used to obtain genomic data from 95 individuals collected at five sites in the SNSM and two in the Colombian Pacific (outgroup). The morphospecies recently diverged in a scenario without gene flow and were recovered as monophyletic. Their phylogenetic relationships were discordant, which is attributed to the presence of incomplete lineage sorting-ILS, which would also explain their shared ancestry among them. The lack of gene flow as well as the recent divergence times given by demography suggests a recent and rapid speciation. However, the reproductive isolation mechanisms that promote or maintain the species boundaries in this group remain unknown and require further investigation. We suggest that this process may have been influenced by the complex topography of the SNSM, traits such as high philopatry, low dispersal ability, and behavioral factors such as habitat preference or to factors related to genetic architecture that influenced the rapid formation of reproductive barriers among populations. Additionally, a pattern of population decline was observed around 200.000 years ago, with recent increases in three morphospecies. Despite the reduction in effective population size, no signs of inbreeding were detected. However, for A. laetissimus , the only species surveyed, the estimated value of N _e and its implications should be interpreted with caution. Ultimately, our findings reveal an evolutionary history shaped by a burst of diversification and abrupt reproductive isolation, highlighting how the resulting endemism and restricted genetic connectivity shape the unique evolutionary trajectory and vulnerability of this threatened montane species.