Complete mitochondrial genome of Ophichthus brevicaudatus reveals novel gene order and phylogenetic relationships of Anguilliformes

Low depth sequencing reveals the critically endangered Batagur kachuga (Red-crowned roofed turtle) mitochondrial genome and its evolutionary implications

The Batagur kachuga (B. kachuga), commonly known as the Red-crowned roofed turtle, is a critically endangered species native to India and its neighboring countries like Bangladesh, and Nepal. The present study is the first report of the complete mitochondrial genome of B. kachuga (16,517 bp) construed via the next-generation sequencing (NGS) approach from eggshell DNA. There are 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs), and one putative control region (CR/D-loop) in the mitogenome. The CR region from the current study reveals conserved TAS, CD, and CSB domains and two AT-rich tandem repeat regions. Most genes are encoded in the heavy strand except the NADH dehydrogenase subunit 6 (ND6) gene and seven tRNA genes. Most PCGs start with the initiation codon ATG, except the COI (Cytochrome Oxidase Subunit-I) gene, which starts with the GTG codon. The present investigation also predicts the distinctive cloverleaf structures of tRNAs except for tRNA-Ser1 and tRNA-Ser-2, which lack a DHU arm. The comparative analysis of Ka/Ks with other 33 species from Order Testudines, in relation to B. kachuga, revealed negative selection in most PCGs, indicating a process of preservation and purification that aids in eliminating undesirable or detrimental substitutes. Phylogenetic analysis of this species has been analysed using the complete mitogenome of 33 turtle species. The maximum likelihood phylogenetic tree strongly supports each family in different clades and also reveals a close relationship between the Pangashura and Batagur genera. Our study suggests the generation of genome-wide molecular data, in terms of mitogenomes, SNPs, and SSRs, is needed to improve the understanding of this species and their phylogenetics and evolutionary relationships, which will help to improve the conservation efforts of this species.

Daphnia sp. (Branchiopoda: Cladocera) Mitochondrial Genome Gene Rearrangement and Phylogenetic Position Within Branchiopoda

In high-altitude (4500 m) freshwater lakes, Daphnia is the apex species and the dominant zooplankton. It frequently dwells in the same lake as the Gammarid. Branchiopoda, a class of Arthropoda, Crustacea, is a relatively primitive group in the subphylum Crustacea, which originated in the Cambrian period of the Paleozoic. The complete mitogenome sequence of Daphnia sp. (Branchiopoda: Cladocera) was sequenced and annotated in this study and deposited in GenBank. The sequence structure of this species was studied by comparing the original sequences with BLAST. In addition, we have also researched the mechanisms of their mitochondrial gene rearrangement by establishing a model. We have used the Bayesian inference [BI] and maximum likelihood [ML] methods to proceed with phylogenetic analysis inference, which generates identical phylogenetic topology that reveals the phylogenetic state of Daphnia. The complete mitogenome of Daphnia sp. shows that it was 15,254 bp in length and included two control regions (CRs) and 37 genes (13 protein-coding genes, 22 tRNAs and two ribosomal RNAs [16S and 12S]). In addition to tRNA-Ser (GCT), other tRNAs have a typical cloverleaf secondary structure. Meanwhile, the mitogenome of Daphnia sp. was clearly rearranged when compared to the mitogenome of typical Daphnia. In a word, we report a newly sequenced mitogenome of Daphnia sp. with a unique rearrangement phenomenon. These results will be helpful for further phylogenetic research and provide a foundation for future studies on the characteristics of the mitochondrial gene arrangement process in Daphnia.

Mitochondrial Genome Characteristics Reveal Evolution of Acanthopsetta nadeshnyi (Jordan and Starks, 1904) and Phylogenetic Relationships

In the present study, the mitochondrial genomic characteristics of Acanthopsetta nadeshnyi have been reported and have depicted the phylogenetic relationship among Pleuronectidae. Combined with a comparative analysis of 13 PCGs, the TN93 model was used to review the neutral evolution and habitat evolution catalysis of the mitogenome to verify the distancing and purification selectivity of the mitogenome in Pleuronectidae. At the same time, a species differentiation and classification model based on mitogenome analysis data was established. This study is expected to provide a new perspective on the phylogenetic relationship and taxonomic status of A. nadeshnyi and lay a foundation for further exploration of environmental and biological evolutionary mechanisms.

Sequence comparison of the mitochondrial genomes of Plesionika species (Caridea: Pandalidae), gene rearrangement and phylogenetic relationships of Caridea

Background

Despite the Caridean shrimps' vast species richness and ecological diversity, controversies persist in their molecular classification. Within Caridea, the Pandalidae family exemplifies significant taxonomic diversity. As of June 25, 2023, GenBank hosts only nine complete mitochondrial genomes (mitogenomes) for this family. The Plesionika genus within Pandalidae is recognized as polyphyletic. To improve our understanding of the mitogenome evolution and phylogenetic relationships of Caridea, this study introduces three novel mitogenome sequences from the Plesionika genus: P.  ortmanni, P. izumiae and P. lophotes.

Methods

The complete mitochondrial genomes of three Plesionika species were sequenced utilizing Illumina's next-generation sequencing (NGS) technology. After assembling and annotating the mitogenomes, we conducted structural analyses to examine circular maps, sequence structure characteristics, base composition, amino acid content, and synonymous codon usage frequency. Additionally, phylogenetic analysis was performed by integrating existing mitogenome sequences of true shrimp available in GenBank.

Results

The complete mitogenomes of the three Plesionika species encompass 37 canonical genes, comprising 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (CR). The lengths of these mitogenomes are as follows: 15,908 bp for P. ortmanni, 16,074 bp for P. izumiae and 15,933 bp for P. lophotes. Our analyses extended to their genomic features and structural functions, detailing base composition, gene arrangement, and codon usage. Additionally, we performed selection pressure analysis on the PCGs of all Pandalidae species available in Genbank, indicating evolutionary purification selection acted on the PCGs across Pandalidae species. Compared with the ancestral Caridea, translocation of two tRNA genes, i.e., trnP or trnT, were found in the two newly sequenced Plesionika species-P. izumiae and P. lophotes. We constructed a phylogenetic tree of Caridea using the sequences of 13 PCGs in mitogenomes. The results revealed that family Pandalidae exhibited robust monophyly, while genus Plesionika appeared to be a polyphyletic group.

Conclusions

Gene rearrangements within the Pandalidae family were observed for the first time. Furthermore, a significant correlation was discovered between phylogenetics of the Caridea clade and arrangement of mitochondrial genes. Our findings offer a detailed exploration of Plesionika mitogenomes, laying a crucial groundwork for subsequent investigations into genetic diversity, phylogenetic evolution, and selective breeding within this genus.

Complete Mitochondrial Genome of Four Peristediidae Fish Species: Genome Characterization and Phylogenetic Analysis

The systematic revision of the family Peristediidae remains an unresolved issue due to their diverse and unique morphology. Despite the popularity of using mitochondrial genome research to comprehensively understand phylogenetic relationships in fish, genetic data for peristediid fish need to be included. Therefore, this study aims to investigate the mitochondrial genomic characteristics and intra-family phylogenetic relationships of Peristediidae by utilizing mitochondrial genome analysis. Therefore, this study aims to investigate the phylogenetic relationship of Peristediidae by utilizing mitochondrial genome analysis. The mitochondrial genome of four species of Peristediidae (Peristedion liorhynchus, Satyrichthys welchi, Satyrichthys rieffeli, and Scalicus amiscus) collected in the East China Sea was studied. The mitochondrial gene sequence lengths of four fish species were 16,533 bp, 16,526 bp, 16,527 bp, and 16,526 bp, respectively. They had the same mitochondrial structure and were all composed of 37 genes and one control region. Most PCGs used ATG as the start codon, and a few used GTG as the start codon. An incomplete stop codon (TA/T) occurred. The AT-skew and GC-skew values of 13 PCGs from four species were negative, and the GC-skew amplitude was greater than that of AT-skew. All cases of D-arm were found in tRNA-Ser (GCT). The Ka/Ks ratio analysis indicated that 13 PCGs were suffering purifying selection. Based on 12 PCGs (excluding ND6) sequences, a phylogenetic tree was constructed using Bayesian inference (BI) and maximum likelihood (ML) methods, providing a further supplement to the scientific classification of Peristediidae fish. According to the results of divergence time, the four species of fish had apparent divergence in the Early Cenozoic, which indicates that the geological events at that time caused the climax of species divergence and evolution.

Comparative Mitogenomic Analysis of Two Snake Eels Reveals Irregular Gene Rearrangement and Phylogenetic Implications of Ophichthidae

The family Ophichthidae has the largest number and the most various species (about 359 valid species) in the order Anguilliformes worldwide. Both morphological and molecular characteristics have been used to assess their taxonomic status. However, due to the ambiguous morphological features, molecular data such as mitochondrial DNA sequences have been implemented for the correct identification and classification of these fishes. In this study, the gene arrangement and structure characteristics of two Ophichthidae mitochondrial genomes were investigated for the first time. The total mitogenome lengths of O. evermanni and O. erabo were 17,759 bp and 17,856 bp, respectively. Comparing with the ancestral mitochondrial gene order, the irregular gene rearrangement happened between ND6 and tRNA-Pro (P) genes with another similar control region emerging between tRNA-Thr (T) and ND6 genes, which could be explained by the tandem duplication and random loss (TDRL) model appropriately. ML phylogenetic tree demonstrated that the family Ophichthidae was monophyletic origin, but genus Ophichthus might be polyphyletic because of the confused cluster relationships among different species.

Mitogenomic Features and Evolution of the Nile River Dominant Tilapiine Species (Perciformes: Cichlidae)

To better understand the diversity and evolution of cichlids, we sequenced, assembled, and annotated the complete mitochondrial genomes of three Nile tilapiine species (Coptodon zillii, Oreochromis niloticus, and Sarotherodon galilaeus) dominating the Nile River waters. Our results showed that the general mitogenomic features were conserved among the Nile tilapiine species. The genome length ranged from 16,436 to 16,631 bp and a total of 37 genes were identified (two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), 13 protein-coding genes (PCGs), and 1 control region). The ND6 was the only CDS that presented a negative AT skew and a positive GC skew. The most extended repeat sequences were in the D-loop followed by the pseudogenes (trnS^GCU). The ND5 showed relatively high substitution rates whereas ATP8 had the lowest substitution rate. The codon usage bias displayed a greater quantity of NNA and NNC at the third position and anti-bias against NNG. The phylogenetic relationship based on the complete mitogenomes and CDS was able to differentiate the three species as previously reported. This study provides new insight into the evolutionary connections between various subfamilies within cichlids while providing new molecular data that can be applied to discriminate between Nile tilapiine species and their populations.

Characterization of the Complete Mitochondrial Genome of the Spotted Catfish Arius maculatus (Thunberg, 1792) and Its Phylogenetic Implications

The spotted catfish, Arius maculatus (Siluriformes), is an important economical aquaculture species inhabiting the Indian Ocean, as well as the western Pacific Ocean. The bioinformatics data in previous studies about the phylogenetic reconstruction of Siluriformes were insufficient and incomplete. In the present study, we presented a newly sequenced A. maculatus mitochondrial genome (mtDNA). The A. maculatus mtDNA was 16,710 bp in length and contained two ribosomal RNA (rRNA) genes, thirteen protein-coding genes (PCGs), twenty-two transfer RNA (tRNA) genes, and one D-loop region. The composition and order of these above genes were similar to those found in most other vertebrates. The relative synonymous codon usage (RSCU) of the 13 PCGs in A. maculatus mtDNA was consistent with that of PCGs in other published Siluriformes mtDNA. Furthermore, the average non-synonymous/synonymous mutation ratio (Ka/Ks) analysis, based on the 13 PCGs of the four Ariidae species, showed a strong purifying selection. Additionally, phylogenetic analysis, according to 13 concatenated PCG nucleotide and amino acid datasets, showed that A. maculatus and Netuma thalassina (Netuma), Occidentarius platypogon (Occidentarius), and Bagre panamensis (Bagre) were clustered as sister clade. The complete mtDNA of A. maculatus provides a helpful dataset for research on the population structure and genetic diversity of Ariidae species.

Novel Gene Rearrangement in the Mitochondrial Genome of Three Garra and Insights Into the Phylogenetic Relationships of Labeoninae

Complete mitochondrial genomes (mitogenomes) can provide valuable information for phylogenetic relationships, gene rearrangement, and molecular evolution. Here, we report the mitochondrial whole genomes of three Garra species and explore the mechanisms of rearrangements that occur in their mitochondrial genomes. The lengths of the mitogenomes’ sequences of Garra dengba, Garra tibetana, and Garra yajiangensis were 16,876, 16,861, and 16,835, respectively. They contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genomes of three Garra species were rearranged compared to other fish mitochondrial genomes. The tRNA-Thr, tRNA-Pro and CR (T-P-CR) genes undergo replication followed by random loss of the tRNA-Thr and tRNA-Pro genes to form tRNA-Thr, CR1, tRNA-Pro and CR2 (T-CR-P-CR). Tandem duplication and random loss best explain this mitochondrial gene rearrangement. These results provide a foundation for future characterization of the mitochondrial gene arrangement of Labeoninae and further phylogenetic studies.

Novel duplication remnant in the first complete mitogenome of Hemitriakis japanica and the unique phylogenetic position of family Triakidae

In this study, we firstly determined the complete mitogenome of the Japanese topeshark (Hemitriakis japonica), which belong to the family Triakidae and was assessed as Endangered A2d on the IUCN Red List in 2021. The mitogenome is 17,301 bp long, has a high AT content (60.0%), and contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, a control region and specially a 594 bp-long non-coding region between Cytb gene and tRNA-Thr gene. The novel non-coding region share high sequence similarity with segments of the former and latter genes, so it was recognized as a duplication remnant. In addition, the Cytb gene and tRNA-Thr gene tandemly duplicated twice while accompanied by being deleted once at least. This is the first report of mitogenomic gene-arrangement in Triakidae. The phylogenetic trees were constructed using Bayesian inference (BI) and maximum likelihood (ML) methods based on the mitogenomic data of 51 shark species and two outgroups. In summary, basing on a novel type of gene rearrangements in houndshark mitogenome, the possibly rearranged process was analyzed and contributed further insight of shark mitogenomes evolution and phylogeny.

The first mitochondrial genome of the genus Exhippolysmata (Decapoda: Caridea: Lysmatidae), with gene rearrangements and phylogenetic associations in Caridea

The complete mitochondrial genome (mitogenome) of animals can provide useful information for evolutionary and phylogenetic analyses. The mitogenome of the genus Exhippolysmata (i.e., Exhippolysmata ensirostris) was sequenced and annotated for the first time, its phylogenetic relationship with selected members from the infraorder Caridea was investigated. The 16,350 bp mitogenome contains the entire set of 37 common genes. The mitogenome composition was highly A + T biased at 64.43% with positive AT skew (0.009) and negative GC skew (- 0.199). All tRNA genes in the E. ensirostris mitogenome had a typical cloverleaf secondary structure, except for trnS1 (AGN), which appeared to lack the dihydrouridine arm. The gene order in the E. ensirostris mitogenome was rearranged compared with those of ancestral decapod taxa, the gene order of trnL2-cox2 changed to cox2-trnL2. The tandem duplication-random loss model is the most likely mechanism for the observed gene rearrangement of E. ensirostris. The ML and BI phylogenetic analyses place all Caridea species into one group with strong bootstrap support. The family Lysmatidae is most closely related to Alpheidae and Palaemonidae. These results will help to better understand the gene rearrangements and evolutionary position of E. ensirostris and lay a foundation for further phylogenetic studies of Caridea.

Novel gene rearrangement in the mitochondrial genome of Muraenesox cinereus and the phylogenetic relationship of Anguilliformes

The structure and gene sequence of the fish mitochondrial genome are generally considered to be conservative. However, two types of gene arrangements are found in the mitochondrial genome of Anguilliformes. In this paper, we report a complete mitogenome of Muraenesox cinereus (Anguilliformes: Muraenesocidae) with rearrangement phenomenon. The total length of the M. cinereus mitogenome was 17,673 bp, and it contained 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNA genes, and two identical control regions (CRs). The mitochondrial genome of M. cinereus was obviously rearranged compared with the mitochondria of typical vertebrates. The genes ND6 and the conjoint trnE were translocated to the location between trnT and trnP, and one of the duplicated CR was translocated to the upstream of the ND6. The tandem duplication and random loss is most suitable for explaining this mitochondrial gene rearrangement. The Anguilliformes phylogenetic tree constructed based on the whole mitochondrial genome well supports Congridae non-monophyly. These results provide a basis for the future Anguilliformes mitochondrial gene arrangement characteristics and further phylogenetic research.

Complete Mitogenomes of Three Carangidae (Perciformes) Fishes: Genome Description and Phylogenetic Considerations

Carangidae are ecologically and economically important marine fish. The complete mitogenomes of three Carangidae species (Alectis indicus, Decapterus tabl, and Alepes djedaba) were sequenced, characterized, and compared with 29 other species of the family Carangidae in this study. The length of the three mitogenomes ranged from 16,530 to 16,610 bp, and the structures included 2 rRNA genes (12S rRNA and 16S rRNA), 1 control region (a non-coding region), 13 protein-coding genes, and 22 tRNA genes. Among the 22 tRNA genes, only tRNA-Ser (GCT) was not folded into a typical cloverleaf secondary structure and had no recognizable DHU stem. The full-length sequences and protein-coding genes (PCGs) of the mitogenomes of the three species all had obvious AT biases. The majority of the AT-skew and GC-skew values of the PCGs among the three species were negative, demonstrating bases T and C were more plentiful than A and G. Analyses of Ka/Ks and overall p-genetic distance demonstrated that ATP8 showed the highest evolutionary rate and COXI/COXII were the most conserved genes in the three species. The phylogenetic tree based on PCGs sequences of mitogenomes using maximum likelihood and Bayesian inference analyses showed that three clades were divided corresponding to the subfamilies Caranginae, Naucratinae, and Trachinotinae. The monophyly of each superfamily was generally well supported. The divergence time analyses showed that Carangidae evolved during three geological periods, the Cretaceous, Paleogene, and Neogene. A. indicus began to differentiate from other species about 27.20 million years ago (Mya) in the early Miocene, while D. tabl (21.25 Mya) and A. djedaba (14.67 Mya) differentiated in the middle Oligocene.

Characterization of the complete mitochondrial DNA sequence of the Ariosoma meeki (Anguilliformes, Congridae)

In this study, the complete mitochondrial genome of Ariosoma meeki was sequenced, assembled and annotated. The circular genome is 16,154 bp in length with nucleotide composition is 28.42% A, 26.53% T, 19.65% G, and 25.40% C and contains 13 protein-coding genes (PCGs), 21 transfer RNA genes (tRNAs), 2 ribosomal RNA unit genes and a large non-coding region (putative control region). To further explore the evolution relationship of the Anguilliformes, we constructed the phylogenetic tree and found that the A. meeki had closer relationship with Ariosoma shiroanago. This study provided the valuable evidence on phylogenetic relationship of the A. meeki at the molecular level and essential resource for further study the molecular phylogenetic, biogeography and adaptive evolution of this lineage.

Comparative Analysis of Complete Mitochondrial Genomes of Three Gerres Fishes (Perciformes: Gerreidae) and Primary Exploration of Their Evolution History

Mitochondrial genome is a powerful molecule marker to explore phylogenetic relationships and reveal molecular evolution in ichthyological studies. Gerres species play significant roles in marine fishery, but its evolution has received little attention. To date, only two Gerres mitochondrial genomes were reported. In the present study, three mitogenomes of Gerres (Gerres filamentosus, Gerres erythrourus, and Gerres decacanthus) were systemically investigated. The lengths of the mitogenome sequences were 16,673, 16,728, and 16,871 bp for G. filamentosus, G. erythrourus, and G. decacanthus, respectively. Most protein-coding genes (PCGs) were initiated with the typical ATG codon and terminated with the TAA codon, and the incomplete termination codon T/TA could be detected in the three species. The majority of AT-skew and GC-skew values of the 13 PCGs among the three species were negative, and the amplitude of the GC-skew was larger than the AT-skew. The genetic distance and Ka/Ks ratio analyses indicated 13 PCGs were suffering purifying selection and the selection pressures were different from certain deep-sea fishes, were which most likely due to the difference in their living environment. The phylogenetic tree was constructed by molecular method (Bayesian Inference (BI) and maximum Likelihood (ML)), providing further supplement to the scientific classification of fish. Three Gerres species were differentiated in late Cretaceous and early Paleogene, and their evolution might link with the geological events that could change their survival environment.