Abstract
We described the complete mitochondrial genome of the Ctenophore Beroe cucumis, which is a circular molecule of 10,487 bp in length. The new mitochondrial genome comprised only 12 genes, making it one of the smallest animals’ mtDNA. Both nucleotide substitution and gene order rearrangements exhibited extreme high evolutionary rate in mitogenomes of Ctenophore. The phylogenetic analysis based on mitogenomics failed to reveal the basal position of Ctenophore within metazoan, owing to the extreme evolutionary rate. Based on the available Ctenophora mitogenomes, we found the optimized primers designed by Geller et al. for DNA barcoding suited for the taxon.
Introduction
Ctenophora (comb jellies) are a small phylum composed of about 200 species (Mills Citation2017; http://www.marinespecies.org/). Although they possess complex nervous and mesoderm-derived muscular systems, the Ctenophora were believed to be the basal animal lineage (Mcfadden et al. Citation2011; Moroz et al. Citation2014). Mitochondrial genomes/genes have been proved to be good tools for studies on molecular evolution and species identification (Geller et al. Citation2013). However, the under-sampled status of Ctenophora’s mitochondrial genomes made the applications of the marker in trouble. Lack of the mitochondrial genomes complicated the primer design of DNA barcode in the Ctenophora. Here, we sequenced a mitogenome belonging to Beroe cucumis.
The sample of B. cucumis was collected in the Jiaozhou Bay (36°06′N, 124°14′E) in August 2016 and preserved under −80 °C in the Center of Deep-Sea Research, Institute of Oceanology, CAS, China (sample code: JZ-20160801BC). The genomic DNA was extracted using Omega plant DNA mini kit (catalog no. D3485-01; Omega Bio-tek, Norcross, GA, USA). The complete mtDNA was obtained by assembling (clc genomic workbench v11) the 10G 150 bp pair-end reads produced by Illumina X10 (). The annotated mtDNA was deposited in the GenBank with the accession number MK361035.
Figure 1. Phylogenetic position of Ctenophora within Metazoan using concatenated proteins in mitochondrial genomes. Numbers at each code are Bayesian posterior probabilities (percentage). GenBank accession numbers of each species were listed in the tree.
The complete mitochondrial genome of B. cucumis is 10,487 bp in length and highly AT-biased. The A + T content (81.51%) is similar to other Ctenophora species (Arafat et al. Citation2018), but higher than typical invertebrate species. The AT-skew and GC-skew of this mtDNA were −0.398 and −0.068 respectively. Like other Ctenophora species, the mitochondrial genome coding genes located on the same strand and lacked all the tRNAs, ATP6/8, and ND6 genes. The PCGs start with ATA (ND2, ND3, and ND5) and ATT (others), and terminated with TAA except for ND2 with TAG. This mitochondrial gene order is identical to that of B. forskalii, but different from other Ctenophora. Gene order of Ctenophora was highly rearranged except for those from the congeneric species as reported before (Arafat et al. Citation2018). The high evolutionary rate (Pett et al. Citation2011) was confirmed by the relative low similarity between the sibling Beroe species based on the blastn comparison results (76% identity for the whole mtDNA and maximum 88% identity for the COX1 gene). Generally, 97% similarity was adopted as the gold standard to distinguish the congeneric species.
We investigated the phylogenetic relationships of ctenophore using the amino acids of mitochondrial genome coding genes. Best partitioning scheme was evaluated in partition finder 2.1. As reported in previous studies, our analysis was unable to resolve the basal position of the ctenophores within metazoans owing to the extreme evolutionary rate (Osigus et al. Citation2013). In our BI (Mr. Bayes) inferred trees, the Ctenophora were positioned as the sister clade of Bilateria. However, the results confirmed the monophyly of Beroe and the high mitochondrial evolutionary rate within Ctenophora.
Finally, we found that the primer-targeted regions of Ctenophora gave a similar level of mismatches compared to other taxa, suggesting the optimized primers for DNA barcoding (Geller et al. Citation2013) suited for the Ctenophora.
Acknowledgements
We thank the faculty of Jiaozhou Bay Marine Ecosystem Research Station, Institute of Oceanology, Chinese Academy of Sciences, for their help of collecting samples.
Disclosure statement
The authors report no conflicts of interest and are responsible for the content and writing of the paper.
Additional information
Funding
References
- Arafat H, Alamaru A, Gissi C, Huchon D. 2018. Extensive mitochondrial gene rearrangements in Ctenophora: insights from benthic Platyctenida. BMC Evol Biol. 18:65.
- Geller J, Meyer C, Parker M, Hawk H. 2013. Redesign of PCR primers for mitochondrial cytochrome c oxidase subunit I for marine invertebrates and application in all-taxa biotic surveys. Mol Ecol Resour. 13:851–861.
- Mcfadden CS, Benayahu Y, Pante E, Thoma JN, Nevarez PA, France SC. 2011. Limitations of mitochondrial gene barcoding in Octocorallia. Mol Ecol Resour. 11:19–31.
- Mills CE. 2017. Phylum Ctenophora: list of all valid species names [Internet Document]. [accessed 2017 Mar 30]. http://faculty.washington.edu/cemills/Ctenolist.html. 1998-Present.
- Moroz LL, Kocot KM, Citarella MR, Dosung S, Norekian TP, Povolotskaya IS, Grigorenko AP, Dailey C, Berezikov E, Buckley KM, et al. 2014. The ctenophore genome and the evolutionary origins of neural systems. Nature. 510:109–114.
- Osigus H, Eitel M, Bernt M, Donath A, Schierwater B. 2013. Mitogenomics at the base of Metazoa. Mol Phylogenet Evol. 69:339–351.
- Pett WF, Ryan J, Pang KC, Mullikin JQ, Martindale MD, Baxevanis A, Lavrov D. 2011. Extreme mitochondrial evolution in the Ctenophore Mnemiopsis leidyi: insights from mtDNA and the nuclear genome. Mitochondr DNA. 22:130–142.