Next Article in Journal
Support Interval for Two-Sample Summary Data-Based Mendelian Randomization
Previous Article in Journal
Peanut AhmTERF1 Regulates Root Growth by Modulating Mitochondrial Abundance
Previous Article in Special Issue
The Complete Mitochondrial Genome of Entemnotrochus rumphii, a Living Fossil for Vetigastropoda (Mollusca: Gastropoda)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
School of Energy and Environmental Engineering, Hebei University of Technology, Tianjin 300401, China
2
Institute of Oceanology, Chinese Academy of Sciences, Qingdao 266071, China
3
University of Chinese Academy of Sciences, Beijing 101400, China
4
Laboratory for Marine Ecology and Environmental Science, Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China
5
Center for Ocean Mega-Science, Chinese Academy of Sciences, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Genes 2023, 14(1), 210; https://doi.org/10.3390/genes14010210
Submission received: 30 December 2022 / Revised: 10 January 2023 / Accepted: 10 January 2023 / Published: 13 January 2023
(This article belongs to the Special Issue Genetic Evolution of Marine Shellfish)

Abstract

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

1. Introduction

The arrangements of mitochondrial genes appear to have changed very little for some phyla of animals. For instance, the most studied arthropods and vertebrates share near-identical gene arrangements within the phyla. Such q consistency of mitogenomic features in the model animals (fruit flies, mouses and frogs, etc.) has led to a ‘textbook description’ of the conservation mitochondrial genome. Yet, more thorough and extensive mitochondrial genome sampling from animals has discovered numerous instances in which such features are absent. The second largest phylum, Mollusca, is particularly rich in exceptions, and the research in recent years has produced a huge list of unexpected departures from the textbook description in Mollusca, including extreme size variations, radical genome rearrangements, frequent gene duplications and losses, the insertion of novel genes and a bizarre inheritance pattern known as “doubly uniparental inheritance” (DUI) [1]. Thus, it is required for us to continue studying the striking variation in the molluscan mitochondrial genome architecture, molecular function and intergenerational transmission, as these can hold great promise for uncovering animal mitochondrial genome-specific biological processes, histories and functions. Molluscan mitochondrial genomes offer a model system for investigating evolution, the variety of functions that mitochondria serve in organismal physiology and the numerous applications of phylogenetic and population biology that can be derived.
The Scaphopoda is a rare and small clade, consisting of nearly 600 described, extant species [2]. They have a common name, tusk-shell, because of their conspicuous ivory-shaped shell, which opens at both the anterior and posterior ends. The anterior opening is wider, with a burrowing foot that extends out. Unlike most other molluscan relatives, they have no gills or ctenidium, instead they use the mantle for respiration. Scaphopods live infaunal in sand or mud to collect foraminifers, which are their main prey, with numerous captacula and grind them up using their large radula apparatus. Scaphopods are ecologically important bioturbators, and in some circumstances, they can be extremely prevalent in benthic marine communities. For example, in a recent biodiversity assessment of the Prince Gustav Channel in the Weddell Sea, Antarctica [3], the most abundant molluscan class was Scaphopoda, with over 11,000 specimens having been collected in just six epibenthic sledge casts (54% of the molluscan diversity sampled). Scaphopods are also evolutionarily important as they are the closest living relatives of Gastropoda and Bivalvia [4,5,6], the two most diverse and economically important classes of mollusks. Yet, the exact placement of scaphopods within Conchifera has been undecided in terms of molluscan phylogeny. Scaphopods have received much less research attention than most other classes of mollusks have. Most of our knowledge about scaphopod diversity and evolution is from species descriptions and morphological studies [7,8], whereas relatively few pieces of work have studied this group from a molecular perspective. To date, only four sets of mitochondrial genome data have been published in this clade, and it is speculated that the scaphopod mitochondrial genomes have the highest degree of gene re-arrangement compared to those of the other remaining classes of mollusc. To verify this and facilitate a more comprehensive understanding of mitochondrial genome evolution in molluscs, here, we present the complete mitogenome of one representative species, the dentaliid scaphopod, P. vernedei.

2. Materials and Methods

2.1. Sampling and DNA Isolation

An adult animal of P. vernedei was sampled in March 2020 from the east China Sea (27°31′32.73″ N, 22°30′39.13″ E). The species identification was carried out by the Department of marine organism taxonomy and phylogeny, IOCAS (Institute of Oceanology, Chinese Academy of Sciences), based on morphology and 16S sequencing. The foot was dissected and preserved in liquid nitrogen. The total genomic DNA was further isolated using CTAB methods.

2.2. Mitochondrial Genome Sequencing and Assembling

The whole mitogenome was sequenced by novogene co. ltd, Tianjing, China. In brief, by using the NexteraXT DNA Library Preparation Kit, we constructed the genomic libraries, with an average length of 350 bp, and performed the sequencing of 150 bp paired-end runs on the Illumina Novaseq 6000 platform (San Diego, CA, USA). To meet the quality control, the raw data were further processed by NGS QC Toolkit v2.3.3 as per manual instructions [9]. Then, the de novo assembling of mitochondrial genome was carried out in SPAdes 3.11.0 [10].

2.3. Annotation and Characterization of Scaphopod Mitochondrial Genome

The annotation of the P. vernedei mitochondrial genome was performed by using MITOS Web Server (http://mitos.bioinf.uni-leipzig.de/index.py, with accession on 12 September 2022). We then used SnapGene Viewer for the manual inspection and corrections of start/stop codons in each gene. The location and length of all of the protein-coding genes (PCGs) were checked under Open Reading Frame Finder (ORFfinder Home, NCBI). The nucleotide composition, composition skew and codon usage of PCGs, as well as the relative synonymous codon usage (RSCU), were analyzed and calculated in PhyloSuite v1.2.0 [11]. We applied the online Proksee tool (https://proksee.ca/, accessed on 27 September 2022) to visualize the mitochondrial genome map. To predict the cloverleaf structures of 22 transfer RNAs, we used tRNAscan-Se [12,13] and MITOS according to the user manual, and we applied ViennaRNA Web Services (http://rna.tbi.univie.ac.at/, accessed on 2 October 2022) to visualize the secondary structures.
To further explore the gene arrangement among the Scaphopoda and other molluscan lineages, we selected 21 complete mitochondrial genomes from multiple lineages (Table 1) and downloaded them from NCBI database with their accession numbers. Of these, 4 species belonged to the Scaphopoda, 3 species belonged to the Protobranchia (presumably the most basal bivalve lineage) and 3 species belonged to vetigastropoda (presumably the most basal gastropod lineages). The Lingulata from Brachiopoda was selected as the outgroups.

2.4. Phylogenetic Placement

To explore the phylogenetic placement of P. vernedei within Scaphopoda, a phylogenetic analysis was carried out as follows. Briefly, 13 protein-coding genes were obtained under PhyloSuite v1.2.2 [11], then, the alignment of their amino acid sequences was performed by Mafft v7 [14]. We used Gblocks 0.91b [15] with default settings to remove those misaligned positions. All of the protein-coding genes were then concatenated into a matrix under PhyloSuite v1.2.2 [11]. We applied ModelFinder [16] with AICc (so-called corrected Akaike Information Criterion, which has a more precise approximation in small samples) to infer a best-fitting partition model. Finally, IQ-Tree v2.1.2 [17] with 1000 ultrafast bootstrap replications was applied to construct a maximum likelihood (ML) tree, which was processed in FigTree v1.4.4 (https://github.com/rambaut/figtree/releases, accessed on 25 December 2022) for visualization.

3. Results and Discussion

3.1. General Features of Mitochondrial Genome

The total length of the circular mitogenome of P. vernedei is 14,519 bp (Figure 1), and this is comparable to those of other four known scaphopod mitogenomes (14,869 of Antalis entails, 14,492 of G. eborea, 13,932 of S. lobatum and 13,790 of Gadilida sp.). Molluscs vary largely in their mitochondrial genomes size. The scaphopods, together with heterobranch gastropods (~13.6–14.1 kb) [18,19,20], are among the smallest mitochondrial genomes that have been reported so far [1], while the largest ones come from Arcidae clams (e.g., the Scapharca kagoshimensis mitogenome is estimated to be ~56.2 kb) [21]. The P. vernedei mitogenome contains 37 mtDNAs that are typical to most metazoans (Table 2 and Figure 1), including thirteen protein coding genes (PCGs), twenty-two tRNAs and two ribosomal RNAs (12S RNA and 16S RNA). Five PCGs (cox1, nad2, cytb, nad3 and atp6) and twelve tRNAs (tRNA-SerTCT, tRNA-Asn, tRNA-His, tRNA-Gly, tRNA-Arg, tRNA-SerTGA, tRNA-Ala, tRNA-LeuTAG, tRNA-LeuTAA, tRNA-Glu, tRNA-Trp and tRNA-Cys) are encoded on the heavy strand, whereas the other eight PCGs, ten two tRNAs and two rRNAs are distributed on the light strand. Among these coding genes, a total of 21 intergenic regions were detected, spanning from 1 to 119 bp, wherein the longest one was situated between 12S RNA and tRNA-Met. In contrast, nine overlaps were detected, ranging from 1 to 18 bp. The longest overlapped region was located between tRNA-Met and 16S RNA (18 bp), which was followed by the cox1 and tRNA-SerTCT overlap (14 bp), as well as the tRNA-His and cox2 overlap (14 bp).
The nucleotide composition was 36.46% for A, 35.45% for T, 13.24% for C and 14.85% for G in P. vernedei, and as a whole, it was biased toward A+T content, at 71.91%. In other Dentaliida species, the A+T contents exceed 70%, indicating that the Dentaliida mitogenomes are more biased toward A+T (Table 3). In addition, P. vernedei has a positive total A-T skew and a negative G-C skew. Among the eleven molluscs plus two outgroups that were analyzed in present study, positive A-T skews were found in six species, while negative G-C skews were found in nine species (Table 3), showing no apparent consistency.

3.2. Protein-Coding Genes (PCGs)

P. vernedei has 13 PCGs in its mitochondrial genome, whose length is 10,842 bp in total, and this accounts for 74.67% of the complete mitogenome. Their size ranges from 159 (atp8) to 1647 bp (nad5), and this is generally similar to most other molluscs, for example, E. rumphii [22]. In P. vernedei, these PCGs use four different triplets as the start codon, including ATG, ATT, TTG and ATA. Such a usage is similar to its sister Dentaliida G. eborea, which also uses ATG, ATT, TTG and ATA as a start codon, but not TTG [23]. In contrast, in other molluscs such as some vetigastropods, ATG is the only start codon used for all of the PCGs [22]. These 13 PCGs contain various types of stop codons; TAA was found in seven PCGs (cox1, cox2, cox3, nad3, nad4, nad4l and nad5), while TAG was found in three PCGs (atp8, cytb and nad6) and T was found in the remaining three PCGs (atp6, nad1 and nad2). Using TAA and TAG as stop codons is widely common in molluscs, and using T is also prevalent in scaphopods such as G. eborea [23].
In terms of the base composition of PCGs, T accounts for most (42.88%), followed by 28.24% of A, 15.22% of C, and 13.66% of G. The A-T content in both the individual PCG and concatenated matrices surpass 60%, suggesting that A-T is more preferred in P. vernedei PCGs. The A-T skews in all of the PCGs are negative, indicating they are biased toward T. With regard to the G-C skew, eleven PCGs show a positive G-C skew, whereas the remaining two PCGs (atp6 and nad6) show a negative one.
As shown in Figure 2, the distribution of relative synonymous codon usage (RSCU) of P. vernedei mitogenome indicates that Leu, Phe and Ile are the top three most frequently utilized amino acids, while UUA-Leu, UCU-Ser2, ACU-Thr and CCU-Pro are the top four most frequently used codons. Among these twenty-two amino acids, nine of them (i.e., Ala, Arg, Gly, Leu1, Pro, Ser1, Ser2, Thr and Val) are encoded by four codons, while the others are encoded by two codons, which is common in other molluscs.
Annotating genes in the molluscan mitochondrial genome is notoriously challenging. While it is easy to identify the protein-encoding genes by searching the similarity of orthologous counterparts, it could be difficult to infer the valid start codon and stop codon. This is partly because molluscs usually have multiple alternatives for both start and end codons [1]. For instance, in addition to ATG, start codons in molluscs also include TTG, ATY, ATA, and GTG (which are responsible for leucine, isoleucine, methionine, and valine encoding, respectively). Each of them has a match to the trnM anticodon for at least two nucleotides, suggesting their dual role in mitochondria, that is, they serve as the transfer RNA for both methionine and, as for protein initiation, formyl-methionine. For the stop codon, in some cases, molluscs do not just terminate at a complete TAA or TAG stop codon but terminate on just a TA or T that is completed to a TAA stop by polyadenylation of the transcript. Another thing that aggravates the difficulty in inferring start and stop codons is the pervasive gene overlapping with the upstream/downstream gene, even when on the same strand. This requires us to take a full consideration of the extent of the ORF evolutionary conservation and manually check the start/stop codon carefully.

3.3. Ribosomal RNAs and Transfer RNAs

There are two ribosomal RNAs: one is 12S ribosomal RNA (rrnS), with a length of 636 bp, and another one is 16S ribosomal RNA (rrnL), with a length of 1200 bp. The former one, rrnS, is located between tRNA-Met and tRNA-Thr, while the latter one, rrnL, is located between tRNA-Val and tRNA-Met. For rrnS, the A-T skew and G-C skew are 0.034 and 0.098, respectively; for rrnL, they are 0.057 and 0.22. Both two ribosomal RNAs are biased toward T and G.
P. vernedei possesses 22 transfer RNA genes, with their size ranging from 59 (tRNA-Arg) to 71 bp (tRNA-Thr). The A-T content in the total transfer RNAs is 75.5%, which also exhibits an obvious A-T bias. Transfer RNAs usually have typical clover-leaf secondary structures, yet in P. vernedei, truncated Di-hydroxy uridine (DHU) arms without paired bases occur in two tRNAs, tRNA-Arg and tRNA-Ser1 (Figure 3). The DHU arm of tRNA-Ser1 is still eleven bases long, but that of the tRNAs Arg has been reduced to five bases. In its close sister species G. eborea, five tRNAs have DHU arms, including tRNA-Ser1, tRNA-Arg, tRNA-Thr and tRNA-Ile [23]. Thus, DHU lost in tRNA-Ser1 and tRNA-Arg can be seen as a synapomorphy in this clade. In addition, for tRNA-Ser1, the DHU loss can be widely founded in other metazoan mitogenomes [24], suggesting it is either an ancestral feature or a convergent consequence.

3.4. Phylogenetic Placement

To understand the placement of P. vernedei and the relationship between the scaphopods, bivalves and gastropods, we analyzed three bivalves, three gastropods, five current available scaphopods and two outgroups (Platynereis and Lingula). The phylogenetic analysis was performed based on the alignment of 13 protein-coding genes (Figure 4A). In the present study, the scaphopods were grouped into a monophyletic clade, within which P. vernedei was placed as a sister to G. eborea with BS = 95. Together with A. entalis, they form a clade with full support, representing Dentaliida. We also recovered the monophyletic Gadilida containing S. lobatum and Gadilida sp. with full support. However, the bootstrap support value between Dentaliida and Gadilida is not high. This is partly because the current matrix contains insufficient phylogenetic signals to resolve this node. We argue that more taxa sampling is needed for revising the scaphopod internal systematics.
The phylogenetic placement of Scaphopoda within Conchifera remains to be a notorious problem. There are at least six competing hypotheses showing discrepancies in scaphopod placement: the Scaphopoda can be placed as sister to (1) Bivalvia in ‘Diasoma’ concept [25,26,27], (2) Cephalopoda in the ‘helcionellid’ concept [28,29,30], (3) Gastropoda [6,31,32], (4) Gastropoda-Cephalopoda in the Haszprunar hypothesis or (5) Pleistomollusca (Bivalvia+ Gastropoda) [4]. Though Cephalopoda and Monoplacophora mitogenomes were not included in present study, we found that Bivalvia and Gastropoda form a clade (Pleistomollusca), and Scaphopoda is a sister to it, with BS = 75. The transcriptome-based phylogenomics results are consistent with the previous findings from large-scale transcriptome-based phylogenomics [4]. However, we hoped that the complete mitochondrial genomes would provide new character sets to help resolve the molluscan relationships, but it is now apparent that molluscan mitogenomes did not provide the useful resolution as they were expected to [2,33]. Even when we were using whole mitogenomes data sets, gastropods can be diphyletic. Mitogenomes can fail to resolve deep relationships of molluscan major lineages, partly because mitogenome has a more accelerated evolutionary rate than the genomes for most metazoans so [34], and more importantly, because molluscan mitogenomes have a surprisingly wide variation in size, rearrangements, duplications/losses and in the introduction of novel genes [1]. Such diverse changes are severally biased to different lineages, making it not a good tool for resolving deep molluscan relationships. We thus argue to use whole genome data to confidently resolve the placement of Scaphopoda in the future.

3.5. Mitochondrial Gene Arrangement

Molluscs break the rules of mitochondrial genomes, as they are one of the few metazoan groups that show extensive gene rearrangement. Such prevalent rearrangements in this phylum can hamper the reconstruction of relationships as discussed above, but on the other hand, such variation occurring among the molluscan lineages is also helpful in understanding the underlying mechanisms of gene reordering. Here, we compared four scaphopod gene arrangements (Gadillida sp. was not included because of its incomplete assembly and annotation, and only these four are currently available and complete) with three protobranch bivalves and three vetigastropods (Figure 4B) because Protobranchia and Vetigastropoda showed deep conservation in the mitochondrial gene order, with only a few steps of rearrangements occurring comparing to those for chitons and other lophotrochozoans, thus, these are supposed to be most approximal to the ancestral states [2].
We found that P. vernedei shared an identical gene order with dentaliid sister G. eborea, yet in another dentaliid A. entails there is an inversion region, namely, trnF-trnK-nad5-trnD-nad4-nad4L. Compared to Protobranchia and Vetigastropoda, the gene order in P. vernedei is highly re-arranged, with only one region (nad6-trnP-nad1) showing the same pattern. Such neighboring relationships between those three blocks can also be seen in Cephalopoda and Monoplacophora [35], indicating that this order is a symplesiomorphy in Conchifera. Surprisingly, in gadilid S. lobatum, the gene order is completely divergent, and it has lost the conserved nad6-trnP-nad1 relationship and shows no conservative region with both the dentaliid groups and other molluscan lineages. This suggests extensive gene rearrangement not only occurred when Scaphopoda split from other lineages, but it also occurred in the split between Gadilida and Dentaliid.

4. Conclusions

Our study presents a complete mitogenome of P. vernedei, whose size was 14,519 bp. It contains typical thirty-seven mtDNAs, including thirteen PCGs, twenty-two tRNAs and two rRNAs. The arrangement of the mitochondrial genes is identical to those of dentaliid G. eborea, but it is quite divergent from those of Gadilida and other molluscan lineages, suggesting that Scaphopoda has gone through substantial gene rearrangement. The extensive rearrangement in scaphopods implies that there has been little conservation, thus making it not a good choice for resolving the placement of Scaphopoda among the molluscs. Given that previous phylogeny analyses using morphological data, transcriptome data, and mitogenome data all leave the scaphopod placement undetermined, we argue for the use of whole genome data, which comprise the most comprehensive genetic information to resolve this problem in the future. Since the scaphopod showed more frequent gene rearrangement than other molluscan lineages did, it provides a model system to study how such substantial rearrangements are in concert with animal physiology and adaptative evolution.

Author Contributions

Y.W. and H.S. conceived of the presented idea. T.Z., Y.W. and H.S. performed sampling, sequencing, and informatics analysis. The first draft of the manuscript was written by T.Z., Y.W. and H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by the Young Elite Scientists Sponsorship Program by CAST (2021QNRC001), the National Science Foundation of China (No. 32002409), the earmarked fund for Modern Agro-industry Technology Research System (CARS-47).

Institutional Review Board Statement

All procedures involving the care and use of animals conformed to the measures for the administration of animal experiment safety of the Institute of Oceanology, Chinese Academy of Science, and were approved by the Institutional Animal Care and Ethics Committee of the Institute of Oceanology, Chinese Academy of Science (Approval Code: Haifa 2020 (37), Approval Date: 17 June 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

The mitochondrial genome of P. vernedei is available from GeneBank under the accession OQ198438.

Acknowledgments

We appreciate the Oceanographic Data Center, IOCAS, for providing a computing server. We also thanks Kevin Kocot from the University of Alabama for providing the annotations of S. lobatum and Gadilida sp.

Conflicts of Interest

The authors have no relevant financial or non-financial interest to disclose.

Abbreviations

BSBootstrap
DHUDi-hydroxy uridine
DUIdoubly uniparental inheritance
MLmaximum likelihood
PCGprotein coding gene
RSCUrelative synonymous codon usage

References

  1. Ghiselli, F.; Gomes-dos-Santos, A.; Adema, C.M.; Lopes-Lima, M.; Sharbrough, J.; Boore, J.L. Molluscan mitochondrial genomes break the rules. Philos. Trans. R. Soc. B 2021, 376, 20200159. [Google Scholar] [CrossRef] [PubMed]
  2. Ponder, W.F.; Lindberg, D.R.; Ponder, J.M. Biology and Evolution of the Mollusca, 1st ed.; CRC Press: Boca Raton, FL, USA, 2020; Volume 1, p. 924. [Google Scholar]
  3. Anderson, M.P.; Fenberg, P.B.; Griffiths, H.J.; Linse, K. Macrobenthic mollusca of the Prince Gustav Channel, eastern Antarctic Peninsula: An area undergoing colonisation. Front. Mar. Sci. 2021, 8, 771369. [Google Scholar] [CrossRef]
  4. Kocot, K.M.; Cannon, J.T.; Todt, C.; Citarella, M.R.; Kohn, A.B.; Meyer, A.; Santos, S.R.; Schander, C.; Moroz, L.L.; Lieb, B. Phylogenomics reveals deep molluscan relationships. Nature 2011, 477, 452–456. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Kocot, K.M. Recent advances and unanswered questions in deep molluscan phylogenetics. Am. Malacol. Bull. 2013, 31, 195–208. [Google Scholar] [CrossRef]
  6. Smith, S.A.; Wilson, N.G.; Goetz, F.E.; Feehery, C.; Andrade, S.; Rouse, G.W.; Giribet, G.; Dunn, C.W. Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature 2011, 480, 364–367. [Google Scholar] [CrossRef]
  7. Reynolds, P.D.; Okusu, A. Phylogenetic relationships among families of the Scaphopoda (Mollusca). Zool. J. Linn. Soc. 1999, 126, 131–154. [Google Scholar] [CrossRef]
  8. Reynolds, P.D. The phylogeny and classification of Scaphopoda (Mollusca): An assessment of current resolution and cladistic reanalysis. Zool. Scr. 1997, 26, 13–21. [Google Scholar] [CrossRef]
  9. Patel, R.K.; Jain, M.J. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data. PLoS ONE 2012, 7, e30619. [Google Scholar] [CrossRef]
  10. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  11. Zhang, D.; Gao, F.; Jakovlić, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef]
  12. Lowe, T.M.; Chan, P.P. tRNAscan-SE On-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, W54–W57. [Google Scholar] [CrossRef] [PubMed]
  13. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef]
  14. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.; Von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [Green Version]
  17. Minh, B.Q.; Schmidt, H.A.; Chernomor, O.; Schrempf, D.; Woodhams, M.D.; Von Haeseler, A.; Lanfear, R. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era. Mol. Biol. Evol. 2020, 37, 1530–1534. [Google Scholar] [CrossRef] [Green Version]
  18. Hatzoglou, E.; Rodakis, G.C.; Lecanidou, R. Complete sequence and gene organization of the mitochondrial genome of the land snail Albinaria coerulea. Genetics 1995, 140, 1353–1366. [Google Scholar] [CrossRef]
  19. Kurabayashi, A.; Ueshima, R. Evolution, Complete sequence of the mitochondrial DNA of the primitive opisthobranch gastropod Pupa strigosa: Systematic implication of the genome organization. Mol. Biol. Evol. 2000, 17, 266–277. [Google Scholar] [CrossRef]
  20. White, T.R.; Conrad, M.M.; Tseng, R.; Balayan, S.; Golding, R.; de Frias Martins, A.M.; Dayrat, B.A. Ten new complete mitochondrial genomes of pulmonates (Mollusca: Gastropoda) and their impact on phylogenetic relationships. BMC Evol. Biol. 2011, 11, 295. [Google Scholar] [CrossRef] [Green Version]
  21. Kong, L.; Li, Y.; Kocot, K.M.; Yang, Y.; Qi, L.; Li, Q.; Halanych, K.M. Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size. Mol. Phylogenet. Evol. 2020, 150, 106857. [Google Scholar] [CrossRef]
  22. Wang, Y.; Ma, P.; Zhang, Z.; Li, C.; Liu, Y.; Chen, Y.; Wang, J.; Wang, H.; Song, H. The Complete Mitochondrial Genome of Entemnotrochus rumphii, a Living Fossil for Vetigastropoda (Mollusca: Gastropoda). Genes 2022, 13, 2061. [Google Scholar] [CrossRef] [PubMed]
  23. Dreyer, H.; Steiner, G. Evolution, The complete sequence and gene organization of the mitochondrial genome of the gadilid scaphopod Siphonondentalium lobatum (Mollusca). Mol. Phylogenet. Evol. 2004, 31, 605–617. [Google Scholar] [CrossRef] [PubMed]
  24. Garey, J.R.; Wolstenholme, D.R. Platyhelminth mitochondrial DNA: Evidence for early evolutionary origin of a tRNAserAGN that contains a dihydrouridine arm replacement loop, and of serine-specifying AGA and AGG codons. J. Mol. Evol. 1989, 28, 374–387. [Google Scholar] [CrossRef] [PubMed]
  25. Runnegar, B.; Pojeta, J., Jr. Molluscan Phylogeny: The Paleontological Viewpoint: The early Paleozoic fossil record shows how living and extinct molluscan classes originated and diversified. Science 1974, 186, 311–317. [Google Scholar] [CrossRef]
  26. Pojeta, J., Jr. The paleontology of rostroconch mollusks and the early history of the Phylum Mollusca. US Geol. Surv. Prof. Pap. 1976, 968, 1–54. [Google Scholar]
  27. De Lacaze Duthiers, H. Histoire de L’organisation, du Développement, des Moeurs et des Rapports Zoologiques du Dentale; Librairie de Victor Masson: Paris, France, 1858. [Google Scholar]
  28. Waller, T.R.; Johnston, P.; Haggart, J. Origin of the molluscan class Bivalvia and a phylogeny of major groups. Bivalves Eon Evol. 1998, 1, 5. [Google Scholar]
  29. Hatschek, B. Lehrbuch der Zoologie, eine morphologische Ubersicht des Thierreiches zur Einfuhrung in das Studium dieser Wissenschaft; Gustav Fischer Verlag: Jena, Germany, 1888. [Google Scholar]
  30. Steiner, G.; Dreyer, H. Molecular phylogeny of Scaphopoda (Mollusca) inferred from 18S rDNA sequences: Support for a Scaphopoda–Cephalopoda clade. Zool. Scr. 2003, 32, 343–356. [Google Scholar] [CrossRef]
  31. Plate, L. Über den Bau und die Verwandtschaftsbeziehungen der Solenoconchen. Zool. Jahrbücher Der Anat. 1892, 5, 301–386. [Google Scholar]
  32. Bronn, H.G. Dr. HG Bronn’s Klassen und Ordnungen des Thier-Reichs, wissenschaftlich dargestellt in Wort und Bild. Dritter Band. Mollusca (Weichthiere); CF Winter’sche-Verlagshandlung: Leipzig, Germany, 1894. [Google Scholar]
  33. Grande, C.; Templado, J.; Zardoya, R. Evolution of gastropod mitochondrial genome arrangements. BMC Ecol. Evol. 2008, 8, 61. [Google Scholar] [CrossRef] [Green Version]
  34. Gissi, C.; Iannelli, F.; Pesole, G. Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species. Heredity 2008, 101, 301–320. [Google Scholar] [CrossRef] [Green Version]
  35. Guerra, D.; Bouvet, K.; Breton, S. Mitochondrial gene order evolution in Mollusca: Inference of the ancestral state from the mtDNA of Chaetopleura apiculata (Polyplacophora, Chaetopleuridae). Mol. Phylogenet. Evol. 2018, 120, 233–239. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The overall mitogenome map of P. vernedei. From the inside to the outside, the innermost scale indicates the length, the inner circle indicates the G-C skew, the middle circle indicates the G-C content, and the outer circle indicates gene arrangement. The G-C skew was centered at zero while the G-C content circus was centered at 50%.
Figure 1. The overall mitogenome map of P. vernedei. From the inside to the outside, the innermost scale indicates the length, the inner circle indicates the G-C skew, the middle circle indicates the G-C content, and the outer circle indicates gene arrangement. The G-C skew was centered at zero while the G-C content circus was centered at 50%.
Genes 14 00210 g001
Figure 2. Distribution of (a) relative synonymous codon usage (RSCU) and (b) amino acids frequency in P. vernedei mitogenome.
Figure 2. Distribution of (a) relative synonymous codon usage (RSCU) and (b) amino acids frequency in P. vernedei mitogenome.
Genes 14 00210 g002
Figure 3. Cloverleaf structures of the 21 tRNA genes of P. vernedei mitochondrial genome.
Figure 3. Cloverleaf structures of the 21 tRNA genes of P. vernedei mitochondrial genome.
Genes 14 00210 g003
Figure 4. Phylogenetic tree (a) and gene orders (b). (a) Three protobranch bivalves, three vetigastropods, five currently available scaphopods, and 2 outgroup species were selected, and the phylogenetic tree was constructed based on the amino acid sequences from 13 protein-coding genes. Bootstrap (BS) support values are shown at each node. (b) Gene order of PCGs among Vetigastropoda, Protobranchia, and Scaphopoda. Gadillida sp. was not included as its assembly and annotation are not complete. The arrangement graph is visualized by iTOL and the boxes with different colors and labels represent different gene blocks.
Figure 4. Phylogenetic tree (a) and gene orders (b). (a) Three protobranch bivalves, three vetigastropods, five currently available scaphopods, and 2 outgroup species were selected, and the phylogenetic tree was constructed based on the amino acid sequences from 13 protein-coding genes. Bootstrap (BS) support values are shown at each node. (b) Gene order of PCGs among Vetigastropoda, Protobranchia, and Scaphopoda. Gadillida sp. was not included as its assembly and annotation are not complete. The arrangement graph is visualized by iTOL and the boxes with different colors and labels represent different gene blocks.
Genes 14 00210 g004
Table 1. Species selected for phylogenetic analysis and gene order comparison.
Table 1. Species selected for phylogenetic analysis and gene order comparison.
PhylumClassSubclassOrderSpeciesAccession Number
MolluscaGastropodaVetigastropodaPleurotomariida Entemnotrochus rumphiiOP354269
LepetellidaHaliotis rufescensKX260956
Fissurella volcanoJN790612
BivalviaProtobranchiaSolemyidaPetrasma pervernicosaKY244080
Solemya velumJQ728447
Nuculida Nucula nucleusEF211991
Scaphopoda Dentaliida Antalis entalisMN098312
G. eboreaAY484748
P. vernedeiThis study
GadilidaSiphonodentalium lobatumAY342055
Gadilida sp.MN104231
AnnelidaPolychaeta ErrantiaPhyllodocida Platynereis dumeriliiAF178678
Brachiopoda Lingulata Lingulida Lingula anatinaKX774482
Table 2. Organization of the coding genes in P. vernedei mitochondrial genome.
Table 2. Organization of the coding genes in P. vernedei mitochondrial genome.
Gene NameLocationLength (bp)CodonAnticodonIntergenic RegionStrand
StartEndStartStop
tRNA-Cys16868 GCA0H
cox16316191557TTGTAA −6H
tRNA-Ser(tct)1606167065 TCT−14H
tRNA-Asn1675174066 GTT4H
nad217612700940ATGT 20H
cytb275538881134ATATAG 54H
tRNA-His3879394062 GTG−10H
cox239274616690ATGTAA −14L
tRNA-Gln4617468064 TTG0L
tRNA-Gly4687475367 TCC6H
cox347565535780ATGTAA 2L
tRNA-Tyr5536560065 GTA0L
tRNA-Arg5599565759 TCG−2H
tRNA-Ser(tga)5657572064 TGA−1H
nad657146199486ATGTAG −7L
tRNA-Pro6200626465 TGG0L
nad162657147883ATGT 0L
atp871487306159ATGTAG 0L
tRNA-Ile7314738168 GAT7L
tRNA-Thr7384745471 TGT2L
12S RNA74568091636 1L
tRNA-Met8211827666 CAT119L
16S RNA825994581200 −18L
tRNA-Val9549961163 TAC90L
tRNA-Ala9621968363 TGC9H
nad3969910,037339TTGTAA 15H
tRNA-Leu(tag)10,03910,10365 TAG1H
tRNA-Leu(taa)10,10510,16965 TAA1H
tRNA-Glu10,17910,24365 TTC9H
tRNA-Trp10,24810,31366 TCA4H
tRNA-Phe10,32210,38968 GAA8L
tRNA-Lys10,39010,45869 TTT0L
nad510,45812,1041647ATTTAA −1L
tRNA-Asp12,14312,20765 GTC38L
nad412,20913,4771269ATTTAA 1L
nad4l13,54913,839291ATTTAA 71L
atp613,85314,519667ATTT 13H
Table 3. Statistics of nucleotide composition in selected molluscan mitogenomes.
Table 3. Statistics of nucleotide composition in selected molluscan mitogenomes.
Species SelectedSize (bp)A%T%G%C%A-T SkewG-C SkewAT Proportion
E. rumphii15,79535.2129.9814.4320.390.08023−0.1711265.19%
H. rufescens16,64635.3924.9313.7525.930.17339−0.3071960.32%
F. volcano17,57525.7235.8326.6911.76−0.164310.3881961.55%
P. pervernicosa16,55434.4531.8412.6421.070.03946−0.2499666.29%
S. velum15,66035.5332.5813.1418.740.04341−0.1756568.11%
N. nucleus13,67132.1331.7815.6620.420.00549−0.1319763.92%
A. entalis14,86935.3040.4112.3911.90−0.067420.0199375.71%
G. eborea14,49236.9837.1412.6713.21−0.00214−0.0208174.12%
P. vernedei14,51936.4635.4513.2414.850.01408−0.0573871.91%
S. lobatum13,93232.0836.2018.9212.81−0.060360.1925868.27%
Gadilida sp.13,79032.0832.8214.3520.74−0.01134−0.1821464.91%
P. dumerilii15,61931.2232.9215.4120.45−0.02646−0.1407164.14%
L. anatina24,87626.2836.5321.2915.90−0.163080.1448362.81%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, T.; Wang, Y.; Song, H. The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei. Genes 2023, 14, 210. https://doi.org/10.3390/genes14010210

AMA Style

Zhang T, Wang Y, Song H. The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei. Genes. 2023; 14(1):210. https://doi.org/10.3390/genes14010210

Chicago/Turabian Style

Zhang, Tianzhe, Yunan Wang, and Hao Song. 2023. "The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei" Genes 14, no. 1: 210. https://doi.org/10.3390/genes14010210

APA Style

Zhang, T., Wang, Y., & Song, H. (2023). The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei. Genes, 14(1), 210. https://doi.org/10.3390/genes14010210

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop