Comparative mitogenome analyses uncover mitogenome features and phylogenetic implications of the subfamily Cobitinae

Whole-genome relaxed selection and molecular constraints in Triplophysa under adapted Qinghai-Tibetan Plateau

High-altitude environments are inhospitable, but Triplophysa, the largest taxon among the three major fish groups in the Qinghai-Tibetan Plateau (QTP), is an exception. However, the evolutionary profiling of the common ancestor and its contribution to the adaptation of existing QTP native species is unclear.We researched the comparative genomics of Triplophysa species and found that the genome-wide genes of Triplophysa and its ancestry have the characteristics of rapid evolution.Moreover, the rapid evolution of the ancestral genes was caused by relaxed selection. Natural selection analysis showed that more ancestral relaxed selection genes were under strongly purifying selection and showed higher expression in QTP endemic Triplophysa species.The change in natural selection might be associated with the adaptation to QTP. It should be noted that SPT5 homolog, DSIF elongation factor subunit (supt5h) experienced relaxed selection in common ancestral populations of Triplophysa but under purifying selection in extant species, which might be related to hypoxia adaptation of QTP. In summary, the extant species in different environments were used to infer the evolutionary profile of the common ancestor and to identify candidate genes based on changes in natural selection. Our work might provide new clues for understanding adaptation to extreme environments.

Comparative morphology and mitogenomics of freshwater mussels Koreosolenaia, Parvasolenaia, and Sinosolenaia (Bivalvia: Unionidae: Gonideinae)

BACKGROUND

Amidst the escalating loss of global biodiversity, freshwater mussels (family Unionidae) have become one of the most imperiled animal groups. Acquiring more biological and phylogenetic information on understudied taxa constitutes a pivotal aspect of conservation biology. Consequently, a comprehensive examination was conducted on Koreosolenaia, Parvasolenaia, and Sinosolenaia from China encompassing morphology, anatomy, distribution, and molecular systematics to provide theoretical support for future species endangerment assessments and biodiversity conservation.

RESULTS

The shell characteristics of Koreosolenaia, Parvasolenaia, and Sinosolenaia were clearly distinct, and the soft-body morphology could also be easily distinguished from each other. The papillae of the incurrent aperture of Sinosolenaia iridinea, Sinosolenaia recognita, and Sinosolenaia oleivora, which were previously described as difficult, exhibited significant variations that could be utilized for species diagnosis. Furthermore, both incurrent and excurrent apertures of the Sinosolenaia species had small cysts on their dorsal surfaces which may be unique to this particular group. Comparative analysis of six mitochondrial genomes (Parvasolenaia rivularis, Koreosolenaia sitgyensis, Sinosolenaia iridinea, Sinosolenaia recognita, Sinosolenaia carinata, and Sinosolenaia oleivora) revealed a completely consistent gene arrangement pattern. Additionally, there was a high consistency in nucleotide base content and skewness, amino acid usage, and relative synonymous codon usage among the six complete mitochondrial genomes. Mitochondrial phylogenomics of these genomes with additional taxa within Gonideinae robustly supported the generic relationships as follows: (Inversidens + ((Microcondylaea + Sinosolenaia) + (Parvasolenaia + (Koreosolenaia + (Ptychorhynchus + (Postolata + Cosmopseudodon)))))).

CONCLUSIONS

The present study provided significant data on the shell morphology and soft-body anatomy of Koreosolenaia, Parvasolenaia, and Sinosolenaia, thereby clarifying the diagnostic characteristics for these challenging taxa. Additionally, we established a robust phylogenetic framework at both the generic and species levels based on mitochondrial genomics.

Complete mitochondrial genome sequence of Nannostomus eques and comparative analysis with Nannostomus beckfordi

The brown pencilfish, Nannostomus eques is a lebiasinid harvested for ornamental purposes; however, its complete mitochondrial genome sequence is still unknown. To enrich the molecular genetic information pertaining to Nannostomus, we present here the first report of the complete mitochondrial genome sequence of Nannostomus eques and compare it with Nannostomus beckfordi. The total lengths of the N. eques and N. beckfordi mitochondrial genomes were 16,673 bp and 16,742 bp, respectively, and there was a double-stranded ring with a heavy chain and a light chain in the whole structures of both. We used PhyloSuite v1.2.1 to construct the maximum likelihood and Bayesian Analysis trees based on tRNAs, rRNAs, and protein-coding genes (PCGs) data and compared them with other Nannostomus species by referring to other studies. Our study found that N. beckfordi has a closer genetic relationship with N. eques than with Lebiasina astrigata, which belongs to the same family, and we discovered some similarities and even rules in Nannostomus species. Our study provides a molecular basis for the conservation and sustainable use of Nannostomus species.

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.

Characterization and phylogenetic implications of newly sequenced mitochondrial genomes of cobitid fish Acantopsis Rungthipae (Boyd, Nithirojpakdee & Page, 2017)

BACKGROUND

Acantopsis rungthipae has significant ornamental and ecological value. This study aimed at structurally characterizing the A. rungthipae mitochondrial genome and elucidate its phylogenetic position in Cobitidae.

METHODS AND RESULTS

High-throughput sequencing technology was used to obtain the complete sequence of the mitochondrial genome of A. rungthipae and reconstruct a Cobitidae phylogenetic tree based on the sequence of 13 protein-coding genes. The entire mitochondrial genome of A. rungthipae was 16,600 bp, containing 22 tRNA genes, 13 protein-coding genes, 2 rRNAs, and 2 non-coding regions (D-loop and OL). The base composition showed a significant AT preference, with the highest A + T content (67.1%) in the D-loop region. Among the protein-coding genes, 12 had ATG as a typical starting codon, while only COXI had GTG as a special starting codon. Twenty-one of the tRNA genes exhibited clover structure, and only tRNA-Ser (GCT) could not fold into a clover structure because of the absence of DHU arms. The phylogenetic tree was reconstructed using the Bayesian and maximum likelihood methods and showed that A. rungthipae and Acantopsis choirorhynchos converged into one branch, and their phylogenetic relationships were relatively close.

CONCLUSIONS

Our findings supplement basic data on the A. rungthipae mitochondrial genome and deepen the understanding of the evolutionary relationships of the genus Acantopsis. Clarifying the evolutionary relationships between different species in Acantopsis lays a solid foundation for subsequent research on fish adaptation and selection pressure.

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.

Comparative mitogenomes provide new insights into phylogeny and taxonomy of the subfamily Xenocyprinae (Cypriniformes: Cyprinidae)

Xenocyprinae is a cyprinid subfamily that not only has a discrete geographic distribution but also has a long history dating to the Early Miocene. However, it is controversial whether systematic classification and some species validity of Xenocyprinae exist, as well as its phylogenetic relationships and evolutionary history. In the present study, we first reviewed the description and taxonomic history of Xenocyprinae, and then the complete mitochondrial genome of Distoechodon compressus, an endemic and locally distributed species belonging to Xenocyprinae, was sequenced and annotated. Finally, all the mitogenomes of Xenocyprinae were compared to reconstruct the phylogenetic relationship and estimate the divergence time. The results showed that the mitogenomes are similar in organization and structure with 16618–16630 bp length from 12 mitogenomes of eight species. Phylogenetic analysis confirmed the monology of Xenocyprinae and illustrated three clades within the Xenocyprinae to consist of ambiguous generic classification. Plagiognathops is a valid genus located at the base of the phylogenetic tree. The genus Xenocypris was originally monophyletic, but X. fangi was excluded. Divergence time estimation revealed that the earliest divergence within Xenocyprinae occurred approximately 12.1 Mya when Plagiognathops separated from the primitive Xenocypris. The main two clades (Xenocypris and (Distoechodon + Pseudobrama + X. fangi)) diverged 10.0 Mya. The major divergence of Xenocyprinae species possibly occurred in the Middle to Late Miocene and Late Pliocene, suggesting that speciation and diversifications could be attributed to the Asian monsoon climate. This study clarifies some controversial issues of systematics and provides essential information on the taxonomy and phylogeny of the subfamily Xenocyprinae.

Comparative Mitogenomics of Two Sympatric Catfishes of Exostoma (Siluriformes: Sisoridae) from the Lower Yarlung Tsangpo River and Its Application for Phylogenetic Consideration

The genus Exostoma is a group of stenotopic and rheophilic glyptosternine catfishes distributed in South and Southeast Asia. So far, comprehensive studies on mitogenomics referring to this genus are very scarce. In this study, we first sequenced and annotated the complete mitochondrial genomes of Exostoma tibetanum and Exostoma tenuicaudatum—two sympatric congeners from the lower Yarlung Tsangpo River, Tibet, China. The mitogenomes of both species contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one light-strand origin of replication, and one control region, with lengths of 16,528 bp and 16,533 bp, respectively. The mitogenome architecture, nucleotide composition, and codon usage of protein-coding genes were almost identical between the two Exostoma species, although some estimated parameters varied. Phylogenetic analysis strongly supported the monophyly of Exostoma in the subfamily Glyptosternae, and Exostoma tibetanum had the closest relationship to Exostoma tenuicaudatum. The divergence time estimation demonstrated that these two species diverged approximately 1.51 Ma during the early Pleistocene, which was speculated to be triggered by the river system changes caused by the uplift of the southeastern Tibetan Plateau. Selection pressure analyses indicated that all protein-coding genes of Exostoma species underwent a strong purifying selection, while minority positive sites from NADH dehydrogenase complex genes were detected. These findings are expected to promote our understanding of the molecular phylogeny of the genus Exostoma and provide valuable mitogenomic resources for the subfamily Glyptosternae

Complete mitochondrial genomes of two catfishes (Siluriformes, Bagridae) and their phylogenetic implications

The mitochondrial genome (mitogenome) has been widely used as a molecular marker to investigate phylogenetic analysis and evolutionary history in fish. However, the study of mitogenomes is still scarce in the family Bagridae. In this study, the mitogenomes of Tachysurusbrachyrhabdion and T.gracilis were sequenced, annotated, and analyzed. The mitogenomes were found to be 16,532 bp and 16,533 bp, respectively, and each contained 37 typical mitochondrial genes, which are 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and a control region. All PCGs begin with the codon ATG, except for the cytochrome c oxidase subunit 1 (COI) gene, while seven PCGs end with an incomplete termination codon. All tRNA genes can fold into their typical cloverleaf secondary structures, except for tRNA^Ser(AGY), which lacks the dihydrouracil arm. The Ka/Ks ratios for all PCGs are far lower than one. Phylogenetic analyses based on Bayesian inference (BI) and maximum likelihood (ML) showed that the two clades in Bagridae excluded Ritarita. The monophyly of Tachysurus supports previous research and the traditional classification that Leiocassis, Pseudobagrus, Pelteobagrus, and Tachysurus belong to one genus (Tachysurus). These findings provide a phylogenetic basis for future phylogenetic and taxonomic studies of Bagridae.

Comparative Mitogenome Analyses of Subgenera and Species Groups in Epeorus (Ephemeroptera: Heptageniidae)

Simple Summary

As one of the most species-rich genera of Ephemeroptera, Epeorus Eaton, 1881, was found to be widely distributed in Holarctic and Oriental regions, and nine subgenera have been reported. Previous phylogenetic studies of Epeorus were mainly focused on morphological characters or several gene fragments. Here, 15 mitogenomes of Epeorus are sequenced and the comparative mitogenome analysis of six subgenera is performed. The gene rearrangement of trnI-trnM-trnQ-NCR-ND2 was first found in the genus. In addition, differences in genetic composition and codon usage between the species with this kind of rearrangement and other Epeorus species were observed. Phylogenetic analyses show that three subgenera, Proepeorus, Belovius and Iron, are not monophyletic groups, and our results imply that gill structures are not always appropriate for the classification of subgenera in Epeorus.

Abstract

Epeorus Eaton, 1881 is a diverse mayfly genus in Heptageniidae comprising more than 100 species which are further divided into nine subgenera and several species groups. However, the classification and the phylogenetic relationships among them are still uncertain. Here, 15 complete mitochondrial genomes of Epeorus were sequenced and compared together with six available ones of same genus in the NCBI database. Based on morphological classification, the 21 mitogenomes were classified into six subgenera (Proepeorus, Epeorus s.str., Belovius, Iron, Caucasiron and Siniron) and four species groups (G1, G2, montanus and longimanus). Among all analyzed mitogenomes, the gene rearrangement of trnI-trnM-trnQ-NCR-ND2 was first found occurring in three species of group G1, whereas the gene block trnI-trnM-trnQ-trnM-ND2 was observed in all other mitogenomes of Epeorus. Furthermore, the genetic composition and codon usage of species in group G1 were also significantly different from all other Epeorus species, except group longimanus. The intergenic spacer between trnA and trnR, which has the stem-loop secondary structure, occurred in all 21 mitogenomes, and the sequences of stems and loops were conserved within species groups. Furthermore, the phylogenetic analyses strongly support the monophyly of all species groups, although three of six recognized subgenera Proepeorus, Belovius, and Iron, were shown as the non-monophyletic groups.

Complete Mitochondrial Genomes of Five Racerunners (Lacertidae: Eremias) and Comparison with Other Lacertids: Insights into the Structure and Evolution of the Control Region

Comparative studies on mitochondrial genomes (mitogenomes) as well as the structure and evolution of the mitochondrial control region are few in the Lacertidae family. Here, the complete mitogenomes of five individuals of Eremias scripta (2 individuals), Eremias nikolskii, Eremias szczerbaki, and Eremias yarkandensis were determined using next-generation sequencing and were compared with other lacertids available in GenBank. The circular mitogenomes comprised the standard set of 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and a long non-coding control region (CR). The extent of purifying selection was less pronounced for the COIII and ND2 genes in comparison with the rest of the PCGs. The codons encoding Leucine (CUN), Threonine, and Isolecucine were the three most frequently present. The secondary structure of rRNA of Lacertidae (herein, E. scripta KZL15 as an example) comprised four domains and 28 helices for 12S rRNA, with six domains and 50 helices for 16S rRNA. Five types and twenty-one subtypes of CR in Lacertidae were described by following the criteria of the presence and position of tandem repeats (TR), termination-associated sequence 1 (TAS1), termination-associated sequence 2 (TAS2), conserved sequence block 1 (CBS1), conserved sequence block 2 (CSB2), and conserved sequence block 3 (CSB3). The compositions of conserved structural elements in four genera, Acanthodactylus, Darevskia, Eremias, and Takydromus, were further explored in detail. The base composition of TAS2 - TATACATTAT in Lacertidae was updated. In addition, the motif "TAGCGGCTTTTTTG" of tandem repeats in Eremias and the motif "GCGGCTT" in Takydromus were presented. Nucleotide lengths between CSB2 and CSB3 remained 35 bp in Eremias and Darevskia. The phylogenetic analyses of Lacertidae recovered the higher-level relationships among the three subfamilies and corroborated a hard polytomy in the Lacertinae phylogeny. The phylogenetic position of E. nikolskii challenged the monophyly of the subgenus Pareremias within Eremias. Some mismatches between the types of CR and their phylogeny demonstrated the complicated evolutionary signals of CR such as convergent evolution. These findings will promote research on the structure and evolution of the CR and highlight the need for more mitogenomes in Lacertidae.

Complete mitogenomes reveal limited genetic variability in the garden dormouse Eliomys quercinus of the Iberian Peninsula

The garden dormouse Eliomys quercinus is a poorly known Western Palearctic species experiencing a global decline. Even though the availability of genetic information is key to assess the driversunderlying demographic changes in wild populations and plan adequate management, data on E. quercinus are still scant. In this study, we reconstructed the complete mitogenomes of four E. quercinus individuals from southern Spain using in–solution enriched libraries, and found evidence of limited genetic variability. We then compared their cytochrome b sequences to those of conspecifics from other countries and supported the divergent but genetically depauperate position of this evolutionarily significant unit (ESU). The information produced will assist future conservation studies on this little–studied rodent.