Fig. 1

Summary of species delimitation results. Shown is a seven-gene phylogeny estimated in BEAST using the partitioned nuclear sequence, a strict molecular clock, and coalescent prior. To the left of the phylogeny are images of the three morphospecies found in this study and a map of the sites where all samples were collected generated using the ggmap package in R with the “get_googlemap()” function. Columns to the right of the phylogeny are the morphospecies for each sample (where available), and the results of four species delimitation methods: the A11 algorithm of BPP, single-threshold GMYC, bGMYC, and multiple-threshold GMYC. If multiple samples were supported by a species delimitation method to represent the same species, that putative species is represented in a unique color in a rectangle separated from other species. Haplotype names for each sample (or group of samples) are shown; where species delimitation methods recovered multiple haplotypes as the same species, the most common haplotype of the group is given first. Note that haplotypes CATUST09 and TUMIG15 were not included in the BPP analyses and so are not depicted in the column for BPP