Evaluating the roles of energetic functional constraints on teleost mitochondrial-encoded protein evolution

Effect of locomotor preference on the evolution of mitochondrial genes in Bovidae

Locomotor preferences and habitat types may drive animal evolution. In this study, we speculated that locomotor preference and habitat type may have diverse influences on Bovidae mitochondrial genes. We used selection pressure and statistical analysis to explore the evolution of mitochondrial DNA (mtDNA) protein-coding genes (PCGs) from diverse locomotor preferences and habitat types. Our study demonstrates that locomotor preference (energy demand) drives the evolution of Bovidae in mtDNA PCGs. The habitat types had no significant effect on the rate of evolution in Bovidae mitochondrial genes. Our study provides deep insight into the adaptation of Bovidae.

Mitochondrial Protein-Coding Gene Expression in the Lizard Sphenomorphus incognitus (Squamata:Scincidae) Responding to Different Temperature Stresses

In the context of global warming, the frequency of severe weather occurrences, such as unexpected cold spells and heat waves, will grow, as well as the intensity of these natural disasters. Lizards, as a large group of reptiles, are ectothermic. Their body temperatures are predominantly regulated by their environment and temperature variations directly impact their behavior and physiological activities. Frequent cold periods and heat waves can affect their biochemistry and physiology, and often their ability to maintain their body temperature. Mitochondria, as the center of energy metabolism, are crucial for maintaining body temperature, regulating metabolic rate, and preventing cellular oxidative damage. Here, we used RT-qPCR technology to investigate the expression patterns and their differences for the 13 mitochondrial PCGs in Sphenomorphus incognitus (Squamata:Scincidae), also known as the brown forest skink, under extreme temperature stress at 4 °C, 8 °C, 34 °C, and 38 °C for 24 h, compared to the control group at 25 °C. In southern China, for lizards, 4 °C is close to lethal, and 8 °C induces hibernation, while 34/38 °C is considered hot and environmentally realistic. Results showed that at a low temperature of 4 °C for 24 h, transcript levels of ATP8, ND1, ND4, COI, and ND4L significantly decreased, to values of 0.52 ± 0.08, 0.65 ± 0.04, 0.68 ± 0.10, 0.28 ± 0.02, and 0.35 ± 0.02, respectively, compared with controls. By contrast, transcript levels of COIII exhibited a significant increase, with a mean value of 1.86 ± 0.21. However, exposure to 8 °C for 24 h did not lead to an increase in transcript levels. Indeed, transcript levels of ATP6, ATP8, ND1, ND3, and ND4 were significantly downregulated, to 0.48 ± 0.11, 0.68 ± 0.07, 0.41 ± 0.08, 0.54 ± 0.10, and 0.52 ± 0.07, respectively, as compared with controls. Exposure to a hot environment of 34 °C for 24 h led to an increase in transcript levels of COI, COII, COIII, ND3, ND5, CYTB, and ATP6, with values that were 3.3 ± 0.24, 2.0 ± 0.2, 2.70 ± 1.06, 1.57 ± 0,08, 1.47 ± 0.13, 1.39 ± 0.56, and 1.86 ± 0.12, respectively, over controls. By contrast, ND4L exhibited a significant decrease (to 0.31 ± 0.01) compared with controls. When exposed to 38 °C, the transcript levels of the 13 PCGs significantly increased, ranging from a 2.04 ± 0.23 increase in ND1 to a 6.30 ± 0.96 rise in ND6. Under two different levels of cold and heat stress, the expression patterns of mitochondrial genes in S. incognitus vary, possibly associated with different strategies employed by this species in response to low and high temperatures, allowing for rapid compensatory adjustments in mitochondrial electron transport chain proteins in response to temperature changes. Furthermore, this underscores once again the significant role of mitochondrial function in determining thermal plasticity in reptiles.

Complete Mitochondrial Genomes of Nedyopus patrioticus: New Insights into the Color Polymorphism of Millipedes

There has been debate about whether individuals with different color phenotypes should have different taxonomic status. In order to determine whether the different color phenotypes of Nedyopus patrioticus require separate taxonomic status or are simply synonyms, here, the complete mitochondrial genomes (mitogenomes) of two different colored N. patrioticus, i.e., red N. patrioticus and white N. patrioticus, are presented. The two mitogenomes were 15,781 bp and 15,798 bp in length, respectively. Each mitogenome contained 13 PCGs, 19 tRNAs, 2 rRNAs, and 1 CR, with a lack of trnI, trnL2, and trnV compared to other Polydesmida species. All genes were located on a single strand in two mitogenomes. Mitochondrial DNA analyses revealed that red N. patrioticus and white N. patrioticus did not show clear evolutionary differences. Furthermore, no significant divergence was discovered by means of base composition analysis. As a result, we suggest that white N. patrioticus might be regarded as a synonym for red N. patrioticus. The current findings confirmed the existence of color polymorphism in N. patrioticus, which provides exciting possibilities for future research. It is necessary to apply a combination of molecular and morphological methods in the taxonomy of millipedes.

Full-length transcriptome reveals the circularly polarized light response-related molecular genetic characteristics of Oratosquilla oratoria

The mantis shrimp is the only animal that can recognize circularly polarized light (CPL), but its molecular genetic characteristics are unclear. Multi-tissue level full-length (FL) transcriptome sequencing of Oratosquilla oratoria, a representative widely distributed mantis shrimp, was performed in the present study. We used comparative transcriptomics to explore the critical genes of O. oratoria selected by CPL and the GNβ gene associated with CPL signal transduction was hypothesized to be positively selected. Furthermore, the FL transcriptomes of O. oratoria compound eyes under five light conditions were sequenced and used to detect alternative splicing (AS). The ASs associated with CPL recognition mainly occurred in the LWS, ARR and TRPC regions. The number of FL transcripts with AS events and annotation information also provided evidence that O. oratoria could recognize LCPL. Additionally, 51 sequences belonging to the LWS, UV and Peropsin gene families were identified based on conserved 7tm domains. The LWS, UV and Peropsin opsins have similar 3D structures with seven domains across the cell membrane and conserved KSLRTPSN, DRY, and QAKK motifs. In conclusion, these results are undoubtedly valuable for perfecting the vision theory of O. oratoria and other mantis shrimp.

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.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Mitogenomic evolutionary rates in bilateria are influenced by parasitic lifestyle and locomotory capacity

The evidence that parasitic animals exhibit elevated mitogenomic evolutionary rates is inconsistent and limited to Arthropoda. Similarly, the evidence that mitogenomic evolution is faster in species with low locomotory capacity is limited to a handful of animal lineages. We hypothesised that these two variables are associated and that locomotory capacity is a major underlying factor driving the elevated rates in parasites. Here, we study the evolutionary rates of mitogenomes of 10,906 bilaterian species classified according to their locomotory capacity and parasitic/free-living life history. In Bilateria, evolutionary rates were by far the highest in endoparasites, much lower in ectoparasites with reduced locomotory capacity and free-living lineages with low locomotory capacity, followed by parasitoids, ectoparasites with high locomotory capacity, and finally micropredatory and free-living lineages. The life history categorisation (parasitism) explained ≈45%, locomotory capacity categorisation explained ≈39%, and together they explained ≈56% of the total variability in evolutionary rates of mitochondrial protein-coding genes in Bilateria. Our findings suggest that these two variables play major roles in calibrating the mitogenomic molecular clock in bilaterian animals.

Natural selection drives the evolution of mitogenomes in Acrossocheilus

The mitochondrial genome plays a crucial role in the balance of energy and heat production in organisms and, thus, may be under natural selection due to its potential role in adaptive divergence and speciation. Here, we studied natural selection on the mitogenome of Acrossocheilus (Cypriniformes: Cyprinidae), a genus of fish that inhabits a broad latitudinal distribution ranging from the tropics and subtropics through temperate regions. Specifically, we used 25 published mitogenome sequences of Acrossocheilus species to investigate phylogenetic relationships in this genus and detected signals of positive selection on 13 protein-coding, mitochondrial genes. We found that relaxed purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, we also found evidence of diversifying selection on some codons, indicating episodes of positive selection. Additionally, we analyzed the mitogenomic data within an environmental modeling framework and found that the Ka/Ks ratio of ATP6 may correlated with a mean diurnal temperature range (p = 0.0449), while the Ka/Ks ratio of COX2 may correlated with precipitation during the driest month (p = 0.00761). These results suggest that the mitogenomes of Acrossocheilus species may be involved in evolutionary adaptations to different habitats. Based on this, we believe that our study provides a new insight into the role of the mitochondrial genome of Acrossocheilus species in adaptation to different environments. During our study, we also discovered several cases of paraphyly and polyphyly among accessions of species and their putative synonyms. Thus, our study suggests that a careful reassessment of the taxonomy of Acrossocheilus is using high-quality molecular data merited.

The complete mitochondrial genome of the roosterfish Nematistius pectoralis Gill 1862: purifying selection in protein coding genes, organization of the control region, and insights into family-level phylogenomic relationships in the recently erected order Carangiformes

The roosterfish Nematistius pectoralis is considered as one of the most magnificent sportfishes worldwide. This study developed the first genomic resource for this trophy-fish that is heavily targeted by the fly-fishing industry. The complete mitochondrial genome of N. pectoralis was assembled using short read sequences and analyzed in detail. The mitochondrial genome of N. pectoralis is 16,537 bp in length and comprises 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (12S and 16S), and 22 transfer RNA genes. A long intergenic space 770 bp in length was assumed to be the D-loop or Control Region (CR). Most of the PCGs and tRNA genes are encoded in the L-strand. All PCGs are under purifying selection and atp8 and nad6 experienced the least selective pressure. All tRNAs exhibit a cloverleaf secondary structure except tRNA-Serine 1 that lacked the D-arm loop. The D-loop of N. pectoralis exhibits three domains commonly described in other fishes; extended terminal associated sequences (ETAS), central, and conserved sequence block (CSB). A ML phylogenetic reconstruction of the newly recognized order Carangiformes based on all 13 mitochondrial PCGs did not support the monophyly of this clade but recognized several families as monophyletic, including Bothidae, Carangidae, Istiophoridae, Latidae, Paralichthyidae, Polynemidae, and Rhombosoleidae. Nematistius pectoralis was sister to a clade composed of Toxotes chatareus (fam. Toxotidae) + Lactarius lactarius (fam. Lactariidae). This genomic resource developed for N. pectoralis will aid in improving our understanding of the population genomics of and strengthen conservation and management strategies in this remarkable trophy-fish.

The Complete Mitochondrial Genomes of Three Sphenomorphinae Species (Squamata: Scincidae) and the Selective Pressure Analysis on Mitochondrial Genomes of Limbless Isopachys gyldenstolpei

In order to adapt to diverse habitats, organisms often evolve corresponding adaptive mechanisms to cope with their survival needs. The species-rich family of Scincidae contains both limbed and limbless species, which differ fundamentally in their locomotor demands, such as relying on the movement of limbs or only body swing to move. Locomotion requires energy, and different types of locomotion have their own energy requirements. Mitochondria are the energy factories of living things, which provide a lot of energy for various physiological activities of organisms. Therefore, mitochondrial genomes could be tools to explore whether the limb loss of skinks are selected by adaptive evolution. Isopachys gyldenstolpei is a typical limbless skink. Here, we report the complete mitochondrial genomes of I. gyldenstolpei, Sphenomorphus indicus, and Tropidophorus hainanus. The latter two species were included as limbed comparator species to the limbless I. gyldenstolpei. The results showed that the full lengths of the mitochondrial genomes of I. gyldenstolpei, S. indicus, and T. hainanus were 17,210, 16,944, and 17,001 bp, respectively. Three mitochondrial genomes have typical circular double-stranded structures similar to other reptiles, including 13 protein-coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and the control region. Three mitochondrial genomes obtained in this study were combined with fifteen mitochondrially complete genomes of Scincidae in the NCBI database; the phylogenetic relationship between limbless I. gyldenstolpei and limbed skinks (S. indicus and T. hainanus) is discussed. Through BI and ML trees, Sphenomorphinae and Mabuyinae were monophyletic, while the paraphyly of Scincinae was also recovered. The limbless skink I. gyldenstolpei is closer to the species of Tropidophorus, which has formed a sister group with (T. hainanus + T. hangman). In the mitochondrial genome adaptations between limbless I. gyldenstolpei and limbed skinks, one positively selected site was found in the branch-site model analysis, which was located in ND2 (at position 28, BEB value = 0.907). Through analyzing the protein structure and function of the selected site, we found it was distributed in mitochondrial protein complex I. Positive selection of some mitochondrial genes in limbless skinks may be related to the requirement of energy to fit in their locomotion. Further research is still needed to confirm this conclusion though.

Molecular Phylogeny and Adaptive Mitochondrial DNA Evolution of Salmonids (Pisces: Salmonidae)

Salmonids are composed of anadromous and freshwater fishes, which is an important model for studying adaptive evolution. Herein, 49 salmonid complete mitochondrial genomes and those of two outgroups were used to infer a robust phylogeny for the family Salmonidae. The BI and RAxML phylogenetic trees based on 13 concatenated mitochondrial protein-coding genes showed well-supported nodes, and topologies were highly congruent. The concatenated 13 mitochondrial protein-coding genes, ND2, ND3, and ND5 genes were shown to have significantly larger dN/dS ratios in anadromous species than in freshwater species of Salmonidae, but the CYTB gene had significantly smaller dN/dS in anadromous species. The FEL analysis identified positively selected sites and negatively selected sites in each mitochondrial protein-coding gene separately. The RELAX program revealed that the ATP8 and CYTB genes supported intensified selection of the anadromous lineages. Our results demonstrated the phylogeny of Salmonidae and explored the mitochondrial DNA evolution pattern between anadromous and freshwater salmonids.

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

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

Insights into the Divergence of Chinese Ips Bark Beetles during Evolutionary Adaptation

Simple Summary

Bark beetle species of the genus Ips are among the major pests of Chinese conifer forests. Based on mitochondrial genome and SNP, we investigated the phylogenetic relationships and evolutionary trends of 19 populations of six Ips species that had serious outbreaks in recent years. Our results demonstrated the relationships between Ips evolution and host plants, pheromones, and altitudinal differences, and provided new insights into the mechanism of adaptive evolution of Ips bark beetles.

Abstract

Many bark beetles of the genus Ips are economically important insect pests that cause severe damage to conifer forests worldwide. In this study, sequencing the mitochondrial genome and restriction site-associated DNA of Ips bark beetles helps us understand their phylogenetic relationships, biogeographic history, and evolution of ecological traits (e.g., pheromones and host plants). Our results show that the same topology in phylogenetic trees constructed in different ways (ML/MP/BI) and with different data (mtDNA/SNP) helps us to clarify the phylogenetic relationships between Chinese Ips bark beetle populations and Euramerican species and their higher order clades; Ips bark beetles are polyphyletic. The structure of the mitochondrial genome of Ips bark beetles is similar and conserved to some extent, especially in the sibling species Ips typographus and Ips nitidus. Genetic differences among Ips species are mainly related to their geographic distribution and different hosts. The evolutionary pattern of aggregation pheromones of Ips species reflects their adaptations to the environment and differences among hosts in their evolutionary process. The evolution of Ips species is closely related to the uplift of the Qinghai-Tibet Plateau and host switching. Our study addresses the evolutionary trend and phylogenetic relationships of Ips bark beetles in China, and also provides a new perspective on the evolution of bark beetles and their relationships with host plants and pheromones.

Genetic Drift and Purifying Selection Shaped Mitochondrial Genome Variation in the High Royal Jelly-Producing Honeybee Strain (Apis mellifera ligustica)

Mitochondrial genomes (mitogenomes) are involved in cellular energy metabolism and have been shown to undergo adaptive evolution in organisms with increased energy-consuming activities. The genetically selected high royal jelly-producing bees (RJBs, Apis mellifera ligustica) in China can produce 10 times more royal jelly, a highly nutritional and functional food, relative to unselected Italian bees (ITBs). To test for potential adaptive evolution of RJB mitochondrial genes, we sequenced mitogenomes from 100 RJBs and 30 ITBs. Haplotype network and phylogenetic analysis indicate that RJBs and ITBs are not reciprocally monophyletic but mainly divided into the RJB- and ITB-dominant sublineages. The RJB-dominant sublineage proportion is 6-fold higher in RJBs (84/100) than in ITBs (4/30), which is mainly attributable to genetic drift rather than positive selection. The RJB-dominant sublineage exhibits a low genetic diversity due to purifying selection. Moreover, mitogenome abundance is not significantly different between RJBs and ITBs, thereby rejecting the association between mitogenome copy number and royal jelly-producing performance. Our findings demonstrate low genetic diversity levels of RJB mitogenomes and reveal genetic drift and purifying selection as potential forces driving RJB mitogenome evolution.

Complete mitochondrial genomes reveal robust phylogenetic signals and evidence of positive selection in horseshoe bats

BACKGROUND

In genus Rhinolophus, species in the Rhinolophus philippinensis and R. macrotis groups are unique because the horseshoe bats in these group have relatively low echolocation frequencies and flight speeds compared with other horseshoe bats with similar body size. The different characteristics among bat species suggest particular evolutionary processes may have occurred in this genus. To study the adaptive evidence in the mitochondrial genomes (mitogenomes) of rhinolophids, especially the mitogenomes of the species with low echolocation frequencies, we sequenced eight mitogenomes and used them for comparative studies of molecular phylogeny and adaptive evolution.

RESULTS

Phylogenetic analysis using whole mitogenome sequences produced robust results and provided phylogenetic signals that were better than those obtained using single genes. The results supported the recent establishment of the separate macrotis group. The signals of adaptive evolution discovered in the Rhinolophus species were tested for some of the codons in two genes (ND2 and ND6) that encode NADH dehydrogenases in oxidative phosphorylation system complex I. These genes have a background of widespread purifying selection. Signals of relaxed purifying selection and positive selection were found in ND2 and ND6, respectively, based on codon models and physicochemical profiles of amino acid replacements. However, no pronounced overlap was found for non-synonymous sites in the mitogenomes of all the species with low echolocation frequencies. A signal of positive selection for ND5 was found in the branch-site model when R. philippinensis was set as the foreground branch.

CONCLUSIONS

The mitogenomes provided robust phylogenetic signals that were much more informative than the signals obtained using single mitochondrial genes. Two mitochondrial genes that encoding proteins in the oxidative phosphorylation system showed some evidence of adaptive evolution in genus Rhinolophus and the positive selection signals were tested for ND5 in R. philippinensis. These results indicate that mitochondrial protein-coding genes were targets of adaptive evolution during the evolution of Rhinolophus species, which might have contributed to a diverse range of acoustic adaptations in this genus.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Sequence and structure comparison of ATP synthase F0 subunits 6 and 8 in notothenioid fish

Mitochondrial changes such as tight coupling of the mitochondria have facilitated sustained oxygen and respiratory activity in haemoglobin-less icefish of the Channichthyidae family. We aimed to characterise features in the sequence and structure of the proteins directly involved in proton transport, which have potential physiological implications. ATP synthase subunit a (ATP6) and subunit 8 (ATP8) are proteins that function as part of the F0 component (proton pump) of the F0F1complex. Both proteins are encoded by the mitochondrial genome and involved in oxidative phosphorylation. To explore mitochondrial sequence variation for ATP6 and ATP8 we analysed sequences from C. gunnari and C. rastrospinosus and compared them with their closely related red-blooded species and eight other vertebrate species. Our comparison of the amino acid sequence of these proteins reveals important differences that could underlie aspects of the unique physiology of the icefish. In this study we find that changes in the sequence of subunit a of the icefish C. gunnari at position 35 where there is a hydrophobic alanine which is not seen in the other notothenioids we analysed. An amino acid change of this type is significant since it may have a structural impact. The biology of the haemoglobin-less icefish is necessarily unique and any insights about these animals will help to generate a better overall understanding of important physiological pathways.

Geographic Variation of Phyllodiaptomus tunguidus Mitogenomes: Genetic Differentiation and Phylogeny

Phyllodiaptomus tunguidus (Copepoda: Calanoida) is largely endemic to and widespread in freshwater in southern China, where it inhabits a complex landscape from lowland to highland across an elevation gradient of 2000m. A deep genetic differentiation can be expected between its most distant geographic populations. Here, we sequenced nine mitogenomes from diverse populations. All mitogenomes contained 37 genes, including 13 protein-coding genes (PCG), two rRNA genes, 22 tRNA genes and one control region. Their base composition, genetic distance and tRNA structure indeed revealed a wide differentiation between mitogenomes. Two P. tunguidus from Guangxi near Vietnam differed from the other seven by up to 10.1%. Their tRNA-Arg had a complete clover-leaf structure, whereas that of the others did not contain an entire dihydrouridine arm. The nine mitogenomes also differed in the length of rRNA. NJ, ML, and Bayesian analyses all split them into two clades, viz. the two P. tunguidus from Guangxi (Clade 1), and the other seven (Clade 2). Both the structure and phylogeny of the mitogenomes suggest that P. tunguidus has complex geographic origin, and its populations in Clade 1 have long lived in isolation from those in Clade 2. They currently reach the level of subspecies or cryptic species. An extensive phylogenetic analysis of Copepoda further verified that Diaptomidae is the most recently diverging family in Calanoida and that P. tunguidus is at the evolutionary apex of the family.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Complete Mitochondrial Genome of Pseudocaranx dentex (Carangidae, Perciformes) Provides Insight into Phylogenetic and Evolutionary Relationship among Carangidae Family

Pseudocaranx dentex (white trevally) which belongs to the Carangidae family, is an important commercial fishery and aquaculture resource in Asia. However, its evolution and population genetics have received little attention which was limited by the mitogenome information absence. Here, we sequenced and annotated the complete mitochondrial genome of P. dentex which was 16,569 bp in length, containing twenty-two tRNAs (transfer RNAs), thirteen PCGs (protein-coding genes), two rRNAs (ribosomal RNAs), and one non-coding region with conservative gene arrangement. The Ka/Ks ratio analysis among Carangidae fishes indicated the PCGs were suffering purify selection and the values were related to the taxonomic status and further influenced by their living habits. Phylogenetic analysis based on the PCGs sequences of mitogenomes among 36 species presented three major clades in Carangidae. According to the phylogenetic tree, we further analyzed the taxonomic confusion of Carangoides equula which was on the same branch with P. dentex but a different branch with Carangoides spp. We inferred Kaiwarinus equula should be the accepted name and belong to the independent Kaiwarinus genus which was the sister genus of Pseudocaranx. This work provides mitochondrial genetic information and verifies the taxonomic status of P. dentex, and further helps to recognize the phylogenetic relationship and evolutionary history of Carangidae.

The first two complete mitogenomes of the order Apodida from deep-sea chemoautotrophic environments: New insights into the gene rearrangement, origin and evolution of the deep-sea sea cucumbers

The deep-sea ecosystem is considered as the largest and most remote biome of the world. It is meaningful and important to elucidate the life origins by exploring the origin and adaptive genetic mechanisms of the large deep-sea organisms. Sea cucumbers (Holothuroidea) are abundant and economically important group of echinoderms, living from the shallow-waters to deep-sea. In this study, we present the mitochondrial genomes of the sea cucumber Chiridota heheva and Chiridota sp. collected from the deep-sea cold seep and hydrothermal vent, respectively. This is the first reported mitochondrial genomes from the order Apodida. The mitochondrial genomes of C. heheva (17,200 bp) and Chiridota sp. (17,199 bp) display novel gene arrangements with the first protein-coding gene rearrangements in the class Holothuroidea. Bases composition analysis showed that the A + T content of deep-sea holothurians were significantly higher than that of the shallow-water groups. We compared the arrangement of genes from the 24 available holothurian mitogenomes and found that the transposition, reverse transposition and tandem-duplication-random-losses (TDRL) may be involved in the evolution of mitochondrial gene arrangements in Holothuroidea. Phylogenetic analysis revealed that the Apodida clustered with Elasipodida, forming two basal deep-sea holothurian clades. The divergence between the deep-sea and shallow-water holothurians was located at 386.93 Mya, during the Late Devonian. Mitochondrial protein-coding genes of deep-sea holothurians underwent relaxed purifying selection. There are 57 positive selected amino acids sites for some mitochondrial genes of the three deep-sea clades, implying they may involve in the adaption of deep-sea sea cucumbers.

The complete mitochondrial genome of the intertidal spider (Desis jiaxiangi) provides novel insights into the adaptive evolution of the mitogenome and the evolution of spiders

BACKGROUND

Although almost all extant spider species live in terrestrial environments, a few species live fully submerged in freshwater or seawater. The intertidal spiders (genus Desis) built silk nests within coral crevices can survive submerged in high tides. The diving bell spider, Argyroneta aquatica, resides in a similar dynamic environment but exclusively in freshwater. Given the pivotal role played by mitochondria in supplying most energy for physiological activity via oxidative phosphorylation and the environment, herein we sequenced the complete mitogenome of Desis jiaxiangi to investigate the adaptive evolution of the aquatic spider mitogenomes and the evolution of spiders.

RESULTS

We assembled a complete mitogenome of the intertidal spider Desis jiaxiangi and performed comparative mitochondrial analyses of data set comprising of Desis jiaxiangi and other 45 previously published spider mitogenome sequences, including that of Argyroneta aquatica. We found a unique transposition of trnL2 and trnN genes in Desis jiaxiangi. Our robust phylogenetic topology clearly deciphered the evolutionary relationships between Desis jiaxiangi and Argyroneta aquatica as well as other spiders. We dated the divergence of Desis jiaxiangi and Argyroneta aquatica to the late Cretaceous at ~ 98 Ma. Our selection analyses detected a positive selection signal in the nd4 gene of the aquatic branch comprising both Desis jiaxiangi and Argyroneta aquatica. Surprisingly, Pirata subpiraticus, Hypochilus thorelli, and Argyroneta aquatica each had a higher Ka/Ks value in the 13 PCGs dataset among 46 taxa with complete mitogenomes, and these three species also showed positive selection signal in the nd6 gene.

CONCLUSIONS

Our finding of the unique transposition of trnL2 and trnN genes indicates that these genes may have experienced rearrangements in the history of intertidal spider evolution. The positive selection signals in the nd4 and nd6 genes might enable a better understanding of the spider metabolic adaptations in relation to different environments. Our construction of a novel mitogenome for the intertidal spider thus sheds light on the evolutionary history of spiders and their mitogenomes.

An Investigation of the Variations in Complete Mitochondrial Genomes of Lingula anatina in the Western Pacific Region

Lingula anatina is a brachiopod widely distributed in the western Pacific region. Even though L. anatina has been targeted for a number of biological studies, there is still limited information on intraspecific genetic variations of L. anatina. In this study, L. anatina specimens were collected from Korea and Vietnam, and complete mitochondrial genome (mitogenome) sequences were analyzed and compared with previous records. The total mitogenomes of L. anatina were 24,875 bp and 25,305 bp in size for Korean and Vietnamese specimens, respectively. Those mitogenomes are extraordinarily longer than the typical mitogenome size for an animal but shorter than the previous record from Yanagawa (Japan) for this species. The gene orders and the sizes of the protein-coding genes are also different from those for the Japanese specimen. Furthermore, the nonsynonymous (Ka) and synonymous (Ks) substitution rates in protein-coding genes (PCGs) were calculated to test the idea of evolutionary rate differences in mitochondrial genomes. The analyses showed relatively low Ka and Ks for the complete mitogenomes from Buan (Korea), Doson (Vietnam) and Yanagawa (Japan). The Ka/Ks ratio was less than 1 in comparisons of three localities, indicating the existence of purifying selection in this species. The phylogenetic analyses showed that L. anatina diverged among localities in the western Pacific region.

Patterns of Natural Selection on Mitochondrial Protein-Coding Genes in Lungless Salamanders: Relaxed Purifying Selection and Presence of Positively Selected Codon Sites in the Family Plethodontidae

There are two distinct lungless groups in caudate amphibians (salamanders and newts) (the family Plethodontidae and the genus Onychodactylus, from the family Hynobiidae). Lunglessness is considered to have evolved in response to environmental and/or ecological adaptation with respect to oxygen requirements. We performed selection analyses on lungless salamanders to elucidate the selective patterns of mitochondrial protein-coding genes associated with lunglessness. The branch model and RELAX analyses revealed the occurrence of relaxed selection (an increase of the dN/dS ratio = ω value) in most mitochondrial protein-coding genes of plethodontid salamander branches but not in those of Onychodactylus. Additional branch model and RELAX analyses indicated that direct-developing plethodontids showed the relaxed pattern for most mitochondrial genes, although metamorphosing plethodontids had fewer relaxed genes. Furthermore, aBSREL analysis detected positively selected codons in three plethodontid branches but not in Onychodactylus. One of these three branches corresponded to the most recent common ancestor, and the others corresponded with the most recent common ancestors of direct-developing branches within Hemidactyliinae. The positive selection of mitochondrial protein-coding genes in Plethodontidae is probably associated with the evolution of direct development.

Not Frozen in the Ice: Large and Dynamic Rearrangements in the Mitochondrial Genomes of the Antarctic Fish

The vertebrate mitochondrial genomes generally present a typical gene order. Exceptions are uncommon and important to study the genetic mechanisms of gene order rearrangements and their consequences on phylogenetic output and mitochondrial function. Antarctic notothenioid fish carry some peculiar rearrangements of the mitochondrial gene order. In this first systematic study of 28 species, we analyzed known and undescribed mitochondrial genome rearrangements for a total of eight different gene orders within the notothenioid fish. Our reconstructions suggest that transpositions, duplications, and inversion of multiple genes are the most likely mechanisms of rearrangement in notothenioid mitochondrial genomes. In Trematominae, we documented an extremely rare inversion of a large genomic segment of 5,300 bp that partially affected the gene compositional bias but not the phylogenetic output. The genomic region delimited by nad5 and trnF, close to the area of the Control Region, was identified as the hot spot of variation in Antarctic fish mitochondrial genomes. Analyzing the sequence of several intergenic spacers and mapping the arrangements on a newly generated phylogeny showed that the entire history of the Antarctic notothenioids is characterized by multiple, relatively rapid, events of disruption of the gene order. We hypothesized that a pre-existing genomic flexibility of the ancestor of the Antarctic notothenioids may have generated a precondition for gene order rearrangement, and the pressure of purifying selection could have worked for a rapid restoration of the mitochondrial functionality and compactness after each event of rearrangement.

A non-exhaustive survey revealed possible genetic similarity in mitochondrial adaptive evolution of marine fish species in the northwestern Pacific

Mitochondrial coding genes involved in the oxidative phosphorylation pathway play vitally important roles in energy production and thermal adaptation. Investigating the underlying molecular mechanism of mitochondrial adaptive evolution is crucial for understanding biodiversity and ecological radiation. In this study, we collated population genetic studies of marine fish species in the northwestern Pacific based on mitochondrial cytochrome b gene sequences, to investigate whether similar patterns could be detected in mitochondrial adaptive evolution. After filtering, nine studies containing eight marine fish species (Ammodytespersonatus, Boleophthalmuspectinirostris, Larimichthyspolyactis, Mugilcephalus, Pampusargenteus, Platycephalus sp.1, Sebastiscusmarmoratus, and Trachidermusfasciatus) belonging to eight different families were retained. Multiple codon-based approaches were used to identify potential sites under selection in each species. By comparison, our results showed that the posterior part of the mitochondrial cytochrome b gene (particularly codon 372 and its neighboring sites) seemed to be involved in the adaptive evolution process, suggesting potential genetic similarity among distantly related species. We also summarized four types of adaptive patterns in the reviewed species, and suggest that the level of genetic differentiation and mitochondrial adaptive evolution might be correlated. Further studies are needed to confirm such relationship by detecting RNA-level evidence and investigating more species and samples.

Genetic analysis of whole mitochondrial genome of Lateolabrax maculatus (Perciformes: Moronidae) indicates the presence of two populations along the Chinese coast

The whole mitochondrial genome ofLateolabrax maculatus(Cuvier, 1828) was used to investigate the reasons for the observed patterns of genetic differentiation among 12 populations in northern and southern China. The haplotype diversity and nucleotide diversity ofL. maculatuswere 0.998 and 0.00169, respectively. Pairwise FSTvalues between populations ranged from 0.001 to 0.429, correlating positively with geographic distance. Genetic structure analysis and haplotype network analysis indicated that these populations were split into two groups, in agreement with geographic segregation and environment. Tajima’s D values, Fu’s Fs tests and Bayesian skyline plot (BSP) indicated that a demographic expansion event may have occurred in the history ofL. maculatus. Through selection pressure analysis, we found evidence of significant negative selection at the ATP6, ND3, Cytb, COX3, COX2 and COX1 genes. In our hypotheses, this study implied that demographic events and selection of local environmental conditions, including temperature, are responsible for population divergence. These findings are a step forward toward the understanding of the genetic basis of differentiation and adaptation, as well as conservation ofL. maculatus.

Mitogenomic Perspectives on the Adaptation to Extreme Alkaline Environment of Amur ide (Leuciscus waleckii)

Amur ide (Leuciscus waleckii, Family Cyprinidae) is widely distributed in Northeast Asia. L. waleckii usually inhabits freshwater environments but can also survive in the Lake Dali Nur, one of the most extreme aquatic environments on the earth, with an alkalinity up to 50 mmol/L (pH 9.6). To investigate mechanisms of mitogenomic evolution underlying adaptation to extreme environments, we determined 30 complete mitogenomes that included Lake Dali Nur (alkaline environment, AL) population and Amur basin (freshwater environment, FW) population. Through phylogenetic and divergence time analysis, we found that AL and FW populations forming distinct two groups which were consistent with geographic divergence (the formation of Lake Dali Nur). In addition, we found that almost of the windows exhibited higher nucleotide diversity in FW population (avg 0.0046) than AL population (avg 0.0012). This result indicated that severe environment selection had remarkably reduced the genetic diversity of mitogenome in AL population and suggested that severe environment selection had remarkably reduced the genetic diversity of mitogenome in the AL population. Compared with the FW population (ω = 0.064), the AL population (ω = 0.092) had a larger mean ω (dN/dS ratios) value for the 13 concatenated mitochondrial protein-coding genes, indicating that the high alkaline tolerated group had accumulated more nonsynonymous mutations. These nonsynonymous mutations had resulted in slightly beneficial amino acid changes that allowed adaption to the severe conditions. This study provides an additional view to decipher the adaptive mitogenome evolution of L. waleckii of the high alkaline environment.

Phylogenetic perspective on the relationships and evolutionary history of the Acipenseriformes

The Acipenseriformes, as one of the earliest extant vertebrates, plays an important role in the evolution of fishes and even the whole vertebrates. Here we collected and analyzed all complete mitochondrial genomes of Acipenseriformes species. Phylogenetic analyses demonstrated that the polytomous branch included Acipenseridae and Polyodontidae formed five clades. The Polyodontidae clade and the Scaphirhynchus clade both were monophyletic group, whereas the Acipenser species and the Huso species both were polyphyletic group. The Bayesian divergence times showed that the origin time for Acipenseriformes was at 318.0 Mya, which was similar to the some previous results of 312.1 Mya, 346.9 Mya and 389.7 Mya. The result was in good consistent with the paleontological data available and the split time of the Pacific and Atlantic Oceans from the Jurassic to the Cretaceous (Laurasia splits in North America and Eurasia). The dN/dS ratios showed the evolutionary rates gradually slow down in five major Acipenseriformes clades from the Clade A (the Pacific sturgeons species) to Clade C (the genus Scaphirhynchus), which was related to the process of geographical formation.

Evolutionary rates of and selective constraints on the mitochondrial genomes of Orthoptera insects with different wing types

Orthoptera is the most diverse order of polyneopterans, and the forewing and hindwing of its members exhibit extremely variability from full length to complete loss in many groups; thus, this order provides a good model for studying the effects of insect flight ability on the evolutionary constraints on and evolutionary rate of the mitochondrial genome. Based on a data set of mitochondrial genomes from 171 species, including 43 newly determined, we reconstructed Orthoptera phylogenetic relationships and estimated the divergence times of this group. The results supported Caelifera and Ensifera as two monophyletic groups, and revealed that Orthoptera originated in the Carboniferous (298.997 Mya). The date of divergence between the suborders Caelifera and Ensifera was 255.705 Mya, in the late Permian. The major lineages of Acrididae seemed to have radiated in the Cenozoic, and the six patterns of rearrangement of 171 Orthoptera mitogenomes mostly occurred in the Cretaceous and Cenozoic. Based on phylogenetic relationships and ancestral state reconstruction, we analysed the evolutionary selection pressure on and evolutionary rate of mitochondrial protein-coding genes (mPCGs). The results indicated that during approximately 300 Mya of evolution, these genes experienced purifying selection to maintain their function. Flightless orthopteran insects accumulated more non-synonymous mutations than flying species and experienced more relaxed evolutionary constraints. The different wing types had different evolutionary rates, and the mean evolutionary rate of Orthoptera mitochondrial mPCGs was 13.554 × 10^-9 subs/s/y. The differences in selection pressures and evolutionary rates observed between the mitochondrial genomes suggested that functional constraints due to locomotion play an important role in the evolution of mitochondrial DNA in orthopteran insects with different wing types.

Complete mitochondrial genomes of eight seahorses and pipefishes (Syngnathiformes: Syngnathidae): insight into the adaptive radiation of syngnathid fishes

Background

The evolution of male pregnancy is the most distinctive characteristic of syngnathids, and their specialized life history traits make syngnathid species excellent model species for many issues in biological evolution. However, the origin of syngnathids and the evolutionary divergence time of different syngnathid species remain poorly resolved. Comprehensive phylogenetic studies of the Syngnathidae will provide critical evidence to elucidate their origin, evolution, and dispersal patterns.

Results

We sequenced the mitochondrial genomes of eight syngnathid species in this study, and the estimated divergence times suggested that syngnathids diverged from other teleosts approximately 48.8 Mya during the Eocene period. Selection analysis showed that many mitochondrial genes of syngnathids exhibited significantly lower Ka/Ks values than those of other teleosts. The two most frequently used codons in syngnathid fishes were different from those in other teleosts, and a greater proportion of the mitochondrial simple sequence repeats (SSRs) were distributed in non-coding sequences in syngnathids compared with other teleosts.

Conclusions

Our study indicated that syngnathid fishes experienced an adaptive radiation process during the early explosion of species. Syngnathid mitochondrial OXPHOS genes appear to exhibit depressed Ka/Ks ratios compared with those of other teleosts, and this may suggest that their mitogenomes have experienced strong selective constraints to eliminate deleterious mutations.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1430-3) contains supplementary material, which is available to authorized users.

Evidence for relaxed selection of mitogenome in rapid-flow cyprinids

BACKGROUND

Hypoxia adaptation is developed in many fish species, which helped them to habitat most of water bodies. However, fishes living under high oxygen concentration may lose this feature. Rapid flows provide high level and stable dissolved oxygen, which facilitate organism's oxygen supply and energy production. Previous studies showed that fish species from rapid-flow habitats exhibited lower hypoxia tolerance compared with fish from intermediate- and slow-flow habitats. Mitochondrial genomes code 13 key components in oxidative phosphorylation pathway; these genes may be under relaxed selection in rapid-flow species.

OBJECTIVES

The primary objectives of this study is to investigate the evolutionary patterns of the 13 mitochondrial OXPHOS genes among nine cyprinids from different water bodies and to test the hypotheses that mitochondrial OXPHOS genes may experience relaxed selection in rapid-flow habitats.

METHODS

We classified nine cyprinid fish species into three groups based on their habitats: rapid-flow, intermediate-flow and slow-flow. To detect relaxed selections, we investigated the 13 protein-coding genes with codon evolution programs RELAX; to estimate evolutionary rates among the cyprinids, free-ratio model in Codeml program was applied; Branch-site models were applied to detect positive selection sites. The polymorphisms of homologous sites were evaluated with PROVEAN program and projected to 3D structure prediction of the proteins using SWISS-MODEL.

RESULTS

We found that nine out of the 13 genes are under relaxed selection in rapid-flow species. Furthermore, dN, dS and dN/dS are relatively increased when compared with those of intermediate-flow species. More amino acid polymorphic sites are presented in rapid-flow species than in intermediate- and slow-flow species. Furthermore, rapid-flow species had more deleterious substitutions than other groups. 3D structure prediction of these proteins and projection of the polymorphic sites indicated that these sites were randomly distributed, suggesting relaxed functional constraints of these proteins in rapid-flow species.

CONCLUSION

Our results suggest that mitochondrial genes are under relaxed selection in rapid-flow cyprinids.

Convergent Evolution of Mitochondrial Genes in Deep-Sea Fishes

Deep seas have extremely harsh conditions including high hydrostatic pressure, total darkness, cold, and little food and oxygen. The adaptations of fishes to deep-sea environment apparently have occurred independently many times. The genetic basis of adaptation for obtaining their energy remains unknown. Mitochondria play a central role in aerobic respiration. Analyses of the available 2,161 complete mitochondrial genomes of 1,042 fishes, including 115 deep-sea species, detect signals of positive selection in mitochondrial genes in nine branches of deep-sea fishes. Aerobic metabolism yields much more energy per unit of source material than anaerobic metabolism. The adaptive evolution of the mtDNA may reflect that aerobic metabolism plays a more important role than anaerobic metabolism in deep-sea fishes, whose energy sources (food) are extremely limited. This strategy maximizes the usage of energy sources. Eleven mitochondrial genes have convergent/parallel amino acid changes between branches of deep-sea fishes. Thus, these amino acid sites may be functionally important in the acquisition of energy, and reflect convergent evolution during their independent invasion of the harsh deep-sea ecological niche.

Relation between mitochondrial DNA hyperdiversity, mutation rate and mitochondrial genome evolution in Melarhaphe neritoides (Gastropoda: Littorinidae) and other Caenogastropoda

Mitochondrial DNA hyperdiversity is primarily caused by high mutation rates (µ) and has potential implications for mitogenome architecture and evolution. In the hyperdiverse mtDNA of Melarhaphe neritoides (Gastropoda: Littorinidae), high mutational pressure generates unusually large amounts of synonymous variation, which is expected to (1) promote changes in synonymous codon usage, (2) reflect selection at synonymous sites, (3) increase mtDNA recombination and gene rearrangement, and (4) be correlated with high mtDNA substitution rates. The mitogenome of M. neritoides was sequenced, compared to closely related littorinids and put in the phylogenetic context of Caenogastropoda, to assess the influence of mtDNA hyperdiversity and high µ on gene content and gene order. Most mitogenome features are in line with the trend in Mollusca, except for the atypical secondary structure of the methionine transfer RNA lacking the TΨC-loop. Therefore, mtDNA hyperdiversity and high µ in M. neritoides do not seem to affect its mitogenome architecture. Synonymous sites are under positive selection, which adds to the growing evidence of non-neutral evolution at synonymous sites. Under such non-neutrality, substitution rate involves neutral and non-neutral substitutions, and high µ is not necessarily associated with high substitution rate, thus explaining that, unlike high µ, a high substitution rate is associated with gene order rearrangement.

Evolutionary divergence of mitochondrial genomes in two Tetranychus species distributed across different climates

There is increasing evidence that mitochondrial genomes (mitogenomes) can be under selection, whereas the selective regimes shaping mitogenome evolution remain largely unclear. To test for mitogenome evolution in relation to the climate adaptation, we explored mtDNA variation in two spider mite (Tetranychus) species that distribute across different climates. We sequenced 26 complete mitogenomes of Tetranychus truncates, which occurs in both warm and cold regions, and nine complete mitogenomes of Tetranychus pueraricola, which is restricted to warm regions. Patterns of evolution in the two species' mitogenomes were compared through a series of d_N /d_S methods and physicochemical profiles of amino acid replacements. We found that: (1) the mitogenomes of both species were under widespread purifying selection; (2) elevated directional adaptive selection was observed in the T. truncatus mitogenome, perhaps linked to the cold climates adaptation of T. truncatus; and (3) the strength of selection varied across genes, and diversifying positive selection detected on ND4 and ATP6 pointed to their crucial roles during adaptation to different climatic conditions. This study gained insight into the mitogenome evolution in relation to the climate adaptation.

Evidence for intra-mitochondrial variation in population genetic structure of Platycephalus sp.1 In the Northwestern Pacific

Understanding evolutionary dynamics of population structure and genetic diversity of marine fish species is of great importance for fishery management and resource conservation. Platycephalus sp.1 is a commercially important marine fish in the Northwestern Pacific. Yet, current understanding of population genetic patterns of Platycephalus sp.1 remains limited. In this study, the Platycephalus sp.1 individuals from four populations in the Northwestern Pacific, which cover large geographic distances and latitudinal gradients, were sampled to examine genetic diversity and population structure by investigating the mitochondrial sequences of Cytochrome b (CYTB, 1141 bp) and control region (CR, 458 bp), and then to assess intra-mitochondrial genetic variability by using a combination of population- and individual-based analytical approaches. Moderate to low levels of genetic diversity were observed in sampled populations by investigating both CR and CYTB datasets. In CR dataset, weak genetic differentiation among populations were observed and population Tokyo Bay was significantly differentiated with other sampled populations. In CYTB dataset, lower genetic differentiations were observed and population Beihai (BH) was significantly differentiated. A non-synonymous substitution was predominantly detected in population BH, and such substitution was tested as a significantly positive selective site in natural selection tests. Based on these findings, we propose that local adaptation influenced by natural selection contributes largely to the intra-mitochondrial variation in Platycephalus sp.1 populations. The present study provides new perspectives on the population genetic structure of Platycephalus sp.1, which could have important implications for the sound management and conservation of this fishery species.

Selective constraints in cold-region wild boars may defuse the effects of small effective population size on molecular evolution of mitogenomes

Spatial range expansion during population colonization is characterized by demographic events that may have significant effects on the efficiency of natural selection. Population genetics suggests that genetic drift brought by small effective population size (Ne) may undermine the efficiency of selection, leading to a faster accumulation of nonsynonymous mutations. However, it is still unknown whether this effect might be balanced or even reversed by strong selective constraints. Here, we used wild boars and local domestic pigs from tropical (Vietnam) and subarctic region (Siberia) as animal model to evaluate the effects of functional constraints and genetic drift on shaping molecular evolution. The likelihood-ratio test revealed that Siberian clade evolved significantly different from Vietnamese clades. Different datasets consistently showed that Siberian wild boars had lower Ka/Ks ratios than Vietnamese samples. The potential role of positive selection for branches with higher Ka/Ks was evaluated using branch-site model comparison. No signal of positive selection was found for the higher Ka/Ks in Vietnamese clades, suggesting the interclade difference was mainly due to the reduction in Ka/Ks for Siberian samples. This conclusion was further confirmed by the result from a larger sample size, among which wild boars from northern Asia (subarctic and nearby region) had lower Ka/Ks than those from southern Asia (temperate and tropical region). The lower Ka/Ks might be due to either stronger functional constraints, which prevent nonsynonymous mutations from accumulating in subarctic wild boars, or larger Ne in Siberian wild boars, which can boost the efficacy of purifying selection to remove functional mutations. The latter possibility was further ruled out by the Bayesian skyline plot analysis, which revealed that historical Ne of Siberian wild boars was smaller than that of Vietnamese wild boars. Altogether, these results suggest stronger functional constraints acting on mitogenomes of subarctic wild boars, which may provide new insights into their local adaptation of cold resistance.

The Complete Mitochondrial Genome of Glyptothorax macromaculatus Provides a Well-Resolved Molecular Phylogeny of the Chinese Sisorid Catfishes

Previous phylogenetic analyses of the Chinese sisorid catfishes have either been poorly resolved or have not included all the 12 sisorid genera. Here, we successfully assembled the first complete mitochondrial genome of the sisorid fish Glyptothorax macromaculatus. Based on this novel mitochondrial genome and previously published mitochondrial genomes in the Sisoridae, we generated maximum likelihood and Bayesian phylogenies. We dated our preferred topology using fossil calibration points. We also tested the protein-coding genes in the mitochondrial genomes of the glyptosternoid fishes for signals of natural selection by comparing the nucleotide substitution rate along the branch ancestral to the glyptosternoid fishes to other branches in our topology. The mitochondrial sequence structure of G. macromaculatus was similar to those known from other vertebrates, with some slight differences. Our sisorid phylogenies were well-resolved and well-supported, with exact congruence between the different phylogenetic methods. This robust phylogeny clarified the relationships among the Chinese sisorid genera and strongly supported the division of the family into three main clades. Interestingly, the glyptosternoid divergence time predicted by our molecular dating analysis coincided with the uplift of the Tibetan Plateau, suggesting that geology may have influenced speciation in the Sisoridae. Among the mitochondrial protein-coding genes, atp8 may have most rapidly evolved, and atp6 may have been subjected to positive selection pressure to adapt to high elevations. In summary, this study provided novel insights into the phylogeny, evolution and high-altitude adaptions of the Chinese sisorid fishes.

Phylogeny, evolution and mitochondrial gene order rearrangement in scale worms (Aphroditiformia, Annelida)

Next-generation sequencing (NGS) has become a powerful tool in phylogenetic and evolutionary studies. Here we applied NGS to recover two ribosomal RNA genes (18S and 28S) from 16 species and 15 mitochondrial genomes from 16 species of scale worms representing six families in the suborder Aphroditiformia (Phyllodocida, Annelida), a complex group of polychaetes characterized by the presence of dorsal elytra or scales. The phylogenetic relationship of the several groups of scale worms remains unresolved due to insufficient taxon sampling and low resolution of individual gene markers. Phylogenetic tree topology based on mitochondrial genomes is comparable with that based on concatenated sequences from two mitochondrial genes (cox1 and 16S) and two ribosomal genes (18S and 28S) genes, but has higher statistical support for several clades. Our analyses show that Aphroditiformia is monophyletic, indicating the presence of elytra is an apomorphic trait. Eulepethidae and Aphroditidae together form the sister group to all other families in this suborder, whereas Acoetidae is sister to Iphionidae. Polynoidae is monophyletic, but within this family the deep-sea subfamilies Branchinotogluminae and Macellicephalinae are paraphyletic. Mitochondrial genomes in most scale-worm families have a conserved gene order, but within Polynoidae there are two novel arrangement patterns in the deep-sea clade. Mitochondrial protein-coding genes in polynoids as a whole have evolved under strong purifying selection, but substitution rates in deep-sea species are much higher than those in shallow-water species, indicating that purifying selection is relaxed in deep-sea polynoids. There are positive selected amino acids for some mitochondrial genes of the deep-sea clade, indicating they may involve in the adaption of deep-sea polynoids. Overall, our study (1) provided more evidence for reconstruction of the phylogeny of Aphroditiformia, (2) provided evidence to refute the assumption that mitochondrial gene order in Errantia is conserved, and (3) indicated that the deep-sea extreme environment may have affected the mitochondrial genome evolution rate and gene order arrangement in Polynoidae.

Limited locomotive ability relaxed selective constraints on molluscs mitochondrial genomes

Mollusca are the second largest phylum in the animal kingdom with different types of locomotion. Some molluscs are poor-migrating, while others are free-moving or fast-swimming. Most of the energy required for locomotion is provided by mitochondria via oxidative phosphorylation. Here, we conduct a comparative genomic analysis of 256 molluscs complete mitochondrial genomes and evaluate the role of energetic functional constraints on the protein-coding genes, providing a new insight into mitochondrial DNA (mtDNA) evolution. The weakly locomotive molluscs, compared to strongly locomotive molluscs, show significantly higher Ka/Ks ratio, which suggest they accumulated more nonsynonymous mutations in mtDNA and have experienced more relaxed evolutionary constraints. Eleven protein-coding genes (CoxI, CoxII, ATP6, Cytb, ND1-6, ND4L) show significant difference for Ka/Ks ratios between the strongly and weakly locomotive groups. The relaxation of selective constraints on Atp8 arise in the common ancestor of bivalves, and the further relaxation occurred in marine bivalves lineage. Our study thus demonstrates that selective constraints relevant to locomotive ability play an essential role in evolution of molluscs mtDNA.

The first complete mitogenome of the South China deep-sea giant isopod Bathynomus sp. (Crustacea: Isopoda: Cirolanidae) allows insights into the early mitogenomic evolution of isopods

In this study, the complete mitochondrial (mt) genome sequence of the South China deep‐sea giant isopod Bathynomus sp. was determined, and this study is the first to explore in detail the mt genome of a deep‐sea member of the order Isopoda. This species belongs to the genus Bathynomus, the members of which are saprophagous residents of the deep‐sea benthic environment; based on their large size, Bathynomus is included in the “supergiant group” of isopods. The mt genome of Bathynomus sp. is 14,965 bp in length and consists of 13 protein‐coding genes, two ribosomal RNA genes, only 18 transfer RNA genes, and a noncoding control region 362 bp in length, which is the smallest control region discovered in Isopoda to date. Although the overall genome organization is typical for metazoans, the mt genome of Bathynomus sp. shows a number of derived characters, such as an inversion of 10 genes when compared to the pancrustacean ground pattern. Rearrangements in some genes (e.g., cob, trnT, nad5, and trnF) are shared by nearly all isopod mt genomes analyzed thus far, and when compared to the putative isopod ground pattern, five rearrangements were found in Bathynomus sp. Two tRNAs exhibit modified secondary structures: The TΨC arm is absent from trnQ, and trnC lacks the DHU. Within the class Malacostraca, trnC arm loss is only found in other isopods. Phylogenetic analysis revealed that Bathynomus sp. (Cymothoida) and Sphaeroma serratum (Sphaeromatidea) form a single clade, although it is unclear whether Cymothoida is monophyletic or paraphyletic. Moreover, the evolutionary rate of Bathynomus sp. (dN/dS [nonsynonymous mutational rate/synonymous mutational rate] = 0.0705) is the slowest measured to date among Cymothoida, which may be associated with its relatively constant deep‐sea environment. Overall, our results may provide useful information for understanding the evolution of deep‐sea Isopoda species.

Comparative mitogenomic analysis of three species of periwinkles: Littorina fabalis, L. obtusata and L. saxatilis

The flat periwinkles, Littorina fabalis and L. obtusata, offer an interesting system for local adaptation and ecological speciation studies. In order to provide genomic resources for these species, we sequenced their mitogenomes together with that of the rough periwinkle L. saxatilis by means of next-generation sequencing technologies. The three mitogenomes present the typical repertoire of 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes and a putative control region. Although the latter could not be fully recovered in flat periwinkles using short-reads due to a highly repetitive fragment, in L. saxatilis this problem was overcome with additional long-reads and we were able to assemble the complete mitogenome. Both gene order and nucleotide composition are similar between the three species as well as compared to other Littorinimorpha. A large variance in divergence was observed across mitochondrial regions, with six- to ten-fold difference between the highest and the lowest divergence rates. Based on nucleotide changes on the whole molecule and assuming a molecular clock, L. fabalis and L. obtusata started to diverge around 0.8 Mya (0.4-1.1 Mya). The evolution of the mitochondrial protein-coding genes in the three Littorina species appears mainly influenced by purifying selection as revealed by phylogenetic tests based on d_N/d_S ratios that did not detect any evidence for positive selection, although some caution is required given the limited power of the dataset and the implemented approaches.

Mitogenomic perspectives on the origin of Tibetan loaches and their adaptation to high altitude

Tibetan loaches are the largest group of Tibetan fishes and are well adapted to the Tibetan Plateau. To investigate the origin of Tibetan loaches and their adaptations to the Tibetan Plateau, we determined 32 complete mitochondrial genomes that included 29 Tibetan loach species, two Barbatula species and Schistura longus. By combining these newly determined sequences with other previously published mitochondrial genomes, we assembled a large mitogenomic data set (11,433 bp) of 96 species in the superfamily Cobitoidea, to investigate the phylogenetic status of the genus Triplophysa. The resulting phylogeny strongly supported that the genus Triplophysa forms a monophyletic group within Nemacheilidae. Our molecular dating time suggests that the lineage leading to the Tibetan loaches and other loaches diverged approximately 23.5 Ma, which falls within the period of recent major uplifts of the Tibetan Plateau in the Early Miocene. Selection analyses revealed that the mitochondrial protein-coding genes of Tibetan loaches have larger ratios of nonsynonymous to synonymous substitutions than do those of non-Tibetan loaches, indicating that Tibetan loaches accumulated more nonsynonymous mutations than non-Tibetan loaches and exhibited rapid evolution. Two positively selected sites were identified in the ATP8 and ND1 genes.

Comparative analysis of complete mitochondrial genomes suggests that relaxed purifying selection is driving high nonsynonymous evolutionary rate of the NADH2 gene in whitefish (Coregonus ssp.)

Several studies have recently reported evidence for positive selection acting on the mitochondrial genome (mitogenome), emphasizing its potential role in adaptive divergence and speciation. In this study we searched 107 full mitogenomes of recently diverged species and lineages of whitefish (Coregonus ssp.) for signals of positive selection. These salmonids show several distinct morphological and ecological differences that may be associated with energetics and therefore potentially positive selection at the mitogenome level. We found that purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, the NADH dehydrogenase 2 gene (ND2) showed a highly elevated dN/dS ratio compared to the other mitochondrial genes, which was significantly higher in whitefish compared to other salmonids. We therefore further examined nonsynonymous evolution in ND2 by (i) mapping amino acid changes to a protein model structure which showed that they were located away from key functional residues of the protein, (ii) locating them in the sequences of other species of fish (Salmonidae, Anguillidae, Scombridae and Percidae) only to find pronounced overlap of nonsynonymous regions. We thus conclude that relaxed purifying selection is driving the evolution of ND2 by affecting mostly regions that have lower functional relevance.