Combined phylogenetic tree based on COI gene fragment reconstructed using maximum likelihood (ML) on RAxML and Bayesian inference (BI) on MrBayes. Support at each node denotes the ML bootstrap value (BV) and BI posterior probability (PP). Nodes highly supported by at least one method (BV ≥ 70 or PP ≥ 0.95) are highlighted by a black dot. Branch lengths were scaled concerning the number of substitutions per site. Each species corresponding to monophyletic clades are represented with different branch colors, species from Japan are highlighted with a colored background: red = H. mogera, azure = H. alba sp. nov., yellow = N. occulta sp. nov., violet = N. insulana sp. nov., light green = N. silvicola sp. nov., blue = N. brevipes, orange = N. terrestris, dark green = N. okinawaensis. The tree is rooted using the species Wraios longiembolus. The right lines denote the results of the species delimitation analysis based on morphological (MO) and molecular (ASAP, mPTP, and bPTP) data.