Molecular evolution of cytochrome C oxidase underlies high-altitude adaptation in the bar-headed goose

Mass Transport: Circulatory System with Emphasis on Nonendothermic Species

Mass transport can be generally defined as movement of material matter. The circulatory system then is a biological example given its role in the movement in transporting gases, nutrients, wastes, and chemical signals. Comparative physiology has a long history of providing new insights and advancing our understanding of circulatory mass transport across a wide array of circulatory systems. Here we focus on circulatory function of nonmodel species. Invertebrates possess diverse convection systems; that at the most complex generate pressures and perform at a level comparable to vertebrates. Many invertebrates actively modulate cardiovascular function using neuronal, neurohormonal, and skeletal muscle activity. In vertebrates, our understanding of cardiac morphology, cardiomyocyte function, and contractile protein regulation by Ca2+ highlights a high degree of conservation, but differences between species exist and are coupled to variable environments and body temperatures. Key regulators of vertebrate cardiac function and systemic blood pressure include the autonomic nervous system, hormones, and ventricular filling. Further chemical factors regulating cardiovascular function include adenosine, natriuretic peptides, arginine vasotocin, endothelin 1, bradykinin, histamine, nitric oxide, and hydrogen sulfide, to name but a few. Diverse vascular morphologies and the regulation of blood flow in the coronary and cerebral circulations are also apparent in nonmammalian species. Dynamic adjustments of cardiovascular function are associated with exercise on land, flying at high altitude, prolonged dives by marine mammals, and unique morphology, such as the giraffe. Future studies should address limits of gas exchange and convective transport, the evolution of high arterial pressure across diverse taxa, and the importance of the cardiovascular system adaptations to extreme environments. © 2017 American Physiological Society. Compr Physiol 7:17-66, 2017.

Rates and patterns of molecular evolution in freshwater versus terrestrial insects

Insect lineages have crossed between terrestrial and aquatic habitats many times, for both immature and adult life stages. We explore patterns in molecular evolutionary rates between 42 sister pairs of related terrestrial and freshwater insect clades using publicly available protein-coding DNA sequence data from the orders Coleoptera, Diptera, Lepidoptera, Hemiptera, Mecoptera, Trichoptera, and Neuroptera. We furthermore test for habitat-associated convergent molecular evolution in the cytochrome c oxidase subunit I (COI) gene in general and at a particular amino acid site previously reported to exhibit habitat-linked convergence within an aquatic beetle group. While ratios of nonsynonymous-to-synonymous substitutions across available loci were higher in terrestrial than freshwater-associated taxa in 26 of 42 lineage pairs, a stronger trend was observed (20 of 31, p_binomial = 0.15, p_Wilcoxon = 0.017) when examining only terrestrial-aquatic pairs including fully aquatic taxa. We did not observe any widespread changes at particular amino acid sites in COI associated with habitat shifts, although there may be general differences in selection regime linked to habitat.

Genetic signals of high-altitude adaptation in amphibians: a comparative transcriptome analysis

BACKGROUND

High-altitude adaptation provides an excellent system for studying how organisms cope with multiple environmental stressors and interacting genetic modifications. To explore the genetic basis of high-altitude adaptation in poikilothermic animals, we acquired transcriptome sequences from a high-altitude population and a low-altitude population of the Asiatic toad (Bufo gargarizans). Transcriptome data from another high-altitude amphibian, Rana kukunoris and its low-altitude relative R. chensiensis, which are from a previous study, were also incorporated into our comparative analysis.

RESULTS

More than 40,000 transcripts were obtained from each transcriptome, and 5107 one-to-one orthologs were identified among the four taxa for comparative analysis. A total of 29 (Bufo) and 33 (Rana) putative positively selected genes were identified for the two high-altitude species, which were mainly concentrated in nutrient metabolism related functions. Using SNP-tagging and F_ST outlier analysis, we further tested 89 other nutrient metabolism related genes for signatures of natural selection, and found that two genes, CAPN2 and ITPR1, were likely under balancing selection. We did not detect any positively selected genes associated with response to hypoxia.

CONCLUSIONS

Amphibians clearly employ different genetic mechanisms for high-altitude adaptation compared to endotherms. Modifications of genes associated with nutrient metabolism feature prominently while genes related to hypoxia tolerance appear to be insignificant. Poikilotherms represent the majority of animal diversity, and we hope that our results will provide useful directions for future studies of amphibians as well as other poikilotherms.

Mitochondrial physiology in the skeletal and cardiac muscles is altered in torrent ducks, Merganetta armata, from high altitudes in the Andes

Torrent ducks inhabit fast-flowing rivers in the Andes from sea level to altitudes up to 4500 m. We examined the mitochondrial physiology that facilitates performance over this altitudinal cline by comparing the respiratory capacities of permeabilized fibers, the activities of 16 key metabolic enzymes and the myoglobin content in muscles between high- and low-altitude populations of this species. Mitochondrial respiratory capacities (assessed using substrates of mitochondrial complexes I, II and/or IV) were higher in highland ducks in the gastrocnemius muscle - the primary muscle used to support swimming and diving - but were similar between populations in the pectoralis muscle and the left ventricle. The heightened respiratory capacity in the gastrocnemius of highland ducks was associated with elevated activities of cytochrome oxidase, phosphofructokinase, pyruvate kinase and malate dehydrogenase (MDH). Although respiratory capacities were similar between populations in the other muscles, highland ducks had elevated activities of ATP synthase, lactate dehydrogenase, MDH, hydroxyacyl CoA dehydrogenase and creatine kinase in the left ventricle, and elevated MDH activity and myoglobin content in the pectoralis. Thus, although there was a significant increase in the oxidative capacity of the gastrocnemius in highland ducks, which correlates with improved performance at high altitudes, the variation in metabolic enzyme activities in other muscles not correlated to respiratory capacity, such as the consistent upregulation of MDH activity, may serve other functions that contribute to success at high altitudes.

Mitonuclear coevolution as the genesis of speciation and the mitochondrial DNA barcode gap

Mitochondrial genes are widely used in taxonomy and systematics because high mutation rates lead to rapid sequence divergence and because such changes have long been assumed to be neutral with respect to function. In particular, the nucleotide sequence of the mitochondrial gene cytochrome c oxidase subunit 1 has been established as a highly effective DNA barcode for diagnosing the species boundaries of animals. Rarely considered in discussions of mitochondrial evolution in the context of systematics, speciation, or DNA barcodes, however, is the genomic architecture of the eukaryotes: Mitochondrial and nuclear genes must function in tight coordination to produce the complexes of the electron transport chain and enable cellular respiration. Coadaptation of these interacting gene products is essential for organism function. I extend the hypothesis that mitonuclear interactions are integral to the process of speciation. To maintain mitonuclear coadaptation, nuclear genes, which code for proteins in mitochondria that cofunction with the products of mitochondrial genes, must coevolve with rapidly changing mitochondrial genes. Mitonuclear coevolution in isolated populations leads to speciation because population-specific mitonuclear coadaptations create between-population mitonuclear incompatibilities and hence barriers to gene flow between populations. In addition, selection for adaptive divergence of products of mitochondrial genes, particularly in response to climate or altitude, can lead to rapid fixation of novel mitochondrial genotypes between populations and consequently to disruption in gene flow between populations as the initiating step in animal speciation. By this model, the defining characteristic of a metazoan species is a coadapted mitonuclear genotype that is incompatible with the coadapted mitochondrial and nuclear genotype of any other population.

The complete mitochondrial genome of Anser indicus (Aves, Anseriformes, Anatidae)

The bar-headed goose is known one of the world's highest-flying birds. In this study, the complete mitochondrial genome of Anser indicus (16,728 bp in length) was sequenced. Similar to the typical mtDNA of other vertebrates, goose mtDNA contained 37 genes (13 protein-coding genes, 2 rRNA genes and 22 tRNA genes) and a non-coding region (D-loop). The characteristics of the mitochondrial genome were analyzed in detail. We deduce that ND5 may be a major gene required for adaptation to high-altitude flight. The complete mitochondrial genome sequence of A. indicus obtained will be useful for phylogenetics, and biological characteristics.

Biogeographic history and high-elevation adaptations inferred from the mitochondrial genome of Glyptosternoid fishes (Sisoridae, Siluriformes) from the southeastern Tibetan Plateau

Background

The distribution of the Chinese Glyptosternoid catfish is limited to the rivers of the Tibetan Plateau and peripheral regions, especially the drainage areas of southeastern Tibet. Therefore, Glyptosternoid fishes are ideal for reconstructing the geological history of the southeastern Tibet drainage patterns and mitochondrial genetic adaptions to high elevations.

Results

Our phylogenetic results support the monophyly of the Sisoridae and the Glyptosternoid fishes. The reconstructed ancestral geographical distribution suggests that the ancestral Glyptosternoids was widely distributed throughout the Brahmaputra drainage in the eastern Himalayas and Tibetan area during the Late Miocene (c. 5.5 Ma). We found that the Glyptosternoid fishes lineage had a higher ratio of nonsynonymous to synonymous substitutions than those found in non-Glyptosternoids. In addition, ω_pss was estimated to be 10.73, which is significantly higher than 1 (p-value 0.0002), in COX1, which indicates positive selection in the common ancestral branch of Glyptosternoid fishes in China. We also found other signatures of positive selection in the branch of specialized species. These results imply mitochondrial genetic adaptation to high elevations in the Glyptosternoids.

Conclusions

We reconstructed a possible scenario for the southeastern Tibetan drainage patterns based on the adaptive geographical distribution of the Chinese Glyptosternoids in this drainage. The Glyptosternoids may have experienced accelerated evolutionary rates in mitochondrial genes that were driven by positive selection to better adapt to the high-elevation environment of the Tibetan Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0516-9) contains supplementary material, which is available to authorized users.

Comparative analysis of the gastrointestinal microbial communities of bar-headed goose (Anser indicus) in different breeding patterns by high-throughput sequencing

The bar-headed goose is currently one of the most popular species for rare birds breeding in China. However, bar-headed geese in captivity display a reduced reproductive rate. The gut microbiome has been shown to influence host factors such as nutrient and energy metabolism, immune homeostasis and reproduction. It is therefore of great scientific and agriculture value to analyze the microbial communities associated with bar-headed geese in order to improve their reproductive rate. Here we describe the first comparative study of the gut microbial communities of bar-headed geese in three different breeding pattern groups by 16SrRNA sequences using the Illumina MiSeq platform. The results showed that Firmicutes predominated (58.33%) among wild bar-headed geese followed by Proteobacteria (30.67%), Actinobacteria (7.33%) and Bacteroidetes (3.33%). In semi-artificial breeding group, Firmicutes was also the most abundant bacteria (62.00%), followed by Bacteroidetes (28.67%), Proteobacteria (4.20%), Actinobacteria (3.27%) and Fusobacteria (1.51%). The microbial communities of artificial breeding group were dominated by Firmicutes (60.67%), Fusobacteria (29.67%) and Proteobacteria (9.33%). Wild bar-headed geese had a significant higher relative abundance of Proteobacteria and Actinobacteria, while semi-artificial breeding bar-headed geese had significantly more Bacteroidetes. The semi-artificial breeding group had the highest microbial community diversity and richness, followed by wild group, and then the artificial breeding group. The marked differences of genus level group-specific microbes create a baseline for future bar-headed goose microbiology research.

How bar-headed geese fly over the Himalayas

Bar-headed geese cross the Himalayas on one of the most iconic high-altitude migrations in the world. Heart rates and metabolic costs of flight increase with elevation and can be near maximal during steep climbs. Their ability to sustain the high oxygen demands of flight in air that is exceedingly oxygen-thin depends on the unique cardiorespiratory physiology of birds in general along with several evolved specializations across the O2 transport cascade.

Mitochondrial divergence between slow- and fast-aging garter snakes

Mitochondrial function has long been hypothesized to be intimately involved in aging processes--either directly through declining efficiency of mitochondrial respiration and ATP production with advancing age, or indirectly, e.g., through increased mitochondrial production of damaging free radicals with age. Yet we lack a comprehensive understanding of the evolution of mitochondrial genotypes and phenotypes across diverse animal models, particularly in species that have extremely labile physiology. Here, we measure mitochondrial genome-types and transcription in ecotypes of garter snakes (Thamnophis elegans) that are adapted to disparate habitats and have diverged in aging rates and lifespans despite residing in close proximity. Using two RNA-seq datasets, we (1) reconstruct the garter snake mitochondrial genome sequence and bioinformatically identify regulatory elements, (2) test for divergence of mitochondrial gene expression between the ecotypes and in response to heat stress, and (3) test for sequence divergence in mitochondrial protein-coding regions in these slow-aging (SA) and fast-aging (FA) naturally occurring ecotypes. At the nucleotide sequence level, we confirmed two (duplicated) mitochondrial control regions one of which contains a glucocorticoid response element (GRE). Gene expression of protein-coding genes was higher in FA snakes relative to SA snakes for most genes, but was neither affected by heat stress nor an interaction between heat stress and ecotype. SA and FA ecotypes had unique mitochondrial haplotypes with amino acid substitutions in both CYTB and ND5. The CYTB amino acid change (Isoleucine → Threonine) was highly segregated between ecotypes. This divergence of mitochondrial haplotypes between SA and FA snakes contrasts with nuclear gene-flow estimates, but correlates with previously reported divergence in mitochondrial function (mitochondrial oxygen consumption, ATP production, and reactive oxygen species consequences).

Evidence of adaptive evolution of alpine pheasants to high-altitude environment from mitogenomic perspective

Adaptive evolutions to high-altitude adaptation have been intensively studied in mammals. However, considering the additional vertebrate groups, new perception regarding selection challenged by high-altitude stress on mitochondrial genome can be gained. To test this hypothesis, we compiled and analyzed the mitochondrial genomes of 5 alpine pheasants and 12 low-altitude species in Phasianidae. The results that evolutionary rates of ATP6 and ND6 showing significant fluctuation among branches when involved with five alpine pheasants revealed both genes might have implications with adapting to highland environment. The radical physico-chemical property changes identified by the modified MM01 model, including composition (C) and equilibrium constant (ionization of COOH) (Pk') in ATP6 and beta-structure tendencies (Pβ), Pk', and long-range non-bonded energy (El) in ND6, suggested that minor overall adjustments in size, protein conformation and relative orientation of reaction interfaces have been optimized to provide the ideal environments for electron transfer, proton translocation and generation of adenosine triphosphate (ATP). Additionally, three unique substitution sites were identified under selection in ND6, which could be potentially important adaptive changes contributing to cellular energy production. Our findings suggested that adaptive evolution may occur in alpine pheasants, which are an important complement to the knowledge of genetic mechanisms against the high-altitude environment in non-mammal animals.

Adaptation of the Mitochondrial Genome in Cephalopods: Enhancing Proton Translocation Channels and the Subunit Interactions

Mitochondrial protein-coding genes (mt genes) encode subunits forming complexes of crucial cellular pathways, including those involved in the vital process of oxidative phosphorylation (OXPHOS). Despite the vital role of the mitochondrial genome (mt genome) in the survival of organisms, little is known with respect to its adaptive implications within marine invertebrates. The molluscan Class Cephalopoda is represented by a marine group of species known to occupy contrasting environments ranging from the intertidal to the deep sea, having distinct metabolic requirements, varied body shapes and highly advanced visual and nervous systems that make them highly competitive and successful worldwide predators. Thus, cephalopods are valuable models for testing natural selection acting on their mitochondrial subunits (mt subunits). Here, we used concatenated mt genes from 17 fully sequenced mt genomes of diverse cephalopod species to generate a robust mitochondrial phylogeny for the Class Cephalopoda. We followed an integrative approach considering several branches of interest–covering cephalopods with distinct morphologies, metabolic rates and habitats–to identify sites under positive selection and localize them in the respective protein alignment and/or tridimensional structure of the mt subunits. Our results revealed significant adaptive variation in several mt subunits involved in the energy production pathway of cephalopods: ND5 and ND6 from Complex I, CYTB from Complex III, COX2 and COX3 from Complex IV, and in ATP8 from Complex V. Furthermore, we identified relevant sites involved in protein-interactions, lining proton translocation channels, as well as disease/deficiencies related sites in the aforementioned complexes. A particular case, revealed by this study, is the involvement of some positively selected sites, found in Octopoda lineage in lining proton translocation channels (site 74 from ND5) and in interactions between subunits (site 507 from ND5) of Complex I.

The effects of temperature on aerobic metabolism: towards a mechanistic understanding of the responses of ectotherms to a changing environment

Because of its profound effects on the rates of biological processes such as aerobic metabolism, environmental temperature plays an important role in shaping the distribution and abundance of species. As temperature increases, the rate of metabolism increases and then rapidly declines at higher temperatures – a response that can be described using a thermal performance curve (TPC). Although the shape of the TPC for aerobic metabolism is often attributed to the competing effects of thermodynamics, which can be described using the Arrhenius equation, and the effects of temperature on protein stability, this account represents an over-simplification of the factors acting even at the level of single proteins. In addition, it cannot adequately account for the effects of temperature on complex multistep processes, such as aerobic metabolism, that rely on mechanisms acting across multiple levels of biological organization. The purpose of this review is to explore our current understanding of the factors that shape the TPC for aerobic metabolism in response to acute changes in temperature, and to highlight areas where this understanding is weak or insufficient. Developing a more strongly grounded mechanistic model to account for the shape of the TPC for aerobic metabolism is crucial because these TPCs are the foundation of several recent attempts to predict the responses of species to climate change, including the metabolic theory of ecology and the hypothesis of oxygen and capacity-limited thermal tolerance.

Mitonuclear Ecology

Eukaryotes were born of a chimeric union between two prokaryotes--the progenitors of the mitochondrial and nuclear genomes. Early in eukaryote evolution, most mitochondrial genes were lost or transferred to the nucleus, but a core set of genes that code exclusively for products associated with the electron transport system remained in the mitochondrion. The products of these mitochondrial genes work in intimate association with the products of nuclear genes to enable oxidative phosphorylation and core energy production. The need for coadaptation, the challenge of cotransmission, and the possibility of genomic conflict between mitochondrial and nuclear genes have profound consequences for the ecology and evolution of eukaryotic life. An emerging interdisciplinary field that I call "mitonuclear ecology" is reassessing core concepts in evolutionary ecology including sexual reproduction, two sexes, sexual selection, adaptation, and speciation in light of the interactions of mitochondrial and nuclear genomes.

High-altitude adaptation of Tibetan chicken from MT-COI and ATP-6 perspective

The problem of hypoxia adaptation in high altitudes is an unsolved brainteaser in the field of life sciences. As one of the best chicken breeds with adaptability to highland environment, the Tibetan chicken, is genetically different from lowland chicken breeds. In order to gain a better understanding of the mechanism of hypoxic adaptability in high altitude, in the present study, we focused on the MT-COI together with ATP-6 gene to explore the regulatory mechanisms for hypoxia adaptability in Tibet chicken. Here, we sequenced MT-COI of 29 Tibetan chickens and 30 Chinese domestic chickens and ATP-6 gene of 28 Tibetan chickens and 29 Chinese domestic chickens. In MT-COI gene, 9 single nucleotide polymorphisms (SNPs) were detected though none of these was a missense mutation, confirming the fact that MT-COI gene is a largely conservative sequence. In ATP-6 gene, 6 single nucleotide polymorphisms (SNPs) were detected and we found a missense mutation (m.9441G > A) in the ATP-6 gene of Tibetan chicken resulting in an amino acid substitution. Due to the critical role of ATP-6 gene in the proton translocation and energy metabolism, we speculated the possibility of this mutation playing an important role in easier energy conversion and metabolism in Tibetan chickens than Chinese domestic chickens so as to better adapt to the harsh environment of the high-altitude areas. The Median-joining profile also suggested that haplotype Ha2 has the ancestral position to the other haplotypes and has significant relationship with high-altitude adaptation in ATP-6 gene. Therefore, we considered that the polymorphism (m.9441G > A) in the ATP-6 gene may affect the specific functions of ATP-6 enzyme relating to high-altitude adaptation of Tibetan chicken and MT-COI gene is a largely conservative sequence.

High-altitude ancestry and hypoxia acclimation have distinct effects on exercise capacity and muscle phenotype in deer mice

The hypoxic and cold environment at high altitudes requires that small mammals sustain high rates of O2 transport for exercise and thermogenesis while facing a diminished O2 availability. We used laboratory-born and -raised deer mice (Peromyscus maniculatus) from highland and lowland populations to determine the interactive effects of ancestry and hypoxia acclimation on exercise performance. Maximal O₂consumption (V̇o(2max)) during exercise in hypoxia increased after hypoxia acclimation (equivalent to the hypoxia at ∼4,300 m elevation for 6-8 wk) and was consistently greater in highlanders than in lowlanders. V̇o(2max) during exercise in normoxia was not affected by ancestry or acclimation. Highlanders also had consistently greater capillarity, oxidative fiber density, and maximal activities of oxidative enzymes (cytochrome c oxidase and citrate synthase) in the gastrocnemius muscle, lower lactate dehydrogenase activity in the gastrocnemius, and greater cytochrome c oxidase activity in the diaphragm. Hypoxia acclimation did not affect any of these muscle traits. The unique gastrocnemius phenotype of highlanders was associated with higher mRNA and protein abundances of peroxisome proliferator-activated receptor γ (PPARγ). Vascular endothelial growth factor (VEGFA) transcript abundance was lower in highlanders, and hypoxia acclimation reduced the expression of numerous genes that regulate angiogenesis and energy metabolism, in contrast to the observed population differences in muscle phenotype. Lowlanders exhibited greater increases in blood hemoglobin content, hematocrit, and wet lung mass (but not dry lung mass) than highlanders after hypoxia acclimation. Genotypic adaptation to high altitude, therefore, improves exercise performance in hypoxia by mechanisms that are at least partially distinct from those underlying hypoxia acclimation.

Oxygen in demand: How oxygen has shaped vertebrate physiology

In response to varying environmental and physiological challenges, vertebrates have evolved complex and often overlapping systems. These systems detect changes in environmental oxygen availability and respond by increasing oxygen supply to the tissues and/or by decreasing oxygen demand at the cellular level. This suite of responses is termed the oxygen transport cascade and is comprised of several components. These components include 1) chemosensory detectors that sense changes in oxygen, carbon dioxide, and pH in the blood, and initiate changes in 2) ventilation and 3) cardiac work, thereby altering the rate of oxygen delivery to, and carbon dioxide clearance from, the tissues. In addition, changes in 4) cellular and systemic metabolism alters tissue-level metabolic demand. Thus the need for oxygen can be managed locally when increasing oxygen supply is not sufficient or possible. Together, these mechanisms provide a spectrum of responses that facilitate the maintenance of systemic oxygen homeostasis in the face of environmental hypoxia or physiological oxygen depletion (i.e. due to exercise or disease). Bill Milsom has dedicated his career to the study of these responses across phylogenies, repeatedly demonstrating the power of applying the comparative approach to physiological questions. The focus of this review is to discuss the anatomy, signalling pathways, and mechanics of each step of the oxygen transport cascade from the perspective of a Milsomite. That is, by taking into account the developmental, physiological, and evolutionary components of questions related to oxygen transport. We also highlight examples of some of the remarkable species that have captured Bill's attention through their unique adaptations in multiple components of the oxygen transport cascade, which allow them to achieve astounding physiological feats. Bill's research examining the oxygen transport cascade has provided important insight and leadership to the study of the diverse suite of adaptations that maintain cellular oxygen content across vertebrate taxa, which underscores the value of the comparative approach to the study of physiological systems.

Comparison of whole mitochondrial genome sequences from two clades of the invasive ascidian, Didemnum vexillum

The mitochondria are the main source of cellular energy production and have an important role in development, fertility, and thermal limitations. Adaptive mitochondrial DNA mutations have the potential to be of great importance in determining aspects of the life history of an organism. Phylogenetic analyses of the globally invasive marine ascidian Didemnum vexillum using the mitochondrial cytochrome c oxidase 1 (COX1) coding region, revealed two distinct clades. Representatives of one clade (denoted by 'B') are geographically restricted to D. vexillum's native region (north-west Pacific Ocean, including Japan), whereas members of the other clade (denoted by 'A') have been introduced and become invasive in temperate coastal areas around the world. Persistence of clade B's restricted distribution may reflect it being inherently less invasive than clade A. To investigate this we sought to determine if the two clades differ significantly in other mitochondrial genes of functional significance, specifically, alterations in amino acids encoded in mitochondrial enzyme subunits. Differences in functional mitochondrial genes could indicate an increased ability for clade A colonies to tolerate a wider range of environmental temperature. Full mitochondrial genomic sequences from D. vexillum clades A and B were obtained and they predict significant sequence differences in genes encoding for enzymes involved in oxidative phosphorylation. Diversity levels were relatively high and showed divergence across almost all genes, with p-distance values between the two clades indicating recent divergence. Both clades showed an excess of rare variants, which is consistent with balancing selection or a recent population expansion. Results presented here will inform future research focusing on examining the functional properties of the corresponding mitochondrial respiration enzymes, of A and B clade enzymes. By comparing closely related taxa that have differing distributions it is possible to identify genes and phenotypes suited to particular environments. The examination of mitochondrial genotypes, and associated enzyme functioning, across populations may aid in our understanding of thermal tolerance and environmental adaptation.

Evolutionary neutrality of mtDNA introgression: evidence from complete mitogenome analysis in roe deer

Introgressive hybridization offers a unique platform for studying the molecular basis of natural selection acting on mitogenomes. Most of the mtDNA protein-coding genes are extremely conserved; however, some of the observed variations have potentially adaptive significance. Here, we evaluated whether the evolution of mtDNA in closely related roe deer species affected by widespread mtDNA introgression is neutral or adaptive. We characterized and compared 16 complete mitogenomes of European (Capreolus capreolus) and Siberian (C. pygargus) roe deer, including four of Siberian origin introgressed into European species. The average sequence divergence of species-specific lineages was estimated at 2.8% and varied across gene classes. Only 21 of 315 fixed differences identified in protein-coding genes represented nonsynonymous changes. Only three of them were determined to have arisen in the C. pygargus lineage since the time to the most recent common ancestor (TMRCA) of both Capreolus species, reflecting a decelerated evolutionary ratio. The almost four-fold higher dN /dS ratio described for the European roe deer lineage is constrained by overall purifying selection, especially pronounced in the ND4 and ND5 genes. We suggest that the highly divergent C. capreolus lineage could have maintained a capability for genomic incorporation of the well-preserved and almost ancestral type of mtDNA present in C. pygargus. Our analyses did not indicate any signs of positive selection for Siberian roe deer mtDNA, suggesting that the present widespread introgression is evolutionarily neutral.

Mitogenomic analyses propose positive selection in mitochondrial genes for high-altitude adaptation in galliform birds

Galliform birds inhabit very diverse habitats, including plateaus that are above 3000 m in altitude. At high altitude, lower temperature and hypoxia are two important factors influencing survival. Mitochondria, as the ultimate oxygen transductor, play an important role in aerobic respiration through oxidative phosphorylation (OXPHOS). We analyzed the mitochondrial genomes of six high-altitude phasianidae birds and sixteen low-altitude relatives in an attempt to determine the role of mitochondrial genes in high-altitude adaptation. We reconstructed the phylogenetic relationships of these phasianidae birds and relatives and found at least four lineages that independently occupied this high-altitude habitat. Selective analyses revealed significant evidence for positive selection in the genes ND2, ND4, and ATP6 in three of the high-altitude lineages. This result strongly suggests that adaptive evolution of mitochondrial genes played a critical role during the independent acclimatization to high altitude by galliform birds.

Positive correlation of trophic level and proportion of sexual taxa of oribatid mites (Acari: Oribatida) in alpine soil systems

We investigated community structure, trophic ecology (using stable isotope ratios; (15)N/(14)N, (13)C/(12)C) and reproductive mode of oribatid mites (Acari, Oribatida) along an altitudinal gradient (2,050-2,900 m) in the Central Alps (Obergurgl, Austria). We hypothesized that (1) the community structure changes with altitude, (2) oribatid mites span over four trophic levels, (3) the proportion of sexual taxa increases with altitude, and (4) the proportion of sexual taxa increases with trophic level, i.e. is positively correlated with the δ(15)N signatures. Oribatid mite community structure changed with altitude indicating that oribatid mites occupy different niches at different altitudes. Oribatid mites spanned over 12 δ(15)N units, i.e. about four trophic levels, which is similar to lowland forest ecosystems. The proportion of sexually reproducing taxa increased from 2,050 to 2,900 m suggesting that limited resource availability at high altitudes favors sexual reproduction. Sexual taxa more frequently occurred higher in the food web indicating that the reproductive mode is related to nutrition of oribatid mites. Generally, oribatid mite community structure changed from being decomposer dominated at lower altitude to being dominated by fungal and lichen feeders, and predators at higher altitude. This supports the view that resources from dead organic material become less available with increasing altitude forcing species to feed on living resources such as fungi, lichens and nematodes. Our findings support the hypothesis that limited resource accessibility (at high altitudes) favors sexually reproducing species whereas ample resource supply (at lower altitudes) favors parthenogenetic species.

Fin whale MDH-1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH-1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH-1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.

The elusive nature of adaptive mitochondrial DNA evolution of an arctic lineage prone to frequent introgression

Mitochondria play a fundamental role in cellular metabolism, being responsible for most of the energy production of the cell in the oxidative phosphorylation (OXPHOS) pathway. Mitochondrial DNA (mtDNA) encodes for key components of this process, but its direct role in adaptation remains far from understood. Hares (Lepus spp.) are privileged models to study the impact of natural selection on mitogenomic evolution because 1) species are adapted to contrasting environments, including arctic, with different metabolic pressures, and 2) mtDNA introgression from arctic into temperate species is widespread. Here, we analyzed the sequences of 11 complete mitogenomes (ten newly obtained) of hares of temperate and arctic origins (including two of arctic origin introgressed into temperate species). The analysis of patterns of codon substitutions along the reconstructed phylogeny showed evidence for positive selection in several codons in genes of the OXPHOS complexes, most notably affecting the arctic lineage. However, using theoretical models, no predictable effect of these differences was found on the structure and physicochemical properties of the encoded proteins, suggesting that the focus of selection may lie on complex interactions with nuclear encoded peptides. Also, a cloverleaf structure was detected in the control region only from the arctic mtDNA lineage, which may influence mtDNA replication and transcription. These results suggest that adaptation impacted the evolution of hare mtDNA and may have influenced the occurrence and consequences of the many reported cases of massive mtDNA introgression. However, the origin of adaptation remains elusive.

Mitochondrial genetics and obesity: evolutionary adaptation and contemporary disease susceptibility

Obesity is a leading risk factor for a variety of metabolic diseases including cardiovascular disease, diabetes, and cancer. Although in its simplest terms, obesity may be thought of as a consequence of excessive caloric intake and sedentary lifestyle, it is also evident that individual propensity for weight gain can vary. The etiology of individual susceptibility to obesity seems to be complex-involving a combination of environmental-genetic interactions. Herein, we suggest that the mitochondrion plays a major role in influencing individual susceptibility to this disease via mitochondrial-nuclear interaction processes and that environmentally influenced selection events for mitochondrial function that conveyed increased reproductive and survival success during the global establishment of human populations during prehistoric times can influence individual susceptibility to weight gain and obesity.

Efficient utilization of aerobic metabolism helps Tibetan locusts conquer hypoxia

Background

Responses to hypoxia have been investigated in many species; however, comparative studies between conspecific geographical populations at different altitudes are rare, especially for invertebrates. The migratory locust, Locusta migratoria, is widely distributed around the world, including on the high-altitude Tibetan Plateau (TP) and the low-altitude North China Plain (NP). TP locusts have inhabited Tibetan Plateau for over 34,000 years and thus probably have evolved superior capacity to cope with hypoxia.

Results

Here we compared the hypoxic responses of TP and NP locusts from morphological, behavioral, and physiological perspectives. We found that TP locusts were more tolerant of extreme hypoxia than NP locusts. To evaluate why TP locusts respond to extreme hypoxia differently from NP locusts, we subjected them to extreme hypoxia and compared their transcriptional responses. We found that the aerobic metabolism was less affected in TP locusts than in NP locusts. RNAi disruption of PDHE1β, an entry gene from glycolysis to TCA cycle, increased the ratio of stupor in TP locusts and decreased the ATP content of TP locusts in hypoxia, confirming that aerobic metabolism is critical for TP locusts to maintain activity in hypoxia.

Conclusions

Our results indicate that TP and NP locusts have undergone divergence in hypoxia tolerance. These findings also indicate that insects can adapt to hypoxic pressure by modulating basic metabolic processes.

Beyond retrograde and anterograde signalling: mitochondrial-nuclear interactions as a means for evolutionary adaptation and contemporary disease susceptibility

Although there is general agreement that most forms of common disease develop as a consequence of a combination of factors, including genetic, environmental and behavioural contributors, the actual mechanistic basis of how these factors initiate or promote diabetes, cancer, neurodegenerative and cardiovascular diseases in some individuals but not in others with seemingly identical risk factor profiles, is not clearly understood. In this respect, consideration of the potential role for mitochondrial genetics, damage and function in influencing common disease susceptibility seems merited, given that the prehistoric challenges were the original factors that moulded cellular function, and these were based upon the mitochondrial-nuclear relationships that were established during evolutionary history. These interactions were probably refined during prehistoric environmental selection events that, at present, are largely absent. Contemporary risk factors such as diet, sedentary lifestyle and increased longevity, which influence our susceptibility to a variety of chronic diseases were not part of the dynamics that defined the processes of mitochondrial-nuclear interaction, and thus cell function. Consequently, the prehistoric challenges that contributed to cell functionality and evolution should be considered when interpreting and designing experimental data and strategies. Although several molecular epidemiological studies have generally supported this notion, studies that probe beyond these associations are required. Such investigation will mark the initial steps for mechanistically addressing the provocative concept that contemporary human disease susceptibility is the result of prehistoric selection events for mitochondrial-nuclear function, which increased the probability for survival and reproductive success during evolution.

Association between MT-CO3 haplotypes and high-altitude adaptation in Tibetan chicken

Genetic mutation in cytochrome c oxidase subunit III gene (MT-CO3) could influence the kinetics of cytochrome c oxidase (COX), which catalyzes oxygen transport capacity in oxidative phosphorylation. However, the potential relationship between MT-CO3 variants and high-altitude adaptation remains poorly understood in Tibetan chicken. Here, we sequenced MT-CO3 gene of 125 Tibetan chickens and 144 Chinese domestic chickens in areas at a low elevation (below 1,000 m). Eight single nucleotide polymorphisms (SNPs) were detected; and five of them (m.10081A>G, m.10115G>A, m.10270G>A, m.10336A>G and m.10447C>T) shared by Tibetan chicken and lowland chicken with the significant difference in their respective allele frequencies. Nine haplotypes (H1-H9) were finally defined. Among them, haplotype H4 was positively associated with high-altitude adaptation whereas haplotypes H6, H7 and H8 had negative association with high-altitude adaptation. The Median-joining profile suggested that haplotype H5 had the ancestral position to the other haplotypes but had no significant relationship with high-altitude adaptation. However, there was only m.10081A>G mutation differed from haplotype H4 and H5. Results also suggested that chickens with A allele at m.10081A>G, had over 2.6 times than those with G allele in the probability of the ability to adapt hypoxia. It suggests that the synonymous mutation m.10081A>G may be a prerequisite for shaping high-altitude adaptation-specific haplotypes.

Genome sequence of ground tit Pseudopodoces humilis and its adaptation to high altitude

Background

The mechanism of high-altitude adaptation has been studied in certain mammals. However, in avian species like the ground tit Pseudopodoces humilis, the adaptation mechanism remains unclear. The phylogeny of the ground tit is also controversial.

Results

Using next generation sequencing technology, we generated and assembled a draft genome sequence of the ground tit. The assembly contained 1.04 Gb of sequence that covered 95.4% of the whole genome and had higher N50 values, at the level of both scaffolds and contigs, than other sequenced avian genomes. About 1.7 million SNPs were detected, 16,998 protein-coding genes were predicted and 7% of the genome was identified as repeat sequences. Comparisons between the ground tit genome and other avian genomes revealed a conserved genome structure and confirmed the phylogeny of ground tit as not belonging to the Corvidae family. Gene family expansion and positively selected gene analysis revealed genes that were related to cardiac function. Our findings contribute to our understanding of the adaptation of this species to extreme environmental living conditions.

Conclusions

Our data and analysis contribute to the study of avian evolutionary history and provide new insights into the adaptation mechanisms to extreme conditions in animals.

Functional modulation of mitochondrial cytochrome c oxidase underlies adaptation to high-altitude hypoxia in a Tibetan migratory locust

Mitochondria are crucial to the hypoxia response of aerobic organisms. However, mitochondrial mechanisms for hypoxia adaptation remain largely unknown. We conducted a comparative study on the mitochondrial hypoxia response and adaptation of the Tibetan Plateau and North China lowland populations of migratory locusts, Locusta migratoria. Compared with lowland locusts, Tibetan locusts presented significantly higher hypoxia tolerance and a better-maintained mitochondrial structure in flight muscles under oxygen partial pressure of 1.6 kPa. The hypoxic treatment inhibited the NADH-linked oxidative phosphorylation (OXPHOS) significantly in both locust populations, but to a less extent in Tibetan locusts. Among the critical components of OXPHOS, only cytochrome c oxidase (COX) exhibited significantly higher activity in Tibetan locusts under normoxia and hypoxia. Pharmacological interventions using NaN(3) confirmed that COX activity inhibition reduced hypoxia tolerance by downregulating OXPHOS in both locust populations. The enhanced COX activity was caused not by protein content, but by elevated catalytic efficiency resulting from the increased ferrocytochrome c affinity of COX and the increased electron transport rate via catalytic redox centres. These findings reveal a novel mechanism that confers mitochondrial robustness against hypoxia by modulating the COX activity, which represents an adaptation to permanent hypoxia in the Tibetan Plateau.

Stepwise colonization of the Andes by ruddy ducks and the evolution of novel β-globin variants

Andean uplift played a key role in Neotropical bird diversification, yet past dispersal and genetic adaptation to high-altitude environments remain little understood. Here we use multilocus population genetics to study population history and historical demographic processes in the ruddy duck (Oxyura jamaicensis), a stiff-tailed diving duck comprising three subspecies distributed from Canada to Tierra del Fuego and inhabiting wetlands from sea level to 4500 m in the Andes. We sequenced the mitochondrial DNA, four autosomal introns and three haemoglobin genes (α(A), α(D), β(A)) and used isolation-with-migration (IM) models to study gene flow between North America and South America, and between the tropical and southern Andes. Our analyses indicated that ruddy ducks dispersed first from North America to the tropical Andes, then from the tropical Andes to the southern Andes. While no nonsynonymous substitutions were found in either α globin gene, three amino acid substitutions were observed in the β(A) globin. Based on phylogenetic reconstruction and power analysis, the first β(A) substitution, found in all Andean individuals, was acquired when ruddy ducks dispersed from low altitude in North America to high altitude in the tropical Andes, whereas the two additional substitutions occurred more recently, when ruddy ducks dispersed from high altitude in the tropical Andes to low altitude in the southern Andes. This stepwise colonization pattern accompanied by polarized β(A) globin amino acid replacements suggest that ruddy ducks first acclimatized or adapted to the Andean highlands and then again to the lowlands. In addition, ruddy ducks colonized the Andean highlands via a less common route as compared to other waterbird species that colonized the Andes northwards from the southern cone of South America.

The paradox of extreme high-altitude migration in bar-headed geese Anser indicus

Bar-headed geese are renowned for migratory flights at extremely high altitudes over the world's tallest mountains, the Himalayas, where partial pressure of oxygen is dramatically reduced while flight costs, in terms of rate of oxygen consumption, are greatly increased. Such a mismatch is paradoxical, and it is not clear why geese might fly higher than is absolutely necessary. In addition, direct empirical measurements of high-altitude flight are lacking. We test whether migrating bar-headed geese actually minimize flight altitude and make use of favourable winds to reduce flight costs. By tracking 91 geese, we show that these birds typically travel through the valleys of the Himalayas and not over the summits. We report maximum flight altitudes of 7290 m and 6540 m for southbound and northbound geese, respectively, but with 95 per cent of locations received from less than 5489 m. Geese travelled along a route that was 112 km longer than the great circle (shortest distance) route, with transit ground speeds suggesting that they rarely profited from tailwinds. Bar-headed geese from these eastern populations generally travel only as high as the terrain beneath them dictates and rarely in profitable winds. Nevertheless, their migration represents an enormous challenge in conditions where humans and other mammals are only able to operate at levels well below their sea-level maxima.

Genetic and phenotypic divergence between low- and high-altitude populations of two recently diverged cinnamon teal subspecies

Spatial variation in the environment can lead to divergent selection between populations occupying different parts of a species' range, and ultimately lead to population divergence. The colonization of new areas can thus facilitate divergence in beneficial traits, yet with little differentiation at neutral genetic markers. We investigated genetic and phenotypic patterns of divergence between low- and high-altitude populations of cinnamon teal inhabiting normoxic and hypoxic regions in the Andes and adjacent lowlands of South America. Cinnamon teal showed strong divergence in body size (PC1; P(ST) = 0.56) and exhibited significant frequency differences in a single nonsynonymous α-hemoglobin amino acid polymorphism (Asn/Ser-α9; F(ST) = 0.60) between environmental extremes, despite considerable admixture of mtDNA and intron loci (F(ST) = 0.004-0.168). Inferences of strong population segregation were further supported by the observation of few mismatched individuals in either environmental extreme. Coalescent analyses indicated that the highlands were most likely colonized from lowland regions but following divergence, gene flow has been asymmetric from the highlands into the lowlands. Multiple selection pressures associated with high-altitude habitats, including cold and hypoxia, have likely shaped morphological and genetic divergence within South American cinnamon teal populations.

Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and nonparallel physiological divergence

Across multiple lakes in North America, lake whitefish (Coregonus clupeaformis) have independently evolved 'dwarf' and 'normal' sympatric species pairs that exhibit pronounced phenotypic and genetic divergence. In particular, traits associated with metabolism have been shown to be highly differentiated between whitefish species. Here, we examine the transcription of genes associated with the five mitochondrial and nuclear genome-encoded oxidative phosphorylation (OXPHOS) complexes, the primary physiological mechanism responsible for the production of ATP, in whitefish species pairs from Cliff Lake and Webster Lake in Maine, USA. We observed OXPHOS gene transcription divergence between dwarf and normal whitefish in each of the two lakes, with the former exhibiting transcription upregulation for genes associated with each of the OXPHOS complexes. We also observed a significant influence of lake on transcription levels for some of the genes, indicating that inter-lake ecological or genetic differences are contributing to variation in OXPHOS gene transcription levels. Together, our results support the hypothesis that metabolic divergence is a critical adaptation involved in whitefish speciation and implicate OXPHOS gene upregulation as a factor involved in meeting the enhanced energetic demands of dwarf whitefish. Further studies are now needed to evaluate the contribution of genetically vs. plasticity driven variation in transcription associated with this critical physiological pathway.

Genomic insights into adaptation to high-altitude environments

Elucidating the molecular genetic basis of adaptive traits is a central goal of evolutionary genetics. The cold, hypoxic conditions of high-altitude habitats impose severe metabolic demands on endothermic vertebrates, and understanding how high-altitude endotherms cope with the combined effects of hypoxia and cold can provide important insights into the process of adaptive evolution. The physiological responses to high-altitude stress have been the subject of over a century of research, and recent advances in genomic technologies have opened up exciting opportunities to explore the molecular genetic basis of adaptive physiological traits. Here, we review recent literature on the use of genomic approaches to study adaptation to high-altitude hypoxia in terrestrial vertebrates, and explore opportunities provided by newly developed technologies to address unanswered questions in high-altitude adaptation at a genomic scale.

Rcf1 and Rcf2, members of the hypoxia-induced gene 1 protein family, are critical components of the mitochondrial cytochrome bc1-cytochrome c oxidase supercomplex

We report that Rcf1 (formerly Aim31), a member of the conserved hypoxia-induced gene 1 (Hig1) protein family, represents a novel component of the yeast cytochrome bc(1)-cytochrome c oxidase (COX) supercomplex. Rcf1 (respiratory supercomplex factor 1) partitions with the COX complex, and evidence that it may act as a bridge to the cytochrome bc(1) complex is presented. Rcf1 interacts with the Cox3 subunit and can do so prior to their assembly into the COX complex. A close proximity of Rcf1 and members of the ADP/ATP carrier (AAC) family was also established. Rcf1 displays overlapping function with another Hig1-related protein, Rcf2 (formerly Aim38), and their joint presence is required for optimal COX enzyme activity and the correct assembly of the cytochrome bc(1)-COX supercomplex. Rcf1 and Rcf2 can independently associate with the cytochrome bc(1)-COX supercomplex, indicating that at least two forms of this supercomplex exist within mitochondria. We provide evidence that the association with the cytochrome bc(1)-COX supercomplex and regulation of the COX complex are a conserved feature of Hig1 family members. Based on our findings, we propose a model where the Hig1 proteins regulate the COX enzyme activity through Cox3 and associated Cox12 protein, in a manner that may be influenced by the neighboring AAC proteins.

Elevated performance: the unique physiology of birds that fly at high altitudes

Birds that fly at high altitudes must support vigorous exercise in oxygen-thin environments. Here I discuss the characteristics that help high fliers sustain the high rates of metabolism needed for flight at elevation. Many traits in the O(2) transport pathway distinguish birds in general from other vertebrates. These include enhanced gas-exchange efficiency in the lungs, maintenance of O(2) delivery and oxygenation in the brain during hypoxia, augmented O(2) diffusion capacity in peripheral tissues and a high aerobic capacity. These traits are not high-altitude adaptations, because they are also characteristic of lowland birds, but are nonetheless important for hypoxia tolerance and exercise capacity. However, unique specializations also appear to have arisen, presumably by high-altitude adaptation, at every step in the O(2) pathway of highland species. The distinctive features of high fliers include an enhanced hypoxic ventilatory response, an effective breathing pattern, larger lungs, haemoglobin with a higher O(2) affinity, further augmentation of O(2) diffusion capacity in the periphery and multiple alterations in the metabolic properties of cardiac and skeletal muscle. These unique specializations improve the uptake, circulation and efficient utilization of O(2) during high-altitude hypoxia. High-altitude birds also have larger wings than their lowland relatives to reduce the metabolic costs of staying aloft in low-density air. High fliers are therefore unique in many ways, but the relative roles of adaptation and plasticity (acclimatization) in high-altitude flight are still unclear. Disentangling these roles will be instrumental if we are to understand the physiological basis of altitudinal range limits and how they might shift in response to climate change.

The mitochondrial paradigm for cardiovascular disease susceptibility and cellular function: a complementary concept to Mendelian genetics

While there is general agreement that cardiovascular disease (CVD) development is influenced by a combination of genetic, environmental, and behavioral contributors, the actual mechanistic basis of how these factors initiate or promote CVD development in some individuals while others with identical risk profiles do not, is not clearly understood. This review considers the potential role for mitochondrial genetics and function in determining CVD susceptibility from the standpoint that the original features that molded cellular function were based upon mitochondrial-nuclear relationships established millions of years ago and were likely refined during prehistoric environmental selection events that today, are largely absent. Consequently, contemporary risk factors that influence our susceptibility to a variety of age-related diseases, including CVD were probably not part of the dynamics that defined the processes of mitochondrial-nuclear interaction, and thus, cell function. In this regard, the selective conditions that contributed to cellular functionality and evolution should be given more consideration when interpreting and designing experimental data and strategies. Finally, future studies that probe beyond epidemiologic associations are required. These studies will serve as the initial steps for addressing the provocative concept that contemporary human disease susceptibility is the result of selection events for mitochondrial function that increased chances for prehistoric human survival and reproductive success.

Cytochrome c oxidase: evolution of control via nuclear subunit addition

According to theory, present eukaryotic cells originated from a beneficial association between two free-living cells. Due to this endosymbiotic event the pre-eukaryotic cell gained access to oxidative phosphorylation (OXPHOS), which produces more than 15 times as much ATP as glycolysis. Because cellular ATP needs fluctuate and OXPHOS both requires and produces entities that can be toxic for eukaryotic cells such as ROS or NADH, we propose that the success of endosymbiosis has largely depended on the regulation of endosymbiont OXPHOS. Several studies have presented cytochrome c oxidase as a key regulator of OXPHOS; for example, COX is the only complex of mammalian OXPHOS with known tissue-specific isoforms of nuclear encoded subunits. We here discuss current knowledge about the origin of nuclear encoded subunits and the appearance of different isozymes promoted by tissue and cellular environments such as hypoxia. We also review evidence for recent selective pressure acting on COX among vertebrates, particularly in primate lineages, and discuss the unique pattern of co-evolution between the nuclear and mitochondrial genomes. Finally, even though the addition of nuclear encoded subunits was a major event in eukaryotic COX evolution, this does not lead to emergence of a more efficient COX, as might be expected from an anthropocentric point of view, for the "higher" organism possessing large brains and muscles. The main function of these subunits appears to be "only" to control the activity of the mitochondrial subunits. We propose that this control function is an as yet under appreciated key point of evolution. Moreover, the importance of regulating energy supply may have caused the addition of subunits encoded by the nucleus in a process comparable to a "domestication scenario" such that the host tends to control more and more tightly the ancestral activity of COX performed by the mtDNA encoded subunits.