the oxidative phosphorylation ( oxphos ) system is unique among fundamental metabolic pathways in animals ( aside from maintenance and gene expression of mitochondria ) in that functional oxphos complexes are composed of subunits encoded by two different genomes .
the different hereditary modes of the mitochondrial ( mt ) and nuclear ( nu ) genomes substantially increase the complexity of maintaining functional interactions among oxphos subunits .
the role of oxphos in the fundamental process of aerobic atp production means that slight perturbations to inter - genomic coordination result in a wide range of human pathologies ( wallace 2010 ) . in addition , minimal divergence of oxphos subunits may result in reproductive incompatibility among closely related populations ( rand et al .
consequently , the successful interaction of oxphos proteins encoded by mt and nu genomes plays a central role in fundamental processes at the cellular , organismal , and population levels of biological organization .
the oxphos system provides a rare opportunity to investigate the coevolution of subunits that must successfully interact in electron transport and atp production , yet are subject to potentially the very different evolutionary forces that govern the mt and nu genomes .
it produces about 80% of energy in the form of atp in almost all eukaryote cells and is the major role - player making mitochondria the  powerhouse of the cell .
 meantime , reactive oxygen species , the by - product of normal metabolism , can damage dna , cell membranes and lipid , decrease bioenergetic efficiency , and thus lead to age , disease , and death in humans ( ballard and melvin 2010 ; wallace 2010 ) .
the oxphos apparatus includes five complexes ( i  v ) , each composed of between 4 and 30 + subunits ( de grassi et al .
complexes i  iv transfer electrons through a series of redox reactions coupled to the transport of protons into the intermembrane space .
complex v employs the resulting electrochemical gradient to phosphorylate adp to atp . like most multimeric protein complexes
an added level of complexity arises from the fact that the composition of four of the five oxphos complexes includes subunits encoded by both mt and nu genomes .
any mismatch between mt and nu genomes likely results in a slowing of electron transfer , which increases reactivity of electron with oxygen , corresponding to a rise in free - radical leak ( lane 2011 ) .
evolution of oxphos must therefore be coordinated among a large number of interacting subunits and among genomes , and compensatory amino acid substitutions may , in some cases , be favored to maintain function .
complex i ( nadh : quinone oxidoreductase ) and complex ii ( succinate dehydrogenase ) receive electrons from reducing equivalents ( nadh and fadh2 ) produced by the krebs cycle .
the electrons are transferred to ubiquinone , which freely diffuses within the mt inner membrane and transfers electrons to complex iii ( cytochrome bc1 ) .
complex iii transfers electrons from ubiquinol to cytochrome c. this small protein is localized on the intermembrane space side of the mt inner membrane .
it transfers electrons to complex iv ( cytochrome c oxidase ) , where oxygen is reduced to water .
coupled with the process of electron transport , complexes i , iii , and iv transport protons to the intermembrane space .
the generated proton gradient across the mt inner membrane drives atp synthesis in complex v ( atp synthase ) ( mitchell 1961 ; saraste 1999 ; boekema and braun 2007 ) . the mt genome has a constant gene content composition in all metazoan species ( with few exceptions ) , including 13 protein - coding genes , 2 rrna genes , and 22 trna genes .
the rrna and trna genes are involved in mt protein synthesis , whereas the 13 protein - coding genes all encode components of oxphos ( shadel and clayton 1997 ; gray et al . 1999 ) .
at least 60 oxphos subunits are encoded in the nu genome . despite the fundamental importance of oxphos to aerobic life ,
mt protein - coding genes are not highly conserved and evolve at rates that is 550 times that of typical nu genes in vertebrates ( lynch 2007 ) .
this increased rate has been observed in a range of animal species , including primates ( brown et al .
1979 ) , galapagos tortoises ( caccone et al . 2004 ) , lice ( johnson et al . 2003 ) , drosophila ( haag - liautard et al . 2008 ) , wasps ( kaltenpoth et al . 2012 ) , and nematodes ( denver et al . 2000 )
. the rapid evolutionary rate of mt genomes may be due to their cell cycle - independent replication of mtdna ( bogenhagen and clayton 1977 ) , increased exposure to mutagenic oxygen radical species ( beal 1996 ) , lack of protective histones , and limited dna repair capacity ( croteau and bohr 1997 ) .
moreover , the mt genome differs from the nu genome in being effectively haploid , maternally inherited , and exhibiting little recombination ( moritz et al . 1987 ; lightowlers et al . 1997 ; ballard and melvin 2010 ) .
haploid maternal inheritance means that the effective population size for mt genes is approximately one - fourth that of nu autosomal genes , increasing the fixation rate of polymorphic sites ( gabriel et al .
in addition , the lack of recombination in animal mitochondria reduces the ability to purge deleterious mutations ( gabriel et al .
incompatibility between mt and nu genomes has been shown to reduce hybrid fitness in copepods ( willett and burton 2001 ) , drosophila ( sackton et al .
nu incompatibility may be an important contributor to the process of speciation ( gershoni et al . 2009 ) and the evolution of two separate sexes ( hadjivasiliou et al .
, amino acid substitutions may be deleterious if they affect domains involved in interaction , yet there is the potential for compensatory changes at interacting amino acid sites to reduce or eliminate any negative effects .
such interacting sites are generally proximal to each other in the three - dimensional protein structure ( pazos and valencia 2008 ) . under the compensatory model ,
substitutions that may destabilize protein structures or inhibit function could be compensated by a subsequent ( or simultaneous ) substitution at an adjacent site ( pollock et al .
, the elevated evolutionary rate of mtdna suggests that deleterious substitutions occur more frequently in mt genes and compensatory changes may then occur in proximal sites of nu - encoded proteins . in bacteria ( sharp and li 1987 ) , drosophila ( comeron and kreitman 1998 ; dunn et al .
2001 ) , and mammals ( wolfe and sharp 1993 ) , the rate of nonsynonymous substitutions ( dn ) is positively correlated with the synonymous substitution rate ( ds ) and both dn and ds are higher in mt than nu genes .
if nonsynonymous substitutions of mt genes drive corresponding nonsynonymous substitutions of nu oxphos genes , we expect to see the acceleration of dn in nu oxphos genes relative to most nu genes not involved in oxphos ( nu non - oxphos ) .
nu coevolution at the molecular level but have been restricted to a few genes or a small number of species ( e.g. , nu genes cyc1 and uqcrfs1 , and mt gene mt - cytb [ willett and burton 2001 ] , or genes involved in oxphos complex iv [ goldberg et al . 2003 ] ) . here
, we describe a broader study of 73 oxphos genes and a comparison between oxphos and 77 non - oxphos housekeeping genes . under a compensatory evolution scenario
, we predicted that dn of mt oxphos > nu oxphos > nu non - oxphos genes .
we tested the compensatory substitution model by comparing the evolutionary rates ( both ds and dn ) of 13 mt protein - coding genes ( mt oxphos genes ) , 60 nu oxphos genes , and 77 non - oxphos genes in 47 vertebrate species , including 7 fishes and 40 mammals .
the seven fish species are phylogenetically disparate , spanning some 250 myr of evolutionary history ( betancur - r et al .
1a ) , whereas most of the mammal lineages emerged within the last 80 myr ( meredith et al .
this study employed genome sequences of 7 teleost fishes and 40 mammals from which we acquired the coding regions of 13 mt oxphos and 60 nu oxphos genes . for comparison
comparators as oxphos genes contribute to critical function and should be highly expressed in virtually all cell types .
these were chosen randomly with respect to genome location and represent a wide range of functions .
they should be broadly representative of functionally important and relatively highly expressed loci in the nu genome .
all sequences were downloaded from the ensembl genome browser ( www.ensembl.org , last accessed september 11 , 2013 ) using a pipeline procedure ( vilella et al .
a phylogenetic gene tree was used to identify orthologous and paralogous sequences so that we used only the former . each gene sequence
a few genes were unavailable in the genome sequences ; however , at least 87% of the genes were present in 33 species
. the majority of the missing genes were not core oxphos subunits whereas the specific genes that were absent varied among species .
1a ) , whereas the phylogeny of mammals was from meredith et al . ( 2011 ) ( fig .
the program codeml implemented in the computer software paml v.4.7 ( yang 2007 ) was used to estimate dn and ds for each gene .
the site - model m0 ( model : 0 , nssites : 0 , and fix - omega : 0 ) was used when estimating dn and ds . for each gene ,
the estimated substitution rates are the sum of dn or ds values from each branch of the tree , so total rate values will vary with the number of taxa in the tree .
the intent here was not to compare the absolute values between fishes and mammals , but instead to examine the patterns among mt oxphos genes , nu oxphos genes , and non - oxphos genes within taxonomic groups ( where the number of branches was identical ) .
codon frequencies were estimated from the average nucleotide frequencies at three codon positions . to avoid problems arising from local optima ( yang and nielsen 1998 ) , replicate runs with three different starting values of  ( 0.3 , 0.9 , and 4.3 ) were used .
different starting  yielded similar , and in many cases identical , ds and dn estimates for each gene among three runs ( supplementary table s2 , supplementary material online ) .
the average ds and dn from the three runs was used in the subsequent analysis .
we tested the hypothesis that there are differences in ds and dn among three groups of genes : 13 mt oxphos genes , 60 nu oxphos genes , and 77 non - oxphos genes .
a shapiro  wilks test was performed to assess the normality of ds and dn for all genes and bartlett s test was performed to evaluate homogeneity of variance .
wallis test followed by post hoc tukey s hsd tests were used to assess differences among the groups .
a wilcoxon signed - rank test was performed to evaluate differences between mt and nu oxphos genes in each of the five oxphos complexes .
each of the oxphos complexes is composed of multiple subunits and each subunit may be more or less important to the core catalytic activity of the complex .
noncore ,  where core genes encode subunits that retain structural , functional , and sequence homology with prokaryotic systems .
for example , unlike human complex i , which consists of 43 subunits , bacterial complex i consists of 14 subunits .
these 14 subunits all have homologs in human complex i , and thus are referred as core subunits .
these include seven nu oxphos genes ( ndufs13 , ndufs78 , and ndufv12 ) and seven mt oxophos genes ( mt - nd16 , mt - nd4l ) ( janssen et al .
similarly , in complex iv , bovine mt oxphos genes ( mt - cox13 ) share high similarity to their bacterial counterparts and can be considered the functional core of the eukaryotic oxidase ( richter and ludwig 2003 ) .
complex ii is only composed of four nu subunits ( sdha - d ) , all of which share high sequence similarity with escherichia coli ( cecchini 2003 ) .
in complex v , six human subunits have bacterial homologs ( atp5a1 , atp5b , atp5c1 , atp5d , atp5f1 , and atp5g3 ) ( yoshida et al . 2001 ) and were categorized as core subunits . because e. coli does not have a homolog of complex iii ( lenn et al .
2008 ) , none of the nu oxphos genes in complex iii were categorized as encoding core subunits .
we estimated the total synonymous and nonsynonymous substitution rates on each tree for fishes ( fig .
1b ) separately using maximum likelihood methods . because neither ds nor dn was normally distributed and each exhibited significant variance heterogeneity , nonparametric statistics were used to compare the different groups of genes .
we found significant differences ( kruskal  wallis ) in ds and dn among the three groups of genes : 13 mt oxphos genes , 60 nu oxphos genes , and 77 non - oxphos genes ( table 1 ) .
the ds of mt oxphos genes was significantly higher than that of nu oxphos and non - oxphos genes in both ( fig .
the ds of mt and nu genes were 16.02  9.64 and 2.02  0.87 , respectively , in fishes and 72.27  26.56 and 8.85  4.15 , respectively , in mammals .
no significant difference of the ds was detected between nu oxphos genes and non - oxphos genes ( fig .
we note that a small number of mt genes had exceptionally high ds values , which may indicate substitutional saturation , leading to violation of the assumptions of the maximum likelihood estimation method .
as predicted , the dn of mt oxphos genes was significantly higher than that of nu oxphos genes , which in turn was significantly higher than that of non - oxphos genes ( fig .
2.comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear oxphos ( nu oxphos ) , and non - oxphos genes in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
whisker - ends are at the 5th and 95th percentiles . * p < 0.05 , * * p < 0.01 , * * * p < 0.001 . 
 table 1statistical summary for the comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear oxphos ( nu oxphos ) , and non - oxphos genes in 7 fishes and 40 mammalsml estimatestaxacomparisonp valuedsfishesmt oxphos > nu oxphos<0.0001dsmammalsmt oxphos > nu oxphos<0.0001dsfishesmt oxphos > non - oxphos<0.0001dsmammalsmt oxphos > non - oxphos<0.0001dsfishesnu oxphos > non - oxphos0.987dsmammalsnu oxphos > non - oxphos0.670dnfishesmt oxphos >
nu oxphos0.037dnmammalsmt oxphos > nu oxphos<0.0001dnfishesmt oxphos > non - oxphos<0.0001dnmammalsmt oxphos > non - oxphos<0.0001dnfishesnu oxphos > non - oxphos0.0007dnmammalsnu oxphos > non - oxphos<0.0001 comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear oxphos ( nu oxphos ) , and non - oxphos genes in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
* p < 0.05 , * * p < 0.01 , * * * p < 0.001 . statistical summary for the comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear oxphos ( nu oxphos ) , and non - oxphos genes in 7 fishes and 40 mammals estimates of ds and dn were assessed for each oxphos complexes .
the ds of mt genes was higher than that of nu genes in all four complexes ( complex ii is composed of subunits only encoded by nu genes ) ( fig . 3a and c ) .
the dn of mt genes was not consistently higher than that of nu genes in all complexes .
for example , the dn of mt genes in complex iii was lower than that of nu genes in fishes ( fig .
3d ) ; the dn of mt genes in complex iv was significantly lower than that of nu genes in fishes ( fig .
3.comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) between mitochondrial ( mt ) and nuclear ( nu ) genes in each complex in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
whisker - ends are at the 5th and 95th percentiles . * p < 0.05 , * * p < 0.01 , * * * p < 0.001 .
comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) between mitochondrial ( mt ) and nuclear ( nu ) genes in each complex in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
whisker - ends are at the 5th and 95th percentiles . * p < 0.05 , * * p < 0.01 , * * * p < 0.001 . because of their different functions and interactions
, oxphos genes may be under different levels of functional constraint , so we separated them into core and noncore groups .
core subunits were defined as those that have bacterial homologs and that participate directly in electron transport or atp synthesis . after partitioning ,
ds of mt oxphos genes was found to be significantly higher than that of all nu gene groups ( fig .
4a and c ) , whereas the dn of noncore oxphos genes was higher than that of core oxphos genes ( fig .
there were no significant differences in dn between nu non - oxphos genes and nu core oxphos genes ( fig .
4.comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear core oxphos ( nu core oxphos ) , nuclear noncore oxphos ( nu noncore oxphos ) , and non - oxphos genes in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
whisker - ends are at the 5th and 95th percentiles . * p < 0.05 , * * p < 0.01 , * * * p < 0.001 . 
 table 2statistical summary for the comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear core oxphos ( nu core oxphos ) , nuclear noncore oxphos ( nu noncore oxphos ) , and non - oxphos genes in 7 fishes and 40 mammalsml estimatestaxacomparisonp valuedsfishesmt oxphos > nu core oxphos<0.0001dsmammalsmt oxphos > nu core oxphos<0.0001dsfishesmt oxphos > nu noncore oxphos<0.0001dsmammalsmt oxphos > nu noncore oxphos<0.0001dsfishesmt oxphos > non - oxphos<0.0001dsmammalsmt oxphos > non - oxphos<0.0001dsfishesnu core oxphos < nu noncore oxphos0.881dsmammalsnu core oxphos > nu noncore oxphos0.9995dsfishesnu core oxphos < non - oxphos0.966dsmammalsnu core oxphos < non - oxphos0.967dsfishesnu noncore oxphos > non - oxphos0.969dsmammalsnu noncore oxphos < non - oxphos0.847dnfishesmt oxphos > nu core oxphos0.002dnmammalsmt oxphos > nu core oxphos<0.0001dnfishesmt oxphos > nu noncore oxphos<0.304dnmammalsmt oxphos > nu noncore oxphos<0.0001dnfishesmt oxphos > non - oxphos<0.0001dnmammalsmt oxphos > non - oxphos<0.0001dnfishesnu core oxphos < nu noncore oxphos0.031dnmammalsnu core oxphos < nu noncore oxphos0.193dnfishesnu core oxphos > non - oxphos0.988dnmammalsnu core oxphos > non - oxphos0.198dnfishesnu noncore oxphos > non - oxphos<0.0001dnmammalsnu noncore oxphos > non - oxphos<0.0001 comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear core oxphos ( nu core oxphos ) , nuclear noncore oxphos ( nu noncore oxphos ) , and non - oxphos genes in 7 fishes ( a , b ) and 40 mammals ( c , d ) .
* p < 0.05 , * * p < 0.01 , * * * p < 0.001 . statistical summary for the comparison of synonymous substitution rate ( ds ) and nonsynonymous substitution rate ( dn ) among mitochondrial oxidative phosphorylation ( mt oxphos ) , nuclear core oxphos ( nu core oxphos ) , nuclear noncore oxphos ( nu noncore oxphos ) , and non - oxphos genes in 7 fishes and 40 mammals
given that mt genome is derived from once free - living bacteria , it is clear that most of the bacterial genes have now been removed to the cell nucleus .
however , what has not been well resolved is why not all genes have been removed to the nucleus ( alberts et al .
the remained genes in mitochondria thus are the resources for any potential problems involved in the interaction between mt 
, we used estimated rates of synonymous and nonsynonymous substitutions among all oxphos genes and a broad genomic sample of housekeeping genes . specifically , we examined the hypothesis of compensatory evolution between mt and nu oxphos genes by comparing their evolutionary rates ( both ds and dn ) relative to non - oxphos genes in two distinct vertebrate lineages .
we found that the mean ds of mt genes was about 79 times higher than that of nu genes .
this magnitude is consistent with the well - established pattern of high ds in animal mt genes ( brown et al .
although the products of nu oxhpos genes are transported into mitochondria where they function , they are expected to show mean ds similar to the non - oxphos genes because the coding sequences reside in the nu genome .
our results show this to be the case and suggest that differences in nonsynonymous rates are due to selective factors rather than different mutation rates . mt 
transplants  interacting partners should result in diminished functional performance , and this disruption should increase as the level of evolutionary divergence increases ( rand et al .
2004 ) . in mouse ( mus musculus domesticus ) cell lines carrying mitochondria from six different murid species spanning 212 myr of divergence ,
a near - linear association between disruption of respiratory chain function and evolutionary distance has been observed ( mckenzie et al .
this phenomenon has also been observed in primates ( kenyon and moraes 1997 ) , copepods ( willett and burton 2001 ) , drosophila ( sackton et al .
introgression experiments of drosophila simulans siii mtdna type into the sympatric population siiii nu background did not show a difference in catalytic properties of mitochondria , indicating that some naturally occurring mutations in mtdna can be accommodated by different nu background and mt 
nu interaction ( pichaud et al . 2012 ) . at the dna sequence level , higher evolutionary rates of mt genes could drive accelerated evolutionary rates of nu genes as compensatory response in the nu genes contributes to the maintenance of function .
several studies of the primate complex iv showed that seven nu oxphos genes : cox4i1 ( wu et al .
2002 ) , cox7a1h ( schmidt et al . 1999 ) , cox5a ( doan et al . 2004 ) , cox8al ( osada et al .
2003 ) , cox6b1 , cox6c , and cox7c ( doan et al . 2004 ) , together with two mt oxphos genes : mt - cox1 ( andrews and easteal 2000 ; wu et al .
2000 ) and mt - cox2 ( adkins and honeycutt 1994 ) , have shown accelerated dn in the lineage leading to hominids relative to other primates . however , in these cases it is not clear whether accelerated dn is due to compensatory changes driven by the higher mt rate , or is due to selection for increased oxphos efficiency associated with increased brain size in the hominid lineage .
this study showed that when all complexes were analyzed together , the dn of 13 mt oxphos genes was significantly higher than that of 60 nu oxphos genes , and both of these were significantly higher than dn of 77 non - oxphos genes .
the apparent acceleration dn of nu oxphos genes is consistent with a compensatory response to the higher dn of mt oxphos genes .
contrary to expectations , our results showed that the dn of mt genes was not always higher than that of nu genes in each complex .
mt genes with higher dn than nu genes were only observed in complexes i and v. in fishes , the dn of mt genes was lower than that of nu genes in complexes iii and iv . in mammals ,
the dn of mt genes was lower than that of nu genes in complex iv , but higher in complex iii .
( 2003 ) found that in the anthropoid lineage , the fast - evolving region ( 12 n - terminal amino acid residues ) of cox8l encodes amino acid sites contacted with mt - cox1 sites in three - dimensional structure .
this suggested structurally mediated cytonuclear coevolution , which was driven by faster evolving mt oxphos genes .
similarly , the only mt gene in complex iii ( mt - cytb ) showed lower dn than its corresponding nu oxphos genes in fishes .
although dn of mt - cytb was higher than nu genes in complex iii in mammals , the difference is not statistically different .
one explanation for the lower dn of mt oxphos genes in complexes iii and iv could be functional constraints on core subunits .
mt - cytb is the only mt gene in complex iii and its structure has been shown to be conserved among vertebrates ( kocher et al .
similarly , mt - cox1 and mt - cox2 have been recognized as the core enzyme catalytic subunits and are conserved from bacteria to bovid ( tsukihara et al .
it therefore seems likely that functional constraints are stronger on core subunits , regardless of which genome encodes them .
the higher dn of nu oxphos genes relative to non - oxphos genes in general is due mainly to higher rates of noncore proteins alone . the more recently derived noncore subunits ( or  accessory oxphos families  )
appear to be important contributors to oxphos assembly and/or stabilization ( ugalde et al . 2004 ) , but do not contribute directly to catalytic activity of electron transport and atp synthesis ( de grassi et al .
thus , maintenance of primary amino acid sequences of the noncore subunits appears to be less important to oxphos function than that of the core subunits and the intensity of selection on these proteins is therefore reduced .
it is possible that different domains within a protein may be under different selective regimes and such domains may have different evolutionary rates .
such differences could be obscured by combining domains and estimating dn for each gene as a whole .
for example , a gene with a small but highly functionally constrained domain combined with large relatively unconserved domains would exhibit high dn despite having an important conserved function . had we obtained high relative dn for core oxphos genes , such a result would have been consistent with a compensatory evolution hypothesis but could have been due to such overestimation of dn . however , a salient result here is the relatively low estimates of dn for several mt and core nu oxphos genes .
it seems unlikely that differential selection on domains could lead to underestimation of dn if the compensatory evolution hypothesis ( which predicts elevated dn values ) was a major factor influence in oxphos evolution .
we can not exclude the possibility of a few compensatory changes in nu oxphos genes as they  keep pace  with changes in mt oxphos genes because assessment of evolutionary rates does not allow examination of individual sites .
neither can we reject compensatory changes in mt genes in response to changes in the nu genome ; indeed , the rapid rate of mt genome evolution could allow these genes to more readily respond to changes in nu components . nonetheless , compensatory changes in either direction do not appear to occur with a frequency sufficient to contribute to differences in evolutionary rates among genes or groups of genes . therefore , despite the necessary interaction of the products of two different genomes , oxphos evolution appears to be driven largely by conventional natural selection for functional efficiency acting on individual subunits , regardless of their genome of origin .
supplementary tables s1 and s2 are available at genome biology and evolution online ( http://www.gbe.oxfordjournals.org/ ) .
 oxidative phosphorylation ( oxphos ) , the major energy - producing pathway in aerobic organisms , includes protein subunits encoded by both mitochondrial ( mt ) and nuclear ( nu ) genomes . 
 how these independent genomes have coevolved is a long - standing question in evolutionary biology . 
 although mt genes evolve faster than most nu genes , maintenance of oxphos structural stability and functional efficiency may involve correlated evolution of mt and nu oxphos genes . 
 the nu oxphos genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt oxphos subunits with which they interact . 
 evolutionary rates of nu oxphos genes should , therefore , be higher than that of nu genes that are not involved in oxphos ( nu non - oxphos ) . 
 we tested the compensatory evolution hypothesis by comparing the evolutionary rates ( synonymous substitution rate ds and nonsynonymous substitution rate dn ) among 13 mt oxphos genes , 60 nu oxphos genes , and 77 nu non - oxphos genes in vertebrates ( 7 fish and 40 mammal species ) . 
 the results from a combined analysis of all oxphos subunits fit the predictions of the hypothesis . 
 however , results from two oxphos complexes did not fit this pattern when analyzed separately . 
 we found that the dn of nu oxphos genes for  core  subunits ( those involved in the major catalytic activity ) was lower than that of  noncore  subunits , whereas there was no significant difference in dn between genes for nu non - oxphos and core subunits . 
 this latter finding suggests that compensatory changes play a minor role in the evolution of oxphos genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits . 
Introduction
Materials and Methods
Results
Discussion
Supplementary Material