Large-scale species delimitation method for hyperdiverse groups

Towards holistic insect monitoring: species discovery, description, identification and traits for all insects

Holistic insect monitoring needs scalable techniques to overcome taxon biases, determine species abundances, and gather functional traits for all species. This requires that we address taxonomic impediments and the paucity of data on abundance, biomass and functional traits. We here outline how these data deficiencies could be addressed at scale. The workflow starts with large-scale barcoding (megabarcoding) of all specimens from mass samples obtained at biomonitoring sites. The barcodes are then used to group the specimens into molecular operational taxonomic units that are subsequently tested/validated as species with a second data source (e.g. morphology). New species are described using barcodes, images and short diagnoses, and abundance data are collected for both new and described species. The specimen images used for species discovery then become the raw material for training artificial intelligence identification algorithms and collecting trait data such as body size, biomass and feeding modes. Additional trait data can be obtained from vouchers by using genomic tools developed by molecular ecologists. Applying this pipeline to a few samples per site will lead to greatly improved insect monitoring regardless of whether the species composition of a sample is determined with images, metabarcoding or megabarcoding. This article is part of the theme issue 'Towards a toolkit for global insect biodiversity monitoring'.

Hapsolutely: a user-friendly tool integrating haplotype phasing, network construction, and haploweb calculation

Motivation

Haplotype networks are a routine approach to visualize relationships among alleles. Such visual analysis of single-locus data is still of importance, especially in species diagnosis and delimitation, where a limited amount of sequence data usually are available and sufficient, along with other datasets in the framework of integrative taxonomy. In diploid organisms, this often requires separating (phasing) sequences with heterozygotic positions, and typically separate programs are required for phasing, reformatting of input files, and haplotype network construction. We therefore developed Hapsolutely, a user-friendly program with an ergonomic graphical user interface that integrates haplotype phasing from single-locus sequences with five approaches for network/genealogy reconstruction.

Results

Among the novel options implemented, Hapsolutely integrates phasing and graphical reconstruction steps of haplotype networks, supports input of species partition data in the common SPART and SPART-XML formats, and calculates and visualizes haplowebs and fields for recombination, thus allowing graphical comparison of allele distribution and allele sharing among subsets for the purpose of species delimitation. The new tool has been specifically developed with a focus on the workflow in alpha-taxonomy, where exploring fields for recombination across alternative species partitions may help species delimitation.

Availability and implementation

Hapsolutely is written in Python, and integrates code from Phase, SeqPHASE, and PopART in C++ and Haxe. Compiled stand-alone executables for MS Windows and Mac OS along with a detailed manual can be downloaded from https://www.itaxotools.org; the source code is openly available on GitHub (https://github.com/iTaxoTools/Hapsolutely).

Coupling DNA barcodes and exon-capture to resolve the phylogeny of Turridae (Gastropoda, Conoidea)

Taxon sampling in most phylogenomic studies is often based on known taxa and/or morphospecies, thus ignoring undescribed diversity and/or cryptic lineages. The family Turridae is a group of venomous snails within the hyperdiverse superfamily Conoidea that includes many undescribed and cryptic species. Therefore 'traditional' taxon sampling could constitute a strong risk of undersampling or oversampling Turridae lineages. To minimize potential biases, we establish a robust sampling strategy, from species delimitation to phylogenomics. More than 3,000 cox-1 "barcode" sequences were used to propose 201 primary species hypotheses, nearly half of them corresponding to species potentially new to science, including several cryptic species. A 110-taxa exon-capture tree, including species representatives of the diversity uncovered with the cox-1 dataset, was build using up to 4,178 loci. Our results show the polyphyly of the genus Gemmula, that is split into up to 10 separate lineages, of which half would not have been detected if the sampling strategy was based only on described species. Our results strongly suggest that the use of blind, exploratory and intensive barcode sampling is necessary to avoid sampling biases in phylogenomic studies.

Temporal Distribution Patterns of Cryptic Brachionus calyciflorus (Rotifera) Species in Relation to Biogeographical Gradient Associated with Latitude

Sympatric distribution and temporal overlap of cryptic zooplankton species pose a challenge to the framework of the niche differentiation theory and the mechanisms allowing competitor coexistence. We applied the methods of phylogenetic analysis, DNA taxonomy, and statistical analysis to study the temporal distribution patterns of the cryptic B. calyciflorus species, an excellent model, in three lakes, and to explore the putative mechanisms for their seasonal succession and temporal overlap. The results showed that in the warm-temperate Lake Yunlong, B. fernandoi and B. calyciflorus s.s. underwent a seasonal succession, which was largely attributed to their differential adaptation to water temperature. In the subtropical Lake Jinghu, B. fernandoi, B. calyciflorus s.s., and B. dorcas exhibited both seasonal succession and temporal overlap. Seasonal successions were largely attributed to their differential adaptation to temperature, and temporal overlap resulted from their differential responses to algal food concentration. In the tropical Lake Jinniu, B. calyciflorus s.s. persisted throughout the year and overlapped with B. dorcas for 5 months. The temporal overlap resulted from their differential responses to copepod predation. These results indicated that the temporal distribution pattern of the cryptic B. calyciforus species and the mechanism that allows competitor coexistence vary with different climate zones.

DNA Barcoding in Species Delimitation: From Genetic Distances to Integrative Taxonomy

Over the past two decades, DNA barcoding has become the most popular exploration approach in molecular taxonomy, whether for identification, discovery, delimitation, or description of species. The present contribution focuses on the utility of DNA barcoding for taxonomic research activities related to species delimitation, emphasizing the following aspects:(1) To what extent DNA barcoding can be a valuable ally for fundamental taxonomic research, (2) its methodological and theoretical limitations, (3) the conceptual background and practical use of pairwise distances between DNA barcode sequences in taxonomy, and (4) the different ways in which DNA barcoding can be combined with complementary means of investigation within a broader integrative framework. In this chapter, we recall and discuss the key conceptual advances that have led to the so-called renaissance of taxonomy, elaborate a detailed glossary for the terms specific to this discipline (see Glossary in Chap. 35 ), and propose a newly designed step-by-step species delimitation protocol starting from DNA barcode data that includes steps from the preliminary elaboration of an optimal sampling strategy to the final decision-making process which potentially leads to nomenclatural changes.

A Glossary of DNA Barcoding Terms

This chapter provides a reference glossary for the protocols in this volume. We have chosen only the very basic terms in the DNA barcode lexicon to include, and provide clear and concise definitions of these terms. We hope the reader finds this glossary useful.

The genetic structuring in pollinating wasps of Ficus hispida in continental Asia

The interaction between figs and fig wasps provides a striking example of obligate brood site pollination mutualism. Monoecious figs, constituting independent radiations in each tropical biome, are present in significant proportions worldwide, but in continental Asia, dioecious figs have diverged into various niches, making the region's assemblage remarkably diverse. However, the reproductive success of figs and fig wasps largely depends on the fig wasp dispersal process. Monoecious fig pollinators in continental Asian tropical rain forests exhibit high gene flow of the plant, while many dioecious fig pollinators have a more restricted gene flow. However, there are limited studies on the genetic structure of dioecious Ficus pollinators that pollinate figs with intermediate gene flow. Here, we used molecular methods to investigate the genetic structure of pollinating wasps of the widely distributed dioecious Ficus hispida in China and Southeast Asia. Sequence data from two gene regions were used: the mitochondrial protein-coding gene cytochrome c oxidase subunit I (COI) and the nuclear 28S genes. Both molecular and morphological results support two fig wasp species at our sampling sites. Our findings suggest that for widely sympatric Ficus species in continental Asia, monoecious figs presenting long gene glow have the fewest fig wasp species, followed by dioecious figs presenting intermediate gene flow, and dioecious figs presenting local gene flow have the most fig wasp species.

Penicillium rhizophilum, a novel species in the section Exilicaulis isolated from the rhizosphere of sugarcane in Southwest Iran

During a survey of species diversity of Penicillium and Talaromyces in sugarcane (Saccharum officinarum) rhizosphere in the Khuzestan province of Iran [1], 195 strains were examined, from which 187 belonged to Penicillium (11 species) and eight to Talaromyces (one species). In the present study, three strains of Penicillium belonging to section Exilicaulis series Restricta, identified as P. restrictum by Ansari et al. [1], were subjected to a phylogenetic study. The multilocus phylogeny of partial β-tubulin, calmodulin and RNA polymerase II second largest subunit genes enabled the recognition of one new phylogenetic species that is here formally described as Penicillium rhizophilum sp. nov. This species is phylogenetically distinct in series Restricta, but it does not show significant morphological differences from other species previously classified in the series. Therefore, we here placed bias on the phylogenetic species concept. The holotype of Penicillium rhizophilum sp. nov. is IRAN 18169F and the ex-type culture is LA30T (=IRAN 4042CT=CBS 149737T).

Abundant, diverse, unknown: Extreme species richness and turnover despite drastic undersampling in two closely placed tropical Malaise traps

Arthropods account for a large proportion of animal biomass and diversity in terrestrial systems, making them crucial organisms in our environments. However, still too little is known about the highly abundant and megadiverse groups that often make up the bulk of collected samples, especially in the tropics. With molecular identification techniques ever more evolving, analysis of arthropod communities has accelerated. In our study, which was conducted within the Global Malaise trap Program (GMP) framework, we operated two closely placed Malaise traps in Padang, Sumatra, for three months. We analyzed the samples by DNA barcoding and sequenced a total of more than 70,000 insect specimens. For sequence clustering, we applied three different delimitation techniques, namely RESL, ASAP, and SpeciesIdentifier, which gave similar results. Despite our (very) limited sampling in time and space, our efforts recovered more than 10,000 BINs, of which the majority are associated with "dark taxa". Further analysis indicates a drastic undersampling of both sampling sites, meaning that the true arthropod diversity at our sampling sites is even higher. Regardless of the close proximity of both Malaise traps (< 360 m), we discovered significantly distinct communities.

Evaluating the taxonomic status of the large sized Tricolia Risso, 1826 in the Northeast Atlantic and Mediterranean Sea

Despite a long history of taxonomic studies on the genus Tricolia Risso, 1826, there is a shortfall on thorough systematic molecular reviews of the taxon from the NE Atlantic and Mediterranean coasts. Aiming to assess the genetic distinctness among morphospecies and the taxonomic status of currently accepted large sized species in these areas, we conducted a molecular phylogenetic analysis of the genus based on one mitochondrial (cox1) and two nuclear (28S and ITS2) markers. Seven Tricolia species were consistently retrieved in the analyses, including a new genetic lineage in the NE Atlantic designated as Tricolia sp. 1. Molecular analyses revealed that only one species, T. azorica, occurs in the NE Atlantic archipelagos. The sister taxa T. pullus (Mediterranean) and T. picta (NE Atlantic) should be classified as distinct species, instead of subspecies of the T. pullus group (sensu Gofas 1982). Tricolia miniata is also a complex of species in the Mediterranean and future studies across the distribution range are necessary to clarify its status.

SNP data reveals the complex and diverse evolutionary history of the blue-ringed octopus genus (Octopodidae: Hapalochlaena) in the Asia-Pacific

The blue-ringed octopus species complex (Hapalochlaena spp.), known to occur from Southern Australia to Japan, currently contains four formally described species (Hapalochlaena maculosa, Hapalochlaena fasciata, Hapalochlaena lunulata and Hapalochlaena nierstraszi). These species are distinguished based on morphological characters (iridescent blue rings and/or lines) along with reproductive strategies. However, the observation of greater morphological diversity than previously captured by the current taxonomic framework indicates that a revision is required. To examine species boundaries within the genus we used mitochondrial (12S rRNA, 16S rRNA, cytochrome c oxidase subunit 1 [COI], cytochrome c oxidase subunit 3 [COIII] and cytochrome b [Cytb]) and genome-wide SNP data (DaRT seq) from specimens collected across its geographic range including variations in depth from 3m to >100m. This investigation indicates substantially greater species diversity present within the genus Hapalochlaena than is currently described. We identified 10,346 SNPs across all locations, which when analysed support a minimum of 11 distinct clades. Bayesian phylogenetic analysis of the mitochondrial COI gene on a more limited sample set dates the diversification of the genus to ∼30mya and corroborates eight of the lineages indicated by the SNP analyses. Furthermore, we demonstrate that the diagnostic lined patterning of H. fasciata found in North Pacific waters and NSW, Australia is polyphyletic and therefore likely the result of convergent evolution. Several "deep water" (> 100m) lineages were also identified in this study with genetic convergence likely to be driven by external selective pressures. Examination of morphological traits, currently being undertaken in a parallel morphological study, is required to describe additional species within the complex.

DNA barcoding of the genus Alburnoides Jeitteles, 1861 (Actinopterygii, Cyprinidae) from Anatolia, Turkey

The present study investigated the ability of DNA barcoding to reliably identify the endemic freshwater species in Turkey, known as biodiversity hotspots. The barcode region (652 bp) of the mitochondrial cytochrome c oxidase subunit I (COI) was used to barcode 153 individuals from 13 morphologically identified species of the genus Alburnoides. Based on the Kimura two-parameter (K2P) evolution model, the average interspecific distance (0.0595) was 31-fold higher than the average intraspecific distance (0.0019). There was a clear-cut barcode gap (0.0158–0.0187) between maximum intraspecific distance (A. tzanevi and A. velioglui) and minimum nearest-neighbour distance (A. freyhofi and A. kurui) for Anatolian Alburnoides species and a common genetic threshold of 0.0158 sequence divergence was defined for species delimitation. The multiple species delimitation methods (ABGD, ASAP, GMYC and bPTP) revealed a total of 11 molecular operational taxonomic units (MOTUs) for 13 morphospecies. Neighbour-joining (NJ), Maximum Likelihood (ML) and Bayesian Inference (BI) tree analysis indicated that all haplotypes were clustered into two major clades, which corresponded to eleven Alburnoides species clusters, with strong bootstrap support. Furthermore, all the specimens clustered in concurrence with the morpho-taxonomic status of the species, except for two species (A. coskuncelebii and A. emineae) that were morphologically differentiated, but showed overlap in variation for COI-based DNA barcode data with other species. Overall, present results identified that COI-based DNA barcoding is effective for species identification and cataloguing of genus Alburnoides in Turkey.

Mediterranean Matters: Revision of the Family Onchidorididae (Mollusca, Nudibranchia) with the Description of a New Genus and a New Species

The Mediterranean Sea hosts a great Nudibranchia diversity and has proved to be particularly intriguing in the case of the family Onchidorididae, a group of dorid nudibranchs that lately increased its diversity with the addition of one recently described Mediterranean species. The Onchidorididae family has a troubled systematic history to date, characterized by uncertainties and genera that are considered valid or not, according to the different authors. This confused taxonomy reflects the lack of a broad and comprehensive view on the phylogenetic relationships occurring between Onchidorididae members, an incorrect interpretation of the diagnostic morphological characters, and a poor knowledge of important biological aspects characterizing the different genera included in the family. To shed some light on the systematics of Onchidorididae, an integrative taxonomic revision was carried out involving morphological, ecological, and molecular analyses on an updated dataset. Mediterranean specimens and species were added to the dataset of the already known Onchidorididae and a new species from the Adriatic Sea (Central Mediterranean Sea) is described here. Furthermore, historical controversies are clarified due to the discovery of new important synapomorphies useful to define genera belonging to the Onchidorididae family and to describe a new genus. Finally, the taxonomic status of all the known Onchidorididae species is investigated and discussed, filling the gap of knowledge on neglected species.

Multilocus phylogeny and species delimitation suggest synonymies of two Lucanus Scopoli, 1763 (Coleoptera, Lucanidae) species names

Phylogenetic relationsships of four nominal Lucanus Scopoli, 1763 species, L.swinhoei Parry, 1874, L.continentalis Zilioli, 1998, L.liuyei Huang & Chen, 2010, and L.wuyishanensis Schenk, 1999, are assessed based on mitochondrial (16S rDNA, COI) and nuclear (28S rDNA, Wingless) genes. The genetic distance is 0.0072 between L.swinhoei and L.continentalis, and 0.0094 between L.wuyishanensis and L.liuyei. Three species-delimitation approaches (ABGD, PTP, and GMYC) consistently showed L.swinhoei + L.continentalis and L.wuyishanensis + L.liuyei as two MOTUs. A new synonymy, L.liuyei = L.wuyishanensis, is proposed. Synonymy of L.swinhoei over L.continentalis is confirmed.

Neither slugs nor snails: a molecular reappraisal of the gastropod family Velutinidae

The systematics of the marine mollusc family Velutinidae has long been neglected by taxonomists, mainly because their often internal and fragile shells offer no morphological characters. Velutinids are usually undersampled owing to their cryptic mantle coloration on the solitary, social or colonial ascidians on which they feed and lay eggs. In this study, we address the worldwide diversity and phylogeny of Velutinidae based on the largest molecular dataset (313 specimens) to date, accounting for &gt; 50% of the currently accepted genera, coupled with morphological and ecological data. Velutinids emerge as a diverse group, encompassing four independent subfamily-level lineages, two of which are newly described herein: Marseniopsinae subfam. nov. and Hainotinae subfam. nov. High diversity was found at genus and species levels, with two newly described genera (Variolipallium gen. nov. and Pacifica gen. nov.) and ≥ 86 species in the assayed dataset, 58 of which are new to science (67%). Velutinidae show a remarkable morphological plasticity in shell morphology, mantle extension and chromatic patterns. This variability is likely to be the result of different selective forces, including habitat, depth and trophic interactions.

DNA barcoding and species delimitation of the genus Oxynoemacheilus (Teleostei: Nemacheilidae) in Anatolia

The DNA barcoding approach was used for the determination of evolutionary relationships and species delimitation of the genus Oxynoemacheilus (Teleostei: Nemacheilidae). The COI barcode region (615 bp amplicon) was used to barcode 444 individuals from 64 morphologically identified species in the genus Oxynoemacheilus and 189 haplotypes were identified. The average of the interspecific p distance (9.59%) was about 21-fold higher than the average intraspecific distance (0.44%). A general genetic threshold of 1.46% sequence divergence was defined for species delimitation. The multiple species delimitation methods (BCM, GMYC, bPTP and TCS) revealed a total of 62 molecular operational taxonomic units for 64 morphospecies with a new loach species from the BuyukMelen River. Neighbour-joining, maximum likelihood and Bayesian inference analyses indicated that all haplotypes were clustered into 62 clades, which corresponded to Oxynoemacheilus species, with strong bootstrap support (≥95%). Furthermore, all samples grouped in concurrence with the taxonomic status of the species except for species groups (O. germencicus-O. cinicus-O. mesudae and O. leontinae-O. namiri) that were showed intraspecific overlap in genetic diversity for COI-based barcodes. In conclusion, our analyses indicate that COI-based barcodes provide reliable species discrimination. Therefore, we currently recommend COI barcodes as the suitable barcode for genus Oxynoemacheilus.

Towards Large-scale Integrative Taxonomy (LIT): resolving the data conundrum for dark taxa

New, rapid, accurate, scalable, and cost-effective species discovery and delimitation methods are needed for tackling “dark taxa,” here defined as groups for which \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$<$\end{document}10\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} of all species are described and the estimated diversity exceeds 1,000 species. Species delimitation for these taxa should be based on multiple data sources (“integrative taxonomy”) but collecting multiple types of data risks impeding a discovery process that is already too slow. We here develop large-scale integrative taxonomy (LIT), an explicit method where preliminary species hypotheses are generated based on inexpensive data that can be obtained quickly and cost-effectively. These hypotheses are then evaluated based on a more expensive type of “validation data” that is only obtained for specimens selected based on objective criteria applied to the preliminary species hypotheses. We here use this approach to sort 18,000 scuttle flies (Diptera: Phoridae) into 315 preliminary species hypotheses based on next-generation sequencing barcode (313 bp) clusters (using objective clustering [OC] with a 3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} threshold). These clusters are then evaluated with morphology as the validation data. We develop quantitative indicators for predicting which barcode clusters are likely to be incongruent with morphospecies by randomly selecting 100 clusters for in-depth validation with morphology. A linear model demonstrates that the best predictors for incongruence between barcode clusters and morphology are maximum p-distance within the cluster and a newly proposed index that measures cluster stability across different clustering thresholds. A test of these indicators using the 215 remaining clusters reveals that these predictors correctly identify all clusters that are incongruent with morphology. In our study, all morphospecies are true or disjoint subsets of the initial barcode clusters so that all incongruence can be eliminated by varying clustering thresholds. This leads to a discussion of when a third data source is needed to resolve incongruent grouping statements. The morphological validation step in our study involved 1,039 specimens (5.8\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} of the total). The formal LIT protocol we propose would only have required the study of 915 (5.1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}: 2.5 specimens per species), as we show that clusters without signatures of incongruence can be validated by only studying two specimens representing the most divergent haplotypes. To test the generality of our results across different barcode clustering techniques, we establish that the levels of incongruence are similar across OC, Automatic Barcode Gap Discovery (ABGD), Poisson Tree Processes (PTP), and Refined Single Linkage (RESL) (used by Barcode of Life Data System to assign Barcode Index Numbers [BINs]). OC and ABGD achieved a maximum congruence score with the morphology of 89\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} while PTP was slightly less effective (84\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document}). RESL could only be tested for a subset of the specimens because the algorithm is not public. BINs based on 277 of the original 1,714 haplotypes were 86\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} congruent with morphology while the values were 89\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for OC, 74\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for PTP, and 72\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\%$\end{document} for ABGD. [Biodiversity discovery; dark taxa; DNA barcodes; integrative taxonomy.]

DNA barcoding for biodiversity assessment: Croatian stoneflies (Insecta: Plecoptera)

Background

The hemi-metabolous aquatic order Plecoptera (stoneflies) constitutes an indispensable part of terrestrial and aquatic food webs due to their specific life cycle and habitat requirements. Stoneflies are considered one of the most sensitive groups to environmental changes in freshwater ecosystems and anthropogenic changes have caused range contraction of many species. Given the critical threat to stoneflies, the study of their distribution, morphological variability and genetic diversity should be one of the priorities in conservation biology. However, some aspects about stoneflies, especially a fully resolved phylogeny and their patterns of distribution are not well known. A study that includes comprehensive field research and combines morphological and molecular identification of stoneflies has not been conducted in Croatia so far. Thus, the major aim of this study was to regenerate a comprehensive and taxonomically well-curated DNA barcode database for Croatian stoneflies, to highlight the morphological variability obtained for several species and to elucidate results in light of recent taxonomy.

Methods

A morphological examination of adult specimens was made using basic characteristics for distinguishing species: terminalia in males and females, head and pronotum patterns, penial morphology, and egg structures. DNA barcoding was applied to many specimens to help circumscribe known species, identify cryptic or yet undescribed species, and to construct a preliminary phylogeny for Croatian stoneflies.

Results

Sequences (658 bp in length) of 74 morphospecies from all families present in Croatia were recovered from 87% of the analysed specimens (355 of 410), with one partial sequence of 605 bp in length for Capnopsis schilleri balcanica Zwick, 1984. A total of 84% morphological species could be unambiguously identified using COI sequences. Species delineation methods confirmed the existence of five deeply divergent genetic lineages, with monophyletic origin, which also differ morphologically from their congeners and represent distinct entities. BIN (Barcode Index Number) assignment and species delineation methods clustered COI sequences into different numbers of operational taxonomic units (OTUs). ASAP delimited 76 putative species and achieved a maximum match score with morphology (97%). ABGD resulted in 62 and mPTP in 61 OTUs, indicating a more conservative approach. Most BINs were congruent with traditionally recognized species. Deep intraspecific genetic divergences in some clades highlighted the need for taxonomic revision in several species-complexes and species-groups. Research has yielded the first molecular characterization of nine species, with most having restricted distributions and confirmed the existence of several species which had been declared extinct regionally.

Molecular phylogeny of the genus Chondrina (Gastropoda, Panpulmonata, Chondrinidae) in the Iberian Peninsula

Chondrina Reichenbach, 1828 is a highly diverse genus of terrestrial molluscs currently including 44 species with about 28 subspecific taxa. It is distributed through North Africa, central and southern Europe, from Portugal in the West to the Caucasus and Asia Minor in the East. Approximately 70% of the species are endemic to the Iberian Peninsula constituting its main center of speciation with 34 species. This genus includes many microendemic taxa, some of them not yet described, confined to limestone habitats (being strictly rock-dwelling species). They are distributed on rocky outcrops up to 2000 m.a.s.l. It is a genus of conical-fusiform snails that differ mainly in shell characters and in the number and position of teeth in their aperture. So far, molecular studies on Chondrina have been based exclusively on the mitochondrial Cytochrome Oxidase subunit I region (COI). These studies gave a first view of the phylogeny of the genus but many inner nodes were not statistically supported. The main objective of the study is to obtain a better understanding of the phylogeny and systematics of the genus Chondrina on the Iberian Peninsula, using multilocus molecular analysis. Partial sequences of the COI and 16S rRNA genes, as well as of the nuclear Internal Transcribed Spacer 1 (ITS1-5.8S) and Internal Transcribed Spacer 2 (5.8S-ITS2-28S) were obtained from individuals of all the extant Chondrina species known from the Iberian Peninsula. In addition to this, the newly obtained COI sequences were combined with those previously published in the GenBank. Phylogenetic relationships were inferred using maximum likelihood and Bayesian methods. The reconstructed phylogenies showed high values of support for more recent branches and basal nodes. Moreover, molecular species delimitation allowed to better definethe studied species and check the presence of new taxa.

Genome Skimming Contributes to Clarifying Species Limits in Paris Section Axiparis (Melanthiaceae)

Paris L. section Axiparis H. Li (Melanthiaceae) is a taxonomically perplexing taxon with considerable confusion regarding species delimitation. Based on the analyses of morphology and geographic distribution of each species currently recognized in the taxon, we propose a revision scheme that reduces the number of species in P. sect. Axiparis from nine to two. To verify this taxonomic proposal, we employed a genome skimming approach to recover the plastid genomes (plastomes) and nuclear ribosomal DNA (nrDNA) regions of 51 individual plants across the nine described species of P. sect. Axiparis by sampling multiple accessions per species. The species boundaries within P. sect. Axiparis were explored using phylogenetic inference and three different sequence-based species delimitation methods (ABGD, mPTP, and SDP). The mutually reinforcing results indicate that there are two species-level taxonomic units in P. sect. Axiparis (Paris forrestii s.l. and P. vaniotii s.l.) that exhibit morphological uniqueness, non-overlapping distribution, genetic distinctiveness, and potential reproductive isolation, providing strong support to the proposed species delimitation scheme. This study confirms that previous morphology-based taxonomy overemphasized intraspecific and minor morphological differences to delineate species boundaries, therefore resulting in an overestimation of the true species diversity of P. sect. Axiparis. The findings clarify species limits and will facilitate robust taxonomic revision in P. sect. Axiparis.

Identifying the Genetic Distance Threshold for Entiminae (Coleoptera: Curculionidae) Species Delimitation via COI Barcodes

The subfamily Entiminae is the largest group in the family Curculionidae, and it has long represented a challenge in traditional and molecular classification. Here, we analyzed intra- and interspecific genetic distances of 621 public COI barcode sequences (658bp) from 39 genera and 110 species of Entiminae, to determine parameters most congruent in retaining established species. We found that the mean intraspecific genetic distance (3.07%) was much smaller than the mean interspecific one (21.96%), but there is a wide range of overlap between intra- and interspecific genetic distances (0.77−18.01%), indicating that there is no consistent, universal barcoding gap. Specifically, DNA barcoding gap analysis for morphospecies revealed that 102 of 110 morphospecies had barcoding gaps, and 9.18% was the optimum threshold of genetic distances for 97 species delimitation. We further confirmed this threshold with barcodes from 27 morphologically identified specimens (including 21 newly reported barcodes) sequenced from five genera and seven species. We also identified thresholds to delimit congeneric species within 14 selected genera (species > 2), which varied from 7.42% (Trichalophus) to 13.48% (Barypeithes). We herein present optimal parameters for species identification in the Entiminae. Our study suggests that despite no universal genetic distance threshold value in subfamily Entiminae, 9.18% is optimal for most species. We recommend a wider sampling of geographic populations to better account for intraspecific distance variation, and that genetic distance thresholds for species delimitation should be refined at the genus level.

Re-evaluating the case for poecilogony in the gastropod Planaxis sulcatus (Cerithioidea, Planaxidae)

Background

Planaxis sulcatus has been touted as a textbook example of poecilogony, with members of this wide-ranging Indo-Pacific marine gastropod said to produce free-swimming veligers as well as brooded juveniles. A recent paper by Wiggering et al. (BMC Evol Biol 20:76, 2020) assessed a mitochondrial gene phylogeny based on partial COI and 16S rRNA sequences for 31 individuals supplemented by observations from the brood pouch of 64 mostly unsequenced individuals. ABGD and bGYMC supported three reciprocally monophyletic clades, with two distributed in the Indo-Pacific, and one restricted to the northern Indian Ocean and Red Sea. Given an apparent lack of correlation between clade membership and morphological differentiation or mode of development, the reported 3.08% maximum K2P model-corrected genetic divergence in COI among all specimens was concluded to represent population structuring. Hence, the hypothesis that phylogenetic structure is evidence of cryptic species was rejected and P. sulcatus was concluded to represent a case of geographic poecilogony.

Results

Our goal was to reassess the case for poecilogony in Planaxis sulcatus with a larger molecular dataset and expanded geographic coverage. We sequenced an additional 55 individuals and included published and unpublished sequence data from other sources, including from Wiggering et al. Our dataset comprised 108 individuals (88 COI, 81 16S rRNA) and included nine countries unrepresented in the previous study. The expanded molecular dataset yielded a maximum K2P model-corrected genetic divergence among all sequenced specimens of 12.09%. The value of 3.08% erroneously reported by Wiggering et al. is the prior maximal distance value that yields a single-species partition in ABGD, and not the maximum K2P intraspecific divergence that can be calculated for the dataset. The bGMYC analysis recognized between two and six subdivisions, while the best-scoring ASAP partitions recognized two, four, or five subdivisions, not all of which were robustly supported in Bayesian and maximum likelihood phylogenetic analyses of the concatenated and single gene datasets. These hypotheses yielded maximum intra-clade genetic distances in COI of 2.56–6.19%, which are more consistent with hypothesized species-level thresholds for marine caenogastropods.

Conclusions

Based on our analyses of a more comprehensive dataset, we conclude that the evidence marshalled by Wiggering et al. in support of Planaxis sulcatus comprising a single widespread, highly variable species with geographic poecilogony is unconvincing and requires further investigation in an integrative taxonomic framework.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01961-7.

Does clustering of DNA barcodes agree with botanical classification directly at high taxonomic levels? Trees in French Guiana as a case study

Characterising biodiversity is one of the main challenges for the coming decades. Most diversity has not been morphologically described and barcoding is now complementing morphological-based taxonomy to further develop inventories. Both approaches have been cross-validated at the level of species and OTUs. However, many known species are not listed in reference databases. One path to speed up inventories using barcoding is to directly identify individuals at coarser taxonomic levels. We therefore studied in barcoding of plants whether morphological-based and molecular-based approaches are in agreement at genus, family and order levels. We used Agglomerative Hierarchical Clustering (with Ward, Complete and Single Linkage) and Stochastic Block Models (SBM), with two dissimilarity measures (Smith-Waterman scores, kmers). The agreement between morphological-based and molecular-based classifications ranges in most of the cases from good to very good at taxonomic levels above species, even though it decreases when taxonomic levels increase, or when using the tetramer-based distance. Agreement is correlated with the entropy of morphological-based classification and with the ratio of the mean within- and mean between-groups dissimilarities. The Ward method globally leads to the best agreement whereas Single Linkage can show poor behaviours. SBM provides a useful tool to test whether or not the dissimilarities are structured by the botanical groups. These results suggest that automatic clustering and group identification at taxonomic levels above species are possible in barcoding.

Diversity and taxonomic revision of tribes Rhipileae and Rhipiliopsideae (Halimedaceae, Chlorophyta) based on molecular and morphological data

Genera and species of the tribes Rhipileae and Rhipiliopsideae are abundant in most coral reef ecosystems worldwide. However, the group has been largely overlooked, and very little genetic data is available to accurately assess its diversity, phylogenetic relationships, and geographic distribution. Our study provided an in-depth reassessment of tribes Rhipileae and Rhipiliopsideae based on a species-rich dataset and the combination of molecular species delimitation, multilocus phylogenetic analyses (tufA, rbcL, and 18S rDNA), and morpho-anatomic observations. Our results revealed an unexpected diversity of 38 morphologically validated species hypotheses, including 20 new species, two of which are described in this paper and one resurrected species (Rhipilia diaphana). Based on our phylogenetic results we proposed to redefine the genera Rhipilia and Rhipiliopsis and described two new genera, Kraftalia gen. nov. (Rhipileae) and Rhipiliospina gen. nov. (Rhipiliopsideae). Finally, we validated Rhipiliella Kraft and included it in the tribe Rhipileae. Although Rhipilia and Rhipiliopsis have a pantropical distribution, none of the species studied here appeared cosmopolitan; instead, they have restricted distributions.

ONTbarcoder and MinION barcodes aid biodiversity discovery and identification by everyone, for everyone

BACKGROUND

DNA barcodes are a useful tool for discovering, understanding, and monitoring biodiversity which are critical tasks at a time of rapid biodiversity loss. However, widespread adoption of barcodes requires cost-effective and simple barcoding methods. We here present a workflow that satisfies these conditions. It was developed via "innovation through subtraction" and thus requires minimal lab equipment, can be learned within days, reduces the barcode sequencing cost to < 10 cents, and allows fast turnaround from specimen to sequence by using the portable MinION sequencer.

RESULTS

We describe how tagged amplicons can be obtained and sequenced with the real-time MinION sequencer in many settings (field stations, biodiversity labs, citizen science labs, schools). We also provide amplicon coverage recommendations that are based on several runs of the latest generation of MinION flow cells ("R10.3") which suggest that each run can generate barcodes for > 10,000 specimens. Next, we present a novel software, ONTbarcoder, which overcomes the bioinformatics challenges posed by MinION reads. The software is compatible with Windows 10, Macintosh, and Linux, has a graphical user interface (GUI), and can generate thousands of barcodes on a standard laptop within hours based on only two input files (FASTQ, demultiplexing file). We document that MinION barcodes are virtually identical to Sanger and Illumina barcodes for the same specimens (> 99.99%) and provide evidence that MinION flow cells and reads have improved rapidly since 2018.

CONCLUSIONS

We propose that barcoding with MinION is the way forward for government agencies, universities, museums, and schools because it combines low consumable and capital cost with scalability. Small projects can use the flow cell dongle ("Flongle") while large projects can rely on MinION flow cells that can be stopped and re-used after collecting sufficient data for a given project.

Solving the Coral Species Delimitation Conundrum

Distinguishing coral species is not only crucial for physiological, ecological and evolutionary studies, but also to enable effective management of threatened reef ecosystems. However, traditional hypotheses that delineate coral species based on morphological traits from the coral skeleton are frequently at odds with tree-based molecular approaches. Additionally, a dearth of species-level molecular markers has made species delimitation particularly challenging in species-rich coral genera, leading to the widespread assumption that inter-specific hybridization might be responsible for this apparent conundrum. Here, we used three lines of evidence - morphology, breeding trials and molecular approaches - to identify species boundaries in a group of ecologically important tabular Acropora corals. In contrast to previous studies, our morphological analyses yielded groups that were congruent with experimental crosses as well as with coalescent-based and allele sharing-based multilocus approaches to species delimitation. Our results suggest that species of the genus Acropora are reproductively isolated and independently evolving units that can be distinguished morphologically. These findings not only pave the way for a taxonomic revision of coral species, but also outline an approach that can provide a solid basis to address species delimitation and provide conservation support to a wide variety of keystone organisms.

Out of taxonomic crypsis: A new trans-arctic cryptic species pair corroborated by phylogenetics and molecular evidence

Cryptic species are a common phenomenon in cosmopolitan marine species. The use of molecular tools has often uncovered cryptic species occupying a fraction of the geographic range of the original morphospecies. Shipworms (Teredinidae) are marine bivalves, living in drift and fixed wood, many of which have a conserved morphology across cosmopolitan distributions. Herein novel and GenBank mitochondrial (cytochrome c oxidase subunit I) and nuclear (18S rRNA) DNA sequences are employed to produce a phylogeny of the Teredinidae and delimit a cryptic species pair in the Psiloteredo megotara complex. The anatomy, biogeography, and ecology of P. megotara, Psiloteredo sp. and Nototeredo edax are compared based on private and historic museum collections and a thorough literature review. Morphological and anatomical characters of P. megotara from the North Atlantic and Psiloteredo sp. from Japan were morphologically indistinguishable, and differ in pallet architecture and soft tissue anatomy from N. edax. The two Psiloteredo species were then delimited as genetically distinct species using four molecular-based methods. Consequently, the Northwest Pacific species, Psiloteredo pentagonalis, first synonymized with N. edax and then with P. megotara, is resurrected. Nototeredo edax, P. megotara and P. pentagonalis are redescribed based upon morphological and molecular characters. Phylogenetic analysis further revealed cryptic species complexes within the cosmopolitan species Bankia carinata and possibly additional cryptic lineages within the cosmopolitan Lyrodus pedicellatus.

A re-analysis of the data in Sharkey et al.'s (2021) minimalist revision reveals that BINs do not deserve names, but BOLD Systems needs a stronger commitment to open science

Halting biodiversity decline is one of the most critical challenges for humanity, but monitoring biodiversity is hampered by taxonomic impediments. One impediment is the large number of undescribed species (here called "dark taxon impediment") whereas another is caused by the large number of superficial species descriptions, that can only be resolved by consulting type specimens ("superficial description impediment"). Recently, Sharkey et al. (2021) proposed to address the dark taxon impediment for Costa Rican braconid wasps by describing 403 species based on COI barcode clusters ("BINs") computed by BOLD Systems. More than 99% of the BINs (387 of 390) were converted into species by assigning binominal names (e.g. BIN "BOLD:ACM9419" becomes Bracon federicomatarritai) and adding a minimal diagnosis (consisting only of a consensus barcode for most species). We here show that many of Sharkey et al.'s species are unstable when the underlying data are analyzed using different species delimitation algorithms. Add the insufficiently informative diagnoses, and many of these species will become the next "superficial description impediment" for braconid taxonomy because they will have to be tested and redescribed after obtaining sufficient evidence for confidently delimiting species. We furthermore show that Sharkey et al.'s approach of using consensus barcodes as diagnoses is not functional because it cannot be applied consistently. Lastly, we reiterate that COI alone is not suitable for delimiting and describing species, and voice concerns over Sharkey et al.'s uncritical use of BINs because they are calculated by a proprietary algorithm (RESL) that uses a mixture of public and private data. We urge authors, reviewers and editors to maintain high standards in taxonomy by only publishing new species that are rigorously delimited with open-access tools and supported by publicly available evidence.

Phylogeography and molecular species delimitation reveal cryptic diversity in Potamolithus (Caenogastropoda: Tateidae) of the southwest basin of the Andes

The species of the genus Potamolithus inhabiting the southwestern basin of the Andes are difficult to distinguish due to small size and similar shell morphology. Only Potamolithus australis and Potamolithus santiagensis have been traditionally recognized in this region, but the occurrence of several morphologically similar undescribed populations could increase the regional richness. Here we delimit described and potentially undescribed cryptic species of the genus using partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Network analysis and diversity indices inferred six highly differentiated haplogroups, many of them sympatric and widespread in the study area. Phylogeographic analyses suggest a scenario of recent diversification and the occurrence of multiple refuges during the successive Pleistocene glaciations. Phylogenetic analysis also recovered six major clades that showed no relationship with physiography. Species delimitation analyses consistently recognized three or four candidate species apart from P. australis and P. santiagensis. Divergence times indicate that speciation of Chilean Potamolithus began at the end of the Pliocene, probably driven by climatic rather than geographic events. Considering the high inter- and intra-basin genetic diversity, conservation efforts should be focused on protecting sympatric taxa in the basins with the highest species richness.

SPART: A versatile and standardized data exchange format for species partition information

A wide range of data types can be used to delimit species and various computer-based tools dedicated to this task are now available. Although these formalized approaches have significantly contributed to increase the objectivity of species delimitation (SD) under different assumptions, they are not routinely used by alpha-taxonomists. One obvious shortcoming is the lack of interoperability among the various independently developed SD programs. Given the frequent incongruences between species partitions inferred by different SD approaches, researchers applying these methods often seek to compare these alternative species partitions to evaluate the robustness of the species boundaries. This procedure is excessively time consuming at present, and the lack of a standard format for species partitions is a major obstacle. Here, we propose a standardized format, SPART, to enable compatibility between different SD tools exporting or importing partitions. This format reports the partitions and describes, for each of them, the assignment of individuals to the "inferred species". The syntax also allows support values to be optionally reported, as well as original trees and the full command lines used in the respective SD analyses. Two variants of this format are proposed, overall using the same terminology but presenting the data either optimized for human readability (matricial SPART) or in a format in which each partition forms a separate block (SPART.XML). ABGD, DELINEATE, GMYC, PTP and TR2 have already been adapted to output SPART files and a new version of LIMES has been developed to import, export, merge and split them.

Diversity, biogeography, evolutionary relationships, and conservation of Eastern Mediterranean freshwater mussels (Bivalvia: Unionidae)

Located at the junction between Europe, Africa, and Asia, with distinct evolutionary origins and varied ecological and geographical settings, together with a marked history of changes in orogeny and configuration of the main river basins, turned the Eastern Mediterranean into a region of high diversity and endemism of freshwater taxa. Freshwater mussels (Bivalvia, Unionidae) from the Western Palearctic have been widely studied in their European range, but little attention has been dedicated to these taxa in the Eastern Mediterranean region and their diversity and phylogeography are still poorly understood. The present study aims to resolve the diversity, biogeography, and evolutionary relationships of the Eastern Mediterranean freshwater mussels. To that end, we performed multiple field surveys, phylogenetic analyses, and a thorough taxonomic revaluation. We reassessed the systematics of all Unionidae species in the region, including newly collected specimens across Turkey, Israel, and Iran, combining COI + 16S + 28S and COI phylogenies with molecular species delineation methods. Phylogeographical patterns were characterized based on published molecular data, newly sequenced specimens, and species distribution data, as well as ancestral range estimations. We reveal that Unionidae species richness in the Eastern Mediterranean is over 70% higher than previously assumed, counting 19 species within two subfamilies, the Unioninae (14) and Gonideinae (5). We propose two new species, Anodonta seddonisp. nov. and Leguminaia anatolicasp. nov. Six additional taxa, Unio delicatusstat. rev., Unio eucirrusstat. rev., Unio huetistat. rev., Unio sesirmensisstat. rev., Unio terminalisstat. rev. removed from the synonymy of Unio tigridis, as well as Unio damascensisstat. rev. removed from the synonymy of Unio crassus, are re-described. The nominal taxa Unio rothi var. komarowi O. Boettger, 1880 and Unio armeniacus Kobelt, 1911 are proposed as new synonyms of Unio bruguierianus, and Anodonta cyrea Drouët, 1881 and Anodonta cilicica Kobelt & Rolle, 1895 as new synonyms of Anodonta anatina. Also, the presence of Unio tumidus in the Maritza River is confirmed. The phylogeographic patterns described here are interpreted concerning major past geological events. Conservation needs and implications are presented, together with populations and species conservation priorities.

Divergence time and species delimitation of microbivalves in the Southern Ocean: the case of Kidderia species

The systematics of Subantarctic and Antarctic near-shore marine benthic invertebrates requires major revision and highlights the necessity to incorporate additional sources of information in the specimen identification chart in the Southern Ocean (SO). In this study, we aim to improve our understanding of the biodiversity of Kidderia (Dall 1876) through molecular and morphological comparisons of Antarctic and Subantarctic taxa. The microbivalves of the genus Kidderia are small brooding organisms that inhabit intertidal and shallow subtidal rocky ecosystems. This genus represents an interesting model to test the vicariance and dispersal hypothesis in the biogeography of the SO. However, the description of Kidderia species relies on a few morphological characters and biogeographic records that raise questions about the true diversity in the group. Here we will define the specimens collected with genetic tools, delimiting their respective boundaries across provinces of the SO, validating the presence of two species of Kidderia. Through the revision of taxonomic issues and species delimitation, it was possible to report that the Antarctic species is Kidderia subquadrata and the species recorded in the Subantarctic islands Diego Ramirez, South Georgia and the Kerguelen Archipelago is Kidderia minuta. The divergence time estimation suggests the origin and diversification of Kidderia lineages are related to historical vicariant processes probably associated with the separation of the continental landmasses close to the late Eocene.

Complete chloroplast genomes of Leptodermis scabrida complex: Comparative genomic analyses and phylogenetic relationships

Leptodermis scabrida complex is one of the important components of genus Leptodermis, which is mainly distributed in the Himalaya Mountains. It includes species of L. gracilis, L. hirsutiflora, L. hirsutiflora var. ciliata, L. kumaonensis, L. pilosa var. acanthoclada and L. scabrida. However, species boundaries and relationships within this complex are unclear based on current morphological and molecular evidence. We sequenced 13 complete chloroplast (cp) genomes representing seven taxa of the complex and two non-Leptodermis scabrida complex taxa. After de novo assembly and annotation, we performed comparative genomic analysis. All cp genomes showed highly conserved structures, and the genome sizes ranged from 154,369 bp to 154,885 bp and possessed the same GC content (37.5%). A total of 113 unique genes were identified in each cp sample, including 79 protein coding genes, 30 tRNAs, and four rRNAs. Repeat sequences and SSRs were detected, showing great similarity among all taxa in this complex. Six highly variable regions, including trnS-trnG, rps2-rpoC2, ndhF, rpl32-ccsA, ccsA-ndhD, and ndhA, were screened as potential molecular markers for phylogenetic reconstruction. Based on a total of 27 complete cp genome sequences, the consistent and robust phylogenetic relationships were firstly constructed and the same species within L. scabrida complex clustered into a group. The divergence time of Leptodermis from ancestral taxa occurred at the middle Eocene, which might be due to geological and climatic changes. The 13 complete cp genome sequences reported will provide new clues for phylogeny elucidation, species identification and evolutionary history speculation of Leptodermis, as well as in Rubiaceae.

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.

Molecular taxonomy confirms that the northeastern Atlantic and Mediterranean Sea harbor a single lancelet, Branchiostoma lanceolatum (Pallas, 1774) (Cephalochordata: Leptocardii: Branchiostomatidae)

Branchiostomatidae (lancelets or amphioxus) comprises about 30 species, several of which are well-established models in evolutionary development. Our zoological and ecological knowledge of the family is nonetheless limited. Despite evident differences can be found among known populations, the taxonomy of Branchiostoma lanceolatum (type species of the genus Branchiostoma) has never been investigated with modern methods through its range in the northeastern Atlantic and Mediterranean Sea. We address this via a multilocus molecular approach and comparing specimens collected from different European populations. Results obtained here confirm the presence of a single species inhabiting the range between the topotypical localities of B. lanceolatum (Atlantic Ocean) and of its junior synonym B. lubricum (Mediterranean Sea), without evincing geographical structure between populations. This suggests that environment most likely drives the characteristics observed in different geographic areas. The long larval phase and the slow mutation rate in lancelets may have played a key role in the evolutionary history of this iconic species.

A multilocus assessment reveals two new synonymies for East Asian Cyclommatus stag beetles (Coleoptera, Lucanidae)

Cyclommatus scutellaris Möllenkamp, 1912, Cyclommatus elsae Kriesche, 1921 and Cyclommatus tamdaoensis Fujita, 2010 are East Asian stag beetle species with long-debated taxonomic relationships due to high intraspecific morphological variability. In this study, we applied multilocus phylogenetic analyses to reassess their relationships. Two mitochondrial genes (16S rDNA, COI) and two nuclear genes (28S rDNA, Wingless) were used to reconstruct the phylogeny through the Bayesian inference (BI) and Maximum Likelihood (ML) methods. Both topologies supported two clades: the clade C. scutellaris was sister to the clade (C. elsae + C. tamdaoensis) with the subclade C. tamdaoensis embedded in the subclade C. elsae. The Kimura 2-parameter (K2P) genetic distance analysis yielded a low mean value (≤0.035) among the three taxa, which was well below the minimum mean value between other Cyclommatus species (≥0.122). We also compared the accuracy and efficiency of two approaches, GMYC and ABGD, in delimitating the three lineages. The result shows that ABGD is a better approach than GMYC. Our molecular data recognizes the three species as different populations of a single species, ranging from Taiwan Island to the continent. Therefore, we propose two new junior synonyms for C. scutellaris: C. tamdaoensis, syn. nov. and C. elsae syn. nov.

Unveiling cryptic diversity among Müllerian co-mimics: insights from the Western Palaearctic Syntomis moths (Lepidoptera: Erebidae: Arctiinae)

Accurate species delimitation is of primary importance in biodiversity assessments and in reconstructing patterns and processes in the diversification of life. However, the discovery of cryptic species in virtually all taxonomic groups unveiled significant gaps in our knowledge of biodiversity. Mimicry complexes are good candidates to source for cryptic species. Indeed, members of mimicry complexes undergo selective pressures on their habitus, which results in strong resemblance even between distantly related species. In this study, we used a multi-locus genetic approach to investigate the presence of cryptic diversity within a group of mimetic day-flying moths whose systematics has long been controversial, the Euro-Anatolian Syntomis. Results showed incongruence between species boundaries and the currently accepted taxonomy of this group. Both mitochondrial and nuclear markers indicate the presence of four, well-distinct genetic lineages. The genetic distance and time of divergence between the Balkan and Italian populations of S. marjana are the same as those found between S. phegea and S. ragazzii, the last two being well-distinct, broadly sympatrically occurring species. The divergence between the two lineages of S. marjana dates back to the Early Pleistocene, which coincided with substantial changes in climatic conditions and vegetation cover in Southern Europe that have likely induced geographic and ecological vicariance. Syntomis populations belonging to the taxa kruegeri (s. str.), albionica and quercii are now considered a separate species from marjana s. str. and are thus distinguished as Syntomis quercii Verity, 1914, bona sp., stat. nov. Our results show that the species richness of mimicry complexes inhabiting temperate regions might still be severely underestimated.

Integrative Taxonomy and Molecular Phylogeny of the Plant-Parasitic Nematode Genus Paratylenchus (Nematoda: Paratylenchinae): Linking Species with Molecular Barcodes

Pin nematodes of the genus Paratylenchus are obligate ectoparasites of a wide variety of plants that are distributed worldwide. In this study, individual morphologically vouchered nematode specimens of fourteen Paratylenchus species, including P. aculentus, P. elachistus, P. goodeyi, P. holdemani, P. idalimus, P. microdorus, P. nanus, P. neoamblycephalus, P. straeleni and P. veruculatus, are unequivocally linked to the D2-D3 of 28S, ITS, 18S rRNA and COI gene sequences. Combined with scanning electron microscopy and a molecular analysis of an additional nine known and thirteen unknown species originating from diverse geographic regions, a total of 92 D2-D3 of 28S, 41 ITS, 57 18S rRNA and 111 COI new gene sequences are presented. Paratylenchus elachistus, P. holdemani and P. neoamblycephalus are recorded for the first time in Belgium and P. idalimus for the first time in Europe. Paratylenchus is an excellent example of an incredibly diverse yet morphologically minimalistic plant-parasitic genus, and this study provides an integrated analysis of all available data, including coalescence-based molecular species delimitation, resulting in an updated Paratylenchus phylogeny and the corrective reassignment of 18 D2-D3 of 28S, 3 ITS, 3 18S rRNA and 25 COI gene sequences that were previously unidentified or incorrectly classified.

An ocean yet to be discovered: increasing systematic knowledge of Indo-Pacific Okenia Menke, 1830 (Nudibranchia:Goniodorididae)

Numerous faunistic and ecological studies have been conducted throughout the Indo-Pacific Ocean to assess its biodiversity. Despite the abundance of research, studies on the species that inhabit the Indo-Pacific are still necessary due to its extent and high species richness. The major species richness of the genus Okenia Menke, 1830 (Nudibranchia, Goniodorididae) is found in the Indo-Pacific Ocean, including 38 of 60 valid species. Nevertheless, this number does not represent the real biodiversity, since at least 20 more species are already reported in field-guides as undescribed species belonging to this genus. The systematics of the genus Okenia are still unclear since it has been the subject of only a few and incomplete studies. In the present paper, we describe five new Okenia species from the coastlines of Japan, Mozambique and Australia: Okenia aurorapapillata sp. nov., Okenia elisae sp. nov., Okenia nakanoae sp. nov., Okenia siderata sp. nov. and Okenia tenuifibrata sp. nov. Moreover, anatomical details not previously described of Okenia atkinsonorum, Okenia barnardi, Okenia cf. echinata, Okenia hallucigenia, Okenia hiroi, Okenia japonica, Okenia pellucida, Okenia pilosa and Okenia rhinorma are provided. New partial sequences of standard markers (COI, 16S rRNA and H3) were obtained and a phylogenetic analysis that included all species with available data was performed. ZOOBANK urn:lsid:zoobank.org:pub:28AE2536-A264–4194–8AE3-C430620572E7

Few and far apart: integrative taxonomy of Australian species of Gladiobela and Pagodibela (Conoidea : Raphitomidae) reveals patterns of wide distributions and low abundance

The deep-sea malacofauna of temperate Australia remains comparatively poorly known. However, a recent influx of DNA-suitable material obtained from a series of deep-sea cruises has facilitated integrative taxonomic study on the Conoidea (Caenogastropoda:Neogastropoda). Building on a recent molecular phylogeny of the conoidean family Raphitomidae, this study focussed on the genera Gladiobela and Pagodibela (both Criscione, Hallan, Puillandre &amp; Fedosov, 2020). We subjected a representative mtDNA cox1 dataset of deep-sea raphitomids to ABGD, which recognised 14 primary species hypotheses (PSHs), 9 of which were converted to secondary species hypotheses (SSHs). Following the additional examination of the shell and hypodermic radula features, as well as consideration of bathymetric and geographic data, seven of these SSHs were recognised as new to science and given full species rank. Subsequently, systematic descriptions are provided herein. Of these, five are attributed to Gladiobela (three of which are endemic to Australia and two more widely distributed) and two are placed in Pagodibela (one endemic to southern Australia and one widespread in the Pacific). The rarity of many ‘turrids’ reported in previous studies is confirmed herein, as particularly indicated by highly disjunct geographic records for two taxa. Additionally, several of the studied taxa exhibit wide Indo-Pacific distributions, suggesting that wide geographic ranges in deep-sea ‘turrids’ may be more common than previously assumed. Finally, impediments to deep-sea ‘turrid’ taxonomy in light of such comparative rarity and unexpectedly wide distributions are discussed.

Deciphering the Taxonomic Delimitation of Ottelia acuminata (Hydrocharitaceae) Using Complete Plastomes as Super-Barcodes

Accurate species delimitation and identification, which is a challenging task in traditional morphology-based taxonomy, is crucial to species conservation. Ottelia acuminata (Hydrocharitaceae) is a severely threatened submerged macrophyte endemic to southwestern China. The taxonomy of O. acuminata, which has long been in dispute, remains unresolved, impeding effective conservation and management practices. Here, we aim to address the long-standing issues concerning species boundary and intraspecific subdivision of O. acuminata using complete plastome sequences as super-barcodes. The taxonomic delimitation of O. acuminata was explored using phylogenetic inference and two independent sequence-based species delimitation schemes: automatic barcode gap discovery (ABGD) and multi-rate Poisson tree processes (mPTP). The reciprocally reinforcing results support the reduction of the closely related congeneric species, O. balansae and O. guanyangensis, as two conspecific varieties of O. acuminata. Within the newly defined O. acuminata, accurate varietal identification can be achieved using plastome super-barcodes. These findings will help inform future decisions regarding conservation, management and restoration of O. acuminata. This case study suggests that the use of plastome super-barcodes can provide a solution for species delimitation and identification in taxonomically difficult plant taxa, thus providing great potential to lessen the challenges of inventorying biodiversity, as well as biologically monitoring and assessing threatened species.

Snails in depth: integrative taxonomy of Famelica, Glaciotomella and Rimosodaphnella (Conoidea: Raphitomidae) from the deep sea of temperate Australia

The deep sea of temperate south-eastern Australia appears to be a ‘hotspot’ for diversity and endemism of conoidean neogastropods of the family Raphitomidae. Following a series of expeditions in the region, a considerable amount of relevant DNA-suitable material has become available. A molecular phylogeny based on this material has facilitated the identification of diagnostic morphological characters, allowing the circumscription of monophyletic genera and the introduction of several new genus-level taxa. Both named and new genera are presently being investigated through integrative taxonomy, with the discovery of a significant number of undescribed species. As part of this ongoing investigation, our study focuses on the genera Famelica Bouchet &amp; Warén, 1980, Glaciotomella Criscione, Hallan, Fedosov &amp; Puillandre, 2020 and Rimosodaphnella Cossmann, 1914. We subjected a comprehensive mitochondrial DNA dataset of representative deep-sea raphitomids to the species delimitation methods ABGD and ASAP that recognised 18 and 15 primary species hypotheses (PSHs) respectively. Following additional evaluation of shell and radular features, and examination of geographic and bathymetric ranges, nine of these PSHs were converted to secondary species hypotheses (SSHs). Four SSHs (two in Famelica and two in Rimosodaphnella) were recognised as new, and formal descriptions are provided herein.