- Research article
- Open Access
The complete mitochondrial genome of Scutopus ventrolineatus (Mollusca: Chaetodermomorpha) supports the Aculifera hypothesis
BMC Evolutionary Biology volume 14, Article number: 197 (2014)
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.
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.
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.
Mollusks are the second largest animal phylum with more than 100000 described extant species and are grouped into eight classes: Solenogastres or Neomeniomorpha, Caudofoveata or Chaetodermomorpha, Polyplacophora, Monoplacophora, Bivalvia, Gastropoda, Cephalopoda, and Scaphopoda -. Morphology-based classifications consider Neomeniomorpha and Chaetodermomorpha as the earliest branching lineages within mollusks because they lack many typical features of mollusks, among which the most conspicuous is the shell . Whether they form a monophyletic group named Aplacophora , or a paraphyletic grade - is still a matter of debate . Some authors - place Polyplacophora (chitons) as sister group of Conchifera (Monoplacophora, Bivalvia, Gastropoda, Cephalopoda, and Scaphopoda), forming the clade Testaria. Alternatively, the Aculifera hypothesis proposes a sister group relationship of Polyplacophora and Aplacophora, suggesting that the aplacophoran morphology was secondarily modified from a chiton-like ancestor ,,,. Phylogenetic relationships within Conchifera are also far from settled due to the highly derived morphologies within each class-level grade, which hinder the discovery of morphological synapomorphies across lineages . The traditional morphology-based hypothesis groups together Bivalvia + Scaphopoda (Diasoma or Loboconcha) and Cephalopoda + Gastropoda (Cyrtosoma or Visceroconcha) ,.
Earlier molecular studies based on partial sequences of one or few genes revealed important phylogenetic inference biases, and failed to recover the monophyly of mollusks and/or of several main lineages within the group . The first studies analyzing the relative phylogenetic position of Monoplacophora , also rendered a surprising result in recovering the group as closely related to Polyplacophora, forming the taxon Serialia, which is in disagreement with most morphological evidence . More recently, attempts to reconstruct the phylogeny of mollusks were based on concatenated matrices spanning many genes. A study  based on 79 ribosomal protein genes recovered the monophyly of the phylum and of all five mollusk classes included in the analyses (Neomeniomorpha, Monoplacophora, and Scaphopoda were missing). However, recovered interclass relationships, although highly supported, were rather unconventional, with Bivalvia and Gastropoda being sister group to Polyplacophora, and this clade being sister group to Cephalopoda and Chaetodermomorpha . The latest phylogenetic studies based on seven housekeeping genes , and on genomic-scale data sets spanning 308  and 1185  genes, respectively, recovered monophyletic Mollusca, Aplacophora, and Aculifera, thus rejecting the Testaria hypothesis. One of these studies included Monoplacophora, which was placed within Conchifera, thus providing no support for the Serialia hypothesis . Interestingly, phylogenomic studies , arrived at highly supported but contrasting conchiferan interclass relationships. While both studies support a basal position of Cephalopoda (+Monoplacophora in ), one favors a clade composed of Gastropoda and Bivalvia (the so-called Pleistomollusca ) whereas the other groups together Gastropoda and Scaphopoda .
During the last decade, complete mitochondrial (mt) genomes have become a standard for phylogenetic reconstruction of animal relationships . Although the number of completely sequenced mollusk mt genomes has increased considerably in the last few years, the majority belong to the most common and economically important mollusk classes i.e., Cephalopoda , Bivalvia , and Gastropoda ,. In addition, there are reported three Polyplacophora, Katharina tunicata and two Sypharochiton species , and two Scaphopoda, Siphonodentalium lobatum and Graptacme eborea mitogenomes (but note that the mt genome of a Chaetodermomorpha, Chaetoderma nitidulum is available at NCBI, although unpublished). The only published study  that has applied thus far whole mt genome data for reconstructing phylogenetic relationships of mollusks failed to recover the monophyly of Mollusca and of many mollusk classes due to long branch attraction (LBA) artifacts. The authors concluded that representatives from all mollusks classes and a denser taxon sampling of most diverse lineages could render a more resolved mollusk interclass phylogeny, and that mt gene order data could become a promising source of phylogenetic information .
In this paper, we present the complete mitochondrial genome of Scutopus ventrolineatus, a representative of the supposedly early-branching Limifossoridae within Chaetodermomorpha. We performed comparative genomic analyses with other available mollusk mt genomes with the specific aim of addressing the evolution of gene order arrangements among the main lineages within the phylum. In addition, we used generated mt sequence data to infer phylogenetic relationships of mollusks, and in particular to test the validity of traditional morphology-based hypotheses that place Chaetodermomorpha at the base of the mollusk tree. Finally, we concatenated mt genomes with publicly available nuclear sequence data in trying to maximize statistical support of the reconstructed mollusk phylogeny.
Results and discussion
Mitochondrial genome organization and structural features
The complete mt genome of Scutopus ventrolineatus was assembled as a 14662 bp circular molecule. Like most metazoan mt genomes, it encodes for 13 protein-coding, 22 tRNA and 2 rRNA genes (Figure 1). The major strand encodes 20 out of the 37 genes (trnF, nad5, trnH, nad4, nad4L, trnS(UCN), cob, nad6, nad1, trnL(UUR), trnL(CUN), rrnS, trnM, trnC, trnQ, trnY, rrnL, trnV, trnG, trnW). Most protein-coding genes start with the codon ATG with the exception of atp8, which begins with GTG. Several genes show complete stop codons, either TAA (as in cox1, atp8, atp6, nad4l, cob, and nad6) or TAG (as in cox2 and nad1). The remaining genes finish with either TA (nad5) or a single T (cox3, nad2, nad4), which presumably become functional stop codons by subsequent polyadenylation of the transcribed messenger RNAs . Four genes overlap with contiguous genes: cox2 with trnD; atp8 with atp6; nad4 with trnH; and nad4L with nad4. The largest noncoding region has 47 bp, and it is located between trnW and trnE genes.
The gene order arrangement of the mt genome of S. ventrolineatus was compared with the mt genome organization in other mollusks. The unpublished mt genome of the chaetodermomorph C. nitidulum is the closest to compare, but it is 40% longer (21008 bp), and shows a duplication of the cox2 gene, as well as large non-coding regions (see Additional file 1: Figure S1). In addition, several tRNA genes (trnV, trnF, trnG, trnW, trnA, trnR, trnI, trnS (UCN)) are reordered and two (trnV and trnS (UCN)) are encoded on opposite strands when compared to the mt genome of S. ventrolineatus (see Additional file 1: Figure S1). Since the mt genome organization of S. ventrolineatus conforms to the mollusk consensus gene order (see below), it would be important to determine the origin of the highly divergent mt genome of C. nitidulum, and discard putative sequencing or assembly errors.
Compared with other mollusks, the mt gene order of S. ventrolineatus is very similar to those of Polyplacophora, Cephalopoda and basal lineages of Gastropoda (Vetigastropoda and Neritimorpha; but not of Patellogastropoda, ) (Figure 2). There are no differences in the relative position of protein-coding genes, but the rrnL gene is translocated in Scutopus mt genome with respect to the above-mentioned lineages (Figure 2). Other differences are found in the relative position of several tRNA genes; in particular trnD, trnE, and trnN are highly rearranged among mollusk lineages (Figure 2). The representative of an early-diverging lineage of Bivalvia (Solemya velum; Solemyoidea) also retains this gene order although a large inversion has occurred affecting a stretch including atp8-atp6-nad5-nad4-nad4L (Figure 2). Hence, our results allow us to propose a consensus ancestral mollusk gene order for mt protein-coding and rRNA genes (Figure 2). Only few translocation and inversion events are required to transform this consensus mollusk gene order into the genome organizations reported for other Spiralia, particularly Phoronida, Brachiopoda, Nemertea, and Entoprocta (Figure 2). The conversion to the genome organization of Annelida requires postulating additional translocations and inversions . In Scaphopoda, the mt gene order is rather different to that of other mollusks, both in Graptacme (Figure 2) and Siphonodentalium. In Graptacme, only the nad6-trnP-nad1 and the nad5-nad4-nad4L arrangements are conserved, whereas gene pairs that are normally conserved in other metazoan mt genomes (e.g. cox1-cox2 and atp8-atp6) are not found (Figure 2). New complete mt genome sequences of more (early-branching) lineages within Scaphopoda are needed in order to determine whether Scaphopoda ancestrally conformed or not to the consensus ancestral mt gene order of mollusks. Within Gastropoda and Bivalvia, highly derived mt genome organizations can be found due to extensive gene rearrangements ,.
Phylogenetic relationships within Mollusca
Despite its popularity, phylogenetic inference based on complete mt genomes has important analytical challenges. In fact, reconstructing the monophyly and phylogenetic relationships of mollusks using solely mt genome data proved to be difficult. A first issue was finding the appropriate outgroup. Most recent phylogenies of spiralians place mollusks within the clade Trochozoa, which also includes annelids, nemerteans, brachiopods, and phoronids . Among available mt genomes of spiralians, we selected representatives from one non-trochozoan spiralian (Entoprocta) and all trochozoan phyla. Another important problem was LBA due to among-lineage rate heterogeneity. Several ingroup lineages, notably Bivalvia, Scaphopoda, and Patellogastropoda (within Gastropoda), exhibited long branches that produced considerable inference biases (not shown), and they had to be excluded in further phylogenetic analyses. At least in several bivalves (families Donacidae, Hyriidae, Margaritiferidae, Mytilidae, Solenidae, Unionidae, and Veneridae), the occurrence of doubly uniparental inheritance (DUI) of the mt genome  has been proposed to cause higher substitution rates and thus faster saturation of the phylogenetic signal, hampering their use in reconstructing deep phylogenies. In addition, Neomeniomorpha and Monoplacophora could not be incorporated into the phylogenetic analyses because complete mt genomes are not available for these lineages thus far. Finally, due to the relatively high substitution rates of mtDNA, saturation is another major issue affecting the reconstruction of deep nodes based on mt markers. Here, we opted to analyze mt data at the amino acid level to reduce the effect of saturation.
The mt data set included 4870 and 2728 positions before and after sites of ambiguous positional homology were discarded. The fact that 44% of the initial length of the alignments was excluded by GBlocks suggests that mt genome amino acid sequence data have complex evolutionary patterns at higher levels of divergence, which would produce phylogenetic reconstruction biases if conflicting positions were not removed. ML (-lnL = 57186.85) and BI (arithmetic mean of the two runs, -lnL = 67264.29) arrived at relatively different trees (Figure 3). In ML, the monophyly of Mollusca is recovered with low bootstrap support whereas in BI, Brachiopoda + Annelida are recovered with maximal posterior probability support as sister group of Chaetodermomorpha + Polyplacophora, rendering Mollusca non-monophyletic. Within mollusks, the clades Aculifera and Conchifera received maximal posterior probability support in BI and low (<50%) bootstrap support in the ML analysis. The phylogenetic analyses failed to recover the monophyly of gastropods due to a LBA artifact (highly supported) between Roboastra + Micromelo (Heterobranchia) and Scaphopoda, both in ML and BI. This clade is the sister group of remaining gastropods with low bootstrap support in ML or of cephalopods with maximal posterior probability support in BI (Figure 3).
Recently, two phylogenies of Mollusca have been reconstructed based on large concatenated data sets of hundreds or thousands of nuclear gene fragments ,. Here, we selected nuclear genes from those matrices that minimized missing data and maximized taxon coverage to resolve the question at hand. The partial deduced amino acid sequences of 13 nuclear ribosomal protein genes (see also ) were aligned for 5 outgroup taxa and main mollusk lineages (19 species). The concatenation of ribosomal proteins yielded an alignment of 2362 positions after removal of ambiguous sites. Phylogenetic analyses under ML (-lnL = 33338.02) and BI (arithmetic mean of the two runs, -lnL = 36866.21) recovered similar topologies that only differed on the relative position of the neritimorph Theodoxus (sister group of Caenogastropoda + Heterobranchia in ML and sister group of Vetigastropoda + Patellogastropoda in BI) and on the closest sister group to mollusks (a nemertean in ML and the annelids in BI) (Figure 4). In both ML and BI, Aplacophora and Aculifera were recovered as monophyletic (as in ,), although with low support (except Aculifera in BI that has a posterior probability of 0.99). In contrast, both phylogenetic analyses failed to recover the monophyly of Conchifera because Cephalopoda was placed as sister group of Aculifera (with low support in ML but with maximal Bayesian posterior probability) (Figure 4). In both ML and BI, Bilvavia was recovered as sister group of Scaphopoda + Gastropoda in agreement with . Interestingly, Scaphopoda and Lottia showed relatively long branches, which clearly classify them as rogue taxa, showing extremely fast evolutionary rates for both mt and nuclear genes; ,. The reconstructed trees based on the nuclear data set improved previous results obtained from ribosomal nuclear protein  and the 18S rRNA gene , which recovered Cephalopoda as the sister group of Chaetodermomorpha or Neomeniomorpha, respectively. However, our reconstructed trees resemble those based on housekeeping genes , that also placed Cephalopoda a sister group to Aculifera. The odd placement of Cephalopoda here and in the above-mentioned studies ,, contradicts the general agreement placing cephalopods within Conchifera . This might be a tree reconstruction artifact produced by LBA biases generated by the inclusion of relatively distantly related spiralian phyla in the outgroup.
Although with a different lineage sampling and species representation, the mt and nuclear data sets rendered rather congruent trees with several nodes in common. Therefore, we concatenated the nuclear and mt data sets. The phylogenetic analyses of the combined (mt and nuclear) data set (5032 sites after removal of ambiguous positions) under ML (-lnL = 101714.93) and BI (arithmetic mean of the two runs, -lnL = 116220.97) arrived at very similar topologies that only differed in the relative position of Bivalvia and Scaphopoda (Figure 5). The monophyly of Aplacophora, Aculifera and Conchifera were recovered both by ML and BI (Figure 5). The corresponding nodes received relatively high support except Conchifera in ML and Aplacophora in BI (Figure 5). In ML and BI, Cephalopoda and Cephalopoda + Scaphopoda were recovered as sister group of the remaining analyzed conchiferan classes, respectively (Figure 5). In ML, Gastropoda was recovered as sister group of Bivalvia + Scaphopoda, whereas in BI Bivalvia was placed nested within Gastropoda, rendering the latter non-monophyletic (Figure 5). These results agree with most recent phylogenomic studies , in the basal position of Cephalopoda within Conchifera, and favor either a close relationship of Gastropoda and Bivalvia to the exclusion of Scaphopoda  or the classical Diasoma hypothesis uniting Bivalvia and Scaphopoda . It is noteworthy that the reconstructed phylogeny based on the combined data set differs from the nuclear-based tree indicating that addition of mt data has a significant (and distinct) contribution to the overall phylogenetic inference. Furthermore, previously encountered LBA artifacts when mt data was analyzed alone (related to Scaphopoda, Bivalvia and Lottia) were ameliorated in the combined analysis.
Reconstructing the monophyly and internal phylogenetic relationships of Mollusca based on molecular data has been challenging over the years. Earlier studies were mostly based on partial gene sequences with insufficient informative characters to reconstruct robust and resolved trees . Recent phylogenomic studies , based on concatenated nuclear genes have produced well-resolved trees. However, comparison of these trees showed some contradicting but highly supported nodes. Another alternative is to investigate phylogenetic information contained in complete mt genomes. A first attempt to use whole mt genome data for clarifying the relative phylogenetic position of mollusks within bilaterians was not satisfying . Our phylogenetic analyses reveal that the use of mt genome data for reconstructing the internal phylogenetic relationships of mollusks is flawed by the existence of high heterogeneity of evolutionary rates among lineages. In particular, the mt genomes of Bivalvia and Lottia introduce detrimental LBA biases. After removing taxa that exhibit long branches in the phylogenetic analyses, mt data are capable of recovering the monophyly of each Aculifera, and Conchifera, in agreement with nuclear data (which additionally recover the monophyly of Aplacophora). Moreover, our analyses indicate that the phylogenetic performance of mt and nuclear data improves when both are combined. At present, the complete mt genomes of several important mollusk lineages (namely Neomeniomorpha and Monoplacophora) are still missing, and the possibility of finding mt genomes with lower substitution rates in Bilvalvia and Patellogastropoda needs to be further explored. It is foreseeable that the addition of these additional mt genomes will improve phylogenetic analyses. In parallel, gapped regions in nuclear genomic data sets will be increasingly reduced allowing in combination with mt data the reconstruction of a robust tree of Mollusca. Such a robust phylogenetic hypothesis has been long-needed as the framework for evolutionary comparative studies within this highly diversified metazoan phylum.
DNA extraction, PCR amplification, cloning and sequencing
Several specimens of Scutopus ventrolineatus Salvini-Plawen, 1968 (Mollusca; Chaetodermomorpha) were collected in March 2010 close to Bergen, in the Norwegian west coast. Total genomic DNA of a single specimen was isolated following standard phenol-chloroform extraction procedures . Two fragments corresponding to partial mt cox1 cox2 genes were PCR -amplified using appropriate universal primer pairs: LCOI 1490, HCO 2198  and cox2F1, cox2R1 , respectively. PCR reactions contained 2.5 μlof 10x Ex Taq Buffer, 2 μlof dNTP Mixture (2.5 mM each), 1.5 μlof each primer, 0.5 μlof template DNA, 0.16 μlEx Taq Hot Start DNA polymerase (5 units/μl; TaKaRa Bio Inc., Otsu, Shiga, Japan), and sterilized distilled water up to 25 μl. The following temperature profile was used: an initial denaturing step at 98°C for 2 min; 35 cycles of denaturing at 98°C for 10 s, annealing at 50°C for 30 s, and extending at 72°C for 1 min; and a final extending step at 72°C for 1 min. PCR products were purified by ethanol precipitation, and sequenced in an automated DNA sequencer ABI PRISM 3700 using the BigDye Terminator v3.1 cycle-sequencing kit (Applied Biosystems; Foster City, CA, USA), and PCR primers.
Newly determined partial sequences of mt cox1 and cox2 were used to design two pairs of specific primers for long range PCR amplification (SVCOX1.F: 5'-TTT TTG ACC CTG CTG GAG GTG GAG AC-3'; SVCOX1.R: 5'-AGA GGG GGG TAT ACA GTC CAC CCA GTC-3'; SVCOX2.F: 5'-TCC CAG CAT TGG GAG TAA AAG CCG AC-3'; SVCOX2.R: 5'-CTC CGC AGA TTT CTG AAC ATT GAC CA-3'). The full mt genome was amplified in two overlapping fragments of 13003 bp and 1771 bp, respectively, using the TaKaRa LA-PCR kit (TaKaRa Bio Inc., Otsu, Shiga, Japan). PCR amplifications were carried out in 50 μlreactions containing 5 μlof 10x LA Buffer II (Mg+2 plus), 8 μlof dNTP Mixture (2.5 mM each), 1 μlof each primer, 0.1 μlof template DNA and 0.5 μlTaq DNA polymerase (5 units/μl). The following temperature profile was used: an initial denaturing step at 94°C for 1 min; 45 cycles of denaturing at 98°C for 10 s, annealing at 57°C for 30 s, and extending at 68°C for 1 min per Kb; and a final extending step at 68°C for 10 min. Sequencing of the two long PCR fragments was achieved with the shotgun technique using the TOPO-Shotgun subcloning Kit (Invitrogen; Life Technologies, Paisley, UK). Random clone libraries were constructed from purified PCR products by shearing them into fragments of 1–3 Kb in size, by repairing fragment ends to form blunt-ends, and by cloning blunt-ended fragments into pCR 4Blunt-TOPO vectors. Clones were sequenced in an automated DNA sequencer ABI PRISM 3700 using the BigDye Terminator v3.1 cycle-sequencing kit (Applied Biosystems; Foster City, CA, USA), and M13 forward and reverse universal primers.
Genome assembly and annotation
The complete mt genome was assembled into a single contig from the shotgun clone sequences using Sequencher v. 5.0.1 (Gene Codes Co.; Ann Arbor, MI, USA). The mt genome was annotated using the MITOS  and DOGMA  webservers. Briefly, protein-coding genes were annotated by identification of their open reading frames and similarity searches against other reported mollusk mt genomes. Ribosomal RNA genes were identified by sequence comparison with other reported mollusk mt genomes, and assumed to extend to the boundaries of adjacent genes. Transfer RNA genes were identified using tRNAscan-SE v. 1.21  and ARWEN v. 1.2 , which can infer cloverleaf secondary structures of the corresponding gene products. The complete mt genome sequence reported in this paper has been deposited at NCBI GenBank under accession number KC757645.
The gene order of the mt genome of S. ventrolineatus was compared to the following mollusk mt genomes: Chaetoderma nitidulum (Chaetodermomorpha; Dreyer and Steiner, unpublished), Katharina tunicata (Polyplacophora; ), Graptacme eborea (Scaphopoda; ), Octopus vulgaris (Cephalopoda; ), Solemya velum (Bivalvia; ), and Haliotis rubra (Gastropoda; ). The following spiralian mt genomes were also included in the gene order comparisons: Paranemertes cf. peregrina (Nemertea; ), Terebratulina retusa (Brachiopoda; ), Loxosomella aloxiata (Entoprocta; ), and Phoronis architecta (Phoronida; ).
Data sets and sequence alignment
Amino acid sequences derived from the 13 mt protein-coding genes were used to assemble the mt data set, which included 14 representatives of the main extant mollusk lineages and 6 species representing several metazoan phyla other than Mollusca (see Additional file 1: Table S1). Similarly, a nuclear data set was constructed with the deduced amino acid sequences of the genes coding for 40S ribosomal proteins S8 and S15, and 60S ribosomal proteins L3, L4, L5, L6, L7, L8, L10a, L16_L10, L17, L18a, and L32. These genes were extracted from Kocot et al. (2011) and were selected among available genes because they minimized missing data. This nuclear data set included five species representing several metazoan phyla other than Mollusca and 19 representatives of the main extant mollusk lineages (see Additional file 1: Table S2). A third data set was constructed combining the mt and nuclear data sets. The three data sets were designed to test specifically the monophylies of Mollusca and Aculifera. Since our phylogenetic analyses were focused on the order level or above, and in order to maximize the completeness of the nuclear and combined data sets, sequences from related species/genera/families/superfamilies (for which there is strong evidence for their monophyly) were merged in some instances (see Additional file 1: Table S2; ). Note that the original study from which nuclear data was extracted  already merged closely related species to maximize gene coverage. The species showing the shortest branches were selected as representatives of the different lineages in the three data sets.
Deduced amino acid sequences of the different mt and nuclear protein-coding genes were downloaded from GenBank and aligned separately using MAFFT v. 7  with default settings. Ambiguously aligned positions were removed using Gblocks, v. 0.19b  with default settings.
Alignment format conversions were performed using the ALTER webserver . For the three analyzed data sets, best-fit partition schemes and models of amino acid replacement were identified using the Akaike information criterion (AIC; ) as implemented in PartitionFinderProtein . For the mt data set, we tested the following a priori partition schemes: (1) all genes combined; (2) genes by functional group (atp, cox, nad, cob); (3) all genes separately except atp8/atp6 and nad4L/nad4, and (4) all genes independently. For the nuclear data set we tested (1) all genes combined, (2) by functional group (40S and 60S genes), and (3) all independent. For the combined data set, we tested all above-mentioned partition schemes. The AIC favored independent gene partitions in the nuclear and the combined data sets, whereas the best partition scheme for the mt data sets was that with all genes analyzed separately except for atp8/atp6 and nad4L/nad4. The resulting best-fit models for each partition are shown in Additional file 1: Table S3.
Phylogenetic relationships were inferred using maximum likelihood (ML) and Bayesian inference (BI). ML analyses were conducted with RAxML v. 7.0.4  using the rapid hill-climbing algorithm. For BI, we used MrBayes v. 3.1.2  running two independent analyses, each consisting in four simultaneous MCMC (Markov chain Monte Carlo) for 10 million generations, sampling every 1,000 generations, and discarding the first 25% generations as burnin (as judged by plots of ML scores and low SD of split frequencies) to prevent sampling before reaching stationarity of Markov chains. Support for internal branches was evaluated by non-parametric bootstrapping  with 1,000 replicates (ML) and by posterior probabilities (BI).
- atp6-8 :
ATP synthase subunits 6-8
- cob :
- cox1-3 :
Cytochrome c oxidase subunits 1-3
Long branch attraction
- nad1-6 4 L :
NADH dehydrogenase subunits 1-6, 4 L
Polymerase chain reaction
- rrnL :
Large ribosomal RNA
- rrnS :
Small ribosomal RNA
- trnX :
Transfer RNA for amino acid X (denoted by the one-letter IUPAC symbol)
Ponder WF, Lindberg DR: Phylogeny and Evolution of the Mollusca. 2008, Univ. California Press, California
Beesley PL, Ross GJB, Wells A: Mollusca: The Southern Synthesis. 1998, CSIRO Publishing, Melbourne
Haszprunar G, Wanninger A: Molluscs. Curr Biol. 2012, 22 (13): R510-R514. 10.1016/j.cub.2012.05.039.
Hyman L: The Invertebrates, vol. 6 (Mollusca I). 1967, McGraw Hill, Boston, MA
Ivanov DL: Origin of Aculifera and problems of monophyly of higher taxa in molluscs. Origin and Evolutionary Radiation of the Mollusca. Edited by: Taylor JD. 1996, Oxford University Press, Oxford, 56-65.
Scheltema A: Aplacophora as progenetic aculiferans and the coelomate origin of mollusks as the sister taxon of Sipuncula. Biol Bull. 1993, 184: 57-78. 10.2307/1542380.
Faller S, Rothe B, Todt C, Schmidt-Rhaesa A, Loesel R: Comparative neuroanatomy of Caudofoveata, Solenogastres, Polyplacophora, and Scaphopoda (Mollusca) and its phylogenetic implications. Zoomorphology. 2012, 131 (2): 149-170. 10.1007/s00435-012-0150-7.
Haszprunar G: Is the Aplacophora monophyletic? a cladistic point of view. Am Malacol Bull. 2000, 15 (2): 115-130.
Haszprunar G, Schander C, Halanych KM: Relationships of the higher molluscan taxa. Towards a Phylogeny of Mollusca. Edited by: Ponder WF, Lindberg DR. 2008, Univ. Calif. Press, Berkeley, 19-32.
Salvini-Plawen L: On the phylogenetic significance of the aplacophoran Mollusca. Iberus. 2003, 21 (1): 67-97.
Todt C: Aplacophoran mollusks—still obscure and difficult?. Am Malacol Bull. 2013, 31 (1): 181-187. 10.4003/006.031.0110.
Salvini-Plawen L: Significance of the placophora for molluscan phylogeny. Venus. 2006, 65: 1-17.
Salvini-Plawen L, Steiner G: Synapomorphies and plesiomorphies in higher classification of Mollusca. Origin and Evolutionary Radiation of the Mollusca. Edited by: Taylor J. 1996, Oxford University Press, Oxford, 29-52.
Waller TR: Origin of the molluscan class Bivalvia and a phylogeny of major groups. Bivalves: An Eon of Evolution. Edited by: Johnston PA, Haggart JW. 1998, University of Calgary Press, Calgary, 1-45.
Sutton MD, Briggs DEG, Siveter DJ, Siveter DJ, Sigwart JD: A Silurian armoured aplacophoran and implications for molluscan phylogeny. Nature. 2012, 490 (7418): 94-97. 10.1038/nature11328.
Scherholz M, Redl E, Wollesen T, Todt C, Wanninger A: Aplacophoran mollusks evolved from ancestors with polyplacophoran-like features. Curr Biol. 2013, 23 (21): 2130-2134. 10.1016/j.cub.2013.08.056.
Kocot KM: Recent advances and unanswered questions in deep molluscan phylogenetics. Am Malacol Bull. 2013, 31 (1): 195-208. 10.4003/006.031.0112.
Passamaneck YJ, Schander C, Halanych KM: Investigation of molluscan phylogeny using large-subunit and small-subunit nuclear rRNA sequences. Mol Phylogenet Evol. 2004, 32 (1): 25-38. 10.1016/j.ympev.2003.12.016.
Wilson N, Rouse G, Giribet G: Assessing the molluscan hypothesis Serialia (Monoplacophora + Polyplacophora) using novel molecular data. Mol Phylogenet Evol. 2010, 54 (1): 187-193. 10.1016/j.ympev.2009.07.028.
Giribet G, Okusu A, Lindgren A, Huff S, Schrodl M, Nishiguchi M: Evidence for a clade composed of molluscs with serially repeated structures: Monoplacophorans are related to chitons. Proc Natl Acad Sci U S A. 2006, 103 (20): 7723-7728. 10.1073/pnas.0602578103.
Wägele J, Letsch H, Klussmann-Kolb A, Mayer C, Misof B, Wägele H: Phylogenetic support values are not necessarily informative: the case of the Serialia hypothesis (a mollusk phylogeny). Front Zool. 2009, 6 (1): 1-15. 10.1186/1742-9994-6-12.
Meyer A, Witek A, Lieb B: Selecting ribosomal protein genes for invertebrate phylogenetic inferences: how many genes to resolve the Mollusca?. Methods Ecol Evol. 2011, 2 (1): 34-42. 10.1111/j.2041-210X.2010.00052.x.
Vinther J, Sperling EA, Briggs DEG, Peterson KJ: A molecular palaeobiological hypothesis for the origin of aplacophoran molluscs and their derivation from chiton-like ancestors. Proc Roy Soc B Biol Sci. 2012, 279 (1732): 1259-1268. 10.1098/rspb.2011.1773.
Kocot KM, Cannon JT, Todt C, Citarella MR, Kohn AB, Meyer A, Santos SR, Schander C, Moroz LL, Lieb B, Halanych KM: Phylogenomics reveals deep molluscan relationships. Nature. 2011, 477 (7365): 452-456. 10.1038/nature10382.
Smith SA, Wilson NG, Goetz FE, Feehery C, Andrade SCS, Rouse GW, Giribet G, Dunn CW: Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature. 2011, 480 (7377): 364-367. 10.1038/nature10526.
Bernt M, Bleidorn C, Braband A, Dambach J, Donath A, Fritzsch G, Golombek A, Hadrys H, Jühling F, Meusemann K, Middendorf M, Misof B, Perseke M, Podsiadlowski L, von Reumont B, Schierwater B, Schlegel M, Stöger M, Simon S, Stadler PF, Stöger I, Struck TH: A comprehensive analysis of bilaterian mitochondrial genomes and phylogeny. Mol Phylogenet Evol. 2013, 69 (2): 352-364. 10.1016/j.ympev.2013.05.002.
S-i Y, Fukuda N, Nakamura M, Aoyama T, Oshima T: Long-term conservation of six duplicated structural genes in cephalopod mitochondrial genomes. Mol Biol Evol. 2004, 21 (11): 2034-2046. 10.1093/molbev/msh227.
Yuan Y, Li Q, Yu H, Kong L: The complete mitochondrial genomes of six heterodont bivalves (Tellinoidea and Solenoidea): variable gene arrangements and phylogenetic implications. PLoS One. 2012, 7 (2): e32353-10.1371/journal.pone.0032353.
Cunha R, Grande C, Zardoya R: Neogastropod phylogenetic relationships based on entire mitochondrial genomes. BMC Evol Biol. 2009, 9 (1): 210-10.1186/1471-2148-9-210.
Medina M, Lal S, VallËs Y, Takaoka TL, Dayrat BA, Boore JL, Gosliner T: Crawling through time: transition of snails to slugs dating back to the Paleozoic, based on mitochondrial phylogenomics. Mar Genom. 2011, 4 (1): 51-59. 10.1016/j.margen.2010.12.006.
Boore J, Brown W: Complete DNA sequence of the mitochondrial genome of the black chiton, Katharina tunicata. Genetics. 1994, 138: 423-443.
Veale AJ, Williams L, Tsai P, Thakur V, Lavery S: The complete mitochondrial genomes of two chiton species(Sypharochiton pelliserpentisandSypharochiton sinclairi)obtained using Illumina next generation sequencing. In Mitochondrial DNA. 2014. doi:10.3109/19401736.2014.905846.,
Dreyer H, Steiner G: The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca). Mol Phylogenet Evol. 2004, 31: 605-617. 10.1016/j.ympev.2003.08.007.
Boore JL, Medina M, Rosenberg LA: Complete sequences of the highly rearranged molluscan mitochondrial genomes of the scaphopod Graptacme eborea and the bivalve Mytilus edulis. Mol Biol Evol. 2004, 21 (8): 1492-1503. 10.1093/molbev/msh090.
Stöger I, Stöger M: Mitogenomics does not resolve deep molluscan relationships (yet?). Mol Phylogenet Evol. 2013, 69 (2): 376-392. 10.1016/j.ympev.2012.11.017.
Ojala D, Montoya J, Attardi G: tRNA punctuation model of RNA processing in human mitochondria. Nature. 1981, 290 (5806): 470-474. 10.1038/290470a0.
Mwinyi A, Meyer A, Bleidorn C, Lieb B, Bartolomaeus T, Podsiadlowski L: Mitochondrial genome sequence and gene order of Sipunculus nudus give additional support for an inclusion of Sipuncula into Annelida. BMC Genomics. 2009, 10 (1): 27-10.1186/1471-2164-10-27.
Dunn CW, Hejnol A, Matus DQ, Pang K, Browne WE, Smith SA, Seaver E, Rouse GW, Obst M, Edgecombe GD, Sørensen MV, Haddock SHD, 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 (7188): 745-749. 10.1038/nature06614.
Breton S, Beaupré HD, Stewart DT, Hoeh WR, Blier PU: The unusual system of doubly uniparental inheritance of mtDNA: isn’t one enough?. Trends Genet. 2007, 23 (9): 465-474. 10.1016/j.tig.2007.05.011.
Meyer A, Todt C, Mikkelsen N, Lieb B: Fast evolving 18S rRNA sequences from Solenogastres (Mollusca) resist standard PCR amplification and give new insights into mollusk substitution rate heterogeneity. BMC Evol Biol. 2010, 10 (1): 70-10.1186/1471-2148-10-70.
Sambrook J, Fritsch E, Maniatis T: Molecular Cloning: A Laboratory Manual. Volume 1, 2, 3. 1989, Cold Spring Harbor Laboratory Press, New York
Folmer O, Black M, Hoeh W, Lutz RA, Vrijenhoek RC: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. 1994, Mol Mar Biol, Biotech
Burger G, Lavrov D, Forget L, Lang B: Sequencing complete mitochondrial and plastid genomes. Nat Protoc. 2007, 2: 603-614. 10.1038/nprot.2007.59.
Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Pütz J, Middendorf M, Stadler PF: MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013, 69 (2): 313-319. 10.1016/j.ympev.2012.08.023.
Wyman S, Jansen R, Boore J: Automatic annotation of organellar genomes with DOGMA. Bioinformatics. 2004, 20 (17): 3252-3255. 10.1093/bioinformatics/bth352.
Schattner P, Brooks A, Lowe T: The tRNAscan-SE, snoscan and snoGPS web servers for the detection of tRNAs and snoRNAs. Nucleic Acids Res. 2005, 33 (suppl 2): W686-W689. 10.1093/nar/gki366.
Laslett D, Canbäck B: ARWEN: a program to detect tRNA genes in metazoan mitochondrial nucleotide sequences. Bioinformatics. 2008, 24 (2): 172-175. 10.1093/bioinformatics/btm573.
Plazzi F, Ribani A, Passamonti M: The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with Conchifera. BMC Genomics. 2013, 14 (1): 409-10.1186/1471-2164-14-409.
Maynard B, Kerr L, McKiernan J, Jansen E, Hanna P: Mitochondrial DNA sequence and gene organization in Australian backup abalone Haliotis rubra (leach). Mar Biotechnol (NY). 2005, 7: 645-658. 10.1007/s10126-005-0013-z.
Chen H-X, Sundberg P, Wu H-Y, Sun S-C: The mitochondrial genomes of two nemerteans, Cephalothrix sp. (Nemertea: Palaeonemertea) and Paranemertes cf. peregrina (Nemertea: Hoplonemertea). Mol Biol Rep. 2011, 38 (7): 4509-4525. 10.1007/s11033-010-0582-4.
Stechmann A, Schlegel M: Analysis of the complete mitochondrial DNA sequence of the brachiopod Terebratulina retusa places Brachiopoda within the protostomes. Proc Roy Soc Lond B Biol Sci. 1999, 266 (1433): 2043-2052. 10.1098/rspb.1999.0885.
S-i Y, Iseto T, Asakawa S, Sasaki T, Shimizu N, Yamagishi A, Oshima T, Hirose E: Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny. Mol Phylogenet Evol. 2008, 47 (2): 612-628. 10.1016/j.ympev.2008.02.013.
Helfenbein KG, Boore JL: The mitochondrial genome of Phoronis architecta—comparisons demonstrate that phoronids are lophotrochozoan protostomes. Mol Biol Evol. 2004, 21 (1): 153-157. 10.1093/molbev/msh011.
Katoh K, Standley DM: MAFFT Multiple Sequence Alignment Software Version 7: improvements in performance and usability. Mol Biol Evol. 2013, 30 (4): 772-780. 10.1093/molbev/mst010.
Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552. 10.1093/oxfordjournals.molbev.a026334.
Glez-Peña D, Gómez-Blanco D, Reboiro-Jato M, Fdez-Riverola F, Posada D: ALTER: program-oriented conversion of DNA and protein alignments. Nucleic Acids Res. 2010, 38 (suppl 2): W14-W18. 10.1093/nar/gkq321.
Akaike H: Information theory and an extension of the maximum likelihood principle. 2nd International Symposium on Information Theory. Edited by: Petrov BN, Csaki F. 1973, Akademiai Kiado, Budapest, Hungary, 267-281.
Lanfear R, Calcott B, Ho SYW, Guindon S: PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012, 29 (6): 1695-1701. 10.1093/molbev/mss020.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22 (21): 2688-2690. 10.1093/bioinformatics/btl446.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.
We are grateful to the late Christoffer Schander (University of Bergen, Norway) for initiating this study. We thank three anonymous reviewers for insightful comments on a previous version of the manuscript. Use of R.V. Hans Brattström for collection of specimens was granted by the Department of Biology, University of Bergen, Norway. We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI). DO was supported by a FPI grant (BES-2008-009562) of the Spanish Ministerio de Economía y Competitividad (MINECO). II was supported by a JAE-pre PhD Fellowship (Consejo Superior de Investigaciones Científicas, CSIC, and the European Social Fund, ESF) and by a Humboldt Research Fellowship for Postdoctoral Researchers (Alexander von Humboldt Foundation) during the last stage of this work. CG is currently a “Ramón y Cajal” postdoctoral fellow supported by the MINECO and the Universidad Autónoma de Madrid, Spain. The work was partly funded by projects of the MINECO to RZ (CGL2007-60954 and CGL2010-18216) and to CG (CGL2011-29916).
The authors declare that they have no competing interests.
RZ conceived the study. CT conducted the sampling and identification of individuals. DO carried out laboratory work. DO, CG and II performed phylogenetic and comparative analyses. All authors contributed to the writing and improving of the manuscript, as well as read and approved the final version.
Electronic supplementary material
Additional file 1: Figure S1.: Comparison of mitochondrial gene orders of Scutopus ventrolineatus and Chaetoderma nitidulum. Table S1. Complete mitochondrial genomes used in the phylogenetic analyses. Table S2. Fragments of nuclear ribosomal proteins used in phylogenetic analyses. Table S3. Best-fit partitions and models selected by Protein Partition Finder. (PDF 569 KB)
About this article
Cite this article
Osca, D., Irisarri, I., Todt, C. et al. The complete mitochondrial genome of Scutopus ventrolineatus (Mollusca: Chaetodermomorpha) supports the Aculifera hypothesis. BMC Evol Biol 14, 197 (2014). https://doi.org/10.1186/s12862-014-0197-9
- Nuclear ribosomal proteins