Disintegrating over space and time: Paraphyly and species delimitation in the Wehrle's Salamander complex

Accurate species delimitation is critical for biodiversity studies. However, species complexes characterized by introgression, high levels of population structure and subtle phenotypic differentiation can be challenging to delimit. Here, we report on a molecular systematic investigation of the woodland salamanders Plethodon wehrlei and Plethodon punctatus, which traditionally have been placed in the Plethodon wehrlei species group. To quantify patterns of genetic variation, we collected genetic samples from throughout the range of both species, including 22 individuals from nine populations of P. punctatus, and 60 individuals from 26 populations of P. wehrlei. From these samples, we sequenced three mtDNA loci (5596 base pairs) and five nuclear loci (3377 base pairs). We inferred time‐calibrated gene trees and species trees using BEAST 2.4.6, and we delimited putative species using a Bayesian implementation of the general mixed Yule‐coalescent model (bGMYC) and STRUCTURE. Finally, we validated putative species using the multispecies coalescent as implemented in Bayesian Phylogenetics and Phylogeography (BPP). We found substantial phylogeographic diversity in P. wehrlei, including multiple geographically cohesive clades and an inferred mitochondrial common ancestor at 11.5 myr (95% HPD: 9.6–13.6 myr) that separated populations formerly assigned to P. dixi from all other populations. We also found that P. punctatus is deeply nested within P. wehrlei, rendering the latter paraphyletic. After discussing the challenges faced by modern species delimitation methods, we recommend retaining P. punctatus because it is ecologically and phenotypically distinct. We further recommend that P. dixi be recognized as a valid species.     

Phylogeny and species delimitation of North European Lumbricillus (Clitellata,Enchytraeidae)

The enchytraeid genus Lumbricillus comprises about 80 described species of clitellate worms, which are up to a few centimetres long, and they mostly inhabit the littoral zone of non‐tropical marine and brackish waters world‐wide. The phylogeny of this genus is poorly studied, but previous work has suggested that Lumbricillus is a non‐monophyletic group. In this study, species boundaries and the phylogeny of this genus is re‐assessed using more than 300 DNA‐barcoded specimens (using COI mtDNA), part of which was also sequenced for two additional mitochondrial and four nuclear molecular markers. Statistical and coalescent based applications were used for the delimitation of a total of 24 species, of which 20 were identified as belonging to 17 described morphospecies; one morphospecies was found to be a complex of four delimited species, and another four delimited species could not be matched with any described species. Furthermore, gene trees, concatenation and multispecies coalescent based species trees were estimated using Bayesian inference. The estimated phylogenies confirm a non‐monophyletic Lumbricillus as L. semifuscus is clearly excluded from the genus. Furthermore, the placement of a monophyletic clade consisting of L. arenarius, L. dubius, and an unidentified species varies between analyses; they are either found as the sister‐group to the genus Grania or as sister‐group to the remaining Lumbricillus, where the latter relationship is supported by the multispecies coalescent, which we consider as the most reliable method.     

A polyphasic approach to the delimitation of diatom species: a case study for the genus Pinnularia (Bacillariophyta)

Diatoms are one of the most abundant and arguably the most species‐rich group of protists. Diatom species delimitation has often been based exclusively on the recognition of morphological discontinuities without investigation of other lines of evidence. Even though DNA sequences and reproductive experiments have revealed several examples of (pseudo)cryptic diversity, our understanding of diatom species boundaries and diversity remains limited. The cosmopolitan pennate raphid diatom genus Pinnularia represents one of the most taxon‐rich diatom genera. In this study, we focused on the delimitation of species in one of the major clades of the genus, the Pinnularia subgibba group, based on 105 strains from a worldwide origin. We compared genetic distances between the sequences of seven molecular markers and selected the most variable pair, the mitochondrial cox1 and nuclear encoded LSU rDNA, to formulate a primary hypothesis on the species limits using three single‐locus automated species delimitation methods. We compared the DNA‐based primary hypotheses with morphology and with other available lines of evidence. The results indicate that our data set comprised 15 species of the P. subgibba group. The vast majority of these taxa have an uncertain taxonomic identity, suggesting that several may be unknown to science and/or members of (pseudo)cryptic species complexes within the P. subgibba group.     

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.     

Testing hypotheses of mitochondrial gene‐tree paraphyly: unravelling mitochondrial capture of the Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) by the Taiwan Scimitar Babbler (Pomatorhinus musicus)

Species‐level paraphyly inferred from mitochondrial gene trees is a prevalent phenomenon in taxonomy and systematics, but there are several potential causes that are not easily explained by currently used methods. This study aimed to test the underlying causes behind the observed paraphyly of Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) via statistical analyses of four mitochondrial (mtDNA) and nine nuclear (nuDNA) genes. Mitochondrial gene trees show paraphyly of P. ruficollis with respect to the Taiwan Scimitar Babbler (Pomatorhinus musicus), but nuclear genealogies support a sister‐group relationship. Predictive coalescent simulations imply several hypothetical explanations, the most likely being mitochondrial capture of P. ruficollis by P. musicus for the observed cyto‐nuclear incongruence. Further approximate Bayesian computation suggests a unidirectional introgression model with substantial level of gene flow from P. ruficollis to P. musicus during their initial divergence during the Late Pleistocene. This specific observation frames several potential causes for incongruent outcomes of mitochondrial and nuclear introgression in general, and on the whole, our results underscore the strength of multiple independent loci for species delimitation and importance of testing hypotheses that explain disparate causes of mitochondrial gene‐tree paraphyly.     

Species delimitation with gene flow: A methodological comparison and population genomics approach to elucidate cryptic species boundaries in Malaysian Torrent Frogs

Accurately delimiting species boundaries is a nontrivial undertaking that can have significant effects on downstream inferences. We compared the efficacy of commonly used species delimitation methods (SDMs) and a population genomics approach based on genomewide single‐nucleotide polymorphisms (SNPs) to assess lineage separation in the Malaysian Torrent Frog Complex currently recognized as a single species (Amolops larutensis). First, we used morphological, mitochondrial DNA and genomewide SNPs to identify putative species boundaries by implementing noncoalescent and coalescent‐based SDMs (mPTP, iBPP, BFD*). We then tested the validity of putative boundaries by estimating spatiotemporal gene flow (fastsimcoal2 , ABBA‐BABA) to assess the extent of genetic isolation among putative species. Our results show that the A. larutensis complex runs the gamut of the speciation continuum from highly divergent, genetically isolated lineages (mean F_st = 0.9) to differentiating populations involving recent gene flow (mean F_st = 0.05; N_m > 5). As expected, SDMs were effective at delimiting divergent lineages in the absence of gene flow but overestimated species in the presence of marked population structure and gene flow. However, using a population genomics approach and the concept of species as separately evolving metapopulation lineages as the only necessary property of a species, we were able to objectively elucidate cryptic species boundaries in the presence of past and present gene flow. This study does not discount the utility of SDMs but highlights the danger of violating model assumptions and the importance of carefully considering methods that appropriately fit the diversification history of a particular system.     

Bayesian Poisson tree processes and multispecies coalescent models shed new light on the diversification of Nawab butterflies in the Solomon Islands (Nymphalidae,Charaxinae, Polyura)

Butterflies of the genus Polyura form a widespread tropical group distributed from Pakistan to Fiji. The rare endemic Polyura epigenes Godman & Salvin, 1888 from the Solomon Islands archipelago represents a case of marked island polymorphism. We sequenced museum specimens of this species across its geographic range to study the phylogeography and genetic differentiation of populations in the archipelago. We used the Bayesian Poisson tree processes and multispecies coalescent models, to study species boundaries. We also estimated divergence times to investigate the biogeographic history of populations. Our molecular species delimitation and nuclear DNA network analyses unambiguously indicate that Malaita populations form an independent metapopulation lineage, as defined in the generalized lineage concept. This lineage, previously ranked as a subspecies, is raised to species rank under the name Polyura bicolor Turlin & Sato, 1995  stat. nov. Divergence time estimates suggest that this lineage split from its sister taxon in the late Pleistocene. At this time, the bathymetric isolation of Malaita from the rest of the archipelago probably prevented gene flow during periods of lower sea level, thereby fostering allopatric speciation. The combination of molecular species delimitation methods, morphological comparisons, and divergence time estimation is useful to study lineage diversification across intricate geographic regions.     

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.     

Cryptic diversity in the Mexican highlands: Thousands of UCE loci help illuminate phylogenetic relationships,species limits and divergence times of montane rattlesnakes (Viperidae: Crotalus)

With the continued adoption of genome‐scale data in evolutionary biology comes the challenge of adequately harnessing the information to make accurate phylogenetic inferences. Coalescent‐based methods of species tree inference have become common, and concatenation has been shown in simulation to perform well, particularly when levels of incomplete lineage sorting are low. However, simulation conditions are often overly simplistic, leaving empiricists with uncertainty regarding analytical tools. We use a large ultraconserved element data set (>3,000 loci) from rattlesnakes of the Crotalus triseriatus group to delimit lineages and estimate species trees using concatenation and several coalescent‐based methods. Unpartitioned and partitioned maximum likelihood and Bayesian analysis of the concatenated matrix yield a topology identical to coalescent analysis of a subset of the data in bpp . ASTRAL analysis on a subset of the more variable loci also results in a tree consistent with concatenation and bpp , whereas the SVDquartets phylogeny differs at additional nodes. The size of the concatenated matrix has a strong effect on species tree inference using SVDquartets , warranting additional investigation on optimal data characteristics for this method. Species delimitation analyses suggest up to 16 unique lineages may be present within the C. triseriatus group, with divergences occurring during the Neogene and Quaternary. Network analyses suggest hybridization within the group is relatively rare. Altogether, our results reaffirm the Mexican highlands as a biodiversity hotspot and suggest that coalescent‐based species tree inference on data subsets can provide a strongly supported species tree consistent with concatenation of all loci with a large amount of missing data.     

Coalescent histories for discordant gene trees and species trees

Given a gene tree and a species tree, a coalescent history is a list of the branches of the species tree on which coalescences in the gene tree take place. Each pair consisting of a gene tree topology and a species tree topology has some number of possible coalescent histories. Here we show that, for each n≥7, there exist a species tree topology S and a gene tree topology G≠S, both with n leaves, for which the number of coalescent histories exceeds the corresponding number of coalescent histories when the species tree topology is S and the gene tree topology is also S. This result has the interpretation that the gene tree topology G discordant with the species tree topology S can be produced by the evolutionary process in more ways than can the gene tree topology that matches the species tree topology, providing further insight into the surprising combinatorial properties of gene trees that arise from their joint consideration with species trees.     

DNA-based species delimitation in algae

Given the problems of species delimitation in algae using morphology or sexual compatibility, molecular data are becoming the standard for delimiting species and testing their traditional boundaries. The idea that species are separately evolving metapopulation lineages, along with theoretical progress in phylogenetic and population genetic analyses, has led to the development of new methods of species delimitation. We review these recent developments in DNA-based species delimitation methods, and discuss how they have changed and continue to change our understanding of algal species boundaries. Although single-locus approaches have proven effective for a first rapid and large-scale assessment of species diversity, species delimitation based on single gene trees falls short due to gene tree–species tree incongruence, caused by confounding processes like incomplete lineage sorting, trans-species polymorphism, hybridization and introgression. Data from unlinked loci and multi-species coalescent methods, which combine principles from phylogenetics and population genetics, may now be able to account for these complicating factors. Several of these methods also provide statistical support regarding species boundaries, which is important because speciation is a process and therefore uncertainty about precise species boundaries is inevitable in recently diverged lineages.     

Species limits in butterflies (Lepidoptera: Nymphalidae): reconciling classical taxonomy with the multispecies coalescent

Species delimitation is at the core of biological sciences. During the last decade, molecular‐based approaches have advanced the field by providing additional sources of evidence to classical, morphology‐based taxonomy. However, taxonomy has not yet fully embraced molecular species delimitation beyond threshold‐based, single‐gene approaches, and taxonomic knowledge is not commonly integrated into multilocus species delimitation models. Here we aim to bridge empirical data (taxonomic and genetic) with recently developed coalescent‐based species delimitation approaches. We use the multispecies coalescent model as implemented in two Bayesian methods (dissect/stacey and bp&p ) to infer species hypotheses. In both cases, we account for phylogenetic uncertainty (by not using any guide tree) and taxonomic uncertainty (by measuring the impact of using a priori taxonomic assignments to specimens). We focus on an entire Neotropical tribe of butterflies, the Haeterini (Nymphalidae: Satyrinae). We contrast divergent taxonomic opinion – splitting, lumping and misclassifying species – in the light of different phenotypic classifications proposed to date. Our results provide a solid background for the recognition of 22 species. The synergistic approach presented here overcomes limitations in both traditional taxonomy (e.g. by recognizing cryptic species) and molecular‐based methods (e.g. by recognizing structured populations, and not raising them to species). Our framework provides a step forward towards standardization and increasing reproducibility of species delimitations.     

Species delimitation and phylogeny in the genus Nasutitermes (Termitidae: Nasutitermitinae) in French Guiana

Species delimitation and identification can be arduous for taxa whose morphologic characters are easily confused, which can hamper global biodiversity assessments and pest species management. Exploratory methods of species delimitation that use DNA sequence as their primary information source to establish group membership and estimate putative species boundaries are useful approaches, complementary to traditional taxonomy. Termites of the genus Nasutitermes make interesting models for the application of such methods. They are dominant in Neotropical primary forests but also represent major agricultural and structural pests. Despite the prevalence, pivotal ecological role and economical impact of this group, the taxonomy of Nasutitermes species mainly depends on unreliable characters of soldier external morphology. Here, we generated robust species hypotheses for 79 Nasutitermes colonies sampled throughout French Guiana without any a priori knowledge of species affiliation. Sequence analysis of the mitochondrial cytochrome oxidase II gene was coupled with exploratory species‐delimitation tools, using the automatic barcode gap discovery method (ABGD) and a generalized mixed Yule‐coalescent model (GMYC) to propose primary species hypotheses (PSHs). PSHs were revaluated using phylogenetic analyses of two more loci (mitochondrial 16S rDNA and nuclear internal transcribed spacer 2) leading to 16 retained secondary species hypotheses (RSSH). Seven RSSHs, represented by 44/79 of the sampled colonies, were morphologically affiliated to species recognized as pests in the Neotropics, where they represent a real invasive pest potential in the context of growing ecosystem anthropization. Multigenic phylogenies based on combined alignments (1426–1784 bp) were also reconstructed to identify ancestral ecological niches and major‐pest lineages, revealing that Guyanese pest species do not form monophyletic groups.     

Among‐tree variability and feedback effects result in different growth responses to climate change at the upper treeline in the Swiss Alps

Upper treeline ecotones are important life form boundaries and particularly sensitive to a warming climate. Changes in growth conditions at these ecotones have wide‐ranging implications for the provision of ecosystem services in densely populated mountain regions like the European Alps. We quantify climate effects on short‐ and long‐term tree growth responses, focusing on among‐tree variability and potential feedback effects. Although among‐tree variability is thought to be substantial, it has not been considered systematically yet in studies on growth–climate relationships. We compiled tree‐ring data including almost 600 trees of major treeline species (Larix decidua, Picea abies, Pinus cembra, and Pinus mugo) from three climate regions of the Swiss Alps. We further acquired tree size distribution data using unmanned aerial vehicles. To account for among‐tree variability, we employed information‐theoretic model selections based on linear mixed‐effects models (LMMs) with flexible choice of monthly temperature effects on growth. We isolated long‐term trends in ring‐width indices (RWI) in interaction with elevation. The LMMs revealed substantial amounts of previously unquantified among‐tree variability, indicating different strategies of single trees regarding when and to what extent to invest assimilates into growth. Furthermore, the LMMs indicated strongly positive temperature effects on growth during short summer periods across all species, and significant contributions of fall (L. decidua) and current year's spring (L. decidua, P. abies). In the longer term, all species showed consistently positive RWI trends at highest elevations, but different patterns with decreasing elevation. L. decidua exhibited even negative RWI trends compared to the highest treeline sites, whereas P. abies, P. cembra, and P. mugo showed steeper or flatter trends with decreasing elevation. This does not only reflect effects of ameliorated climate conditions on tree growth over time, but also reveals first signs of long‐suspected negative and positive feedback of climate change on stand dynamics at treeline.     

Estimating the temporal and spatial extent of gene flow among sympatric lizard populations (genus Sceloporus) in the southern Mexican highlands

Interspecific gene flow is pervasive throughout the tree of life. Although detecting gene flow between populations has been facilitated by new analytical approaches, determining the timing and geography of hybridization has remained difficult, particularly for historical gene flow. A geographically explicit phylogenetic approach is needed to determine the overlap of ancestral populations. In this study, we performed population genetic analyses, species delimitation, simulations and a recently developed approach of species tree diffusion to infer the phylogeographic history, timing and geographic extent of gene flow in lizards of the Sceloporus spinosus group. The two species in this group, S. spinosus and S. horridus, are distributed in eastern and western portions of Mexico, respectively, but populations of these species are sympatric in the southern Mexican highlands. We generated data consisting of three mitochondrial genes and eight nuclear loci for 148 and 68 individuals, respectively. We delimited six lineages in this group, but found strong evidence of mito‐nuclear discordance in sympatric populations of S. spinosus and S. horridus owing to mitochondrial introgression. We used coalescent simulations to differentiate ancestral gene flow from secondary contact, but found mixed support for these two models. Bayesian phylogeography indicated more than 60% range overlap between ancestral S. spinosus and S. horridus populations since the time of their divergence. Isolation–migration analyses, however, revealed near‐zero levels of gene flow between these ancestral populations. Interpreting results from both simulations and empirical data indicate that despite a long history of sympatry among these two species, gene flow in this group has only recently occurred.     

Resolving species boundaries in the critically imperiled freshwater mussel species,Fusconaia mitchelli (Bivalvia: Unionidae)

Species are a fundamental unit of biology, and defining accurate species boundaries is integral to effective conservation and management of imperiled taxa. Freshwater mussels (Bivalvia: Unionidae) are among the most imperiled groups of organisms in North America, yet species boundaries remain uncertain for many taxa. The False Spike, Fusconaia mitchelli (Simpson in Dall, 1895), is a freshwater mussel considered to be endemic to central Texas (Brazos, Colorado, and Guadalupe drainages). Recent research revealed significant intraspecific genetic variation between geographically separated populations of F. mitchelli, which could be indicative of speciation; however, small sample sizes for several of the populations precluded formal taxonomic revision. Here, we increase taxon sampling and use multilocus DNA sequence data and traditional morphometrics to re‐evaluate species boundaries in F. mitchelli. We sequenced three loci: the protein‐coding mitochondrial DNA genes cytochrome c oxidase subunit 1 and NADH dehydrogenase 1, and the nuclear internal transcribed spacer 1. Phylogenetic analyses depicted deep genetic divergence between F. mitchelli in the Guadalupe and those in the Brazos and Colorado drainages, which was further supported by available biogeographic information. Morphometric analyses and coalescent‐based species delimitation models integrating both DNA sequence and morphological data provided strong support for the divergence observed between the two geographically isolated clades of F. mitchelli. Based on these results, we revise taxonomy accordingly by elevating the junior synonym Fusconaia iheringi (Wright, 1898) to represent the Brazos and Colorado populations and restrict the distribution of F. mitchelli to the Guadalupe River drainage. Our findings may impact pending management decisions to protect F. mitchelli under the U.S. Endangered Species Act.     

Molecular systematics,species limits,and diversification of the genus Dendrocolaptes (Aves: Furnariidae): Insights on biotic exchanges between dry and humid forest types in the Neotropics

The true diversity and interspecific limits in the Neotropical endemic avian genus Dendrocolaptes (Furnariidae) remain a highly controversial subject, with previous genus‐wide assessments, based mostly on morphological characters, producing poorly resolved phylogenies. The lack of well‐resolved, robust, and taxonomically densely sampled phylogenies for Dendrocolaptes prevents reliable inferences on the genus’ actual species diversity and evolutionary history. Here, we analyzed 2,741 base pairs of mitochondrial and nuclear genes from 43 specimens belonging to all species and the majority of subspecies described for Dendrocolaptes to evaluate species limits and reconstruct its diversification through time. Our phylogenies recovered a monophyletic Dendrocolaptes, with two main highly supported internal clades corresponding to the D. certhia and D. picumnus species complexes. Also, our analyses supported the monophyly of most Dendrocolaptes species recognized today, except D. picumnus, which was consistently recovered as paraphyletic with respect to D. hoffmannsi. A coalescent‐based test supported a total of 15 different lineages in Dendrocolaptes and indicated that the number of currently accepted species within the genus may be greatly underestimated. Particularly relevant, when combined with previous analyses based on plumage characters, comparative high levels of genetic differentiation and coalescent analyses support the recognition of D. picumnus transfasciatus as a full species that is already under threat. Ancestral area reconstructions suggest that diversification in Dendrocolaptes was centered in lowland Amazonia, with several independent dispersal events leading to differentiation into different adjacent dry and high elevation forest types throughout the Neotropics, mainly during the Middle and Late Pleistocene.     

Assessing the diversity and distribution of Cephalothrix species (Nemertea: Palaeonemertea) in European waters by comparing different species delimitation methods

Soft‐bodied marine taxa, like ribbon worms (Nemertea), often lack clear diagnostic morphological characters impeding traditional species delimitation. Therefore, recent studies concentrated on molecular genetic methods to solve taxonomic issues. Different delimitation methods were employed to explore species boundaries and the presence of cryptic species. However, the performance of the different delimitation methods needs to be tested. A particularly promising nemertean genus in this regard is the palaeonemertean genus Cephalothrix that is commonly found in European waters. In order to gain information on the number and distribution of European cephalotrichids and to test different tree‐based and non‐tree‐based delimitation methods, we analyzed a dataset comprising the barcoding region of the mitochondrial cytochrome c oxidase subunit I (COI) of 215 European Cephalothrix specimens, of which 78 were collected for this study. Our results show the presence of 12–13 European lineages of which several can be assigned to known European species. Analyzing a second dataset comprising 74 additional sequences from the Pacific and the Atlantic Oceans helped identify some of the unassigned European specimens. One resulting clade seems to represent a non‐native introduced Cephalothrix species, while another has never been recorded from Europe before. In our analysis, especially the tree‐based methods and the phylogenetic analysis proved to be a useful tool when delimiting species. It remains unclear whether the different identified clades result from cryptic speciation or from a high genetic variability of the COI gene.