Mitogenome Sequence Research Articles

The complete mitochondrial genome of Hemiclepsis yangtzenensis (Clitellata: Glossiphoniidae)

We present the complete mitochondrial genome sequence of a recently described new leech species named Hemiclepsis yangtzenensis Yang & Bolotov 2021 collected in central China. The mitochondrial genome is 14,984 bp in length and consists of 13 protein-coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes, all of which are encoded on a single strand. It exhibited a strong A + T bias of 72.87%. There is a large non-coding region (614 bp) located between the tRNA-Arg and tRNA-His genes, wherein we identified 40 short dispersed repeats, 13–22 bp long, 8 of which were direct, 20 inverted, and 12 palindromic. Phylogenetic analysis of 20 Hirudinea mitogenome sequences resolved monophyletic Glossiphoniidae, and H. yangtzenensis formed a sister lineage with Glossiphonia concolor.

Morphometrical characterization of the Atlantic mudskipper species (Periophthalmus barbarus) (Linnaeus, 1766) (Perciformes; Gobiiae) from Abonema in Port Harcourt, Rivers State, Nigeria

Mudskippers are borderline organisms whose biology, systematics and evolution has attracted research attention for so long. There are controversies in their systematics especially because the morphological classifies are plastic and can be dependent on life stage, living conditions and the expertise of the scientist. This paper sort to characterize the Nigerian mudskippers using firstly morphometry, then DNA barcoding and the complete mitogenome sequences. Ten parameters were measured on One hundred specimens and the result was subjected to exploratory data analysis. The Nigerian mudskipper is in size and shape not different from other species found around the world. The standard deviation of all the measurements range from low (0.24-055) for six of the parameters, medium for two (0.56-0.86) and high for the remaining two (0.87-1.17). Low and medium standard deviation is an indication of very low variability in the parameters with the data clustering around the mean closely. Scatter plotting indicated five of the specimen that can be classified as outliers. Principal component analysis indicated two namely PC1 and PC2 as being responsible for 68% of the variance observed. Principal component plots 1 and 2 indicated head width, head depth and snout length loaded positively while only body depth loaded positively for the PC2. The Nigerian species were also amphidromous. The essence of this initial analysis was to see if body measurements will indicate any significant grouping as we plan to select samples for molecular analysis.

The Carniolan Honeybee from Slovenia-A Complete and Annotated Mitochondrial Genome with Comparisons to Closely Related Apis mellifera Subspecies.

Simple SummaryThe western honeybee, Apis mellifera, is a globally distributed bee species with many recognised subspecies, one of which is Apis mellifera carnica, the Carniolan honeybee. Apis m. carnica is native to southern Central Europe and parts of the Balkans, with the locus classicus in Slovenia. It is also widely popular with beekeepers in parts of Central and Northern Europe and other parts of the world, including the USA, Canada, and even New Zealand. In Slovenia, A. m. carnica is protected, with measures to conserve the subspecies’ autochthonous domestic population in place. Such efforts depend heavily upon genomic and phylogenetic information. In this study, we sequenced and annotated the mitochondrial genome of a specimen from Slovenia and compared the obtained data with a previously published sample of the A. m. carnica from Austria and the closely related Italian honeybee A. m. ligustica. We found several features unique to the new mitochondrial genome. We also phylogenetically analyzed the relationship between our sequence and the selected available A. mellifera mitochondrial sequences. The acquired position of the sequenced A. m. carnica from Slovenia on the phylogenetic tree brings new evidence for close relationships among C and O lineages and reflects their recent historical matrilinear ancestry.The complete mitochondrial genome of the Carniolan honeybee (Apis mellifera carnica) from Slovenia, a homeland of this subspecies, was acquired in two contigs from WGS data and annotated. The newly obtained mitochondrial genome is a circular closed loop of 16,447 bp. It comprises 37 genes (13 protein coding genes, 22 tRNA genes, and 2 rRNA genes) and an AT-rich control region. The order of the tRNA genes resembles the order characteristic of A. mellifera. The mitogenomic sequence of A. m. carnica from Slovenia contains 44 uniquely coded sites in comparison to the closely related subspecies A. m. ligustica and to A. m. carnica from Austria. Furthermore, 24 differences were recognised in comparison between A. m. carnica and A. m. ligustica subspecies. Among them, there are three SNPs that affect translation in the nd2, nd4, and cox2 genes, respectively. The phylogenetic placement of A. m. carnica from Slovenia within C lineage deviates from the expected position and changes the perspective on relationship between C and O lineages. The results of this study represent a valuable addition to the information available in the phylogenomic studies of A. mellifera—a pollinator species of worldwide importance. Such genomic information is essential for this local subspecies’ conservation and preservation as well as its breeding and selection.

Five new mitogenomes sequences of Calidridine sandpipers (Aves: Charadriiformes) and comparative mitogenomics of genus Calidris.

BackgroundThe genus Calidris (Charadriiformes, Scolopacidae) includes shorebirds known as dunlin, knots, and sanderlings. The relationships between species nested within Calidris, including Eurynorynchus, Limicola and Aphriza, are not well-resolved.MethodsSamples were collected from Xiaoyangkou, Rudong County, Jiangsu Province, China. Mitogenomes were sequenced using the Illumina Novaseq 6000 platform for PE 2 × 150 bp sequencing, and then checked for PCR products. Protein-coding genes were determined using an Open Reading Frame Finder. tRNAscan-SE, MITOS, and ARWEN were used to confirm tRNA and rRNA annotations. Bioinformatic analyses were conducted using DnaSP 5.1 and MEGA X. Phylogenic trees were constructed using maximum likelihood and Bayesian analyses.ResultsWe sequenced and annotated the mitogenome of five species and obtained four complete mitogenomes and one nearly complete mitogenome. Circular mitogenomes displayed moderate size variation, with a mean length of 16,747 bp, ranging from 16,642 to 16,791 bp. The mitogenome encoded a control region and a typical set of 37 genes containing two rRNA genes, 13 protein-coding genes, and 22 tRNA genes. There were four start codons, four stop codons, and one incomplete stop codon (T–). The nucleotide composition was consistently AT-biased. The average uncorrected pairwise distances revealed heterogeneity in the evolutionary rate for each gene; the COIII had a slow evolutionary rate, whereas the ATP8 gene had a fast rate. dN/dS analysis indicated that the protein-coding genes were under purifying selection. The genetic distances between species showed that the greatest genetic distance was between Eurynorhynchus pygmeus and Limicola falcinellus (22.5%), and the shortest was between E. pygmeus and Calidris ruficollis (12.8%). Phylogenetic trees revealed that Calidris is not a monophyletic genus, as species from the genera Eurynorynchus and Limicola were nested within Calidris. The molecular data obtained in this study are valuable for research on the taxonomy, population genetics, and evolution of birds in the genus Calidris.

Comparative Analyses of Complete Peronosporaceae (Oomycota) Mitogenome Sequences-Insights into Structural Evolution and Phylogeny.

Members of the Peronosporaceae (Oomycota, Chromista), which currently consists of 25 genera and approximately 1,000 recognized species, are responsible for disease on a wide range of plant hosts. Molecular phylogenetic analyses over the last two decades have improved our understanding of evolutionary relationships within Peronosporaceae. To date, 16 numbered and three named clades have been recognized; it is clear from these studies that the current taxonomy does not reflect evolutionary relationships. Whole organelle genome sequences are an increasingly important source of phylogenetic information, and in this study, we present comparative and phylogenetic analyses of mitogenome sequences from 15 of the 19 currently recognized clades of Peronosporaceae, including 44 newly assembled sequences. Our analyses suggest strong conservation of mitogenome size and gene content across Peronosporaceae but, as previous studies have suggested, limited conservation of synteny. Specifically, we identified 28 distinct syntenies amongst the 71 examined isolates. Moreover, 19 of the isolates contained inverted or direct repeats, suggesting repeated sequences may be more common than previously thought. In terms of phylogenetic relationships, our analyses of 34 concatenated mitochondrial gene sequences resulted in a topology that was broadly consistent with previous studies. However, unlike previous studies concatenated mitochondrial sequences provided strong support for higher-level relationships within the family.

Characterization of the mitochondrial genomes of Bradysia hygida, Phytosciara flavipes and Trichosia splendens (Diptera: Sciaridae) and novel insights on the control region of sciarid mitogenomes.

Sciarids, also called "fungus gnats" are small, almost entirely dark-coloured insects. Sciarid larvae feed on different substrates and can infest agricultural crops and mushroom nurseries, causing economic losses. Of the 2174 Diptera mitogenome sequences currently available in GenBank, only eight are from the Sciaridae family, none of which are complete circular molecules. Here we describe the mitogenome sequences of three sciarid species: Phytosciara flavipes, Trichosia splendens and Bradysia hygida and provide novel insights on the control region of sciarid mitogenomes. The assembled mitogenomes range from 16,062 bp in P. flavipes to 17,095 bp in B. hygida. All 13 protein coding genes, 22 tRNAs and 2 rRNAs characteristic of insect mitogenomes were identified, but the sequence of the control region could not be determined. Experimental results suggest that the B. hygida control region is about 21 kb long resulting in a 37 kb long mitogenome which constitutes the largest insect mitochondrial genome described so far. Phylogenetic analysis using all Bibionomorpha mitogenome sequences available in GenBank strongly supports the Sciaridae monophyly and led to the identification of species and subfamily specific gene rearrangements. Our study extends the knowledge of this large and diverse insect family that includes agricultural pest species.

Three complete mitochondrial genomes of shortfin mako sharks, Isurus oxyrinchus, from the Atlantic and Pacific Oceans

We present complete mitogenome sequences of three shortfin mako sharks (Isurus oxyrinchus) sampled from the western Pacific, and eastern and western Atlantic oceans. Mitogenome sequence lengths ranged between 16,699 bp and 16,702 bp, and all three mitogenomes contained one non-coding control region, two rRNA genes, 22 tRNA genes, and 13 protein-coding genes. Comparative assessment of five mitogenomes from globally distributed shortfin makos (the current three and two previously published mitogenomes) yielded 98.4% identity, with the protein-coding genes ATP8, ATP6, and ND5 as the most variable regions (sequence identities of 96.4%, 96.5%, and 97.6%, respectively). These mitogenome sequences contribute resources for assessing the genetic population dynamics of this endangered oceanic apex predator.

The complete mitochondrial genome of Sinopodisma hengshanica (Orthoptera: Acrididae) and its phylogenetic implication

The complete mitochondrial genome (mitogenome) of Sinopodisma hengshanica (Orthoptera: Acrididae) was firstly determined and analyzed in the present study. Assembled mitogenome sequence of S. hengshanica is 15,623 bp in size, containing 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), as well as one A + T-rich region. Its gene component and arrangement are identical with other Acrididae species. The overall nucleotide composition is as follows: A (42.89%), G (10.34%), T (33.07%), and C (13.70%). All PCGs are initiated by typical ATN codons, and terminated with harboring complete stop codons TAA and TAG. Furthermore, phylogenetic trees were reconstructed based on 13 PCGs to validate the taxonomic status of S. hengshanica, exhibiting a close relationship with Sinopodisma rostellocerca.

Mitogenome Analyses of Water Buffalo: Closeness of the Genetic Architecture of River Buffalo and Wild Buffalo (Bubalus arnee) Excludes Swamp Buffalo

Water buffalo species belongs to genus Bubalus and involves the wild and the domestic types. Domestic water buffaloes are reclassified into river and swamp types. Here, we have reconstructed the phylogenetic relationship between members of water buffalo species based on mitogenome sequences. The phylogenetic tree included lowland anoa, mountain anoain addition to wild, river and swamp water buffalo types.Results indicated the paraphyly of water buffalo species since lowland anoa and mountain anoa were clustered with river and wild water buffalo in a separate clade, leaving swamp buffalo at the outside. Results also indicated close genetic architecture between the mitogenomes of river and wild water buffalo types. These findings may indicate that river buffalo and swamp buffalo belong to two distinct species and that river buffalo and notswamp buffalo is a descendant of wild water buffalo.Mitogenome analysis provided more information about the taxonomic status of water buffalo types and verified their relationship with the members of genus Bubalus.

Comparative genomic analysis of vertebrate mitochondrial reveals a differential of rearrangements rate between taxonomic class

Vertebrate mitochondrial genomes have been extensively studied for genetic and evolutionary purposes, these are normally believed to be extremely conserved, however, different cases of gene rearrangements have been reported. To verify the level of rearrangement and the mitogenome evolution, we performed a comparative genomic analysis of the 2831 vertebrate mitochondrial genomes representing 12 classes available in the NCBI database. Using a combination of bioinformatics methods, we determined there is a high number of errors in the annotation of mitochondrial genes, especially in tRNAs. We determined there is a large variation in the proportion of rearrangements per gene and per taxonomic class, with higher values observed in Actinopteri, Amphibia and Reptilia. We highlight that these are results for currently available vertebrate sequences, so an increase in sequence representativeness in some groups may alter the rearrangement rates, so in a few years it would be interesting to see if these rates are maintained or altered with the new mitogenome sequences. In addition, within each vertebrate class, different patterns in rearrangement proportion with distinct hotspots in the mitochondrial genome were found. We also determined that there are eleven convergence events in gene rearrangement, nine of which are new reports to the scientific community.

Characterization of the first mitogenomes of the smallest fish in the world, Paedocypris progenetica, from peat swamp of Peninsular Malaysia, Selangor, and Perak.

The two complete mitochondrial genomes (mitogenomes) of Paedocypris progenetica, the smallest fish in the world which belonged to the Cyprinidae family, were sequenced and assembled. The circular DNA molecules of mitogenomes P1-P. progenetica and S3-P. progenetica were 16,827 and 16,616 bp in length, respectively, and encoded 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region. The gene arrangements of P. progenetica were identical to those of other Paedocypris species. BLAST and phylogenetic analyses revealed variations in the mitogenome sequences of two Paedocypris species from Perak and Selangor. The circular DNA molecule of P. progenetica yield a standard vertebrate gene arrangement and an overall nucleotide composition of A 33.0%, T 27.2%, C 23.5%, and G 15.5%. The overall AT content of this species was consistent with that of other species in other genera. The negative GC-skew and positive AT-skew of the control region in P. progenetica indicated rich genetic variability and AT nucleotide bias, respectively. The results of this study provide genomic variation information and enhance the understanding of the mitogenome of P. progenetica. They could later deliver highly valuable new insight into data for phylogenetic analysis and population genetics.

Variation in the Mitochondrial Genome of the Chagas Disease Vector Triatoma infestans (Hemiptera: Reduviidae).

Chagas' disease is transmitted mainly by members of the subfamily Triatominae (Hemiptera: Reduviidae). Among them, Triatoma infestans (Klug) is the main vector of the disease in Southern Cone of Latin America. In order to contribute to knowledge of the genetic variation between triatomine vectors, in the present study, we analyzed the intraspecific and interspecific variations of the seven mitogenomes available from Triatominae. In addition, in order to examine their evolutionary relationships with others species of Reduviidae and to estimate the divergence time of the main lineages, we constructed phylogenetic trees including mitogenome sequences of 30 species from Reduviidae. Comparative analysis between mitochondrial DNA sequences from two specimens of T. infestans revealed a total of 54 variable sites. Triatoma infestans, Triatoma dimidiata (Latreille), Triatoma rubrofasciata (De Geer), Triatoma migrans (Breddin), Rhodnius pictipes (Stål), and Panstrongylus rufotuberculatus (Champion) present similar mitogenome organization and the length differences observed among these species are primarily caused by variations in control region (CR) and intergenic spacers (IGS). The relative synonymous codon usage values (RSCU) were similar in the six species of Triatominae, and in agreement with the observed in other insects, a biased use of A and C nucleotides in the majority strand was detected. The monophyly of five subfamilies was strongly supported (Phymatinae, Peiratinae, Triatominae, Stenopodainae, and Harpactorinae), while the sampled species of Reduviinae were grouped with one specie from the Salyavatinae subfamily. The oldest subfamily is Phymatinae at 100.3 Mya (99.6-102.2 Mya) and the youngest is Triatominae and Stenopodainae at 52.6 Mya (42.5-63.7 Mya). The estimated diversification time for the Triatominae subfamily agrees with the Andean uplift geological event. An analysis with more mitogenomes from more Triatominae species would be necessary to provide sufficient evidence to support this finding.

Comparative mitogenomics of freshwater snails of the genus Bulinus, obligatory vectors of Schistosoma haematobium, causative agent of human urogenital schistosomiasis

Among the snail genera most responsible for vectoring human-infecting schistosomes, Bulinus, Biomphalaria, and Oncomelania, the former is in many respects the most important. Bulinid snails host the most common human blood fluke, Schistosoma haematobium, responsible for approximately two-thirds of the estimated 237 million cases of schistosomiasis. They also support transmission of schistosomes to millions of domestic and wild animals. Nonetheless, our basic knowledge of the 37 Bulinus species remains incomplete, especially with respect to genome information, even including mitogenome sequences. We determined complete mitogenome sequences for Bulinus truncatus, B. nasutus, and B. ugandae, and three representatives of B. globosus from eastern, central, and western Kenya. A difference of the location of tRNA-Asp was found between mitogenomes from the three species of the Bulinus africanus group and B. truncatus. Phylogenetic analysis using partial cox1 sequences suggests that B. globosus is a complex comprised of multiple species. We also highlight the status of B. ugandae as a distinct species with unusual interactions with the S. haematobium group parasites deserving of additional investigation. We provide sequence data for potential development of genetic markers for specific or intraspecific Bulinus studies, help elucidate the relationships among Bulinus species, and suggest ways in which mitogenomes may help understand the complex interactions between Schistosoma and Bulinus snails and their relatives.

Highly diversified mitochondrial genomes provide new evidence for interordinal relationships in the Arachnida.

Arachnida is an exceptionally diverse class in the Arthropoda, consisting of 20 orders and playing crucial roles in the terrestrial ecosystems. However, their interordinal relationships have been debated for over a century. Rearranged or highly rearranged mitochondrial genomes (mitogenomes) were consistently found in this class, but their various extent in different lineages and efficiency for resolving arachnid phylogenies are unclear. Here, we reconstructed phylogenetic trees using mitogenome sequences of 290 arachnid species to decipher interordinal relationships as well as diversification through time. Our results recovered monophyly of ten orders (i.e. Amblypygi, Araneae, Ixodida, Mesostigmata, Opiliones, Pseudoscorpiones, Ricinulei, Sarcoptiformes, Scorpiones and Solifugae), while rejecting monophyly of the Trombidiformes due to the unstable position of the Eriophyoidea. The monophyly of Acari (subclass) was rejected, possibly due to the long-branch attraction of the Pseudoscorpiones. The monophyly of Arachnida was further rejected because the Xiphosura nested within arachnid orders with unstable positions. Mitogenomes that are highly rearranged in mites but less rearranged or conserved in the remaining lineages point to their exceptional diversification in mite orders; however, shared derived mitochondrial (mt) gene clusters were found within superfamilies rather than interorders, confusing phylogenetic signals in arachnid interordinal relationships. Molecular dating results show that arachnid orders have ancient origins, ranging from the Ordovician to the Carboniferous, yet have significantly diversified since the Cretaceous in orders Araneae, Mesostigmata, Sarcoptiformes, and Trombidiformes. By summarizing previously resolved key positions of some orders, we propose a plausible arachnid tree of life. Our results underline a more precise framework for interordinal phylogeny in the Arachnida and provide new insights into their ancient evolution.

Maternal Lineages of Gepids from Transylvania.

According to the written historical sources, the Gepids were a Germanic tribe that settled in the Carpathian Basin during the Migration Period. They were allies of the Huns, and an independent Gepid Kingdom arose after the collapse of the Hun Empire. In this period, the Carpathian Basin was characterized by so-called row-grave cemeteries. Due to the scarcity of historical and archaeological data, we have a poor knowledge of the origin and composition of these barbarian populations, and this is still a subject of debate. To better understand the genetic legacy of migration period societies, we obtained 46 full mitogenome sequences from three Gepid cemeteries located in Transylvania, Romania. The studied samples represent the Classical Gepidic period and illustrate the genetic make-up of this group from the late 5th and early 6th centuries AD, which is characterized by cultural markers associated with the Gepid culture in Transylvania. The genetic structure of the Gepid people is explored for the first time, providing new insights into the genetic makeup of this archaic group. The retrieved genetic data showed mainly the presence of Northwestern European mitochondrial ancient lineages in the Gepid group and all population genetic analyses reiterated the same genetic structure, showing that early ancient mitogenomes from Europe were the major contributors to the Gepid maternal genetic pool.

Characterization of the complete mitochondrial genome of Laomedia astacina De Haan, 1841 (Decapoda: Laomediidae) and its phylogenetic implications

In this study, we determined the complete mitochondrial genome (mitogenome) of Laomedia astacina De Haan, 1841 using next-generation sequencing technology. The total length of the mitogenome sequence of L. astacina is 14,795 base pairs, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, and two ribosomal RNA genes. The overall composition of the mitogenome is estimated to be 35.3% A, 38.0% T, 13.8% C, and 12.9% G, indicating that the L. astacina mitogenome is rich in A + T (73.3%). The phylogenetic relationships of 13 decapod species were constructed based on the 13 PCGs by the maximum-likelihood approach using IQtree software.

The complete mitochondrial genome of Philus pallescens Bates, 1866 (Coleoptera: Vesperidae) and its phylogenetic implications

We sequenced and assembled the complete mitochondrial genome of Philus pallescens from Madou, Tainan County, Taiwan. The complete mitogenome of P. pallescens is 15,750 bp long, and contains 13 protein-coding, 22 tRNA and two rDNA genes. Nucleotide compositions of the mitogenome of P. pallescens are A: 38.08%, T: 32.25%, C: 18.67%, and G: 11.00%. The AT and GC skewness of the mitogenome sequence were 0.0828 and −0.25845 respectively, showing the genome composition skewed toward adenine and cytosine. The phylogenetic position of Chrysomelidae is sister to all the other families in the superfamily Chrysomeloidea. The results indicate that Chrysomeloidea Cerambycidae is not a monophyletic group. Philus is phylogenetically close to Spiniphilus. Vesperidae is monophyletic and sister to Disteniidae. Mitogenomic data from this study will provide useful information for further studies on the population genetics, speciation, and pest management of P. pallescens.

The complete mitochondrial genome of a satin moth: Leucoma chrysoscela (Lepidoptera: Erebidae, Lymantriinae)

We describe the mitogenome sequence of Leucoma chrysoscela (Collenette, 1934) collected in the Longtan National Forest Park, which is located in the southeast of China. The assembled mitogenome is 15,508 bp in length and consists 13 protein coding genes, 22 transfer-RNA genes, 2 ribosomal-RNA genes, and one A + T rich region. The most common start codon for 13 PCGs is ATT and the most common termination codon is TAA. The overall G + C content was only 20.45% in the heavy strand. The result of phylogenetic analysis shows that the relationship of L. chrysoscela is close to the species in the same subfamily Lymantriinae.

The phylogeny of the Anderson's White-bellied Rat (Niviventer andersoni) based on complete mitochondrial genomes.

The phylogenetic structure of the genus Niviventer has been studied based on several individual mitochondrial and nuclear genes, but the results seem to be inconsistent. In order to clarify the phylogeny of Niviventer, we sequenced the complete mitochondrial genome of white‐bellied rat (Niviventer andersoni of the family Muridae) by next‐generation sequencing. The 16,291 bp mitochondrial genome consists of 22 transfer RNA genes, 13 protein‐coding genes (PCGs), two ribosomal RNA genes, and one noncoding control region (D‐Loop). Phylogenetic analyses of the nucleotide sequences of all 13 PCGs, PCGs minus ND6, and the entire mitogenome sequence except for the D‐loop revealed well‐resolved topologies supporting that N. andersoni was clustered with N. excelsior forming a sister division with N. confucianus, which statistically rejected the hypothesis based on the tree of cytochrome b (cytb) gene that N. confucianus is sister to N. fulvescens. Our research provides the first annotated complete mitochondrial genome of N. andersoni, extending the understanding about taxonomy and mitogenomic evolution of the genus Niviventer.

Mitogenomics of the suborder Cottoidei (Teleostei: Perciformes): Improved assemblies, mitogenome features, phylogeny, and ecological implications.

We determined the mitogenome of Cyclopterus lumpus using a hybrid sequencing approach, and another four closely related species in the Liparidae based on available next-generation sequence data. We found that the mitogenome of C. lumpus was 17,266bp in length, where the length and organisation were comparable to those reported for cottoids. However, we found a GC-homopolymer region in the intergenic space between tRNALeu2 and ND1 in liparids and cyclopterids. Phylogenetic reconstruction confirmed the monophyly of infraorders and firmly supported a sister-group relationship between Cyclopteridae and Liparidae. Purifying selection was the predominant force in the evolution of cottoid mitogenomes. There was significant evidence of relaxed selective pressures along the lineage of deep-sea fish, while selection was intensified in the freshwater lineage. Overall, our analysis provides a necessary expansion in the availability of mitogenomic sequences and sheds light on mitogenomic adaptation in Cottoidei fish inhabiting different aquatic environments.