“…Recent empirical and theoretical studies have shown the limitations of Bayesian species delimitation approaches based on the multi-species coalescent (MSC) model (Rannala & Yang, 2013;Yang & Rannala, 2010), pointing out that such methods are not able to statistically distinguish genetic structure due to population isolation from true species boundaries and tend to over-split genetically differentiated populations rather than capturing species divergence (Huang, 2018;Leaché et al, 2018;Leaché, Zhu, Rannala, & Yang, 2019;Sukumaran & Knowles, 2017). The problem of species overestimation can be particularly exacerbated when using vast genome-wide data due to the high power of a large number of markers (hundreds to tens of thousands) to detect fine-grain population structure (Hime et al, 2016;Noguerales et al, 2018;Sukumaran & Knowles, 2017). In this line, recent empirical studies have found that geographically continuous intraspecific populations presenting genotypic differences simply resulted from isolation by distance can be delimited as distinct species when using hundreds of loci (Huang, 2018; Pyron, Hsieh, Lemmon, Lemmon, & Hendry, 2016).…”