2.4.1 Phylogenetic tree
Neighbor-Joining (NJ) method in MEGA 7.0 software was used for tree construction (Satiou et al., 1987) with Kimura-2-Parameter (K-2-P) genetic distances (Kimura et al., 1980). Meanwhile, to evaluate the accuracy of DNA barcoding for classifying species as well as the evolutionary relationship of species, we downloaded the sequencing data of relevant species from BOLD, including those from Germany, United States, Mexico, Canada, Norway, Italy, Finland, Japan, Iran, etc. (Raupach et al., 2014; Stein et al., 2013) (Part I, Supporting information). A total of 1000 replicates were utilized to detect the reliability of each branch of the tree in order to obtain bootstrap values.