The necessity of DNA taxonomy to reveal cryptic diversity and spatial distribution of meiofauna, with a focus on Nemertea

Circumtropical distribution and cryptic species of the meiofaunal enteropneust Meioglossus (Harrimaniidae, Hemichordata)

Hemichordata has always played a central role in evolutionary studies of Chordata due to their close phylogenetic affinity and shared morphological characteristics. Hemichordates had no meiofaunal representatives until the surprising discovery of a microscopic, paedomorphic enteropneust Meioglossus psammophilus (Harrimaniidae, Hemichordata) from the Caribbean in 2012. No additional species have been described since, questioning the broader distribution and significance of this genus. However, being less than a millimeter long and superficially resembling an early juvenile acorn worm, Meioglossus may easily be overlooked in both macrofauna and meiofauna surveys. We here present the discovery of 11 additional populations of Meioglossus from shallow subtropical and tropical coralline sands of the Caribbean Sea, Red Sea, Indian Ocean, and East China Sea. These geographically separated populations show identical morphology but differ genetically. Our phylogenetic reconstructions include four gene markers and support the monophyly of Meioglossus. Species delineation analyses revealed eight new cryptic species, which we herein describe using DNA taxonomy. This study reveals a broad circumtropical distribution, supporting the validity and ecological importance of this enigmatic meiobenthic genus. The high cryptic diversity and apparent morphological stasis of Meioglossus may exemplify a potentially common evolutionary 'dead-end' scenario, where groups with highly miniaturized and simplified body plan lose their ability to diversify morphologically.

Sampling multiple life stages significantly increases estimates of marine biodiversity

Biodiversity assessments are critical for setting conservation priorities, understanding ecosystem function and establishing a baseline to monitor change. Surveys of marine biodiversity that rely almost entirely on sampling adult organisms underestimate diversity because they tend to be limited to habitat types and individuals that can be easily surveyed. Many marine animals have planktonic larvae that can be sampled from the water column at shallow depths. This life stage often is overlooked in surveys but can be used to relatively rapidly document diversity, especially for the many species that are rare or live cryptically as adults. Using DNA barcode data from samples of nemertean worms collected in three biogeographical regions—Northeastern Pacific, the Caribbean Sea and Eastern Tropical Pacific—we found that most species were collected as either benthic adults or planktonic larvae but seldom in both stages. Randomization tests show that this deficit of operational taxonomic units collected as both adults and larvae is extremely unlikely if larvae and adults were drawn from the same pool of species. This effect persists even in well-studied faunas. These results suggest that sampling planktonic larvae offers access to a different subset of species and thus significantly increases estimates of biodiversity compared to sampling adults alone. Spanish abstract is available in the electronic supplementary material.

Kulikovia alborostrata and Kulikovia fulva comb. nov. (Nemertea: Heteronemertea) are Sister Species with Prezygotic Isolating Barriers

The heteronemertean Kulikovia alborostrata (Takakura, 1898) was originally described as Lineus alborostratus based on material from Misaki, Japan. Although this species was regarded as consisting of two color variants, purple and brown-yellow, the identity of these variants has never been examined based on topotypes. In this study, we performed a multi-locus phylogeny reconstruction, species delimitation analyses, and cross-fertilization experiments to examine the species status of Takakura's original taxon concept consisting of these color variants. Our results suggest that the purple type is identical to Lineus alborostratus Takakura, 1898 auct. (currently Kulikovia alborostrata), whereas the brown-yellow type is conspecific with Lineus fulvus Iwata, 1954, originally established from Hokkaido. These two species appear to have a sister-taxon relationship and are reproductively isolated from each other by prezygotic mechanisms involving gamete incompatibility, minimally separated with 2.8% (16S rRNA) and 14.4% (COI) uncorrected p-distances. We propose that the purple type be considered as representing the true identity of the nominal species Lineus alborostratus (currently assigned to the genus Kulikovia) to maintain the common usage of the name. Although Takakura's type material is not extant, we consider that neotypification is unnecessary in this case because no taxonomic/nomenclatural confusion persists. We also propose to transfer Lineus fulvus to yield Kulikovia fulva comb. nov.

Incomplete lineage sorting and ancient admixture, and speciation without morphological change in ghost-worm cryptic species

Morphologically similar species, that is cryptic species, may be similar or quasi-similar owing to the deceleration of morphological evolution and stasis. While the factors underlying the deceleration of morphological evolution or stasis in cryptic species remain unknown, decades of research in the field of paleontology on punctuated equilibrium have originated clear hypotheses. Species are expected to remain morphologically identical in scenarios of shared genetic variation, such as hybridization and incomplete lineage sorting, or in scenarios where bottlenecks reduce genetic variation and constrain the evolution of morphology. Here, focusing on three morphologically similar Stygocapitella species, we employ a whole-genome amplification method (WGA) coupled with double-digestion restriction-site associated DNA sequencing (ddRAD) to reconstruct the evolutionary history of the species complex. We explore population structure, use population-level statistics to determine the degree of connectivity between populations and species, and determine the most likely demographic scenarios which generally reject for recent hybridization. We find that the combination of WGA and ddRAD allowed us to obtain genomic-level data from microscopic eukaryotes (∼1 millimetre) opening up opportunities for those working with population genomics and phylogenomics in such taxa. The three species share genetic variance, likely from incomplete lineage sorting and ancient admixture. We speculate that the degree of shared variation might underlie morphological similarity in the Atlantic species complex.

A Bioinformatics Tutorial for Comparative Development Genomics in Diverse Meiofauna

Miniaturization, which is a common feature in animals, is particularly manifest in meiofauna-animals sharing peculiar phenotypic features that evolved as adaptations to the highly specialized aquatic interstitial habitat. While revealing much about the extreme phyletic diversity of meiofauna, the genome structure of meiofaunal species could also characterize the phenotype of ancestral states as well as explain the origin and evolution of miniaturization. Here, we present a practical bioinformatics tutorial for genome assembly, genome comparison, and characterization of Hox clusters in meiofaunal species.

ASAP: assemble species by automatic partitioning

Here, we describe Assemble Species by Automatic Partitioning (ASAP), a new method to build species partitions from single locus sequence alignments (i.e., barcode data sets). ASAP is efficient enough to split data sets as large 10⁴ sequences into putative species in several minutes. Although grounded in evolutionary theory, ASAP is the implementation of a hierarchical clustering algorithm that only uses pairwise genetic distances, avoiding the computational burden of phylogenetic reconstruction. Importantly, ASAP proposes species partitions ranked by a new scoring system that uses no biological prior insight of intraspecific diversity. ASAP is a stand-alone program that can be used either through a graphical web-interface or that can be downloaded and compiled for local usage. We have assessed its power along with three others programs (ABGD, PTP and GMYC) on 10 real COI barcode data sets representing various degrees of challenge (from small and easy cases to large and complicated data sets). We also used Monte-Carlo simulations of a multispecies coalescent framework to assess the strengths and weaknesses of ASAP and the other programs. Through these analyses, we demonstrate that ASAP has the potential to become a major tool for taxonomists as it proposes rapidly in a full graphical exploratory interface relevant species hypothesis as a first step of the integrative taxonomy process.

DNA Barcoding for Delimitation of Putative Mexican Marine Nematodes Species

Nematode biodiversity is mostly unknown; while about 20,000 nematode species have been described, estimates for species diversity range from 0.1 to 100 million. The study of nematode diversity, like that of meiofaunal organisms in general, has been mostly based on morphology-based taxonomy, a time-consuming and costly task that requires well-trained specialists. This work represents the first study on the taxonomy of Mexican nematodes that integrates morphological and molecular data. We added eleven new records to the Mexican Caribbean nematode species list: Anticomidae sp.1, Catanema sp.1, Enoploides gryphus, Eurystomina sp.1, Haliplectus bickneri, Metachromadora sp.1, Odontophora bermudensis, Oncholaimus sp.1, Onyx litorale, Proplatycoma fleurdelis, and Pontonema cf. simile. We improved the COI database with 57 new sequences from 20 morphotypes. All COI sequences obtained in this work are new entries for the international genetic databases GenBank and BOLD. Among the studied sites, we report the most extensive species record (12 species) at Cozumel. DNA barcoding and species delineation methods supported the occurrence of 20 evolutionary independent entities and confirmed the high taxonomic resolution of the COI gene. Different approaches provided consistent results: ABGD and mPTP methods disentangled 20 entities, whereas Barcode Index Numbers (BINs) recovered 22 genetic species. Results support DNA barcoding being an efficient, fast, and low-cost method to integrate into morphological observations in order to address taxonomical shortfalls in meiofaunal organisms.

Broad North Atlantic distribution of a meiobenthic annelid - against all odds

DNA barcoding and population genetic studies have revealed an unforeseen hidden diversity of cryptic species among microscopic marine benthos, otherwise exhibiting highly similar and simple morphologies. This has led to a paradigm shift, rejecting cosmopolitism of marine meiofauna until genetically proven and challenging the "Everything is Everywhere, but the environment selects" hypothesis that claims ubiquitous distribution of microscopic organisms. With phylogenetic and species delimitation analyses of worldwide genetic samples of the meiofaunal family Dinophilidae (Annelida) we here resolve three genera within the family and showcase an exceptionally broad, boreal, North Atlantic distribution of a single microscopic marine species with no obvious means of dispersal besides vicariance. With its endobenthic lifestyle, small size, limited migratory powers and lack of pelagic larvae, the broad distribution of Dinophilus vorticoides seems to constitute a "meiofaunal paradox". This species feasts in the biofilm among sand grains, but also on macroalgae and ice within which it can likely survive long-distance rafting dispersal due to its varying lifecycle stages; eggs encapsulated in cocoons and dormant encystment stages. Though often neglected and possibly underestimated among marine microscopic species, dormancy may be a highly significant factor for explaining wide distribution patterns and a key to solving this meiofaunal paradox.

Cephalotrichella echinicola, sp. nov. (Palaeonemertea : Cephalotrichellidae), a new nemertean associated with sea urchins from Nha Trang Bay (South China Sea)

About 50 nemertean species have been reported to live in symbiotic relationships with other invertebrates, but only two hoplonemertean species are associated with echinoderms (starfish). The palaeonemertean Cephalotrichella echinicola, sp. nov. is described from samples collected in Nha Trang Bay, Vietnam, South China Sea. The species is the first known nemertean associated with sea urchins (Metalia sternalis and M. spatagus), living on both the oral and the aboral surfaces of the host and freely moving among its spines. The internal morphology of the new species is described based on histological sections and confocal laser scanning microscopy with phalloidin and antibody labelling. Sequences of three nuclear gene markers (18S rRNA, 28S rRNA, and histone H3) and two mitochondrial gene markers (16S rRNA and COI) were compared with those of other palaeonemertean species, and a phylogenetic analysis suggested that C. echinicola is closest to the free-living Cephalotrichella alba Gibson & Sundberg, 1992. Both the morphological data and the phylogenetic analysis provide additional evidence for distinguishing the families Cephalotrichidae and Cephalotrichellidae and support the rejection of Cephalotrichella as a junior synonym of Cephalothrix. A new diagnosis of the genus Cephalotrichella is given. http://zoobank.org/urn:lsid:zoobank.org:pub:091B5D56-71B2-4F4C-9AD8-F666B4610DE2

Hidden diversity: Phylogeography of genus Ototyphlonemertes Diesing, 1863 (Ototyphlonemertidae: Hoplonemertea) reveals cryptic species and high diversity in Chilean populations

Ototyphlonemertes is a cosmopolitan genus of meiofaunal nemerteans. Their morphological characters are insufficient to reliably identify and delimit species. Consequently, some of the species are considered cosmopolitan despite anticipated low dispersion capability of the adults and a short planktonic larval phase. Indeed, recent studies show that some species actually comprise cryptic species, and populations are connected by stochastic events of long-distance dispersion. Based solely on morphological traits, a Lactea and a Pallida morph of Ototyphlonemertes are recognized here from collections at eight and five locations respectively along the Chilean coast. To assess the phylogeographic patterns of their populations, two mitochondrial markers (COI and COX3) of 162 specimens of Lactea and 25 of Pallida were sequenced. Final sequences are 605bp and 362bp for COI and COX3, respectively. Results from phylogenetic and haplotype network analyses suggest that the Lactea morph comprises up to three independent evolutionary units (one with only COX3 sequences). A COI gene tree including other previously published Ototyphlonemertes sequences groups the Chilean Lactea with other Lactea, while the Chilean Pallida is grouped with other Pallida. Different structuring and gene flow patterns found for the four groups support the hypothesis that these are four independent evolutionary entities with different ecological, dispersal and demographical characteristics.

Genetic diversity of marine oligochaetous clitellates in selected areas of the South Atlantic as revealed by DNA barcoding

Marine oligochaetous clitellates are poorly investigated in the South Atlantic Ocean, especially along the east coast of South America. Closely related species are often difficult to distinguish based on morphology. The lack of specialists and modern identification guides have been pointed out as the main reasons for the scarcity of studies in the South Atlantic Ocean as a whole. To increase the knowledge of this group in the South Atlantic, the genetic diversity of a sample of marine oligochaetous clitellates from Brazil, South Africa and Antarctica was assessed by the Automatic Barcode Gap Discovery (ABGD) and the generalised mixed Yule coalescent (GMYC) approaches. In total, 80 cytochrome c oxidase subunit I (COI) sequences were obtained, each with ~658bp, estimated to represent 32 distinct putative species. ABGD established a barcoding gap between 3% and 14% divergence for uncorrected p-distances and the estimates of GMYC were largely concordant. All the clusters or putative species were genetically associated with previously known species or genera. This study thus confirms the adequacy of the COI barcoding approach combined with a genetic divergence threshold at the order of 10% for marine oligochaetous clitellates.

Species diversity in the marine microturbellarian Astrotorhynchus bifidus sensu lato (Platyhelminthes: Rhabdocoela) from the Northeast Pacific Ocean

Increasing evidence suggests that many widespread species of meiofauna are in fact regional complexes of (pseudo-)cryptic species. This knowledge has challenged the 'Everything is Everywhere' hypothesis and also partly explains the meiofauna paradox of widespread nominal species with limited dispersal abilities. Here, we investigated species diversity within the marine microturbellarian Astrotorhynchus bifidus sensu lato in the Northeast Pacific Ocean. We used a multiple-evidence approach combining multi-gene (18S, 28S, COI) phylogenetic analyses, several single-gene and multi-gene species delimitation methods, haplotype networks and conventional taxonomy to designate Primary Species Hypotheses (PSHs). This included the development of rhabdocoel-specific COI barcode primers, which also have the potential to aid in species identification and delimitation in other rhabdocoels. Secondary Species Hypotheses (SSHs) corresponding to morphospecies and pseudo-cryptic species were then proposed based on the minimum consensus of different PSHs. Our results showed that (a) there are at least five species in the A. bifidus complex in the Northeast Pacific Ocean, four of which can be diagnosed based on stylet morphology, (b) the A. bifidus complex is a mixture of sympatric and allopatric species with regional and/or subglobal distributions, (c) sympatry occurs on local (sample sites), regional (Northeastern Pacific) and subglobal (Northern Atlantic, Arctic, Northeastern Pacific) scales. Mechanisms for this co-occurrence are still poorly understood, but we hypothesize they could include habitat differentiation (spatial and/or seasonal) and life history characteristics such as sexual selection and dispersal abilities. Our results also suggest the need for improved sampling and exploration of molecular markers to accurately map gene flow and broaden our understanding of species diversity and distribution of microturbellarians in particular and meiofauna in general.

Molecular species delimitation in the genus Eumerus (Diptera: Syrphidae)

Eumerus is one of the most diverse genera of hoverfly worldwide. Species delimitation within genus is considered to be difficult due to: (a) lack of an efficient key; (b) non-defined taxonomical status of a large number of species; and (c) blurred nomenclature. Here, we present the first molecular study to delimit species of the genus by using a fragment of the mitochondrial cytochrome-c oxidase subunit I gene (COI) gene. We assessed 75 specimens assigned to 28 taxa originating from two biogeographic zones: 22 from the western Palaearctic and six from the Afrotropical region. Two datasets were generated based on different sequence lengths to explore the significance of availability of more polymorphic sites for species delimitation; dataset A with a total length of 647 bp and dataset B with 746 bp. Various tree inference approaches and Poisson tree processes models were applied to evaluate the putative 'taxonomical' vs. 'molecular' taxa clusters. All analyses resulted in high taxonomic resolution and clear species delimitation for both the dataset lengths. Furthermore, we revealed a high number of mitochondrial haplotypes and high intraspecific variability. We report two major monophyletic clades, and seven 'molecular' groups of taxa formed, which are congruent with morphology-based taxonomy. Our results support the use of the mitochondrial COI gene in species diagnosis of Eumerus.

Phylogeography of the Wheat Stem Sawfly, Cephus cinctus Norton (Hymenoptera: Cephidae): Implications for Pest Management

The wheat stem sawfly, Cephus cinctus Norton (Hymenoptera: Cephidae), is a key pest of wheat in the northern Great Plains of North America, and damage resulting from this species has recently expanded southward. Current pest management practices are inadequate and uncertainty regarding geographic origin, as well as limited data on population structure and dynamics across North America impede progress towards more informed management. We examined the genetic divergence between samples collected in North America and northeastern Asia, the assumed native range of C. cinctus using two mitochondrial regions (COI and 16S). Subsequently, we characterized the structure of genetic diversity in the main wheat producing areas in North America using a combination of mtDNA marker and microsatellites in samples collected both in wheat fields and in grasses in wildlands. The strong genetic divergence observed between North American samples and Asian congeners, in particular the synonimized C. hyalinatus, did not support the hypothesis of a recent American colonization by C. cinctus. Furthermore, the relatively high genetic diversity both with mtDNA and microsatellite markers offered additional evidence in favor of the native American origin of this pest. The genetic diversity of North American populations is structured into three genetic clusters and these are highly correlated with geography. Regarding the recent southern outbreaks in North America, the results tend to exclude the hypothesis of recent movement of damaging wheat stem sawfly populations from the northern area. The shift in host plant use by local populations appears to be the most likely scenario. Finally, the significance of these findings is discussed in the context of pest management.

DNA Taxonomy of Paranemertes (Nemertea: Hoplonemertea) with Spirally Fluted Stylets

Of the 14 nominal species that are now or have ever been assigned to the genus Paranemertes Coe, 1901 , four have been reported to have stylets with a spirally fluted or braided appearance. Although differentiation in color patterns has been documented among species/populations, these nemerteans share similar external characters. Using the sequence datasets of mitochondrial cytochrome c oxidase subunit I (COI), 16S rRNA, and nuclear 28S rRNA genes of specimens from 14 localities of Canada, USA, Russia, Japan, and China, we analyzed the genetic differentiation and reconstructed the phylogenetic trees for these nemerteans. In conjunction with the external characters, we discuss their taxonomy and species delimitation. An analysis based on COI dataset showed high genetic variations among populations and even among worms from the same geographic area. The analyzed 111 individuals were assigned into seven networks by statistical parsimony analysis. The inter-network uncorrected p-distances ranged from 0.044 to 0.172 and the mean intra-network uncorrected p-distances varied from 0.001 to 0.005. With the exception of two networks that contain specimens from the East China Sea, all networks were well-supported by the results of Bayesian and neighbor-joining analyses on the COI data. Phylogenetic trees based on 16S rRNA and 28S rRNA datasets were basically similar to the COI trees, but specimens in some networks were merged into larger clades. Present molecular analyses support the validity of P. sanjuanensis and the synonymization of P. cylindracea with P. peregrina. Nemerteans previously recorded as P. peregrina may contain several species and sympatric speciation might have been occurred in this nemertean group.

Integrative Taxonomy of the Micrura alaskensis Coe, 1901 Species Complex (Nemertea: Heteronemertea), with Descriptions of a New Genus Maculaura gen. nov. and Four New Species from the NE Pacific

Micrura alaskensis Coe, 1901 is a common intertidal heteronemertean known from eastern and northwest Pacific (Alaska to Ensenada, Mexico and Akkeshi, Japan, respectively). It is an emerging model system in developmental biology research. We present evidence from morphology of the adults, gametes, and sequences of cytochrome c oxidase subunit I and 16S rRNA genes that it is not one, but a complex of five, cryptic species. All five of these species co-occur at least in part of their geographic range (e.g. southern Oregon). Preliminary cross-hybridization experiments suggest that at least some of these species are reproductively isolated. The five species share characteristics of adult morphology (e.g. accessory buccal glands) and at least four are known to possess a unique larval morphotype--pilidium maculosum. We propose that these characters define a new genus, Maculaura gen. nov., which contains the following five species: Maculaura alaskensis comb. nov., Maculaura aquilonia sp. nov., Maculaura cerebrosa sp. nov., Maculaura oregonensis sp. nov., and Maculaura magna sp. nov. It is unclear which of the five species Coe originally encountered and described. We chose to retain the name "alaskensis" for the species that current researchers know as "Micrura alaskensis", although, presently, it is only known from Washington and Oregon, and has not been collected from Alaska. Maculaura aquilonia sp. nov. is the only member of the genus we have encountered in Alaska, and we show that it also occurs in the Sea of Okhotsk, Russia.

Trilobodrilus itoi sp. nov., with a Re-Description of T. nipponicus (Annelida: Dinophilidae) and a Molecular Phylogeny of the Genus

The marine interstitial annelid Trilobodrilus itoi sp. nov., the sixth member of the genus, is described on the basis of specimens collected intertidally at Ishikari Beach, Hokkaido, Japan; this is the second species in the genus described from the Pacific Rim. In addition, T. nipponicus Uchida and Okuda, 1943 is re-described based on fresh topotypic material from Akkeshi, Hokkaido, Japan. From both species, we determined sequences of the nuclear 18S and 28S rRNA genes, and the mitochondrial cytochrome c oxidase subunit I (COI) gene. Molecular phylogenetic trees based on concatenated sequences of the three genes showed that T. itoi and T. nipponicus form a clade, which was the sister group to a clade containing the two European congeners T. axi Westheide, 1967 and T. heideri Remane, 1925. The Kimura two-parameter distance for COI was 22.5-22.7% between T. itoi and T. nipponicus, comparable with interspecific values in other polychaete genera. We assessed the taxonomic utility of epidermal inclusions and found that the known six species can be classified into three groups.

Thirty-Five Years of Nemertean (Nemertea) Research--Past, Present, and Future

Developments in nemertean research over the last 35+ years are reviewed from a systematist's perspective. Nemertean systematics and classification, until fairly recently, was not based on explicit phylogenetic hypotheses, but rather on subjective assessment of "important characters". The first cladistic analyses appeared in the 1980s and were criticized at the time by leading researchers in nemertean systematics for not taking into account convergent evolution in ribbon worm morphology. The first molecular study involving the phylum Nemertea appeared in 1992, followed by reports later in the 1990s and early 2000s. Molecular information is now commonplace in nemertean research, and has changed our understanding of evolutionary relationships within the phylum, as well as our view on species and intraspecific variation. Challenges in nemertean systematics and taxonomy are discussed, with special emphasis on future species descriptions, and how to deal with a number of species names that in all likelihood never will be encountered again. Suggestions for how to deal with these challenges are discussed.