Delimiting Species-Poor Data Sets using Single Molecular Markers: A Study of Barcode Gaps, Haplowebs and GMYC

Genetic and morphological evidence support the specific status of the endemic Ericaandevalensis (Ericales, Ericaceae)

Assessing the taxonomic status of closely related taxa is crucial in plant systematics and can have important implications for conservation and human plant use. Ericaandevalensis Cabezudo & Rivera is a metallophyte endemic species from highly metal-polluted soils of SW Iberian Peninsula, an area with a mining history going back more than 5,000 years. Ericaandevalensis is closely related to Ericamackayana Bab., a northern Iberian species also present in western Ireland. The status of E.andevalensis as a species or subspecies subordinated to E.mackayana is subject to debate. Here, we assessed the genetic and phenotypic relationship between both species, including the population structure of E.andevalensis. We used high throughput sequencing to determine genome-wide Single Nucleotide Polymorphisms (SNPs), and morphometric analyses from 35 reproductive and vegetative traits. The morphological analysis showed at least eight characters that can discriminate the two species, from which ovary hairiness and the size of leaf glandular hairs were the most informative. Genetic analyses showed that each species formed a monophyletic cluster with full support, separated by an interspecific genetic distance >4-fold higher than intra-specific distance. Population genetic analyses of E.andevalensis shows that populations are highly structured, with the Portuguese one as the most isolated and less variable. These results support the recognition of E.andevalensis as a distinct species with a highly constrained ecological requirements and a narrow geographic distribution, but with a limited gene flow between populations. We discuss the implications of these outcomes in conservation policies and potential uses of E.andevalensis such as decontamination of polluted soils.

Exploring Barbronia species diversity and phylogenetic relationship within Suborder Erpobdelliformes (Clitellata: Annelida)

Background

Barbronia, a genus of freshwater macrophagous leeches, belongs to Erpobdelliformes (Salifidae: Clitellata: Annelida), and B. weberi, a well-known leech within this genus, has a worldwide distribution. However, the systematics of Barbronia have not yet been adequately investigated, primarily due to a few molecular markers, and only 20 Barbronia sequences available in the GenBank database. This gap significantly limits our understanding of the Barbronia species identification, as well as the phylogenetic placement of the genus Barbronia within Salifidae.

Methods

Next-generation sequencing (NGS) was used to simultaneously capture the entire mitochondrial genome and the full-length 18S/28S rDNA sequences. The species boundary of Barbronia species was estimated using bGMYC and bPTP methods, based on all available Barbronia COI sequences. Uncorrected COI p-distance was calculated in MEGA. A molecular data matrix consisting of four loci (COI, 12S, 18S, and 28S rDNA) for outgroups (three Haemopis leeches) and 49 erpobdellid leeches, representing eight genera within the Suborder Erpobdelliformes was aligned using MAFFT and LocARNA. This matrix was used to reconstruct the phylogenetic relationship of Barbronia via Bayesian inference (BI) and the maximum likelihood (ML) method.

Results

The full lengths of the mitochondrial genome, 18S and 28S rDNAs of B. cf. gwalagwalensis, are 14847 bp, 1876 bp 1876 bp, and 2863 bp, respectively. Both bGMYC and bPTP results based on COI data are generally congruent, suggesting that the previously proposed taxa (B. arcana, B. weberi formosana, and B. wuttkei or Erpobdella wuttkei) are synonyms of B. weberi. The specimens listed in the B. gwalagwalensis group, however, are split into at least two Primary Species Hypotheses (PSHs). The p-distance of the first PSH is less than 1.3% but increased to 4.5% when including the secondary PSH (i.e., B. cf. gwalagwalensis). In comparison, the interspecific p-distance between the B. weberi group and the B. gwalagwalensis group ranged from 6.4% to 8.7%, and the intraspecific p-distance within the B. weberi group is less than 0.8%. Considering the species delimitation results and the sufficient large p-distance, the specimen sampled in China is treated as B. cf. gwalagwalensis. The monophyly of the four Erpobdelliformes families Salifidae, Orobdellidae, Gastrostomobdellidae sensu stricto and Erpobdellidae is well supported in ML and BI analysis based on a data of four markers. Within the Salifidae, a well-supported Barbronia is closely related to a clade containing Odontobdella and Mimobdella, and these three genera are sister to a clade consisted of Salifa and Linta. According to the results of this study, the strategy of simultaneous obtaining both whole mitochondria and nuclear markers from extensively sampled Salifids species using NGS is expected to fathom both the species diversity of B. gwalagwalensis and the evolutionary relationship of Salifidae.

Genetic Diversity and DNA Barcoding of Thrips in Bangladesh

Thrips are economically important pests, and some species transmit plant viruses that are widely distributed and can damage vegetables and cash crops. Although few studies on thrips species have been conducted in Bangladesh, the variation and genetic diversity of thrips species remain unknown. In this study, we collected thrips samples from 16 geographical locations throughout the country and determined the nucleotide sequences of the mitochondrial cytochrome c oxidase subunit 1 (mtCOI) gene in 207 thrips individuals. Phylogenetic analysis revealed ten genera (Thrips, Haplothrips, Megalothrips, Scirtothrips, Frankliniella, Dendrothripoides, Astrothrips, Microcephalothrips, Ayyaria, and Bathrips) and 19 species of thrips to inhabit Bangladesh. Among these, ten species had not been previously reported in Bangladesh. Intraspecific genetic variation was diverse for each species. Notably, Thrips palmi was the most genetically diverse species, containing 14 haplotypes. The Mantel test revealed no correlation between genetic and geographical distances. This study revealed that thrips species are expanding their host ranges and geographical distributions, which provides valuable insights into monitoring the diversity of and control strategies for these pests.

DNA barcodes provide insights into the diversity and biogeography of the non-biting midge Polypedilum (Diptera, Chironomidae) in South America

South America, particularly within its tropical belt, is renowned for its unparalleled high levels of species richness, surpassing other major biomes. Certain neotropical areas harbor fragmented knowledge of insect diversity and face imminent threats from biodiversity loss and climate change. Hence, there is an urgent need for rapid estimation methods to complement slower traditional taxonomic approaches. A variety of algorithms for delimiting species through single-locus DNA barcodes have been developed and applied for rapid species diversity estimates across diverse taxa. However, tree-based and distance-based methods may yield different group assignments, leading to potential overestimation or underestimation of putative species. Here, we investigate the performance of different DNA-based species delimitation approaches to rapidly estimate the diversity of Polypedilum (Chironomidae, Diptera) in South America. Additionally, we test the hypothesis that significant differences exist in the community structure of Polypedilum fauna between South America and its neighboring regions, particularly the Nearctic. Our analysis encompasses a dataset of 1492 specimens from 598 locations worldwide, with a specific focus on South America. Within this region, we analyzed a subset of 247 specimens reported from 37 locations. Using various methods including the Barcode Index Number (BIN), Bayesian Poisson tree processes (bPTP), multi-rate Poisson tree processes (mPTP), single-rate Poisson tree processes (sPTP), and generalized mixed Yule coalescent (sGMYC), we identify molecular operational taxonomic units (MOTUs) ranging from 267 to 520. Our results indicate that the sGMYC method is the most suitable for estimating putative species in our dataset, resulting in the identification of 75 species in the Neotropical region, particularly in South America. Notably, this region exhibited higher species richness in comparison to the Palearctic and Oriental realms. Additionally, our findings suggest potential differences in species composition of Polypedilum fauna between the Neotropical and the adjacent Nearctic realms, highlighting high levels of endemism and species richness in the first. These results support our hypothesis that there are substantial differences exist in species composition between the Polypedilum fauna in South America and the neighboring regions.

Speciation on the Roof of the World: Parallel Fast Evolution of Cryptic Mole Vole Species in the Pamir-Alay-Tien Shan Region

Speciation is not always accompanied by morphological changes; numerous cryptic closely related species were revealed using genetic methods. In natural populations of Ellobius tancrei (2n = 54-30) and E. alaicus (2n = 52-48) of the Pamir-Alay and Tien Shan, the chromosomal variability due to Robertsonian translocations has been revealed. Here, by comprehensive genetic analysis (karyological analyses as well as sequencing of mitochondrial genes, cytb and COI, and nuclear genes, XIST and IRBP) of E. alaicus and E. tancrei samples from the Inner Tien Shan, the Alay Valley, and the Pamir-Alay, we demonstrated fast and independent diversification of these species. We described an incompletely consistent polymorphism of the mitochondrial and nuclear markers, which arose presumably because of habitat fragmentation in the highlands, rapid karyotype changes, and hybridization of different intraspecific varieties and species. The most intriguing results are a low level of genetic distances calculated from mitochondrial and nuclear genes between some phylogenetic lines of E. tancrei and E. alaicus, as well significant species-specific chromosome variability in both species. The chromosomal rearrangements are what most clearly define species specificity and provide further diversification. The "mosaicism" and inconsistency in polymorphism patterns are evidence of rapid speciation in these mammals.

Surrounded by Kindred: Spermophilus major Hybridization with Other Spermophilus Species in Space and Time

Among the numerous described cases of hybridization in mammals, the most intriguing are (a) cases of introgressive hybridization deeply affecting the evolutionary history of species, and (b) models involving not a pair of species but a multi-species complex. Therefore, the hybridization history of the russet ground squirrel Spermophilus major, whose range has repeatedly changed due to climatic fluctuations and now borders the ranges of four related species, is of great interest. The main aims of this study were to determine the direction and intensity of gene introgression, the spatial depth of the infiltration of extraneous genes into the S. major range, and to refine the hypothesis of the hybridogenic replacement of mitochondrial genomes in the studied group. Using phylogenetic analysis of the variability of mitochondrial (CR, cytb) and nuclear (SmcY, BGN, PRKCI, c-myc, i6p53) markers, we determined the contribution of neighboring species to the S. major genome. We showed that 36% of S. major individuals had extraneous alleles. All peripheral species that were in contact with S. major contributed towards its genetic variability. We also proposed a hypothesis for the sequence and localization of serial hybridization events. Our assessment of the S. major genome implications of introgression highlights the importance of implementing conservation measures to protect this species.

Comparing the Efficiency of Single-Locus Species Delimitation Methods within Trochoidea (Gastropoda: Vetigastropoda)

In the context of diminishing global biodiversity, the validity and practicality of species delimitation methods for the identification of many neglected and undescribed biodiverse species have been paid increasing attention. DNA sequence-based species delimitation methods are mainly classified into two categories, namely, distance-based and tree-based methods, and have been widely adopted in many studies. In the present study, we performed three distance-based (ad hoc threshold, ABGD, and ASAP) and four tree-based (sGMYC, mGMYC, PTP, and mPTP) analyses based on Trochoidea COI data and analyzed the discordance between them. Moreover, we also observed the performance of these methods at different taxonomic ranks (the genus, subfamily, and family ranks). The results suggested that the distance-based approach is generally superior to the tree-based approach, with the ASAP method being the most efficient. In terms of phylogenetic methods, the single threshold version performed better than the multiple threshold version of GMYC, and PTP showed higher efficiency than mPTP in delimiting species. Additionally, GMYC was found to be significantly influenced by taxonomic rank, showing poorer efficiency in datasets at the genus level than at higher levels. Finally, our results highlighted that cryptic diversity within Trochoidea (Mollusca: Vetigastropoda) might be underestimated, which provides quantitative evidence for excavating the cryptic lineages of these species.

A DNA barcode reference library for endemic Ponto-Caspian amphipods

The Ponto-Caspian region is an endemicity hotspot that harbours several crustacean radiations, among which amphipods are the most diverse. These poorly known species are severely threatened in their native range, while at the same time they are invading European inland waters with significant ecological consequences. A proper taxonomic knowledge of this fauna is paramount for its conservation within the native region and monitoring outside of it. Here, we assemble a DNA barcode reference library for nearly 60% of all known Ponto-Caspian amphipod species. We use several methods to define molecular operational taxonomic units (MOTUs), based on two mitochondrial markers (COI and 16S), and assess their congruence with current species-level taxonomy based on morphology. Depending on the method, we find that 54–69% of species had congruent morpho-molecular boundaries. The cases of incongruence resulted from lumping distinct morphospecies into a single MOTU (7–27%), splitting a morphospecies into several MOTUs (4–28%), or both (4–11%). MOTUs defined by distance-based methods without a priori divergence thresholds showed the highest congruence with morphological taxonomy. These results indicate that DNA barcoding is valuable for clarifying the diversity of Ponto-Caspian amphipods, but reveals that extensive work is needed to resolve taxonomic uncertainties. Our study advances the DNA barcode reference library for the European aquatic biota, paving the way towards improved taxonomic knowledge needed to enhance monitoring and conservation efforts.

Lack of Statistical Rigor in DNA Barcoding Likely Invalidates the Presence of a True Species' Barcode Gap

DNA barcoding has been largely successful in satisfactorily exposing levels of standing genetic diversity for a wide range of taxonomic groups through the employment of only one or a few universal gene markers. However, sufficient coverage of geographically-broad intra-specific haplotype variation within genomic databases like the Barcode of Life Data Systems (BOLD) and GenBank remains relatively sparse. As reference sequence libraries continue to grow exponentially in size, there is now the need to identify novel ways of meaningfully analyzing vast amounts of available DNA barcode data. This is an important issue to address promptly for the routine tasks of specimen identification and species discovery, which have seen broad adoption in areas as diverse as regulatory forensics and resource conservation. Here, it is demonstrated that the interpretation of DNA barcoding data is lacking in statistical rigor. To highlight this, focus is set specifically on one key concept that has become a household name in the field: the DNA barcode gap. Arguments outlined herein specifically center on DNA barcoding in animal taxa and stem from three angles: (1) the improper allocation of specimen sampling effort necessary to capture adequate levels of within-species genetic variation, (2) failing to properly visualize intra-specific and interspecific genetic distances, and (3) the inconsistent, inappropriate use, or absence of statistical inferential procedures in DNA barcoding gap analyses. Furthermore, simple statistical solutions are outlined which can greatly propel the use of DNA barcoding as a tool to irrefutably match unknowns to knowns on the basis of the barcoding gap with a high degree of confidence. Proposed methods examined herein are illustrated through application to DNA barcode sequence data from Canadian Pacific fish species as a case study.

Uncertainties and risks in delimiting species of Cambeva (Siluriformes: Trichomycteridae) with single-locus methods and geographically restricted data

Abstract Cambeva contains species with complex taxonomy or poorly delimitated in terms of morphology and geopraphic distribution. We conducted an extensive review of Cambeva populations from coastal drainages of Southern to Southeastern Brazil to evaluate species geographic limits with an integrative analysis including morphological and molecular data (COI). We test if two single-locus methods, Bayesian Poisson Tree Processes (bPTP) and Generalized Mixed Yule Coalescent (GMYC), are efficient to delimit species boundaries in Cambeva by the comparison with the diagnosable morphological units. Using GMYC, we also evaluated the combination of tree and molecular clock priors to reconstruct the input phylogeny and assessed how well the implemented model fitted our empirical data. Eleven species were identified using a morphological diagnosability criterion: Cambeva balios, C. barbosae, C. botuvera, C. cubataonis, C. davisi, C. guaraquessaba, C. iheringi, C. tupinamba, and C. zonata and two treated as undescribed species. In contrast with previous knowledge, many of them have wider distribution and high intraspecific variation. Species delimitation based on single-locus demonstrated incongruences between the methods and strongly differed from the morphological delimitation. These disagreements and the violation of the GMYC model suggest that a single-locus data is insufficient to delimit Cambeva species and the failure may be attributable to events of mitochondrial introgression and incomplete lineage sorting.

Hidden diversity revealed in the freshwater snails, Bythinella and Pseudamnicola, in the Island of Crete

Crete with its complex geomorphological history is the island with the highest number of endemism observed in animal and plant taxa throughout the Aegean archipelago. While other groups of organisms within Crete are well-studied, the freshwater gastropod fauna still remains poorly investigated. Bythinella and Pseudamnicola, are 2 genera of freshwater springsnails, both present on the island, inhabiting springs and other freshwater habitats. Here, we conduct a comprehensive study on the distribution of the different genetic lineages of the 2 gastropod genera in order to assess the mode of their differentiation on the island and infer the actual number of species present in the island. Towards these aims, sequence data from the mitochondrial gene were used and analyzed within a phylogenetic framework. For Bythinella, our results strongly support at least 5 delineated Bythinella spp. inhabiting Crete, which correspond to the already described species from previous studies with the addition of a new one. Bythinella analyses reveal an old time-frame of differentiation with vicariant phenomena being more likely the main drivers shaping the present-day distribution of the genus' genetic lineages. For Pseudamnicola, our data indicate the presence of at least 2 delineated Pseudamnicola spp. with a differentiation more consistent to an isolation-by-distance pattern of a relatively recent origin. Dispersion processes followed by isolation of the populations and/or recent speciation, seem to be the underlying process for the current distribution of Pseudamnicola lineages.

Type specimen sequencing, multilocus analyses, and species delimitation methods recognize the cosmopolitan Corallina berteroi and establish the northern Japanese C. yendoi sp. nov. (Corallinaceae, Rhodophyta)

A partial rbcL sequence of the lectotype specimen of Corallina berteroi shows that it is the earliest available name for C. ferreyrae. Multilocus species delimitation analyses (ABGD, SPN, GMYC, bPTP, and BPP) using independent or concatenated COI, psbA, and rbcL sequences recognized one, two, or three species in this complex, but only with weak support for each species hypothesis. Conservatively, we recognize a single worldwide species in this complex of what appears to be multiple, evolving populations. Included in this species, besides C. ferreyrae, are C. caespitosa, the morphologically distinct C. melobesioides, and, based on a partial rbcL sequence of the holotype specimen, C. pinnatifolia. Corallina berteroi, not C. officinalis, is the cosmopolitan temperate species found thus far in the NE Atlantic, Mediterranean Sea, warm temperate NW Atlantic and NE Pacific, cold temperate SW Atlantic (Falkland Islands), cold and warm temperate SE Pacific, NW Pacific and southern Australia. Also proposed is C. yendoi sp. nov. from Hokkaido, Japan, which was recognized as distinct by 10 of the 13 species discrimination analyses, including the multilocus BPP.

Blue mussels of the Mytilus edulis species complex from South America: The application of species delimitation models to DNA sequence variation

Smooth-shelled blue mussels, Mytilus spp., have a worldwide antitropical distribution and are ecologically and economically important. Mussels of the Mytilus edulis species complex have been the focus of numerous taxonomic and biogeographical studies, in particular in the Northern hemisphere, but the taxonomic classification of mussels from South America remains unclear. The present study analysed 348 mussels from 20 sites in Argentina, Chile, Uruguay and the Falkland Islands on the Atlantic and Pacific coasts of South America. We sequenced two mitochondrial locus, Cytochrome c Oxidase subunit I (625 bp) and 16S rDNA (443 bp), and one nuclear gene, ribosomal 18S rDNA (1770 bp). Mitochondrial and nuclear loci were analysed separately and in combination using maximum likelihood and Bayesian inference methods to identify the combination of the most informative dataset and model. Species delimitation using five different models (GMYC single, bGMYC, PTP, bPTP and BPP) revealed that the Mytilus edulis complex in South America is represented by three species: native M. chilensis, M. edulis, and introduced Northern Hemisphere M. galloprovincialis. However, all models failed to delimit the putative species Mytilus platensis. In contrast, however, broad spatial scale genetic structure in South America using Geneland software to analyse COI sequence variation revealed a group of native mussels (putatively M. platensis) in central Argentina and the Falkland Islands. We discuss the scope of species delimitation methods and the use of nuclear and mitochondrial genetic data to the recognition of species within the Mytilus edulis complex at regional and global scales.

The effect of landscape and human settlement on the genetic differentiation and presence of Paragonimus species in Mesoamerica

Foodborne diseases are a neglected research area, and despite the existence of many tools for diagnosis and genetic studies, very little is known about the effect of the landscape on the genetic diversity and presence of parasites. One of these foodborne disease is paragonimiasis, caused by trematodes of the genus Paragonimus, which is responsible for a high number of infections in humans and wild animals. The main Paragonimus sp reported in Mesoamerica is Paragonimus mexicanus, yet there are doubts about its correct identification as a unique species throughout the region. This, together with a lack of detailed knowledge about their ecology, evolution and differentiation, may complicate the implementation of control strategies across the Mesoamerican region. We had the goal of delimiting the species of P. mexicanus found throughout Mesoamerica and determining the effect of landscape and geology on the diversity and presence of the parasite. We found support for the delimitation of five genetic groups. The genetic differentiation among these groups was positively affected by elevation and the isolation of river basins, while the parasite's presence was affected negatively only by the presence of human settlements. These results suggest that areas with lower elevation, connected rivers basins, and an absence of human settlements have low genetic differentiation and high P. mexicanus presence, which may increase the risk of Paragonimus infection. These demonstrate the importance of accurate species delimitation and consideration of the effect of landscape on Paragonimus in the proposal of adequate control strategies. However, other landscape variables cannot be discarded, including temperature, rainfall regime, and spatial scale (local, landscape and regional). These additional variables were not explored here, and should be considered in future studies.

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.

A taxonomist's nightmare - Cryptic diversity in Caribbean intertidal arthropods (Arachnida, Acari, Oribatida)

There has been a long controversy about what defines a species and how to delimitate them which resulted in the existence of more than two dozen different species concepts. Recent research on so-called "cryptic species" heated up this debate as some scientists argue that these cryptic species are only a result of incompatible species concepts. While this may be true, we should keep in mind that all concepts are nothing more than human constructs and that the phenomenon of high phenotypic similarity despite reproductive isolation is real. To investigate and understand this phenomenon it is important to classify and name cryptic species as it allows to communicate them with other fields of science that use Linnaean binomials. To provide a common framework for the description of cryptic species, we propose a possible protocol of how to formally name and describe these taxa in practice. The most important point of this protocol is to explain which species concept was used to delimitate the cryptic taxon. As a model, we present the case of the allegedly widespread Caribbean intertidal mite Thalassozetes barbara, which in fact consists of seven phenotypically very similar but genetically distinct species. All species are island or short-range endemics with poor dispersal abilities that have evolved in geographic isolation. Stabilizing selection caused by the extreme conditions of the intertidal environment is suggested to be responsible for the morphological stasis of this cryptic species complex.

DNA barcoding of Austrian snow scorpionflies (Mecoptera, Boreidae) reveals potential cryptic diversity in Boreus westwoodi

Background

Snow scorpionflies (genus Boreus) belong to a family of Mecoptera, Boreidae, that has been vastly neglected by entomological researchers due to their shift in seasonality to the winter months. Their activity during this time is regarded as a strategy for predator avoidance and regular sightings on snow fields suggest that this also facilitates dispersal. However, many aspects about snow scorpionflies, especially systematics, taxonomy, distribution of species, phylogenetics and phylogeography have remained fairly unexplored until today. In this study, we fill some of these gaps by generating a reference DNA barcode database for Austrian snow scorpionflies in the frame of the Austrian Barcode of Life initiative and by characterising morphological diversity in the study region.

Methods

Initial species assignment of all 67 specimens was based on male morphological characters previously reported to differ between Boreus species and, for females, the shape of the ovipositor. DNA barcoding of the mitochondrial cytochrome c oxidase subunit 1 (COI) gene was carried out for all 67 samples and served as a basis for BIN assignment, genetic distance calculations, as well as alternative species delimitation analyses (ABGD, GMYC, bGMYC, bPTP) and a statistical parsimony network to infer phylogenetic relationships among individual samples/sampling sites.

Results

Morphological investigations suggested the presence of both Boreus hyemalis and Boreus westwoodi in Austria. DNA barcoding also separated the two species, but resulted in several divergent clades, the paraphyly of B. westwoodi in Austria, and high levels of phylogeographic structure on a small geographic scale. Even though the different molecular species delimitation methods disagreed on the exact number of species, they unequivocally suggested the presence of more than the traditionally recognized two Boreus species in Austria, thus indicating potential cryptic species within the genus Boreus in general and especially in B. westwoodi.

DNA Barcoding of Marine Mollusks Associated with Corallina officinalis Turfs in Southern Istria (Adriatic Sea)

Presence of mollusk assemblages was studied within red coralligenous algae Corallina officinalis L. along the southern Istrian coast. C. officinalis turfs can be considered a biodiversity reservoir, as they shelter numerous invertebrate species. The aim of this study was to identify mollusk species within these settlements using DNA barcoding as a method for detailed identification of mollusks. Nine locations and 18 localities with algal coverage range above 90% were chosen at four research areas. From 54 collected samples of C. officinalis turfs, a total of 46 mollusk species were identified. Molecular methods helped identify 16 gastropod, 14 bivalve and one polyplacophoran species. COI sequences for two bivalve species (Musculus cf. costulatus (Risso, 1826) and Gregariella semigranata (Reeve, 1858)) and seven gastropod species (Megastomia winfriedi Peñas & Rolán, 1999, Eatonina sp. Thiele, 1912, Eatonina cossurae (Calcara, 1841), Crisilla cf. maculata (Monterosato, 1869), Alvania cf. carinata (da Costa, 1778), Vitreolina antiflexa (Monterosato, 1884) and Odostomia plicata (Montagu, 1803)) represent new BINs in BOLD database. This study contributes to new findings related to the high biodiversity of mollusks associated with widespread C. officinalis settlements along the southern coastal area of Istria.

Use of species delimitation approaches to tackle the cryptic diversity of an assemblage of high Andean butterflies (Lepidoptera: Papilionoidea)

Cryptic biological diversity has generated ambiguity in taxonomic and evolutionary studies. Single-locus methods and other approaches for species delimitation are useful for addressing this challenge, enabling the practical processing of large numbers of samples for identification and inventory purposes. This study analyzed an assemblage of high Andean butterflies using DNA barcoding and compared the identifications based on the current morphological taxonomy with three methods of species delimitation (automatic barcode gap discovery, generalized mixed Yule coalescent model, and Poisson tree processes). Sixteen potential cryptic species were recognized using these three methods, representing a net richness increase of 11.3% in the assemblage. A well-studied taxon of the genus Vanessa, which has a wide geographical distribution, appeared with the potential cryptic species that had a higher genetic differentiation at the local level than at the continental level. The analyses were useful for identifying the potential cryptic species in Pedaliodes and Forsterinaria complexes, which also show differentiation along altitudinal and latitudinal gradients. This genetic assessment of an entire assemblage of high Andean butterflies (Papilionoidea) provides baseline information for future research in a region characterized by high rates of endemism and population isolation.

Addressing the Linnean shortfall in a cryptic species complex

Biodiversity is being lost at an alarming rate, but the rate of this loss is likely to be underestimated as a result of a deficit in taxonomic knowledge (i.e. the Linnean shortfall). This knowledge gap is more extensive for morphologically indistinct taxa. The advancement of molecular techniques and delimitation methods has facilitated the identification of such cryptic species, but a majority of these taxa remain undescribed. To investigate the effects of taxonomic uncertainty on understanding of biodiversity, we applied the general lineage concept of species to an amphipod species complex, the Gammaruslacustris lineage that occupies springs of the northern Chihuahuan Desert, which is emerging in contemporary times. We investigated species boundaries using a validation-based approach and examined genetic structure of the lineage using a suite of microsatellite markers to identify independently evolving metapopulations. Our results show that each spring contains a genetically distinct population that is geographically isolated from other springs, suggesting evolutionary independence and status as separate species. Additionally, we observed subtle interspecific morphological variation among the putative species. We used multiple lines of evidence to formally describe four new species (Gammarus langi sp. nov., G. percalacustris sp. nov., G. colei sp. nov. and G. malpaisensis sp. nov.) endemic to the northern Chihuahuan Desert. Cryptic speciation is likely to be high in other aquatic taxa within these ecosystems, and across arid landscapes throughout North America and elsewhere, suggesting that the magnitude of the Linnean shortfall is currently underestimated in desert springs worldwide.

New insights into the genetic diversity of the stone crayfish: taxonomic and conservation implications

Background

Austropotamobius torrentium is a freshwater crayfish species native to central and south-eastern Europe, with an intricate evolutionary history and the highest genetic diversity recorded in the northern-central Dinarides (NCD). Its populations are facing declines, both in number and size across its entire range. By extanding current knowledge on the genetic diversity of this species, we aim to assist conservation programmes. Multigene phylogenetic analyses were performed using different divergence time estimates based on mitochondrial and, for the first time, nuclear DNA markers on the largest data set analysed so far. In order to reassess taxonomic relationships within this species we applied several species delimitation methods and studied the meristic characters with the intention of finding features that would clearly separate stone crayfish belonging to different phylogroups.

Results

Our results confirmed the existence of high genetic diversity within A. torrentium, maintained in divergent phylogroups which have their own evolutionary dynamics. A new phylogroup in the Kordun region belonging to NCD has also been discovered. Due to the incongruence between implemented species delimitation approaches and the lack of any morphological characters conserved within lineages, we are of the opinion that phylogroups recovered on mitochondrial and nuclear DNA are cryptic subspecies and distinct evolutionary significant units.

Conclusions

Geographically and genetically isolated phylogroups represent the evolutionary legacy of A. torrentium and are highly relevant for conservation due to their evolutionary distinctiveness and restricted distribution.

Population genetics of the brooding coral Seriatopora hystrix reveals patterns of strong genetic differentiation in the Western Indian Ocean

Coral reefs provide essential goods and services but are degrading at an alarming rate due to local and global anthropogenic stressors. The main limitation that prevents the implementation of adequate conservation measures is that connectivity and genetic structure of populations are poorly known. Here, the genetic diversity and connectivity of the brooding scleractinian coral Seriatopora hystrix were assessed at two scales by genotyping ten microsatellite markers for 356 individual colonies. S. hystrix showed high differentiation, both at large scale between the Red Sea and the Western Indian Ocean (WIO), and at smaller scale along the coast of East Africa. As such high levels of differentiation might indicate the presence of more than one species, a haploweb analysis was conducted with the nuclear marker ITS2, confirming that the Red Sea populations are genetically distinct from the WIO ones. Based on microsatellite analyses three groups could be distinguished within the WIO: (1) northern Madagascar, (2) south-west Madagascar together with one site in northern Mozambique (Nacala) and (3) all other sites in northern Mozambique, Tanzania and Kenya. These patterns of restricted connectivity could be explained by the short pelagic larval duration of S. hystrix, and/or by oceanographic factors, such as eddies in the Mozambique Channel (causing larval retention in northern Madagascar but facilitating dispersal from northern Mozambique towards south-west Madagascar). This study provides an additional line of evidence supporting the conservation priority status of the Northern Mozambique Channel and should inform coral reef management decisions in the region.

Cryptic Diversity in the Monotypic Neotropical Micromoth Genus Angelabella (Lepidoptera: Gracillariidae) in the Peru-Chile Desert

Simple Summary

The Neotropical Region harbors a highly diverse and poorly known fauna of leaf miners of the micromoth family Gracillariidae (Lepidoptera). Angelabella is a genus of Gracillariidae whose geographic range is restricted to a few valleys of the arid environments of the Peru-Chile desert. Only one species is currently included in this genus. The aims of this study were to explore the geographic range, determine the spatial distribution of mitochondrial lineages, and test lineage conspecificity hypotheses in Angelabella. The spatial distribution of genetic diversity indicated four spatial clusters, three of which are north of the previously known geographic range. These groups were defined as different species by four species delimitation methods. These results suggest that Angelabella harbors at least four morphologically cryptic species with restricted, not overlapping geographic ranges. This study shows that adequate single locus sequence analysis can be useful to discover surprising biodiversity patterns in underexplored environments, providing the base to plan further studies involving little-known organisms.

Abstract

Angelabella (Lepidoptera: Gracillariidae: Oecophyllembiinae) is considered a monotypic Neotropical genus of leaf miner micromoths known only from a few valleys of the arid environments of the Peru-Chile desert, particularly the southernmost part of Peru and northernmost part of Chile (type locality), where natural populations of its primary host plant occur. The geographic distribution of potential host plants provides a scenario for a wider range for this micromoth genus. The aims of this study were to explore the geographic range of Angelabella, determine the spatial distribution of mitochondrial lineages, and test lineage conspecificity hypotheses. The spatial distribution of genetic diversity indicated the presence of four spatial clusters, three of which are north of the previously known geographic range. Genetic distances were 0.2–0.8% and 3.6–8.3% (K2P) between haplotypes of the same and different spatial clusters, respectively. Phylogenetic relationships indicated reciprocal monophyly among the four spatial clusters, suggesting that allopatric differentiation processes have governed the recent history of Angelabella in these arid environments. These groups were defined as different species by four species delimitation methods, suggesting that Angelabella is not a monotypic genus, but harbors at least four morphologically cryptic allopatric species with restricted geographic ranges, including the type species and three candidate species.

A little frog leaps a long way: compounded colonizations of the Indian Subcontinent discovered in the tiny Oriental frog genus Microhyla (Amphibia: Microhylidae)

Frogs of the genus Microhyla include some of the world's smallest amphibians and represent the largest radiation of Asian microhylids, currently encompassing 50 species, distributed across the Oriental biogeographic region. The genus Microhyla remains one of the taxonomically most challenging groups of Asian frogs and was found to be paraphyletic with respect to large-sized fossorial Glyphoglossus. In this study we present a time-calibrated phylogeny for frogs in the genus Microhyla, and discuss taxonomy, historical biogeography, and morphological evolution of these frogs. Our updated phylogeny of the genus with nearly complete taxon sampling includes 48 nominal Microhyla species and several undescribed candidate species. Phylogenetic analyses of 3,207 bp of combined mtDNA and nuDNA data recovered three well-supported groups: the Glyphoglossus clade, Southeast Asian Microhyla II clade (includes M. annectens species group), and a diverse Microhyla I clade including all other species. Within the largest major clade of Microhyla are seven well-supported subclades that we identify as the M. achatina, M. fissipes, M. berdmorei, M. superciliaris, M. ornata, M. butleri, and M. palmipes species groups. The phylogenetic position of 12 poorly known Microhyla species is clarified for the first time. These phylogenetic results, along with molecular clock and ancestral area analyses, show the Microhyla-Glyphoglossus assemblage to have originated in Southeast Asia in the middle Eocene just after the first hypothesized land connections between the Indian Plate and the Asian mainland. While Glyphoglossus and Microhyla II remained within their ancestral ranges, Microhyla I expanded its distribution generally east to west, colonizing and diversifying through the Cenozoic. The Indian Subcontinent was colonized by members of five Microhyla species groups independently, starting with the end Oligocene-early Miocene that coincides with an onset of seasonally dry climates in South Asia. Body size evolution modeling suggests that four groups of Microhyla have independently achieved extreme miniaturization with adult body size below 15 mm. Three of the five smallest Microhyla species are obligate phytotelm-breeders and we argue that their peculiar reproductive biology may be a factor involved in miniaturization. Body size increases in Microhyla-Glyphoglossus seem to be associated with a burrowing adaptation to seasonally dry habitats. Species delimitation analyses suggest a vast underestimation of species richness and diversity in Microhyla and reveal 15-33 undescribed species. We revalidate M. nepenthicola, synonymize M. pulverata with M. marmorata, and provide insights on taxonomic statuses of a number of poorly known species. Further integrative studies, combining evidence from phylogeny, morphology, advertisement calls, and behavior will result in a better systematic understanding of this morphologically cryptic radiation of Asian frogs.

Cryptic diversity and species boundaries within the Paragalago zanzibaricus species complex

The recently described genus Paragalago is a complex of several nocturnal and morphologically cryptic species distributed in the forests of eastern Africa. Species diversity within this genus has been mainly described using species-specific differences in their loud calls. However, molecular data are still lacking for this group and species boundaries remain unclear. In this study, we explore species diversity within the zanzibaricus-complex using a combination of mitochondrial and nuclear data and comparing multiple species delimitation methods. Our results consistently support the existence of three independent lineages, P. cocos, P. zanzibaricus, and P. granti, confirming previous hypotheses based on vocal data. We conclude that these three lineages represent valid cryptic species and we hypothesize that speciation within this complex was characterized by cycles of forest expansion and contraction in the Plio-Pleistocene.

The tree that hides the forest: cryptic diversity and phylogenetic relationships in the Palaearctic vector Obsoletus/Scoticus Complex (Diptera: Ceratopogonidae) at the European level

Background

Culicoides obsoletus is an abundant and widely distributed Holarctic biting midge species, involved in the transmission of bluetongue virus (BTV) and Schmallenberg virus (SBV) to wild and domestic ruminants. Females of this vector species are often reported jointly with two morphologically very close species, C. scoticus and C. montanus, forming the Obsoletus/Scoticus Complex. Recently, cryptic diversity within C. obsoletus was reported in geographically distant sites. Clear delineation of species and characterization of genetic variability is mandatory to revise their taxonomic status and assess the vector role of each taxonomic entity. Our objectives were to characterize and map the cryptic diversity within the Obsoletus/Scoticus Complex.

Methods

Portion of the cox1 mitochondrial gene of 3763 individuals belonging to the Obsoletus/Scoticus Complex was sequenced. Populations from 20 countries along a Palaearctic Mediterranean transect covering Scandinavia to Canary islands (North to South) and Canary islands to Turkey (West to East) were included. Genetic diversity based on cox1 barcoding was supported by 16S rDNA mitochondrial gene sequences and a gene coding for ribosomal 28S rDNA. Species delimitation using a multi-marker methodology was used to revise the current taxonomic scheme of the Obsoletus/Scoticus Complex.

Results

Our analysis showed the existence of three phylogenetic clades (C. obsoletus clade O2, C. obsoletus clade dark and one not yet named and identified) within C. obsoletus. These analyses also revealed two intra-specific clades within C. scoticus and raised questions about the taxonomic status of C. montanus.

Conclusions

To our knowledge, our study provides the first genetic characterization of the Obsoletus/Scoticus Complex on a large geographical scale and allows a revision of the current taxonomic classification for an important group of vector species of livestock viruses in the Palaearctic region.

Molecular phylogeny and species delimitation of the genus Dicerapanorpa (Mecoptera: Panorpidae)

Given that species is the fundamental unit in systematic biology, rigorous species delimitation is crucial for taxonomic studies, yet routine species delimitation remains an ongoing challenge in the taxonomic practice of insects. The two-horned scorpionfly Dicerapanorpa is a small genus in Panorpidae (Mecoptera) endemic to the Qinling-Bashan and Hengduan mountains, a biodiversity hotspot. However, species of Dicerapanorpa are difficult to delineate owing to marked intraspecific variation and interspecific similarity. Here, we investigate the diversity and species boundaries of Dicerapanorpa using an integrative approach based on DNA barcoding, morphological, geometric morphometric and molecular phylogenetic analyses. This integrative analyses confirmed the 13 described species of Dicerapanorpa and revealed three new species: Dicerapanorpa lativalva sp. nov., Dicerapanorpa hualongshana sp. nov. and Dicerapanorpa minshana sp. nov. Most molecular operational taxonomic units are in congruence with morphological clusters. Possible reasons for several discordances in Dicerapanorpa are tentatively discussed.

DNA barcoding of the rodent genus Oligoryzomys (Cricetidae: Sigmodontinae): mitogenomic-anchored database and identification of nuclear mitochondrial translocations (Numts)

DNA barcoding has become a standard method for species identification in taxonomically complex groups. An important step of the barcoding process is the construction of a library of voucher-based material that was properly identified by independent methods, free of inaccurate identification, and paralogs. We provide here a cytochrome oxidase I (mt-Co1) DNA barcode database for species of the genus Oligoryzomys, based on type material and karyotyped specimens, and anchored on the mitochondrial genome of one species of Oligoryzomys, O. stramineus. To evaluate the taxonomic determination of new COI sequences, we assessed species intra/interspecific genetic distances (barcode gap), performed the General Mixed Yule Coalescent method (GMYC) for lineages' delimitation, and identified diagnostic nucleotides for each species of Oligoryzomys. Phylogenetic analyses of Oligoryzomys were performed on 2 datasets including 14 of the 23 recognized species of this genus: a mt-Co1 only matrix, and a concatenated matrix including mt-Co1, cytochrome b (mt-Cytb), and intron 7 of the nuclear fibrinogen beta chain gene (i7Fgb). We recovered nuclear-mitochondrial translocated (Numts) pseudogenes on our samples and identified several published sequences that are cases of Numts. We analyzed the rate of non-synonymous and synonymous substitution, which were higher in Numts in comparison to mtDNA sequences. GMYC delimitations and DNA barcode gap results highlight the need for further work that integrate molecular, karyotypic, and morphological analyses, as well as additional sampling, to tackle persistent problems in the taxonomy of Oligoryzomys.

Cryptic Diversity, but to What Extent? Discordance Between Single-Locus Species Delimitation Methods Within Mainland Anoles (Squamata: Dactyloidae) of Northern Central America

Single-locus molecular barcoding is a useful method for identifying overlooked and undescribed biodiversity, providing the groundwork for further systematic study and taxonomic investigation. A variety of methods for delimiting species from barcoding libraries have been developed and applied, allowing for rapid estimates of species diversity in a broad range of taxa. However, tree-based and distance-based analyses can infer different group assignments, potentially over- or underestimating the number of putative species groups. Here, we explore diversity of mainland species of anole lizards from the Chortís Block biogeographical province of northern Central America using a DNA barcoding approach, generating and analyzing cytochrome oxidase subunit I (COI) sequences for over 400 samples assignable to 33 of 38 (86.8%) native and one introduced mainland species. We subsequently tested the effects different models of nucleotide substitution, different species-delimitation algorithms, and reducing our dataset had on species delimitation estimates. We performed of two distance-based (ABGD, RESL) and three tree-based (bPTP, mPTP, GMYC) analyses on both the full dataset and a dataset consisting only of unique halotypes. From 34 nominal taxa, analyses of the full dataset recovered between 34 and 64 operational taxonomic units (OTUs), while analyses of the reduced dataset inferred between 36 and 59. Reassigning individuals to either mPTP-inferred or ABGD clustered (7.2% threshold) groups improved the detection of a barcoding gap across three different models of nucleotide substitution, removing overlap between intra- and interspecific distances. Our results highlight the underestimated diversity of mainland Chortís Block anoles, but the lack of congruence between analyses demonstrates the importance of considering multiple analytical methods when dealing with single-locus datasets. We recommend future studies consider the effects of different models of nucleotide substitution on proposed barcoding gaps, as well as the effect reducing a dataset to unique haplotypes may have on proposed diversity estimates.

Molecular phylogeographic analyses and species delimitations reveal that Leopoldamys edwardsi (Rodentia: Muridae) is a species complex

Leopoldamys edwardsi is a species with wide distribution ranges in southern China but is not discussed in studies on geographic variation and species differentiation. We used 2 mitochondrial (Cytb, CO1) and 3 nuclear (GHR, IRBP and RAG1) genes to clarify species phylogeography and geographical differentiation. Maximum likelihood (ML) and Bayesian phylogenetic inference (BI) trees consistently indicated that L. edwardsi is a species complex containing 3 main lineages with high Kimura-2-parameter (K2P) divergences (i.e. lineages L_N , L_S and L_HN ) found in the northern and southern China and Hainan Island, respectively. The 3 species delimitation methods, automated barcoding gap discovery, Bayesian poisson tree process analysis and Bayesian phylogenetics and phylogeography, consistently supported the existence of cryptic species. Divergence times among the main lineages were inferred to be during the Pleistocene, with L_HN /L_S split at 1.33 Ma and L_N /(L_HN +L_S ) at 2.61 Ma; the diversifications of L. edwardsi complex might be caused by the rapid uplifts of Tibetan Plateau, paleoclimate change and complex topography. The divergence between L_HN and L_S was probably related to the separation of Hainan Island from the mainland via the formation of the Qiongzhou Strait. Lineages L_N and (L_S +L_HN ) likely diverged due to the Wuyi-Nanling mountain range forming a dispersal barrier. Our results suggested that L. edwardsi complex contains at least 3 distinct species: L_HN represents L. hainanensis, endemic to Hainan Island and previously considered as a subspecies L. e. hainanensis; L_S represents a cryptic species distributed throughout the southern Chinese continent; and L_N represents the nominotypical species L. edwardsi.

Cox1 barcoding versus multilocus species delimitation: validation of two mite species with contrasting effective population sizes

BACKGROUND

The cox1-barcoding approach is currently extensively used for high-throughput species delimitation and discovery. However, this method has several limitations, particularly when organisms have large effective population sizes. Paradoxically, most common, abundant, and widely distributed species may be misclassified by this technique.

RESULTS

We conducted species delimitation analyses for two host-specific lineages of scab mites of the genus Caparinia, having small population sizes. Cox1 divergence between these lineages was high (7.4-7.8%) while that of nuclear genes was low (0.06-0.53%). This system was contrasted with the medically important American house dust mite, Dermatophagoides farinae, a globally distributed species with very large population size. This species has two distinct, sympatric cox1 lineages with 4.2% divergence. We tested several species delimitation algorithms PTP, GMYC, ABGD, BPP, STACEY and PHRAPL, which inferred different species boundaries for these entities. Notably, STACEY recovered the Caparinia lineages as two species and D. farinae as a single species. BPP agreed with these results when the prior on ancestral effective population sizes was set to expected values, although delimitation of Caparinia was still equivocal. No other cox1 species delimitation algorithms inferred D. farinae as a single species, despite the fact that the nuclear CPW2 gene shows some evidence for introgression between the cox1 groups. This indicates that the cox1-barcoding approach may result in excessive species splitting.

CONCLUSIONS

Our research highlights the importance of using nuclear genes and demographic characteristics to infer species boundaries rather than relying on a single-gene barcoding approach, particularly for putative species having large effective population sizes.

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Uncovering hidden specific diversity of Andean glassfrogs of the Centrolene buckleyi species complex (Anura: Centrolenidae)

The glassfrog Centrolene buckleyi has been recognized as a species complex. Herein, using coalescence-based species delimitation methods, we evaluate the specific diversity within this taxon. Four coalescence approaches (generalized mixed Yule coalescents, Bayesian general mixed Yule-coalescent, Poisson tree processes, and Bayesian Poisson tree processes) were consistent with the delimitation results, identifying four lineages within what is currently recognized as C. buckleyi. We propose three new candidate species that should be tested with nuclear markers, morphological, and behavioral data. In the meantime, for conservation purposes, candidate species should be considered evolutionary significant units, in light of observed population crashes in the C. buckleyi species complex. Finally, our results support the validity of C. venezuelense, formerly considered as a subspecies of C. buckleyi.

Comparison of Methods for Molecular Species Delimitation Across a Range of Speciation Scenarios

Species are fundamental units in biological research and can be defined on the basis of various operational criteria. There has been growing use of molecular approaches for species delimitation. Among the most widely used methods, the generalized mixed Yule-coalescent (GMYC) and Poisson tree processes (PTP) were designed for the analysis of single-locus data but are often applied to concatenations of multilocus data. In contrast, the Bayesian multispecies coalescent approach in the software Bayesian Phylogenetics and Phylogeography (BPP) explicitly models the evolution of multilocus data. In this study, we compare the performance of GMYC, PTP, and BPP using synthetic data generated by simulation under various speciation scenarios. We show that in the absence of gene flow, the main factor influencing the performance of these methods is the ratio of population size to divergence time, while number of loci and sample size per species have smaller effects. Given appropriate priors and correct guide trees, BPP shows lower rates of species overestimation and underestimation, and is generally robust to various potential confounding factors except high levels of gene flow. The single-threshold GMYC and the best strategy that we identified in PTP generally perform well for scenarios involving more than a single putative species when gene flow is absent, but PTP outperforms GMYC when fewer species are involved. Both methods are more sensitive than BPP to the effects of gene flow and potential confounding factors. Case studies of bears and bees further validate some of the findings from our simulation study, and reveal the importance of using an informed starting point for molecular species delimitation. Our results highlight the key factors affecting the performance of molecular species delimitation, with potential benefits for using these methods within an integrative taxonomic framework.

The hitchhiker's guide to single-locus species delimitation

Molecular approaches to species delimitation are increasingly used to ascertain the number of species in a sample prior to taxonomic, ecological or physiological studies. Although multilocus approaches are gaining fast in popularity, single-gene methods still predominate in the literature. However, available simulation benchmarks of these methods focus exclusively on species-poor samples and/or tree-based approaches: as a result, travellers in the land of single-locus species delimitation lack a comprehensive "hitchhiker's guide" highlighting the sweet spots and dangers on their road. To fill this gap, we compared the performances of distance-based (ABGD, "automatic barcode gap discovery"), allele sharing-based (haplowebs) and tree-based approaches (GMYC, "generalized mixed Yule-coalescent" and PTP, "Poisson tree processes") to detect interspecific boundaries in samples of 6, 60 and 120 simulated species with various speciation rates, effective population sizes, mutation rates and sampling patterns. We found that all approaches performed poorly when population sizes and speciation rates were large, with haplowebs yielding best results followed by ABGD then tree-based approaches. The latter's error type was mostly oversplitting, whereas ABGD chiefly overlumped and haplowebs leaned either way depending on simulation parameters: such widely divergent error patterns suggest that, if all three types of methods agree, then the resulting delimitation is probably correct. Perfect congruence being quite rare, travellers in search of a one-size-fit-all approach to single-locus species delimitation should forget it; however, our hitchhiker's guide raises hope that such species delimitation's Holy Grail may be found in the relatively uncharted nearby land of multilocus species delimitation.

Formal description of sequence-based voucherless Fungi: promises and pitfalls, and how to resolve them

There is urgent need for a formal nomenclature of sequence-based, voucherless Fungi, given that environmental sequencing has accumulated more than one billion fungal ITS reads in the Sequence Read Archive, about 1,000 times as many as fungal ITS sequences in GenBank. These unnamed Fungi could help to bridge the gap between 115,000 to 140,000 currently accepted and 2.2 to 3.8 million predicted species, a gap that cannot realistically be filled using specimen or culture-based inventories. The Code never aimed at placing restrictions on the nature of characters chosen for taxonomy, and the requirement for physical types is now becoming a constraint on the advancement of science. We elaborate on the promises and pitfalls of sequence-based nomenclature and provide potential solutions to major concerns of the mycological community. Types of sequence-based taxa, which by default lack a physical specimen or culture, could be designated in four alternative ways: (1) the underlying sample ('bag' type), (2) the DNA extract, (3) fluorescent in situ hybridization (FISH), or (4) the type sequence itself. Only (4) would require changes to the Code and the latter would be the most straightforward approach, complying with three of the five principal functions of types better than physical specimens. A fifth way, representation of the sequence in an illustration, has been ruled as unacceptable in the Code. Potential flaws in sequence data are analogous to flaws in physical types, and artifacts are manageable if a stringent analytical approach is applied. Conceptual errors such as homoplasy, intragenomic variation, gene duplication, hybridization, and horizontal gene transfer, apply to all molecular approaches and cannot be used as a specific argument against sequence-based nomenclature. The potential impact of these phenomena is manageable, as phylogenetic species delimitation has worked satisfactorily in Fungi. The most serious shortcoming of sequence-based nomenclature is the likelihood of parallel classifications, either by describing taxa that already have names based on physical types, or by using different markers to delimit species within the same lineage. The probability of inadvertently establishing sequence-based species that have names available is between 20.4 % and 1.5 % depending on the number of globally predicted fungal species. This compares favourably to a historical error rate of about 30 % based on physical types, and this rate could be reduced to practically zero by adding specific provisions to this approach in the Code. To avoid parallel classifications based on different markers, sequence-based nomenclature should be limited to a single marker, preferably the fungal ITS barcoding marker; this is possible since sequence-based nomenclature does not aim at accurate species delimitation but at naming lineages to generate a reference database, independent of whether these lineages represent species, closely related species complexes, or infraspecies. We argue that clustering methods are inappropriate for sequence-based nomenclature; this approach must instead use phylogenetic methods based on multiple alignments, combined with quantitative species recognition methods. We outline strategies to obtain higher-level phylogenies for ITS-based, voucherless species, including phylogenetic binning, 'hijacking' species delimitation methods, and temporal banding. We conclude that voucherless, sequence-based nomenclature is not a threat to specimen and culture-based fungal taxonomy, but a complementary approach capable of substantially closing the gap between known and predicted fungal diversity, an approach that requires careful work and high skill levels.

A hyper-diverse genus of acanthocephalans revealed by tree-based and non-tree-based species delimitation methods: Ten cryptic species of Neoechinorhynchus in Middle American freshwater fishes

The genus Neoechinorhynchus represents a hyper-diverse group of acanthocephalans, parasites of fresh and brackish water fish and freshwater turtles, with approximately 116 species described worldwide. Forty-nine species have been recorded in the Americas, nine of them in Middle America. Even though species delimitation methods using DNA sequences have been rarely used for parasitic helminths, the genetic library for species of Neoechinorhynchus has grown in the past few years, enhancing the possibility of using these methods for inferring evolutionary relationships and for establishing more robust species boundaries. In this study, we used non-tree-based and tree-based methods through a coalescent approach to explore the species limits of specimens of Neoechinorhynchus collected in 57 localities across Middle America. We sequenced a large number of individuals to build a comprehensive dataset for three genes: the mitochondrial cytochrome c oxidase subunit I (352 individuals), the internal transcribed spacers (330 individuals), and the D2 + D3 domains of the large subunit (278 individuals). Several species delimitation methods were implemented, i.e., Automatic Barcode Gap Discovery (ABGD), General Mixed Yule-Coalescent Model (GMYC), Bayesian species delimitation (BPP) and species tree (∗BEAST). Additionally, we conducted a detailed morphological study of the diagnostic traits associated with the proboscis of 184 males and 169 females. Overall, our analyses allowed us to validate nine nominal species of Neoechinorhynchus and to identify 10 additional genetic lineages herein regarded as candidate species. This unexpected genetic diversity and the lack of reliable morphological traits show that the genus Neoechinorhynchus includes a group of cryptic species, at least in Middle America.

Quantification of complex modular architecture in plants

Morphometrics, the assignment of quantities to biological shapes, is a powerful tool to address taxonomic, evolutionary, functional and developmental questions. We propose a novel method for shape quantification of complex modular architecture in thalloid plants, whose extremely reduced morphologies, combined with the lack of a formal framework for thallus description, have long rendered taxonomic and evolutionary studies extremely challenging. Using graph theory, thalli are described as hierarchical series of nodes and edges, allowing for accurate, homologous and repeatable measurements of widths, lengths and angles. The computer program MorphoSnake was developed to extract the skeleton and contours of a thallus and automatically acquire, at each level of organization, width, length, angle and sinuosity measurements. Through the quantification of leaf architecture in Hymenophyllum ferns (Polypodiopsida) and a fully worked example of integrative taxonomy in the taxonomically challenging thalloid liverwort genus Riccardia, we show that MorphoSnake is applicable to all ramified plants. This new possibility of acquiring large numbers of quantitative traits in plants with complex modular architectures opens new perspectives of applications, from the development of rapid species identification tools to evolutionary analyses of adaptive plasticity.

Explaining large mitochondrial sequence differences within a population sample

Mitochondrial DNA sequence is frequently used to infer species' boundaries, as divergence is relatively rapid when populations are reproductively isolated. However, the shared history of a non-recombining gene naturally leads to correlation of pairwise differences, resulting in mtDNA clusters that might be mistaken for evidence of multiple species. There are four distinct processes that can explain high levels of mtDNA sequence difference within a single sample. Here, we examine one case in detail as an exemplar to distinguish among competing hypotheses. Within our sample of tree wētā (Hemideina crassidens; Orthoptera), we found multiple mtDNA haplotypes for a protein-coding region (cytb/ND1) that differed by a maximum of 7.9%. From sequencing the whole mitochondrial genome of two representative individuals, we found evidence of constraining selection. Heterozygotes were as common as expected under random mating at five nuclear loci. Morphological traits and nuclear markers did not resolve the mtDNA groupings of individuals. We concluded that the large differences found among our sample of mtDNA sequences were simply owing to a large population size over an extended period of time allowing an equilibrium between mutation and drift to retain a great deal of genetic diversity within a single species.