Fig. 1 The automatic partition results of dataset with ABGD method

Fig. 2 Lineages-through-time plot based on the time calibrated tree obtained from all 121 haplotypes. The sharp increase in branching rate, corresponding to the transition from interspecies to intraspecies branching events.

Fig. 3 Species richness estimation results with different methods. Pink colour represents ABGD-based method. Green colour represents GMYC-based method. Blue colour represents morphology-based method. Black error bar indicates confidence interval.

Fig. 4 Rank-abundance curves and fittings of five models. a-1)Species rank-abundance curves based on morphological method, a-2) Species rank-abundance curves based on ABGD method, a-3) Species rank-abundance curves based on GMYC method; b-1) Species abundance models based on morphological method, b-2) Species abundance models based on ABGD method, b-3) Species abundance models based on GMYC method. Black points represent the real abundance data. Black curve represents broken stick model. Red curve represents preemption model. Green curve represents log normal model. Dark blue curve represents Zipf model. Light blue curve represents Zipf-Mandelbrot model.