Skip to main content


Singleton molecular species delimitation based on COI-5P barcode sequences revealed high cryptic/undescribed diversity for Chinese katydids (Orthoptera: Tettigoniidae)



DNA barcoding has been developed as a useful tool for species discrimination. Several sequence-based species delimitation methods, such as Barcode Index Number (BIN), REfined Single Linkage (RESL), Automatic Barcode Gap Discovery (ABGD), a Java program uses an explicit, determinate algorithm to define Molecular Operational Taxonomic Unit (jMOTU), Generalized Mixed Yule Coalescent (GMYC), and Bayesian implementation of the Poisson Tree Processes model (bPTP), were used. Our aim was to estimate Chinese katydid biodiversity using standard DNA barcode cytochrome c oxidase subunit I (COI-5P) sequences.


Detection of a barcoding gap by similarity-based analyses and clustering-base analyses indicated that 131 identified morphological species (morphospecies) were assigned to 196 BINs and were divided into four categories: (i) MATCH (83/131 = 64.89%), morphospecies were a perfect match between morphospecies and BINs (including 61 concordant BINs and 22 singleton BINs); (ii) MERGE (14/131 = 10.69%), morphospecies shared its unique BIN with other species; (iii) SPLIT (33/131 = 25.19%, when 22 singleton species were excluded, it rose to 33/109 = 30.28%), morphospecies were placed in more than one BIN; (iv) MIXTURE (4/131 = 5.34%), morphospecies showed a more complex partition involving both a merge and a split. Neighbor-joining (NJ) analyses showed that nearly all BINs and most morphospecies formed monophyletic cluster with little variation. The molecular operational taxonomic units (MOTUs) were defined considering only the more inclusive clades found by at least four of seven species delimitation methods. Our results robustly supported 61 of 109 (55.96%) morphospecies represented by more than one specimen, 159 of 213 (74.65%) concordant BINs, and 3 of 8 (37.5%) discordant BINs.


Molecular species delimitation analyses generated a larger number of MOTUs compared with morphospecies. If these MOTU splits are proven to be true, Chinese katydids probably contain a seemingly large proportion of cryptic/undescribed taxa. Future amplification of additional molecular markers, particularly from the nuclear DNA, may be especially useful for specimens that were identified here as problematic taxa.


Taxonomic ambiguities and uncertainties are frequently generated due to cryptic or hidden species [1]. Species identification based on morphological characters requires experienced taxonomists [2]. Recently, DNA barcode have been recommended for the insect biodiversity evaluation [3]. DNA barcoding employs a single or a few standardized, highly variable and easily amplified DNA fragments for species identification [4, 5]. The 5′ portion of mitochondrial cytochrome c oxidase subunit I (COI-5P) has become the standard insect barcoding marker. Numerous studies rely on COI-5P as the only molecular information for insect species delimitation and identification [6,7,8,9]. DNA barcodes not only substantially improve the accuracy of species identifications but also accelerate the study of taxonomically difficult and hyperdiverse taxon. The introduction of DNA barcoding as an auxiliary method in taxonomy has many benefits. Firstly, DNA barcodes can lead to an easy assignment of specimens of certain life stages (e.g. eggs, larvae, nymphs or pupae) to known species [10]. Second, DNA barcoding requires a solid taxonomic background to use as a reference [11]. DNA barcodes can also add scientific value to standard museum specimens, as the information they contain is revealed through molecular analyses [12]. Third, DNA barcodes can accelerate the discovery of cryptic/undescribed species and have been incorporated into the new species description for several zoological groups [13,14,15,16].

Recently, several similarity-based species delimitation approaches, e.g. Barcode Index Number (BIN), REfined Single Linkage algorithm (RESL), Automatic Barcode Gap Discovery (ABGD) and a Java program uses an explicit, determinate algorithm to define Molecular Operational Taxonomic Unit (jMOTU); and clustering-based approaches, e.g. Generalized Mixed Yule Coalescent (GMYC), Bayesian implementation of the Poisson Tree Processes model (bPTP) have highlighted extensive inconsistency in morphological taxonomy [17, 18]. ABGD analysis generate MOTUs based on features in sequence distance distributions [19]. RESL employs single linkage clustering as a tool for the preliminary assignment of records into one MOTU [20]. BIN system developed within the Barcode of Life Data (BOLD, system to register the OTUs delineated by RESL [20]. jMOTU, a Java program uses an explicit, determinate algorithm to define MOTU [21]. GMYC and bPTP based on quite different models. GMYC applies single or multiple time thresholds to delimit species in a Maximum likelihood context, using ultrametric trees [22, 23]. bPTP is similar to GMYC, but used substitution calibrated trees [24]. Either single approach may lead to over- and/or underestimating species diversity. Here, we carried out different approaches for species delineation using DNA sequence data in order to have more robust results. The distinction between intraspecific and interspecific genetic divergence is critical for DNA barcoding; greater intraspecific divergence produces a greater likelihood of overlap with interspecific divergence [11]. For similarity-based approaches, the distance cut-off used for the determination of MOTUs is important, but arbitrary. Even no one threshold captures all species concepts or operational criteria [12]. Relaxed clustering-based methods that permit larger divergences within cohesive clusters may give even greater utility to similarity-based approaches [25]. The choice of a species delimitation method from molecular data has a considerable effect on estimated species entities and, thus, also on species richness estimates [26]. BOLD system [27] provides a unique environment for sharing data across projects; it not only supports all phases of the analytical pathway, from specimen collection to a tightly validated barcode library, but has already integrated many analysis tools. The BOLD system assigns a BIN for all barcode records. The BIN system will also help us to focus on those taxa that share the same BIN or split amongst multiple BINs [28].

Single-locus species delimitation methods have become popular due to the adoption of the DNA barcoding paradigm [17]. When delimiting putative species based on single-locus data, researchers should consider using both clustering- and similarity-based methods to account for the shortcomings of different methods [17, 29, 30]. Previously reported cases of high-failure rates in using traditional morphospecies definitions were largely resolved upon using MOTUs instead of traditionally described morphospecies, which suggested that some morphospecies may require taxonomic revision [25].

DNA barcoding has been used to document both grasshoppers [10, 31,32,33,34] and katydids [35, 36]. Unfortunately, these studies included only a few species and a limited number of specimens from each species, which prevents the rigorous assessment of species boundaries among closely related lineages and for calculation of intraspecific distances. Recently, Hawlitschek et al. (2017) presented a large-scale DNA barcode data set that includes 748 COI sequences from 127 species of Central European crickets, katydids and grasshoppers [37].

The Katydid diversity is rich but woefully underexplored in China. Both Mecopoda elongata and Gampsocleis gratiosa have a long history as singing pets in China. Researchers have reported that numerous species belong to the family Tettigoniidae Krauss, 1902 [38], but only a very limited COI-5P barcode records were available in the GenBank and BOLD systems. Barcode-based species identification relies on a comparison of its DNA barcode with those of determined individuals. To be effective, species-level assignments require a reference sequences database which represents all known species [39].

As different methods may yield inconformity conclusions [40], the accurate species identification and/or delimitation requires further integrative analysis. This study represents the first large-scale barcoding study of the family Tettigoniidae. Different molecular species delimitation methods such as BIN, RESL, jMOTU, ABGD, GMYC, and bPTP were applied to an unexplored Chinese katydid fauna. The main aims were (i) to present the largest species-level barcoding study for the Chinese katydids to date and then characterize the range of genetic divergence; (ii) to evaluate the correspondence between the identified morphospecies and the defined barcode groupings using molecular species delimitation methods; (iii) to infer species diversity and compare it to traditionally identified morphospecies sorting; and (iv) to test for the existence of a hidden diversity among otherwise well-defined taxa. Here, we define independently barcode groupings that remain morphologically indistinguishable from each other as cryptic or hidden species.

Materials and methods

Specimen collection and morphological identification

We collected 2576 katydid specimens throughout China. Specimens were fixed in absolute ethanol and were transferred to − 20 °C storage prior to genomic DNA extraction. Whenever possible, more than one location was sampled for each species to survey for intraspecific variation. Due to the polymorphism and remarkably wide distribution, a broader sampling (n ≥ 10 specimens) was particularly intended for 39 morphospecies, namely, Conanalus axinus, Conanalus pieli, Conocephalus bambusanus, Conocephalus gladiatus, Conocephalus longipennis, Conocephalus maculatus, Conocephalus melaenus, Deflorita deflorita, Ducetia japonica, Ducetia spina, Elimaea cheni, Euconocephalus pallidus, Euxiphidiopsis capricercus, Gampsocleis gratiosa, Gampsocleis sedakovii, Gampsocleis sinensis, Gampsocleis ussuriensis, Hemielimaea chinensis, Hexacentrus japonicas, Hexacentrus unicolor, Isopsera denticulate, Kuwayamaea brachyptera, Mecopoda niponensis, Orthelimaea trapzialis, Parapsyra nigrovittata, Parapsyra notabilis, Phaneroptera falcate, Pseudorhynchus concisus, Ruidocollaris truncatolobata, Ruspolia dubia, Ruspolia jezoensis, Ruspolia lineosa, Sinochlora szechwanensis, Tettigonia chinensis, Xizicus fascipes, Xizicus howardi, Xizicus kweichowensis, Xizicus spathulatus, and Xizicus szechwanensis.

Specimens were sorted morphologically and were taxonomically identified at least to the subfamily-level. For ease of information management, all unidentified specimens without a scientific name were assigned to an interim species (hereafter noted as BIN-species) by BINs provided by BOLD systems, and were noted by the generic/subfamily name plus sp.1, sp.2 and so on. For example, Atlanticus spp. were assigned to 19 BINs, so we noted them as Atlanticus sp1, Atlanticus sp2, Atlanticus sp3, etc. The specimen described as ‘undescribed genus’ were identified up to the subfamily level. All voucher specimens were preserved at the College of Life Sciences, Hebei University. More precise taxonomic determinations have been added for some specimens since their initial identification, and further taxonomic detail will be added to the BOLD systems as work progresses after publication.

DNA extraction, amplification and sequencing

DNA extraction from the leg muscle tissue of each specimen was carried out using the TIANamp Genomic DNA Kit (Tiangen Biotech, Beijing, China), following the manufacturer’s instructions. The universal primer pairs LCO1490/HCO2198 [41] were used to amplify and sequence the animal barcode region. PCRs were performed with 96-well plates. The reaction master mix consisted of 2400 μL 2 × Premix Taq™, 480 μL each primer (10 μmol/L) and 1152 μL water. This mixture was prepared for each plate, and each well contained 1 × Premix Taq™, 1 μmol each primer, 3~5 ng genomic DNA. The PCR profile was comprised of an initial denaturation step of 2 min at 95 °C, and 35 cycles of 30 s at 94 °C, 40 s at 50 °C and 1 min at 72 °C, with a final extension of 7 min at 72 °C.

Amplicons were checked through a 1% agarose gel and bi-directional sequencing was performed at GENEWIZ (Tianjin, China). Sequences were manually edited and assembled into a consensus sequence using SeqMan Pro [42]. Consensus sequences, specimen collection data, specimen images and sequence trace files were uploaded to the Barcode of Life Data System (BOLD) and are available to the public domain as part of the project DNA Barcoding to Katydids from China (DBKC).

Data analysis

We constructed our primary dataset from BOLD systems, including all public records that were geographically limited in China, and with a length ≥ 600 bp. To reduce computational requirements, we divided our entire dataset into four subsets by a subfamily or consonant subfamilies based on the four monophyletic lineages that were recognized by our previous mitogenomic Bayesian inference analysis with the site-heterogeneous CAT-GTR model [43]. The DS-DBCHL dataset is composed of 596 barcode sequences from Conocephalinae, Hexacentrinae, and Lipotactinae. The DS-DBMEC dataset is composed of 376 barcode sequences from Meconematinae. The DS-DBPPM dataset is composed of 993 barcode sequences from Phaneropterinae, Pseudophyllinae, and Mecopodinae. The DS-DBTB dataset is composed of 200 barcode sequences from Tettigoniinae and Bradyporinae.

Sequence analysis module of BOLD systems

BIN was used as a registry for the records on the BOLD systems [27], which provided a means of confirming the concordance between barcode sequence clusters and species designations [20, 28]. Cluster sequence analysis using RESL was independent of the BIN registry of BOLD systems [20]. Genetic distances were calculated and summarized using the “Distance Summary” and “Barcode Gap Analysis” tools on BOLD systems [27]. Barcode gap analysis provides the distribution of distances within each species and the distance to the nearest neighbor (NN) of each species. Species are tested for the presence of the barcode gap. The NN distance is the genetic distance between a species and its closest congeneric relative. All sequences (> 600 bp) were aligned using MUSCLE [44], phylogeny model used was the Kimura 2-parameter (K2P) [45], and pairwise deletion of missing data was done. The correlation between the maximum intraspecific variation (K2P) against record count and maximum geographic extent (km) of sampled individuals was determined for each species sampled from more than one sites.

Similarity-based methods: jMOTUs and ABGD

In performing the additional “Sequence analysis module” of BOLD, we also applied two similarity-based methods to generate MOTUs. Each of the four datasets was analyzed with jMOTU_define 2.04 [21, 46] using different cut-offs (from 0 to 25 bp). ABGD analysis was performed to sort the sequences into hypothetical species that are based on the barcode gap, which can be observed whenever the divergence among organisms that belonging to the same species is smaller than the divergence among organisms from different species [19].

Clustering-based methods: GMYC and bPTP

Previous studies suggest that an additional bias may be introduced for clustering-based methods when duplicate haplotypes are not removed [17]. Prior to species delimitation analyses with two clustering-based methods, we applied DAMBE [47] to remove duplicate haplotypes. A total of 530 unique haplotypes from the DS-DBPPM dataset, 147 haplotypes from DS-DBTB, 158 haplotypes from DS-DBMEC, and 390 haplotypes from DS-DBCHL were included in the further analyses. Ultrametric trees were estimated with BEAST v1.8.3 [48] using a Yule speciation prior and an uncorrelated lognormal relaxed clock. The best-fitting substitution models were selected under the Bayesian Information Criteria (BIC), as was implemented in jModelTest 2.1.7 [49]. Each of the four datasets was analyzed for 200 million iterations with the first 10% discarded as burn-in. Posterior probabilities (PP) were estimated under a sampling frequency of every 10,000 steps. Tracer v.1.6 ( was used to determine when the analyses became stable and to check whether the effective sample size (ESS) values were greater than 200, as recommended by Drummond et al. (2007). The consensus trees obtained before the Markov chain reached stable and convergent likelihood values were discarded as burn-in with TreeAnnotator v.1.7. The resulting ultrametric trees were used for both single-threshold GMYC (sGMYC) [22] and multiple-threshold GMYC (mGMYC) [23] analyses using the Splits [50] and Ape [51] libraries.

The coalescent clustering-based method (bPTP) was performed using the online server ( and the Bayesian Inference trees from MrBayes 3.2 [52]. We ran bPTP analyses for 500,000 MCMC generations with a thinning of 500 and a burn-in of 0.1. Convergence of the MCMC chain was assessed as recommended by Zhang et al. (2013). Outgroups were pruned before conducting bPTP analyses to avoid bias that may arise if some of the outgroup taxa were too distantly related to the ingroup taxa [24].

Comparison morphospecies and MOTU of species delimitation method outputs

All NJ-K2P trees of unique COI-5P haplotypes were performed using MEGA v7.0 [53]. The results of the species delimitation methods were summarized on the NJ-K2P trees with a midpoint root. Four different taxonomic scenarios between morphospecies (equated with BIN-species for unidentified specimens) and MOTU clustering methods outputs can be distinguished: (i) ‘MATCH’, whereby the members of a species were placed in one MOTU that had no other members; (ii) ‘MERGE’, whereby the members of a species were placed in one MOTU together with members from another species; (iii) ‘SPLIT’, whereby the members of a species were assigned to more than one MOTU that had no other specimens from another species; and (iv) ‘MIXTURE’, whereby each species show a more complex partition involving both ‘MERGE’ and ‘SPLIT’. To further compare the results of different species delimitation methods, we also employed the adjusted Wallace coefficients analysis [54] to quantify MOTUs agreement with Linnaean species labels or among MOTUs from different species delimitation methods through the website Comparing Partitions ( [55]. Here, we excluded singletons and only discussed species or MOTUs that are represented by more than one specimen. Finally, MOTUs were defined considering only the clades represent groups of barcodes recovered in at least four of the seven species delimitation methods [20].


DNA was extracted from 2576 Chinese katydid specimens, of which 2131 specimens (82.73%) were successfully sequenced for COI-5P barcode. All records were removed that less than 600 bp, contained contaminants, had stop codons, flagged as misidentifications or errors. In summary, we generated 2131 original COI-5P sequences from 131 morphospecies, including 528 specimens that were identified to genus level and one specimen that was placed at the subfamily level. The remaining unidentified lineages were represented using BINs, because they were either unable to be reliably identified based on the available reference materials or they were still undescribed. Additional 34 published COI-5P barcode sequences (Additional file 1: Table S1) were retrieved from GenBank in the BOLD system and included for further analyses. The entire dataset containing 2165 Chinese katydids COI-5P barcode sequences comprised 1225 distinct haplotypes, and represented 60 genera, 9 subfamilies of the family Tettigoniidae Krauss, 1902, including Bradyporinae (n = 1), Conocephalinae (n = 490), Hexacentrinae (n = 102), Lipotactinae (n = 4), Meconematinae (n = 376), Mecopodinae (n = 53), Phaneropterinae (n = 861), Pseudophyllinae (n = 79), and Tettigoniinae (n = 199). Nearly a third of the morphospecies (39/131 = 29.77%) included 10 or more specimens. The number of barcode sequences per species varied from 1 up to 142 in the commonly occurring Ducetia japonica, dispersed throughout China.

Intra- and interspecific genetic divergences

The intra- and interspecific genetic divergences within different taxonomic ranks are detailed in Table 1. The intraspecific divergence for 109 morphospecies represented by more than one specimen averaged 1.54% (ranging from 0 to 27.45%). However, the intraspecific divergence for 77 unidentified BIN-species represented by more than one specimen averaged 0.39% (ranging from 0 to 2.81%). The identified morphospecies, which were assigned to more than one BIN, were the major cause of higher intraspecificity. The mean interspecific divergence for identified morphospecies and unidentified BIN-species at genus level were 15.35% (ranging from 0 to 28.74%) and 12.86% (ranging from 1.07 to 28.07%), respectively. The mean interspecific divergence for identified morphospecies and unidentified BIN-species at family level were 22.29% (ranging from 2.16 to 32.88%) and 21.67% (ranging from 3.27 to 32.78%), respectively. The normalized mean intraspecific and minimum interspecific distance were 1.40 ± 0.02 and 0% for 109 identified morphospecies, in contrast to 0.45 ± 0.01 and 1.07% for 77 unidentified BIN-species, respectively (Table 2).

Table 1 Summary of K2P distances between barcode sequences at each taxonomic level
Table 2 Normalized divergence statistics

For the entire dataset, the distance of 22 morphospecies and 16 BIN-species to its NN was smaller than 2%, in which the distance of 16 of 22 morphospecies to NN was less than maximum intraspecific distance (Additional file 2: Table S2). Meanwhile, there are 13 morphospecies and one BIN-species, the distance to NN was larger than 2%, but less than the maximum intraspecific distance (Table 3). Deep intraspecific divergences (> 2%) overlapping with the distance to NN were detected in 21 morphospecies and one BIN-species, namely, Conanalus robustus (9.25%), Conocephalus bidentatus (18.85%), Conocephalus longipennis (17.81%), Euconocephalus pallidus (3.45%), Euxiphidiopsis capricercus (27.45%), Gampsocleis gratiosa (4.91%), Gampsocleis sedakovii (4.9%), Gampsocleis ussuriensis (3.78%), Hexacentrus japonicus (4.4%), Mecopoda niponensis (7.97%), Phyllomimus sinicus (8.09%), Pseudorhynchus pyrgocoryphus (10.36%), Ruidocollaris truncatolobata (6.59%), Ruspolia dubia (3.27%), Ruspolia yunnana (3.78%), Sinochlora szechwanensis (7.92%), Xiphidiopsis autumnalis (12.14%), Xiphidiopsis bituberculata (2.65%), Xiphidiopsis gurneyi (13.02%), Xizicus spathulatus (9.43%), Xizicus howardi (7.79%), and Ruidocollaris sp. 7 (2.81%). The linear regressions analysis indicate the maximum intraspecific variation (K2P) was significantly correlated with record count (P < 0.001) and maximum geographical extent of sampled individuals (P < 0.001), but had limited explanatory power (Adjusted R-square = 0.080 and 0.142) (Fig. 1).

Table 3 Highlighted records which the distance to nearest neighbor (NN) less than 2% or Max intra-specific distance
Fig. 1

Maximum intraspecific variation (K2P) against record count (a) and maximum geographic extent (b) of sampled individuals. a linear regression, y = 0.0007x + 0.0174, Adjusted R-square = 0.142, P < 0.001; b linear regression, y = 2E-05x + 0.0131, Adjusted R-square = 0.080, P < 0.001

Comparing BIN-species and morphospecies identifications

The BOLD-implemented refined single-linkage algorithm that provided the BIN assignments used a 2.2% p-distance seed threshold but then refined the groupings for individual BINs and neighboring clusters based on the level of continuity in the distribution of genetic divergences among sequences [20, 56]. In total, 1635 sequences of 131 morphospecies were assigned to 196 BINs, including 52 singleton BINs, 136 concordant BINs, and 8 discordant BINs. Only one barcode of Ducetia japonica RBTC523–16 was not assigned to any BIN, because it contained more than 1% Ns. It was excluded from the analyses of species delimitation except for NJ-K2P tree. These cases of discrepancy are discussed in more detail in subsequent sections. The unidentified 529 specimens without formal (binomial) names were assigned to 150 BIN-species, including 77 concordant BINs and nearly half (n = 73) singleton BINs (Table 4). On average, we identified a different BIN for every 6.25 specimens, and 44 of 346 BINs (12.72%) included not less than 10 specimens, in which only two BINs (BOLD:ACE7214 and ADB5001) included more than 100 specimens.

Table 4 MOTU number of Chinese Katydids inferred by each single-locus species delimitation methods

There are four different taxonomic scenarios between BINs and morphospecies, including MATCH, SPLIT, MERGE, or MIXTURE [20]. Approximately 64.9% (n = 84) of the identified morphospecies were MATCH, and 1:1 corresponded with BINs (including 63 concordant BINs and 21 singleton BINs) (Table 5). The discordances between morphospecies and BINs included the members of multiple morphospecies pooled into one BIN (MERGE), a morphospecies split into more than one BIN (SPLIT), or both (MIXTURE). We found 18 morphospecies (6 pairs, two triads) that shared BINs (Table 6), in which 14 morphospecies were MERGE and shared its unique BIN with other species. The remaining 4 morphospecies were MIXTURE, and were also assigned to an additional BIN, including Conocephalus longipennis (2), Euconocephalus pallidus (1), Ruspolia dubia (3), and Xizicus howardi (2). Excluded singleton species, with a relatively high percentage of morphospecies, were SPLIT, and were assigned to more than one BIN (34/109 = 31.19%), namely, Conanalus pieli (7BINs), Conanalus robustus (3BINs), Conocephalus bidentatus (3BINs), Conocephalus longipennis (3BINs), Conocephalus maculatus (11BINs), Ducetia japonica (3BINs), Elimaea nautica (3BINs), Euconocephalus pallidus (2BINs), Euxiphidiopsis capricercus (2BINs), Euxiphidiopsis spathulata (2BINs), Gampsocleis carinata (3BINs), Gampsocleis gratiosa (5BINs), Hexacentrus japonicus (4BINs), Isopsera denticulata (2BINs), Kuwayamaea brachyptera (4BINs), Mecopoda niponensis (4BINs), Orthelimaea trapzialis (2BINs), Phyllomimus detersus (2BINs), Phyllomimus sinicus (6BINs), Pseudorhynchus pyrgocoryphus (2BINs), Ruidocollaris truncatolobata (2BINs), Ruspolia dubia (4BINs), Ruspolia lineosa (2BINs), Ruspolia yunnana (3BINs), Sinochlora szechwanensis (2BINs), Tettigonia chinensis (2BINs), Xiphidiopsis autumnalis (2BINs), Xiphidiopsis gurneyi (4BINs), Xizicus fascipes (2BINs), Xizicus howardi (3BINs), Xizicus kweichowensis (2BINs), Xizicus magnus (3BINs), Xizicus spathulatus (3BINs), and Xizicus szechwanensis (2BINs) (Table 7).

Table 5 Morphospecies-BIN perfect matches in this study
Table 6 Discordance BINs report for different morphospecies were assigned to one BIN
Table 7 Summary of the 37 morphospecies were assigned to more than one BIN

Monophyletic morphospecies or BIN-species recovered by NJ-K2P trees

The NJ-K2P trees based on COI-5P haplotypes are shown in Additional files 3, 4, 5 and 6. The BIN 2.2% seed threshold was calibrated against morphological species using a selected groups of taxa: bees, butterflies and moths, fish, and birds [20]. The DBCHL dataset included 596 COI-5P barcode sequences, and was assigned to 77 BINs, including 13 singleton BINs, 61 concordant BINs, and 3 discordant BINs. NJ analysis with 390 COI-5P haplotypes sequences showed that all BIN-species represented by more than one specimen formed a monophyletic clade, except for BOLD:ACH8981, ACN8107 and ADE4977 (Additional file 3). The 579 sequences represented 40 identified morphospecies, in which 39 species included more than one specimen. 30 identified morphospecies formed monophyletic clusters. In addition, three specimens of Conanalus brevicaudus shared a unique COI-5P haplotype. The members of Ruspolia dubia, R. jezoensis, and R. liangshanensis were grouped jointly and formed a larger monophyletic clade with a low divergence (Additional file 3). The BOLD:ADE4977 includes a triad of species, including all members of Ruspolia jezoensis (n = 10), R. liangshanensis (n = 5) and most of the R. dubia (n = 38). The remaining members of Ruspolia dubia were assigned three BINs, BOLD:ADE5391 (n = 5), ADE5392 (n = 1), and ACD5503 (n = 6). The members of Euconocephalus pallidus were split into two closely related clades (Additional file 3), and clade B1 was formed by two specimens of E. pallidus and all specimens of E. nasutus. The singleton Pseudorhynchus sp. is nested within the P. pyrgocoryphus clade (Additional file 3). Conanalus robustus was split into two relatively distant clades (Additional file 3). Clade D1 corresponded to BOLD:ADB9302, and the members of clade D2 with high divergences were assigned to two BINs (BOLD:ADB9301 and ADB9303). Two specimens identified as Conocephalus japonicus is nested within the C. longlpennis clade (Additional file 3). The widely distributed morphospecies Ruspolia lineosa was monophyletic but contained two deeply subclusters. Almost all species delimitation analyses suggested that R. lineosa split into two putative species, R. lineosa BOLD:ACD5256 and ACD5257. Only ABGD suggested that R. lineosa to be a distinct species. Both sGMYC and mGMYC subsplit R. lineosa BOLD:ACD5257 into two parts.

The DBMEC dataset included 376 COI-5P barcode sequences and was assigned to 56 BINs, which included 24 singleton BINs, 28 concordant BINs, and 4 discordant BINs. NJ analysis with 158 COI-5P haplotypes sequences showed that all BIN-species with more than one sequence formed a monophyletic clade (Additional file 4). The 341 sequences represented 40 identified morphospecies, in which 29 species included more than one specimen. 22 identified morphospecies revealed nonoverlapping monophyletic clusters. The remaining 7 morphospecies were not monophyletic. Euxiphidiopsis capricercus was split into two reciprocally monophyletic clusters (Additional file 4) and corresponded to two BINs (BOLD:ADB5001 and ADE2467). Both Xiphidiopsis gurneyi and X. autumnalis were split into two reciprocally monophyletic clusters, and four clusters were grouped jointly as a separate clade (Additional file 4). The singleton species identified as Xiphidiopsis maculatus was embedded into the Xizicus spathulatus clade (Additional file 4). The members of Xizicus howardi were recovered in three reciprocally monophyletic clusters (Additional file 4) and corresponded to three BINs (BOLD:ACD5539, ADB5688 and ADE3141). Xizicus concavilaminus and X. kulingensis were grouped jointly (Additional file 4), and shared a single BIN (BOLD:ADB3332). Xizicus tinkhami and X. laminatus were grouped jointly (Additional file 4) and shared a single BIN (BOLD:ADB5868). Xiphidiopsis bituberculata and X. minorincisus were grouped jointly (Additional file 4) and shared a single BIN (BOLD:ADB3697). Due to their small size, the active dispersal abilities of Meconematinae katydids were highly limited.

The DBPPM dataset included 993 COI-5P barcode sequences and was assigned to 181 BINs, which included 76 singleton BINs, and 105 concordant BINs. Only one barcode (Ducetia japonica RBTC523–16) without BIN records corresponded to sequences that did not fulfill the barcode compliance standards. NJ analysis with 530 COI-5P haplotypes sequences showed that all BIN-species with more than one sequence formed a monophyletic clade (Additional file 5). The 551 sequences represented 43 identified morphospecies, in which 33 species included more than one specimen. Twenty-seven identified morphospecies were revealed in nonoverlapping monophyletic clusters. Ruidocollaris truncatolobata was split into two relatively distant clusters (Additional file 5) and corresponded to two BINs (BOLD:ACD6433 and ACD7529). One specimen that was identified as Mecopoda sp. is nested within the clade of M. niponensis (Additional file 5). Sinochlora szechwanensis was split into two relatively distant clades (Additional file 5). Clade C1 was formed exclusively by specimens from Yuexi, Anhui and was closely related to the species Sinochlora longifissa. Clade C2 was formed by the remaining specimens, and was closely related to Sinochlora sp2 (BOLD:ADB3463). Two specimens that were identified as Phyllomimus sp14 (BOLD:ADB3808) and Phyllomimus sp15 (BOLD:ADB6425) are nested within the P. sinicus clade (Additional file 5).

The DBTB dataset included 200 COI-5P barcode sequences, and was assigned to 32 BINs, which included 12 singleton BINs, 19 concordant BINs, and 1 discordant BINs. NJ analysis with 147 COI-5P haplotypes sequences showed that all BIN-species with more than one sequence formed a monophyletic clade, except for BOLD: ADA6837 and ADB3445 (Additional file 6). A total of 8 identified morphospecies were represented by 164 sequences, including 7 species that were represented by more than one specimen. The NJ analysis based on K2P distances revealed nonoverlapping clusters for 4 identified morphospecies, Chizuella bonneti, Gampsocleis carinata, G. gratiosa, and Tettigonia chinensis. In contrast, Gampsocleis sedakovii, G. sinensis, and G. ussuriensis were grouped jointly. The remaining 36 sequences were provisionally assigned into 19 putative species based on BINs.

Concordance among MOTUs from similarity-based species delimitation methods

Because of its strong taxonomic performance and speed, RESL was adopted to generate MOTUs for the barcode sequences on BOLD systems [27]. The results of RESL analyses generated 349 MOTUs, which had only small discrepancies in comparison with the BINs (Table 4). For the DBCHL dataset, RESL analysis generated 71 MOTUs. The differences between BINs and RESL were that (i) four BINs (BOLD:ADE4977 representing Ruspolia dubia, R. jezoensis, R. liangshanensis, and ADE5391, ADE5392, ACD5503 representing R. dubia) were pooled in one MOTU; (ii) two BINs (BOLD:ACH8981 and ADE5243 representing Ruspolia yunnana) were pooled in one MOTU; (iii) two BINs (BOLD:ACD6726 representing Euconocephalus pallidus and E. nasutus, and BOLD:AAP6087 representing E. pallidus) were pooled in one MOTU. For the DBMEC dataset, RESL analysis generated 56 MOTUs. For the DBPPM dataset, RESL analysis generated 180 MOTUs. Hemielimaea omeishanica (BOLD:ACD5212) were split into two MOTUs. Ruidocollaris sp7 (BOLD:ADB6075) were split into four MOTUs. Meanwhile, RESL recovered Ducetia japonica as one MOTU, which was assigned to three BINs (BOLD:ACD7324, ADB6191, ACE7214). Two singleton BINs, BOLD:ADB3808 representing Phyllomimus sp14 and ADB6425 representing Phyllomimus sp15, were pooled into one MOTU. For the DBTB dataset, RESL analysis generated 42 MOTUs. The members of BOLD:AAY1322 (representing Gampsocleis sedakovii, G. sinensis, G. ussuriensis) were split into 10 MOTUs, and members of BOLD:ADA6837 (representing G. gratiosa) were split into three MOTUs.

A 2% divergence criterion has been proposed as a general rule-of-thumb for species boundaries with COI-5P [6]. The results of jMOTU analyses at different cutoffs (from 0 to 25 bp) are shown in Fig. 2. A total of 318 MOTUs were determined by a 13 bp (~ 2%) distance cut-off, including 61 MOTUs of DBCHL dataset, 57 MOTUs of DBMEC dataset, 169 MOTUs of DBPPM dataset, and 31 MOTUs of DBTB dataset. ABGD was only based on similarity among sequences, without considering the phylogenetic relationships [19]. The perfect match of ABGD approaching between the initial and the recursive partitions occurred at nucleotide divergence values of 2.15%. The ABGD analyses generated the most conservative results (Table 4) and inferred 255 MOTUs, including 42 MOTUs of DBCHL dataset, 53 MOTUs of DBMEC dataset, 140 MOTUs of DBPPM dataset, and 20 MOTUs of DBTB dataset (Fig. 3).

Fig. 2

Correspondence between the genetically putative molecular species (MOTUs) number and the cut-off value (bp) generated by jMOTU

Fig. 3

Correspondence between the genetically distinct MOTUs number and prior maximal distance (%) by ABGD based on K2P model

Concordance among MOTUs from clustering-based species delimitation methods

Both sGMYC and mGMYC coalescence-based clustering of the specimens were partitioned in the data far more than in all of the other methods (Table 4). The mGMYC analysis was by far the most sensitive of the methods compared, inferring a total of 397 entities, which was slightly higher than sGMYC (n = 382). For the DBCHL dataset, the sGMYC analysis identified 87 ML entities (95% confidence interval = 77–95): 71 ML clusters (95% confidence interval = 63–76) and 16 singletons. The mGMYC analysis identified 5 independent switches between speciation and coalescent processes, resulting in 94 ML entities (95% confidence interval = 85–105): 73 ML clusters (95% confidence interval = 67–77) and 21 singletons. For the DBMEC dataset, the sGMYC analysis identified 66 ML entities (95% confidence interval = 48–73): 40 ML clusters (95% confidence interval = 31–40) and 26 singletons. The mGMYC analysis identified 4 independent switches between speciation and coalescent processes, resulting in 60 ML entities (95% confidence interval = 51–76): 36 ML clusters (95% confidence interval = 31–38) and 24 singletons. For the DBPPM dataset, the sGMYC analysis identified 197 ML entities (95% confidence interval = 151–219): 112 ML clusters (95% confidence interval = 93–114) and 85 singletons. The mGMYC analysis identified 6 independent switches between speciation and coalescent processes, resulting in 206 ML entities (95% confidence interval = 179–234): 106 ML clusters (95% confidence interval = 97–113) and 100 singletons. For the DBTB dataset, the sGMYC analysis identified 32 ML entities (95% confidence interval = 20–43): 24 ML clusters (95% confidence interval = 16–29) and 8 singletons. The mGMYC analysis identified 4 independent switches between speciation and coalescent processes, resulting in 37 ML entities (95% confidence interval = 29–59): 23 ML clusters (95% confidence interval = 22–27) and 14 singletons. bPTP analyses inferred 312 MOTUs with wide confidence intervals from MCMC analyses, including 55 MOTUs of DBCHL dataset, 56 MOTUs of DBMEC dataset, 178 MOTUs of DBPPM dataset, and 23 MOTUs of DBTB dataset.

Previous studies have found that in some cases, GMYC could lead to an overestimation of the number of species [57, 58]. GMYC requires prior construction of a ultrametric tree, which does not necessarily reflect the real divergence between species [30]. Alternatively, PTP estimates branching processes using the expected number of substitutions (vs. time in GMYC) and thus utilizes a nonultrametric phylogenetic tree as input [17]. Moreover, in contrast to GMYC, bPTP appeared less sensitive to the sampling regime [17].

Bidirectional concordance among species delimitation methods with adjusted Wallace coefficients

The adjusted Wallace coefficients were used to compare the bidirectional concordance among species delimitation methods for the identified katydid specimens dataset (Table 8). There were markedly directional results in discriminatory power between the molecular species delimitation methods and morphospecies or across molecular species delimitation methods. For example, the adjusted Wallace coefficient value (0.954) from BIN to morphology meant two specimens within a BIN had a 95.4% chance of belonging to the same morphospecies. In contrast, two specimens within a morphospecies had only a 70.9% chance of belonging to the same BIN. Overall, molecular species delimitation methods had a strong explanatory ability (0.872–0.969) for morphospecies; in contrast, the morphospecies had a generally low explanatory ability (0.391–0.995) for molecular species delimitation results. Both sGMYC and mGMYC were less concordant with morphospecies than other molecular species delimitation results. GMYC inferred a substantially unrealistically high number of katydid MOTUs (Table 4). The morphospecies was best able to explain the results of both ABGD (0.955) and bPTP (0.940). They generally exhibited a more modest ability to explain the molecular results in comparison with all other methods (0.996–1 for ABGD and 0.990–0.998 for bPTP to be explained by all other molecular species delimitation results). The differences among the remaining methods (e.g., BIN, jMOTU, and RESL) in their concordance to the current taxonomy were modest.

Table 8 Bidirectional concordance among clustering methods for identified katydids specimens dataset using Adjusted Wallace and 95% CI


Species was defined as lineages that evolve separately from each other [59]. Determining the species boundaries is one of the central debates in biology [18]. DNA barcoding was widely used for species identification and/or delimitation. Recent research on Central European Orthoptera found that ninety-three of these 122 species (76.2%, including all Ensifera) could be reliably identified using DNA barcodes [37]. In European diving beetles, 36% of multiply sampled species were nonmonophyletic [60]. Our study provides barcode data for 131 identified morphospecies and 148 unidentified BIN-species of Chinese katydids. There was a perfect correspondence between BIN membership and morphospecies in 83 cases, while another 34 species split into more than one BIN. more than one species merges as one BIN or in a combination of merges and splits. The maximum intraspecific distance is less than 3% in 74 of the 109 identified morphospecies (67.89%) represented by multiple individuals. Maximum intraspecific distance is less than 3% in all unidentified BIN-species. Our results revealed a much higher diversity in Chinese katydids than the current taxonomy suggests. There are more katydid species to be described and cryptic lineages within currently recognized species.

COI-5P barcode and BIN sharing

The causes for barcode and BIN sharing in closely related species include imperfect taxonomy [32], nonfunctional nuclear-encoded mitochondrial pseudogenes (Numts), hybridization, and incomplete lineage sorting [37, 61, 62]. Previous studies have found evidence for frequent hybridization across orthopteran closely related species, such as the genus Chorthippus [63], Aglaothorax [64], Tetrix [65]. COI-5P barcodes sharing was found in two cases, Gampsocleis sedakovii (GHF077_16) vs. G. sinensis (GHF074_16), as well as G. sinensis (RBTC480_16, RBTC1222_16, RBTC1209_16, RBTC1193_16) vs. G. ussuriensis (GHF028_16). The barcodes of G. sedakovii, G. sinensis and G. ussuriensis pooled into one discordance BIN (BOLD:AAY1322), but was not supported by our other analyses. The morphological high similarity between G. sinensis and G. ussuriensis. Another possible reason is that G. ussuriensis might in fact synonymised to G. sinensis. Meanwhile, G. sedakovii and G. ussuriensis occur in sympatry over large parts of their distribution ranges. Hybridization in sympatry has resulted the transfer barcodes from G. sedakovii to G. sinensis and/or G. ussuriensis, causing COI-5P barcode and BIN sharing.

Hausmann et al. (2013) suggested that cases of BIN sharing among allopatric, slightly divergent genetic clusters represent recently separated lineages that have recently speciated or are still undergoing genetic differentiation [28]. The bPTP analysis indicated four Ruspolia species, R. dubia, R. jezoensis, R. liangshanensis, and R. yunnana, pooled into one MOTU. Meanwhile, both RESL and jMOTU analysis indicated R. dubia, R. jezoensis, and R. liangshanensis pooled into one MOTU. Consistent results has also been previously observed in which with regard to R. jezoensis synonymised to R. dubia, and R. liangshanensis may be recently separated from R. dubia [36]. Conocephalus japonicus is nested within the C. longipennis cluster. These species formed very recently indeed, and young species (incomplete lineage sorting) remain within its sister species’ coalescent lead to BIN sharing.

In the last few years, some locally distributed Xizicus species have been described. Our analyses found X. concavilaminus and X. kulingensis pooled into one MOTU (BOLD:ADB3332). X. laminatus and X. tinkhami pooled into one MOTU (BOLD:ADB5868) except for GMYC analyses. Xizicus rehni and X. howardi share BOLD:ACD5539, but still exhibit subclusters that separate species at a very low distance. Three discordance BINs (BOLD: ACD5539, ADB3332, and ADB5868) between Xizicus species were supported by most analyses. This phenomenon might reflect their relatively recent split or the current taxonomy of Xizicus too detailed.

Morphospecies split into more than one BIN

Our results demonstrate the existence of more than one separate lineage in several katydids with wide geographic distribution range. Some species split into more than one BIN may referred to sister clusters on the barcode trees may representing true potential cryptic diversity [31]. Many species with wide geographic distribution range were placed in either a single BIN or a few, but Conocephalus maculatus was outliers, being assigned to 11 BINs. 10 of 11 BINs in C. maculatus represented by two or more specimens. Four BINs of C. maculatus co-occur at different sites across China, such as BOLD:ACN8107 from Hainan, Yunnan, Guangxi; BOLD:ADB6356 from Xizang, Yunnan; BOLD:ACD2116 from Zhejiang, Jiangxi; BOLD:ADB6002 from Xizang, Yunnan. C. maculatus is one of the most widespread species of genus Conocephalus, and exhibit one monophyletic cluster. C. maculatus is very likely to represent a species complex. Previous research found specimens of one Canadian spiders Tetragnatha versicolor were assigned to 20 BINs [9].

The maximum intraspecific distances possessed up to 18.85% in Conocephalus bidentatus. All barcodes of C. bidentatus exhibit one monophyletic cluster, and clearly distinct subclusters reflected by three BINs. Three BINs within C. bidentatus reflect geographic clustering with BOLD:ADB6577 from Fujiang and Zhejiang, BOLD:ADB9596 from Sichuang and BOLD:ADC0531 from Guizhou. Our analyses support C. bidentatus split into three MOTUs except for ABGD and mGMYC treating ADB6577 and ADB9596 as one MOTU.

Xiphidiopsis autumnalis and X. gurneyi form one monophyletic cluster. The two BINs within X. autumnalis reflect geographic clustering with ADE1666 from Hainan, ADE1667 from Guangxi. Meanwhile, the four BINs (BOLD:ADB7052, ADE1668, ADE1669, and ADE1670) within X. gurneyi reflect different sampled localities. Xizicus howardi split into three BINs: ACD5539 from Guangxi, Henan, Hubei, and Zhejiang, ACD5688 from Zhejiang, and ADE3141 from Zhejiang, presumably with a species status. Two or more species are cryptic if they are morphologically similar, biologically distinct, and misclassified as a single species [66]. Cryptic species complexes, in which the component taxa have not diverged morphologically too much, are very difficult to identify, and their discovery is frequently a matter of chance [67]. This phenomenon may suggest possible cryptic species or there are more species to be described.

Conflicts among different species delimitation approaches are very common. Both Gampsocleis carinata and G. gratiosa were split into two MOTUs except for ABGD analysis. The two BINs (BOLD:ADA6038 and ADA5568) of Gampsocleis carinata were predominantly matched to most other clustering methods, which potentially implies detectable intraspecific diversity within Gampsocleis carinata and G. gratiosa, or the probable existence of more than one species. The species Tettigonia chinensis was supported by three clustering methods (ABGD, mGMYC and bPTP). Meanwhile, it was split into two BINs (BOLD:ACD6622 and ACD6623) and was supported by both RESL and jMOTU results.

DNA-based species delimitation may be compromised by limited sampling effort and species rarity, including “singleton” representatives of species, which hampers estimates of intra- versus interspecies evolutionary processes [68]. A broader intraspecific sampling is a critical step for increasing the success of species identification, and a special effort was made to achieve this aim [69]. Previous studies demonstrated broader geographical sampling decreases the barcoding gap between species and hence reduces the accuracy of DNA barcoding [60]. Ducetia japonica have been found are distributed over a huge area extending from Pakistan in the West to the Solomon Islands in the East and from Northern China in the North to northern Australia in the South [70]. We have a broader sampling of Ducetia japonica distribution in China. Our results demonstrate the existence of three separate lineages (BOLD:ACD7324, ADB6191, and ACE7214) in D. japonica. The different song types indicated clearly that D. japonica as presently understood is not a homogeneous, extremely widespread species, but a complex of several distinct species [70].

The nonfunctional nuclear-encoded mitochondrial pseudogenes (Numts) are a potential source of barcoding error [31]. Previous studies have showed that Numts had devastating influence on DNA barcoding results and were very hard to detect [71,72,73]. Our analyses supported Euxiphidiopsis capricercus split into two BINs (BOLD:ADB5001 and ADE2467) except for GMYC. Only one specimen (HLXX121–16) correspond to BOLD:ADE2467, and differed from the remaining E. capricercus specimens by 27.45%. E. capricercus HLXX121–16 as the NN of Gampsocleis gratiosa, the distance only 2.49%. Note that excluding E. capricercus BOLD:ADE2467 results the intraspecific distances (1.23%) significant decrease. BOLD:ADE2467 far from the Euxiphidiopsis cluster on the NJ-K2P tree. The sequence divergence within E. capricercus and geographical distance is correlated significantly. E. capricercus HLXX121–16 (BOLD:ADE2467).


No universal barcode gap was observed in our four data sets. There are moderately variable results from different delimitation methods. Our research supported the contention of Ortiz and Francke (2016) contention that combining evidence from multiple delimitation methods obtains better-supported results [18]. To diminish the probability of species under- or overestimation solutions, we determined separate MOTUs that were recovered in at least four of the seven species delimitation analyses. Excluding singletons (22 identified species, 125 BINs), we recognized 62 robustly supported identified morphospecies, and 166 BIN-species. The GMYC model exhibited a characteristic “overestimating” solutions.

If most MOTU splits detected in this study reflect cryptic/undescribed taxa, the true species count for Chinese katydids could be a large proportion higher than currently recognized. Moreover, only less than 20% species (50 of 279) were represented by not less than 10 specimens, and expanded sample sizes might reveal more barcode splits. Here, we refrain from taxonomic descriptions, as this requires a thorough morphological and taxonomic study for each putative taxon. It is also important to note that there could be noise in our results, potentially due to considerable unidentified specimens. Nevertheless, our results support COI-5P efficacy for rapid delimitation of katydid species and for indicating likely cryptic/undescribed species for further exploration.


  1. 1.

    Lobo J, Teixeira MAL, Borges LMS, Ferreira MSG, Hollatz C, Gomes PT, Sousa R, Ravara A, Costa MH, Costa FO. Starting a DNA barcode reference library for shallow water polychaetes from the southern European Atlantic coast. Mol Ecol Resour. 2016;16(1):298–313.

  2. 2.

    Fiser Pecnikar Z, Buzan EV. 20 years since the introduction of DNA barcoding: from theory to application. J Appl Genet. 2014;55(1):43–52.

  3. 3.

    Montagna M, Mereghetti V, Lencioni V, Rossaro B. Integrated taxonomy and DNA barcoding of alpine midges (Diptera: Chironomidae). PLoS One. 2016;11(3):e0149673.

  4. 4.

    Hollingsworth PM. Refining the DNA barcode for land plants. Proc Natl Acad Sci U S A. 2011;108(49):19451–2.

  5. 5.

    Vijayan K, Tsou CH. DNA barcoding in plants: taxonomy in a new perspective. Curr Sci India. 2010;99(11):1530–41.

  6. 6.

    Hebert PD, Ratnasingham S, deWaard JR. Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proceedings Biological sciences / The Royal Society. 2003;270(Suppl 1):S96–9.

  7. 7.

    Sundberg P, Kvist S, Strand M. Evaluating the utility of single-locus DNA barcoding for the identification of ribbon worms (phylum Nemertea). PLoS One. 2016;11(5):e0155541.

  8. 8.

    Vasconcelos R, Montero-Mendieta S, Simo-Riudalbas M, Sindaco R, Santos X, Fasola M, Llorente G, Razzetti E, Carranza S. Unexpectedly high levels of cryptic diversity uncovered by a complete DNA barcoding of reptiles of the Socotra archipelago. PLoS One. 2016;11(3):e0149985.

  9. 9.

    Blagoev GA, deWaard JR, Ratnasingham S, deWaard SL, Lu LQ, Robertson J, Telfer AC, Hebert PDN. Untangling taxonomy: a DNA barcode reference library for Canadian spiders. Mol Ecol Resour. 2016;16(1):325–41.

  10. 10.

    Yassin A, Amedegnato C, Cruaud C, Veuille M. Molecular taxonomy and species delimitation in Andean Schistocerca (Orthoptera: Acrididae). Mol Phylogenet Evol. 2009;53(2):404–11.

  11. 11.

    Pramual P, Adler PH. DNA barcoding of tropical black flies (Diptera: Simuliidae) of Thailand. Mol Ecol Resour. 2014;14:262–71.

  12. 12.

    Alex Smith M, Fernandez-Triana JL, Eveleigh E, Gomez J, Guclu C, Hallwachs W, Hebert PD, Hrcek J, Huber JT, Janzen D, et al. DNA barcoding and the taxonomy of Microgastrinae wasps (hymenoptera, Braconidae): impacts after 8 years and nearly 20 000 sequences. Mol Ecol Resour. 2013;13(2):168–76.

  13. 13.

    Arrigoni R, Berumen ML, Chen CA, Terraneo TI, Baird AH, Payri C, Benzoni F. Species delimitation in the reef coral genera Echinophyllia and Oxypora (Scleractinia, Lobophylliidae) with a description of two new species. Mol Phylogenet Evol. 2016;105:146–59.

  14. 14.

    Fernandez-Triana JL. Eight new species and an annotated checklist of Microgastrinae (hymenoptera, Braconidae) from Canada and Alaska. ZooKeys. 2010;63:1–53.

  15. 15.

    Fu Z, Toda MJ, Li NN, Zhang YP, Gao JJ. A new genus of anthophilous drosophilids, Impatiophila (Diptera, Drosophilidae): morphology, DNA barcoding and molecular phylogeny, with descriptions of thirty-nine new species. Zootaxa. 2016;4120(1):1–100.

  16. 16.

    Seidel M. Morphology and DNA barcoding reveal a new species of Eudicella from East Africa (Coleoptera: Scarabaeidae: Cetoniinae). Zootaxa. 2016;4137(4):535–44.

  17. 17.

    Blair C, Bryson RW Jr. Cryptic diversity and discordance in single-locus species delimitation methods within horned lizards (Phrynosomatidae: Phrynosoma). Mol Ecol Resour. 2017;17:1168–82.

  18. 18.

    Ortiz D, Francke OF. Two DNA barcodes and morphology for multi-method species delimitation in Bonnetina tarantulas (Araneae: Theraphosidae). Mol Phylogenet Evol. 2016;101:176–93.

  19. 19.

    Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, automatic barcode gap discovery for primary species delimitation. Mol Ecol. 2012;21(8):1864–77.

  20. 20.

    Ratnasingham S, Hebert PD. A DNA-based registry for all animal species: the barcode index number (BIN) system. PLoS One. 2013;8(7):e66213.

  21. 21.

    Jones M, Ghoorah A, Blaxter M. jMOTU and Taxonerator: turning DNA barcode sequences into annotated operational taxonomic units. PLoS One. 2011;6(4):e19259.

  22. 22.

    Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler AP. Sequence-based species delimitation for the DNA taxonomy of undescribed insects. Syst Biol. 2006;55(4):595–609.

  23. 23.

    Monaghan MT, Wild R, Elliot M, Fujisawa T, Balke M, Inward DJ, Lees DC, Ranaivosolo R, Eggleton P, Barraclough TG, et al. Accelerated species inventory on Madagascar using coalescent-based models of species delineation. Syst Biol. 2009;58(3):298–311.

  24. 24.

    Zhang J, Kapli P, Pavlidis P, Stamatakis A. A general species delimitation method with applications to phylogenetic placements. Bioinformatics. 2013;29(22):2869–76.

  25. 25.

    Renaud AK, Savage J, Adamowicz SJ. DNA barcoding of northern nearctic muscidae (Diptera) reveals high correspondence between morphological and molecular species limits. BMC Ecol. 2012;12:24.

  26. 26.

    Papadopoulou A, Cardoso A, Gomez-Zurita J. Diversity and diversification of Eumolpinae (Coleoptera: Chrysomelidae) in New Caledonia. Zool J Linn Soc-Lond. 2013;168(3):473–95.

  27. 27.

    Ratnasingham S, Hebert PD. Bold: the barcode of life data system ( Mol Ecol Notes 2007, 7(3):355–364.

  28. 28.

    Hausmann A, Godfray HC, Huemer P, Mutanen M, Rougerie R, van Nieukerken EJ, Ratnasingham S, Hebert PD. Genetic patterns in European geometrid moths revealed by the barcode index number (BIN) system. PLoS One. 2013;8(12):e84518.

  29. 29.

    Hamilton CA, Formanowicz DR, Bond JE. Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in north American tarantulas. PLoS One. 2011;6(10):e26207.

  30. 30.

    Puillandre N, Modica MV, Zhang Y, Sirovich L, Boisselier MC, Cruaud C, Holford M, Samadi S. Large-scale species delimitation method for hyperdiverse groups. Mol Ecol. 2012;21(11):2671–91.

  31. 31.

    Zhao L, Lin LL, Zheng ZM. DNA barcoding reveals polymorphism in the pygmy grasshopper Tetrix bolivari (Orthoptera, Tetrigidae). ZooKeys. 2016;582:111–20.

  32. 32.

    Huang J, Zhang A, Mao S, Huang Y. DNA barcoding and species boundary delimitation of selected species of Chinese Acridoidea (Orthoptera: Caelifera). PLoS One. 2013;8(12):e82400.

  33. 33.

    Lehmann AW, Devriese H, Tumbrinck J, Skejo J, Lehmann GUC, Hochkirch A. The importance of validated alpha taxonomy for phylogenetic and DNA barcoding studies: a comment on species identification of pygmy grasshoppers (Orthoptera, Tetrigidae). ZooKeys. 2017;679:139–44.

  34. 34.

    Jesúsbonilla VSD, Barrientoslozano L, Zaldívarriverón A. Sequence-based species delineation and molecular phylogenetics of the transitional Nearctic–Neotropical grasshopper genus Taeniopoda (Orthoptera, Romaleidae). Syst Biodivers. 2017;15(6):600–17.

  35. 35.

    Guo HF, Guan B, Shi FM, Zhou ZJ. DNA barcoding of genus Hexacentrus in China reveals cryptic diversity within Hexacentrus japonicus (Orthoptera, Tettigoniidae). ZooKeys. 2016;596:53–63.

  36. 36.

    Zhou ZJ, Li RL, Huang DW, Shi FM. Molecular identification supports most traditional morphological species of Ruspolia (Orthoptera: Conocephalinae). Invertebr Syst. 2012;26(5–6):451–6.

  37. 37.

    Hawlitschek O, Moriniere J, Lehmann GUC, Lehmann AW, Kropf M, Dunz A, Glaw F, Detcharoen M, Schmidt S, Hausmann A, et al. DNA barcoding of crickets, katydids and grasshoppers (Orthoptera) from Central Europe with focus on Austria, Germany and Switzerland. Mol Ecol Resour. 2017;17(5):1037–53.

  38. 38.

    Cigliano MM, Braun H, Eades DC, Otte D: Orthoptera Species File. Version 5.0/5.0. [11/25/2018]. <>.

  39. 39.

    Virgilio M, Backeljau T, Nevado B, De Meyer M. Comparative performances of DNA barcoding across insect orders. BMC bioinformatics. 2010;11:206.

  40. 40.

    Carstens BC, Pelletier TA, Reid NM, Satler JD. How to fail at species delimitation. Mol Ecol. 2013;22(17):4369–83.

  41. 41.

    Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.

  42. 42.

    Burland TG. DNASTAR's Lasergene sequence analysis software. Methods Mol Biol. 2000;132:71–91.

  43. 43.

    Zhou Z, Zhao L, Liu N, Guo H, Guan B, Di J, Shi F. Towards a higher-level Ensifera phylogeny inferred from mitogenome sequences. Mol Phylogenet Evol. 2017;108:22–33.

  44. 44.

    Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC bioinformatics. 2004;5:113.

  45. 45.

    Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20.

  46. 46.

    Floyd R, Blaxter ML: Available from: <>. 2006.

  47. 47.

    Xia XH. DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol. 2013;30(7):1720–8.

  48. 48.

    Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.

  49. 49.

    Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.

  50. 50.

    Ezard T, Fujisawa T, Barraclough TG: SPLITS: Species Limits by threshold statistics. R Package Version 1.0–11. <>. 2009.

  51. 51.

    Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.

  52. 52.

    Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Hohna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42.

  53. 53.

    Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016; 33(7):1870–4.

  54. 54.

    Wallace DL. A method for comparing two hierarchical clustering’s: comment. J Am Stat Assoc. 1983;78(383):569–76.

  55. 55.

    Severiano A, Pinto FR, Ramirez M, Carriço JA. Adjusted Wallace coefficient as a measure of congruence between typing methods. J Clin Microbiol. 2011;49(11):3997–4000.

  56. 56.

    Young RG, Abbott CL, Therriault TW, Adamowicz SJ. Barcode-based species delimitation in the marine realm: a test using Hexanauplia (Multicrustacea: Thecostraca and Copepoda). Genome / National Research Council Canada = Genome / Conseil national de recherches Canada. 2016:1–14.

  57. 57.

    Talavera G, Dinca V, Vila R. Factors affecting species delimitations with the GMYC model: insights from a butterfly survey. Methods Ecol Evol. 2013;4(12):1101–10.

  58. 58.

    Lohse K. Can mtDNA barcodes be used to delimit species? A response to Pons et al. (2006). Syst Biol. 2009;58(4):439–41.

  59. 59.

    De Queiroz K. Species concepts and species delimitation. Syst Biol. 2007;56(6):879–86.

  60. 60.

    Bergsten J, Bilton DT, Fujisawa T, Elliott M, Monaghan MT, Balke M, Hendrich L, Geijer J, Herrmann J, Foster GN, et al. The effect of geographical scale of sampling on DNA barcoding. Syst Biol. 2012;61(5):851–69.

  61. 61.

    Funk DJ, Omland KE. Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu Rev Ecol Evol S. 2003;34:397–423.

  62. 62.

    Ermakov OA, Simonov E, Surin VL, Titov SV, Brandler OV, Ivanova NV, Borisenko AV. Implications of hybridization, NUMTs, and overlooked diversity for DNA barcoding of Eurasian ground squirrels. PLoS One. 2015;10(1):e0117201.

  63. 63.

    Rohde K, Hau Y, Weyer J, Hochkirch A. Wide prevalence of hybridization in two sympatric grasshopper species may be shaped by their relative abundances. BMC Evol Biol. 2015;15:191.

  64. 64.

    Cole JA. Reinforcement and a cline in mating behaviour evolve in response to secondary contact and hybridization in shield-back katydids (Orthoptera: Tettigoniidae). J Evol Biol. 2016;29(9):1652–66.

  65. 65.

    Hochkirch A, Lemke I. Asymmetric mate choice, hybridization, and hybrid fitness in two sympatric grasshopper species. Behav Ecol Sociobiol. 2011;65(8):1637–45.

  66. 66.

    Bickford D, Lohman DJ, Sodhi NS, Ng PK, Meier R, Winker K, Ingram KK, Das I. Cryptic species as a window on diversity and conservation. Trends Ecol Evol. 2007;22(3):148–55.

  67. 67.

    Furman A, Postawa T, Oztunc T, Coraman E. Cryptic diversity of the bent-wing bat, Miniopterus schreibersii (Chiroptera: Vespertilionidae), in Asia minor. BMC Evol Biol. 2010;10:121.

  68. 68.

    Ahrens D, Fujisawa T, Krammer HJ, Eberle J, Fabrizi S, Vogler AP. Rarity and incomplete sampling in DNA-based species delimitation. Syst Biol. 2016;65(3):478–94.

  69. 69.

    TL SA, Chauveau O, Eggers L, de Souza-Chies TT. Species discrimination in Sisyrinchium (Iridaceae): assessment of DNA barcodes in a taxonomically challenging genus. Mol Ecol Resour. 2013;14(2):324–35.

  70. 70.

    Heller KG, Ingrisch S, Liu CX, Shi FM, Hemp C, Warchalowska-Sliwa E, Rentz DCF. Complex songs and cryptic ethospecies: the case of the Ducetia japonica group (Orthoptera: Tettigonioidea: Phaneropteridae: Phaneropterinae). Zool J Linn Soc-Lond. 2017;181(2):286–307.

  71. 71.

    Moulton MJ, Song H, Whiting MF. Assessing the effects of primer specificity on eliminating numt coamplification in DNA barcoding: a case study from Orthoptera (Arthropoda: Insecta). Mol Ecol Resour. 2010;10(4):615–27.

  72. 72.

    Song H, Buhay JE, Whiting MF, Crandall KA. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proc Natl Acad Sci U S A. 2008;105(36):13486–91.

  73. 73.

    Jordal BH, Kambestad M. DNA barcoding of bark and ambrosia beetles reveals excessive NUMTs and consistent east-west divergence across Palearctic forests. Mol Ecol Resour. 2014;14(1):7–17.

Download references


We would like to thank Guanglin Xie, Zhilin Chen, Ping Wang, Baojie Du, and Qiong Song for their assistance in collection of the katydid specimens.


This work financially supported by the National Natural Science Foundation of China (No. 31471985). The authors declare that the funding body has no role in the design of the study and the collection, analysis, and interpretation of data and in writing the manuscript.

Availability of data and materials

New sequences from this study are stored at BOLD systems under the project DNA Barcoding to Katydids from China (DBKC).

Author information

ZZ and SF contributed to the design, data analysis, writing and general discussion of the taxonomic part. GH, HL, CJ, and CX took lead on the sampling, general laboratory work. All authors approved the final version of the manuscript.

Correspondence to Zhijun Zhou or Fuming Shi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Additional 34 COI-5P sequences that were mined from GenBank. (XLSX 23 kb)

Additional file 2:

Table S2. The distance within-species and to its nearest neighbor (NN). (XLSX 21 kb)

Additional file 3:

Comparison of the species delimitation results of Chinese katydids based on an analysis of 390 unique COI-5P haplotypes of the DBCHL dataset. A midpoint-rooted NJ-K2P tree was implemented in MEGA 7.0. Terminals were labeled with Sequence/Process ID, Species identifications, plus BIN. * indicated a haplotype representing more than one specimen. ** indicated a haplotype shared by more than one species. On the right: summary of putative species delimitation drawn by BINs, RESL, jMOTU, ABGD, sGMYC, mGMYC and bPTP (one column per method). Black codes represented putative MOTUs defined by at least four of the seven species delimitation methods. Grey codes represented MOTUs defined by less than four of the seven species delimitation methods. Other color codes for each column represented clustering together as a single MOTU. (TIF 8640 kb)

Additional file 4:

Comparison of the species delimitation results of Chinese katydids based on an analysis of 158 unique COI-5P haplotypes of the DBMEC dataset. A midpoint-rooted NJ-K2P tree was implemented in MEGA 7.0. Terminals were labeled with Sequence/Process ID, Species identifications, plus BIN. * indicated a haplotype representing more than one specimen. ** indicated a haplotype shared by more than one species. On the right: summary of putative species delimitation drawn by BINs, RESL, jMOTU, ABGD, sGMYC, mGMYC and bPTP (one column per method). Black codes represented putative MOTUs defined by at least four of the seven species delimitation methods. Grey codes represented MOTUs defined by less than four of the seven species delimitation methods. Other color codes for each column represented clustering together as a single MOTU. (TIF 3270 kb)

Additional file 5:

Comparison of the species delimitation results of Chinese katydids based on an analysis of 530 unique COI-5P haplotypes of the DBPPM dataset. A midpoint-rooted NJ-K2P tree was implemented in MEGA 7.0. Terminals were labeled with Sequence/Process ID, Species identifications, plus BIN. * indicated a haplotype representing more than one specimen. ** indicated a haplotype shared by more than one species. On the right: summary of putative species delimitation drawn by BINs, RESL, jMOTU, ABGD, sGMYC, mGMYC and bPTP (one column per method). Black codes represented putative MOTUs defined by at least four of the seven species delimitation methods. Grey codes represented MOTUs defined by less than four of the seven species delimitation methods. Other color codes for each column represented clustering together as a single MOTU (TIF 1150 kb)

Additional file 6:

Comparison of the species delimitation results of Chinese katydids based on an analysis of 147 unique COI-5P haplotypes of the DBTB dataset. A midpoint-rooted NJ-K2P tree was implemented in MEGA 7.0. Terminals were labeled with Sequence/Process ID, Species identifications, plus BIN. * indicated a haplotype representing more than one specimen. ** indicated a haplotype shared by more than one species. On the right: summary of putative species delimitation drawn by BINs, RESL, jMOTU, ABGD, sGMYC, mGMYC and bPTP (one column per method). Black codes represented putative MOTUs defined by at least four of the seven species delimitation methods. Grey codes represented MOTUs defined by less than four of the seven species delimitation methods. Other color codes for each column represented clustering together as a single MOTU. (TIF 3020 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Zhou, Z., Guo, H., Han, L. et al. Singleton molecular species delimitation based on COI-5P barcode sequences revealed high cryptic/undescribed diversity for Chinese katydids (Orthoptera: Tettigoniidae). BMC Evol Biol 19, 79 (2019).

Download citation


  • Katydids
  • DNA barcoding
  • Species delimitation
  • Cryptic species
  • Clustering-based method, Similarity-based method