A comparison of DNA-based methods for delimiting species in a Cretan land snail radiation reveals shortcomings of exclusively molecular taxonomy

Differentiation in the Trochulus hispidus complex and related taxa (Pulmonata: Hygromiidae): morphology, ecology and their relation to phylogeography

In this study we investigated the morphology and ecology of representatives of the taxonomically ambiguous genus Trochulus. The main focus was on the T. hispidus complex, which comprises several genetically highly divergent mitochondrial clades, as determined in a parallel molecular genetic study. We analysed shell morphology and anatomical traits and asked whether the clades are differentiated in these characters. In addition, the related species T. oreinos and T. striolatus were investigated and compared with the T. hispidus complex. Finally, we compared the ecological requirements of the taxa. Among the genetic clades of the T. hispidus complex there was no clear morphological differentiation and geographic populations could not be distinguished based on their morphology. The investigated characters of the genital anatomy did not allow discrimination of any of the T. hispidus clades and were not even diagnostic for the group as a whole. The morphotype of T. sericeus is present in all clades and thus cannot be assigned to a genetic group or any specific population. Thus, our morphological data do not provide evidence that any of the mitochondrial T. hispidus clades represent separate species. Concerning interspecific delimitation, the T. hispidus complex was clearly differentiated from T. striolatus and T. oreinos by shell morphological and anatomical characters, e.g. sculpture of shell surface and details of the penis. Finally, the habitat of T. oreinos is different from those of the other two species. In contrast to the lack of correspondence between genetic and morphological differentiation within the T. hispidus complex, related species display intraspecific morphological differentiation corresponding with mitochondrial clades: within T. striolatus there was a slight morphological differentiation between the subspecies T. s. striolatus, T. s. juvavensis and T. s. danubialis. The two subspecies of T. oreinos could be discriminated by a small but consistent difference in the cross-section of the penis. The unequal levels of intraspecific differentiation are caused by different evolutionary histories as a consequence of disparities in ecological demands, dispersal ability and use of glacial refugia: both the T. hispidus complex and T. striolatus are fast-spreading, euryoecious organisms which are able to (re-)colonize habitats and survive under different climate conditions. While the T. hispidus complex probably survived the Pleistocene in several glacial refugia, for T. striolatus one glacial refugium is suggested. Trochulus oreinos differs from the other taxa, as it is a slow disperser with a narrow ecological niche. We suggest that its subspecies spent at least the last glaciation in or close to the presently inhabited areas.

DNA barcode-based delineation of putative species: efficient start for taxonomic workflows

The analysis of DNA barcode sequences with varying techniques for cluster recognition provides an efficient approach for recognizing putative species (operational taxonomic units, OTUs). This approach accelerates and improves taxonomic workflows by exposing cryptic species and decreasing the risk of synonymy. This study tested the congruence of OTUs resulting from the application of three analytical methods (ABGD, BIN, GMYC) to sequence data for Australian hypertrophine moths. OTUs supported by all three approaches were viewed as robust, but 20% of the OTUs were only recognized by one or two of the methods. These OTUs were examined for three criteria to clarify their status. Monophyly and diagnostic nucleotides were both uninformative, but information on ranges was useful as sympatric sister OTUs were viewed as distinct, while allopatric OTUs were merged. This approach revealed 124 OTUs of Hypertrophinae, a more than twofold increase from the currently recognized 51 species. Because this analytical protocol is both fast and repeatable, it provides a valuable tool for establishing a basic understanding of species boundaries that can be validated with subsequent studies.

Paraphyly and budding speciation in the hairy snail (Pulmonata, Hygromiidae)

Delimitation of species is often complicated by discordance of morphological and genetic data. This may be caused by the existence of cryptic or polymorphic species. The latter case is particularly true for certain snail species showing an exceptionally high intraspecific genetic diversity. The present investigation deals with the Trochulus hispidus complex, which has a complicated taxonomy. Our analyses of the COI sequence revealed that individuals showing a T. hispidus phenotype are distributed in nine highly differentiated mitochondrial clades (showing p-distances up to 19%). The results of a parallel morphometric investigation did not reveal any differentiation between these clades, although the overall variability is quite high. The phylogenetic analyses based on 12S, 16S and COI sequences show that the T. hispidus complex is paraphyletic with respect to several other morphologically well-defined Trochulus species (T. clandestinus, T. villosus, T. villosulus and T. striolatus) which form well-supported monophyletic groups. The nc marker sequence (5.8S-ITS2-28S) shows only a clear separation of T. o. oreinos and T. o. scheerpeltzi, and a weakly supported separation of T. clandestinus, whereas all other species and the clades of the T. hispidus complex appear within one homogeneous group. The paraphyly of the T. hispidus complex reflects its complicated history, which was probably driven by geographic isolation in different glacial refugia and budding speciation. At our present state of knowledge, it cannot be excluded that several cryptic species are embedded within the T. hispidus complex. However, the lack of morphological differentiation of the T. hispidus mitochondrial clades does not provide any hints in this direction. Thus, we currently do not recommend any taxonomic changes. The results of the current investigation exemplify the limitations of barcoding attempts in highly diverse species such as T. hispidus.

Revisiting species delimitation within the genus Oxystele using DNA barcoding approach

The genus Oxystele, a member of the highly diverse marine gastropod superfamily Trochoidea, is endemic to southern Africa. Members of the genus include some of the most abundant molluscs on southern African shores and are important components of littoral biodiversity in rocky intertidal habitats. Species delimitation within the genus is still controversial, especially regarding the complex O. impervia / O. variegata. Here, we assessed species boundaries within the genus using DNA barcoding and phylogenetic tree reconstruction. We analysed 56 specimens using the mitochondrial gene COI. Our analysis delimits five molecular operational taxonomic units (MOTUs), and distinguishes O. impervia from O. variegata. However, we reveal important discrepancies between MOTUs and morphology-based species identification and discuss alternative hypotheses that can account for this. Finally, we indicate the need for future study that includes additional genes, and the combination of both morphology and genetic techniques (e.g. AFLP or microsatellites) to get deeper insight into species delimitation within the genus.

A general species delimitation method with applications to phylogenetic placements

MOTIVATION

Sequence-based methods to delimit species are central to DNA taxonomy, microbial community surveys and DNA metabarcoding studies. Current approaches either rely on simple sequence similarity thresholds (OTU-picking) or on complex and compute-intensive evolutionary models. The OTU-picking methods scale well on large datasets, but the results are highly sensitive to the similarity threshold. Coalescent-based species delimitation approaches often rely on Bayesian statistics and Markov Chain Monte Carlo sampling, and can therefore only be applied to small datasets.

RESULTS

We introduce the Poisson tree processes (PTP) model to infer putative species boundaries on a given phylogenetic input tree. We also integrate PTP with our evolutionary placement algorithm (EPA-PTP) to count the number of species in phylogenetic placements. We compare our approaches with popular OTU-picking methods and the General Mixed Yule Coalescent (GMYC) model. For de novo species delimitation, the stand-alone PTP model generally outperforms GYMC as well as OTU-picking methods when evolutionary distances between species are small. PTP neither requires an ultrametric input tree nor a sequence similarity threshold as input. In the open reference species delimitation approach, EPA-PTP yields more accurate results than de novo species delimitation methods. Finally, EPA-PTP scales on large datasets because it relies on the parallel implementations of the EPA and RAxML, thereby allowing to delimit species in high-throughput sequencing data.

AVAILABILITY AND IMPLEMENTATION

The code is freely available at www.exelixis-lab.org/software.html. .

How to describe a cryptic species? Practical challenges of molecular taxonomy

BACKGROUND

Molecular methods of species delineation are rapidly developing and widely considered as fast and efficient means to discover species and face the 'taxonomic impediment' in times of biodiversity crisis. So far, however, this form of DNA taxonomy frequently remains incomplete, lacking the final step of formal species description, thus enhancing rather than reducing impediments in taxonomy. DNA sequence information contributes valuable diagnostic characters and -at least for cryptic species - could even serve as the backbone of a taxonomic description. To this end solutions for a number of practical problems must be found, including a way in which molecular data can be presented to fulfill the formal requirements every description must meet. Multi-gene barcoding and a combined molecular species delineation approach recently revealed a radiation of at least 12 more or less cryptic species in the marine meiofaunal slug genus Pontohedyle (Acochlidia, Heterobranchia). All identified candidate species are well delimited by a consensus across different methods based on mitochondrial and nuclear markers.

RESULTS

The detailed microanatomical redescription of Pontohedyle verrucosa provided in the present paper does not reveal reliable characters for diagnosing even the two major clades identified within the genus on molecular data. We thus characterize three previously valid Pontohedyle species based on four genetic markers (mitochondrial cytochrome c oxidase subunit I, 16S rRNA, nuclear 28S and 18S rRNA) and formally describe nine cryptic new species (P. kepii sp. nov., P. joni sp. nov., P. neridae sp. nov., P. liliae sp. nov., P. wiggi sp. nov., P. wenzli sp. nov., P. peteryalli sp. nov., P. martynovi sp. nov., P. yurihookeri sp. nov.) applying molecular taxonomy, based on diagnostic nucleotides in DNA sequences of the four markers. Due to the minute size of the animals, entire specimens were used for extraction, consequently the holotype is a voucher of extracted DNA ('DNA-type'). We used the Character Attribute Organization System (CAOS) to determine diagnostic nucleotides, explore the dependence on input data and data processing, and aim for maximum traceability in our diagnoses for future research. Challenges, pitfalls and necessary considerations for applied DNA taxonomy are critically evaluated.

CONCLUSIONS

To describe cryptic species traditional lines of evidence in taxonomy need to be modified. DNA sequence information, for example, could even serve as the backbone of a taxonomic description. The present contribution demonstrates that few adaptations are needed to integrate into traditional taxonomy novel diagnoses based on molecular data. The taxonomic community is encouraged to join the discussion and develop a quality standard for molecular taxonomy, ideally in the form of an automated final step in molecular species delineation procedures.

Morphological and genetic evidence for multiple evolutionary distinct lineages in the endangered and commercially exploited red lined torpedo barbs endemic to the Western Ghats of India

Red lined torpedo barbs (RLTBS) (Cyprinidae: Puntius) endemic to the Western Ghats Hotspot of India, are popular and highly priced freshwater aquarium fishes. Two decades of indiscriminate exploitation for the pet trade, restricted range, fragmented populations and continuing decline in quality of habitats has resulted in their 'Endangered' listing. Here, we tested whether the isolated RLTB populations demonstrated considerable variation qualifying to be considered as distinct conservation targets. Multivariate morphometric analysis using 24 size-adjusted characters delineated all allopatric populations. Similarly, the species-tree highlighted a phylogeny with 12 distinct RLTB lineages corresponding to each of the different riverine populations. However, coalescence-based methods using mitochondrial DNA markers identified only eight evolutionarily distinct lineages. Divergence time analysis points to recent separation of the populations, owing to the geographical isolation, more than 5 million years ago, after the lineages were split into two ancestral stocks in the Paleocene, on north and south of a major geographical gap in the Western Ghats. Our results revealing the existence of eight evolutionarily distinct RLTB lineages calls for the re-determination of conservation targets for these cryptic and endangered taxa.

New metrics for comparison of taxonomies reveal striking discrepancies among species delimitation methods in Madascincus lizards

Delimiting and describing species is fundamental to numerous biological disciplines such as evolution, macroecology, and conservation. Delimiting species as independent evolutionary lineages may and often does yield different outcomes depending on the species criteria applied, but methods should be chosen that minimize the inference of objectively erroneous species limits. Several protocols exploit single-gene or multi-gene coalescence statistics, assignment tests or other rationales related to nuclear DNA (nDNA) allele sharing to automatically delimit species. We apply seven different species delimitation protocols to a taxonomically confusing group of Malagasy lizards (Madascincus), and compare the resulting taxonomies with two newly developed metrics: the Taxonomic index of congruence C tax which quantifies the congruence between two taxonomies, and the Relative taxonomic resolving power index R tax which quantifies the potential of an approach to capture a high number of species boundaries. The protocols differed in the total number of species proposed, between 9 and 34, and were also highly incongruent in placing species boundaries. The Generalized Mixed Yule-Coalescent approach captured the highest number of potential species boundaries but many of these were clearly contradicted by extensive nDNA admixture between sympatric mitochondrial DNA (mtDNA) haplotype lineages. Delimiting species as phenotypically diagnosable mtDNA clades failed to detect two cryptic species that are unambiguous due to a lack of nDNA gene flow despite sympatry. We also consider the high number of species boundaries and their placement by multi-gene Bayesian species delimitation as poorly reliable whereas the Bayesian assignment test approach provided a species delimitation highly congruent with integrative taxonomic practice. The present study illustrates the trade-off in taxonomy between reliability (favored by conservative approaches) and resolving power (favored by inflationist approaches). Quantifying excessive splitting is more difficult than quantifying excessive lumping, suggesting a priority for conservative taxonomies in which errors are more liable to be detected and corrected by subsequent studies.

Continuing fragmentation of a widespread species by geographical barriers as initial step in a land snail radiation on crete

The phylogeographic structure of the land snail Xerocrassa mesostena on Crete inferred from AFLP markers and mitochondrial cox1 sequences can be explained by three mechanisms: gene flow restriction, population expansion and leptokurtic dispersal. Gene flow restriction by geographic barriers caused subdivision of the gene pool into distinct clusters. Population expansion was probably facilitated by deforestation of Crete in the postglacial. Newly available areas were colonized by leptokurtic dispersal, i.e. slow active expansion resulting in isolation by distance within the clusters and occasional long distance dispersal events that resulted in departures from the isolation by distance model. Less than one percent of the AFLP markers show correlations with environmental variables. Random phylogeographic breaks in the distribution of the mitochondrial haplotype groups indicate that single locus markers, especially mitochondrial DNA, might result in a misleading picture of the phylogeographic structure of a species. Restriction of gene flow between metapopulations caused by geographical barriers can interact with sexual selection resulting in the differentiation of these metapopulations into separate species without noticeable ecological differentiation. Evidence for gene flow between parapatrically distributed evolutionary units representing different stages of the speciation process suggests that the ongoing process of fragmentation of the X. mesostena complex might be an example for parapatric speciation. The lack of ecological differentiation between these units confirms theoretical predictions that divergent selection for local adaptation is not required for rapid speciation.

Exploring species level taxonomy and species delimitation methods in the facultatively self-fertilizing land snail genus Rumina (gastropoda: pulmonata)

Delimiting species in facultatively selfing taxa is a challenging problem of which the terrestrial pulmonate snail genus Rumina is a good example. These snails have a mixed breeding system and show a high degree of shell and color variation. Three nominal species (R. decollata, R. saharica and R. paivae) and two color morphs within R. decollata (dark and light) are currently recognized. The present study aims at evaluating to what extent these entities reflect evolutionary diverging taxonomic units, rather than fixed polymorphisms due to sustained selfing. Therefore, a phylogenetic analysis of nuclear (ITS1, ITS2) and mitochondrial DNA (COI, CytB, 12S rDNA, 16S rDNA) sequences was performed. Putative species in Rumina, inferred from the mitochondrial DNA phylogeny, were compared with those proposed on the basis of the COI gene by (1) DNA barcoding gap analysis, (2) Automatic Barcode Gap Discovery, (3) the species delimitation plug-in of the Geneious software, (4) the Genealogical Sorting Index, and (5) the General Mixed Yule Coalescent model. It is shown that these methods produce a variety of different species hypotheses and as such one may wonder to what extent species delimitation methods are really useful. With respect to Rumina, the data suggest at least seven species, one corresponding to R. saharica and six that are currently grouped under the name R. decollata. The species-level status of R. paivae is rejected.

Barcoding against a paradox? Combined molecular species delineations reveal multiple cryptic lineages in elusive meiofaunal sea slugs

BACKGROUND

Many marine meiofaunal species are reported to have wide distributions, which creates a paradox considering their hypothesized low dispersal abilities. Correlated with this paradox is an especially high taxonomic deficit for meiofauna, partly related to a lower taxonomic effort and partly to a high number of putative cryptic species. Molecular-based species delineation and barcoding approaches have been advocated for meiofaunal biodiversity assessments to speed up description processes and uncover cryptic lineages. However, these approaches show sensitivity to sampling coverage (taxonomic and geographic) and the success rate has never been explored on mesopsammic Mollusca.

RESULTS

We collected the meiofaunal sea-slug Pontohedyle (Acochlidia, Heterobranchia) from 28 localities worldwide. With a traditional morphological approach, all specimens fall into two morphospecies. However, with a multi-marker genetic approach, we reveal multiple lineages that are reciprocally monophyletic on single and concatenated gene trees in phylogenetic analyses. These lineages are largely concordant with geographical and oceanographic parameters, leading to our primary species hypothesis (PSH). In parallel, we apply four independent methods of molecular based species delineation: General Mixed Yule Coalescent model (GMYC), statistical parsimony, Bayesian Species Delineation (BPP) and Automatic Barcode Gap Discovery (ABGD). The secondary species hypothesis (SSH) is gained by relying only on uncontradicted results of the different approaches ('minimum consensus approach'), resulting in the discovery of a radiation of (at least) 12 mainly cryptic species, 9 of them new to science, some sympatric and some allopatric with respect to ocean boundaries. However, the meiofaunal paradox still persists in some Pontohedyle species identified here with wide coastal and trans-archipelago distributions.

CONCLUSIONS

Our study confirms extensive, morphologically cryptic diversity among meiofauna and accentuates the taxonomic deficit that characterizes meiofauna research. We observe for Pontohedyle slugs a high degree of morphological simplicity and uniformity, which we expect might be a general rule for meiofauna. To tackle cryptic diversity in little explored and hard-to-sample invertebrate taxa, at present, a combined approach seems most promising, such as multi-marker-barcoding (i.e., molecular systematics using mitochondrial and nuclear markers and the criterion of reciprocal monophyly) combined with a minimum consensus approach across independent methods of molecular species delineation to define candidate species.