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

Genetic diversity of diamondback moth, Plutella xylostella L. (Lepidoptera: Plutellidae) populations in Türkiye

BACKGROUND

Diamondback moth (DBM), Plutella xylostella (L.) (Lepidoptera: Plutellidae), is an important worldwide pest of plants belonging to the Brassicaceae family. In this study, we investigated genetic diversity of DBM populations in Brassicaceae production areas in Türkiye using the partial mtDNA CO1 gene region.

METHODS

We determined 43 samples from 11 different populations for haplotype variations using the partial mitochondrial DNA sequences a 684 bp fragment of the CO1 gene.

RESULTS

The results indicated that, the average haplotype diversity (Hd) was determined as 0.962 and nucleotide diversity (π) was determined as 0.557%. In neutrality tests, negative values were obtained in Tajima's D and Fu' Fs tests (Fu' Fs=-0.40, Tajima's D=-0.01). Tajima's D test was not found significant (p > 0.05). Fst value among DBM population estimates ranged from 0 to 0.631. Barcode gap distance was determined as 1.6%, but the intraspecies of genetic distance were found to be 0.15%.

CONCLUSION

In conclusion, the presented study provided detailed and fundamental information about the genetic diversity of DBM populations in Türkiye. Further studies are needed to develop alternative pest management strategies for DBM populations integrating genetic approaches.

DNA barcoding of the genus Gampsocleis (Orthoptera, Tettigoniidae) from China

DNA barcoding is a useful addition to the traditional morphology-based taxonomy. A ca. 650 bp fragment of the 5' end of mitochondrial cytochrome c oxidase subunit I (hereafter COI-5P) DNA barcoding was sued as a practical tool for Gampsocleis species identification. DNA barcodes from 889 specimens belonging to 8 putative Gampsocleis species was analyzed, including 687 newly generated DNA barcodes. These barcode sequences were clustered/grouped into Operational Taxonomic Units (OTUs) using the criteria of five algorithms, namely Barcode Index Number (BIN) System, Assemble Species by Automatic Partitioning (ASAP), 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). The Taxon ID Tree grouped sequences of morphospecies and almost all MOTUs in distinct nonoverlapping clusters. Both long- and short-winged Gampsocleis species are reciprocally monophyletic in the Taxon ID Tree. In BOLD, 889 barcode sequences are assigned to 17 BINs. The algorithms ASAP, jMOTU, bPTP and GMYC clustered the barcode sequences into 6, 13, 10, and 23 MOTUs, respectively. BIN, ASAP, and bPTP algorithm placed three long-winged species, G. sedakovii, G. sinensis and G. ussuriensis within the same MOTU. All species delimitation algorithms split two short-winged species,G. fletcheri and G. gratiosa into at least two MOTUs each, except for ASAP algorithm. More detailed molecular and morphological integrative studies are required to clarify the status of these MOTUs in the future.

Molecular Species Delimitation Using COI Barcodes of Mealybugs (Hemiptera: Pseudococcidae) from Coffee Plants in Espírito Santo, Brazil

Mealybugs are insects belonging to the family Pseudococcidae. This family includes many plant-pest species with similar morphologies, which may lead to errors in mealybug identification and delimitation. In the present study, we employed molecular-species-delimitation approaches based on distance (ASAP) and coalescence (GMYC and mPTP) methods to identify mealybugs collected from coffee and other plant hosts in the states of Espírito Santo, Bahia, Minas Gerais, and Pernambuco, Brazil. We obtained 171 new COI sequences, and 565 from the BOLD Systems database, representing 26 candidate species of Pseudococcidae. The MOTUs estimated were not congruent across different methods (ASAP-25; GMYC-30; mPTP-22). Misidentifications were revealed in the sequences from the BOLD Systems database involving Phenacoccus solani × Ph. solenopsis, Ph. tucumanus × Ph. baccharidis, and Planacoccus citri × Pl. minor species. Ten mealybug species were collected from coffee plants in Espírito Santo. Due to the incorrect labeling of the species sequences, the COI barcode library of the dataset from the database needs to be carefully analyzed to avoid the misidentification of species. The systematics and taxonomy of mealybugs may be improved by integrative taxonomy which may facilitate the integrated pest management of these pests.

DNA Barcoding of Culicoides Latreille (Diptera: Ceratopogonidae) From Thailand Reveals Taxonomic Inconsistencies and Novel Diversity Among Reported Sequences

Recent focus on Culicoides species diversity in Thailand was prompted by a need to identify vectors responsible for the transmission of African Horse Sickness in that country. To assist rapid genetic identification of species, we sampled mitochondrial cytochrome c oxidase subunit I (COI) DNA barcodes (N = 78) from 40 species of Culicoides biting midge from Thailand, including 17 species for which DNA barcodes were previously unavailable. The DNA barcodes were assigned to 39 Barcode Identification Numbers (BINs) representing terminal genetic clusters at the Barcode of Life Data systems (BOLD). BINs assisted with comparisons to published conspecific DNA barcodes and allowed partial barcodes obtained from seven specimens to be associated with BINs by their similarity. Some taxonomic issues were revealed and attributed to the possible misidentification of earlier reported specimens as well as a potential synonymy of C. elbeli Wirth & Hubert and C. menglaensis Chu & Liu. Comparison with published BINs also revealed genetic evidence of divergent population processes and or potentially cryptic species in 16 described taxa, flagged by their high levels of COI sequence difference among conspecifics. We recommend the BOLD BIN system to researchers preparing DNA barcodes of vouchered species for public release. This will alert them to taxonomic incongruencies between their records and publicly released DNA barcodes, and also flag genetically deep and potentially novel diversity in described species.

A DNA barcode library for katydids, cave crickets, and leaf-rolling crickets (Tettigoniidae, Rhaphidophoridae and Gryllacrididae) from Zhejiang Province, China

Barcode libraries are generally assembled with two main objectives in mind: specimen identification and species discovery/delimitation. In this study, the standard COI barcode region was sequenced from 681 specimens belonging to katydids (Tettigoniidae), cave crickets (Rhaphidophoridae), and leaf-rolling crickets (Gryllacrididae) from Zhejiang Province, China. Of these, four COI-5P sequences were excluded from subsequent analyses because they were likely NUMTs (nuclear mitochondrial pseudogenes). The final dataset consisted of 677 barcode sequences representing 90 putative species-level taxa. Automated cluster delineation using the Barcode of Life Data System (BOLD) revealed 118 BINs (Barcodes Index Numbers). Among these 90 species-level taxa, 68 corresponded with morphospecies, while the remaining 22 were identified based on reverse taxonomy using BIN assignment. Thirteen of these morphospecies were represented by a single barcode (so-called singletons), and each of 19 morphospecies were split into more than one BIN. The consensus delimitation scheme yielded 55 Molecular Operational Taxonomic Units (MOTUs). Only four morphospecies (I _max > DNN) failed to be recovered as monophyletic clades (i.e., Elimaeaterminalis, Phyllomimusklapperichi, Sinochloraszechwanensis and Xizicushowardi), so it is speculated that these may be species complexes. Therefore, the diversity of katydids, cave crickets, and leaf-rolling crickets in Zhejiang Province is probably slightly higher than what current taxonomy would suggest.

The efficiency of the COI gene as a DNA barcode and an overview of Orthoptera (Caelifera and Ensifera) sequences in the BOLD System

Orthoptera, among the oldest and most numerous insect lineages, is an excellent model for evolutionary studies but has numerous taxonomic problems. To mitigate these issues, the cytochrome c oxidase subunit I (COI), standardized with the DNA barcode for Metazoa, is increasingly used for specimen identification and species delimitation. We tested the performance of COI as a DNA barcode in Orthoptera, using two analyses based on intra- and interspecific distances, barcode gap and Probability of Correct Identification (PCI); and estimated species richness through Automatic Barcode Gap Discovery (ABGD) and Assemble Species by Automatic Partitioning (ASAP). We filtered all sequences of Orthoptera available in Barcode of Life Data System (BOLD) and used 11,605 COI sequences, covering 1,132 species, 226 genera, and 18 families. The overall average PCI was 73.86%. For 82.2% of genera the barcode gap boxplots were classified as good or intermediate, indicating that COI can be effective as a DNA barcode in Orthoptera, although with varying efficiency depending on the need for more information. ABGD and ASAP inferred species richness similar to labels informed by BOLD for the suborders Caelifera and Ensifera. The representation of Orthoptera in the BOLD database and the results of these analyses are discussed.

Phylogeny and Integrative Taxonomy of the Genera Gymnaetoides and Pseudotachycines (Orthoptera: Rhaphidophoridae)

The genera Gymnaetoides and Pseudotachycines are endemic to China and are morphologically homogeneous. The few available diagnostic characters make species identification particularly challenging. Species cannot be classified according to the given generic diagnosis, and phylogenetic analyses have not been reported. Here, we reconstruct the phylogeny using Bayesian inference and maximum likelihood and employ four approaches to delimit species. The results suggest that both Gymnaetoides and Pseudotachycines are paraphyletic. Therefore, we revise their taxonomy based on the combination of morphological characters and molecular data. A new genus Homotachycines Zhu & Shi gen. nov. is erected, and six new combinations are proposed. Species delimitation identifies 15 new species and one new subspecies: Gymnaetoides huangshanensis, G. petalus, G. yangmingensis, G. lushanensis, Pseudotachycines procerus, P. procerus guizhouensis, P. zhengi, P. nephrus, P. sagittus, P. fengyangshanensis, Homotachycines triangulus, H. quadratus, H. baokangensis, H. fusus, H. concavus, and H. qinlingensis sp. nov. Moreover, we find that the shapes of the dorsal lateral lobes and the dorsal median lobe of the male genitalia are also important characters for identifying these genera and that the shapes of the dorsal and lateral sclerites of the male genitalia are suitable for the classifications of species.

DNA Barcoding and Phylogeny of Acari Species Based on ITS and COI Markers

Acari harbor numerous minute species of agricultural economic importance, mainly Tetranychidae and Phytoseiidae. Great efforts have been established by means of recovering morphological, molecular, and phylogenetic traits for species identification. Traditional identification still relies on external diagnostic characters, which are limited and usually exhibit large phenotypic plasticity within the species, rendering them useless for species delimitation and identification. We decided to increase the number of sequences of the Acari mitochondrial COI (Cytochrome C oxidase I) marker and ITS nuclear ribosomal DNA region for species identification in Tetranychidae and Phytoseiidae. The molecular data allow us to establish species boundaries and phylogenetic relationships among several clades of Acari, mainly Tetranychidae and Phytoseiidae. Sequence comparisons between complete COI and the Acari mitochondrial COI, ITS1-5,8S-ITS2, and ITS2 among all Acari sequences have demonstrated that the selected regions, even small, gave enough informative positions for both species’ identification and phylogenetic studies. Analyses of both DNA regions have unveiled their use as species identification characters, with special emphasis on Acari mitochondrial COI for Tetranychidae and Phytoseiidae species in comparison with the Folmer fragment, which has been universally used as a barcode marker. We demonstrated that the Acari mitochondrial COI region is also a suitable marker to establish a barcode dataset for Acari identification. Our phylogenetic analyses are congruent with other recent works, showing that Acari is a monophyletic group, of which Astigmata, Ixodida, Mesostigmata, Oribatida, and Prostigmata are also monophyletic.

A DNA barcode survey of insect biodiversity in Pakistan

Although Pakistan has rich biodiversity, many groups are poorly known, particularly insects. To address this gap, we employed DNA barcoding to survey its insect diversity. Specimens obtained through diverse collecting methods at 1,858 sites across Pakistan from 2010-2019 were examined for sequence variation in the 658 bp barcode region of the cytochrome c oxidase 1 (COI) gene. Sequences from nearly 49,000 specimens were assigned to 6,590 Barcode Index Numbers (BINs), a proxy for species, and most (88%) also possessed a representative image on the Barcode of Life Data System (BOLD). By coupling morphological inspections with barcode matches on BOLD, every BIN was assigned to an order (19) and most (99.8%) were placed to a family (362). However, just 40% of the BINs were assigned to a genus (1,375) and 21% to a species (1,364). Five orders (Coleoptera, Diptera, Hemiptera, Hymenoptera, Lepidoptera) accounted for 92% of the specimens and BINs. More than half of the BINs (59%) are so far only known from Pakistan, but others have also been reported from Bangladesh (13%), India (12%), and China (8%). Representing the first DNA barcode survey of the insect fauna in any South Asian country, this study provides the foundation for a complete inventory of the insect fauna in Pakistan while also contributing to the global DNA barcode reference library.

A practical approach to revise the Oxynoemacheilus bergianus species group (Teleostei: Nemacheilidae)

The Oxynoemacheilus bergianus species group is revised based on tree topology (ML, NJ, MP), distance (K2P and ASAP) and Poisson tree process analyses of DNA barcode data tested against morphometric and morphological characters including colour patterns. The O. bergianus species group is distinguished from other Oxynoemacheilus groups based on morphological characters: its constituent species have a slender caudal peduncle, a suborbital flap in the male, a mottled or blotched colour pattern, and lack bold, black spots on the caudal-fin base. It is also supported as a monophyletic unit in our molecular analysis.                 The O. bergianus group includes 10 molecular clades following congruently well-supported NJ, MP and ML based entities. Species described as O. bergianus, O. banarescui, O. erdali, O. fatsaensis, O. samanticus, and O. simavicus from Turkey, O. lenkoranensis from Azerbaijan, and O. longipinnis and O. parvinae from Iran belong to this species group. The group includes also four unnamed molecular clades. We were unable to detect external differences between any of the molecular clades in colour pattern or any morphometric or morphological characters examined. In the 10 molecular clades in the O. bergianus species group, the intraclade K2P distance ranges from 0.0–1.8% while the distances between molecular clades ranges from 0.6–5.9%. To resolve the species diversity of this group, we also analysed the intraspecific and interspecific variability in the K2P distance of DNA barcode data from 53 other Oxynoemacheilus species. Here, the intraspecific variability ranges from 0.0–2.4% while the interspecific K2P distance ranges from 1.2–20.8%. In the O. bergianus species group, only four groups are detected by the mPTP species delimitation approach distinguished by a K2P distance of 2.9% or more. We treat these four groups as valid species, corresponding to O. banarescui, O. bergianus, O. fatsaensis, and O. simavicus. Oxynoemacheilus samanticus from the Kızılırmak and Seyhan drainages, O. lenkoranensis from the Caspian basin, O. erdali from the Euphrates, and O. longipinnis and O. parvinae from the Tigris drainage are treated as synonyms of O. bergianus. Fishes from an unnamed molecular clade from the upper Tigris, and from a second unnamed clade from the upper Euphrates, are both identified as O. bergianus. Oxynoemacheilus bergianus might be a junior synonym of O. bergi from the Kura. The distribution range of O. simavicus, described from the Simav drainage in the Marmara basin, is expanded to the east and two molecularly differentiated population groups occur in the Sakarya drainage, the Büyük Melen River and potentially in other adjacent coastal streams. Oxynoemacheilus fatsaensis, described from the coastal stream Elekçi in northern Anatolia, is also widespread in the Yeşilırmak drainage. Morphological characters proposed to distinguish O. fatsaensis from the other species of the O. bergianus group could not be confirmed by our data on fishes from the Yeşilırmak. This study also discusses the theoretical background, our reasons for conducting this revision in the way we did, and what the alternatives would be.  

The dark side of pseudoscorpion diversity: The German Barcode of Life campaign reveals high levels of undocumented diversity in European false scorpions

DNA barcoding is particularly useful for identification and species delimitation in taxa with conserved morphology. Pseudoscorpions are arachnids with high prevalence of morphological crypsis. Here, we present the first comprehensive DNA barcode library for Central European Pseudoscorpiones, covering 70% of the German pseudoscorpion fauna (35 out of 50 species). For 21 species, we provide the first publicly available COI barcodes, including the rare Anthrenochernes stellae Lohmander, a species protected by the FFH Habitats Directive. The pattern of intraspecific COI variation and interspecific COI variation (i.e., presence of a barcode gap) generally allows application of the DNA barcoding approach, but revision of current taxonomic designations is indicated in several taxa. Sequences of 36 morphospecies were assigned to 74 BINs (barcode index numbers). This unusually high number of intraspecific BINs can be explained by the presence of overlooked cryptic species and by the accelerated substitution rate in the mitochondrial genome of pseudoscorpions, as known from previous studies. Therefore, BINs may not be an appropriate proxy for species numbers in pseudoscorpions, while partitions built with the ASAP algorithm (Assemble Species by Automatic Partitioning) correspond well with putative species. ASAP delineated 51 taxonomic units from our data, an increase of 42% compared with the present taxonomy. The Neobisium carcionoides complex, currently considered a polymorphic species, represents an outstanding example of cryptic diversity: 154 sequences from our dataset were allocated to 23 BINs and 12 ASAP units.

A dearth of data: fitting parasitoids into ecological networks

Studying parasitoids can provide insights into global diversity estimates, climate change impacts, and agroecosystem service provision. However, this potential remains largely untapped due to a lack of data on how parasitoids interact with other organisms. Ecological networks are a useful tool for studying and exploiting the impacts of parasitoids, but their construction is hindered by the magnitude of undescribed parasitoid species, a sparse knowledge of host ranges, and an under-representation of parasitoids within DNA-barcode databases (we estimate <5% have a barcode). Here, we advocate the use of DNA metabarcoding to construct the host-parasitoid component of multilayer networks. While the incorporation of parasitoids into network-based analyses has far ranging applications, we focus on its potential for assessing ecosystem service provision within agroecosystems.

Gene Flow Increases Phylogenetic Structure and Inflates Cryptic Species Estimations: A Case Study on Widespread Philippine Puddle Frogs (Occidozyga laevis)

In cryptic amphibian complexes, there is a growing trend to equate high levels of genetic structure with hidden cryptic species diversity. Typically, phylogenetic structure and distance-based approaches are used to demonstrate the distinctness of clades and justify the recognition of new cryptic species. However, this approach does not account for gene flow, spatial, and environmental processes that can obfuscate phylogenetic inference and bias species delimitation. As a case study, we sequenced genome-wide exons and introns to evince the processes that underlie the diversification of Philippine Puddle Frogs-a group that is widespread, phenotypically conserved, and exhibits high levels of geographically-based genetic structure. We showed that widely adopted tree- and distance-based approaches inferred up to 20 species, compared to genomic analyses that inferred an optimal number of five distinct genetic groups. Using a suite of clustering, admixture, and phylogenetic network analyses, we demonstrate extensive admixture among the five groups and elucidate two specific ways in which gene flow can cause overestimations of species diversity: (1) admixed populations can be inferred as distinct lineages characterized by long branches in phylograms; and (2) admixed lineages can appear to be genetically divergent, even from their parental populations when simple measures of genetic distance are used. We demonstrate that the relationship between mitochondrial and genome-wide nuclear p-distances is decoupled in admixed clades, leading to erroneous estimates of genetic distances and, consequently, species diversity. Additionally, genetic distance was also biased by spatial and environmental processes. Overall, we showed that high levels of genetic diversity in Philippine Puddle Frogs predominantly comprise metapopulation lineages that arose through complex patterns of admixture, isolation-by-distance, and isolation-by-environment as opposed to species divergence. Our findings suggest that speciation may not be the major process underlying the high levels of hidden diversity observed in many taxonomic groups and that widely-adopted tree- and distance-based methods overestimate species diversity in the presence of gene flow.

DNA barcoding of selected short-horned grasshoppers (Orthoptera: Acrididae) from Indian Himalayan region

In the context of Indian zoogeography, the DNA barcode data of short-horned grasshoppers (family Acrididae) are limited in global databases. Hence, the present study was aimed to collect selected Acridid species from the Indian Himalayan regions and generate DNA barcode data to enrich the global database. The estimated K2P genetic distances, Bayesian analysis (BA) topology and multiple species delimitation methods (ABGD, bPTP, and GMYC) clearly discriminate all the studied species. Based on high genetic distance (7.5%), multiple clades, and more than one molecular operational taxonomic unit, the present study elucidates the allopatric speciation and presence of possible cryptic diversity of Oxya japonica within India, China, and Russia. The present study suggests the collection of multiple specimens from different geographical locations and the generation of more DNA barcode data would facilitate the actual diversity of this insect group.

Assembling a DNA barcode reference library for the spiders (Arachnida: Araneae) of Pakistan

Morphological study of 1,795 spiders from sites across Pakistan placed these specimens in 27 families and 202 putative species. COI sequences >400 bp recovered from 1,782 specimens were analyzed using neighbor-joining trees, Bayesian inference, barcode gap, and Barcode Index Numbers (BINs). Specimens of 109 morphological species were assigned to 123 BINs with ten species showing BIN splits, while 93 interim species included representatives of 98 BINs. Maximum conspecific divergences ranged from 0-5.3% while congeneric distances varied from 2.8-23.2%. Excepting one species pair (Oxyopes azhari-Oxyopes oryzae), the maximum intraspecific distance was always less than the nearest-neighbor (NN) distance. Intraspecific divergence values were not significantly correlated with geographic distance. Most (75%) BINs detected in this study were new to science, while those shared with other nations mainly derived from India. The discovery of many new, potentially endemic species and the low level of BIN overlap with other nations highlight the importance of constructing regional DNA barcode reference libraries.