Cost-efficient PCR based DNA barcoding of marine invertebrate specimens with NovaSeq amplicon sequencing

BACKGROUND

The marine environment harbors high biodiversity; however, it is poorly understood. Nucleotide sequence data of all marine organisms should be accumulated before natural and/or anthropogenic environmental changes jeopardize the marine environment. In this study, we report a cost-effective and easy DNA barcoding method. This method can be readily adopted without using library preparation kits. It includes multiplex PCR of short targets, indexing PCR, and outsourcing to a sequencing service using the NovaSeq system.

METHODS AND RESULTS

We targeted four mitochondrial genes [cytochrome c oxidase subunit I (COI), COIII, 16S rRNA (16S), and 12S rRNA (12S)] and three nuclear genes [18S rRNA (18S), 28S rRNA (28S), internal transcribed spacer 2 (ITS2)] in 95 marine invertebrate specimens, which were primarily annelids. The primers, including adapters and indices for NovaSeq sequencing, were newly designed. Two PCR runs were conducted. The 1st PCR amplified specific loci with universal primers and the 2nd added sequencing adapters and indices to the 1st PCR products. The gene sequences obtained from the FASTQ files were subjected to BLAST search and phylogenetic analyses. One run using 95 specimens yielded sequences averaging 2816 bp per specimen for a total length of six loci. Nuclear genes were more successfully assembled compared with mitochondrial genes. A weak but significantly negative correlation was observed between the average length of each locus and success rate of the assembly. Some of the sequences were almost identical to the sequences obtained from specimens collected far from Japan, indicating the presence of potentially invasive species identified for the first time.

CONCLUSIONS

We obtained gene sequences efficiently using next-generation sequencing rather than Sanger sequencing. Although this method requires further optimization to increase the success rate for some loci, it is used as a first step to select specimens for further analyses by determining the specific loci of the targets.

Scaphopoda is the sister taxon to Bivalvia: Evidence of ancient incomplete lineage sorting

The almost simultaneous emergence of major animal phyla during the early Cambrian shaped modern animal biodiversity. Reconstructing evolutionary relationships among such closely spaced branches in the animal tree of life has proven to be a major challenge, hindering understanding of early animal evolution and the fossil record. This is particularly true in the species-rich and highly varied Mollusca where dramatic inconsistency among paleontological, morphological, and molecular evidence has led to a long-standing debate about the group's phylogeny and the nature of dozens of enigmatic fossil taxa. A critical step needed to overcome this issue is to supplement available genomic data, which is plentiful for well-studied lineages, with genomes from rare but key lineages, such as Scaphopoda. Here, by presenting chromosome-level genomes from both extant scaphopod orders and leveraging complete genomes spanning Mollusca, we provide strong support for Scaphopoda as the sister taxon of Bivalvia, revitalizing the morphology-based Diasoma hypothesis originally proposed 50 years ago. Our molecular clock analysis confidently dates the split between Bivalvia and Scaphopoda at ~520 Ma, prompting a reinterpretation of controversial laterally compressed Early Cambrian fossils, including Anabarella, Watsonella, and Mellopegma, as stem diasomes. Moreover, we show that incongruence in the phylogenetic placement of Scaphopoda in previous phylogenomic studies was due to ancient incomplete lineage sorting (ILS) that occurred during the rapid radiation of Conchifera. Our findings highlight the need to consider ILS as a potential source of error in deep phylogeny reconstruction, especially in the context of the unique nature of the Cambrian Explosion.

Fusceulimoides kohtsukai gen. et sp. nov., a Minute Eulimid Gastropod Parasitic on the Little Brittle Star Ophiactis savignyi in Central Japan

Gastropods of over a dozen genera in the family Eulimidae have been identified as parasites of brittle stars, and many more remain to be discovered and described for a comprehensive understanding of the evolutionary history of their host-parasite relationships. In this study, we describe Fusceulimoides kohtsukai gen. et sp. nov., parasitic on the little brittle star, Ophiactis savignyi (Ophiactidae), in Kanagawa, central Japan. The new genus is distinguished from other eulimid genera by the combination of following seven conchological characters: (1) a very small size of up to 1.7 mm high, (2) a colorless translucent appearance, (3) a conical pupiform shape with a paucispiral protoconch, (4) slightly convex teleoconch whorls, (5) a remarkably large body whorl occupying 65-70% of the total shell height, (6) a broad, somewhat squarish and laterally expanded aperture with a strongly curved outer lip, and (7) a developed parietal callus without an indentation or depression in the umbilical area. A multi-locus molecular phylogeny revealed its distant relationship to Hemiliostraca + Pyramidelloides, a previously known clade of ophiuroid parasites, thereby suggesting multiple origins of this host-parasite association in Eulimidae.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

Molecular phylogenetic assessment of Spirobranchus kraussii-complex (Annelida: Serpulidae) from the Japanese Archipelago

Spirobranchus kraussii (Annelida: Serpulidae) was recognized as being widely distributed both in the Pacific and Atlantic Oceans. However, the sampling records far from its type locality (South Africa) have been questioned. Actually, recent molecular phylogenetic studies showed that S. kraussii contains genetically distinct species. In this study, we performed molecular phylogenetic analyses of S. cf. kraussii collected from Japan using the nucleotide sequences of a mitochondrial gene and two nuclear genes. Three lineages were recovered within Spirobranchus kraussii-complex in Japan, and one (Spirobranchus sp. 6) showed moderate genetic difference (approximately 4%) in the mitochondrial cytb gene sequence from Spirobranchus sp. 1, an undescribed sequenced species from Honshu Island, Japan. However, the nucleotide sequences of the 18S rRNA gene and ITS2 region were nearly indistinguishable. The other lineage was clearly distinct from the other previously sequenced species and is thus considered to be another distinct species of this species complex (Spirobranchus sp. 5). Although detailed morphological assessment of these lineages is necessary to define their taxonomic status, the present study provided further implications for the species diversity within the S. kraussii-complex.

Complete mitochondrial genomes of two scaphopod molluscs

Complete mitochondrial genomes were determined for two scaphopod molluscs: the dentaliid Antalis entalis and an unidentified Antarctic gadilid. Both genomes are complete except, in Gadilida sp. indet., a short stretch of nad5 was undetermined and trnR could not be annotated. Organization of the Gadilida sp. genome is nearly identical to that previously reported for the gadilid Siphonodentalium whereas trnK, nad5, trnD, nad4, and nad4l are transposed to the opposite strand in the previously published Graptacme genome relative to that of Antalis. Phylogenetic analysis of the 13 protein-coding and 2 rRNA genes recovered Scaphopoda, Gadilida, and Dentaliida monophyletic with maximal support.

Staggered Hox expression is more widespread among molluscs than previously appreciated

Hox genes are expressed along the anterior-posterior body axis in a colinear fashion in the majority of bilaterians. Contrary to polyplacophorans, a group of aculiferan molluscs with conserved ancestral molluscan features, gastropods and cephalopods deviate from this pattern by expressing Hox genes in distinct morphological structures and not in a staggered fashion. Among conchiferans, scaphopods exhibit many similarities with gastropods, cephalopods and bivalves, however, the molecular developmental underpinnings of these similar traits remain unknown. We investigated Hox gene expression in developmental stages of the scaphopod Antalis entalis to elucidate whether these genes are involved in patterning morphological traits shared by their kin conchiferans. Scaphopod Hox genes are predominantly expressed in the foot and mantle but also in the central nervous system. Surprisingly, the scaphopod mid-stage trochophore exhibits a near-to staggered expression of all nine Hox genes identified. Temporal colinearity was not found and early-stage and late-stage trochophores, as well as postmetamorphic individuals, do not show any apparent traces of staggered expression. In these stages, Hox genes are expressed in distinct morphological structures such as the cerebral and pedal ganglia and in the shell field of early-stage trochophores. Interestingly, a re-evaluation of previously published data on early-stage cephalopod embryos and of the gastropod pre-torsional veliger shows that these developmental stages exhibit traces of staggered Hox expression. Considering our results and all gene expression and genomic data available for molluscs as well as other bilaterians, we suggest a last common molluscan ancestor with colinear Hox expression in predominantly ectodermal tissues along the anterior-posterior axis. Subsequently, certain Hox genes have been co-opted into the patterning process of distinct structures (apical organ or prototroch) in conchiferans.

Ecological innovations in the Cambrian and the origins of the crown group phyla

Simulation studies of the early origins of the modern phyla in the fossil record, and the rapid diversification that led to them, show that these are inevitable outcomes of rapid and long-lasting radiations. Recent advances in Cambrian stratigraphy have revealed a more precise picture of the early bilaterian radiation taking place during the earliest Terreneuvian Series, although several ambiguities remain. The early period is dominated by various tubes and a moderately diverse trace fossil record, with the classical ‘Tommotian’ small shelly biota beginning to appear some millions of years after the base of the Cambrian at ca 541 Ma. The body fossil record of the earliest period contains a few representatives of known groups, but most of the record is of uncertain affinity. Early trace fossils can be assigned to ecdysozoans, but deuterostome and even spiralian trace and body fossils are less clearly represented. One way of explaining the relative lack of clear spiralian fossils until about 536 Ma is to assign the various lowest Cambrian tubes to various stem-group lophotrochozoans, with the implication that the groundplan of the lophotrochozoans included a U-shaped gut and a sessile habit. The implication of this view would be that the vagrant lifestyle of annelids, nemerteans and molluscs would be independently derived from such a sessile ancestor, with potentially important implications for the homology of their sensory and nervous systems.

The ParaHox gene Gsx patterns the apical organ and central nervous system but not the foregut in scaphopod and cephalopod mollusks

Background

It has been hypothesized that the ParaHox gene Gsx patterned the foregut of the last common bilaterian ancestor. This notion was corroborated by Gsx expression in three out of four lophotrochozoan species, several ecdysozoans, and some deuterostomes. Remarkably, Gsx is also expressed in the bilaterian anterior-most central nervous system (CNS) and the gastropod and annelid apical organ. To infer whether these findings are consistent with other mollusks or even lophotrochozoans, we investigated Gsx expression in developmental stages of representatives of two other molluscan classes, the scaphopod Antalis entalis and the cephalopod Idiosepius notoides.

Results

Gsx is not expressed in the developing digestive tract of Antalis entalis and Idiosepius notoides. Instead, it is expressed in cells of the apical organ in the scaphopod trochophore and in two cells adjacent to this organ. Late-stage trochophores express Aen-Gsx in cells of the developing cerebral and pedal ganglia and in cells close to the pavilion, mantle, and foot. In postmetamorphic specimens, Aen-Gsx is expressed in the cerebral and pedal ganglia, the foot, and the nascent captacula. In early squid embryos, Ino-Gsx is expressed in the cerebral, palliovisceral, and optic ganglia. In late-stage embryos, Ino-Gsx is additionally expressed close to the eyes and in the supraesophageal and posterior subesophageal masses and optic lobes. Developmental stages close to hatching express Ino-Gsx only close to the eyes.

Conclusions

Our results suggest that Gsx expression in the foregut might not be a plesiomorphic trait of the Lophotrochozoa as insinuated previously. Since neither ecdysozoans nor deuterostomes express Gsx in their gut, a role in gut formation in the last common bilaterian ancestor appears unlikely. Gsx is consistently expressed in the bilaterian anterior-most CNS and the apical organ of lophotrochozoan larvae, suggesting a recruitment of Gsx into the formation of this organ in the Lophotrochozoa. The cephalopod posterior subesophageal mass and optic ganglia and the scaphopod pedal ganglia also express Gsx. In summary, Gsx expression only appears to be conserved in the anterior-most brain region during evolution. Accordingly, Gsx appears to have been recruited into the formation of other expression domains, e.g., the apical organ or the foregut, in some lophotrochozoans.

The continuing debate on deep molluscan phylogeny: evidence for Serialia (Mollusca, Monoplacophora + Polyplacophora)

Molluscs are a diverse animal phylum with a formidable fossil record. Although there is little doubt about the monophyly of the eight extant classes, relationships between these groups are controversial. We analysed a comprehensive multilocus molecular data set for molluscs, the first to include multiple species from all classes, including five monoplacophorans in both extant families. Our analyses of five markers resolve two major clades: the first includes gastropods and bivalves sister to Serialia (monoplacophorans and chitons), and the second comprises scaphopods sister to aplacophorans and cephalopods. Traditional groupings such as Testaria, Aculifera, and Conchifera are rejected by our data with significant Approximately Unbiased (AU) test values. A new molecular clock indicates that molluscs had a terminal Precambrian origin with rapid divergence of all eight extant classes in the Cambrian. The recovery of Serialia as a derived, Late Cambrian clade is potentially in line with the stratigraphic chronology of morphologically heterogeneous early mollusc fossils. Serialia is in conflict with traditional molluscan classifications and recent phylogenomic data. Yet our hypothesis, as others from molecular data, implies frequent molluscan shell and body transformations by heterochronic shifts in development and multiple convergent adaptations, leading to the variable shells and body plans in extant lineages.

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

Comparative morphology among representatives of main taxa of Scaphopoda and basal protobranch Bivalvia (Mollusca)

This study deals with detailed morphology and anatomy of 4 species of Scaphopoda and 5 species of protobranch Bivalvia. Both classes are traditionally grouped in the taxon Diasoma, which has been questioned by different methodologies, such as molecular and developmental. This study is developed under a phylogenetic methodology with the main concern in performing it in an intelligible and testable methodology. The analyzed Scaphopoda species came from the Brazilian coast and belong to the family Dentaliidae [(1) Coccodentalium carduus; (2) Paradentalium disparile] and Gadiliidae; [(3) Polyschides noronhensis, n. sp. from Fernando de Noronha Archipelago; (4) Gadila braziliensis]. These species represent the main branches of the class Scaphopoda. From protobranch bivalves, representatives of the families Solemyidae [(5) Solemya occidentalis, from Florida; S. notialis, n. sp. from S.E. Brazil], Nuculanidae [(6) Propeleda carpentieri from Florida], and Nuculidae [(7) Ennucula puelcha, from south Brazil] are included. These species represent the main branches of the basal Bivalvia. The descriptions on the anatomy of S. occidentalis and of P. carpentieri are published elsewhere. The remaining are included here, for which a complete taxonomical treatment is performed. Beyond these species, representatives of other taxa are operationally included as part of the ingroup (indices are then shared with them), as a procedure to test the morphological monophyly of Diasoma. These taxa are: two lamellibranch bivalves [(8) Barbatia - Arcidae; (9) Serratina - Tellinidae; both published elsewhere;, and Propilidium (10) Patellogastropoda, and (11) Nautilus, basal Cephalopoda, based on basal taxa. The effective outgroups are (12) Neopilina (Monoplacophora) and (13) Hanleya (Polyplacophora). The phylogenetic analysis based on morphology revealed that the taxon Diasoma is supported by 14 synapomorphies, and is separated from Cyrtosoma (Gastropoda + Cephalopoda). Although they are not the main goal of this paper, the taxa Scaphopoda and Bivalvia are supported by 8 and by 7 synapomorphies respectively. The taxon Protobranchia resulted paraphyletic. Both scaphopod orders resulted monophyletic. The obtained cladogram is: ((((Coccodentalium carduus - Paradentalium disparile) (Polyschides noronhensis - Gadila brasiliensis)) ((Solemya occidentalis - S. notialis) (Propeleda carpenteri (Ennucula puelcha (Barbatia cancellaria - Serratina capsoides))))) (Propilidium curumim - Nautilus pompilius - Lolliguncula brevis)).

Evidence for a clade composed of molluscs with serially repeated structures: monoplacophorans are related to chitons

Monoplacophorans are among the rarest members of the phylum Mollusca. Previously only known from fossils since the Cambrian, the first living monoplacophoran was discovered during the famous second Galathea deep-sea expedition. The anatomy of these molluscs shocked the zoological community for presenting serially repeated gills, nephridia, and eight sets of dorsoventral pedal retractor muscles. Seriality of organs in supposedly independent molluscan lineages, i.e., in chitons and the deep-sea living fossil monoplacophorans, was assumed to be a relic of ancestral molluscan segmentation and was commonly accepted to support a direct relationship with annelids. We were able to obtain one specimen of a monoplacophoran Antarctic deep-sea species for molecular study. The first molecular data on monoplacophorans, analyzed together with the largest data set of molluscs ever assembled, clearly illustrate that monoplacophorans and chitons form a clade. This "Serialia" concept may revolutionize molluscan systematics and may have important implications for metazoan evolution as it allows for new interpretations for primitive segmentation in molluscs.

Complete DNA sequence of the mitochondrial genome of the sea-slug, Aplysia californica: conservation of the gene order in Euthyneura

We have sequenced and characterized the complete mitochondrial genome of the sea slug, Aplysia californica, an important model organism in experimental biology and a representative of Anaspidea (Opisthobranchia, Gastropoda). The mitochondrial genome of Aplysia is in the small end of the observed sizes of animal mitochondrial genomes (14,117 bp, NCBI Accession No. NC_005827). The Aplysia genome, like most other mitochondrial genomes, encodes genes for 2 ribosomal subunit RNAs (small and large rRNAs), 22 tRNAs, and 13 protein subunits (cytochrome c oxidase subunits 1-3, cytochrome b apoenzyme, ATP synthase subunits 6 and 8, and NADH dehydrogenase subunits 1-6 and 4L). The gene order is virtually identical between opisthobranchs and pulmonates, with the majority of differences arising from tRNA translocations. In contrast, the gene order from representatives of basal gastropods and other molluscan classes is significantly different from opisthobranchs and pulmonates. The Aplysia genome was compared to all other published molluscan mitochondrial genomes and phylogenetic analyses were carried out using a concatenated protein alignment. Phylogenetic analyses using maximum likelihood based analyses of the well aligned regions of the protein sequences support both monophyly of Euthyneura (a group including both the pulmonates and opisthobranchs) and Opisthobranchia (as a more derived group). The Aplysia mitochondrial genome sequenced here will serve as an important platform in both comparative and neurobiological studies using this model organism.

Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences

The Mollusca represent one of the most morphologically diverse animal phyla, prompting a variety of hypotheses on relationships between the major lineages within the phylum based upon morphological, developmental, and paleontological data. Analyses of small-ribosomal RNA (SSU rRNA) gene sequence have provided limited resolution of higher-level relationships within the Mollusca. Recent analyses suggest large-subunit (LSU) rRNA gene sequences are useful in resolving deep-level metazoan relationships, particularly when combined with SSU sequence. To this end, LSU (approximately 3.5 kb in length) and SSU (approximately 2 kb) sequences were collected for 33 taxa representing the major lineages within the Mollusca to improve resolution of intraphyletic relationships. Although the LSU and combined LSU+SSU datasets appear to hold potential for resolving branching order within the recognized molluscan classes, low bootstrap support was found for relationships between the major lineages within the Mollusca. LSU+SSU sequences also showed significant levels of rate heterogeneity between molluscan lineages. The Polyplacophora, Gastropoda, and Cephalopoda were each recovered as monophyletic clades with the LSU+SSU dataset. While the Bivalvia were not recovered as monophyletic clade in analyses of the SSU, LSU, or LSU+SSU, the Shimodaira-Hasegawa test showed that likelihood scores for these results did not differ significantly from topologies where the Bivalvia were monophyletic. Analyses of LSU sequences strongly contradict the widely accepted Diasoma hypotheses that bivalves and scaphopods are closely related to one another. The data are consistent with recent morphological and SSU analyses suggesting scaphopods are more closely related to gastropods and cephalopods than to bivalves. The dataset also presents the first published DNA sequences from a neomeniomorph aplacophoran, a group considered critical to our understanding of the origin and early radiation of the Mollusca.

The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca)

Comparisons of mitochondrial gene sequences and gene arrangements can be informative for reconstructing high-level phylogenetic relationships. We determined the complete sequence of the mitochondrial genome of Siphonodentalium lobatum, (Mollusca, Scaphopoda). With only 13,932 bases, it is the shortest molluscan mitochondrial genome reported so far. The genome contains the usual 13 protein-coding genes, two rRNA and 22 tRNA genes. The ATPase subunit 8 gene is exceptionally short. Several transfer RNAs show truncated TpsiC arms or DHU arms. The gene arrangement of S. lobatum is markedly different from all other known molluscan mitochondrial genomes and shows low similarity even to an unpublished gene order of a dentaliid scaphopod. Phylogenetic analyses of all available complete molluscan mitochondrial genomes based on amino acid sequences of 11 protein-coding genes yield trees with low support for the basal branches. None of the traditionally accepted molluscan taxa and phylogenies are recovered in all analyses, except for the euthyneuran Gastropoda. S. lobatum appears as the sister taxon to two of the three bivalve species. We conclude that the deep molluscan phylogeny is probably beyond the resolution of mitochondrial protein sequences. Moreover, assessing the phylogenetic signal in gene order data requires a much larger taxon sample than is currently available, given the exceptional diversity of this character set in the Mollusca.