Structure and variation of the mitochondrial genome of fishes

Complete Mitochondrial Genomes of Endemic Charrs of Salvelinus: Implications for Taxonomy and Species Identification

Gene fragments on which databases usually rely are not always suited for identification of endemic species of disputable origin. Clarification of the unresolved issues requires more comprehensive mitochondrial genome-based phylogenetics to enhance understanding of the relationship and delineation of species as a basis for their conservation. This study focuses on endemic charrs (Salvelinus, Salmonidae). The complete mitochondrial genomes (mitogenomes) of Longhead charr Salvelinus kronocius and Schmidt's charr Salvelinus schmidti from Lake Kronotskoye (Kamchatka Peninsula, Russia) were sequenced and compared with the mitogenomes of a putative ancestor species, Northern Dolly Varden Salvelinus malma malma. We analysed 36 mitogenomes, including six ones sequenced in this study, 26 charr sequences published in GenBank, and four outgroup sequences. We obtained additional evidence that the lacustrine charrs belong to the Beringian lineage with very low level of mitogenome divergence. Our findings demonstrate that the mitogenome sequencing serves as a valuable tool in clarifying the taxonomic status of narrow-range endemics or rare charrs and may be useful for future genetic investigations of charrs and other groups with high level of diversity toward proper assessment of the status of the controversial taxa. Expanding genetic resources is necessary for species identification and biodiversity monitoring, as well as for conservation measures aimed at controlling adverse natural and anthoropogenic impacts on fish population systems.

Comparative analysis of the complete mitochondrial genomes of Firmiana danxiaensis and F. kwangsiensis (Malvaceae), two endangered Firmiana species in China

MAIN CONCLUSION

We reported the mitogenomes of F. danxiaensis and F. kwangsiensis for the first time. Mitogenome size and structure differ significantly between them. Firmiana danxiaensis and F. kwangsiensis belong to the Firmiana genus and are distributed in the Danxia and Karst regions of southern China. Both species have been designated as endangered. Currently, the chloroplast genomes of F. danxiaensis and F. kwangsiensis have been sequenced, but the mitochondrial genome (mitogenome) of these two species has not been reported. To further understand the mitogenome characteristics, evolution, and phylogeny of F. danxiaensis and F. kwangsiensis, we assembled the mitogenomes of these two species based on a combination of Illumina and Nanopore sequencing methods. The mitogenome of F. danxiaensis exhibits a branching structure consisting of nine circular molecules with a total length of 938,890 bp, while the F. kwangsiensis has a circular structure with a length of 736,334 bp. Compared to F. kwangsiensis, F. danxiaensis has more tRNA genes, SSRs, tandem repeats, and dispersed repeats, while the codon use patterns are similar in these two species. There were 24 and 23 homologous sequences between mitogenome and chloroplast genome of F. danxiaensis and F. kwangsiensis, accounting for 0.37% and 0.49% of the mitogenome, respectively. In addition, the Ka/Ks ratio and the nucleic acid diversity analysis revealed that most of the mitochondria protein-coding genes in F. danxiaensis and F. kwangsiensis are highly conserved and may have undergone purifying selection. Furthermore, the collinear and comparative analysis showed that extensive genomic rearrangement events existed among the Malvaceae species. Lastly, a phylogenetic tree based on shared mitochondrial PCGs of 29 species revealed that F. danxiaensis and F. kwangsiensis form a sister group with high support values. Overall, the current study reports two mitogenomes (F. danxiaensis and F. kwangsiensis) in the Firmiana genus for the first time, which will help enhance comprehension of the mitogenome evolutionary patterns within Firmiana and promote the evolutionary and comparative genomic analyses within Malvaceae species.

Mitogenomic Characterization and Phylogenetic Insights of the Ornamental Sail-Fin Molly (Poecilia velifera) in Non-Native Indonesian Waters

The ornamental fish Poecilia velifera (Sail-fin molly, Poeciliidae) has spread widely to various non-native ecosystems around the world, far from its native habitat in the Yucatan Peninsula, Mexico. Despite the availability of some partial mitochondrial and nuclear genetic information, the complete mitogenomic structure and its variation remain unknown for this species, which is essential for a comprehensive genetic characterization and detailed phylogenetic investigation. This study applied next-generation sequencing to generate the de novo mitogenome of morphologically identified P. velifera from a non-native brackish water ecosystem in Banten Province, Indonesia. The resulting mitogenome was 16,627 bp in length and encompassed 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and a non-coding control region (CR). The result enhances our understanding of the genetic makeup of P. velifera compared to its congeners. Furthermore, the identified nucleotide variations within the conserved blocks of the CR region could provide insights into the functional role of this non-coding region. Bayesian phylogenetic inference using concatenated PCGs distinguished P. velifera from its congeners and showed monophyletic clustering of Poecilia in the family Poeciliidae, consistent with earlier evolutionary hypotheses. This first mitogenome of P. velifera paves the way for using multiple mitochondrial markers in species identification and understanding population structure in the near future. In addition, looking into the genetic evidence of this ornamental species in a non-native ecosystem, the study emphasizes the importance of strict quarantine regulations to protect Indonesia's native fish species.

The complete mitogenome of Amazonian Brachyplatystoma filamentosum and the evolutionary history of body size in the order Siluriformes

The order Siluriformes (catfish) is one of the largest groups of fish. Diversity in the body size among its species, which range from a few centimeters to 4 meters, makes Siluriformes an interesting group to investigate the body size evolution. Here, we present the complete mitogenome of Brachyplatystoma filamentosum (Piraíba), the largest Amazonian catfish, to explore the evolutionary history of Siluriformes and their body size dynamics. The Piraíba's mtDNA is 16,566 bp long, with a GC content of 42.21% and a D-loop of 911 bp. Phylogenetic analysis was conducted using protein-coding sequences, tRNAs, and rRNAs from mtDNA of Piraíba and 137 other Siluriformes species. Time-calibrated maximum likelihood trees estimated the origin of the order Siluriformes to be ~118.4 Ma, with the Loricarioidei suborder diversifying first, followed by Diplomystoidei and Siluroidei. The Siluroidei suborder experienced rapid expansion around 94.1 Ma. Evolutionary dynamics revealed 16 positive and 11 negative directional body size changes in Siluriformes, with no global trend toward larger or smaller sizes, and with Piraíba showing a significant size increase (5.65 times over 40.8 Ma). We discuss how biological, ecological and environmental factors could have shaped the evolution of body size in this group.

Complete mitochondrial genome of Inimicus didactylus (Pallas, 1769) and its phylogenetic analysis

Inimicus didactylus is a venomous fish belonging to the family Synanceiidae, which is closely related to the true stonefishes (Synanceiidae, Synanceia spp.). Here, we employed high-throughput sequencing technology to assemble and annotate the complete mitochondrial genome of I. didactylus. The mitochondrial genome, with a length of 16,670 base pairs, exhibits A-T bias in its base composition (58.19%) and consists of 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Phylogenetic analysis indicates a close relationship between I. didactylus and I. japonicus. This study provides a genomic resource that enhances ecology and evolution researches on Synanceiidae species.

The First Complete Mitochondrial Genomes for the Genus Dianema (Siluriformes: Callichthyidae): Dianema longibarbis and D. urostriatum

BACKGROUND/OBJECTIVES

To date, no information is available on the complete mitochondrial genome of the genus Dianema (Siluriformes: Callichthyidae), a callichthyid catfish. In this study, we report on two complete mitochondrial genome sequences of Dianema longibarbis Cope, 1872, and Dianema urostriatum Miranda Ribeiro, 1912, the only two recognized species within the genus Dianema.

METHODS

DNA sequencing was performed using the HiSeq platform to obtain their complete mitogenomes. To confirm phylogenetic distance, two phylogenetic trees were established using maximum-likelihood and Bayesian inference methods with all concatenated protein-coding sequences (PCGs) and two ribosomal RNA (rRNA) genes from the D. longibarbis and D. urostriatum mitogenomes, along with 32 mitogenomes retrieved from Siluriformes.

RESULTS

The complete mitogenomes of D. longibarbis and D. urostriatum are 16,493 and 16,495 base pairs in length, respectively. Their nucleotide compositions are 31.79% A, 27.53% T, 25.86% C, and 14.82% G for D. longibarbis, and 31.69% A, 27.04% T, 26.36% C, and 14.91% G for D. urostriatum. Both mitogenomes contain 13 PCGs, 22 transfer RNA (tRNA) genes, and two rRNA genes. Phylogenetic results based on all PCGs and two rRNAs genes confirm D. longibarbis as a sister species to D. urostriatum in the subfamily Callichthyinae.

CONCLUSIONS

In contrast to the extensive mitochondrial studies on species in the Corydoradinae, species in the Callichthyinae have been largely understudied. This study provides valuable insights into genetic diversity and evolutionary complexity by presenting the first mitochondrial genome analysis of two Dianema species, a genus within the Callichthyinae.

Mitogenome of Neolissochilus pnar, the largest cavernicolous species of Mahseer

The study of the mitogenome of Neolissochilus pnar, the world's largest cave fish, uncovered its structural features, gene content and evolutionary dynamics within mahseer. Its mitogenome is of 16,440 base pairs, resembling those of the teleost species and exhibits a high degree of conservation in genes arrangement. It comprises 37 mitochondrial genes, including 13 protein-coding genes (PCGs), 22 tRNA genes (tRNAs), 2 rRNA genes (rRNAs) and a control region. Notably, the distribution of genes on the L- and H-strands is consistent with that of the typical teleost. The study reveals the lengths and variations in PCGs in mahseer species, displaying a range from 164 to 11,404 bp. The tRNA and rRNA genes and the control region also demonstrate conservation among the species. A robust phylogenetic analysis, employing Bayesian and ASAP methods, supports the classification of N. pnar within the Neolissochilus genus and validates the taxonomic status of this species. Selection pressure analyses indicate positive selection in seven genes: COII , COIII, Cytb, ND1, ND2, ND5 and ND6. These findings suggest the dynamic nature of mitochondrial evolution in mahseer species. The purifying selection preserve essential mitochondrial functions, and additionally, the specific sites in ND5 and ND6 genes undergo episodic positive or diversifying selection, likely in response to environmental changes or selective pressures. In conclusion, this research enriches our understanding of N. pnar visa-vis other mahseers' mitogenomes, pointing to its possible mitogenome evolution to adaptation to cave environment.

The mitochondrial genomes of Iberian freshwater and diadromous fishes

The Iberian Peninsula, in southwestern Europe, is home to a distinctive freshwater fish fauna, predominantly composed of endemic species. This is a consequence of the prolonged isolation from western Europe caused by the Pyrenees, the diverse geological and climatic gradients, and the isolation of river basins. Freshwater and diadromous fishes have diversified in the Iberian Peninsula and include 72 currently recognized native species, 50 of which are endemic to the region. Habitat loss and degradation, the introduction of invasive species, and climate change have placed Iberian freshwater and diadromous fishes among the most threatened groups of vertebrates, with some species on the brink of extinction. Here, we present 60 new complete mitochondrial genome assemblies out of the 109 freshwater and diadromous fish species found in the Iberian Peninsula, including the mitogenomes of 37 endemics. These resources are crucial for characterising the mitochondrial evolution of species, reconstructing phylogeny and paleogeography, advancing species identification, delineation, and monitoring, and ultimately supporting conservation planning.

Comprehensive analysis of 111 Pleuronectiformes mitochondrial genomes: insights into structure, conservation, variation and evolution

BACKGROUND

Pleuronectiformes, also known as flatfish, are important model and economic animals. However, a comprehensive genome survey of their important organelles, mitochondria, has been limited. Therefore, we aim to analyze the genomic structure, codon preference, nucleotide diversity, selective pressure and repeat sequences, as well as reconstruct the phylogenetic relationship using the mitochondrial genomes of 111 flatfish species.

RESULTS

Our analysis revealed a conserved gene content of protein-coding genes and rRNA genes, but varying numbers of tRNA genes and control regions across species. Various gene rearrangements were found in flatfish species, especially for the rearrangement of nad5-nad6-cytb block in Samaridae family, the swapping rearrangement of nad6 and cytb gene in Bothidae family, as well as the control region translocation and tRNA-Gln gene inversion in the subfamily Cynoglossinae, suggesting their unique evolutionary history and/or functional benefit. Codon usage showed obvious biases, with adenine being the most frequent nucleotide at the third codon position. Nucleotide diversity and selective pressure analysis suggested that different protein-coding genes underwent varying degrees of evolutionary pressure, with cytb and cox genes being the most conserved ones. Phylogenetic analysis using both whole mitogenome information and concatenated independently aligned protein-coding genes largely mirrored the taxonomic classification of the species, but showed different phylogeny. The identification of simple sequence repeats and various long repetitive sequences provided additional complexity of genome organization and offered markers for evolutionary studies and breeding practices.

CONCLUSIONS

This study represents a significant step forward in our comprehension of the flatfish mitochondrial genomes, providing valuable insights into the structure, conservation and variation within flatfish mitogenomes, with implications for understanding their evolutionary history, functional genomics and fisheries management. Future research can delve deeper into conservation biology, evolutionary biology and functional usages of variations.

Characterization of Complete Mitochondrial Genome and Phylogeny of Three Echeneidae Species

Species of the family Echeneidae are renowned for their capacity to adhere to various hosts using a sucking disc. This study aimed to examine the mitochondrial genome characteristics of three fish species (Echeneis naucrates, Remora albescens, and Remora remora) within the family Echeneidae and determine their phylogenetic relationships. The findings revealed that the mitochondrial genome lengths of the three species were 16,611 bp, 16,648 bp, and 16,623 bp, respectively, containing 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a D-loop region. Most PCGs utilized ATG as the initiation codon, while only cox I used the GTG as the initiation codon. Additionally, seven genes employed incomplete termination codons (T and TA). The majority of PCGs in the three species displayed negative AT-skew and GC-skew values, with the GC-skew amplitude being greater than the AT-skew. The Ka/Ks ratios of the 13 PCGs did not exceed 1, demonstrating these species had been subjected to purification selection. Furthermore, only tRNA-Ser (GCT) lacked the D arm, while other tRNAs exhibited a typical cloverleaf secondary structure. Bayesian inference (BI) and maximum likelihood (ML) methods were utilized to construct a phylogenetic tree of the three species based on the 13 PCGs. Remora remora was identified as a distinct group, while R. osteochir and R. brachyptera were classified as sister taxa. This study contributes to the mitochondrial genome database of the family Echeneidae and provides a solid foundation for further systematic classification research in this fish group.

Complete Mitochondrial Genome of Niphon spinosus (Perciformes: Niphonidae): Genome Characterization and Phylogenetic Analysis

The species Niphon spinosus (Cuvier, 1829) is the only representative of the family Niphonidae and the genus Niphon, and its taxonomic history is complicated; it is still unclear in a phylogenetic sense. In this study, we report the complete mitochondrial genome of N. spinosus (OP391482), which was determined to be 16,503 bp long with biased A + T contents (53.8%) using next-generation technology. The typical set of 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes, and one control region (D-loop) are included in the mitochondrial genome. The H-strand encoded 28 genes (14 tRNA, 2 rRNA, and 12 PCGs), and D-loop, whereas the L-strand encoded the remaining 9 genes (8 tRNA and ND6). Its nucleotide composition, gene arrangement, codon usage patterns, and tRNA secondary structures are identical with other members of the Percoidei suborder. Furthermore, we reconstructed phylogenetic trees based on the 13 PCGs. The resulting phylogenetic trees showed N. spinosus placing as a separate lineage within the family Niphonidae, its close relationship to Trachinus draco (Trachinidae), and the clustering of major subfamilies like Luciopercinae and Percinae of the Percoidei suborder. These findings will contribute to future studies on the evolutionary history, population genetics, molecular taxonomy, and phylogeny of N. spinosus and related species.

Assembly of Mitochondrial Genomes Using Nanopore Long-Read Technology in Three Sea Chubs (Teleostei: Kyphosidae)

Complete mitochondrial genomes have become markers of choice to explore phylogenetic relationships at multiple taxonomic levels and they are often assembled using whole genome short-read sequencing. Herein, using three species of sea chubs as an example, we explored the accuracy of mitochondrial chromosomes assembled using Oxford Nanopore Technology (ONT) Kit 14 R10.4.1 long reads at different sequencing depths (high, low and very low or genome skimming) by comparing them to 'gold' standard reference mitochondrial genomes assembled using Illumina NovaSeq short reads. In two species of sea chubs, Girella nigricans and Kyphosus azureus, ONT long-read assembled mitochondrial genomes at high sequencing depths (> 25× whole [nuclear] genome) were identical to their respective short-read assembled mitochondrial genomes. Not a single 'homopolymer insertion', 'homopolymer deletion', 'simple substitution', 'single insertion', 'short insertion', 'single deletion' or 'short deletion' were detected in the long-read assembled mitochondrial genomes after aligning each one of them to their short-read counterparts. In turn, in a third species, Medialuna californiensis, a 25× sequencing depth long-read assembled mitochondrial genome was 14 nucleotides longer than its short-read counterpart. The difference in total length between the latter two assemblies was due to the presence of a short motif 14 bp long that was repeated (twice) in the long read but not in the short-read assembly. Read subsampling at a sequencing depth of 1× resulted in the assembly of partial or complete mitochondrial genomes with numerous errors, including, among others, simple indels, and indels at homopolymer regions. At 3× and 5× subsampling, genomes were identical (perfect) or almost identical (quasiperfect, 99.5% over 16,500 bp) to their respective Illumina assemblies. The newly assembled mitochondrial genomes exhibit identical gene composition and organisation compared with cofamilial species and a phylomitogenomic analysis based on translated protein-coding genes suggested that the family Kyphosidae is not monophyletic. The same analysis detected possible cases of misidentification of mitochondrial genomes deposited in GenBank. This study demonstrates that perfect (complete and fully accurate) or quasiperfect (complete but with a single or a very few errors) mitochondrial genomes can be assembled at high (> 25×) and low (3-5×) but not very low (1×, genome skimming) sequencing depths using ONT long reads and the latest ONT chemistries (Kit 14 and R10.4.1 flowcells with SUP basecalling). The newly assembled and annotated mitochondrial genomes can be used as a reference in environmental DNA studies focusing on bioprospecting and biomonitoring of these and other coastal species experiencing environmental insult. Given the small size of the sequencing device and low cost, we argue that ONT technology has the potential to improve access to high-throughput sequencing technologies in low- and moderate-income countries.

The Complete Mitochondrial Genome of the Luciocyprinus langsoni (Cypriniformes: Cyprinidae): Characterization, Phylogeny, and Genetic Diversity Analysis

Background: Luciocyprinus langsoni is a species belonging to the Cyprinidae family. The objective of this study is to gain a comprehensive understanding of its evolutionary history and genetic characteristics. Methods: The complete mitochondrial genome of L. langsoni was determined using overlapping PCR. A phylogenetic analysis was conducted based on 13 protein-coding genes from 48 species. A population genetic diversity analysis using the COI gene and a selection analysis of 13 protein-coding genes were also performed. Results: The mitogenome is 16,586 base pairs long and consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and two control regions. It has a high adenine-thymine (A + T) content. The phylogenetic analysis confirms the placement of L. langsoni within the subfamily Cyprininae. The population genetic diversity analysis reveals low variability in the Hechi Longjiang population. The selection analysis shows that all 13 protein-coding genes have evolved under purifying selection with Ka/Ks ratios below 1. Conclusions: These results enhance our understanding of L. langsoni's evolutionary history and lay a genetic foundation for future studies in population genetics and phylogenetics.

Complete mitochondrial genome of Lepidocephalichthysberdmorei and its phylogenetic status within the family Cobitidae (Cypriniformes)

In this study, the complete mitochondrial genome of Lepidocephalichthysberdmorei was first determined by the primer walking sequence method. The complete mitochondrial genome was 16,574 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA (rRNA) genes, and a control region (D-loop). The gene arrangement pattern was identical to that of other teleosts. The overall base composition was 29.9% A, 28.5% T, 25.5% C, and 16.1% G, with an A+T bias of 58.4%. Furthermore, phylogenetic analyses were conducted based on 13 PCGs from the mitochondrial genomes of 18 cobitid species using with three different methods (Neighbor-joining, Maximum likelihood, and Bayesian inference). All methods consistently showed that the four species of the genus Lepidocephalichthys form a monophyletic group. This study would provide effective molecular information for the Lepidocephalichthys species as well as novel genetic marker for the study of species identification.

Morphometric and phylogenetic analysis of a commercial fish Leiognathusequula (Teleostei, Leiognathidae)

The complete mitochondrial genome and phylogenetic analysis for Leiognathusequula from the South China Sea was performed using shallow genome skimming. For accurate species identification and redescription, morphometric and meristic characters were examined and compared with previous descriptions. To facilitate the identification of species and to enable comparison with the mitochondrial genome phylogeny, molecular comparisons were conducted using three mitochondrial genes: large ribosomal RNA (16S rRNA), cytochrome c oxidase subunit 1 (COX1), and NADH dehydrogenase (ND5). The mitogenome (16 398 bp) comprised 38 mitochondrial genes, similar to most bony fishes: 13 protein-coding genes (PCGs), 2 rRNA and 22 transfer RNA genes, and 1 non-coding region. The complete mitogenome comprised 30.7% A, 24.2% T, 15.0% G, and 30.1% C. The A+T content (54.9%) was higher than the G+C content (45.1%). All PCGs started with the typical ATG codon, except COX1, which started with GTG. Seven PCGs ended with incomplete stop codons (TA or T). The Ka/Ks ratios of all PCGs were < 1, indicating purifying selection. The phylogenies of Leiognathidae, both based on the amino acid sequences of the 13 PCGs and the single genes 16S RNA and COX1, were monophyletic with high nodal support (> 75%). L.brevirostris (NC 026232) is believed to be a Nuchequula species, whereas L.ruconius (NC 057225) is not classified under Leiognathus in the NCBI database. The phylogenetic trees divided the Leiognathidae family into three clades. The mitogenome phylogeny suggested that the Leiognathidae and Chaetodontidae clades are sister groups. These findings provide important genetic data for population genetics research and a phylogenetic analysis of Leiognathidae.

Partial and complete sequence of small and large subunit ribosomal RNA genes, tRNA-Val gene in some species of family Labridae

Background

Mitochondrial genomes play a key role in molecular biology research by providing essential information about evolutionary links, population history, and genetic diversity.

Aim

The aim of this investigation was to produce a partial sequence of 12S rRNA and 16S rRNA genes, as well as a complete sequence of tRNA-Val gene in some species of family Labridae.

Methods

Five species of labrid fishes (Oxycheilinus digramma, Cheilio inermis, Epibulus insidiator, Coris aygula, and Gomphosus caeruleus) belonging to Family Labridae were collected from the Red Sea, thereafter, taken to a laboratory for morphological identification in accordance with. Using forward and reverse primers, genome DNA was amplified through polymerase chain reaction.

Results

The tRNA-Val gene's entire sequence, the 12S rRNA gene's partial sequence, and the 16S rRNA gene's partial sequence were all submitted to GenBank/NCBI with accession numbers (PP962382.1-PP962386.1). The sequences' outcomes showed that the average A + T values were higher than the C + G values.

Conclusion

The partial sequences of 12S RNA and 16S RNA, and the whole sequence of the tRNA-Val gene, were arranged so that, the 12S RNA and 16S RNA have been distinguished by the tRNA-Val gene.

Insights into the structural features and phylogenetic implications of the complete mitochondrial genome of Fasin rainbow fish (Melanotaenia fasinensis)

The Fasin rainbow fish, scientifically named Melanotaenia fasinensis, is highly prized by aquarium enthusiasts for its vibrant colors and adaptability to artificial aquatic environments. This species is endemic to the karst landscape of the Bird's Head region in Papua, Indonesia, and belongs to the family Melanotaeniidae. Discovered relatively recently in 2010, this species was designated as endangered by the International Union for Conservation of Nature (IUCN) in 2021. However, there is currently insufficient data regarding its phylogenetic position. To address this gap, our study employed next-generation sequencing (NGS) to analyze the entire mitochondrial genome of M. fasinensis. The mitochondrial genome comprises 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes, with a total length of 16,731 base pairs. The base composition of the mitogenome revealed percentages of 27.76% adenine (A), 27.34% thymine (T), 16.15% guanine (G), and 28.75% cytosine (C) residues. Our phylogenetic analysis based on sequence data indicated that all species of the Melanotaeniidae family clustered together on the same branch. Furthermore, the intergeneric and interspecific taxonomic positions were explicit and clear. Phylogenetically, Melanotaeniidae were more closely related to the family Isonidae than to the family Atherinomorus. The phylogenetic position of M. fasinensis was relatively basal within the genus Melanotaenia. This study provides valuable molecular insights for further exploration of the phylogeography and evolutionary history of M. fasinensis and other members of the genus Melanotaenia.

The complete mitochondrial genome of Chilean Jack Mackerel, Trachurus murphyi Nichols, 1920 (Teleostei, Carangidae)

The Chilean jack mackerel (Trachurus murphyi, Nichols, 1920) is a pelagic fish of high fishery interest, with a global capture of 828,358 mt in 2021. We report the complete mitochondrial genome of T. murphyi (Teleostei, Carangidae), collected in Bahia de Zapallar, Chile (32°33'02.59'' S, 71°27'55.38'' W). The size of the mitogenome is 16,561 bp (H-strand composition: 25.9% A, 16.7% C, 29.8% G, and 27.5% T). The mitogenome has the classic vertebrate gene content of 13 protein-coding, two rRNA, and 22 tRNA genes, as found in Carangidae and other Teleostei families. Phylogenetic analysis using mitochondrial genomes of 22 related species revealed that T. murphyi formed a well-supported monophyletic group with the other Trachurus species, being T. simmetricus its closest relative. Sequencing the mitochondrial genome from T. murphyi is the first step in developing traceability tools based on DNA analysis to enforce fishing quotas and trace the processed food and foodstuff containing Chilean jack mackerel following the objective of the South Pacific Regional Fisheries Management Organization (SPRFMO).

Elucidating the Mitogenomic Blueprint of Pomadasys perotaei from the Eastern Atlantic: Characterization and Matrilineal Phylogenetic Insights into Haemulid Grunts (Teleostei: Lutjaniformes)

The parrot grunt fish, Pomadasys perotaei, has a limited distribution in the Eastern Atlantic Ocean and is an important species in marine capture fisheries across several West African countries. Despite its ecological and economic significance, the mitogenomic information for this species is lacking. This study utilized next-generation sequencing to generate the de novo mitogenome of P. perotaei from Eastern Atlantic. The resulting mitogenome is 16,691 base pairs and includes 13 protein-coding genes (PCGs), 22 transfer RNAs, two ribosomal RNAs, and an AT-rich control region (CR). Most of the PCGs exhibit nonsynonymous (Ka) and synonymous (Ks) substitution rates of less than '1', indicating strong negative selection across haemulid fishes. The control region of Pomadasys species contains four conserved domains, as seen in other teleost's, with polymorphic nucleotides that can be used to study population structures through the amplification of short mitochondrial gene fragments. Additionally, Bayesian phylogenetic analysis based on PCGs revealed a non-monophyletic clustering pattern of Pomadasys within the haemulid matrilineal tree. Overall, the structural characterization and phylogenetic analysis enhance our understanding of the genetic composition and evolutionary history of Pomadasys species from the Indo-West Pacific and Eastern Atlantic Oceans.

Mitogenomic Architecture of Atlantic Emperor Lethrinus atlanticus (Actinopterygii: Spariformes): Insights into the Lineage Diversification in Atlantic Ocean

The evolutionary history of emperors, particularly in the Atlantic and Indo-West Pacific Oceans, remains largely unmapped. This study explores the maternal lineage evolution of Lethrinids by examining the complete mitogenome of Lethrinus atlanticus, which is endemic to the Eastern Atlantic Ocean. Utilizing advanced next-generation sequencing, we found that the mitogenome spans 16,789 base pairs and encompasses 37 genes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and an AT-rich control region (CR). Our analysis indicates a preference for AT base pairs in the L. atlanticus mitogenome (53.10%). Most PCGs begin with the ATG codon, except for COI, which starts with GTG. Relative synonymous codon usage reveals high frequencies for alanine, leucine, proline, serine, and threonine. The ratio of nonsynonymous to synonymous substitutions suggests strong negative selection across all PCGs in Lethrinus species. Most transfer RNAs exhibit typical cloverleaf structures, with the exception of tRNA-serine (GCT), which lacks a dihydrouracil stem. Comparative analysis of conserved sequence blocks across the CRs of three Lethrinus species shows notable differences in length and nucleotide composition. Phylogenetic analysis using concatenated PCGs clearly distinguishes all Lethrinus species, including L. atlanticus, and sheds light on the evolutionary relationships among Spariformes species. The estimated divergence time of approximately 20.67 million years between L. atlanticus and its Indo-West Pacific relatives provides insights into their historical separation and colonization during the late Oligocene. The distribution of Lethrinids may be influenced by ocean currents and ecological factors, potentially leading to their speciation across the Eastern Atlantic and Indo-West Pacific. This study enhances our understanding of the genetic diversity and phylogenetic relationships within Lethrinus species. Further exploration of other emperor fish mitogenomes and comprehensive genomic data could provide vital insights into their genetic makeup, evolutionary history, and environmental adaptability in marine ecosystems globally.

Comparative mitochondrial genomics analysis of selected species of Schizothoracinae sub family to explore the differences at mitochondrial DNA level

A comprehensive analysis of the whole mitochondrial genomes of the Schizothoracinae subfamily of the family Cyprinidae has been revealed for the first time. The species analyzed include Schizothorax niger, Schizothorax esocinus, Schizothorax labiatus and Schizothorax plagoistomus. The total mitochondrial DNA (mtDNA) length was determined to be 16585 bp, 16583 bp, 16582 bp and 16576 bp, respectively with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and 2 non-coding area genes. The combined mean base compositions of the four species were as follows: A: 29.91 % T: 25.47 % G: 17.65 % C 27.01 %. The range of the GC content is 45-44 %, respectively. All protein coding genes (PCGs) commenced with the typical ATG codon, except for the cytochrome c oxidase subunit 1 (COX1) gene with GTG. The analysis of vital amino acid biosynthesis genes (COX1, ATPase 6, ATPase 8) in four different species revealed no significant differences. All 13 PCGs had Ka/Ks ratios that were all lesser than one, demonstrating purifying selection on those molecules. These tRNA genes were predicted to fold into the typical cloverleaf secondary structures with normal base pairing and ranged in size from 66 to 75 nucleotides. Additionally, the phylogenetic tree analysis revealed that S. esocinus species that was most alike to S. labiatus. This study provides critical data for phylogenetic analysis of the Schizothoracinae subfamily, which will help to resolve taxonomic difficulties and identify evolutionary links. Detailed mtDNA data are an invaluable resource for studying genetic diversity, population structure, and gene flow. Understanding genetic makeup can help inform conservation plans, identify unique populations, and track genetic variation to ensure effective preservation.

Unraveling the Mitogenomic Characteristics and Phylogenetic Implications of Leuciscus merzbacheri (Zugmayer, 1912), an Endangered Fish in the Junggar Basin of Xinjiang, Northwest China

BACKGROUND

Leuciscus merzbacheri is a rare and endangered fish in Xinjiang, China. As a representative species of the fauna in the Junggar Basin, it is of high economic and scientific value. The genetic data are still limited, and the mitochondrial genomic characteristics remain unexplored.

METHODS

A high-throughput sequencing method was used to obtain the complete mitogenome of L. merzbacheri.

RESULTS

The full length of the circular DNA was 16,609 bp, and it consisted of 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs and 2 non-coding regions. The overall nucleotide compositions of both the mitogenome and PCGs showed an obvious AT preference with percentages of 54.20% and 53.60%, respectively. Three commonly used amino acids were Leu (16.43%), Ala (8.95%) and Thr (7.85%) in turn. All tRNAs could form the typical clover structures excluding tRNA-Ser AGY. The presumed secondary structures of two rRNAs contained several stem-loop domains, and the structure of 12S rRNA seemed to be more stable than that of 16S rRNA. Extended termination sequence regions (ETASs), central conserved regions (CSB-F, CSB-E and CSB-D), and conserved sequence regions (CSB-1, CSB-2 and CSB-3) were identified in the control region. The phylogenetic tree showed that L. merzbacheri was recovered with strong supports as a sister to the other members of the genus. The location in the outermost branch implied that it might be a relatively ancient species among its congeners.

CONCLUSIONS

This study would complement the genetic data on L. merzbacheri and contribute to a better understanding of molecular evolution in Leuciscus as well.

Endemic Radiation of African Moonfish, Selene dorsalis (Gill 1863), in the Eastern Atlantic: Mitogenomic Characterization and Phylogenetic Implications of Carangids (Teleostei: Carangiformes)

This study offers an in-depth analysis of the mitochondrial genome of Selene dorsalis (Gill 1863), a species native to the Eastern Atlantic Ocean. The circular mitochondrial DNA molecule measures 16,541 base pairs and comprises 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, two ribosomal RNA genes, and a control region (CR). The nucleotide composition exhibits a notable adenine-thymine (AT) bias, accounting for 53.13%, which aligns with other species in the Carangidae family. Most PCGs initiate with the ATG codon, with the exception of Cytochrome C oxidase subunit I, which starts with GTG. Analysis of relative synonymous codon usage reveals that leucine and serine are the most prevalent amino acids in the mitochondrial genome of S. dorsalis and its congeners (S. vomer and S. setapinnis). All tRNAs display the typical cloverleaf structure, though tRNA Serine (S1) lacks a dihydrouracil arm. Pairwise comparisons of synonymous and nonsynonymous substitutions for all PCGs yielded values below '1', indicating strong purifying selection. The CR spans 847 bp, representing 5.12% of the mitochondrial genome, and is characterized by high AT content (62.81%). It is situated between tRNA-Pro (TGG) and tRNA-Phe (GAA). The CR contains conserved sequence blocks, with CSB-1 being the longest at 22 bp and CSB-D the shortest at 18 bp. Phylogenetic analysis, using Bayesian and Maximum-likelihood trees constructed from concatenated PCGs across 72 species, successfully differentiates S. dorsalis from other carangids. This study also explores how ocean currents and gyres might influence lineage diversification and parapatric speciation of Selene species between the Atlantic and Pacific Oceans. These results highlight the importance of the mitochondrial genome in elucidating the structural organization and evolutionary dynamics of S. dorsalis and its relatives within marine ecosystems.

Long-read de novo genome assembly of Gulf toadfish (Opsanus beta)

BACKGROUND

The family Batrachoididae are a group of ecologically important teleost fishes with unique life histories, behavior, and physiology that has made them popular model organisms. Batrachoididae remain understudied in the realm of genomics, with only four reference genome assemblies available for the family, with three being highly fragmented and not up to current assembly standards. Among these is the Gulf toadfish, Opsanus beta, a model organism for serotonin physiology which has recently been bred in captivity.

RESULTS

Here we present a new, de novo genome and transcriptome assemblies for the Gulf toadfish using PacBio long read technology. The genome size of the final assembly is 2.1 gigabases, which is among the largest teleost genomes. This new assembly improves significantly upon the currently available reference for Opsanus beta with a final scaffold count of 62, of which 23 are chromosome scale, an N50 of 98,402,768, and a BUSCO completeness score of 97.3%. Annotation with ab initio and transcriptome-based methods generated 41,076 gene models. The genome is highly repetitive, with ~ 70% of the genome composed of simple repeats and transposable elements. Satellite DNA analysis identified potential telomeric and centromeric regions.

CONCLUSIONS

This improved assembly represents a valuable resource for future research using this important model organism and to teleost genomics more broadly.

The expansion of invasive species to the East: new sites of the bullheads (genus Ameiurus Rafinesque 1820) in Ukraine with morphological and genetic identification

This study confirms the extended distribution of two invasive species of the genus Ameiurus in Ukraine. Specifically, A. melas is recorded for the first time in the Southern Buh basin and A. nebulosus has expanded further eastward within the Dnipro basin. Material collected in 2019 and 2022 was identified by morphological features and confirmed by molecular genetic analysis. The most reliable morphological characters for distinguishing these two species include anal-fin membrane pigmentation (light or black), gill raker count (fewer or more than 16), and serrations on the pectoral-fin spine (well-developed along the full length or small, absent near the tip). The analysis of the cytochrome oxidase subunit I barcoding marker identified all samples from the Dnipro Basin (Tnia and Velykyi Luh localities) as A. nebulosus, while all specimens from the Vinnytsia region within the Southern Buh basin (Sotskoho and Vyshenske lakes) were attributed to A. melas. The maximum-likelihood analysis revealed clearly separated clades with high bootstrap support (>75%), strongly supporting the presence of the two separate species. This study suggests the potential for further eastward expansion of both species within Ukraine: A. nebulosus in the northern direction and A. melas in the southern direction.

Exploring the mitochondrial genomes and phylogenetic relationships of trans-Andean Bryconidae species (Actinopterygii: Ostariophysi: Characiformes)

Comparative mitogenomics and its evolutionary relationships within Bryconidae remains largely unexplored. To bridge this gap, this study assembled 15 mitogenomes from 11 Bryconidae species, including five newly sequenced. Salminus mitogenomes, exceeding 17,700 bp, exhibited the largest size, contrasting with a median size of 16,848 bp in the remaining species (Brycon and Chilobrycon). These mitogenomes encode 37 typical mitochondrial genes, including 13 protein-coding, 2 ribosomal RNA, and 22 transfer RNA genes, and exhibit the conserved gene arrangement found in most fish species. Phylogenetic relationships, based on the maximum-likelihood method, revealed that the trans-Andean species (found in northwestern South America) clustered into two main sister clades. One clade comprised the trans-Andean species from the Pacific slope, Brycon chagrensis and Chilobrycon deuterodon. The other clade grouped the trans-Andean species from the Magdalena-Cauca Basin Brycon moorei and Salminus affinis, with their respective cis-Andean congeners (found in eastern South America), with Brycon rubricauda as its sister clade. Since the current members of Brycon are split in three separated lineages, the systematic classification of Bryconidae requires further examination. This study provides novel insights into mitogenome characteristics and evolutionary pathways within Bryconidae, standing as crucial information for prospective phylogenetic and taxonomic studies, molecular ecology, and provides a valuable resource for environmental DNA applications.

Complete Mitochondrial Genome of Clupeonella abrau (Maliatsky, 1930) (Clupeiformes), an Endemic Freshwater Fish from the Lake Abrau (Russia)

The Abrau sprat (tyulka or sardelka) Clupeonella abrau (Maliatsky, 1930) is an endemic fish of the Lake Abrau (Krasnodar Krai, Russia). The full C. abrau mitochondrial genome (16 650 bp) showed a gene arrangement conserved in Clupeidae and 98.8% similarity with the mitochondrial genome of the related species Black and Caspian Sea sprat C. cultriventris from the Black Sea. The COX1 gene sequence was additionally studied in a museum specimen collected in the Lake Abrau in 1938. Variability in modern Abrau sprat COX1 gene locus was estimated at approximately 0.15%, the difference between C. abrau and C. cultriventris was 1.2%, and the difference between the museum and modern C. abrau specimens from the Lake Abrau was 0.92%. The study confirmed that the Abrau sprat is present in the fish community and is capable of reproducing in the lake. Various scenarios were proposed to explain colonization of the Lake Abrau by C. abrau.

The complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909)

Holocephali is a subclass of chondrichthyans with ample geographic distribution in marine ecosystems. Holocephalan species are organized into three families: Callorhinchidae, Chimaeridae, and Rhinochimaeridae. Despite the critical ecological and evolutionary importance, genomic information from holocephalans is still scarce, particularly from rhinochimaerids. The present study provides the first complete mitogenome of the Atlantic longnose chimaera Rhinochimaera atlantica (Holt & Byrne, 1909). The whole mitogenome was sequenced from an R. atlantica specimen, collected on the Porcupine Bank (NE Atlantic), by Illumina high-throughput sequencing. The R. atlantica mitogenome has 17,852 nucleotides with 13 protein-coding genes, 22 transfer RNA, and two ribosomal RNA genes. Nine of these genes are in the complementary strand. This mitogenome has a GC content of 41.5% and an AT content of 58.5%. The phylogenetic reconstruction provided here, using all the available complete and partial Holocephali mitogenomes, places R. atlantica in the Rhinochimaeridae family, as expected. This genomic resource will be useful in the genomic characterization of this species.

Venlafaxine exposure alters mitochondrial respiration and mitomiR abundance in zebrafish brains

Wastewater treatment plant (WWTP) effluent often releases pharmaceuticals like venlafaxine (a serotonin-norephinephrine reuptake inhibitor antidepressant) to freshwater ecosystems at levels causing adverse metabolic effects on fish. Changes to fish metabolism can be regulated by epigenetic mechanisms like microRNA (small RNA molecules that regulate mRNA translation), including regulating mitochondrial mRNAs. Nuclear-encoded microRNAs regulate mitochondrial gene expression in mammals, and have predicted effects in fish. We aimed to identify whether venlafaxine exposure changed mitochondrial respiration and resulted in differentially abundant mitochondrial microRNA (mitomiRs) in zebrafish brains. In vitro exposure of brain homogenate to below environmentally relevant concentrations of venlafaxine (<1 µg/L) caused a decrease in mitochondrial respiration, although this was not driven by changes to mitochondrial Complex I or II function. To identify whether these effects occur in vivo, zebrafish were exposed to 1 µg/L venlafaxine for 0, 1, 6, 12, 24, and 96 h. In vivo, venlafaxine exposure had no significant effects on brain mitochondrial respiration; however, select mitomiRs (dre-miR-301a-5p, dre-miR-301b-3p, and dre-miR-301c-3p) were also measured, because they were bioinformatically predicted to regulate mitochondrial cytochrome c oxidase subunit I (COI) abundance. These mitomiRs were differentially regulated based on venlafaxine exposure (with miR-301c-3p abundance differing during the day and miR-301b-3p being lower in exposed fish at night), and with respect to sex and time sampled. Overall, the results demonstrated that in vitro venlafaxine exposure to zebrafish brain caused a decrease in mitochondrial respiration, but these effects were not seen after acute in vivo exposure. Results may have differed because in vivo exposure allows for fish to mitigate effects through mechanisms that could include mitomiR regulation, and because fish were only acutely exposed. Environ Toxicol Chem 2024;43:1569-1582. © 2024 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Complete mitogenome and intra-family comparative mitogenomics showed distinct position of Pama Croaker Otolithoides pama

The Pama Croaker, Otolithoides pama, is an economically important fish species in Bangladesh. Intra-family similarities in morphology and typical barcode sequences of cox1 create ambiguities in its identification. Therefore, morphology and the complete mitochondrial genome of O. pama, and comparative mitogenomics within the family Sciaenidae have been studied. Extracted genomic DNA was subjected to Illumina-based short read sequencing for De-Novo mitogenome assembly. The complete mitogenome of O. pama (Accession: OQ784575.1) was 16,513 bp, with strong AC biasness and strand asymmetry. Relative synonymous codon usage (RSCU) among 13 protein-coding genes (PCGs) of O. pama was also analyzed. The studied mitogenomes including O. pama exhibited consistent sizes and gene orders, except for the genus Johnius which possessed notably longer mitogenomes with unique gene rearrangements. Different genetic distance metrics across 30 species of Sciaenidae family demonstrated 12S rRNA and the control region (CR) as the most conserved and variable regions, respectively, while most of the PCGs undergone a purifying selection. Different phylogenetic trees were congruent with one another, where O. pama was distinctly placed. This study would contribute to distinguishing closely related fish species of Sciaenidae family and can be instrumental in conserving the genetic diversity of O. pama.

Genomic analysis and phylogenetic characterization of Himalayan snow trout, Schizothorax esocinus based on mitochondrial protein-coding genes

BACKGROUND

Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes.

METHODS

Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis.

RESULTS

The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%.

CONCLUSION

The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species.

Global patterns of nuclear and mitochondrial genetic diversity in marine fishes

Genetic diversity is a fundamental component of biodiversity. Examination of global patterns of genetic diversity can help highlight mechanisms underlying species diversity, though a recurring challenge has been that patterns may vary by molecular marker. Here, we compiled 6862 observations of genetic diversity from 492 species of marine fish and tested among hypotheses for diversity gradients: the founder effect hypothesis, the kinetic energy hypothesis, and the productivity-diversity hypothesis. We fit generalized linear mixed effect models (GLMMs) and explored the extent to which various macroecological drivers (latitude, longitude, temperature (SST), and chlorophyll-a concentration) explained variation in genetic diversity. We found that mitochondrial genetic diversity followed geographic gradients similar to those of species diversity, being highest near the Equator, particularly in the Coral Triangle, while nuclear genetic diversity did not follow clear geographic patterns. Despite these differences, all genetic diversity metrics were correlated with chlorophyll-a concentration, while mitochondrial diversity was also positively associated with SST. Our results provide support for the kinetic energy hypothesis, which predicts that elevated mutation rates at higher temperatures increase mitochondrial but not necessarily nuclear diversity, and the productivity-diversity hypothesis, which posits that resource-rich regions support larger populations with greater genetic diversity. Overall, these findings reveal how environmental variables can influence mutation rates and genetic drift in the ocean, caution against using mitochondrial macrogenetic patterns as proxies for whole-genome diversity, and aid in defining global gradients of genetic diversity.

Complete mitochondrial genome of critically endangered catfish Hemibagrus punctatus (Jerdon, 1849) and comparative analysis for insights into the phylogeny of hemibagrids through mitogenomic approach

BACKGROUND

Hemibagrus punctatus (Jerdon, 1849) is a critically endangered bagrid catfish endemic to the Western Ghats of India, whose population is declining due to anthropogenic activities. The current study aims to compare the mitogenome of H. punctatus with that of other Bagrid catfishes and provide insights into their evolutionary relationships.

METHODS AND RESULTS

Samples were collected from Hemmige Karnataka, India. In the present study, the mitogenome of H. punctatus was successfully assembled, and its phylogenetic relationships with other Bagridae species were studied. The total genomic DNA of samples was extracted following the phenol-chloroform isoamyl alcohol method. Samples were sequenced, and the Illumina paired-end reads were assembled to a contig length of 16,517 bp. The mitochondrial genome was annotated using MitoFish and MitoAnnotator (Iwasaki et al., 2013). A robust phylogenetic analysis employing NJ (Maximum composite likelihood) and ASAP methods supports the classification of H. punctatus within the Bagridae family, which validates the taxonomic status of this species. In conclusion, this research enriches our understanding of H. punctatus mitogenome, shedding light on its evolutionary dynamics within the Bagridae family and contributing to the broader knowledge of mitochondrial genes in the context of evolutionary biology.

CONCLUSIONS

The study's findings contribute to a better understanding of the mitogenome of H. punctatus and provide insights into the evolutionary relationships within other Hemibagrids.

Molecular Identification of Acrossocheilus jishouensis (Teleostei: Cyprinidae) and Its Complete Mitochondrial Genome

Molecular identification, such as DNA barcoding, is a useful tool that is widely applied in distinguishing species. To identify the cyprinid Acrossocheilus jishouensis, which was previously known to be restricted to only its type locality, we conducted molecular identification of this species based on 23 samples in five localities. Molecular identification based on the mitochondrial COI gene sequence showed that the morphologically similar samples from the five populations were all A. jishouensis, as the mean genetic distances between populations were very small (0.1-1.6%); thus, the distribution of this species was substantially expanded. The whole mitochondrial genome of one sample was also assembled, which was 16,594 bp in length and consisted of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. All PCGs began with ATG except the COI gene, which started with GTG; seven PCGs used the complete stop codon TAA, while four terminated in T(AA) and two ended with TAG. The overall base composition reflected a higher proportion of A+T than G+C and a positive AT-skew and negative GC-skew pattern except for the opposite in ND6. Phylogenetic relationships inferred using BI and ML methods revealed that both Acrossocheilus and Onychostoma were nonmonophyletic, which indicated that the traditional diagnoses between these two genera need to be assessed further. The results of this study not only expanded the known distribution ranges of A. jishouensis, but also provided a valuable data resource for future molecular and evolutionary studies of Acrossocheilus and other cyprinids in Barbinae.

Complete mitogenome of Lepidopygopsis typus, an evolutionarily-distinct, endangered cyprinid fish from the Western Ghats Biodiversity Hotspot: Phylogenetic relationships and implications for conservation

The mitogenome of Lepidopygopsis typus, an evolutionarily distinct, endangered, cyprinid fish from the Western Ghats Biodiversity Hotspot, was characterized. Total length of the mitogenome was 16,729 bp, and comprised of 13 protein coding, 22 tRNA and two rRNA genes. Thirteen protein coding genes contributed a total nucleotide length of 11,408 bp, which coded for a total of 3794 amino acids and accounting for 68.2 % of the mitogenome. A maximum likelihood phylogenetic tree based on mitogenomes of members of cyprinid subfamilies Torinae, Schizothoracinae and Schizopygopsinae, recovered L. typus in a well-resolved clade of Torinae with members of the Western Ghats endemic genus Hypselobarbus as sister taxa. Selection analysis suggested that the branch of L. typus and Hypselobarbus spp. was under intensified selection with reference to other members of Torinae. There were 19 codons under diversifying selection in L. typus, which could be the result of positive selection for adapting to high-altitude, upstream tributaries of Periyar River, where the species has a restricted distribution. With respect to Hypselobarbus spp., L. typus was under relaxed selection with about 68 % of codons experiencing neutral evolution. Restricted distribution, low population size and relaxation in selection can likely trigger extinction in L. typus, and therefore urgent conservation and monitoring plans are required to secure the future of this evolutionary distinct, and globally endangered species.

Mitochondrial genome structure and composition in 70 fishes: a key resource for fisheries management in the South Atlantic

BACKGROUND

Phylogenetic gaps of public databases of reference sequences are a major obstacle for comparative genomics and management of marine resources, particularly in the Global South, where economically important fisheries and conservation flagship species often lack closely-related references. We applied target-enrichment to obtain complete mitochondrial genomes of marine ichthyofauna from the Brazilian coast selected based on economic significance, conservation status and lack of phylogenetically-close references. These included sardines (Dorosomatidae, Alosidae), mackerels (Scombridae) croakers (Sciaenidae), groupers (Epinephelidae) and snappers (Lutjanidae).

RESULTS

Custom baits were designed to enrich mitochondrial DNA across a broad phylogenetic range of fishes. Sequencing generated approximately 100k reads per sample, which were assembled in a total of 70 complete mitochondrial genomes and include fifty-two new additions to GenBank, including five species with no previous mitochondrial data. Departures from the typical gene content and order occurred in only three taxa and mostly involved tRNA gene duplications. Start-codons for all genes, except Cytochrome C Oxidase subunit I (COI), were consistently ATG, whilst a wide range of stop-codons deviated from the prevailing TAA. Phylogenetic analysis confirmed assembly accuracy and revealed signs of cryptic diversification within the Mullus genus. Lineage delimitation methods using Sardinella aurita and S. brasiliensis mitochondrial genomes support a single Operational Taxonomic Unit.

CONCLUSIONS

Target enrichment was highly efficient, providing complete novel mitochondrial genomes with little sequencing effort. These sequences are deposited in public databases to enable subsequent studies in population genetics and adaptation of Latin American fish species and serve as a vital resource for conservation and management programs that rely on molecular data for species and genus-level identification.

Characterization of the Complete Mitochondrial Genome of Schizothorax kozlovi (Cypriniformes, Cyprinidae, Schizothorax) and Insights into the Phylogenetic Relationships of Schizothorax

Schizothorax kozlovi is an endemic and vulnerable fish species found in the upper Yangtze River in China. Over the past few years, the population resources of S. kozlovi have been nearly completely depleted owing to multiple contributing threats. While the complete mitochondrial genomes serve as important molecular markers for phylogenetic and genetic studies, the mitochondrial genome of S. kozlovi has still received little attention. In this study, we analyzed the characterization of the mitochondrial genome of S. kozlovi and investigated the phylogenetic relationships of Schizothorax. The complete mitochondrial genome of S. kozlovi was 16,585 bp in length, which contained thirty-seven genes (thirteen protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), twenty-two transfer RNA genes (tRNAs)) and two non-coding regions for the origin of light strand (OL) and the control region (CR). There were nine overlapping regions and seventeen intergenic spacers regions in the mitochondrial genome. The genome also showed a bias towards A + T content (55.01%) and had a positive AT-skew (0.08) and a negative GC-skew (-0.20). All the PCGs employed the ATG or GTG as the start codon and TAA, TAG, or single T as the stop codon. Additionally, all of the tRNAs displayed a typical cloverleaf secondary structure, except trnS1 which lacked the D arm. The phylogenetic analysis, based on the maximum likelihood (ML) and Bayesian inference (BI) methods, revealed that the topologies of the phylogenetic tree divided the Schizothorax into four clades and did not support the classification of Schizothorax based on morphology. The phylogenetic status of S. kozlovi was closely related to that of S. chongi. The present study provides valuable genomic information for S. kozlovi and new insights in phylogenetic relationships of Schizothorax. These data could also offer fundamental references and guidelines for the management and conservation of S. kozlovi and other species of Schizothorax.

Comparative mitochondrial genome brings insights to slight variation in gene proportion and large intergenic spacer and phylogenetic relationship of mudskipper species

Fish mitochondrial genome have been largely studied worldwide for evolutionary and other genetic purposes and the structure and gene organization are commonly conservative. However, several studies have demonstrated that this scenario may present variations in some taxa, showing differentiation on the gene rearrangement. In this study, the complete mitogenome of terrestrial fish Boleophthalmus dussumieri was generated and compared with other species of the Exudercidae fishes. The newly complete mitogenome generated is circular and 16,685 bp of length, and it contained 13 protein-coding genes (PCGs), two ribosomal RNA (rRNAs), 22 transfer RNA genes (tRNAs), and one control region (CR), with high conservative structure, like other Mudskippers. Most of the PCG showed similar codon usage bias. The gene length was found to be different specially for the CR, 12S rRNA gene and ND5 gene in some taxon. All the Boleophthalmus species showed a gene duplication in the CR, except for B. dussumieri, and they presented a long intergenic spacer specially on the tRNA-Pro/ OH Tandem duplication/random loss (TDRL) and dimer-mitogenome and nonrandom loss (DMNL) are suitable to explain the mitogenome rearrangement observed in this study. The phylogenetic analysis well supported the monophyly of all mudskipper species and the analysis positioned the Periophthalmus clade as the most basal of the terrestrial fishes. This finding provides basis and brings insights for gene variation, gene rearrangements and replications showing evidence for variety of mitochondrial structure diversity within mudskippers.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae)

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

Adaptive evolution of mitochondrial genomes in Triplophysa cavefishes

Extreme conditions in caves pose survival challenges for cave dwellers, who gradually develop adaptive survival features. Cavefishes are one of the most successful animals among cave dwellers. Triplophysa cavefishes are an important group of cavefishes, and they show remarkable adaptability to the extreme environments of caves. However, there is a limited understanding of their adaptation mechanisms. In this study, eight complete mitochondrial genomes of Triplophysa cavefishes were newly obtained, and their genomic characteristics, including the base composition, base bias, and codon usage, were analyzed. Phylogenetic analysis was carried out based on 13 mitochondrial protein-coding genes from 44 Nemacheilidae species. This showed that Triplophysa cavefishes and non-cavefishes separate into two reciprocally monophyletic clades, suggesting a single origin of the cave phenotype. Positive selection analysis strongly suggested that the selection pressure in cavefishes is higher than that in non-cavefishes. Furthermore, the ND5 gene in cavefishes showed evidence of positive selection, which suggests that the gene may play an important role in the adaptation of cavefishes to the cave environment. Protein structure analysis of the ND5 subunit implied that the sites of positive selection in cavefishes might allow them to acquire lower ND5 protein stability, compared to that in non-cavefishes, which might help the accumulation of nonsynonymous (mildly deleterious) mutations. Together, our study revealed the genetic signatures of cave adaptation in Triplophysa cavefishes from the perspective of energy metabolism.

Four new complete mitochondrial genomes of Gobioninae fishes (Teleostei: Cyprinidae) and their phylogenetic implications

The subfamily Gobioninae is one of the most diverse fish groups within Cyprinidae. Their taxonomy and phylogenetic relationships are not completely resolved. In this study, the complete mitochondrial genomes (mitogenome) of four Gobioninae species (Microphysogobio elongatus, Microphysogobio chinssuensis, Gobio rivuloides and Rhinogobio nasutus) were sequenced and compared. The mitogenomes of four species ranges from 16603 bp to 16609 bp in length, consisting of 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. Most PCGs had significant codon usage bias. Except for the tRNASer (GCT), all the nucleotide substitutions of tRNA loops higher than the stems could fold into a stable secondary structure. The nucleotide compositions of Gobioninae mitogenome were biased toward A/T, and NAD4 was subjected to low purification selection and had a faster evolution rate among 13 PCGs. Bayesian inference and maximum likelihood phylogenetic analyses showed the consistent results. The four sequenced species clustered together with their congener species. However, more samples and mitogenome data are needed to untangle the phylogenetic relationships among genera Microphysogobio, Romanogobio, Hugobio, Biwia and Platysmacheilus.

Description of eight new mitochondrial genomes for the genus Neoarius and phylogenetic considerations for the family Ariidae (Siluriformes)

The genus Neoarius, known as marine catfish, is a group of the family Ariidae, composed of 10 species found in Oceania. None of the species in this genus have their mitochondrial genome described, which is highly valuable in phylogenetic and molecular evolution studies. For the present work, eight species from the Neoarius genus were selected: Neoarius utarus, Neoarius midgleyi, Neoarius graeffei, Neoarius leptaspis, Neoarius berenyi, Neoarius paucus, Neoarius pectoralis, and Neoarius aff. graeffei. DNA sequences of the eight species were obtained through the NCBI Sequence Read Archive (SRA) database, and the mitochondrial genomes were assembled using the NOVOplasty tool on the Galaxy platform, subsequently annotated with the MitoAnnotator tool. We then utilized the protein-coding genes from the mitogenomes to estimate the phylogenetic relationships within the group, including seven additional mitogenomes available in the NCBI. In all species, the mitochondrial genomes presented 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 D-loop.

Mitochondrial replication's role in vertebrate mtDNA strand asymmetry

Mitogenomes are defined as compact and structurally stable over aeons. This perception results from a vertebrate-centric vision, where few types of mtDNA rearrangements are described. Here, we bring a new light to the involvement of mitochondrial replication in the strand asymmetry of the vertebrate mtDNA. Using several species of deep-sea hatchetfish (Sternoptychidae) displaying distinct mtDNA structural arrangements, we unravel the inversion of the coding direction of protein-coding genes (PCGs). This unexpected change is coupled with a strand asymmetry nucleotide composition reversal and is shown to be directly related to the strand location of the Control Region (CR). An analysis of the fourfold redundant sites of the PCGs (greater than 6000 vertebrates), revealed the rarity of this phenomenon, found in nine fish species (five deep-sea hatchetfish). Curiously, in Antarctic notothenioid fishes (Trematominae), where a single PCG inversion (the only other record in fish) is coupled with the inversion of the CR, the standard asymmetry is disrupted for the remaining PCGs but not yet reversed, suggesting a transitory state. Our results hint that a relaxation of the classic vertebrate mitochondrial structural stasis promotes disruption of the natural balance of asymmetry of the mtDNA. These findings support the long-lasting hypothesis that replication is the main molecular mechanism promoting the strand-specific compositional bias of this unique and indispensable molecule.

The complete mitochondrial genome of the Tanichthys flavianalis (Cypriniformes: Tanichthyidae) and its phylogeny

In the present study, we have sequenced, annotated and characterized the complete mitogenome of Tanichthys flavianalis Li, Liao & Shen, 2022, which is an endemic minnow to Hainan Province, China. The entire mitogenome size is 16,545 bp, containing 37 genes coding for 13 mitochondrial proteins, 22 tRNAs, 2 rRNAs, and one putative control region. It has the typical fish mitochondrial gene arrangement, though with 9 gene overlaps and 17 intergenic spacers scattering throughout the mitogenome, and a high A + T bias (60.8%) comparing to other Cyprinids. ATG is used as the start codon in most mitochondrial protein coding genes (PCGs) except for cox1 which uses GTG, whereas the canonical TAA/TAG are prevailing stop codons except for cox2, cox3, nad4 and cytb. Phylogenetic analyses with other reported species from Tanichyidae were conducted using IQ-TREE and MrBayes, based on the concatenated 13 PCGs. The results clustered T. flavinalis with the newly described T. albiventris Li Bohlen & Liao, 2022, and proved the reciprocally monophyletic relationship between Chinese and Vietnamese Tanichthys.

Genetic Diversity in Natural Populations of the Near-Threatened Species Lignobrycon myersi (Characiformes, Triportheidae): Implications for Species Conservation

The river basins of Brazil contain a highly diverse ichthyofauna of remarkable endemism, including several threatened species. Accordingly, Lignobrycon myersi is a fish species distributed only in a few rivers from the state of Bahia, northeastern Brazil. Since this species is classified as Near Threatened and is poorly studied, efforts to understand the genetic structure of populations and putative cryptic forms should help define efficient strategies of management and conservation. Herein, the molecular identification and the population genetic diversity of specimens of L. myersi across their range (Almada, Contas, and Cachoeira river basins) were assessed using mitochondrial markers (16S rDNA and D-Loop, respectively). The inferences based on phylogenetics, genetic distance, and species delimitation methods invariably identified all samples as L. myersi. In addition, sequencing of D-loop fragments revealed significant haplotype diversity and a considerable level of population genetic structure. Despite their geographic isolation, these data suggested that populations from Almada and Contas rivers represent a single evolutionary lineage that could be managed as a whole. In contrast, the population from Cachoeira River was highly differentiated from the others and should be managed separately as a unique and endemic unit, particularly focused on the conservation of native habitats.

The New Mitochondrial Genome of Hemiculterella wui (Cypriniformes, Xenocyprididae): Sequence, Structure, and Phylogenetic Analyses

Hemiculterella wui is an endemic small freshwater fish, distributed in the Pearl River system and Qiantang River, China. In this study, we identified and annotated the complete mitochondrial genome sequence of H. wui. The mitochondrial genome was 16,619 bp in length and contained 13 protein coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. The nucleotide composition of the mitochondrial genome was 29.9% A, 25.3% T, 27.4% C, and 17.5% G, respectively. Most PCGs used the ATG start codon, except COI and ATPase 8 started with the GTG start codon. Five PCGs used the TAA termination codon and ATPase 8 ended with the TAG stop codon, and the remaining seven genes used two incomplete stop codons (T and TA). Most of the tRNA genes showed classical cloverleaf secondary structures, except that tRNASer(AGY) lacked the dihydrouracil loop. The average Ka/Ks value of the ATPase 8 gene was the highest, while the average Ka/Ks value of the COI gene was the lowest. Phylogenetic analyses showed that H. wui has a very close relationship with Pseudohemiculter dispar and H. sauvagei. This study will provide a valuable basis for further studies of taxonomy and phylogenetic analyses in H. wui and Xenocyprididae.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

Trends of Eurasian Perch (Perca fluviatilis) mtDNA ATP6 Region Genetic Diversity within the Hydro-Systems of the Eastern Part of the Baltic Sea in the Anthropocene

The intraspecific genetic diversity of freshwater fish inhabiting hydro-systems of the macrogeographic area spreading from the Black to Baltic Seas requires comprehensive investigation from fundamental and practical perspectives. The current study focused on the involvement of the mtDNA ATP6 region in the adaptability and microevolution of Perca fluviatilis within phylogeographic and anthropogenic contexts. We sequenced a 627 bp fragment encompassing the ATP6 region and used it for genetic analysis of 193 perch caught in Latvia, Lithuania, Belarus, and Ukraine, representing natural and anthropogenically impacted populations. We evaluated patterns of intraspecific genetic diversity in the ATP6 region and phylogeographic trends within the studied area compared with previously established D-loop trends. Evaluation of ATP6 coding sequence variability revealed that among 13 newly detected haplotypes, only two were caused by non-synonymous substitutions of amino acids of the protein. PCoA revealed three genetic groups (I-III) based on the ATP6 region that encompassed four previously described genetic groups established based on the mtDNA D-loop. The two mtDNA regions (D-loop and ATP6) have microevolved at least partially independently. Prolonged anthropogenic impacts may generate new point mutations at the ATP6 locus, but this phenomenon could be mainly concealed by natural selection and reparation processes.

Mitophylogeny of Pangasiid Catfishes and its Taxonomic Implications for Pangasiidae and the Suborder Siluroidei

Pangasiidae (catfish order: Siluriformes) comprises 30 valid catfish species in four genera: Pangasius, Pangasianodon, Helicophagus, and Pseudolais. Their systematics are frequently revised due to the addition of newly described species. Although Pangasiidae is known to be a monophyletic family, the generic and phylogenetic relationships among the taxa are poorly resolved. This study characterized three newly obtained complete mitogenomes of Mekong River catfishes from Vietnam (Pangasius mekongensis, Pangasius krempfi, and Pangasianodon hypophthalmus), as well as the inter-and intrafamilial relationships of the Pangasiidae and catfish families in Siluroidei. The genomic features of their mitogenomes were similar to those of previously reported pangasiids, including all regulatory elements, extended terminal associated sequences (ETAS), and conserved sequence blocks (CSBs) (CSB-1, CSB-2, CSB-3, and CSBs, A to F) in the control region. A comprehensive phylogeny constructed from datasets of multiple 13 PCG sequences from 117 complete mitogenomes of 32 recognized siluriform families established Pangasiidae as monophyletic and a sister group of Austroglanididae. The [Pangasiidae + Austroglanididae] + (Ictaluridae + Cranoglanididae) + Ariidae] clade is a sister to the "Big Africa" major clade of Siluriformes. Furthermore, both phylogenies constructed from the single barcodes (83 partial cox1 and 80 partial cytB, respectively) clearly indicate genus relationships within Pangasiidae. Pangasianodon was monophyletic and a sister to the (Pangasius + Helicophagus + Pseudolais) group. Within the genus Pangasius, P. mekongensis was placed as a sister taxon to P. pangasius. Pangasius sanitwongsei was found to be related to and grouped with Pangasianodon, but in single-gene phylogenies, it was assigned to the Pangasius + Helicophagus + Pseudolais group. The datasets in this study are useful for studying pangasiid systematics, taxonomy and evolution.

Comparative Mitogenome Analyses Uncover Mitogenome Features and Phylogenetic Implications of the Reef Fish Family Holocentridae (Holocentriformes)

To understand the molecular mechanisms and adaptive strategies of holocentrid fish, we sequenced the mitogenome of eight species within the family Holocentridae and compared them with six other holocentrid species. The mitogenomes were found to be 16,507-16,639 bp in length and to encode 37 typical mitochondrial genes, including 13 PCGs, two ribosomal RNAs, and 22 transfer RNA genes. Structurally, the gene arrangement, base composition, codon usage, tRNA size, and putative secondary structures were comparable between species. Of the 13 PCGs, nad6 was the most specific gene that exhibited negative AT-skews and positive GC-skews. Most of the genes begin with the standard codon ATG, except cox1, which begins with the codon GTG. By examining their phylogeny, Sargocentron and Neoniphon were verified to be closely related and to belong to the same subfamily Holocentrinae, while Myripristis and Ostichthys belong to the other subfamily Myripristinae. The subfamilies were clearly distinguished by high-confidence-supported clades, which provide evidence to explain the differences in morphology and feeding habits between the two subfamilies. Selection pressure analysis indicated that all PCGs were subject to purifying selection. Overall, our study provides valuable insight into the habiting behavior, evolution, and ecological roles of these important marine fish.