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

Complete Mitochondrial Genomes and Phylogenetic Analysis of Genus Henricia (Asteroidea: Spinulosida: Echinasteridae)

The genus Henricia is known to have intraspecific morphological variations, making species identification difficult. Therefore, molecular phylogeny analysis based on genetic characteristics is valuable for species identification. We present complete mitochondrial genomic sequences of Henricia longispina aleutica, H. reniossa, and H. sanguinolenta for the first time in this study. This study will make a significant contribution to our understanding of Henricia species and its relationships within the class Asteroidea. Lengths of mitochondrial genomes of the three species are 16,217, 16,223, and 16,194 bp, respectively, with a circular form. These genomes contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, and a D-loop. The gene order and direction aligned with other asteroid species. Phylogenetic relationship analysis showed that our Henricia species were in a monophyletic clade with other Henricia species and in a large clade with species (Echinaster brasiliensis) from the same family. These findings provide valuable insight into understanding the phylogenetic relationships of species in the genus Henricia.

Phylogenetic relationships and adaptation in deep-sea carideans revealed by mitogenomes

The deep-sea environment is characterized by extreme and inhospitable conditions, including oxygen depletion, low temperatures, high pressure, absence of light, and limited food availability. Mitochondria and mitogenomes play a crudial role in aerobic respiration to generate energy for eukaryotes. Here, using the Illumina Hiseq 4000 platform, we performed mitogenome sequencing for five deep-sea caridean species: Lebbeus shinkaiae, Lebbeus Formosus, Glyphocrangon regalis, Heterocarpus dorsalis, and Heterocarpus laevigatus, and five deep-sea caridean mitogenomes were assembled and identified. Each of the five mitogenomes contained 13 protein-coding genes, 2 rRNAs and 22 tRNAs. Specific elements, such as tandem repeats and AT-rich sequences, were observed in the control regions of Lebbeus formosus and Lebbeus shinkaiae, potentially take a role in regulating mitochondrial genome replication and transcription. The gene order of all obtained mitogenomes follows caridean ancestral type organization. Phylogenetic analysis shows a robustly supported phylogenetic tree for the infraorder Caridea. The monophyly of the families included in this study was strongly supported. This study supports the monophyly of Oplophoroidea, but rejects the monophyletic status of Nematocarcinoidea, Crangonoidea, and Alpheoidea. At the genus level, Plesionika is polyphyletic and Rimicaris is paraphyletic in our analysis. Furthermore, Paralebbeus may be considered invalid and synonymous with Lebbeus. Positive selection analysis reveals evidence for adaptive changes in the mitogenome of different deep-sea caridean lineages. Nine residues located in cox1, cox3, atp6, nad1, nad2, nad4, nad5, nad6 and cytb were determined to have undergone positive selection. Mitogenome of different deep-sea lineages experienced different positive selection, and the lineage represented by Alvinocarididae living in deep-sea hydrothermal vents experienced the strongest positive selection. This study provides valuable insights into the adaptive evolution of deep-sea shrimps at the mitochondrial, highlighting the mitogenomic strategy that contribute to their unique adaptations in the deep-sea environment.

Systematics of deep-sea starfish order Brisingida (Echinodermata: Asteroidea), with a revised classification and assessments of morphological characters

Brisingida Fisher 1928 is one of the seven currently recognised starfish orders, and one of the least known taxa as being exclusive deep-sea inhabitants. Modern deep-sea expeditions revealed their common occurrences in various deep-sea settings including seamounts, basins and hydrothermal vent peripheral, underlining the necessity of clarifying their global diversity and phylogeny. In this study, we present a comprehensive molecular phylogeny of Brisingida which encompasses the highest taxonomic diversity to date. DNA sequences (COI, 16S, 12S and 28S) were obtained from 225 specimens collected in the global ocean, identified as 58 species spanning 15 of the 17 extant genera. Phylogenetic relationship was inferred using both maximum likelihood and Bayesian inference methods, revealing polyphyletic families and genera and indicating nonnegligible bias in prior morphology-based systematics. Based on the new phylogeny, a novel classification of the order, consisting of 5 families and 17 genera, is proposed. Families Odinellidae, Brisingasteridae and Novodiniidae (sensu Clark and Mah, 2001) were resurrected to encompass the genera Odinella, Brisingaster and Novodinia. Brisingidae and Freyellidae were revised to include 11 and 3 genera, respectively. A new genus and species, two new subgenera and seven new combinations are described and a key to each genus and family is provided. Transformations of morphological traits were evaluated under the present phylogenetic hypothesis. A series of paedomorphic characters were found in many genera and species, which led to a high degree of homoplasy across phylogenetically distant groups. Our results provide new insights in the phylogeny and ontogeny of the order, and highlight the necessity to evaluate character convergence under sound phylogenetic hypothesis.

The Complete Mitochondrial Genomes of Two Rock Scallops (Bivalvia: Spondylidae) Indicate Extensive Gene Rearrangements and Adaptive Evolution Compared with Pectinidae

Different from the diverse family Pectinidae, the Spondylidae is a small group with a single genus that shares the sedentary life habit of cementing themselves to the substrate. However, little information related to the genetic diversity of Spondylidae has been reported. In the present study, the complete mitochondrial genomes of Spondylus versicolor and S. spinosus were sequenced and compared with those of pectinids. The mtDNA of S. versicolor and S. spinosus show similar patterns with respect to genome size, AT content, AT skew, GC skew, and codon usage, and their mitogenomic sizes are longer than most pectinid species. The mtDNA of S. spinosus is 27,566 bp in length, encoding 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes, while an additional tRNA-Met was found in the mtDNA of S. versicolor, which is 28,600 bp in length. The monophylies of Spondylidae and Pectinidae were well supported, but the internal relationships within Pectinidae remain unresolved due to the paraphyly of the genus Mimachlamy and the controversial position of the tribe Aequipectinini. The gene orders of S. versicolor and S. spinosus are almost identical but differ greatly from species of the Pectinidae, indicating extensive gene rearrangements compared with Pectinidae. Positive selection analysis revealed evidence of adaptive evolution in the branch of Spondylidae. The present study could provide important information with which to understand the evolutionary progress of the diverse and economically significant marine bivalve Pectinoidea.

A computational study on mitogenome-encoded proteins of Pavo cristatus and Pavo muticus identifies key genetic variations with functional implications

BACKGROUND

The Pavo cristatus population, native to the Indian subcontinent, is thriving well in India. However, the Pavo muticus population, native to the tropical forests of Southeast Asia, has reduced drastically and has been categorised as an endangered group. To understand the probable genetic factors associated with the decline of P. muticus, we compared the mitogenome-encoded proteins (13 proteins) between these two species.

RESULTS

Our data revealed that the most frequent variant between these two species was mtND1, which had an alteration in 9.57% residues, followed by mtND5 and mtATP6. We extended our study on the rest of the proteins and observed that cytochrome c oxidase subunits 1, 2, and 3 do not have any change. The 3-dimensional structure of all 13 proteins was modeled using the Phyre2 programme. Our data show that most of the proteins are alpha helical, and the variations observed in P. muticus reside on the surface of the respective proteins. The effect of variation on protein function was also predicted, and our results show that amino acid substitution in mtND1 at 14 sites could be deleterious. Similarly, destabilising changes were observed in mtND1, 2, 3, 4, 5, and 6 and mtATP6-8 due to amino acid substitution in P. muticus. Furthermore, protein disorder scores were considerably altered in mtND1, 2, and 5 of P. muticus.

CONCLUSIONS

The results presented here strongly suggest that variations in mitogenome-encoded proteins of P. cristatus and P. muticus may alter their structure and functions. Subsequently, these variations could alter energy production and may correlate with the decline in the population of P. muticus.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Strategy of micro-environmental adaptation to cold seep among different brittle stars’ colonization

Diffusing fluid from methane seepage in cold seep field creates zones with physicochemical gradients and divergent ecosystems like the mussel beds and clam beds. Three species of brittle stars (Ophiuroidea) were discovered in the Haima cold seep fields, of which Ophiophthalmus serratus and Histampica haimaensis were found on top of or within mussel beds and clam beds, whereas Amphiura sp. was only collected from muds in the clam bed assemblage. Here, we evaluated the genetic signatures of micro-environmental adaptation of brittle stars to cold seep through the comparison of mitogenomes. This study provided two complete mitogenome sequences of O. serratus and Amphiura sp. and compared with those of H. haimaensis and other non-seep species. We found that the split events of the seep and non-seep species were as ancient as the Cretaceous period (∼148–98 Mya). O. serratus and H. haimaensis display rapid residue mutation and mitogenome rearrangements compared to their shallow or deep-sea relatives, in contrast, Amphiura sp. only show medium, regardless of nucleotide mutation rate or mitogenome rearrangement, which may correlate with their adaptation to one or two micro-ecosystems. Furthermore, we identified 10 positively selected residues in ND4 in the Amphiura sp. lineage, suggesting important roles of the dehydrogenase complex in Amphiura sp. adaptive to the cold seep environment. Our results shed light on the different evolutionary strategies during colonization in different micro-environments.

The First Complete Mitochondrial Genome of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) and Phylogenetic Relationships within Infraorder Brachyura

Characterizing the complete mitochondrial genome (mitogenome) of an organism is useful for genomic studies in taxonomy and evolution. The mitogenomic characteristics of Eucrate crenata (Decapoda: Brachyura: Goneplacidae) have never been studied. The present study decodes the first mitogenome of E. crenata by high-throughput sequencing (HTS). The length of the mitogenome is 15,597 bp, and it contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL), and 22 transfer RNA genes. There are 14 and 23 genes observed on the heavy and light strands, respectively. E. crenata possesses a trnH-cac translocation, with the trnH-cac shifted between trnE-gaa and trnF-ttc instead of the usual location between nad5 and nad4 in decapods. Phylogenetic analyses based on the current dataset of 33 Brachyuran mitogenomes indicate that E. crenata. is closely related to Ashtoret lunaris of Matutidae. The similar codon usage and rearrangements in the two species provide evidence for their close phylogenetic relationship. Positive selection analysis showed that one residue located in cox1 was identified as a positively selected site with high BEB value (>95%), indicating that this gene was under positive selection pressure. This study is the first complete mitogenome record for the family Goneplacidae, and the results obtained may improve the understanding of the phylogeny of Goneplacidae in Brachyura.

Mitogenomics provides new insights into the phylogenetic relationships and evolutionary history of deep-sea sea stars (Asteroidea)

The deep sea (> 200 m) is considered as the largest and most remote biome, which characterized by low temperatures, low oxygen level, scarce food, constant darkness, and high hydrostatic pressure. The sea stars (class Asteroidea) are ecologically important and diverse echinoderms in all of the world’s oceans, occurring from the intertidal to the abyssal zone (to about 6000 m). To date, the phylogeny of the sea stars and the relationships of deep-sea and shallow water groups have not yet been fully resolved. Here, we recovered five mitochondrial genomes of deep-sea asteroids. The A+T content of the mtDNA in deep-sea asteroids were significantly higher than that of the shallow-water groups. The gene orders of the five new mitogenomes were identical to that of other asteroids. The phylogenetic analysis showed that the orders Valvatida, Paxillosida, Forcipulatida are paraphyletic. Velatida was the sister order of all the others and then the cladeValvatida-Spinulosida-Paxillosida-Notomyotida versus Forcipulatida-Brisingida. Deep-sea asteroids were nested in different lineages, instead of a well-supported clade. The tropical Western Pacific was suggested as the original area of asteroids, and the temperate water was initially colonized with asteroids by the migration events from the tropical and cold water. The time-calibrated phylogeny showed that Asteroidea originated during Devonian-Carboniferous boundary and the major lineages of Asteroidea originated during Permian–Triassic boundary. The divergence between the deep-sea and shallow-water asteroids coincided approximately with the Triassic-Jurassic extinction. Total 29 positively selected sites were detected in fifteen mitochondrial genes of five deep-sea lineages, implying a link between deep-sea adaption and mitochondrial molecular biology in asteroids.

The complete mitogenome of Phymorhynchus sp. (Neogastropoda, Conoidea, Raphitomidae) provides insights into the deep-sea adaptive evolution of Conoidea

The deep-sea environment is characterized by darkness, hypoxia, and high hydrostatic pressure. Mitochondria play a vital role in energy metabolism; thus, they may endure the selection process during the adaptive evolution of deep-sea organisms. In the present study, the mitogenome of Phymorhynchus sp. from the Haima methane seep was completely assembled and characterized. This mitogenome is 16,681 bp in length and contains 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of most sequenced conoidean gastropods. Some special elements, such as tandem repeat sequences and AT-rich sequences, which are involved in the regulation of the replication and transcription of the mitogenome, were observed in the control region. Phylogenetic analysis revealed that Conoidea is divided into two separate clades with high nodal support. Positive selection analysis revealed evidence of adaptive changes in the mitogenomes of deep-sea conoidean gastropods. Eight residues located in atp6, cox1, cytb, nad1, nad4, and nad5 were determined to have undergone positive selection. This study explores the adaptive evolution of deep-sea conoidean gastropods and provides valuable clues at the mitochondrial level regarding the exceptional adaptive ability of organisms in deep-sea environments.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

The complete mitochondrial genome of a cold seep gastropod Phymorhynchus buccinoides (Neogastropoda: Conoidea: Raphitomidae)

Phymorhynchus is a genus of deep-sea snails that are most distributed in hydrothermal vent or cold seep environments. In this study, we presented the complete mitochondrial genome of P. buccinoides, a cold seep snail from the South China Sea. It is the first mitochondrial genome of a cold seep member of the superfamily Conoidea. The mitochondrial genome is 15,764 bp in length, and contains 13 protein-coding genes (PCGs), 2 rRNA genes, and 22 tRNA genes. These genes are encoded on the positive strand, except for 8 tRNA genes that are encoded on the negative strand. The start codon ATG and 3 types of stop codons, TAA, TAG and the truncated termination codon T, are used in the 13 PCGs. All 13 PCGs in the 26 species of Conoidea share the same gene order, while several tRNA genes have been translocated. Phylogenetic analysis revealed that P. buccinoides clustered with Typhlosyrinx sp., Eubela sp., and Phymorhynchus sp., forming the Raphitomidae clade, with high support values. Positive selection analysis showed that a residue located in atp6 (18 S) was identified as the positively selected site with high posterior probabilities, suggesting potential adaption to the cold seep environment. Overall, our data will provide a useful resource on the evolutionary adaptation of cold seep snails for future studies.

The complete mitochondrial genome of Calyptogena marissinica (Heterodonta: Veneroida: Vesicomyidae): Insight into the deep-sea adaptive evolution of vesicomyids

The deep-sea chemosynthetic environment is one of the most extreme environments on the Earth, with low oxygen, high hydrostatic pressure and high levels of toxic substances. Species of the family Vesicomyidae are among the dominant chemosymbiotic bivalves found in this harsh habitat. Mitochondria play a vital role in oxygen usage and energy metabolism; thus, they may be under selection during the adaptive evolution of deep-sea vesicomyids. In this study, the mitochondrial genome (mitogenome) of the vesicomyid bivalve Calyptogena marissinica was sequenced with Illumina sequencing. The mitogenome of C. marissinica is 17,374 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes (rrnS and rrnL) and 22 transfer RNA genes. All of these genes are encoded on the heavy strand. Some special elements, such as tandem repeat sequences, “G(A)_nT” motifs and AT-rich sequences, were observed in the control region of the C. marissinica mitogenome, which is involved in the regulation of replication and transcription of the mitogenome and may be helpful in adjusting the mitochondrial energy metabolism of organisms to adapt to the deep-sea chemosynthetic environment. The gene arrangement of protein-coding genes was identical to that of other sequenced vesicomyids. Phylogenetic analyses clustered C. marissinica with previously reported vesicomyid bivalves with high support values. Positive selection analysis revealed evidence of adaptive change in the mitogenome of Vesicomyidae. Ten potentially important adaptive residues were identified, which were located in cox1, cox3, cob, nad2, nad4 and nad5. Overall, this study sheds light on the mitogenomic adaptation of vesicomyid bivalves that inhabit the deep-sea chemosynthetic environment.