Characterization of mitochondrial genome of sea cucumber Stichopus horrens: a novel gene arrangement in Holothuroidea

A novel deep-benthic sea cucumber species of Benthodytes (Holothuroidea, Elasipodida, Psychropotidae) and its comprehensive mitochondrial genome sequencing and evolutionary analysis

BACKGROUND

The holothurians, commonly known as sea cucumbers, are marine organisms that possess significant dietary, nutritional, and medicinal value. However, the National Center for Biotechnology Information (NCBI) currently possesses only approximately 70 complete mitochondrial genome datasets of Holothurioidea, which poses limitations on conducting comprehensive research on their genetic resources and evolutionary patterns. In this study, a novel species of sea cucumber belonging to the genus Benthodytes, was discovered in the western Pacific Ocean. The genomic DNA of the novel sea cucumber was extracted, sequenced, assembled and subjected to thorough analysis.

RESULTS

The mtDNA of Benthodytes sp. Gxx-2023 (GenBank No. OR992091) exhibits a circular structure spanning 17,386 bp, comprising of 13 protein-coding genes (PCGs), 24 non-coding RNAs (2 rRNA genes and 22 tRNA genes), along with two putative control regions measuring 882 bp and 1153 bp, respectively. It exhibits a high AT% content and negative AT-skew, which distinguishing it from the majority of sea cucumbers in terms of environmental adaptability evolution. The mitochondrial gene homology between Gxx-2023 and other sea cucumbers is significantly low, with less than 91% similarity to Benthodytes marianensis, which exhibits the highest level of homology. Additionally, its homology with other sea cucumbers is below 80%. The mitogenome of this species exhibits a unique pattern in terms of start and stop codons, featuring only two types of start codons (ATG and ATT) and three types of stop codons including the incomplete T. Notably, the abundance of AT in the Second position of the codons surpasses that of the First and Third position. The gene arrangement of PCGs exhibits a relatively conserved pattern, while there exists substantial variability in tRNA. Evolutionary analysis revealed that it formed a distinct cluster with B. marianensis and exhibited relatively distant phylogenetic relationships with other sea cucumbers.

CONCLUSIONS

These findings contribute to the taxonomic diversity of sea cucumbers in the Elasipodida order, thereby holding significant implications for the conservation of biological genetic resources, evolutionary advancements, and the exploration of novel sea cucumber resources.

Complete mitochondrial genome of a sea cucumber, Euapta godeffroyi (Echinodermata, Holothuroidea, Apodida, Synaptidae)

We determined the complete mitochondrial genome sequence of a holothurian Euapta godeffroyi belonging to the order Apodida. The complete mitogenome of E. godeffroyi was 16,410 bp in length and consisted of 13 protein-coding genes (PCGs), two ribosomal RNA genes, and 22 transfer RNA genes. The orders of PCGs and rRNAs did not match those of any recorded holothurian mitogenomes. The maximum likelihood (ML) phylogenetic tree placed E. godeffroyi as the sister group to chiridotid species and supported the monophyly of the order Apodida.

The Complete Mitochondrial Genome of Stichopus naso (Aspidochirotida: Stichopodidae: Stichopus) and Its Phylogenetic Position

The mitochondrial genome is widely used to study the molecular evolution of and perform phylogenetic analyses on animals. In this study, the complete mitochondrial genome (mitogenome) of Stichopus naso was sequenced. The mitogenome was 16,239 bp in length and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The genome composition showed positive AT-skew (0.023) and negative GC-skew (−0.158). The order of the mitochondrial genes was consistent with those from the Stichopus and Isostichopus species, whereas it was different from those of other species of Aspidochirotida. The phylogenetic analysis, based on the nucleotide sequences of 13 PCGs through the methods of Bayesian inference (BI) and maximum likelihood (ML), indicated that S. naso has close relationships with S. horrens and S. monotuberculatus, and belongs to a member of Stichopodidae. Our study provides a reference mitogenome for further molecular evolution studies and phylogenetic research on sea cucumbers.

The first two complete mitogenomes of the order Apodida from deep-sea chemoautotrophic environments: New insights into the gene rearrangement, origin and evolution of the deep-sea sea cucumbers

The deep-sea ecosystem is considered as the largest and most remote biome of the world. It is meaningful and important to elucidate the life origins by exploring the origin and adaptive genetic mechanisms of the large deep-sea organisms. Sea cucumbers (Holothuroidea) are abundant and economically important group of echinoderms, living from the shallow-waters to deep-sea. In this study, we present the mitochondrial genomes of the sea cucumber Chiridota heheva and Chiridota sp. collected from the deep-sea cold seep and hydrothermal vent, respectively. This is the first reported mitochondrial genomes from the order Apodida. The mitochondrial genomes of C. heheva (17,200 bp) and Chiridota sp. (17,199 bp) display novel gene arrangements with the first protein-coding gene rearrangements in the class Holothuroidea. Bases composition analysis showed that the A + T content of deep-sea holothurians were significantly higher than that of the shallow-water groups. We compared the arrangement of genes from the 24 available holothurian mitogenomes and found that the transposition, reverse transposition and tandem-duplication-random-losses (TDRL) may be involved in the evolution of mitochondrial gene arrangements in Holothuroidea. Phylogenetic analysis revealed that the Apodida clustered with Elasipodida, forming two basal deep-sea holothurian clades. The divergence between the deep-sea and shallow-water holothurians was located at 386.93 Mya, during the Late Devonian. Mitochondrial protein-coding genes of deep-sea holothurians underwent relaxed purifying selection. There are 57 positive selected amino acids sites for some mitochondrial genes of the three deep-sea clades, implying they may involve in the adaption of deep-sea sea cucumbers.

Comparison of the complete mitochondrial genome of Phyllophorus liuwutiensis (Echinodermata: Holothuroidea: Phyllophoridae) to that of other sea cucumbers

Sea cucumber species are abundant (>1400 species) and widely distributed globally. mtDNA sequencing is frequently used to identify the phylogenetic and evolutionary relationships among species. However, there are no reports on the mitochondrial genome of Phyllophorus liuwutiensis. Here, we performed mtDNA sequencing of P. liuwutiensis to examine its phylogenetic relationships with other echinoderms. Its mitochondrial genome (15 969 bp) contains 37 coding genes, including 13 protein‐coding genes, 22 tRNA genes and 2 rRNA genes. Except for one protein‐coding gene (nad6) and five tRNA genes encoded on the negative strand, all other genes were encoded on the positive strand. The mitochondrial bases of P. liuwutiensis were composed of 29.55% T, 22.16% C, 35.64% A and 12.64% G. The putative control region was 703 bp in length. Seven overlapping regions (1–10 bp) were found. The noncoding region between the genes ranged from 1 to 130 bp in length. One putative control region has been found in the P. liuwutiensis mitogenome. All of the tRNA genes were predicted to fold into a cloverleaf structure. In addition, we compared the gene arrangements of six echinoderms, revealing that the gene order of P. liuwutiensis was a new arrangement.

The complete mitochondrial genome sequence of Stichopus variegatus (Echinodermata: Holothuroidea) and phylogenetic studies of Echinodermata

At present, there exist some confusing issues on the species classification and phylogeny in Echinodermata. In this study, we first determined and described the complete mitochondrial genome of Stichopus variegatus. The complete mitogenome sequence had a circular mapping molecular with the total length of 16,315 bp and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and a putative control region. To further validate the newly determined sequences, phylogenetic trees involving all the Holothuroidea and other Echinodermata species available in GenBank Database were constructed. These results would be used for the species identification and further phylogenetic studies of Echinodermata.

The complete mitochondrial genome of a tropical sea cucumber, Holothuria leucospilota

The mitochondrial genome of Holothuria leucospilota was 15,906 bp in length, containing 13 protein-coding genes (PCGs), 22 tRNA genes, and two rRNA genes. There were four initiation codons (ATG, ATT, ATC, and ATA) for the PCGs, and the termination codon of most PCGs was TAA, except for nad4 (TAG) and nad6 (TAG). Only one PCG (nad6) and five tRNA genes (tRNA^Ser(UCN), tRNA^Gln, tRNA^Ala, RNA^Val, and tRNA^Asp) were encoded on the light chain; the other genes were encoded on the heavy chain. H. leucospilota was most closely related to Holothuria scabra in a phylogenetic tree.

Complete mitochondrial genome of Holothuria leucospilata (Holothuroidea, Holothuriidae) and phylogenetic analysis

The complete Holothuria leucospilata mitochondrial genome was determined and analyzed in this work. It had a circular mapping molecular with a total length of 15,904 bp and contained 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 putative control region. Phylogenetic analysis showed that H. leucospilata clustered together with Holothuria scabra and Holothuria forskali. The complete mitochondrial genome provided in this work would be used for elucidation of Holothuroidea conservation genetics and evolutionary relationships.

Complete mitochondrial genome of Benthodytes marianensis (Holothuroidea: Elasipodida: Psychropotidae): Insight into deep sea adaptation in the sea cucumber

Complete mitochondrial genomes play important roles in studying genome evolution, phylogenetic relationships, and species identification. Sea cucumbers (Holothuroidea) are ecologically important and diverse members, living from the shallow waters to the hadal trench. In this study, we present the mitochondrial genome sequence of the sea cucumber Benthodytes marianensis collected from the Mariana Trench. To our knowledge, this is the first reported mitochondrial genome from the genus Benthodytes. This complete mitochondrial genome is 17567 bp in length and consists of 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Most of these genes are coded on the positive strand except for one protein-coding gene (nad6) and five tRNA genes which are coded on the negative strand. Two putative control regions (CRs) have been found in the B. marianensis mitogenome. We compared the order of genes from the 10 available holothurian mitogenomes and found a novel gene arrangement in B. marianensis. Phylogenetic analysis revealed that B. marianensis clustered with Peniagone sp. YYH-2013, forming the deep-sea Elasipodida clade. Positive selection analysis showed that eleven residues (24 S, 45 S, 185 S, 201 G, 211 F and 313 N in nad2; 108 S, 114 S, 322 C, 400 T and 427 S in nad4) were positively selected sites with high posterior probabilities. We predict that nad2 and nad4 may be the important candidate genes for the further investigation of the adaptation of B. marianensis to the deep-sea environment.

The first complete mitochondrial genome of the Mariana Trench Freyastera benthophila (Asteroidea: Brisingida: Brisingidae) allows insights into the deep-sea adaptive evolution of Brisingida

Starfish (phylum Echinodermata) are ecologically important and diverse members of marine ecosystems in all of the world's oceans, from the shallow water to the hadal zone. The deep sea is recognized as an extremely harsh environment on earth. In this study, we present the mitochondrial genome sequence of Mariana Trench starfish Freyastera benthophila, and this study is the first to explore in detail the mitochondrial genome of a deep-sea member of the order Brisingida. Similar to other starfish, it contained 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes (duplication of two tRNAs: trnL and trnS). Twenty-two of these genes are encoded on the positive strand, while the other 15 are encoded on the negative strand. The gene arrangement was identical to those of sequenced starfish. Phylogenetic analysis showed the deep-sea Brisingida as a sister taxon to the traditional members of the Asteriidae. Positive selection analysis indicated that five residues (8 N and 16 I in atp8, 47 D and 196 V in nad2, 599 N in nad5) were positively selected sites with high posterior probabilities. Compared these features with shallow sea starfish, we predict that variation specifically in atp8, nad2, and nad5 may play an important role in F. benthophila's adaptation to deep-sea environment.

Complete mitochondrial genome of Parastichopus californicus (Aspidochirotida: Stichopodidae)

In this study, we first determined and described the complete mitogenome sequence of Parastichopus californicus, which was 16 727 bp in length. The mitochondrial genome had the canonical mitochondrial gene content, including 13 protein-coding genes, 2 rRNA genes and 22 tRNA genes. The overall base composition of the heavy-strand was 31.5% A, 18 % G, 20.6% C and 29.9% T, with a high A + T content of 61.4%. ML phylogenetic tree indicated that P. californicus and P. nigripunctatus were clustered in one branch belonging to the genus Parastichopus. This conclusion was identical to the former result by the methods of morphological taxonomy.

Complete mitochondrial genome of the sandfish Holothuria scabra (Holothuroidea, Holothuriidae)

The complete mitochondrial genome (mitogenome) sequence of Holothuria scabra, an economically and ecologically important tropical sea cucumber, was first sequenced and annotated. The mitochondrial DNA is 15,779 bp in length and contains 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and a 456 bp putative control region, of which gene order is identical to the echinoderm ground pattern. Comparative analyses between H. scabra and other holothurians revealed three new findings: (1) the mitogenome of H. scabra is highly compact having five regions with overlapping genes and least intergenic nucleotides among the sequenced holothurians, only accounting for 3.58% of its mitogenome; (2) the genus Holothuria mitogenomes show a pattern of high interspecies divergence at the 13 PCGs, and the genetic p-distance reaches 25.68% between H. scabra and H. forskali; (3) the incomplete stop codon T of cox2 shared with H. forskali may be a common feature in the genus Holothuria.

Recombination and evolution of duplicate control regions in the mitochondrial genome of the Asian big-headed turtle, Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3' end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.