Skip to main content

The adaptive evolution of the mammalian mitochondrial genome



The mitochondria produce up to 95% of a eukaryotic cell's energy through oxidative phosphorylation. The proteins involved in this vital process are under high functional constraints. However, metabolic requirements vary across species, potentially modifying selective pressures. We evaluate the adaptive evolution of 12 protein-coding mitochondrial genes in 41 placental mammalian species by assessing amino acid sequence variation and exploring the functional implications of observed variation in secondary and tertiary protein structures.


Wide variation in the properties of amino acids were observed at functionally important regions of cytochrome b in species with more-specialized metabolic requirements (such as adaptation to low energy diet or large body size, such as in elephant, dugong, sloth, and pangolin, and adaptation to unusual oxygen requirements, for example diving in cetaceans, flying in bats, and living at high altitudes in alpacas). Signatures of adaptive variation in the NADH dehydrogenase complex were restricted to the loop regions of the transmembrane units which likely function as protons pumps. Evidence of adaptive variation in the cytochrome c oxidase complex was observed mostly at the interface between the mitochondrial and nuclear-encoded subunits, perhaps evidence of co-evolution. The ATP8 subunit, which has an important role in the assembly of F0, exhibited the highest signal of adaptive variation. ATP6, which has an essential role in rotor performance, showed a high adaptive variation in predicted loop areas.


Our study provides insight into the adaptive evolution of the mtDNA genome in mammals and its implications for the molecular mechanism of oxidative phosphorylation. We present a framework for future experimental characterization of the impact of specific mutations in the function, physiology, and interactions of the mtDNA encoded proteins involved in oxidative phosphorylation.


Mitochondrial DNA (mtDNA) has long been treated as an ideal marker because of its convenience for reconstruction of gene genealogy and population history inference. However, the selective neutrality assumption of the mtDNA is simplistic because variation in mitochondrial protein-coding genes involved in oxidative phosphorylation (responsible for the production of up to 95% of the energy of eukaryotic cells) can directly influence metabolic performance. Because of the importance of this biochemical pathway, evaluating selective pressures acting on mtDNA proteins could provide key insight into the adaptive evolution of the mtDNA genome as has been suggested by recent empirical evidence (e.g. Ruiz-Pesini et al. 2004 [1], Moyer et al. 2005 [2], Bazin et al. 2006 [3]). As amino acid changes cause inefficiencies in the electron transfer chain system, oxidative phosphorylation produces reactive oxygen molecules, causing oxidative damage to mitochondrial and cellular proteins, lipids and nucleic acids, and eventually interrupting the production of mitochondrial energy.

Mutations in mitochondrial-encoded genes can influence the production of reactive oxygen species in mice [4] and have been implicated in a large number of human and mouse diseases [5]. Some amino acid changes may also improve aerobic capacity and adaptation to new thermal environments [611]. Furthermore, mutations in mitochondrial genes have been implicated in exercise intolerance in humans [12] (see Table 1). Metabolic capacity varies widely among mammalian species [13], and variation in the elements of the oxidative phosphorylation pathway have been linked to different life history traits and environmental adaptations [1416]. In addition, functional interactions between mitochondrial- and nuclear-encoded proteins may result in co-adaptation to maintain or improve metabolic fitness [9, 17].

Table 1 Examples of mutations in mitochondrial genes that cause exercise intolerance in humans.

The mtDNA comprises a closed circular DNA strand that encodes the following proteins involved in the oxidative phosphorylation machinery: seven subunits of the NADH dehydrogenase or NADH ubiquinone oxidoreductase complex (ND: ND1, 2, 3, 4, 4L, 5 and 6), the cytochrome b subunit of the ubiquinol cytochrome c oxidoreductase or cytochrome bc1 complex (CytB), three subunits of the cytochrome c oxidase complex (COX), and two subunits of ATP synthase (ATPase: ATP6 and ATP8) [18] (Figure 1A). These key components of four of the five complexes involved in oxidative phosphorylation, combine with other subunits encoded by the nuclear DNA genome [18] (Figure 1B). In the oxidative phosphorylation, reducing equivalents resulting from the oxidation of nutrients such as glucose, are carried by a series of molecules (the electron transport chain), which have increasing standard reduction potentials. The resulting free energy is transformed into a proton gradient that is used by ATP synthase to produce ATP (adenosine triphosphate) from ADP (adenosine diphosphate) [18].

Figure 1

The mammalian mitochondrial genome and its protein-coding gene repertoire involved in the oxidative phosphorylation pathway. (A) Schematic representation of genes within mammalian mitochondrial genome (~7,000 bp). Genes on the outer circle are transcribed from the light-strand. Location of the tRNAs (red boxes) conform to the canonical placental mammalian arrangement. Abbreviations: HSP, putative heavy-strand promoter; OHR, origin of heavy-strand replication; OLR, origin of light-strand replication. (B) Simplified view of the mitochondrial oxidative phosphorylation machinery. Complexes I (NADH dehydrogenase) and II (succinate dehydrogenase) receive electrons from either NADH or FADH2. Electrons are then carried between complexes by the carrier molecules coenzyme Q/ubiquinone (UQ) and cytochrome c (CYC). The potential energy of these electron transfer events is used to pump protons against the gradient, from the mitochondrial matrix into the intermembrane space [Complexes I and III (cytochrome bc1) and IV (cytochrome c oxidase)]. ATP synthesis by Complex V (ATP synthase) is driven by the proton gradient, and occurs in the mitochondrial matrix. MM: mitochondria matrix; IM: intermembrane space.

Although mtDNA sequence data have been extensively used in mammalian phylogenetics [1924], relatively little attention has been devoted to the study of molecular adaptation of mitochondrial encoded proteins. In this work we combine molecular evolution analyses with crystallographic and secondary structure prediction analyses to explore how mitochondrial genetic variation may be linked to the diverse metabolic patterns of 41 mammalian species (Table 2) from each of the four major clades of the placental mammals. Since most of the cell's ATP production results from mitochondrial oxidative phosphorylation, we also use estimates of metabolic rates based on oxygen (O2) consumption under aerobic conditions as a surrogate for physiological or metabolic differences among the species in our analyses.

Table 2 The 41 mammalian species in study.

Results and Discussion

The mitogenomic phylogenetic reconstruction (Figure 2) obtained with Bayesian inference methods [25, 26] is mostly congruent with previous comprehensive analyses from nuclear and mtDNA genes [27] where mammalian species form four major groups: Laurasiatheria, Euarchontoglires, Xenarthra, and Afrotheria. The superordinal tree is resolved with posterior probabilities greater than 0.99. The relative positioning of the orders is in agreement with that obtained using Maximum Likelihood [see Additional file 1, Fig. S2], supporting the previous topology.

Figure 2

Mammalian mitogenomic phylogenetic tree. Consensus phylogenetic tree (50 percent majority rule) constructed from the combined set of MCMC runs resulting topologies in MrBayes. Numbers give the percentage of posterior probability support for each clade averaged over the runs.

Significant physicochemical amino acid changes among residues in mitochondrial protein coding genes were identified by the algorithm implemented in TreeSAAP [28], which compares the observed distribution of physicochemical changes inferred from a phylogenetic tree with an expected distribution based on the assumption of completely random amino acid replacement expected under the condition of selective neutrality (see details in the Methods section). There are more modifications in the number of radical amino acid property changes in the tips (80%) relative to the interior branches (20%) of the mammalian tree (Figure 3). The interior branches with the highest number of radical amino acid property changes are those that lead to the house mouse and the rat (node 78 to 80; a total of 55 changes across all proteins), the branch leading to the elephant shrews (node 45 to 46; 43 changes across all proteins), the branch leading to the guinea pig and the greater cane rat (node 78 to 79; 21 changes across all proteins), and the branch uniting the humpback whale and white-beaked dolphin (node 64 to 65; 17 changes across all proteins). There was an average of 10 changes per interior branch. The proteins with the highest average number of changing properties per site are ND and ATPase (0.7 and 0.8 average changes per site, respectively), while CytB and COX were the lowest (0.3 and 0.2 changes per site, respectively) (Figure 4). Noteworthy, a distinction arises between loop areas and transmembrane domains for all proteins complexes. Loop areas are, in general, more affected by positive selection than transmembrane domains (0.52 property changes showing strong positive selection per site vs 0.39 for the latter).

Figure 3

Radical physicochemical amino acid changes varying across the mammalian mitogenomic phylogenetic tree. Representation of the standardized number of strong positively selected amino acid properties across mitochondrial protein-coding genes varying within the branches of the mammalian mitochondrial tree (see Methods for details).

Figure 4

Radical physicochemical amino acid changes among residues in mammalian mitochondrial protein-coding genes. Number of strong positively selected amino acid properties in the mammalian mtDNA protein-coding genes of the oxidative phosphorylation chain.

Amino acid properties with signals of strong positive selection accumulated at a rate roughly equivalent to the mutation rate of the gene itself (i.e. mutation rate of ATP > ND > CytB > COX; see Lopez et al. 1997 [29]) (Figure 5; [see Additional file 1, Fig. S3]). This correlation is more pronounced for the protein-coding genes with higher mutation rates, such as ATPase and ND, as is most apparent in analyses of the variation along the interior branches, where 20% of the radical amino acid property changes occur [see Additional file 1, Fig. S3]. The best correlation between overall mutation rate and number of sites with radical amino acid changes was observed for NDs, while the existence of several outliers for ATPase slightly reduced the strength of the correlation [see Additional file 1, Fig. S3].

Figure 5

Correlation between amino acid property variation and genetic distance. Correlation between the number of positively selected amino acid properties and the branch length (genetic distance) in mammalian NDs, CytB, COXs, and ATPs protein-coding genes.

The biochemical complexity of the oxidative phosphorylation processes precludes a clear discussion on the functional implications of the amino acid properties that are under selection (Figures 6 and 7). Negative selection dominates the categories of radical changes. Moderate changes (categories 1 and 2) characterize most of the positive selection detected. Properties under positive destabilizing selection (the power to be at the N-terminal, refractive index, long-range non-bonded energy, coil tendencies, compressibility, turn tendencies and the power to be at the C-terminal) will interfere both at a chemical and structural level. However, since we are considering events as diverse as protein-protein interactions, molecular oxygen and proton diffusion and electron-transfer, it is not possible to establish a direct correlation with the previously referred properties. Discussion on varying amino acid properties will therefore be made specifically for those sites which can be mapped on available structural data.

Figure 6

Amino acid properties under positive (blue) and negative (red) selection in mammalian mitochondrial protein-coding genes. Conservative changes correspond to conservative categories 1 and 2 and radical changes to categories 7 and 8 (P ≤ 0.001). (C: chemical; S: structural; O: other [1])

Figure 7

Amino acid properties under positive (blue) and negative (red) selection in transmembrane and loop regions. Conservative changes correspond to conservative categories 1 and 2 and radical changes to categories 7 and 8 (P ≤ 0.001). (C: chemical; S: structural; O: other [1])

NADH dehydrogenase

The sites with the highest number (between 17 and 18) of strong positively selected changes in amino acid properties are within ND2, ND4 and ND5 [see Additional file 1, Fig. S4]. These subunits show a high average in the number of such changes per residue (1.1, 0.6 and 0.8, respectively, for an overall average of 0.5) (Figure 4). NADH dehydrogenase is the first (Figure 1B) and the largest enzyme complex in the respiratory chain. It receives electrons from the oxidation of NADH and provides electrons for reduction of quinone to quinol. This is coupled to the translocation of four protons across the inner membrane, generating an electrochemical proton gradient (Figure 1B). Complex I is an L-shaped complex (as shown in low resolution electron microscopy analysis [30]) that contains all seven mtDNA encoded subunits in a membrane-embedded arm. The subunits of the "peripheral arm" (that projects into the mitochondrial matrix) are encoded by the nuclear DNA genome. Forty-six different units have been identified in complex I from bovine heart mitochondria [31], although the 14 bacterial subunits are the minimum needed for sufficient energy transduction by complex I [32]. The electron transfer chain events occur in the peripheral arm subunits and the proton pumping occurs in the membrane-embedded module. ND1 and ND2 have been suggested to be located at the junction between the peripheral and the membrane arm [32], while ND4 and ND5 should occur at the distal end of the latter [31]. ND2, ND4, and ND5 are suggested to be the actual proton pumping devices because of their sequence homology with a class of Na+/H+ antiporters [32]. The overall variation in these subunits, which have already been assigned some function, is larger than the observed variation in the subunits with still unknown functions. Mutations in these subunits may interfere with the efficiency of the proton-pumping process. This could either occur through chemical changes that hinder/improve the proton translocation or by disrupting/improving the long-range redox-linked conformational changes that are suggested to occur in order to, for example, activate ND5 which is far away from the electron transfer events [32]. The number of TM domains predicted in this study (Figure 8) for ND2, ND4, and ND5 is similar to that present in their bacterial counterparts (which are 11, 12, and 17, respectively) [33]. Figure 8 shows that the sites with higher variation (20 < number of changes < 110) are located only in the loop regions, suggesting that functional constraints are acting upon the TM domains, which would be consistent with their putative proton-pumping function.

Figure 8

Amino acid property variation in NADH dehydrogenase subunits ND2/3/5. Topological assignment of the sites that present a high number of radically changing properties under positive-destabilizing selection in three subunits of the NADH dehydrogenase complex that are suggested to be proton pumping devices. The transmembrane domain average prediction is shown in grey (for details see Material and Methods section). MM: mitochondria matrix; IM: intermembrane space.

Cytochrome b

CytB is an extremely conserved protein, reflecting its fundamental role in energy production in the mitochondria. It catalyses reversible electron transfer from ubiquinol to cytochrome c coupled to proton translocation (Q-cycle [34, 35]). A quinol molecule at the Q0 site donates an electron to cytochrome c via the iron-sulfur protein (ISP) and cytochrome c1 (Figure 9). A second electron passes sequentially through the bL and bH heme ending up in a quinone/semiquinone radical at the Qi site. In a complete Q cycle, two quinol molecules are the Q0 site and one molecule of quinol is regenerated while four protons are translocated across the membrane. Available x-ray data showing cytochrome bc1 inhibitors bound to both Qi and Q0 sites suggest that these mutations have functional consequences. The high degree of conservation among CytB sequences (Figure 4) made it difficult to find evidence of positive selection using TreeSAAP. Sites detected by TreeSAAP were mainly located at the interface between mitochondrial and nuclear-encoded subunits (Figure 9). However, insight into the evolution of this protein was obtained by inspection of the amino acid substitutions observed in the sequence alignment of the various mammalian species.

Figure 9

Three-dimensional representation of relevant variable amino acids in mammalian cytochrome b. Illustrative representation of some of the amino acid variable sites (see Figure 10) located in relevant functional spots of the bovine CytB structure (pdb code: 1PPJ [66]). The prosthetic groups in CytB are represented in black, and the Qi and QP bound inhibitors in orange (ant: antimycin; stig: stigmatellin). Mutations in sites shown in red have been related to exercise intolerance in humans (see Table 1). Mutations in sites shown in pink are under strong positive selection according to the TreeSAAP analysis (over 5 positively selected properties). In site 277, the alanine present in the bovine structure is shown as ball and stick and the van der Waals surface for the arginine found in dugong is also depicted. MM: mitochondria matrix; IM: intermembrane space.

The African savanna elephant Loxodonta Africana has two atypical amino acid replacements at positions 16 and 260 (Figure 10) that are conserved in other members of the Elephantidae family (the African forest elephant Loxodonta cyclotis, the Asian elephant Elephas maximus, and both the extinct woolly mammoth Mammuthus primigenius and the American mastodon Mammut americanum; [see Additional file 1, Fig. S5]). Another uncommon amino acid occurs in most of the Elephantidae group at position 266, but not in the American mastodon, suggesting that this mutation originated less than 24 million years ago, after the divergence of the American mastodon from the other four proboscidean species analyzed [36] [see Additional file 1, Fig. S5]). At the N-terminal site 16, adjacent to the Qi pocket, there is a positively charged residue (lysine) that will affect the binding of ligands at this site (Figure 9). Residues 260 and 266 are located in loop ef, which contributes to the formation of the ISP binding crater. Changes in shape complementarity between the ISP and CytB have been shown to be important for the control mechanism of ISP conformational change, an important feature of cytochrome bc1 mechanism [37]. The aspartate on position 260, a conserved asparagine spot (Figure 10), may affect both the interaction with the ISP (protein recognition event) and other electron transfer events, as it adds an extra negative charge to this interface (Figure 9). Also on loop ef, the shift from the conserved proline on site 266 to a leucine will alter the rigidity of the secondary structure, with implications in the protein-protein interface contacts (Figure 9).

Figure 10

Amino acid variation in functional sites of mammalian cytochrome b. Amino acid variation in the selected sites of CytB presented in Figure 9 across all the mammalian species surveyed in this study.

The cetaceans have an alanine in position 266 (Figure 10), which will have a similar consequence to that referred above. The biochemistry involved in aquatic mammals is still a matter of debate, but it has been found that mitochondria of seals can survive ischemia for much longer than terrestrial mammals [38], suggesting changes in function or regulation. The humpback whale has another interesting mutation. An arginine in site 110 (Figure 10) changes the local net charge, and being located just outside the coordination sphere of heme bH, such a difference can interfere with the electron transfer events (Figure 9). Such alteration is also seen in two of the Chiroptera species, the Ryukyu flying fox and the New Zealand long-tailed bat (Figure 11). Curiously, the sloth, one of the mammals showing a lower metabolic rate than expected given its body size [14], has a positively charged residue on site 111 (a lysine instead of the highly conserved glutamate), a yet more radical change in net electric charge at this location (Figure 10).

Figure 11

Three-dimensional representation of variable amino acid in mammalian cytochrome c oxidase subunits I, II and III. Sites within COX subunits showing a high number of amino acid properties positive selected [see Additional file 1, Fig. S4] mapped into the bovine structure (pdb code: 1V54 [67]). Side-chains with over 5 amino acid properties positive selected are presented as black spheres. Sites showing particular mutational trends are shown in white (see text for further details). Mutations in sites shown in red have been related to exercise intolerance in humans (see Table 1). Prosthetic groups are presented as grey spheres. The side chains of the hydrophobic loop in COXII that is involved in cytochrome c docking and electron transfer are depicted as sticks together with their van der Waals surface. The nuclear encoded subunits of Complex IV are shown as in silver. MM: mitochondria matrix; IM: intermembrane space.

Another region that shows radical amino acid variation is helix cd2, close to the hinge region of ISP (Figure 9). Mutations in the hinge region were shown to have drastic consequences in the catalytic activity of ISP, by hindering the conformational changes that are required for cytochrome bc1 function [39]. Two species have a proline residue on site 158: the greater cane rat, a rodent with spiny fur on the back, and the pangolin, the scaly-anteater, which has a low metabolic rate because of the combination of an invertebrate diet and a large body size [14]. Sites 16, 159, 162 and 263 show an elevated number of amino acid changes (Figure 10), suggestive of adaptive relevance in many mammalian species. Site 277 presents a highly conserved alanine residue. It is placed in the middle of helix F1, within the Q0 pocket (Figure 9). Two peculiar species (dugong and alpaca) with distinct metabolic requirements show very radical amino acid changes at this site. The dugong is an aquatic mammal that is more closely related to elephants than to other marine mammals [40]. It is sometimes referred to as a sea cow because of its strict sea-grass diet, combining several interesting features from the metabolic point of view: a large body size, a low energy diet and aquatic environment adaptation. It has an arginine residue in position 277, which will not only cause extra steric hindrance because of the size of the side chain relative to alanine (in Figure 9 the van der Waals surface of an arginine side chain is presented; it clearly overlaps the stigmatellin binding position), but will also change the binding mode of the ligand, as it is positively charged. Finally, an important change was detected in the alpaca, a domesticated breed of South American camel-like ungulates that lives at an altitude of 3500 to 5000 meters above sea-level presenting metabolic adaptations to the low O2 environment [41]. This species has a proline at site 277, which will drastically alter the local secondary structure, disrupting the alpha helix and therefore changing the shape of the Q0 pocket. Such a mutation is not present in the closest relatives of the alpaca [see Additional file 1, Fig. S5]. Curiously, the Old World members of the Camelidae family (the dromedary Camelus dromedarius and the bactrian camel C. bactrianus) have an aspartate in position 16, instead of the asparagine exhibited by the four species that inhabit South America (alpaca Lama pacos, guanaco L. guanicoe, llama L. glama, and vicuna Vicugna vicugna) [see Additional file 1, Fig. S5]. Several mutations in human CytB have been related to exercise intolerance [12] (Table 1; Figure 9), all of which have similar chemical effects including a change in the net charge around the heme groups and in the binding pockets and disruption of local secondary structures close to the substrate binding areas.

Cytochrome c Oxidase

Complex IV is the terminal electron transfer chain complex that catalyses the electron oxidations of four consecutive reduced cytochrome c molecules and the concomitant reduction of one O2 molecule to water [42]. COXI and II subunits are directly involved in the electron transfer and proton translocation processes, while the other 11 subunits are thought to have regulatory roles. Results obtained with TreeSAAP were mapped on an available 3D structure for COX (Figure 11). COXI and II are the most conserved of all the 12 subunits analyzed (Figure 3; Figure 4; [see Additional file 1, Fig. S4]). Recent studies on taxa that had clearly experienced significant changes in their metabolic needs detected positive selection acting in COXI and COXII, namely in carnivorous plants [7] and in high-performance fish [6]. As demonstrated in high-performance fish COXII subunits, the highly variable sites in mammals are located at the interface between mitochondrial and nuclear-encoded subunits (Table 3), suggesting either an unknown biological role or the occurrence of compensatory mutations in the nuclear subunits (co-evolution). The known functional spots in the complex were found to be very conserved, namely the portion of COXII directly involved in the electron transfer from cytochrome c to the binuclear CuA center (hydrophobic loop consisting of His102-Tyr105 depicted in Figure 11), the negative patch that surrounds it (Asp119, Glu132, Asp139, Glu157, Asp158) [43], the metal binding sites, and the proposed proton/O2/water pathways [4446]. COXIII had more variable sites, which is expected due to the fact that it has no associated redox cofactors. Also, COXIII has been shown to be unnecessary in some bacterial organisms [45]. Some other interesting changes in COX included the K216E mutation (charge inversion) in COXII observed only in cetaceans, and Y182H in COXIII only seen in the hippopotamus and the dugong (Figure 11). Site 185 in COXIII shows a big variation in amino acid type across all species. Finally, the mutation M29K in COXII has been shown to cause exercise intolerance in humans (Table 1; Figure 11).

Table 3 Variable amino acid sites located in COXI/II/III. These sites located on the mtDNA encoded subunits of COX (subunits I, II and III) show significant amino acid properties variation and are in contact with the nuclear encoded subunits of Complex IV (subunits IV, VB, VIA, VIB, VIIA, VIIB, VIIC and VIII).

ATP synthase

The proton-gradient that results from H+ pumping into the intermembrane space, is used by ATPase to synthesize ATP. The proton channel is located in the membrane sector (F0) which is connected to the catalytic component (F1), located on the matrix side of the membrane (Figure 12A). The latter, when separated from the membrane, behaves as a soluble ATPase. Large cooperative conformational changes occur in order to couple the passage of protons through the membrane arm and the production of ATP [47, 48]. In the proposed mechanism for E. coli ATPase, protons that have accumulated in the periplasm enter the assembly via subunit a (corresponding to ATP6 in yeast [49]). One proton binds between two c subunits (corresponding to ATP9 subunits in yeast [49]). In order for the proton to reach the exit channel, the c subunits (in a total of 10 in E. coli), that are arranged as a cylinder, have to rotate, releasing the proton after 10 steps of proton binding. This rotation movement involves the γ and ε subunits, that remain fixed to the top of one set of c subunits. The rotation of γ within the α/β subunits induces conformational changes that release ATP from the alternating catalytic cycles. The ε subunit (homologous to mitochondrial subunit δ and IF1 regulatory protein) is responsible for determining whether complex V acts as a synthase or catalyses the reverse reaction (pumping protons from the cytoplasm/matrix to the periplasm/intermembrane space) at the expenses of ATP hydrolysis. Subunits b and δ (equivalent to mitochondrial OSCP) keep the α/β subunits in a fixed position. The conservation of ATP6 reflects its key role in the coupling of the proton flow with the rotation of the c subunits: as for the ND complex, the sites with higher variation are located only in the predicted loop regions (Figure 12B).

Figure 12

Rotary model for E. coli F 1 F 0 ATPase and variation in the mammalian ATP6 subunit. A) Rotary model for E. coli F1F0 ATPase (see text for details); B) Topological assignment of the sites that present a high number of strong positively selected amino acid properties under positive-destabilizing selection in ATP6 (corresponds to the a subunit in E. coli). The transmembrane domains location is shown in grey (for details see Material and Methods section). The dark grey domain was only predicted by one of the three methods used.

The ATP8 gene encodes a core subunit of the F0 component of ATPase. In its absence, the ATPase in yeast contains no ATP6 subunit, which suggests an important role in the assembly of F0 [50]. Nevertheless, ATP8 subunit has some highly variable sites (Figure 12B), presenting the higher average of radically changing amino acid properties per residue, suggesting some variation of its regulatory role across species.

Oxidative phosphorylation vs metabolic rates

Some mutations in mitochondrial genes are responsible for severe phenotypic effects related to metabolic capacity, such as exercise intolerance in humans [11] (see Table 1). However, the scaling of metabolic rates in relation to the ATP yield through oxidative phosphorylation is not straightforward, as the total metabolic rate is influenced by multiple characteristics varying across species. The scaling of metabolic rates by body size has been discussed for more than a century [13] and remains controversial [13, 5154], although generally metabolic rate increases as body mass increases (see this trend in our dataset: Figure 13 and Table 2). However, there are numerous exceptions such as for large tropical ant and termite predators (e.g. the sloth, the aardvark, some pangolins, tamanduas and armadillos), for which basal rate of metabolism decreases when their size increases [14]. Furthermore, variation in size is accompanied by changes in the metabolic rates inherent to cells (in vitro studies showed that these decline with increasing body mass; [55]). Some of the approaches to scale the whole body basal metabolic rate take into account the scaling of individual organ masses and metabolic rates [56]. However, variation can occur in the activity of oxidative enzymes, such as in the case of acclimatization to altitude by lamas and alpacas [41].

Figure 13

Basal metabolic rate vs log 10 (body mass). The basal metabolic rate (BMR) is presented in three categories (low, intermediate and high). When the BMR value was not available for the species in study, either that of a close relative (indicated in parenthesis) or the average values of several closely related species (AVE) were used (together with the corresponding body mass value). The increase in BMR with the body mass can be easily observed by comparing the average values for each category (dotted lines).

The variation between the metabolic rates in different species is a consequence of multiple factors, including the need to maintain body temperature, the number of mitochondria and the volume densities and/or cristae surface, and the fact that relative organ mass and organ metabolite rate varies interspecifically. For example, in reptiles, the lower metabolic rates, compared to mammals, are due to a combination of smaller internal organs, lower mitochondrial volume and cristae surface densities [55].

The scaling of metabolic rates is thus an intricate issue, and even recent multiple-cause models [52] are flawed [57]. Adding to the complexity of interspecies metabolic rates analysis is the random accumulation of variation in the coding sequences of proteins directly involved in energy production and differential selective pressures that arise as mutations affecting mitochondrial ATP production.


We present a mammalian phylogeny based on variation in protein-coding mtDNA genes among 41 representative species. Sequence analyses were complemented with functional analyses to assess the potential importance of mutations leading to radical changes in the physicochemical properties of the amino acids. Most of the mtDNA protein-coding genes were extremely conserved, reflecting their vital role in oxidative phosphorylation. However, much of the observed variation had plausible adaptive significance.

The ND2, ND4, and ND5 complex I genes showed higher than average adaptive variation, with all of the variable sites located in the assessed loop regions of these putative protons pumps (3D structural data are needed to further confirm these interpretations and to measure the functional implications).

The available high resolution 3D structure of CytB facilitated interpretation of the functional implications of mutations occurring at portions of the protein which resulted in extreme amino acid properties variation in species with peculiar metabolic requirements (such as adaptation to low energy diet vs large body size, namely in elephant, dugong, sloth, and pangolin; and adaptation to extreme O2 requirements, i.e. diving in cetaceans, flying in bats, and high altitudes resistance in alpacas). The adaptive variation in COX was restricted mostly to the interface between mitochondrial and nuclear-encoded subunits, suggesting either co-evolution or some influence in the regulatory role of the latter. Among the ATPase subunits, ATP8 which has an important role in the assembly of F0, showed the highest amount of adaptive variation in this analysis. ATP6, which has an essential role in the ATPase rotor performance, showed a high adaptive variation in predicted loop areas. Interpretation of possible functional roles of these changes is limited, however, by the lack of experimental and structural data for these genes.

Our study provides insight into the adaptive evolution of the mtDNA genome in mammals, which may have facilitated the successful radiation and diversification of mammalian species into different environments and habits. The evidence of positive selection acting in important functional regions of the various mammalian mtDNA proteins provides the framework for future experimental characterization of the impact of specific mutations in the function, physiology, and interactions of the mtDNA encoded proteins involved in the oxidative phosphorylation.


Phylogenetic analyses

A mammalian mitogenomic phylogeny was constructed using 12 of the 13 protein-coding genes of the mtDNA genome of 41 species representative of all mammalian orders (Table 2). The ND6 gene was excluded because it is encoded by the light-strand which has a significantly different base composition from the heavy-chain [58]. Gaps and ambiguous sites adjacent to gaps were removed, resulting in a total alignment of 10,587 nucleotides (3,529 amino acids). The third codon position was excluded from the phylogenetic analysis (7,058 nucleotides were used) because of observed nucleotide saturation [see Additional file 1, Fig. S1].

Bayesian inference methods with Markov chain Monte Carlo (MCMC) sampling were used in MrBayes [25, 26] to assess phylogenetic relationships among the species. We used a General-Time-Reversible substitution model [59] with the invariant site plus gamma options (five categories) after determining the optimal model of sequence substitution with MrModeltest 2.2 [60]. One cold and four incrementally heated chains were run for 2,000,000 generations with chains I = 2, 3, 4, and 5 incrementally heated with heat being 1/(1+ [i-1]T) and T = 0.2. Trees were sampled every 100 generations from the last 1,000,000 generated (well after the chain reached stationarity) and 10,000 trees were used for inferring Bayesian posterior probability. The burn-in fraction performance was evaluated using the program Tracer v1.4 Bayesian methods have been successfully applied to estimation of the tree topology of placental animals using both mitochondrial and nuclear data [27]. A maximum likelihood phylogenetic tree was constructed in PAUP 4.0b10 [61] after determining the optimal model of sequence substitution (TVM+I+G) with Modeltest 3.04 [62].

Adaptive evolution analyses

Selection in protein-coding genes is generally assessed by estimating ω, the ratio between nonsynonymous and synonymous substitution rates (dN/dS) [63]. However, this statistical approach for detecting molecular adaptation is largely biased against even moderately conservative proteins as it does not allow the possibility that adaptation may come in the form of very few amino acid changes. Thus, significant physicochemical amino acid changes among residues in mitochondrial protein coding genes were identified by the algorithm implemented in TreeSAAP [28], which compares the observed distribution of physicochemical changes inferred from a phylogenetic tree with an expected distribution based on the assumption of completely random amino acid replacement expected under the condition of selective neutrality. The evaluation of the magnitude of property change at nonsynonymous residues and their location on a protein 3D-strcuture may provide important insight into the structural and functional consequences of the substitutions [64]. Eight magnitude categories (1 to 8) represent one-step nucleotide changes in a codon and rank the correspondent variation in a property scale of the coded amino acid. Categories 1 to 3 indicate small variation in the amino acid characteristics while categories 6 to 8 represent the most radical substitutions. By accounting for the property changes across the data set, a set of relative frequencies changes for each category is obtained allowing to test the null hypothesis under the assumption of neutral conditions [65]. The categories for which the observed numbers of amino acid replacements in the data set is significantly different from the null model (z-scores > 1.645; P < 0.05) are considered as being potentially affected by selective pressures [65]. Here we focus on amino acid differences that correspond to radical physicochemical variation (positive-destabilizing selection) and are expected to be linked with significant changes in function. TreeSAAP categorizes each amino acid site by positive and negatively destabilizing using 31 properties (henceforth amino acid positions will be referred as sites). To detect strong directional selective pressure, only changes corresponding to categories 7 and 8 (the 2 most radical property changes categories) at the P ≤ 0.001 level were considered. The total number of changes per site is the sum of those occurring in each branch of the phylogeny. The number of changes in amino acid properties was standardized relatively to the overall size of the protein when comparing different complexes (weight factor = total number of amino acids in the complex/total number of amino acids in ATPase, which is the smallest protein complex).

Protein structure analyses

The functional relevance of the amino acid mutations was discussed in the context of existing three-dimensional (3D) structures of mtDNA encoded proteins (CytB [PDB:1PPJ] [66]; COX [PDB:1V54] [67]). For those proteins with unknown 3D structures (ND and ATPase), topologies for transmembrane (TM) subunits were predicted using hidden Markov model (HMM) based servers for topology prediction of transmembrane proteins [6871]. The algorithm used by the program PRODIV-TMHMM [70] has proven to be very reliable at predicting 3D topologies as it incorporates evolutionary information from multiple sequence alignments and assigns amino acid residues to different TM regions according to their properties. However, since even homologous sequences from the same protein family can have inverted topologies [72], some caution is necessary when using topologies predicted by these automated approaches. We have therefore delineated putative TM domains by integrating the results from PRODIV-TMHMM with two other reliable HMM based methods [73], TMHMM [68, 69] and HMMTOP [71]. Graphic representations of the 3D structures were created with the program VMD [74].



mitochondrial DNA


adenosine triphosphate


adenosine diphosphate


nicotinamide adenine dinucleotide


NADH dehydrogenase


cytochrome b


cytochrome c oxidase complex


ATP synthase.


  1. 1.

    Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio V, Wallace DC: Effects of Purifying and Adaptive Selection on Regional Variation in Human mtDNA. Science. 2004, 303 (5655): 223-226.

    PubMed  Article  Google Scholar 

  2. 2.

    Moyer GR, Winemeller KO, McPhee MV, Turner TF: Historical demography, selection, and coalescence of mitochondrial and nuclear genes in Prochilodus species of northern South America. Evolution. 2005, 59 (3): 599-610.

    PubMed  Article  Google Scholar 

  3. 3.

    Bazin E, Glemin S, Galtier N: Population Size Does Not Influence Mitochondrial Genetic Diversity in Animals. Science. 2006, 312 (5773): 570-572.

    PubMed  Article  Google Scholar 

  4. 4.

    Moreno-Loshuertos R, Acin-Perez R, Fernandez-Silva P, Movilla N, Perez-Martos A, Rodriguez de Cord S, Gallardo ME, Enriquez JA: Differences in reactive oxygen species production explain the phenotypes associated with common mouse mitochondrial DNA variants. Nat Genet. 2006, 38 (11): 1261-1268.

    PubMed  Article  Google Scholar 

  5. 5.

    Wallace DC: Mitochondrial Diseases in Man and Mouse. Science. 1999, 283 (5407): 1482-1488.

    PubMed  Article  Google Scholar 

  6. 6.

    Dalziel AC, Moyes CD, Fredriksson E, Lougheed SC: Molecular evolution of cytochrome c oxidase in high-performance fish (teleostei: Scombroidei). J Mol Evol. 2006, 62 (3): 319-331.

    PubMed  Article  Google Scholar 

  7. 7.

    Jobson RW, Nielsen R, Laakkonen L, Wikstrom M, Albert VA: Adaptive evolution of cytochrome c oxidase: Infrastructure for a carnivorous plant radiation. PNAS. 2004, 101 (52): 18064-18068.

    PubMed  Article  Google Scholar 

  8. 8.

    Grossman LI, Wildman DE, Schmidt TR, Goodman M: Accelerated evolution of the electron transport chain in anthropoid primates. Trends Genet. 2004, 20 (11): 578-585.

    PubMed  Article  Google Scholar 

  9. 9.

    Blier PU, Dufresne F, Burton RS: Natural selection and the evolution of mtDNA-encoded peptides: evidence for intergenomic co-adaptation. Trends Genet. 2001, 17 (7): 400-406.

    PubMed  Article  Google Scholar 

  10. 10.

    Ballard JW, Whitlock MC: The incomplete natural history of mitochondria. Mol Ecol. 2004, 13 (4): 729-744.

    PubMed  Article  Google Scholar 

  11. 11.

    Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, Brandon M, Easley K, Chen E, Brown MD, Sukernik RI, Olckers A, Wallace DC: Natural selection shaped regional mtDNA variation in humans. PNAS. 2003, 100 (1): 171-176.

    PubMed  Article  Google Scholar 

  12. 12.

    Rankinen T, Bray MS, Hagberg JM, Perusse L, Roth SM, Wolfarth B, Bouchard C: The human gene map for performance and health-related fitness phenotypes: the 2005 update. Med Sci Sports Exerc. 2006, 38 (11): 1863-1888.

    PubMed  Article  Google Scholar 

  13. 13.

    Suarez RK, Darveau CA, Childress JJ: Metabolic scaling: a many-splendoured thing. Comp Biochem Physiol B Biochem Mol Biol. 2004, 139 (3): 531-541.

    PubMed  Article  Google Scholar 

  14. 14.

    McNab BK: Energy constraints on carnivore diet. Nature. 2000, 407 (6804): 584-

    PubMed  Article  Google Scholar 

  15. 15.

    Arends A, McNab BK: The comparative energetics of 'caviomorph' rodents. Comp Biochem Physiol A Mol Integr Physiol. 2001, 130 (1): 105-122.

    PubMed  Article  Google Scholar 

  16. 16.

    Bonaccorso FJ, McNab BK: Standard energetics of leaf-nosed bats (Hipposideridae): its relationship to intermittent- and protracted-foraging tactics in bats and birds. J Comp Physiol [B]. 2003, 173 (1): 43-53.

    Google Scholar 

  17. 17.

    Rand DM, Haney RA, Fry AJ: Cytonuclear coevolution: the genomics of cooperation. Trends Ecol Evol. 2004, 19 (12): 645-653.

    PubMed  Article  Google Scholar 

  18. 18.

    Saraste M: Oxidative Phosphorylation at the fin de siècle. Science. 1999, 283 (5407): 1488-1493.

    PubMed  Article  Google Scholar 

  19. 19.

    Lin YH, Waddell PJ, Penny D: Pika and vole mitochondrial genomes increase support for both rodent monophyly and glires. Gene. 2002, 294 (1-2): 119-129.

    PubMed  Article  Google Scholar 

  20. 20.

    Reyes A, Gissi C, Catzeflis F, Nevo E, Pesole G, Saccone C: Congruent Mammalian Trees from Mitochondrial and Nuclear Genes Using Bayesian Methods. Mol Biol Evol. 2004, 21 (2): 397-403.

    PubMed  Article  Google Scholar 

  21. 21.

    Phillips MJ, Penny D: The root of the mammalian tree inferred from whole mitochondrial genomes. Mol Phylogenet Evol. 2003, 28 (2): 171-185.

    PubMed  Article  Google Scholar 

  22. 22.

    Arnason U, Adegoke JA, Bodin K, Born EW, Esa YB, Gullberg A, Nilsson M, Short RV, Xu X, Janke A: Mammalian mitogenomic relationships and the root of the eutherian tree. PNAS. 2002, 99 (12): 8151-8156.

    PubMed  Article  Google Scholar 

  23. 23.

    Murata Y, Nikaido M, Sasaki T, Cao Y, Fukumoto Y, Hasegawa M, Okada N: Afrotherian phylogeny as inferred from complete mitochondrial genomes. Mol Phylogenet Evol. 2003, 28 (2): 253-260.

    PubMed  Article  Google Scholar 

  24. 24.

    Broughton RE, Milam JE, Roe BA: The Complete Sequence of the Zebrafish (Danio rerio) Mitochondrial Genome and Evolutionary Patterns in Vertebrate Mitochondrial DNA. Genome Res. 2001, 11 (11): 1958-1967.

    PubMed  Google Scholar 

  25. 25.

    Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP: Bayesian Inference of Phylogeny and Its Impact on Evolutionary Biology. Science. 2001, 294 (5550): 2310-2314.

    PubMed  Article  Google Scholar 

  26. 26.

    Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574.

    PubMed  Article  Google Scholar 

  27. 27.

    Murphy WJ, Eizirik E, O'Brien SJ, Madsen O, Scally M, Douady CJ, Teeling E, Ryder OA, Stanhope MJ, de Jong WW, Springer MS: Resolution of the Early Placental Mammal Radiation Using Bayesian Phylogenetics. Science. 2001, 294 (5550): 2348-2351.

    PubMed  Article  Google Scholar 

  28. 28.

    Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA: TreeSAAP: Selection on Amino Acid Properties using phylogenetic trees. Bioinformatics. 2003, 19 (5): 671-672.

    PubMed  Article  Google Scholar 

  29. 29.

    Lopez JV, Culver M, Stephens JC, Johnson WE, O'Brien SJ: Rates of nuclear and cytoplasmic mitochondrial DNA sequence divergence in mammals. Mol Biol Evol. 1997, 14 (3): 277-286.

    PubMed  Article  Google Scholar 

  30. 30.

    Radermacher M, Ruiz T, Clason T, Benjamin S, Brandt U, Zickermann V: The three-dimensional structure of complex I from Yarrowia lipolytica: a highly dynamic enzyme. J Struct Biol. 2006, 154 (3): 269-279.

    PubMed  Article  Google Scholar 

  31. 31.

    Sazanov LA, Walker JE: Cryo-electron crystallography of two sub-complexes of bovine complex I reveals the relationship between the membrane and peripheral arms. J Mol Biol. 2000, 302 (2): 455-464.

    PubMed  Article  Google Scholar 

  32. 32.

    Brandt U: Energy converting NADH:quinone oxidoreductase (complex I). Annu Rev Biochem. 2006, 75: 69-92.

    PubMed  Article  Google Scholar 

  33. 33.

    Mathiesen C, Hagerhall C: Transmembrane topology of the NuoL, M and N subunits of NADH:quinone oxidoreductase and their homologues among membrane-bound hydrogenases and bona fide antiporters. Biochim Biophys Acta. 2002, 1556 (2-3): 121-132.

    PubMed  Article  Google Scholar 

  34. 34.

    Trumpower BL: The protonmotive Q cycle. Energy transduction by coupling of proton translocation to electron transfer by the cytochrome bc1 complex. J Biol Chem. 1990, 265 (20): 11409-11412.

    PubMed  Google Scholar 

  35. 35.

    Mitchell P: Possible molecular mechanisms of the protonmotive function of cytochrome systems. J Theor Biol. 1976, 62 (2): 327-367.

    PubMed  Article  Google Scholar 

  36. 36.

    Rohland N, Malaspinas AS, Pollack JL, Slatkin M, Matheus P, Hofreiter M: Proboscidean Mitogenomics: Chronology and Mode of Elephant Evolution Using Mastodon as Outgroup. PLoS Biol. 2007, 5 (8): e207-

    PubMed  Article  Google Scholar 

  37. 37.

    Esser L, Gong X, Yang S, Yu L, Yu CA, Xia D: Surface-modulated motion switch: Capture and release of iron-sulfur protein in the cytochrome bc1 complex. PNAS. 2006, 103 (35): 13045-13050.

    PubMed  Article  Google Scholar 

  38. 38.

    Castellini MA, Castellini JM: Defining the limits of diving biochemistry in marine mammals. Comp Biochem Physiol B Biochem Mol Biol. 2004, 139 (3): 509-518.

    PubMed  Article  Google Scholar 

  39. 39.

    Darrouzet E, Valkova-Valchanova M, Moser CC, Dutton PL, Daldal F: Uncovering the [2Fe2S] domain movement in cytochrome bc1 and its implications for energy conversion. PNAS. 2000, 97 (9): 4567-4572.

    PubMed  Article  Google Scholar 

  40. 40.

    Shoshani J: Understanding proboscidean evolution: a formidable task. Trends in Ecology & Evolution. 1998, 13 (12): 480-487.

    Article  Google Scholar 

  41. 41.

    Hochachka PW, Stanley C, Merkt J, Sumar-Kalinowski J: Metabolic meaning of elevated levels of oxidative enxymes in high altitude adapted animals: an interpretive hypotheses. Respiration Physiology. 1982, 52: 303-313.

    Article  Google Scholar 

  42. 42.

    Belevich I, Verkhovsky MI, Wikstrom M: Proton-coupled electron transfer drives the proton pump of cytochrome c oxidase. Nature. 2006, 440 (7085): 829-832.

    PubMed  Article  Google Scholar 

  43. 43.

    Overholtzer MH, Yakowec PS, Cameron V: The Effect of Amino Acid Substitutions in the Conserved Aromatic Region of Subunit II of Cytochrome c Oxidase in Saccharomyces cerevisiae. J Biol Chem. 1996, 271 (13): 7719-7724.

    PubMed  Article  Google Scholar 

  44. 44.

    Hofacker I, Schulten K: Oxygen and proton pathways in cytochrome c oxidase. Proteins. 1998, 30 (1): 100-107.

    PubMed  Article  Google Scholar 

  45. 45.

    Richter OM, Ludwig B: Cytochrome c oxidase--structure, function, and physiology of a redox-driven molecular machine. Rev Physiol Biochem Pharmacol. 2003, 147: 47-74.

    PubMed  Google Scholar 

  46. 46.

    Schmidt B, McCracken J, Ferguson-Miller S: A discrete water exit pathway in the membrane protein cytochrome c oxidase. PNAS. 2003, 100 (26): 15539-15542.

    PubMed  Article  Google Scholar 

  47. 47.

    Fillingame RH, Angevine CM, Dmitriev OY: Mechanics of coupling proton movements to c-ring rotation in ATP synthase. FEBS Lett. 2003, 555 (1): 29-34.

    PubMed  Article  Google Scholar 

  48. 48.

    Leyva JA, Bianchet MA, Amzel LM: Understanding ATP synthesis: structure and mechanism of the F1-ATPase (Review). Mol Membr Biol. 2003, 20 (1): 27-33.

    PubMed  Article  Google Scholar 

  49. 49.

    Mueller DM: Partial assembly of the yeast mitochondrial ATP synthase. J Bioenerg Biomembr. 2000, 32 (4): 391-400.

    PubMed  Article  Google Scholar 

  50. 50.

    Hadikusumo RG, Meltzer S, Choo WM, Jean-Francois MJ, Linnane AW, Marzuki S: The definition of mitochondrial H+ ATPase assembly defects in mit- mutants of Saccharomyces cerevisiae with a monoclonal antibody to the enzyme complex as an assembly probe. Biochim Biophys Acta. 1988, 933 (1): 212-222.

    PubMed  Article  Google Scholar 

  51. 51.

    Weibel ER: Physiology: the pitfalls of power laws. Nature. 2002, 417 (6885): 131-132.

    PubMed  Article  Google Scholar 

  52. 52.

    Darveau CA, Suarez RK, Andrews RD, Hochachka PW: Allometric cascade as a unifying principle of body mass effects on metabolism. Nature. 2002, 417 (6885): 166-170.

    PubMed  Article  Google Scholar 

  53. 53.

    White CR, Seymour RS: Mammalian basal metabolic rate is proportional to body mass2/3. PNAS. 2003, 100 (7): 4046-4049.

    PubMed  Article  Google Scholar 

  54. 54.

    West GB, Woodruff WH, Brown JH: Allometric scaling of metabolic rate from molecules and mitochondria to cells and mammals. PNAS. 2002, 99 (90001): 2473-2478.

    PubMed  Article  Google Scholar 

  55. 55.

    Else PL, Hulbert AJ: An allometric comparison of the mitochondria of mammalian and reptilian tissues: the implications for the evolution of endothermy. J Comp Physiol [B]. 1985, 156 (1): 3-11.

    Article  Google Scholar 

  56. 56.

    Wang ZM, O'Connor TP, Heshka S, Heymsfield SB: The Reconstruction of Kleiber's Law at the Organ-Tissue Level. J Nutr. 2001, 131 (11): 2967-2970.

    PubMed  Google Scholar 

  57. 57.

    West GB, Savage VM, Gillooly J, Enquist BJ, Woodruff WH, Brown JH: Physiology: Why does metabolic rate scale with body size?. Nature. 2003, 421 (6924): 713; discussion 714-

    PubMed  Article  Google Scholar 

  58. 58.

    Gibson A, Gowri-Shankar V, Higgs PG, Rattray M: A Comprehensive Analysis of Mammalian Mitochondrial Genome Base Composition and Improved Phylogenetic Methods. Mol Biol Evol. 2005, 22 (2): 251-264.

    PubMed  Article  Google Scholar 

  59. 59.

    Saccone C, Lanave C, Pesole G, Preparata G: Influence of base composition on quantitative estimates of gene evolution. Methods Enzymol. 1990, 183: 570-583.

    PubMed  Article  Google Scholar 

  60. 60.

    Nylander JAA: MrModeltest v2. 2004, Evolutionary Biology Centre, Uppsala University

    Google Scholar 

  61. 61.

    Swofford DL: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 1998, Sunderland, Massachussetts , Sinauer Associates

    Google Scholar 

  62. 62.

    Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. 1998, 14 (9): 817-818.

    Google Scholar 

  63. 63.

    Yang Z: Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol. 1998, 15 (5): 568-573.

    PubMed  Article  Google Scholar 

  64. 64.

    McClellan DA, Palfreyman EJ, Smith MJ, Moss JL, Christensen RG, Sailsbery JK: Physicochemical Evolution and Molecular Adaptation of the Cetacean and Artiodactyl Cytochrome b Proteins. Mol Biol Evol. 2005, 22 (3): 437-455.

    PubMed  Article  Google Scholar 

  65. 65.

    McClellan DA, McCracken KG: Estimating the Influence of Selection on the Variable Amino Acid Sites of the Cytochrome b Protein Functional Domains. 2001, 18 (6): 917-925.

    Google Scholar 

  66. 66.

    Huang LS, Cobessi D, Tung EY, Berry EA: Binding of the respiratory chain inhibitor antimycin to the mitochondrial bc1 complex: a new crystal structure reveals an altered intramolecular hydrogen-bonding pattern. J Mol Biol. 2005, 351 (3): 573-597.

    PubMed  Article  Google Scholar 

  67. 67.

    Tsukihara T, Shimokata K, Katayama Y, Shimada H, Muramoto K, Aoyama H, Mochizuki M, Shinzawa-Itoh K, Yamashita E, Yao M, Ishimura Y, Yoshikawa S: The low-spin heme of cytochrome c oxidase as the driving element of the proton-pumping process. PNAS. 2003, 100 (26): 15304-15309.

    PubMed  Article  Google Scholar 

  68. 68.

    Krogh A, Larsson B, von Heijne G, Sonnhammer EL: Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001, 305 (3): 567-580.

    PubMed  Article  Google Scholar 

  69. 69.

    Sonnhammer EL, von Heijne G, Krogh A: A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998, 6: 175-182.

    PubMed  Google Scholar 

  70. 70.

    Viklund H, Elofsson A: Best {alpha}-helical transmembrane protein topology predictions are achieved using hidden Markov models and evolutionary information. Protein Sci. 2004, 13 (7): 1908-1917.

    PubMed  Article  Google Scholar 

  71. 71.

    Tusnady GE, Simon I: Principles governing amino acid composition of integral membrane proteins: application to topology prediction. J Mol Biol. 1998, 283 (2): 489-506.

    PubMed  Article  Google Scholar 

  72. 72.

    Saaf A, Johansson M, Wallin E, von Heijne G: Divergent evolution of membrane protein topology: The Escherichia coli RnfA and RnfE homologues. PNAS. 1999, 96 (15): 8540-8544.

    PubMed  Article  Google Scholar 

  73. 73.

    Nilsson J, Persson B, von Heijne G: Consensus predictions of membrane protein topology. FEBS Lett. 2000, 486 (3): 267-269.

    PubMed  Article  Google Scholar 

  74. 74.

    Humphrey W, Dalke A, Schulten K: VMD: visual molecular dynamics. J Mol Graph. 1996, 14 (1): 33-8, 27-8.

    PubMed  Article  Google Scholar 

  75. 75.

    Hinttala R, Smeets R, Moilanen JS, Ugalde C, Uusimaa J, Smeitink JAM, Majamaa K: Analysis of mitochondrial DNA sequences in patients with isolated or combined oxidative phosphorylation system deficiency. J Med Genet. 2006, 43 (11): 881-886.

    PubMed  Article  Google Scholar 

  76. 76.

    Simon DK, Friedman J, Breakefield XO, Jankovic J, Brin MF, Provias J, Bressman SB, Charness ME, Tarsy D, Johns DR, Tarnopolsky MA: A heteroplasmic mitochondrial complex I gene mutation in adult-onset dystonia. Neurogenetics. 2003, 4 (4): 199-205.

    PubMed  Article  Google Scholar 

  77. 77.

    Andreu AL, Tanji K, Bruno C, Hadjigeorgiou GM, Sue CM, Jay C, Ohnishi T, Shanske S, Bonilla E, DiMauro S: Exercise intolerance due to a nonsense mutation in the mtDNA ND4 gene. Ann Neurol. 1999, 45 (6): 820-823.

    PubMed  Article  Google Scholar 

  78. 78.

    Karadimas CL, Greenstein P, Sue CM, Joseph JT, Tanji K, Haller RG, Taivassalo T, Davidson MM, Shanske S, Bonilla E, DiMauro S: Recurrent myoglobinuria due to a nonsense mutation in the COX I gene of mitochondrial DNA. Neurology. 2000, 55 (5): 644-649.

    PubMed  Article  Google Scholar 

  79. 79.

    Rahman S, Taanman JW, Cooper JM, Nelson I, Hargreaves I, Meunier B, Hanna MG, Garcia JJ, Capaldi RA, Lake BD, Leonard JV, Schapira AH: A missense mutation of cytochrome oxidase subunit II causes defective assembly and myopathy. Am J Hum Genet. 1999, 65 (4): 1030-1039.

    PubMed  Article  Google Scholar 

  80. 80.

    Hanna MG, Nelson IP, Rahman S, Lane RJ, Land J, Heales S, Cooper MJ, Schapira AH, Morgan-Hughes JA, Wood NW: Cytochrome c oxidase deficiency associated with the first stop-codon point mutation in human mtDNA. Am J Hum Genet. 1998, 63 (1): 29-36.

    PubMed  Article  Google Scholar 

  81. 81.

    Horvath R, Scharfe C, Hoeltzenbein M, Do BH, Schroder C, Warzok R, Vogelgesang S, Lochmuller H, Muller-Hocker J, Gerbitz KD, Oefner PJ, Jaksch M: Childhood onset mitochondrial myopathy and lactic acidosis caused by a stop mutation in the mitochondrial cytochrome c oxidase III gene. J Med Genet. 2002, 39 (11): 812-816.

    PubMed  Article  Google Scholar 

  82. 82.

    Dumoulin R, Sagnol I, Ferlin T, Bozon D, Stepien G, Mousson B: A novel gly290asp mitochondrial cytochrome b mutation linked to a complex III deficiency in progressive exercise intolerance. Mol Cell Probes. 1996, 10 (5): 389-391.

    PubMed  Article  Google Scholar 

  83. 83.

    Andreu AL, Hanna MG, Reichmann H, Bruno C, Penn AS, Tanji K, Pallotti F, Iwata S, Bonilla E, Lach B, Morgan-Hughes J, DiMauro S, Shanske S, Sue CM, Pulkes T, Siddiqui A, Clark JB, Land J, Iwata M, Schaefer J: Exercise Intolerance Due to Mutations in the Cytochrome b Gene of Mitochondrial DNA. N Engl J Med. 1999, 341 (14): 1037-1044.

    PubMed  Article  Google Scholar 

  84. 84.

    Legros F, Chatzoglou E, Frachon P, Ogier De Baulny H, Laforet P, Jardel C, Godinot C, Lombes A: Functional characterization of novel mutations in the human cytochrome b gene. Eur J Hum Genet. 2001, 9 (7): 510-518.

    PubMed  Article  Google Scholar 

  85. 85.

    Andreu AL, Checcarelli N, Iwata S, Shanske S, DiMauro S: A missense mutation in the mitochondrial cytochrome b gene in a revisited case with histiocytoid cardiomyopathy. Pediatr Res. 2000, 48 (3): 311-314.

    PubMed  Article  Google Scholar 

  86. 86.

    Schuelke M, Krude H, Finckh B, Mayatepek E, Janssen A, Schmelz M, Trefz F, Trijbels F, Smeitink J: Septo-optic dysplasia associated with a new mitochondrial cytochrome b mutation. Ann Neurol. 2002, 51 (3): 388-392.

    PubMed  Article  Google Scholar 

  87. 87.

    AnAge Database. []

Download references


This work was supported in part by the Project POCTI/BSE/47559/2002 and PTDC/BIA-BDE/69144/2006 from the Portuguese Foundation for Science and Technology (FCT) and by the National Geographic Society Grant 7483-03. RF is funded by FCT fellowship SFRH/BPD/26769/2006. This project has been funded in part with federal funds from the National Cancer Institute, National Institutes of Health, under contract N01-CO-12400. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government. Comments made by three anonymous referees improved a previous version of this manuscript.

Author information



Corresponding author

Correspondence to Agostinho Antunes.

Additional information

Authors' contributions

RF performed all phylogenetic, evolutionary and structure-function analyses and drafted the manuscript.

WEJ participated in the drafting and coordination of the study.

SJB participated in the drafting and coordination of the study.

MJR participated in the drafting and coordination of the study.

AA participated in the design, genetic analyses, drafting and coordination of the study.

All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Supporting information. Saturation plots (S1), maximum likelihood tree (S2), number of strong positively selected amino acid properties plotted against the branch length (S3), number of strong positively selected amino acid properties across all sites in each mammalian mtDNA-encoded subunit of the respiratory chain (S4) and the sites indicated in Figure 10 in four elephant relatives and 6 alpaca relatives (S5). (PDF 2 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

da Fonseca, R.R., Johnson, W.E., O'Brien, S.J. et al. The adaptive evolution of the mammalian mitochondrial genome. BMC Genomics 9, 119 (2008).

Download citation


  • Humpback Whale
  • Amino Acid Property
  • Interior Branch
  • Mitochondrial Protein Code Gene
  • Local Secondary Structure