Exploring species boundaries with multiple genetic loci using empirical data from non‐biting midges

Over the past decade, molecular approaches to species delimitation have seen rapid development. However, species delimitation based on a single locus, for example, DNA barcodes, can lead to inaccurate results in cases of recent speciation and incomplete lineage sorting. Here, we compare the performance of Automatic Barcode Gap Discovery (ABGD), Bayesian Poisson tree processes (PTP), networks, generalized mixed Yule coalescent (GMYC) and Bayesian phylogenetics and phylogeography (BPP) models to delineate cryptic species previously detected by DNA barcodes within Tanytarsus (Diptera: Chironomidae) non‐biting midges. We compare the results from analyses of one mitochondrial (cytochrome c oxidase subunit I [COI]) and three nuclear (alanyl‐tRNA synthetase 1 [AATS1], carbamoyl phosphate synthetase 1 [CAD1] and 6‐phosphogluconate dehydrogenase [PGD]) protein‐coding genes. Our results show that species delimitation based on multiple nuclear DNA markers is largely concordant with morphological variation and delimitations using a single locus, for example, the COI barcode. However, ABGD, GMYC, PTP and network models led to conflicting results based on a single locus and delineate species differently than morphology. Results from BPP analyses on multiple loci correspond best with current morphological species concept. In total, 10 lineages of the Tanytarsus curticornis species complex were uncovered. Excluding a Norwegian population of Tanytarsus brundini which might have undergone recent hybridization, this suggests six hitherto unrecognized species new to science. Five distinct species are well supported in the Tanytarsus heusdensis species complex, including two species new to science.     

Species‐delimitation and phylogenetic analyses of some cosmopolitan species of Hypnea (Rhodophyta) reveal synonyms and misapplied names to H. cervicornis,including a new species from Brazil

Hypnea has an intricate nomenclatural history due to a wide pantropical distribution and considerable morphological variation. Recent molecular studies have provided further clarification on the systematics of the genus; however, species of uncertain affinities remain due to flawed taxonomic identification. Detailed analyses coupled with literature review indicated a strong relationship among H. aspera, H. cervicornis, H. flexicaulis, and H. tenuis, suggesting a need for further taxonomic studies. Here, we analyzed sequences from two molecular markers (COI‐5P and rbcL) and performed several DNA‐based delimitation methods (mBGD, ABGD, SPN, PTP and GMYC). These molecular approaches were contrasted with morphological and phylogenetic evidence from type specimens and/or topotype collections of related species under a conservative approach. Our results demonstrate that H. aspera and H. flexicaulis represent heterotypic synonyms of H. cervicornis and indicate the existence of a misidentified Hypnea species, widely distributed on the Brazilian coast, described here as a new species: H. brasiliensis. Finally, inconsistencies observed among our results based on six different species delimitation methods evidence the need for adequate sampling and marker choice for different methods.     

Development of high‐resolution DNA barcodes for Dioscorea species discrimination and phylogenetic analysis

The genus Dioscorea is widely distributed in tropical and subtropical regions, and is economically important in terms of food supply and pharmaceutical applications. However, DNA barcodes are relatively unsuccessful in discriminating between Dioscorea species, with the highest discrimination rate (23.26%) derived from matK sequences. In this study, we compared genic and intergenic regions of three Dioscorea chloroplast genomes and found that the density of SNPs and indels in intergenic sites was about twice and seven times higher than that of SNPs and indels in the genic regions, respectively. A total of 52 primer pairs covering highly variable regions were designed and seven pairs of primers had 80%–100% PCR success rate. PCR amplicons of 73 Dioscorea individuals and assembled sequences of 47 Dioscorea SRAs were used for estimating intraspecific and interspecific divergence for the seven loci: The rpoB‐trnC locus had the highest interspecific divergence. Automatic barcoding gap discovery (ABGD), Poisson tree processes (PTP), and generalized mixed Yule coalescence (GMYC) analysis were applied for species delimitation based on the seven loci and successfully identified the majority of species, except for species in the Enantiophyllum section. Phylogenetic analysis of 51 Dioscorea individuals (28 species) showed that most individuals belonging to the same species tended to cluster in the same group. Our results suggest that the variable loci derived from comparative analysis of plastid genome sequences could be good DNA barcode candidates for taxonomic analysis and species delimitation.     

Phylogeny and species reassessment of Hyalopterus (Aphididae,Aphidinae)

The broadly distributed genus Hyalopterus currently comprises three formally recognized species that are highly similar morphologically and hence difficult to be identified with certainty. This group has undergone multiple revisions in the past century, but none of these has assessed species from Asia, which has hampered our understanding of the species diversity within this genus. Based on a comprehensive data set from morphological data and host-associated data, and by coalescent-based delimitation approaches, the Hyalopterus species boundaries, distribution and diversity were clarified here to further reveal the composition of the species. Two single-locus (ML-GMYC and mPTP) and two multilocus (BPP and STACEY) delimitation methods were conducted based on extensive sampling. Then, the phylogenetic relationships and morphological divergence were assessed. Our data strongly supported that the number of recognized species in Hyalopterus had likely been underestimated. The phylogenetic analyses recovered four major clades, which corresponded to distinct host-plant preferences. Also, the morphological analyses showed significant differentiation for only one of the newly recognized candidate species uncovered by the delimitation approaches, suggesting the existence of at least two independent evolutionary lineages within Hyalopterus arundiniformis, which showed different patterns of host association. Moreover, based on our data, the taxonomic misidentification of H. arundiniformis in China was corrected here. This study lays the groundwork for the thorough taxonomic revision of Hyalopterus and for future evolutionary studies and underlines the importance of an integrated framework for species determination.     

An empirical comparison of character‐based and coalescent‐based approaches to species delimitation in a young avian complex

The process of discovering species is a fundamental responsibility of systematics. Recently, there has been a growing interest in coalescent‐based methods of species delimitation aimed at objectively identifying species early in the divergence process. However, few empirical studies have compared these new methods with character‐based approaches for discovering species. In this study, we applied both a character‐based and a coalescent‐based approaches to delimit species in a closely related avian complex, the light‐vented/Taiwan bulbul (Pycnonotus sinensis/Pycnonotus taivanus). Population aggregation analyses of plumage, mitochondrial and 13 nuclear intron character data sets produced conflicting species hypotheses with plumage data suggesting three species, mitochondrial data suggesting two species, and nuclear intron data suggesting one species. Such conflict is expected among recently diverged species, and by integrating all sources of data, we delimited three species verified with independently congruent character evidence as well as a more weakly supported fourth species identified by a single character. Attempts to validate species hypothesis using Bayesian Phylogenetics and Phylogeography (BPP), a coalescent‐based method of species delimitation, revealed several issues that can seemingly affect statistical support for species recognition. We found that θ priors had a dramatic impact on speciation probabilities, with lower values consistently favouring splitting and higher values consistently favouring lumping. More resolved guide trees also resulted in overall higher speciation probabilities. Finally, we found suggestive evidence that BPP is sensitive to the divergent effects of nonrandom mating caused by intraspecific processes such as isolation‐with‐distance, and therefore, BPP may not be a conservative method for delimiting independently evolving population lineages. Based on these concerns, we questioned the reliability of BPP results and based our conclusions about species limits exclusively on character data.     

A Good Compromise: Rapid and Robust Species Proxies for Inventorying Biodiversity Hotspots Using the Terebridae (Gastropoda: Conoidea)

Devising a reproducible approach for species delimitation of hyperdiverse groups is an ongoing challenge in evolutionary biology. Speciation processes combine modes of passive and adaptive trait divergence requiring an integrative taxonomy approach to accurately generate robust species hypotheses. However, in light of the rapid decline of diversity on Earth, complete integrative approaches may not be practical in certain species-rich environments. As an alternative, we applied a two-step strategy combining ABGD (Automated Barcode Gap Discovery) and Klee diagrams, to balance speed and accuracy in producing primary species hypotheses (PSHs). Specifically, an ABGD/Klee approach was used for species delimitation in the Terebridae, a neurotoxin-producing marine snail family included in the Conoidea. Delimitation of species boundaries is problematic in the Conoidea, as traditional taxonomic approaches are hampered by the high levels of variation, convergence and morphological plasticity of shell characters. We used ABGD to analyze gaps in the distribution of pairwise distances of 454 COI sequences attributed to 87 morphospecies and obtained 98 to 125 Primary Species Hypotheses (PSHs). The PSH partitions were subsequently visualized as a Klee diagram color map, allowing easy detection of the incongruences that were further evaluated individually with two other species delimitation models, General Mixed Yule Coalescent (GMYC) and Poisson Tree Processes (PTP). GMYC and PTP results confirmed the presence of 17 putative cryptic terebrid species in our dataset. The consensus of GMYC, PTP, and ABGD/Klee findings suggest the combination of ABGD and Klee diagrams is an effective approach for rapidly proposing primary species proxies in hyperdiverse groups and a reliable first step for macroscopic biodiversity assessment.     

Multilocus species delimitation in Mesoamerican Scaptotrigona stingless bees (Apidae: Meliponini) supports the existence of cryptic species

Accelerating taxonomic knowledge and making accurate species identifications are critically important given the current biodiversity crisis, particularly in biodiversity hotspots such as Mesoamerica. Objective species delimitation that reduces investigator‐driven bias is fundamental to the establishment of appropriate conservation strategies, above all in managed species. Previous morphological and molecular studies on three managed stingless bee species of the genus Scaptotrigona distributed in Mexico (S. mexicana, S. pectoralis and S. hellwegeri) suggested that both S. mexicana and S. hellwegeri are cryptic species complexes. Herein we tested species delimitation by analysing sequence information of five markers (two mitochondrial: cox1 and 16S, and three nuclear: ITS1, EF1‐α, ArgK) within a Bayesian coalescent framework to test the putative species. We obtained two different hypotheses using a Generalized Mixed Yule Coalescent (GMYC) model: four (cox1) and six (16S) species. After the species validation step with the Bayesian species‐delimitation analysis (BPP), we suggest that only S. mexicana is a complex of two species with different distribution (along the Pacific and the Atlantic coasts, respectively). We highly recommend avoiding colony exchange between geographical regions in order to conserve the genetic integrity of both taxa.     

Mitochondrial divergence suggests unexpected high species diversity in the opsariichthine fishes (Teleostei: Cyprinidae) and the revalidation of Opsariichthys macrolepis

Opsariichthine (sensu Oceanologi Et Limnologia Sinica, 1982, 13, 293–298) is a cyprinid group consisting of five genera and endemic to East Asia. Previous studies suggested that there may be many possible cryptic species in this group, but this has not been confirmed. In this study, using mitochondrial cyt b sequences on 1,388 samples and 739 haplotypes, we showed very high species diversity within this group. The results showed that phylogenetic relationships of the opsariichthine group were as ([Nipponocypris‐Parazacco‐Candidia] + [Zacco + Opsariichthys]), and there were multiple deep lineages within several species, flagging putative cryptic species. When a 3% genetic distance was used as a threshold for species delimitation, 35 haplogroups were found, nine haplogroups in Candidia‐Parazacco‐Nipponocypris group, six haplogroups in the Zacco group, and 20 haplogroups in the Opsariichthys group. We consider all of them to be putative until determination of distinct species based on the tree topology, geographic distributions, or a combination of both. In addition, two kinds of species delimitation tools, ABGD and PTP, were applied to construct molecular operational taxonomic units (MOTUs). The ABGD method revealed nine MOTUs in Candidia‐Parazacco‐Nipponocypris group, two MOTUs in the Zacco group, and 17 MOTUs in the Opsariichthys group. And the PTP method revealed 10 MOTUs in Candidia‐Parazacco‐Nipponocypris group, 10 MOTUs in the Zacco group, and 29 MOTUs in the Opsariichthys group. Therefore, there should be more species in the opsariichthine group than presently described. Based on the molecular data and morphological characteristics, we proposed Opsariichthys macrolepis as a valid species and described its morphological diagnostic characters.     

Integration of conflict into integrative taxonomy: fitting hybridization in species delimitation of Mesocarabus (Coleoptera: Carabidae)

In species differentiation, characters may not diverge synchronously, and there are also processes that shuffle character states in lineages descendant from a common ancestor. Species are thus expected to show some degree of incongruence among characters; therefore, taxonomic delimitation can benefit from integrative approaches and objective strategies that account for character conflict. We illustrate the potential of exploiting conflict for species delimitation in a study case of ground beetles of the subgenus Carabus (Mesocarabus), where traditional taxonomy does not accurately delimit species. The molecular phylogenies of four mitochondrial and three nuclear genes, cladistic analysis of the aedeagus, ecological niche divergence and morphometry of pronotal shape in more than 500 specimens of Mesocarabus show that these characters are not fully congruent. For these data, a three‐step operational strategy is proposed for species delimitation by (i) delineating candidate species based on the integration of incongruence among conclusive lines of evidence, (ii) corroborating candidate species with inconclusive lines of evidence and (iii) refining a final species proposal based on an integrated characterization of candidate species based on the evolutionary analysis of incongruence. This procedure provided a general understanding of the reticulate process of hybridization and introgression acting on Mesocarabus and generated the hypothesis of seven Mesocarabus species, including two putative hybrid lineages. Our work emphasizes the importance of incorporating critical analyses of character and phylogenetic conflict to infer both the evolutionary history and species boundaries through an integrative taxonomic approach.     

Algorithmic single‐locus species delimitation: effects of sampling effort,variation and nonmonophyly in four methods and 1870 species of beetles

The vast number of undescribed species and the fast rate of biodiversity loss call for new approaches to speed up alpha taxonomy. A plethora of methods for delimiting species or operational taxonomic units (OTUs) based on sequence data have been published in recent years. We test the ability of four delimitation methods (BIN, ABGD, GMYC, PTP) to reproduce established species boundaries on a carefully curated DNA barcode data set of 1870 North European beetle species. We also explore how sampling effort, intraspecific variation, nearest neighbour divergence and nonmonophyly affect the OTU delimitations. All methods produced approximately 90% identity between species and OTUs. The effects of variation and sampling differed between methods. ABGD was sensitive to singleton sequences, while GMYC showed tendencies for oversplitting. The best fit between species and OTUs was achieved using simple rules to find consensus between discordant OTU delimitations. Using several approaches simultaneously allows the methods to compensate for each other's weaknesses. Barcode‐based OTU‐picking is an efficient way to delimit putative species from large data sets where the use of more sophisticated methods based on multilocus or genomic data is not feasible.     

Population and phylogenomic decomposition via genotyping‐by‐sequencing in Australian Pelargonium

Species delimitation has seen a paradigm shift as increasing accessibility of genomic‐scale data enables separation of lineages with convergent morphological traits and the merging of recently diverged ecotypes that have distinguishing characteristics. We inferred the process of lineage formation among Australian species in the widespread and highly variable genus Pelargonium by combining phylogenomic and population genomic analyses along with breeding system studies and character analysis. Phylogenomic analysis and population genetic clustering supported seven of the eight currently described species but provided little evidence for differences in genetic structure within the most widely distributed group that containing P. australe. In contrast, morphometric analysis detected three deep lineages within Australian Pelargonium; with P. australe consisting of five previously unrecognized entities occupying separate geographic ranges. The genomic approach enabled elucidation of parallel evolution in some traits formerly used to delineate species, as well as identification of ecotypic morphological differentiation within recognized species. Highly variable morphology and trait convergence each contribute to the discordance between phylogenomic relationships and morphological taxonomy. Data suggest that genetic divergence among species within the Australian Pelargonium may result from allopatric speciation while morphological differentiation within and among species may be more strongly driven by environmental differences.     

Taxonomy of the long‐tailed mouse Oligoryzomys destructor (Sigmodontinae: Oryzomyini) with the designation of neotypes for Hesperomys destructor Tschudi, 1844 and Hesperomys melanostoma Tschudi, 1844

The genus Oligoryzomys, distributed from southern South America to southern North America, is the most diverse of the tribe Oryzomyini of sigmodontine rodents. Even when 22 species are currently recognized, species boundaries are unclear for several forms. The species Oligoryzomys destructor is one of the least studied species of the genus and is the one with the largest distribution along the Andes (from southern Colombia to northern Bolivia). The species was described without the selection of a holotype and indication of its type locality. In addition, several taxa are regarded as synonyms of O. destructor. These facts are relevant because previous analysis of DNA sequences has shown that O. destructor represents a species complex. Herein, in addition to test the phylogenetic position of O. destructor within the genus Oligoryzomys, we assess patterns of morphological and molecular variation of O. destructor and its associated nominal forms aimed to assess the boundaries of the species. As part of the study, we selected neotypes for Hesperomys destructor and H. melanostoma. At the light of our results, we recognized O. destructor as a species with two subspecies, O. d. destructor and O. d. spodiurus. Also, we discuss the role of Andean rivers, and their different permeability, as allopatric barriers molding the structure of O. destructor.     

A genomic evaluation of taxonomic trends through time in coast horned lizards (genus Phrynosoma)

Determining the boundaries between species and deciding when to describe new species are challenging practices that are particularly difficult in groups with high levels of geographic variation. The coast horned lizards (Phrynosoma blainvillii, Phrynosoma cerroense and P. coronatum) have an extensive geographic distribution spanning many distinctive ecological regions ranging from northern California to the Cape Region of Baja California, Mexico, and populations differ substantially with respect to external morphology across much of this range. The number of taxa recognized in the group has been reevaluated by herpetologists over 20 times during the last 180 years, and typically without the aid of explicit species delimitation methods, resulting in a turbulent taxonomy containing anywhere from one to seven taxa. In this study, we evaluate taxonomic trends through time by ranking 15 of these species delimitation models (SDMs) using coalescent analyses of nuclear loci and SNPs in a Bayesian model comparison framework. Species delimitation models containing more species were generally favoured by Bayesian model selection; however, several three‐species models outperformed some four‐ and five‐species SDMs, and the top‐ranked model, which contained five species, outperformed all SDMs containing six species. Model performance peaked in the 1950s based on marginal likelihoods estimated from nuclear loci and SNPs. Not surprisingly, SDMs based on genetic data outperformed morphological taxonomies when using genetic data alone to evaluate models. The de novo estimation of population structure favours a three‐population model that matches the currently recognized integrative taxonomy containing three species. We discuss why Bayesian model selection might favour models containing more species, and why recognizing more than three species might be warranted.     

Integrative taxonomy of the primitively segmented spider genus Ganthela (Araneae: Mesothelae: Liphistiidae): DNA barcoding gap agrees with morphology

Species delimitation is difficult for taxa in which the morphological characters are poorly known because of the rarity of adult morphs or sexes, and in cryptic species. In primitively segmented spiders, family Liphistiidae, males are often unknown, and female genital morphology – usually species‐specific in spiders – exhibits considerable intraspecific variation. Here, we report on an integrative taxonomic study of the liphistiid genus Ganthela Xu & Kuntner, 2015, endemic to south‐east China, where males are only available for two of the seven morphological species (two known and five undescribed). We obtained DNA barcodes (cytochrome c oxidase subunit I gene, COI) for 51 newly collected specimens of six morphological species and analysed them using five species‐delimitation methods: DNA barcoding gap, species delimitation plugin [P ID(Liberal)], automatic barcode gap discovery (ABGD), generalized mixed Yule‐coalescent model (GMYC), and statistical parsimony (SP). Whereas the first three agreed with the morphology, GMYC and SP indicate several additional species. We used the consensus results to delimit and diagnose six Ganthela species, which in addition to the type species Ganthela yundingensis Xu, 2015, completes the revision of the genus. Although multi‐locus phylogenetic approaches may be needed for complex taxonomic delimitations, our results indicate that even single‐locus analyses based on the COI barcodes, if integrated with morphological and geographical data, may provide sufficiently reliable species delimitation. © 2015 The Linnean Society of London     

HOST SWITCHING PROMOTES DIVERSITY IN HOST‐SPECIALIZED MYCOPARASITIC FUNGI: UNCOUPLED EVOLUTION IN THE BIATOROPSIS‐USNEA SYSTEM

Fungal mycoparasitism—fungi parasitizing other fungi—is a common lifestyle in some basal lineages of the basidiomycetes, particularly within the Tremellales. Relatively nonaggressive mycoparasitic fungi of this group are in general highly host specific, suggesting cospeciation as a plausible speciation mode in these associations. Species delimitation in the Tremellales is often challenging because morphological characters are scant. Host specificity is therefore a great aid to discriminate between species but appropriate species delimitation methods that account for actual diversity are needed to identify both specialist and generalist taxa and avoid inflating or underestimating diversity. We use the Biatoropsis‐Usnea system to study factors inducing parasite diversification. We employ morphological, ecological, and molecular data—methods including genealogical concordance phylogenetic species recognition (GCPSR) and the general mixed Yule‐coalescent (GMYC) model—to assess the diversity of fungi currently assigned to Biatoropsis usnearum. The degree of cospeciation in this association is assessed with two cophylogeny analysis tools (ParaFit and Jane 4.0). Biatoropsis constitutes a species complex formed by at least seven different independent lineages and host switching is a prominent force driving speciation, particularly in host specialists. Combining ITS and nLSU is recommended as barcode system in tremellalean fungi.     

Integrative taxonomy uncovers hidden species diversity in the rheophilic genus Potamometra (Hemiptera: Gerridae)

The genus Potamometra Bianchi, 1896 represents big rheophilic semi-aquatic bugs that typically inhabit middle-altitude mountainous streams. Here, we integrated molecular and morphological data for delimiting species boundaries and understanding the evolutionary history of the genus Potamometra. Twenty-seven complete mitochondrial genomes of Potamometra were sequenced, with samples representing most of the known geographically distributed locations around the Sichuan Basin. The results of different species delimitation methods (ABGD, bPTP, GMYC and BPP) based on the monolocus or multilocus data strongly supported the existence of two cryptic new species (Potamometra anderseni Zheng, Ye & Bu, sp. nov. and Potamometra zhengi Zheng, Ye & Bu, sp. nov.) although more entities were found in the tree-based delimitation methods. The two new species were successfully validated using morphological characters within a detailed taxonomic framework. Phylogenetic analyses supported the reciprocal monophyly of the seven highly node-supported clades, which were matched with the five known species and two new taxa. A novel gene arrangement pattern that two trnF (trnF1 and trnF2) genes separated by an intergenic spacer (IGS) were found in all the species except the sister group of Potamometra berezowskii Bianchi, 1896 and Potamometra linnavuorii Chen, Nieser & Bu, 2016. This gene rearrangement event could be explained by the tandem duplication and random loss (TDRL) model. Our study emphasized that the combination of molecular sequence data, morphological characters and mitochondrial structural information could improve the accuracy of species delimitation.