Research Article |
Corresponding author: Joong-Ki Park ( jkpark@ewha.ac.kr ) Academic editor: Andrew Davinack
© 2022 Yucheol Lee, Joong-Ki Park.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Lee Y, Park J-K (2022) Complete mitochondrial genome of Conus lischkeanus Weinkauff, 1875 (Neogastropoda, Conidae) and phylogenetic implications of the evolutionary diversification of dietary types of Conus species. ZooKeys 1088: 173-185. https://doi.org/10.3897/zookeys.1088.78990
|
The family Conidae, commonly known as cone snails, is one of the most intriguing gastropod groups owing to their diverse array of feeding behaviors (diets) and toxin peptides (conotoxins). Conus lischkeanus Weinkauff, 1875 is a worm-hunting species widely distributed from Africa to the Northwest Pacific. In this study, we report the mitochondrial genome sequence of C. lischkeanus and inferred its phylogenetic relationship with other Conus species. Its mitochondrial genome is a circular DNA molecule (16,120 bp in size) composed of 37 genes: 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. Phylogenetic analyses of concatenated nucleotide sequences of 13 PCGs and two ribosomal RNA genes showed that C. lischkeanus belongs to the subgenus Lividoconus group, which is grouped with species of the subgenus Virgiconus, and a member of the largest assemblage of worm-hunting (vermivorous) species at the most basal position in this group. Mitochondrial genome phylogeny supports the previous hypothesis that the ancestral diet of cone snails was worm-hunting, and that other dietary types (molluscivous or piscivorous) have secondarily evolved multiple times from different origins. This new, complete mitochondrial genome information provides valuable insights into the mitochondrial genome diversity and molecular phylogeny of Conus species.
Cone snail, dietary type evolution, Lividoconus
The genus Conus Linnaeus, 1758 is a well-known predatory gastropod group that produces venomous peptides, called conotoxins, to capture prey and defend against predators (
The implementation of new sequencing technologies (e.g., next-generation sequencing; NGS) and various bioinformatics tools has allowed mitochondrial genome sequencing to be markedly easier, cost-effective, and widely used for studying phylogeny in various metazoan groups, including Conus species (
Conus lischkeanus specimen was collected from Moonseom, Jeju Island, Korea, preserved in 95% ethanol solution, and deposited in the Marine Mollusk Resource Bank of Korea (MMRBK; voucher specimen no. MMRBK6746) in Seoul, Korea. The specimen was morphologically identified based on shell characters, which include a conical last whorl covered with yellow-brown periostracum and an angular shoulder. Total genomic DNA was extracted from the foot tissue using an E.Z.N.A. mollusc DNA kit (Omega Bio-tek, Norcross, GA, USA) following the manufacturer’s instructions.
Whole-genome sequencing libraries were prepared using the MGIEasy DNA library prep kit (BGI, Shenzhen, China) according to the manufacturer’s instructions and quantified using the QuantiFluor ssDNA System (Promega Corporation, Madison, WI, USA). Sequencing was conducted on the MGISEQ-2000 system with 150 base-pair reads. A total of 48,608,637 raw reads were obtained, and adapter-trimmed using a skewer program (
To determine the relationship between C. lischkeanus and other Conus species, phylogenetic analyses were performed for the nucleotide sequences of 13 protein-coding genes (PCGs) and two rRNA genes from 39 complete or nearly complete mitochondrial genomes of the family Conidae (Table
Complete mitochondrial genomes used for phylogenetic analysis in this study.
Family | Species | Diet | GenBank | Reference |
---|---|---|---|---|
Conidae | Conus victoriae | Molluscivorous | — |
|
Conus gloriamaris | Molluscivorous | KU996360 | — | |
Conus textile | Molluscivorous | DQ862058 |
|
|
Conus episcopatus | Molluscivorous | — |
|
|
Conus marmoreus | Molluscivorous | — |
|
|
Conus nobilis | Molluscivorous | KX263253 |
|
|
Conus ermineus | Piscivorous | KY864977 |
|
|
Conus tulipa | Piscivorous | KR006970 |
|
|
Conus consors | Piscivorous | KF887950 |
|
|
Conus striatus | Piscivorous | KX156937 |
|
|
Conus betulinus | Vermivorous | — |
|
|
Conus sponsalis | Vermivorous | — |
|
|
Conus arenatus | Vermivorous | — |
|
|
Conus goudeyi | Vermivorous | KY864975 |
|
|
Conus ebraeus | Vermivorous | — |
|
|
Conus coronatus | Vermivorous | — |
|
|
Conus miliaris | Vermivorous | — |
|
|
Conus pseudonivifer | Vermivorous | KY864969 |
|
|
Conus venulatus | Vermivorous | KX263250 |
|
|
Conus ateralbus | Vermivorous | KY864970 |
|
|
Conus byssinus | Vermivorous | KY864973 |
|
|
Conus pulcher | Vermivorous | KY864972 |
|
|
Conus genuanus | Vermivorous | KY864974 |
|
|
Conus hybridus | Vermivorous | KX263252 |
|
|
Conus guanche | Vermivorous | KY801847 |
|
|
Conus ventricosus | Vermivorous | KX263251 |
|
|
Conus miruchae | Vermivorous | KY864971 |
|
|
Conus borgesi | Vermivorous | EU827198 |
|
|
Conus infinitus | Vermivorous | KY864967 |
|
|
Conus spurius | Vermivorous | KY864976 |
|
|
Conus virgo | Vermivorous | — |
|
|
Conus quercinus | Vermivorous | KY609509 |
|
|
Conus lischkeanus | Vermivorous | OL632021 | This study | |
Conus lividus | Vermivorous | — |
|
|
Conus tabidus | Vermivorous | KY864968 |
|
|
Conus lenavati | Vermivorous | — |
|
|
Conus tribblei | Vermivorous | KT199301 |
|
|
Conus imperialis | Vermivorous | — |
|
|
Conus capitaneus | Vermivorous | KX155573 |
|
|
Conasprella wakayamaensis | Vermivorous | KX263254 |
|
|
Californiconus californicus | All | KX263249 |
|
|
Profundiconus teramachii | Vermivorous | KX263256 |
|
|
Borsoniidae | Tomopleura sp. | — | KX263259 |
|
Conus lischkeanus is widely distributed from East Africa to the western Pacific (
Gene | Start | Stop | Strand direction | Length (bp) | Codon (start) | Codon (stop) | Overlapping regions | Intergenic spacers |
---|---|---|---|---|---|---|---|---|
cox1 | 1 | 1,548 | H | 1,548 | ATG | TAA | 1 | 166 |
cox2 | 1,715 | 2,401 | H | 687 | ATG | TAA | — | — |
tRNA-Asp (trnD) (gtc) | 2,402 | 2,468 | H | 67 | — | — | ||
atp8 | 2,469 | 2,630 | H | 162 | ATG | TAA | — | 6 |
atp6 | 2,637 | 3,359 | H | 723 | ATG | TAA | — | 11 |
tRNA-Met (trnM) (cat) | 3,371 | 3,438 | L | 68 | — | 12 | ||
tRNA-Tyr (trnY) (gta) | 3,451 | 3,516 | L | 66 | — | 1 | ||
tRNA-Cys (trnC) (gca) | 3,518 | 3,582 | L | 65 | — | — | ||
tRNA-Trp (trnW) (tca) | 3,583 | 3,648 | L | 66 | — | — | ||
tRNA-Gln (trnQ) (ttg) | 3,646 | 3,711 | L | 66 | 3 | 24 | ||
tRNA-Gly (trnG) (tcc) | 3,736 | 3,802 | L | 67 | — | 35 | ||
tRNA-Glu (trnE) (ttc) | 3,838 | 3,902 | L | 65 | — | — | ||
small subunit rRNA (rrnS) | 3,903 | 4,854 | H | 952 | — | — | ||
tRNA-Val (trnV) (tac) | 4,855 | 4,921 | H | 67 | — | — | ||
large subunit rRNA (rrnL) | 4,922 | 6,296 | H | 1,375 | — | — | ||
tRNA-Leu1 (trnL1) (tag) | 6,297 | 6,366 | H | 70 | — | 6 | ||
tRNA-Leu2 (trnL2) (taa) | 6,373 | 6,441 | H | 69 | — | — | ||
nad1 | 6,442 | 7,383 | H | 942 | ATG | TAG | — | 16 |
tRNA-Pro (trnP) (tgg) | 7,400 | 7,468 | H | 69 | — | — | ||
nad6 | 7,469 | 7,975 | H | 507 | ATG | TAA | — | 13 |
cob | 7,989 | 9,128 | H | 1,140 | ATG | TAA | — | 11 |
tRNA-Ser2 (trnS2) (tga) | 9,140 | 9,204 | H | 65 | — | 16 | ||
tRNA-Thr (trnT) (tgt) | 9,221 | 9,289 | L | 69 | — | 22 | ||
nad4L | 9,312 | 9,608 | H | 297 | ATG | TAG | — | — |
nad4 | 9,602 | 10,984 | H | 1,383 | ATG | TAG | 7 | — |
tRNA-His (trnH) (gtg) | 10,984 | 11,049 | H | 66 | 1 | — | ||
nad5 | 11,050 | 12,765 | H | 1,716 | ATG | TAA | — | — |
tRNA-Phe (trnF) (gaa) | 12,765 | 12,829 | H | 65 | 1 | — | ||
D-loop | 12,830 | 13,415 | H | 586 | — | — | ||
cox3 | 13,416 | 14,195 | H | 780 | ATG | TAA | — | 34 |
tRNA-Lys (trnK) (ttt) | 14,230 | 14,298 | H | 69 | — | 9 | ||
tRNA-Ala (trnA) (tgc) | 14,308 | 14,374 | H | 67 | — | 22 | ||
tRNA-Arg (trnR) (tcg) | 14,397 | 14,465 | H | 69 | — | 11 | ||
tRNA-Asn (trnN) (gtt) | 14,477 | 14,545 | H | 69 | — | 12 | ||
tRNA-Ile (trnI) (gat) | 14,558 | 14,626 | H | 69 | — | 5 | ||
nad3 | 14,632 | 14,985 | H | 354 | ATG | TAA | — | 15 |
tRNA-Ser1 (trnS1) (gct) | 15,001 | 15,068 | H | 68 | — | — | ||
nad2 | 15,069 | 1 | H | 1,053 | ATG | TAA | — | — |
The lengths of 13 PCGs of C. lischkeanus mitochondria range from 162 bp (atp8) to 1,716 bp (nad5) and contain 3,751 codons, excluding termination codons. The base composition of PCGs is 26.3% A, 39.2% T, 17.5% G, and 17.0% C, and the overall AT content was 65.5%, which is very similar to that of the entire mitochondrial genome sequence (AT content of 66.1%; Table
Nucleotide composition of the mitochondrial genome of Conus lischkeanus.
Nucleotide sequence | Length (bp) | A (%) | C (%) | G (%) | T (%) | A+T (%) | G+C (%) |
---|---|---|---|---|---|---|---|
Entire sequence | 16,120 | 29.0 | 16.3 | 17.6 | 37.1 | 66.1 | 33.9 |
Protein coding sequence | 11,292 | 26.3 | 17.0 | 17.5 | 39.2 | 65.5 | 34.5 |
Codon position* | |||||||
1st | 3,751 | 26.9 | 17.2 | 24.7 | 31.2 | 58.1 | 41.9 |
2nd | 3,751 | 18.3 | 20.9 | 16.6 | 44.2 | 62.5 | 37.5 |
3rd | 3,751 | 33.4 | 13.1 | 11.4 | 42.1 | 75.5 | 24.5 |
Ribosomal RNA gene sequence | 2,327 | 35.5 | 14.4 | 18.4 | 31.6 | 67.2 | 32.8 |
Transfer RNA gene sequence | 1,481 | 34.0 | 16.2 | 17.7 | 32.1 | 66.1 | 33.9 |
D-loop region sequence | 586 | 31.1 | 15.4 | 17.6 | 35.8 | 67.1 | 32.9 |
Twenty-two tRNA genes were found in the mitochondrial genome of C. lischkeanus. The length of tRNA genes range from 65 bp (trnC, trnE, trnS2, and trnF) to 70 bp (trnL1) (Table
Phylogenetic analysis using ML and BI methods yield similar results with respect to the tree topology, as shown in Fig.
Phylogenetic relationships of the genus Conus based on concatenated nucleotide sequences (13 protein coding genes plus two rRNA genes). Numbers above branches are statistical support values for ML (bootstrap values, > 70)/BI (posterior probability values, > 0.7). *: determined in this study.
This work was supported by National Marine Biodiversity Institute of Korea (2022M01100) and the Basic Science Research Program of the National Research Foundation of Korea (NRF) funded by the Ministry of Education (no. 2020R1I1A1A01074213; 2020R1A2C2005393).