Species delineation
The BPP analyses conducted using the 3RAD data set significantly supported the species model that considered five evolutionary lineages (PP > 0.95; Table S3), the four putative species proposed by Aguilar Velasco et al. (2016), E. ruidum spp. 1-4, and the specimens from Pinotepa and the state of Guerrero that these authors regarded as E. ruidum sp. 2x3 and as part of E. ruidum sp. 3, respectively. In contrast, when we assigned the samples according to their geographic proximity, our results widely varied (Table S4).
The Bayes factor comparisons of the two tested species delimitation models using the phased UCE data consistently favored a five species model, with E. ruidum spp. 1, 2, 3 and 4, as well as the specimens from Pinotepa and Guerrero, Mexico, each representing separate species (Table S4).