- Open Access
Complete mitochondrial genome sequences of two parasitic/commensal nemerteans, Gononemertes parasita and Nemertopsis tetraclitophila (Nemertea: Hoplonemertea)
Parasites & Vectorsvolume 7, Article number: 273 (2014)
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.
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.
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.
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.
The phylum Nemertea (ribbon worm) includes about 1280 named species . Most of them are free-living in marine, freshwater and terrestrial habitats, but there are about 50 species reported to be associated with other animals; host organisms include poriferans, cnidarians, bivalves, echiurans, crustaceans, echinoderms and ascidians. The position of Nemertea among metazoans was traditionally considered to be close to the acoelomate Platyhelminthes, but comparative ultrastructure studies and molecular phylogenetic analyses during recent decades have supported it to be a member of the Lophotrochozoa [2–6]. The phylogenetic relationship of the phylum is still unsettled in parts, and conclusions may be dependent on different markers and analytical methods [7–9]. A recent analysis based on four nuclear and two mitochondrial loci further suggested that an expanded taxon sampling at family and generic level was required for getting a better understanding of nemertean affinities .
To date, there are 12 complete or near-complete nemertean mitogenome sequences available in GenBank. From these, we can infer some interesting patterns in terms of genome organization. For instance, Palaeonemertea and Heteronemertea bear larger mitogenomes than the more recently diverged hoplonemertean taxon Distromatonemertea. The gene arrangement within the phylum is not conserved, but generally stable within each of the three major groups (Palaeo-, Hetero- and Hoplonemertea). Nevertheless, a fuller understanding of the evolutionary patterns of nemertean mitogenome evolution requires denser taxon sampling, particularly of taxa that have adopted unusual lifestyles, such as Malacobdella and Carcinonemertes. In the present study, we determined the first complete mitogenome sequences of two parasitic/commensal nemerteans, Gononemertes parasita Bergendal, 1900 and Nemertopsis tetraclitophila Gibson, 1990, which taxonomically belong to Monostilifera (a group that contains most known symbiotic nemerteans). G. parasita lives in the branchial chamber of some ascidians in European waters , whereas N. tetraclitophila has been recorded from the mantle cavity of the balanomorph barnacle Tetraclita squamosa (Bruguiére, 1789) in Hong Kong, China . Worms of both species seem to be firmly associated with a host, and possess some adaptive features that might be related to none free-living lifestyle, e.g., the greater number of gonads than most free-living monostiliferans; the absence of a proboscis apparatus (G. parasita) [11, 12]. Mostly based on reproductive adaptations, Roe has argued that G. parasita and another Nemertopis species living in barnacles (Nemertopis quadripunctata (Quoy & Gaimard, 1833), which feeds on the eggs of the barnacles and possibly on the barnacles themselves) should be regarded as parasites . However, the ecology, particularly the feeding biology, of G. parasita and N. tetraclitophila has not been well understood. Therefore, the two species are cautiously mentioned as “parasitic/commensal” in the present paper.
Specimens and DNA extraction
Gononemertes parasita was collected from the branchial chamber of the sea squirt Ascidia obliqua Alder, 1863 near Tjärnö, Sweden. Nemertopsis tetraclitophila was collected from the mantle cavity of the barnacle Tetraclita squamosa in Shenzhen, China. For either species, total DNA was extracted from a single specimen using the Genomic DNA Extraction Kit (OMEGA) following the manufacturer’s instructions and stored at −20°C.
PCR amplification and sequencing
Small fragments such as cox1, rrnS-rrnL, cob and cox3 were amplified with universal primers, and then specific primers were designed for the amplification of long fragments (Additional file 1: Table S1). All PCR reactions were carried out in a reaction volume of 25 μl containing 12.5 μl Premix Taq (LA version 2.0) (TaKaRa Clone Tech), 0.5 μl each primer, 0.5 μl DNA template and 11 μl distilled H2O. The PCR amplifications were performed under the following conditions: 4 min at 94°C, followed by 35 cycles of 30 s at 94°C, 30 s at 48–50°C (according to primers), 1–10 min (according to the length of products) at 72°C, followed by a 10 min elongation. The PCR products were separated by agarose gel electrophoresis and purified using DNA gel extraction kit (OMEGA). The purified PCR products were ligated into pEASY-T1 vector (Transgen, China) and sequenced by primer walking on an ABI 3730 Sequencer.
Genome assembly and annotation
All the sequences were compared with other nemerteans to prevent contaminations from a host or bacteria. The obtained fragments of mitogenomes were assembled with Codoncode Aligner 5.0.1. Identification of protein-coding genes and rRNA genes was performed by BLAST searches (http://www.ncbi.nlm.nih.gov/BLAST) and by alignment to known hoplonemertean mitogenomes. Most tRNA genes were identified by tRNAscan-SE 1.21 , and additional tRNA genes were inferred with RNAfold . The mitogenome was visualized using CGView . The nucleotide composition and codon usage were calculated with DAMBE 5 . Multiple alignments of genes were generated by Clustal X  with default settings and amino acid translation was carried out using MEGA 5.0 . The full mitogenome sequences of Gononemertes parasita [KF572481] and Nemertopsis tetraclitophila [KF572482] were submitted to GenBank and compared with Cephalothrix hongkongiensis [NC_012821], Cephalothrix sp. [NC_014869], Iwatanemertes piperata [KF719984], Lineus viridis [NC_012889], Lineus alborostratus [NC_018356], Nectonemertes cf. mirabilis [NC_017874], Amphiporus formidabilis [KC710979], Emplectonema gracile [NC_016952], Paranemertes cf. peregrina [NC_014865], Zygeupolia rubens [NC_017877], Prosadenoporus spectaculum [KC710980] and Nipponnemertes punctatula [KC710981].
Phylogenetic analyses of the 14 available nemertean mitogenomes were carried out as follows: i) nucleotide-level analysis of protein-coding genes, with 3rd codon position removed; ii) nucleotide-level analysis of protein-coding genes, with 3rd codon position removed, rRNA and tRNA genes, iii) amino acid-level analysis of protein-coding genes. The saturation test was carried out based on the transition and transversion substitutions vs. the Tamura-Nei (TN93) distance of three codon positions by DAMBE 5 , and the third codon position which tended to be saturated (the transition and transversion substitution values do not increase as the genetic distance increase) was not used in phylogenetic analyses. The outgroups Katharina tunicata [NC_001636] and Terebratulina retusa [NC_000941] were selected based on their close relationships with Nemertea in previous studies [20, 21]. All datasets were aligned with Clustal X with default settings . Poorly aligned positions were excluded using Gblocks Version 0.91b  allowing less strict flanking positions and other default parameters. For nucleotide sequences MODELTEST  and MRMODELTEST , and for amino acid sequences ProtTest 2.4  were used to select the best-fit substitution models (the model parameters were estimated when the concatenated nucleotides/amino acids were treated as a single partition). Based on the Akaike Information Criterion (AIC), the best-fit model for nucleotides was the GTR + I + G and for amino acid sequences was the MtRev + G + F. The ML analysis was performed with PHYML 3.0 program (http://www.atgc-montpellier.fr/phyml/)  with 100 bootstrap replicates. Bayesian inference was conducted using MrBayes version 3.1.2 . Four Monte Carlo Markov chains (MCMC) were run for 1,000,000 generations, sampling every 100 generations. The first 2500 trees were omitted as burn-in. To ensure convergence, the run was not ended until the average standard deviation of split frequencies reached <0.01 and the PSRF values were close to 1 for all parameters. To investigate the contribution of different genes, the nucleotide data matrix containing the 1st and 2nd codon positions, rRNA and tRNA sequences was subjected to a heuristic parsimony analysis (i.e. hsearch addseq = random nreps = 1000 swap = TBR multrees = yes start = stepwise) in PAUP* 4.0  and TreeRot.v3  was used to calculate the partitioned Bremer support (PBS) values [30, 31] of each gene partition on the tree nodes.
Results and discussion
Genome organization and base composition
As observed in the previously determined Hoplonemertea mitogenomes, both of the present mitogenomes also include 13 protein-coding genes, two rRNAs and 22 tRNAs genes, all encoded on the coding strand except for trnP and trnT (Figure 1 and Table 1). The gene orders are identical to previously published hoplonemertean mitogenomes without exceptions. There are several overlaps throughout the two mitogenomes, for example, the 8-bp overlaps between nad6 and cob (Table 1).
The nucleotide composition of the coding strand is biased toward T and A in these two mitogenomes, as is the case in most metazoan mitogenomes . The A + T content of the coding strands in G. parasita and N. tetraclitophila is 68.8% and 71.2% respectively, which falls within the range of the previously sequenced nemertean mitogenomes (from 64.7% in Lineus alborostratus to 75.7% in Cephalothrix sp.) (Table 2). The A + T biased composition is particularly pronounced at the third codon position of the protein-coding genes (75.4% and 82.5%, respectively). The coding strands bear several poly-T stretches with the longest one being 20 Ts in G. parasita and 33 Ts in N. tetraclitophila, which have proved to be detrimental to PCR amplification [33, 34]. Among lophotrochozoans, AT- and GC skews always show high inter- or intra-phylum variation, which might affect phylogenetic analyses . The nucleotide skewness for the coding strands of N. tetraclitophila (AT-skew = −0.41, GC-skew = 0.34) and G. parasita (AT-skew = −0.46, GC-skew = 0.28) is biased toward T and G. A similar trend has been observed in other Nemertea mitogenomes (Figure 2): the negative AT-skew ranges from −0.46 (G. parasita) to −0.27 (Cephalothrix sp., C. hognkongiensis and L. alborostratus) and the GC-skew is always positive varying from 0.18 (Cephalothrix sp.) to 0.44 (N. punctatula). It is noteworthy that the nucleotide skews of the mNCRs are different among species (Figure 2), which reflects the relatively higher variability of mNCR. The mNCR of G. parasita is biased toward T and C, which is contrary to other hoplonemerteans.
The canonical start codons ATG and GTG are used in most protein-coding genes of the G. parasita and N. tetraclitophila mitogenomes. An exceptional case is the nad3 gene of N. tetraclitophila, which was inferred to be initiated by the ATT codon (Table 1), and its length (343 bp) is shorter than that of other Monostilifera species (354 bp). Nonstandard initiation codons were also inferred in previously sequenced nemertean mitogenomes, e.g., the cox1 (TCT) of Cephalothrix sp. and C. hongkongiensis. The majority of the protein-coding genes appear to use the stop codons TAA or TAG, except that the nad5, nad3 and nad2 genes in N. tetraclitophila and the nad2 gene in G. parasita use a single T as the termination codon, most of which are adjacent to a protein-coding gene and occasionally a tRNA gene (Table 1). The incomplete termination codon T has been proposed to be converted into the complete stop codon TAA through polyadenylation during posttranscriptional mRNA processing . The overall length of protein-coding genes in the known nemertean mitogenomes varies from 11066 to 11268 bp. The protein-coding genes in seven Monostilifera mitogenomes are shorter than that in the other nemerteans (Figure 3A). The two present mitogenomes do not exhibit apparent length change compared to other hoplonemertean mitogenomes, unlike in some parasitic insects whose protein-coding gene sizes are significantly smaller than those of free-living ones .
The overall nucleotide composition of 13 protein-coding genes in G. parasita and N. tetraclitophila mitogenomes are AT biased (68.1% and 70.9%, respectively). For both species, the third codon position has a considerably higher AT content (75.4% and 82.5%, respectively) than the first and second codon positions and the lowest content of C (Table 2). According to the analysis of relative synonymous codon usage (RSCU), the two- and four-fold degenerate codons prefer the one ending with T (Additional file 2: Table S2), for example, GCT (2.811) is more frequently used than the other three codons (0.297-0.486) for Ala. Corresponding to the high percentage of T in both mitogenomes, the most frequently used codon is TTT (17.6% and 16.3%, respectively), and Phe is the most frequently used amino acid (19.2% and 16.8%, respectively) (Additional file 2: Table S2). The other preferred amino acids in both species are Leu, Val, Gly and Ser, all of which might be associated with transmembrane functions. Similar codon usage and amino acid composition patterns have been observed in previously sequenced Nemertea mitogenomes .
Ribosomal and transfer RNA genes
The ribosomal RNA genes (rrnL and rrnS) are located at the same location as in other nemertean mitogenomes, separated by trnV. The rrnL gene is 1,108 bp in G. parasita and 1,090 bp in N. tetraclitophila, and the A + T contents are 72.7% and 73.6%, respectively. The rrnS gene is 837 bp and 773 bp, and the A + T content is 70.9% and 72.8%, respectively (Table 2). At the 5′ end of rrnS gene in G. parasita, there is a region of 58 bp (TGTTTATTGGTATATTTTGATAAGTACTTTTAGTTTTATTCTATTTTTTTTCTTGTTT), which can neither be aligned with other nemertean rrnS sequences nor does it show any similarity with any remaining parts of the mitogenome, making the rrnS gene in G. parasita the longest among enoplan mitogenomes (Figure 3B). This insertion is also one major reason that G. parasita bears the largest mitogenome within Distromatonemertea. Except for rrnS of G. parasita, the rRNA genes of monostiferans are apparently shorter than that of other nemerteans (Figure 3B).
A + T contents in the tRNA genes is slightly lower than in the remainder of the mitogenomes. The anticodons of 22 tRNAs in both mitogenomes are the same as in other hoplonemerteans. All tRNA genes can be folded into conventional cloverleaf-like structures, except for trnS1(UCN) and trnS2(AGN) of G. parasita, and trnS2(AGN) of N. tetraclitophila. The structures of trnS2 of both species conform to the secondary structure achieved for known hoplonemertean mitogenomes, all lacking a DHU-arm which is replaced by a DHU-loop [38, 39]. trnS1 of G. parasita was inferred to be 59 bp, which makes it one of the shortest known tRNA genes of nemerteans. It has a 5-T DHU-loop instead of a DHU-arm. Uncanonical secondary structures of tRNA genes occur frequently during animal evolution .
There are a total of 265 bp and 250 bp non-coding nucleotides throughout the mitogenomes of G. parasita and N. tetraclitophila, accounting for 1.8% and 1.7% of the whole mitogenomes, respectively. The mNCRs are 120 bp and 137 bp, respectively, both located between trnW and trnS2. The A + T content (70.8% and 73.7%) of both mNCRs is slightly higher compared with the whole coding strands, but not as high as that of the third codon position (Table 2). Besides poly-T/C/G stretches, the two mNCRs have a similarity of 33%, which reflects a rapid evolutionary rate. Tandem repeats like those in Amphiporus formidabilis and Nipponnemertes punctatula are not detected. In both, N. tetraclitophila and G. parasita, the mNCRs have the potential to fold into hairpin-like structures at the 5′ end (not shown), which might be involved in the beginning of replication and transcription . The second longest mNCRs in the mitogenomes of N. tetraclitophila and G. parasita are both located between cox1 and trnW (33 bp and 28 bp, respectively), in agreement with other Monostilifera species .
The concatenated datasets for amino acid and nucleotide sequences of the 13 protein-coding genes (excluding the 3rd codon position) yielded 3,056 and 6,721 aligned sites, respectively. The third dataset (comprising 8,962 nucleotide sites) was constructed by adding informative rRNA and tRNA gene sites to the above nucleotide dataset, which can help avoid directional migration resulting from only using protein-coding genes . According to the Partitioned Bremer support (PBS) analysis , the rRNA and tRNA sequences contribute 17.7% and 11.7% (Table 3) of phylogenetic signal, respectively, making them promising for phylogenetic analysis.
Based on these three datasets, ML and BI analyses yielded identical tree topologies (Figure 4). All of them support the hypothesis that Hoplonemertea has a closer relationship with Heteronemertea than with Palaeonemertea, represented here by two Cephalothrix species that form the earliest divergent clade with high bootstrap values and posterior probabilities. As documented in previous studies [7, 44], Polystilifera (Nectonemertes cf. mirabilis) is the sister group to Monostilifera; Nipponnemertes is sister to the other monostiliferans which make up the group Distromatonemertea . The two present taxa, G. parasita and N. tetraclitophila, exhibit early divergent positions in the analyzed Distromatonemertea species. A recent analysis based on data of six genes also placed G. parasita in a basal Distromatonemertea clade containing mostly symbiotic and terrestrial species , whereas it was placed at a different position in the phylogenetic analysis of cox1 and 18S rRNA sequences . No similar species of the genus Nemertopsis have been studied in previous phylogenetic analyses. The position of the congeneric free-living species, Nemertopsis bivittata, was more or less different in previous analyses [8, 9] and seems to be different from the placement of N. tetraclitophila in the present study, which calls for further studies about the interrelationships within the genus Nemertopsis.
The complete mitochondrial genomes of Gononemertes parasita and Nemertopsis tetraclitophila, both of which possess some morphological characteristics adaptive to their lifestyle, are 14742 bp and 14597 bp, respectively. They are identical to the previously published mitogenomes of free-living hoplonemerteans in gene content and gene order, and have similar patterns in nucleotide richness and skewness. The length of whole genomes, as well as protein-coding genes and ribosomal RNA genes, is relatively conservative within Distromatonemertea and shorter (with the exception of the rrnS of G. parasita) than that of the other nemerteans. As in other hoplonemerteans, the coding strands of the present two mitogenomes bear some poly-T stretches; the tRNA genes usually exhibit cloverleaf-like structure except for trnS; the major non-coding regions exhibit AT-rich and hairpin-like structures that may be involved in transcription and replication. Some differences are found between the present mitogenomes and other hoplonemertean mitogenomes. For example, in G. parasita the mNCR is biased toward T and C (contrary to that in other hoplonemerteans) and the rrnS gene has a unique 58-bp insertion at 5′ end, and in N. tetraclitophila the nad3 gene starts with the ATT codon (ATG in other hoplonemerteans). However, we cannot conclude that these differences are related to their special lifestyle, because similar variations may also exist among free-living nemerteans and available mitogenomic data of nemerteans are stilled limited. Phylogenetic analyses show that both G. parasita and N. tetraclitophila are early divergent within the analyzed Distromatonemertea species.
- atp6 and atp8:
ATP synthase subunits 6 and 8
- cytb :
- cox1-3 Cytochrome c:
oxidase subunits I-III
- nad1-6 and nad4L:
NADH dehydrogenase subunits 1–6 and 4 L
- rrnL and rrnS:
The large and small subunits of ribosomal RNA
- trnX :
Transfer RNA molecules with the one-letter code of corresponding amino acid
Major non-coding region
Polymerase chain reaction
Kajihara H, Chernyshev AV, Sun S-C, Sundberg P, Crandall FB: Checklist of nemertean genera and species published between 1995–2007. Species Div. 2008, 13: 245-274.
Turbeville JM, Field KG, Raff RA: Phylogenetic position of phylum Nemertini, inferred from 18S rRNA sequences: molecular data as a test of morphological character homology. Mol Biol Evol. 1992, 9: 235-249.
Passamaneck Y, Halanych KM: Lophotrochozoan phylogeny assessed with LSU and SSU data: evidence of lophophorate polyphyly. Mol Phylogenet Evol. 2006, 40: 20-28. 10.1016/j.ympev.2006.02.001.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sørensen MV, Haddock SH, Schmidt-Rhaesa A, Okusu A, Kristensen RM, Wheeler WC, Martindale MQ, Giribet G: Broad phylogenomic sampling improves resolution of the animal tree of life. Nature. 2008, 452: 745-749. 10.1038/nature06614.
Podsiadlowski L, Braband A, Struck TH, von Dohren J, Bartolomaeus T: Phylogeny and mitochondrial gene order variation in Lophotrochozoa in the light of new mitogenomic data from Nemertea. BMC Genomics. 2009, 10: 364-10.1186/1471-2164-10-364.
Turbeville JM: An ultrastructural analysis of coelomogenesis in the hoplonemertine Prosorhochmus americanus and the polychaete Magelona sp. J Morphol. 1986, 187: 51-60. 10.1002/jmor.1051870105.
Sundberg P, Turbeville JM, Lindh S: Phylogenetic relationships among higher Nemertean (Nemertea) Taxa inferred from 18S rDNA sequences. Mol Phylogenet Evol. 2001, 20: 327-334. 10.1006/mpev.2001.0982.
Thollesson M, Norenburg JL: Ribbon worm relationships: a phylogeny of the phylum Nemertea. Proc Biol Sci. 2003, 270: 407-415. 10.1098/rspb.2002.2254.
Andrade SCS, Strand M, Schwartz M, Chen HX, Kajihara H, Von Döhren J, Sun SC, Junoy J, Thiele M, Norenburg JL, Turbeville JM, Giribet G, Sundberg P: Disentangling ribbon worm relationships: multi-locus analysis supports traditional classification of the phylum Nemertea. Cladistics. 2012, 28: 141-159. 10.1111/j.1096-0031.2011.00376.x.
Strand M, Samuelsson H, Sundberg P: Nationalnyckeln till Sveriges flora och fauna. Stjärnmaskar − slemmaskar. Sipuncula − Nemertea. 2010, Uppsala: ArtDatabanken, SLU
Gibson R: The macrobenthic nemertean fauna of Hong Kong. Proceedings of the Second International Marine Biological Workshop: the Marine Flora and Fauna of Hong Kong and Southern China, 1986. Volume 1. Edited by: Morton B. 1990, Hong Kong: Hong Kong University Press, 33-212.
Bergendal D: Über ein Paar sehr eigenthümliche nordische Nemertinen. Zoologischer Anzeiger. 1900, 23: 313-328.
Roe P: Ecological implications of the reproductive biology of symbiotic nemerteans. Hydrobiologia. 1988, 156: 13-22. 10.1007/BF00027973.
Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: 955-964. 10.1093/nar/25.5.955.
Denman RB: Using RNAFOLD to predict the activity of small catalytic RNAs. Biotechniques. 1993, 15: 1090-1095.
Grant JR, Stothard P: The CGView server: a comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36: W181-W184. 10.1093/nar/gkn179.
Xia X: DAMBE5: a comprehensive software package for data analysis in molecular biology and evolution. Mol Biol Evol. 2013, 30: 1720-1728. 10.1093/molbev/mst064.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Golombek A, Tobergte S, Nesnidal MP, Purschke G, Struck TH: Mitochondrial genomes to the rescue - diurodrilidae in the myzostomid trap. Mol Phylogenet Evol. 2013, 68: 312-326. 10.1016/j.ympev.2013.03.026.
Stöger I, Schrodl M: Mitogenomics does not resolve deep molluscan relationships (yet?). Mol Phylogenet Evol. 2013, 69: 376-392. 10.1016/j.ympev.2012.11.017.
Talavera G, Castresana J: Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007, 56: 564-577. 10.1080/10635150701472164.
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
Nylander JAA: MrModeltest v2. Program distributed by the author. 2004, Evolutionary Biology Centre, Uppsala University, [http://www.abc.se/]
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O: New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010, 59: 307-321. 10.1093/sysbio/syq010.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.
Swofford DL: PAUP* – Phylogenetic Analysis Using Parsimony (*and Other Methods). Ver. 4.0. [Computer software and manual]. Sunderland, MA: Sinauer Associates. Book PAUP* – Phylogenetic Analysis Using Parsimony (*and Other Methods). Ver. 4.0. [Computer software and manual]. 1999, Sunderland, MA: Sinauer Associates
TreeRot version 3. [http://people.bu.edu/msoren/TreeRot.html]
Bremer K: The limits of amino acid sequence data in angiosperm phylogenetic reconstruction. Evolution. 1988, 42: 795-803. 10.2307/2408870.
Bremer K: Branch support and tree stability. Cladistics. 1994, 10: 295-304. 10.1111/j.1096-0031.1994.tb00179.x.
Hassanin A, Leger N, Deutsch J: Evidence for multiple reversals of asymmetric mutational constraints during the evolution of the mitochondrial genome of metazoa, and consequences for phylogenetic inferences. Syst Biol. 2005, 54: 277-298. 10.1080/10635150590947843.
Riepsamen AH, Gibson T, Rowe J, Chitwood DJ, Subbotin SA, Dowton M: Poly(T) variation in heteroderid nematode mitochondrial genomes is predominantly an artefact of amplification. J Mol Evol. 2011, 72: 182-192. 10.1007/s00239-010-9414-3.
Gibson T, Farrugia D, Barrett J, Chitwood DJ, Rowe J, Subbotin S, Dowton M: The mitochondrial genome of the soybean cyst nematode, Heterodera glycines. Genome. 2011, 54: 565-574. 10.1139/g11-024.
Nesnidal MP, Helmkampf M, Bruchhaus I, Hausdorf B: The complete mitochondrial genome of Flustra foliacea (Ectoprocta, Cheilostomata) - compositional bias affects phylogenetic analyses of lophotrochozoan relationships. BMC Genomics. 2011, 12: 572-10.1186/1471-2164-12-572.
Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290: 470-474. 10.1038/290470a0.
McMahon DP, Hayward A, Kathirithamby J: The mitochondrial genome of the ‘twisted-wing parasite’ Mengenilla australiensis (Insecta, Strepsiptera): a comparative study. BMC Genomics. 2009, 10: 603-10.1186/1471-2164-10-603.
Chen HX, Sun SC, Sundberg P, Ren WC, Norenburg JL: 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. BMC Genomics. 2012, 13: 139-10.1186/1471-2164-13-139.
Chen HX, Sundberg P, Wu HY, Sun SC: The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea). Mol Biol Rep. 2011, 38: 4509-4525. 10.1007/s11033-010-0582-4.
Gasser RB, Jabbar A, Mohandas N, Hoglund J, Hall RS, Littlewood DT, Jex AR: Assessment of the genetic relationship between Dictyocaulus species from Bos taurus and Cervus elaphus using complete mitochondrial genomic datasets. Parasit Vectors. 2012, 5: 241-10.1186/1756-3305-5-241.
Sun WY, Sun SC: A description of the complete mitochondrial genomes of Amphiporus formidabilis. Mol Biol Rep, in press: Prosadenoporus spectaculum and Nipponnemertes punctatula (Nemertea: Hoplonemertea: Monostilifera), doi:10.1007/s11033-014-3438-5
Arnason E, Rand DM: Heteroplasmy of short tandem repeats in mitochondrial DNA of Atlantic cod, Gadus morhua. Genetics. 1992, 132: 211-220.
Kumazawa Y, Nishida M: Sequence evolution of mitochondrial tRNA genes and deep-branch animal phylogenetics. J Mol Evol. 1993, 37: 380-398.
Sundberg P, Chernyshev AV, Kajihara H, Kanneby T, Strand M: Character-matrix based descriptions of two new nemertean (Nemertea) species. Zool J Linn Soc. 2009, 157: 264-294. 10.1111/j.1096-3642.2008.00514.x.
Junoy J, Andrade SCS, Giribet G: Phylogenetic placement of a new hoplonemertean species commensal on ascidians. Invertebr Syst. 2010, 24: 616-629.
This work is supported by the National Natural Science Foundation of China (31172046; 30970333). We are grateful to Mr. Hai-Lin Shen for assistance in collecting specimens.
The authors declare that they have no competing interests.
S-CS conceived and designed the study. W-YS and D-LX performed the experiments. W-YS, D-LX and WS analyzed the data. MS, PS and S-CS collected and identified specimens. WY-S drafted the manuscript. S-C S, H-XC and PS revised the manuscript. All authors read and approved the final manuscript.