ND6 gene "lost" and found: evolution of mitochondrial gene rearrangement in Antarctic notothenioids

Mitoribosome structure with cofactors and modifications reveals mechanism of ligand binding and interactions with L1 stalk

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNAVal. The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed us to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transitions in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide a description of the structure and function of the human mitoribosome.

Characterization of the mitochondrial Huso huso genome and new aspects of its organization in the presence of tandem repeats in 12S rRNA

BACKGROUND

The sturgeon group has been economically significant worldwide due to caviar production. Sturgeons consist of 27 species in the world. Mitogenome data could be used to infer genetic diversity and investigate the evolutionary history of sturgeons. A limited number of complete mitogenomes in this family were sequenced. Here, we annotated the mitochondrial Huso huso genome, which revealed new aspects of this species.

RESULTS

In this species, the mitochondrial genome consisted of 13 genes encoding proteins, 22tRNA and 2rRNA, and two non-coding regions that followed other vertebrates. In addition, H. huso had a pseudo-tRNA-Glu between ND6 and Cytb and a 52-nucleotide tandem repeat with two replications in 12S rRNA. This duplication event is probably related to the slipped strand during replication, which could remain in the strand due to mispairing during replication. Furthermore, an 82 bp repeat sequence with three replications was observed in the D-loop control region, which is usually visible in different species. Regulatory elements were also seen in the control region of the mitochondrial genome, which included termination sequences and conserved regulatory blocks. Genomic compounds showed the highest conservation in rRNA and tRNA, while protein-encoded genes and nonencoded regions had the highest divergence. The mitochondrial genome was phylogenetically assayed using 12 protein-encoding genes.

CONCLUSIONS

In H. huso sequencing, we identified a distinct genome organization relative to other species that have never been reported. In recent years, along with the advancement in sequencing identified more genome rearrangements. However, it is an essential aspect of researching the evolution of the mitochondrial genome that needs to be recognized.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

Structure of mitoribosome reveals mechanism of mRNA binding, tRNA interactions with L1 stalk, roles of cofactors and rRNA modifications

The mitoribosome translates mitochondrial mRNAs and regulates energy conversion that is a signature of aerobic life forms. We present a 2.2 Å resolution structure of human mitoribosome together with validated mitoribosomal RNA (rRNA) modifications, including aminoacylated CP-tRNA Val . The structure shows how mitoribosomal proteins stabilise binding of mRNA and tRNA helping to align it in the decoding center, whereas the GDP-bound mS29 stabilizes intersubunit communication. Comparison between different states, with respect to tRNA position, allowed to characterize a non-canonical L1 stalk, and molecular dynamics simulations revealed how it facilitates tRNA transition in a way that does not require interactions with rRNA. We also report functionally important polyamines that are depleted when cells are subjected to an antibiotic treatment. The structural, biochemical, and computational data illuminate the principal functional components of the translation mechanism in mitochondria and provide the most complete description so far of the structure and function of the human mitoribosome.

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

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

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

Simple Summary

Duplicated control regions have been reported several times in the tree frog family Rhacophoridae, and previous studies have mostly relied on sequence analysis to reconstruct their evolution. This is the first study to employ a phylogenetic method to demonstrate the existence of concerted and parallel evolution succinctly and intuitively in the duplicated control regions of the family Rhacophoridae. Phylogenetic relationships were also used to illustrate the parallel evolution of ATP8 loss of function in the genus Polypedates. In general, this study elucidated the evolutionary patterns and pathways of mitochondrial gene rearrangement of the family Rhacophoridae from a phylogenetic perspective, which aids in understanding the evolutionary history of this fascinating tree frog taxon from a molecular evolution standpoint.

Abstract

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

A Unique Mitochondrial Gene Block Inversion in Antarctic Trematomin Fishes: A Cautionary Tale

Many Antarctic notothenioid fishes have major rearrangements in their mitochondrial (mt) genomes. Here, we report the complete mt genomes of 3 trematomin notothenioids: the bald notothen (Trematomus (Pagothenia) borchgrevinki), the spotted notothen (T. nicolai), and the emerald notothen (T. bernacchii). The 3 mt genomes were sequenced using next-generation Illumina technology, and the assemblies verified by Sanger sequencing. When compared with the canonical mt gene order of the Antarctic silverfish (Pleuragramma antarctica), we found a large gene inversion in the 3 trematomin mt genomes that included tRNA^Ile, ND1, tRNA^Leu2, 16S, tRNA^Val, 12S, tRNA^Phe, and the control region. The trematomin mt genomes contained 3 intergenic spacers, which are thought to be the remnants of previous gene and control region duplications. All control regions included the characteristic conserved regulatory sequence motifs. Although short-read next-generation DNA sequencing technology has allowed the rapid and cost-effective sequencing of a large number of complete mt genomes, it is essential in all cases to verify the assembly in order to prevent the publication and use of erroneous data.

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.

Foster thy young: enhanced prediction of orphan genes in assembled genomes

Proteins encoded by newly-emerged genes ('orphan genes') share no sequence similarity with proteins in any other species. They provide organisms with a reservoir of genetic elements to quickly respond to changing selection pressures. Here, we systematically assess the ability of five gene prediction pipelines to accurately predict genes in genomes according to phylostratal origin. BRAKER and MAKER are existing, popular ab initio tools that infer gene structures by machine learning. Direct Inference is an evidence-based pipeline we developed to predict gene structures from alignments of RNA-Seq data. The BIND pipeline integrates ab initio predictions of BRAKER and Direct inference; MIND combines Direct Inference and MAKER predictions. We use highly-curated Arabidopsis and yeast annotations as gold-standard benchmarks, and cross-validate in rice. Each pipeline under-predicts orphan genes (as few as 11 percent, under one prediction scenario). Increasing RNA-Seq diversity greatly improves prediction efficacy. The combined methods (BIND and MIND) yield best predictions overall, BIND identifying 68% of annotated orphan genes, 99% of ancient genes, and give the highest sensitivity score regardless dataset in Arabidopsis. We provide a light weight, flexible, reproducible, and well-documented solution to improve gene prediction.

Different gene rearrangements of the genus Dardanus (Anomura: Diogenidae) and insights into the phylogeny of Paguroidea

Complete mitochondrial genomes (mitogenomes) can provide useful information for phylogenetic relationships, gene rearrangement, and molecular evolution. In this study, the complete mitogenomes of two hermit crabs, Dardanus arrosor and Dardanus aspersus, were sequenced for the first time and compared with other published mitogenomes of Paguroidea. Each of the two mitogenomes contains an entire set of 37 genes and a putative control region, but they display different gene arrangements. The different arrangements of the two mitogenomes might be the result of transposition, reversal, and tandem duplication/random loss events from the ancestral pancrustacean pattern. Genome sequence similarity analysis reveals the gene rearrangement in 15 Paguroidea mitogenomes. After synteny analysis between the 15 Paguroidea mitogenomes, an obvious rearranged region is found in D. aspersus mitogenome. Across the 13 protein-coding genes (PCGs) tested, COI has the least and ND6 has the largest genetic distances among the 15 hermit crabs, indicating varied evolution rates of PCGs. In addition, the dN/dS ratio analysis shows that all PCGs are evolving under purifying selection. The phylogenetic analyses based on both gene order and sequence data present the monophyly of three families (Paguridae, Coenobitidae, and Pylochelidae) and the paraphyly of the family Diogenidae. Meanwhile, the phylogenetic tree based on the nucleotide sequences of 13 PCGs shows that two Dardanus species formed a sister group with five Coenobitidae species. These findings help to better understand the gene rearrangement and phylogeny of Paguroidea, as well as provide new insights into the usefulness of mitochondrial gene order as a phylogenetic marker.

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.

Characterization of the mitochondrial genomes of two toads, Anaxyrus americanus (Anura: Bufonidae) and Bufotes pewzowi (Anura: Bufonidae), with phylogenetic and selection pressure analyses

Mitogenomes are useful in analyzing phylogenetic relationships and also appear to influence energy metabolism, thermoregulation and osmoregulation. Much evidence has accumulated for positive selection acting on mitochondrial genes associated with environmental adaptation. Hence, the mitogenome is a likely target for environmental selection. The family Bufonidae (true toads) has only nine complete and four partial mitogenomes published compared to the 610 known species of this family. More mitogenomes are needed in order to obtain a clearer understanding of the phylogenetic relationships within Bufonidae that are currently controversial. To date, no mitogenomes have been reported from the genera Anaxyrus and Bufotes. Anaxyrus americanus can live in low temperature environments and Bufotes pewzowi can live in high salinity environments. We sequenced the mitogenomes of these two species to discuss the phylogenetic relationships within Bufonidae and the selection pressures experienced by specimens living in low temperature or saline environments. Like other toads, the circular mitogenomes of both species contained the typical 37 genes. Anaxyrus americanus had the highest A+T content of the complete mitogenome among the Bufonidae. In addition, A. americanus showed a negative AT-skew in the control region, whereas Bufotes pewzowi showed a positive AT-skew. Additionally, both toad species had unique molecular features in common: an ND1 gene that uses TTG as the start codon, an extra unpaired adenine (A) in the anticodon arm of trnS (AGY), and the loss of the DHU loop in trnC. The monophyly of Bufonidae was corroborated by both BI and ML trees. An analysis of selective pressure based on the 13 protein coding genes was conducted using the EasyCodeML program. In the branch model analysis, we found two branches of A. americanus and Bufotes pewzowi that were under negative selection. Additionally, we found two positively selected sites (at positions 115 and 119, BEB value > 0.90) in the ND6 protein in the site model analysis. The residue D (119) was located only in A. americanus and may be related to adaptive evolution in low temperature environments. However, there was no evidence of a positively selected site in Bufotes pewzowi in this study.

Characterization of complete mitochondrial genome of Pogonophryne albipinna (Perciformes: Artedidraconidae)

The complete mitochondrial genome of Pogonophryne albipinna was determined by the MiSeq platform, which was the first report in the family Artedidraconidae. The circular form of its mitochondrial genome was 17,086 bp, which contained the canonical eukaryotic 37 genes. The gene orders of P. albipinna was identical to the other icefish species, in which there was additional non-coding region and translocation of ND6 gene. Except for ATP6 gene and COI (GTG), 11 genes begin with the typical start codon, while incomplete stop codons (T– –) were identified in COII, ND4, and CytB. Phylogenetic tree with the currently known mitogenomes in suborder Notothenioidei showed that P. albipinna was located distinctly from those in Bathydraconidae and Nototheniidae forming a unique cluster as Artedidraconidae. The first complete mitochondrial genome of P. albipinna would be the fundamental data to understand the evolutional relationship of icefish species in the Antarctic Oceans.

Mechanisms of gene rearrangement in 13 bothids based on comparison with a newly completed mitogenome of the threespot flounder, Grammatobothus polyophthalmus (Pleuronectiformes: Bothidae)

Background

The mitogenomes of 12 teleost fish of the bothid family (order Pleuronectiformes) indicated that the genomic-scale rearrangements characterized in previous work. A novel mechanism of genomic rearrangement called the Dimer-Mitogenome and Non-Random Loss (DMNL) model was used to account for the rearrangement found in one of these bothids, Crossorhombus azureus.

Results

The 18,170 bp mitogenome of G. polyophthalmus contains 37 genes, two control regions (CRs), and the origin of replication of the L-strand (O_L). This mitogenome is characterized by genomic-scale rearrangements: genes located on the L-strand are grouped in an 8-gene cluster (Q-A-C-Y-S₁-ND6-E-P) that does not include tRNA-N; genes found on the H-strand are grouped together (F-12S … CytB-T) except for tRNA-D that was translocated inside the 8-gene L-strand cluster. Compared to non-rearranged mitogenomes of teleost fishes, gene organization in the mitogenome of G. polyophthalmus and in that of the other 12 bothids characterized thus far is very similar. These rearrangements could be sorted into four types (Type I, II, III and IV), differing in the particular combination of the CR, tRNA-D gene and 8-gene cluster and the shuffling of tRNA-V. The DMNL model was used to account for all but one gene rearrangement found in all 13 bothid mitogenomes. Translocation of tRNA-D most likely occurred after the DMNL process in 10 bothid mitogenomes and could have occurred either before or after DMNL in the three other species. During the DMNL process, the tRNA-N gene was retained rather than the expected tRNA-N′ gene. tRNA-N appears to assist in or act as O_L function when the O_L secondary structure could not be formed from intergenic sequences. A striking finding was that each of the non-transcribed genes has degenerated to a shorter intergenic spacer during the DMNL process. These findings highlight a rare phenomenon in teleost fish.

Conclusions

This result provides significant evidence to support the existence of dynamic dimeric mitogenomes and the DMNL model as the mechanism of gene rearrangement in bothid mitogenomes, which not only promotes the understanding of mitogenome structural diversity, but also sheds light on mechanisms of mitochondrial genome rearrangement and replication.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.

Comparative mitogenomics of the Decapoda reveals evolutionary heterogeneity in architecture and composition

The emergence of cost-effective and rapid sequencing approaches has resulted in an exponential rise in the number of mitogenomes on public databases in recent years, providing greater opportunity for undertaking large-scale comparative genomic and systematic research. Nonetheless, current datasets predominately come from small and disconnected studies on a limited number of related species, introducing sampling biases and impeding research of broad taxonomic relevance. This study contributes 21 crustacean mitogenomes from several under-represented decapod infraorders including Polychelida and Stenopodidea, which are used in combination with 225 mitogenomes available on NCBI to investigate decapod mitogenome diversity and phylogeny. An overview of mitochondrial gene orders (MGOs) reveals a high level of genomic variability within the Decapoda, with a large number of MGOs deviating from the ancestral arthropod ground pattern and unevenly distributed among infraorders. Despite the substantial morphological and ecological variation among decapods, there was limited evidence for correlations between gene rearrangement events and species ecology or lineage specific nucleotide substitution rates. Within a phylogenetic context, predicted scenarios of rearrangements show some MGOs to be informative synapomorphies for some taxonomic groups providing strong independent support for phylogenetic relationships. Additional comparisons for a range of mitogenomic features including nucleotide composition, strand asymmetry, unassigned regions and codon usage indicate several clade-specific trends that are of evolutionary and ecological interest.

The complete mitogenome of the splendid japalure Japalura splendida (Squamata, Agamidae)

The complete mitogenome of Japalura splendida (16,673 bp in length) is determined and analyzed in this study. It contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and one non-coding regions. All the genes in J. splendida are distributed on the H-strand, except for the ND6 gene and seven tRNA genes which are encoded on the L-strand. The phylogenetic tree suggests that J. splendida and Japalura flaviceps formed a sister group and reveals the order ((((Acanthosaura lepidogaster, Acanthosaura armata), ((J. splendida, J. flaviceps), Pseudocalotes microlepis))), Calotes versicolor) with substantial support for the monophyly.

Large-scale mitochondrial gene rearrangements in the hermit crab Pagurus nigrofascia and phylogenetic analysis of the Anomura

Complete mitochondrial genome (mitogenome) provides important information for better understanding of gene rearrangement, molecular evolution and phylogenetic analysis. Currently, only a few Paguridae mitogenomes have been reported. Herein, we described the complete mitogenome of hermit crab Pagurus nigrofascia. The total length was 15,423 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNA genes, as well as an AT-rich region. The genome composition was highly A + T biased (71.4%), and exhibited a negative AT-skew (-0.006) and GC-skew (-0.138). Eight tRNA genes, two PCGs and an AT-rich region found to be rearranged with respect to the pancrustacean ground pattern gene order. Duplication-random loss and recombination model were adopted to explain the large-scale gene rearrangement events. Two phylogenetic trees of Anomura involving 12 families were constructed. The results showed that all Paguridae species were clustered into one clade except Pagurus longicarpus, which for the first time imposed raises doubt about the morphological taxonomy of this species. Furthermore, the present study found that higher- level phylogenetic relationships within Anomura were controversial, compared with the previous studies. Our results help to better understand gene rearrangements and the evolutionary status of P. nigrofascia and lay foundation for further phylogenetic study of Anomura.

New Insight into Parrots' Mitogenomes Indicates That Their Ancestor Contained a Duplicated Region

Mitochondrial genomes of vertebrates are generally thought to evolve under strong selection for size reduction and gene order conservation. Therefore, a growing number of mitogenomes with duplicated regions changes our view on the genome evolution. Among Aves, order Psittaciformes (parrots) is especially noteworthy because of its large morphological, ecological, and taxonomical diversity, which offers an opportunity to study genome evolution in various aspects. Former analyses showed that tandem duplications comprising the control region with adjacent genes are restricted to several lineages in which the duplication occurred independently. However, using an appropriate polymerase chain reaction strategy, we demonstrate that early diverged parrot groups contain mitogenomes with the duplicated region. These findings together with mapping duplication data from other mitogenomes onto parrot phylogeny indicate that the duplication was an ancestral state for Psittaciformes. The state was inherited by main parrot groups and was lost several times in some lineages. The duplicated regions were subjected to concerted evolution with a frequency higher than the rate of speciation. The duplicated control regions may provide a selective advantage due to a more efficient initiation of replication or transcription and a larger number of replicating genomes per organelle, which may lead to a more effective energy production by mitochondria. The mitogenomic duplications were associated with phenotypic features and parrots with the duplicated region can live longer, show larger body mass as well as predispositions to a more active flight. The results have wider implications on the presence of duplications and their evolution in mitogenomes of other avian groups.

Raising orphans from a metadata morass: A researcher's guide to re-use of public 'omics data

More than 15 petabases of raw RNAseq data is now accessible through public repositories. Acquisition of other 'omics data types is expanding, though most lack a centralized archival repository. Data-reuse provides tremendous opportunity to extract new knowledge from existing experiments, and offers a unique opportunity for robust, multi-'omics analyses by merging metadata (information about experimental design, biological samples, protocols) and data from multiple experiments. We illustrate how predictive research can be accelerated by meta-analysis with a study of orphan (species-specific) genes. Computational predictions are critical to infer orphan function because their coding sequences provide very few clues. The metadata in public databases is often confusing; a test case with Zea mays mRNA seq data reveals a high proportion of missing, misleading or incomplete metadata. This metadata morass significantly diminishes the insight that can be extracted from these data. We provide tips for data submitters and users, including specific recommendations to improve metadata quality by more use of controlled vocabulary and by metadata reviews. Finally, we advocate for a unified, straightforward metadata submission and retrieval system.

Landscape genomics: natural selection drives the evolution of mitogenome in penguins

Background

Mitochondria play a key role in the balance of energy and heat production, and therefore the mitochondrial genome is under natural selection by environmental temperature and food availability, since starvation can generate more efficient coupling of energy production. However, selection over mitochondrial DNA (mtDNA) genes has usually been evaluated at the population level. We sequenced by NGS 12 mitogenomes and with four published genomes, assessed genetic variation in ten penguin species distributed from the equator to Antarctica. Signatures of selection of 13 mitochondrial protein-coding genes were evaluated by comparing among species within and among genera (Spheniscus, Pygoscelis, Eudyptula, Eudyptes and Aptenodytes). The genetic data were correlated with environmental data obtained through remote sensing (sea surface temperature [SST], chlorophyll levels [Chl] and a combination of SST and Chl [COM]) through the distribution of these species.

Results

We identified the complete mtDNA genomes of several penguin species, including ND6 and 8 tRNAs on the light strand and 12 protein coding genes, 14 tRNAs and two rRNAs positioned on the heavy strand. The highest diversity was found in NADH dehydrogenase genes and the lowest in COX genes. The lowest evolutionary divergence among species was between Humboldt (Spheniscus humboldti) and Galapagos (S. mendiculus) penguins (0.004), while the highest was observed between little penguin (Eudyptula minor) and Adélie penguin (Pygoscelis adeliae) (0.097). We identified a signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, which is consistent with the hypothesis that purifying selection is constraining mitogenome evolution to maintain Oxidative phosphorylation (OXPHOS) proteins and functionality. Pairwise species maximum-likelihood analyses of selection at codon sites suggest positive selection has occurred on ATP8 (Fixed-Effects Likelihood, FEL) and ND4 (Single Likelihood Ancestral Counting, SLAC) in all penguins. In contrast, COX1 had a signature of strong negative selection. ND4 Ka/Ks ratios were highly correlated with SST (Mantel, p-value: 0.0001; GLM, p-value: 0.00001) and thus may be related to climate adaptation throughout penguin speciation.

Conclusions

These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Such knowledge may be particularly useful for developing predictive models of how these species may respond to severe climatic changes in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4424-9) contains supplementary material, which is available to authorized users.

Isolation and characterization of the mitochondrial genome of Gymnodraco acuticeps (Perciformes: Bathydraconidae) with phylogenetic consideration

Gymnodraco acuticepsis is an Antarctic fish living in the Southern Ocean. Until now, studies on G. acuticeps are still limited. As an Antarctic fish, obtaining and characterization of the mitochondrial genome of G. acuticeps will be important for elucidation of the mechanism of cold-adapting evolution in mitochondrion. In this study, we first isolated and characterized the mitochondrial genome sequence of G. acuticeps with 15,987 bp in length. It contained of 34 genes (12 protein-coding genes, 20 transfer RNA genes, 2 ribosomal RNA genes) and a partial putative control region. Gene organization and nucleotide composition of obtained mito-genome were similar to those of other Antarctic fish. Twenty-eight genes were encoded by heavy strand, while six genes were encoded by light strand. Further, the phylogenetic tree, which based on 12 protein-coding genes, revealed that the G. acuticeps was genetically closest to species Parachaenichthys charcoti among 18 species. We hope this work would be helpful for the population genetics and molecular evolution studies.

Structure and variation of the mitochondrial genome of fishes

Background

The mitochondrial (mt) genome has been used as an effective tool for phylogenetic and population genetic analyses in vertebrates. However, the structure and variability of the vertebrate mt genome are not well understood. A potential strategy for improving our understanding is to conduct a comprehensive comparative study of large mt genome data. The aim of this study was to characterize the structure and variability of the fish mt genome through comparative analysis of large datasets.

Results

An analysis of the secondary structure of proteins for 250 fish species (248 ray-finned and 2 cartilaginous fishes) illustrated that cytochrome c oxidase subunits (COI, COII, and COIII) and a cytochrome bc1 complex subunit (Cyt b) had substantial amino acid conservation. Among the four proteins, COI was the most conserved, as more than half of all amino acid sites were invariable among the 250 species. Our models identified 43 and 58 stems within 12S rRNA and 16S rRNA, respectively, with larger numbers than proposed previously for vertebrates. The models also identified 149 and 319 invariable sites in 12S rRNA and 16S rRNA, respectively, in all fishes. In particular, the present result verified that a region corresponding to the peptidyl transferase center in prokaryotic 23S rRNA, which is homologous to mt 16S rRNA, is also conserved in fish mt 16S rRNA. Concerning the gene order, we found 35 variations (in 32 families) that deviated from the common gene order in vertebrates. These gene rearrangements were mostly observed in the area spanning the ND5 gene to the control region as well as two tRNA gene cluster regions (IQM and WANCY regions). Although many of such gene rearrangements were unique to a specific taxon, some were shared polyphyletically between distantly related species.

Conclusions

Through a large-scale comparative analysis of 250 fish species mt genomes, we elucidated various structural aspects of the fish mt genome and the encoded genes. The present results will be important for understanding functions of the mt genome and developing programs for nucleotide sequence analysis. This study demonstrated the significance of extensive comparisons for understanding the structure of the mt genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3054-y) contains supplementary material, which is available to authorized users.

Integrated studies of organismal plasticity through physiological and transcriptomic approaches: examples from marine polar regions

Transcriptomic methods are now widely used in functional genomic research. The vast amount of information received from these studies comes along with the challenge of developing a precise picture of the functional consequences and the characteristic regulatory mechanisms. Here we assess recent studies in marine species and their adaptation to polar (and seasonal) cold and explore how they have been able to draw reliable conclusions from transcriptomic patterns on functional consequences in the organisms. Our analysis indicates that the interpretation of transcriptomic data suffers from insufficient understanding of the consequences for whole organism performance and fitness and comes with the risk of supporting only preliminary and superficial statements.We propose that the functional understanding of transcriptomic data may be improved by their tighter integration into overarching physiological concepts that support the more specific interpretation of the 'omics' data and, at the same time, can be developed further through embedding the transcriptomic phenomena observed. Such possibilities have not been fully exploited.In the context of thermal adaptation and limitation, we explore preliminary evidence that the concept of oxygen and capacity limited thermal tolerance (OCLTT) may provide sufficient complexity to guide the integration of such data and the development of associated functional hypotheses. At the same time, we identify a lack of methodological approaches linking genes and function to higher levels of integration, in terms of organism and ecosystem functioning, at temporal and geographical scales, to support more reliable conclusions and be predictive with respect to the effects of global changes.

Antarctic notothenioid fish: what are the future consequences of 'losses' and 'gains' acquired during long-term evolution at cold and stable temperatures?

Antarctic notothenioids dominate the fish fauna of the Southern Ocean. Evolution for millions of years at cold and stable temperatures has led to the acquisition of numerous biochemical traits that allow these fishes to thrive in sub-zero waters. The gain of antifreeze glycoproteins has afforded notothenioids the ability to avert freezing and survive at temperatures often hovering near the freezing point of seawater. Additionally, possession of cold-adapted proteins and membranes permits them to sustain appropriate metabolic rates at exceptionally low body temperatures. The notothenioid genome is also distinguished by the disappearance of traits in some species, losses that might prove costly in a warmer environment. Perhaps the best-illustrated example is the lack of expression of hemoglobin in white-blooded icefishes from the family Channichthyidae. Loss of key elements of the cellular stress response, notably the heat shock response, has also been observed. Along with their attainment of cold tolerance, notothenioids have developed an extreme stenothermy and many species perish at temperatures only a few degrees above their habitat temperatures. Thus, in light of today's rapidly changing climate, it is critical to evaluate how these extreme stenotherms will respond to rising ocean temperatures. It is conceivable that the remarkable cold specialization of notothenioids may ultimately leave them vulnerable to future thermal increases and threaten their fitness and survival. Within this context, our review provides a current summary of the biochemical losses and gains that are known for notothenioids and examines these cold-adapted traits with a focus on processes underlying thermal tolerance and acclimation capacity.

The complete mitochondrial genome sequence and gene organization of Trematomus bernacchii (Perciformes: Nototheniidae) with phylogenetic consideration

In this study, the complete mitochondrial DNA sequence was obtained from Trematomus bernacchii. The full length of mitogenome was 16 018 bp, including 20 transfer RNA genes, a 16S ribosomal RNA genes and an incomplete 12S rRNA, 12 protein-coding genes and an incomplete ND6 protein-coding gene, and 2 non-coding regions: control region (D-loop) and origin of light-strand (L-strand) replication (OL). Most genes were located on the heavy-strand (H-strand), while ND6 and seven tRNA genes were located on the L-strand. T. bernacchii’s mitogenome was comprised of 24.52% for A, 25.44% for C, 20.38% for G and 29.66% for T, with a higher A + T content (54.18%). From the NJ phylogenetic tree, it revealed that T. bernacchii’s was genetically closest to species Dissostichus eleginoides among seven species within suborder Notothenioidei. This study could lay the foundation for the future studies in species identification, taxonomic status, molecular systematics, stock evaluation and conservation genetics for T. bernacchii’s.

Complete mitochondrial genome of the Antarctic icefish, Chaenocephalus aceratus (Perciforms, Channichthyidae)

The complete sequence of the mitochondrial genome of the Antarctic icefish Chaenocephalus aceratus was determined to be 17311 bp in length, and to contain 13 protein coding genes (PCGs), 22 tRNA genes and 2 rRNA genes. The total A+T content is 52.8%. The notothenioid-exclusive ND6/tRNAGlu translocation was observed in the mitogenome of C. aceratus. Generally, the order and contents of the other genes are identical with those of other fishes. Antarctic icefishes, the only vertebrates which do not have hemoglobins, have evolved to survive sub-freezing temperature. Therefore, the whole mitogenome sequences of C. aceratus will provide the insights into resolving the evolutionary history of icefish.

Complete mitochondrial genome of the Antarctic silverfish, Pleuragramma antarcticum (Perciformes, Nototheniidae)

The complete sequence of the mitochondrial genome of Pleuragramma antarcticum was obtained from the assembled genome data. The complete sequence was determined to be 18,460 bp in length, containing 13 protein coding genes, 22 tRNA genes, 2 rRNA genes, 2 control regions and a repeat. The general order and contents of the genes are identical with those of other fishes. The Antarctic notothenioid-specific translocation of ND6/tRNA(Glu) exists in the mitogenome. Pleuragramma antarcticum is the only pelagic species which dominates Antarctic fish fauna by more than 90%. Thus, the whole mitogenome sequences of P. antarcticum will provide the basis to understand the deep phylogeny of the Antarctic species.

Recombination and evolution of duplicate control regions in the mitochondrial genome of the Asian big-headed turtle, Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3' end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.

A comprehensive description and evolutionary analysis of 22 grouper (perciformes, epinephelidae) mitochondrial genomes with emphasis on two novel genome organizations

Groupers of the family Epinephelidae are a diverse and economically valuable group of reef fishes. To investigate the evolution of their mitochondrial genomes we characterized and compared these genomes among 22 species, 17 newly sequenced. Among these fishes we identified three distinct genome organizations, two of them never previously reported in vertebrates. In 19 of these species, mitochondrial genomes followed the typical vertebrate canonical organization with 13 protein-coding genes, 22 tRNAs, two rRNAs, and a non-coding control region. Differing from this, members of genus Variola have an extra tRNA-Ile between tRNA-Val and 16S rRNA. Evidence suggests that this evolved from tRNA-Val via a duplication event due to slipped strand mispairing during replication. Additionally, Cephalopholisargus has an extra tRNA-Asp in the midst of the control region, likely resulting from long-range duplication of the canonical tRNA-Asp through illicit priming of mitochondrial replication by tRNAs. Along with their gene contents, we characterized the regulatory elements of these mitochondrial genomes' control regions, including putative termination-associated sequences and conserved sequence blocks. Looking at the mitochondrial genomic constituents, rRNA and tRNA are the most conserved, followed by protein-coding genes, and non-coding regions are the most divergent. Divergence rates vary among the protein-coding genes, and the three cytochrome oxidase subunits (COI, II, III) are the most conserved, while NADH dehydrogenase subunit 6 (ND6) and the ATP synthase subunit 8 (ATP8) are the most divergent. We then tested the phylogenetic utility of this new mt genome data using 12 protein-coding genes of 48 species from the suborder Percoidei. From this, we provide further support for the elevation of the subfamily Epinephelinae to family Epinephelidae, the resurrection of the genus Hyporthodus, and the combination of the monotypic genera Anyperodon and Cromileptes to genus Epinephelus, and Aethaloperca to genus Cephalopholis.

Mitochondrial acclimation capacities to ocean warming and acidification are limited in the antarctic Nototheniid Fish, Notothenia rossii and Lepidonotothen squamifrons

Antarctic notothenioid fish are characterized by their evolutionary adaptation to the cold, thermostable Southern Ocean, which is associated with unique physiological adaptations to withstand the cold and reduce energetic requirements but also entails limited compensation capacities to environmental change. This study compares the capacities of mitochondrial acclimation to ocean warming and acidification between the Antarctic nototheniid Notothenia rossii and the sub-Antarctic Lepidonotothen squamifrons, which share a similar ecology, but different habitat temperatures. After acclimation of L. squamifrons to 9°C and N. rossii to 7°C (normocapnic/hypercapnic, 0.2 kPa CO2/2000 ppm CO2) for 4-6 weeks, we compared the capacities of their mitochondrial respiratory complexes I (CI) and II (CII), their P/O ratios (phosphorylation efficiency), proton leak capacities and mitochondrial membrane fatty acid compositions. Our results reveal reduced CII respiration rates in warm-acclimated L. squamifrons and cold hypercapnia-acclimated N. rossii. Generally, L. squamifrons displayed a greater ability to increase CI contribution during acute warming and after warm-acclimation than N. rossii. Membrane unsaturation was not altered by warm or hypercapnia-acclimation in both species, but membrane fatty acids of warm-acclimated L. squamifrons were less saturated than in warm normocapnia-/hypercapnia-acclimated N. rossii. Proton leak capacities were not affected by warm or hypercapnia-acclimation of N. rossii. We conclude that an acclimatory response of mitochondrial capacities may include higher thermal plasticity of CI supported by enhanced utilization of anaplerotic substrates (via oxidative decarboxylation reactions) feeding into the citrate cycle. L. squamifrons possesses higher relative CI plasticities than N. rossii, which may facilitate the usage of energy efficient NADH-related substrates under conditions of elevated energy demand, possibly induced by ocean warming and acidification. The observed adjustments of electron transport system complexes with a higher flux through CI under warming and acidification suggest a metabolic acclimation potential of the sub-Antarctic L. squamifrons, but only limited acclimation capacities for N. rossii.

Cloning and characterization of a Δ9-desaturase gene of the Antarctic fish Chionodraco hamatus and Trematomus bernacchii

Chionodraco hamatus and Trematomus bernacchii are perciforms, members of the fish suborder Notothenioidei that live in the Antarctic Ocean and experience very cold and persistent environmental temperature. These fish have biochemical and molecular features that allow them to live at these extreme cold temperatures. Fine tuning of the level of unsaturated fatty acids content in membrane is a key mechanism of living organisms to adapt to cold and high temperatures. Desaturases are key enzymes that synthesize unsaturated fatty acyl-CoAs from saturated fatty acids. We cloned and sequenced a Δ(9)-desaturase gene and its cDNA of C. hamatus, and the cDNA of T. bernacchii. The coded proteins are virtually identical and share homology to other Δ(9)-desaturase fish sequences. These proteins contain, in the first trans-membrane domain, two cysteine residues that may form a disulfur bond present in the corresponding membrane region of Δ(9)-desaturase proteins of other Antarctic fish but not in Eleginops maclovinus that experiences higher environmental temperatures and in all other Δ(9)-desaturase genes of mammals present in data bases. C. hamatus Δ(9)-desaturase gene complements a Saccharomyces cerevisiae mutant lacking Δ(9)-desaturase (Ole1) gene. Analysis of sequence homology of the trans-membrane domains of Δ(9)-desaturase and the cytoplasmic region of the same proteins of Antarctic fish, non-Antarctic fish and mammals suggest that the significant differences found in the homologous sequences of the first trans-membrane domain may be due to the specific lipid content of their membrane.

Replication-Associated Mutational Pressure (RMP) Governs Strand-Biased Compositional Asymmetry (SCA) and Gene Organization in Animal Mitochondrial Genomes

The nucleotide composition of the light (L-) and heavy (H-) strands of animal mitochondrial genomes is known to exhibit strand-biased compositional asymmetry (SCA). One of the possibilities is the existence of a replication-associated mutational pressure (RMP) that may introduce characteristic nucleotide changes among mitochondrial genomes of different animal lineages. Here, we discuss the influence of RMP on nucleotide and amino acid compositions as well as gene organization. Among animal mitochondrial genomes, RMP may represent the major force that compels the evolution of mitochondrial protein-coding genes, coupled with other process-based selective pressures, such as on components of translation machinery— tRNAs and their anticodons. Through comparative analyses of sequenced mitochondrial genomes among diverse animal lineages and literature reviews, we suggest a strong RMP effect, observed among invertebrate mitochondrial genes as compared to those of vertebrates, that is either a result of positive selection on the invertebrate or a relaxed selective pressure on the vertebrate mitochondrial genes.

Multiple independent origins of mitochondrial control region duplications in the order Psittaciformes

Mitochondrial genomes are generally thought to be under selection for compactness, due to their small size, consistent gene content, and a lack of introns or intergenic spacers. As more animal mitochondrial genomes are fully sequenced, rearrangements and partial duplications are being identified with increasing frequency, particularly in birds (Class Aves). In this study, we investigate the evolutionary history of mitochondrial control region states within the avian order Psittaciformes (parrots and cockatoos). To this aim, we reconstructed a comprehensive multi-locus phylogeny of parrots, used PCR of three diagnostic fragments to classify the mitochondrial control region state as single or duplicated, and mapped these states onto the phylogeny. We further sequenced 44 selected species to validate these inferences of control region state. Ancestral state reconstruction using a range of weighting schemes identified six independent origins of mitochondrial control region duplications within Psittaciformes. Analysis of sequence data showed that varying levels of mitochondrial gene and tRNA homology and degradation were present within a given clade exhibiting duplications. Levels of divergence between control regions within an individual varied from 0-10.9% with the differences occurring mainly between 51 and 225 nucleotides 3' of the goose hairpin in domain I. Further investigations into the fates of duplicated mitochondrial genes, the potential costs and benefits of having a second control region, and the complex relationship between evolutionary rates, selection, and time since duplication are needed to fully explain these patterns in the mitochondrial genome.

Mitochondrial function in Antarctic nototheniids with ND6 translocation

Fish of the suborder Notothenioidei have successfully radiated into the Southern Ocean and today comprise the dominant fish sub-order in Antarctic waters in terms of biomass and species abundance. During evolution in the cold and stable Antarctic climate, the Antarctic lineage of notothenioids developed several unique physiological adaptations, which make them extremely vulnerable to the rapid warming of Antarctic waters currently observed. Only recently, a further phenomenon exclusive to notothenioid fish was reported: the translocation of the mitochondrial gene encoding the NADH Dehydrogenase subunit 6 (ND6), an indispensable part of complex I in the mitochondrial electron transport system.This study investigated the potential physiological consequences of ND6 translocation for the function and thermal sensitivity of the electron transport system in isolated liver mitochondria of the two nototheniid species Notothenia coriiceps and Notothenia rossii, with special attention to the contributions of complex I (NADH DH) and complex II (Succinate DH) to oxidative phosphorylation. Furthermore, enzymatic activities of NADH:Cytochrome c Oxidoreductase and Cytochrome C Oxidase were measured in membrane-enriched tissue extracts.During acute thermal challenge (0-15°C), capacities of mitochondrial respiration and enzymatic function in the liver could only be increased until 9°C. Mitochondrial complex I (NADH Dehydrogenase) was fully functional but displayed a higher thermal sensitivity than the other complexes of the electron transport system, which may specifically result from its unique amino acid composition, revealing a lower degree of stability in notothenioids in general. We interpret the translocation of ND6 as functionally neutral but the change in amino acid sequence as adaptive and supportive of cold stenothermy in Antarctic nototheniids. From these findings, an enhanced sensitivity to ocean warming can be deduced for Antarctic notothenioid fish.

A phylogeny of the temperate seabasses (Moronidae) characterized by a translocation of the mt-nd6 gene

The entire mitochondrial genome of the striped bass Morone saxatilis was sequenced together with the mitochondrial (mt) control regions of the white bass Morone chrysops, white perch Morone americana, yellow bass Morone mississippiensis, spotted seabass Dicentrarchus punctatus, European seabass Dicentrarchus labrax and the Japanese seabass Lateolabrax japonicus. The resultant 17 580 base pair circular genome of M. saxatilis contains 38 genes (13 proteins, 23 transfer RNAs and two ribosomal RNAs) and a control region bordered by the proline and phenylalanine mitochondrial tRNAs. Gene arrangement was similar to other vertebrates, except that the mt-nd6 gene was found within the control region rather than the canonical position between the mt-nd5 and mt-cyb genes. This translocation was found in all the Morone and Dicentrarchus species studied without functional copies or pseudogenes in the ancestral position. In L. japonicus, the mt-nd6 gene was found in the canonical position without evidence of an mt-nd6 gene in the control region. A Bayesian analysis of these and published mt-nd6 sequences from 45 other Perciformes grouped the Morone and Dicentrarchus species monophyletically with a probability of 1·00 with respect to L. japonicus and all other perciforms, and placed the Dicentrarchus species in the basal position. These data reinforce current placement of L. japonicus outside the Moronidae and provide a clear evolutionary character to define this family. The phylogeny of the Moronidae presented here also supports the hypothesis of an anadromous common ancestor to this family that gave rise to the North American estuarine and freshwater species. A series of tandem repeats previously reported in M. saxatilis was found in the control region of all Morone species between the mt-nd6 and mt-rnr1 genes, but not in either Dicentrarchus species, which reinforces the continued use of these two separate genera.