The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea)

Ribbon worms (phylum Nemertea) from Bodega Bay, California: A largely undescribed diversity

The diversity of nemerteans along the Pacific coast of the United States is regarded as well characterized, but there remain many cryptic, undescribed, and "orphan" species (those known only in their larval form). Recent sampling of nemerteans in Oregon and Washington has begun to fill in these taxonomic gaps, but nemertean diversity in California has received relatively little attention over the past 60 years. During the summers of 2019 and 2020, nemertean specimens were collected from 20 locations in the Bodega Bay region of northern California, USA, including rocky intertidal shores, sandy beaches, mudflats, and other habitats. Based on morphological assessment and DNA sequence analysis (partial Cytochrome Oxidase I and 16S rRNA genes), our surveys identified 34 nemertean species. Only 13 of these (38%) can be confidently assigned to described species. Another 11 represent species that are new to science, including members of the genera Riserius, Nipponnemertes, Poseidonemertes, Zygonemertes, Nemertellina, Oerstedia, and three species of uncertain affiliation. The remaining ten species include undescribed or cryptic species of uncertain status that have been found previously along the Pacific Coast of North America. Our surveys also document extensions of known geographic ranges for multiple species, including the first records in California of Antarctonemertesphyllospadicola, Cephalothrixhermaphroditica, and Maculauraoregonensis. This is the first report of the genus Nemertellina in the northeast Pacific and Riserius in California. Overall, our findings highlight how much remains to be learned about the diversity and distribution of nemerteans in the northeast Pacific.

Characterization of the complete mitochondrial genome of Quasilineus sinicus Gibson, 1990 (Nemertea: Heteronemertea) and its phylogenetic implications

In this study, we sequenced and characterized the complete mitochondrial genome (mitogenome) of Quasilineus sinicus Gibson, 1990 (Heteronemertea, Nemertea) using Illumina sequencing technology. The circular mitogenome was 16,358 bp in length and comprised 22 transfer RNA genes, 13 protein-coding genes, and two ribosomal RNA genes. Its overall base composition included 20.82% A, 41.06% T, 26.68% G, and 11.44% C; in fact, the mitogenome had a high A + T content of 61.88%. Furthermore, our phylogenetic analysis demonstrated that Paleonemertea, Pilidiophora, and Hoplonemertea were monophyletic groups, and Q. sinicus was most closely related to Iwatanemertes piperata.

Mitochondrial gene order of the freshwater bryozoan Cristatella mucedo retains ancestral lophotrochozoan features

Bryozoans are aquatic colonial suspension-feeders abundant in many marine and freshwater benthic communities. At the same time, the phylum is under studied on both morphological and molecular levels, and its position on the metazoan tree of life is still disputed. Bryozoa include the exclusively marine Stenolaemata, predominantly marine Gymnolaemata and exclusively freshwater Phylactolaemata. Here we report the mitochondrial genome of the phylactolaemate bryozoan Cristatella mucedo. This species has the largest (21,008 bp) of all currently known bryozoan mitogenomes, containing a typical metazoan gene compendium as well as a number of non-coding regions, three of which are longer than 1500 bp. The trnS1/trnG/nad3 region is presumably duplicated in this species. Comparative analysis of the gene order in C. mucedo and another phylactolaemate bryozoan, Pectinatella magnifica, confirmed their close relationships, and revealed a stronger similarity to mitogenomes of phoronids and other lophotrochozoan species than to marine bryozoans, indicating the ancestral nature of their gene arrangement. We suggest that the ancestral gene order underwent substantial changes in different bryozoan cladesshowing mosaic distribution of conservative gene blocks regardless of their phylogenetic position. Altogether, our results support the early divergence of Phylactolaemata from the rest of Bryozoa.

Comparative analysis of complete mitochondrial genomes with Cerithioidea and molecular phylogeny of the freshwater snail, Semisulcospira gottschei (Caenogastropoda, Cerithioidea)

In this study, we determined the complete mitochondrial genome (mt genome) of Semisulcospira gottschei for the first time and then compared it with the mt genome of species belonging to Cerithioidea. The mt genome consists of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and non-coding region with a total length of 16,101 bp. The type of constitutive genes and the direction of the coding strand which appeared in the mt genome were the same as the ones observed in Cerithioidea except for the tRNA-Q and tRNA-R positions. The S. gottschei mt genome had a non-coding region with an AT-rich loop between tRNA-F and tRNA-C regions. In regard to molecular phylogeny, two types of analysis were performed to confirm the introgressive hybridization of S. gottschei and to identify the phylogenetic location among the species in Caenogastropoda. As a result, S. gottschei used in this study belonged to the same clade as other non-introgressed S. gottschei. As for the molecular phylogenic analysis of species belonging to Caenogastropoda, S. gottschei was found to be the closest to S. coreana taxonomically and to be included in Cerithioidea.

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.

Mitochondrial genome of Micrura bella (Nemertea: Heteronemertea), the largest mitochondrial genome known to phylum Nemertea

The complete mitochondrial genome (mitogenome) of Micrura bella was sequenced and analyzed. Being the largest mitogenome known to phylum Nemertea, the genome is 16 847 bp in length. It encodes 37 genes typical to metazoan mitogenomes and has the same gene arrangement with the other Heteronemertea mitogenomes sequenced to date. The genome has the maximal number of non-coding nucleotides (2037 bp at 25 sites) in Nemertea mitogenomes, among which two large non-coding regions were found (507 and 508 bp, respectively).

The complete mitochondrial genome of Scutopus ventrolineatus (Mollusca: Chaetodermomorpha) supports the Aculifera hypothesis

Background

With more than 100000 living species, mollusks are the second most diverse metazoan phylum. The current taxonomic classification of mollusks recognizes eight classes (Neomeniomorpha, Chaetodermomorpha, Polyplacophora, Monoplacophora, Cephalopoda, Gastropoda, Bivalvia, and Scaphopoda) that exhibit very distinct body plans. In the past, phylogenetic relationships among mollusk classes have been contentious due to the lack of indisputable morphological synapomorphies. Fortunately, recent phylogenetic analyses based on multi-gene data sets are rendering promising results. In this regard, mitochondrial genomes have been widely used to reconstruct deep phylogenies. For mollusks, complete mitochondrial genomes are mostly available for gastropods, bivalves, and cephalopods, whereas other less-diverse lineages have few or none reported.

Results

The complete DNA sequence (14662 bp) of the mitochondrial genome of the chaetodermomorph Scutopus ventrolineatus Salvini-Plawen, 1968 was determined. Compared with other mollusks, the relative position of protein-coding genes in the mitochondrial genome of S. ventrolineatus is very similar to those reported for Polyplacophora, Cephalopoda and early-diverging lineages of Bivalvia and Gastropoda (Vetigastropoda and Neritimorpha; but not Patellogastropoda). The reconstructed phylogenetic tree based on combined mitochondrial and nuclear sequence data recovered monophyletic Aplacophora, Aculifera, and Conchifera. Within the latter, Cephalopoda was the sister group of Gastropoda and Bivalvia + Scaphopoda.

Conclusions

Phylogenetic analyses of mitochondrial sequences showed strong among-lineage rate heterogeneity that produced long-branch attraction biases. Removal of long branches (namely those of bivalves and patellogastropods) ameliorated but not fully resolved the problem. Best results in terms of statistical support were achieved when mitochondrial and nuclear sequence data were concatenated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0197-9) contains supplementary material, which is available to authorized users.

Complete mitochondrial genome sequences of two parasitic/commensal nemerteans, Gononemertes parasita and Nemertopsis tetraclitophila (Nemertea: Hoplonemertea)

BACKGROUND

Most nemerteans (phylum Nemertea) are free-living, but about 50 species are known to be firmly associated with other marine invertebrates. For example, Gononemertes parasita is associated with ascidians, and Nemertopsis tetraclitophila with barnacles. There are 12 complete or near-complete mitochondrial genome (mitogenome) sequences of nemerteans available in GenBank, but no mitogenomes of none free-living nemerteans have been determined so far. In the present paper complete mitogenomes of the above two parasitic/commensal nemerteans are reported.

METHODS

The complete mitochondrial genomes (mitogenome) of G. parasita and N. tetraclitophila were amplified by conventional and long PCR. Phylogenetic analyses of maximum likelihood (ML) and Bayesian inference (BI) were performed with both concatenated nucleotide and amino acid sequences.

RESULTS

Complete mitogenomes of G. parasita and N. tetraclitophila are 14742 bp and 14597 bp in size, respectively, which are within the range of published Hoplonemertea mitogenomes. Their gene orders are identical to that of published Hoplonemertea mitogenomes, but different from those of Palaeo- and Heteronemertea species. All the coding genes, as well as major non-coding regions (mNCRs), are AT rich, which is especially pronounced at the third codon position. The AT/GC skew pattern of the coding strand is the same among nemertean mitogenomes, but is variable in the mNCRs. Some slight differences are found between mitogenomes of the present species and other hoplonemerteans: in G. parasita the mNCR is biased toward T and C (contrary to other hoplonemerteans) and the rrnS gene has a unique 58-bp insertion at the 5' end; in N. tetraclitophila the nad3 gene starts with the ATT codon (ATG in other hoplonemerteans). Phylogenetic analyses of the nucleotide and amino acid datasets show early divergent positions of G. parasita and N. tetraclitophila within the analyzed Distromatonemertea species, and provide strong support for the close relationship between Hoplonemertea and Heteronemertea.

CONCLUSION

Gene order is highly conserved within the order Monostilifera, particularly within the Distromatonemertea, and the special lifestyle of G. parasita and N. tetraclitophila does not bring significant variations to the overall structures of their mitogenomes in comparison with free-living hoplonemerteans.

A description of the complete mitochondrial genomes of Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula (Nemertea: Hoplonemertea: Monostilifera)

We sequenced the complete mitochondrial genomes (mitogenomes) of three Hoplonemertea species, Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula, which are 14,616, 14,655 and 15,354 bp in length, respectively. Each of the three circular mitogenomes consists of 37 typical genes and some non-coding regions. The nucleotide composition of the coding strand is biased toward T, almost a half of total nucleotides in these mitogenomes. There are many poly-T tracts across these mitogenomes, which exhibit T-number variation within different clones of protein-coding genes, mainly resulting from false PCR amplification. The major non-coding regions have tandem repeat motifs and hairpin-like structures that may be associated with the initiation of replication or transcription. Data published to date for nemerteans show that Palaeonemertea species usually bear the largest mitogenomes, while representatives in the more recently derived Distromatonemertea clade bear the smallest ones; and that the gene arrangement of mitogenomes seems to be variable within the phylum Nemertea, but stable within either of Heteronemertea and Hoplonemertea.

Taxonomic identity of a tetrodotoxin-accumulating ribbon-worm Cephalothrix simula (Nemertea: Palaeonemertea): a species artificially introduced from the Pacific to Europe

We compared the anatomy of the holotype of the palaeonemertean Cephalothrix simula ( Iwata, 1952 ) with that of the holotypes of Cephalothrix hongkongiensis Sundberg, Gibson and Olsson, 2003 and Cephalothrix fasciculus ( Iwata, 1952 ), as well as additional specimens from Fukue (type locality of C. simula) and Hiroshima, Japan. While there was no major morphological discordance between these specimens, we found discrepancies between the actual morphology and some statements in the original description of C. simula with respect to supposedly species-specific characters. Our observation indicates that these three species cannot be discriminated by the anatomical characters so far used to distinguish congeners. For objectivity of scientific names, topogenetypes of the mitochondrial cytochrome c oxidase subunit I (COI) sequences are designated for C. simula, C. hongkongiensis, and C. fasciculus. Analysis of COI sequence showed that the Hiroshima population can be identified as C. simula, which has been found in previous studies from Trieste, Italy, and also from both the Mediterranean and Atlantic coasts of the Iberian Peninsula, indicating an artificial introduction via (1) ballast water, (2) ship-fouling communities, or (3) the commercially cultured oyster Crassostrea gigas ( Thunberg, 1793 ) brought from Japan to France in 1970s. Cephalothrix simula is known to be toxic, as it contains large amounts of tetrodotoxin (TTX). We report here that the grass puffer Takifugu niphobles ( Jordan and Snyder, 1901 )-also known to contain TTX- consumes C. simula. We suggest that the puffer may be able to accumulate TTX by eating C. simula.

Mutation and selection cause codon usage and bias in mitochondrial genomes of ribbon worms (Nemertea)

The phenomenon of codon usage bias is known to exist in many genomes and it is mainly determined by mutation and selection. To understand the patterns of codon usage in nemertean mitochondrial genomes, we use bioinformatic approaches to analyze the protein-coding sequences of eight nemertean species. Neutrality analysis did not find a significant correlation between GC12 and GC3. ENc-plot showed a few genes on or close to the expected curve, but the majority of points with low-ENc values are below it. ENc-plot suggested that mutational bias plays a major role in shaping codon usage. The Parity Rule 2 plot (PR2) analysis showed that GC and AT were not used proportionally and we propose that codons containing A or U at third position are used preferentially in nemertean species, regardless of whether corresponding tRNAs are encoded in the mitochondrial DNA. Context-dependent analysis indicated that the nucleotide at the second codon position slightly affects synonymous codon choices. These results suggested that mutational and selection forces are probably acting to codon usage bias in nemertean mitochondrial genomes.

The complete mitochondrial genome of Pseudochauhanea macrorchis (Monogenea: Chauhaneidae) revealed a highly repetitive region and a gene rearrangement hot spot in Polyopisthocotylea

The complete mitochondrial genome of Pseudochauhanea macrorchis was determined and compared with other monogenean mitochondrial genomes from GenBank. The circular genome was 15,031 bp in length and encoded 36 genes (12 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs) typically found in flatworms. Structures of the mitochondrial genome were mostly concordant with that known for Microcotyle sebastis and Polylabris halichoeres, but also contained two noted features-a gene rearrangement hot spot and the highly repetitive region (HRR) in major non-coding region (NCR). The gene rearrangement hot spot located between the cox3 and nad5 genes, including a cluster of tRNA genes, nad6 gene and one major NCR. The HRR seemed to be a unique feature of the polyopisthocotylean mitochondrial genomes. In conclusion, the present study provided new molecular data for future studies of the comparative mitochondrial genomics and also served as a resource of markers for the studies of species populations and monogenean phylogenetics.

A comparative study of nemertean complete mitochondrial genomes, including two new ones for Nectonemertes cf. mirabilis and Zygeupolia rubens, may elucidate the fundamental pattern for the phylum Nemertea

BACKGROUND

The mitochondrial genome is important for studying genome evolution as well as reconstructing the phylogeny of organisms. Complete mitochondrial genome sequences have been reported for more than 2200 metazoans, mainly vertebrates and arthropods. To date, from a total of about 1275 described nemertean species, only three complete and two partial mitochondrial DNA sequences from nemerteans have been published. Here, we report the entire mitochondrial genomes for two more nemertean species: Nectonemertes cf. mirabilis and Zygeupolia rubens.

RESULTS

The sizes of the entire mitochondrial genomes are 15365 bp for N. cf. mirabilis and 15513 bp for Z. rubens. Each circular genome contains 37 genes and an AT-rich non-coding region, and overall nucleotide composition is AT-rich. In both species, there is significant strand asymmetry in the distribution of nucleotides, with the coding strand being richer in T than A and in G than C. The AT-rich non-coding regions of the two genomes have some repeat sequences and stem-loop structures, both of which may be associated with the initiation of replication or transcription. The 22 tRNAs show variable substitution patterns in nemerteans, with higher sequence conservation in genes located on the H strand. Gene arrangement of N. cf. mirabilis is identical to that of Paranemertes cf. peregrina, both of which are Hoplonemertea, while that of Z. rubens is the same as in Lineus viridis, both of which are Heteronemertea. Comparison of the gene arrangements and phylogenomic analysis based on concatenated nucleotide sequences of the 12 mitochondrial protein-coding genes revealed that species with closer relationships share more identical gene blocks.

CONCLUSION

The two new mitochondrial genomes share many features, including gene contents, with other known nemertean mitochondrial genomes. The tRNA families display a composite substitution pathway. Gene order comparison to the proposed ground pattern of Bilateria and some lophotrochozoans suggests that the nemertean ancestral mitochondrial gene order most closely resembles the heteronemertean type. Phylogenetic analysis proposes a sister-group relationship between Hetero- and Hoplonemertea, which supports one of two recent alternative hypotheses of nemertean phylogeny.