- Research article
- Open Access
Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance
BMC Evolutionary Biology volume 15, Article number: 173 (2015)
Oxidative phosphorylation (OXPHOS) is the primary source of ATP in eukaryotes and serves as a mechanistic link between variation in genotypes and energetic phenotypes. While several physiological and anatomical factors may lead to increased aerobic capacity, variation in OXPHOS proteins may influence OXPHOS efficiency and facilitate adaptation in organisms with varied energy demands. Although there is evidence that natural selection acts on OXPHOS genes, the focus has been on detection of directional (positive) selection on specific phylogenetic branches where traits that increase energetic demands appear to have evolved. We examined patterns of selection in a broader evolutionary context, i.e., on multiple lineages of fishes with extreme high and low aerobic performance.
We found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. Positive selection is not consistently associated with high performance taxa and appears to be strongest on lineages containing low performance taxa. In contrast, within high performance lineages, purifying (negative) selection appears to predominate.
We provide evidence that selection on OXPHOS varies in both form and intensity within and among lineages through evolutionary time. These results provide evidence for fluctuating selection on OXPHOS associated with divergence in aerobic performance. However, in contrast to previous studies, positive selection was strongest on low performance taxa suggesting that adaptation of OXPHOS involves many factors beyond enhancing ATP production in high performance taxa. The broader pattern indicates a complex interplay between organismal adaptations, ATP demand, and OXPHOS function.
Physiological processes may serve as mechanistic links between genotypes and organismal phenotypes. Accordingly, adaptations in genes of energy metabolism pathways may facilitate the evolution of organismal structures and life habits with diverse energy requirements. Proteins encoded by the mitochondrial genome serve as core subunits of the oxidative phosphorylation (OXPHOS) system, the primary source of ATP in eukaryotic cells. Consequently, organismal traits with differing ATP demands may be influenced by adaptation in mitochondrial genomes [1–5]. Our understanding of patterns and rates of adaptive change in mitochondrial genomes remains limited despite extensive use of mitochondrial genes as molecular markers in evolutionary studies  and the important role of mitochondria in many human pathologies . Although there is ample evidence of natural selection acting on mitochondrial genes [8–15], the functional significance of adaptive mitochondrial change is rarely known.
Evidence of positive selection on OXPHOS genes has been associated with evolution of a variety of energetically demanding characteristics (reviewed in ), including origin of large brains in anthropoid primates , powered flight in bats , and adaptation to cold environment in polar bears . These suggest a significant role for OXPHOS in organismal adaptation, and because divergence among lineages often involves traits with different energy usage, OXPHOS evolution may be an important factor in the diversification of life.
OXPHOS functional efficiency may be particularly important to energy intensive processes such as locomotion. Variation in locomotive performance of fishes is among the most extreme among vertebrates, ranging from largely sedentary filter feeders and sit-and-wait predators to highly migratory species and active pelagic foragers. For example, seahorses and flounders spend much of their time nearly motionless, while tunas and marlins are “high-performance” swimmers that exhibit high aerobic metabolism and prolonged fast swimming [20, 21]. Highly-active fish taxa exhibit many morphological and physiological adaptations that enhance swimming performance (reviewed in ). Such adaptations include modifications of body shape and hydrodynamics [23–25], swimming form and mechanical kinematics [26–28], muscle composition [29–31], metabolic rates , heart volume and aerobic capacity [32–34], and mitochondrial structure and concentration . However, much less is known about the molecular adaptations that influence organismal performance.
At the molecular level, expression levels of proteins directly involved in energy metabolism may be increased in highly mobile fish species. Tunas and marlins, which have higher cruising speeds than other active fishes , have been shown to have elevated myoglobin  and higher concentrations of metabolic enzymes in heart and skeletal muscle (reviewed in ). Elevated activities of citrate synthase (which catalyzes the first reaction of the Krebs cycle), carnitine-palmitoyl transferase, and 3-hydroxy-o-acyl-CoA dehydrogenase (rate-limiting enzymes in fatty acid oxidation), reflect the increased aerobic metabolic potential of scombrid fishes . In addition, OXPHOS genes are differentially expressed between morphs of lake whitefish that differ in activity levels in foraging behavior .
Alternatively, or in addition, to variation in gene expression and post-transcriptional modification, divergent energy demands may lead to adaptive evolution in the structure of specific OXPHOS proteins. Variation in OXPHOS proteins could influence the efficiency of ATP production by affecting how tightly electron transport and proton pumping are coupled in the electron transport chain. Modifications in the structure of OXPHOS complexes I–IV caused by amino acid substitutions in constituent proteins could affect “slip reactions”, resulting in more or less protons pumped by the electron-transport-chain for each electron pair transferred (the H+/2e ratio) (reviewed in ). Alternatively, substitutions in proteins of ATP synthase could modify the amount of ATP made for each proton driven through it (the H+/ATP ratio) . These phenomena are consistent with prior reports of positive selection in the cytochrome c oxidase genes of the high performance billfishes [39, 40], the MT-ND2 and MT-ND5 genes of some migratory Pacific salmon , and MT-ND2, MT-ND4 and MT-ND5 genes of pelagic Atlantic herring .
Here, we investigate patterns of adaptation on OXPHOS genes in diverse group of fish taxa with different swimming performance. We hypothesized that positive natural selection would affect OXPHOS efficiency among divergent lineages with long-term differences in ATP demand. Thus molecular adaptation was predicted to be associated with locomotor intensity (speed, duration and frequency). We examined evidence for positive selection on all mitochondrial OXPHOS genes from six fish groups that can be classified into three different swimming performance categories based on general locomotion patterns . Tunas and billfishes represent high performance swimmers, mackerels and jacks represent moderate (or moderate-high) performance swimmers, and flatfishes and seahorses + pipefishes represent low performance swimmers. Pelagic fishes exhibit highly aerobic locomotion with greater endurance than sedentary fishes, and among the pelagic fishes, tunas and billfishes may maintain the highest speeds for the longest duration . Seahorses and pipefishes exhibit much lower frequency, duration, and speed of locomotion, while flatfishes exhibit lower frequency, and sustained speed. These factors as well as their unusual swimming forms and kinematics suggest that much less of their total energy budget is devoted to locomotion, so we categorized them as low-performance swimmers. Recent phylogenetic analyses [43–45] indicate that these taxa are arranged into two monophyletic groups each containing representatives of all three performance classes (Fig. 1). Two previous studies (39, 40) examined positive selection in some high-performance fishes that are included in this study. We extended these studies by including all 13 mitochondrial OXPHOS genes and examined the evidence of positive selection on multiple phylogenetic branches, including six independent fish lineages as well as some ancestral branches. In contrast to previous studies, this approach allowed us to make mitochondrial genome-wide assessments of selection across taxa exhibiting a wide range of locomotor performance. We also inferred functional significance from the position of positively selected amino acid sites in the 3-dimensional structure of specific enzyme complexes. Our results indicate that selection on OXPHOS genes is indeed associated with divergent swimming habits among these fishes. However, selection is heterogeneous over evolutionary time and positive selection is not strictly associated with high performance taxa. Thus we provide new insights on the evolution of swimming diversity in fishes and the adaptation of OXPHOS genes relative to organismal energetic performance on a broad phylogenetic scale.
We acquired sequences of 13 mitochondrial protein genes (all of them encode subunits in OXPHOS) from the representative species of each fish group from National Center for Biotechnology Information (NCBI) GenBank (Additional file 1: Table S1). To complement the mitochondrial genes, we added 7 nuclear genes which are commonly used in fish phylogenies (obtained from NCBI GenBank) or were developed as part of the Fish Tree of Life project . These genes are rag1 (recombination activating protein 1), rag2 (recombination activating protein 2), rhodopsin, tmo4c4 (anonymous, see ), zic1 (zinc finger protein 1), myh6 (myosin, heavy chain six), and btbd7 (BTB domain containing seven). For species without available sequences, we designed primers (Additional file 1: Table S2) and amplified and sequenced specific genes via standard polymerase chain reaction (PCR). PCR product was sequenced in both forward and reverse directions using ABI BigDye terminator chemistry and an ABI Prism 3130 XL Genetic Analyzer. Sequences have been deposited in GenBank (Additional file 1: Table S3). Particular genes for a few species could not be amplified. In such cases, sequences were obtained from congeners, yielding 19 “chimaeric” individuals. The 20 genes (13 mitochondrial and 7 nuclear) were concatenated for phylogenetic analyses .
Maximum likelihood trees (with 1000 bootstrap replicates) were estimated in the program RAxML v.7.0.4 . Bayesian phylogenetic analysis was performed with MrBayes v.3.1 [50, 51]. Swimming performance states were reconstructed from extant taxa with parsimony and likelihood using outgroups from  in Mesquite .
Analysis of patterns of natural selection
Sequences of 13 mitochondrial OXPHOS genes were used to examine patterns of natural selection. We first ran a series of random sites models (M0, M1a, M2a, M3, M7, M8a and M8) implemented in Codeml of PAML v.4.7 . Likelihood ratio tests (LRTs) were conducted on the likelihood values produced by specific pairs of models: M3 vs. M0, M2a vs. M1a, M8 vs. M7, and M8 vs. M8a. Next we ran branch-site models on designated branches with three different starting ω values (0.1, 1.0, and 4.3). We had no prior knowledge of which branch(es) could have experienced positive natural selection (except for branches leading to tunas and billfishes); therefore, we designated one branch (from b1-25 in Fig. 2) as the foreground branch in each test. LRTs were performed to determine if the more complex model A is significantly better than the null model. In the branches that showed evidence of positive selection, we used Bayes empirical Bayes (BEB) to calculate the probability of amino acid sites are under positive selection.
The 13 mitochondrial gene sequences were also analyzed using TreeSAAP , which measures the selective influences on 31 structural and biochemical amino acid properties, and performs goodness-of-fit and categorical statistical tests. We used a sliding window size of 20 amino acids (Additional file 1: Figure S2). Because of the large number of sites under positive selection, we only show sites with significance when its p <0.001 for clarity. Then we used a sliding window size of 1 and one amino acid site is defined as being positively selected as long as one of the 31 properties showed significance.
Mapping positively selected amino acid sites onto 3-dimensional (3D) crystal protein structure
The MT-ND genes that encode subunits in OXPHOS complex I are considered highly conserved between prokaryotes and vertebrates . We inferred the potential function of positively selected amino acid sites belonging to MT-ND subunits by mapping them onto the 3D crystal structure of Thermus thermophiles (PDB ID: 4HEA), a Gram negative eubacterium . We mapped the positively selected sites of MT-CYB onto chicken bc1 complex C chain (PDB ID: 1BCC) crystal structure because a study suggests the structure of the catfish MT-CYB protein resembles that of chickens . Similarly, we mapped the positively selected sites onto bovine MT-CO1-3 subunits 3D structure (PDB ID: 1OCC) because of the highly conserved structure from bacteria to bovine [58, 59]. All the mapping was conducted using Geneious pro (v.7.0).
Phylogeny and ancestral state reconstruction
Both maximum likelihood and Bayesian analyses generated the same topology (Fig. 2) regardless of partition scheme with relatively high bootstrap support and posterior probabilities (see Additional file 1: Figure S1). The recent comprehensive phylogenetic analysis of  resolved the phylogeny of bony fishes to the level of taxonomic family, while a study by  focused on relationships of several pelagic fish groups. Our results are consistent with both of these studies for the taxa in common to each. Specifically, we recovered two monophyletic groups, each containing taxa with low, moderate, and high swimming performance. In one group, the tunas and mackerels (Scombriformes) are sister to the seahorses and pipefishes (Syngnathiformes), while the series Carangaria contains the billfishes (Istiophoriformes), jacks (Carangiformes), and flatfishes (Pleuronectiformes)  (see also Phylogenetic Classification of Bony Fishes--Version 3: www.deepfin.org).
Swimming performance states were reconstructed from extant taxa with parsimony and likelihood using outgroups from . These results of reconstructed swimming performance suggest that the ancestors of each major group were low performance swimmers (Fig. 1), although incomplete taxon sampling of related groups could mislead the reconstructions. Even if the ancestors were moderate performance swimmers, it is clear that high performance swimming evolved independently in the tunas and billfishes.
Analysis of positive selection
Positive selection under site models cannot identify particular branches where positive selection has occurred, but they can detect positive selection among sites where it occurs as long-term trends or among multiple branches separated on the tree. A significant difference was detected between model M0 and M3, which suggests ω is variable among sites (Table 1). We used three model pairs to test for positive selection: M1a vs. M2a, M7 vs. M8, and M8a vs. M8. Neither M2a-M1a nor M8a-M8 showed significant difference. However, M8-M7 showed a significant difference, suggesting some positive selection signal on certain sites somewhere on the tree.
Table 2 lists the resulting likelihood values, likelihood ratio tests (LRTs), and estimated model parameters for each branch examined under branch-site models. We performed a false discovery rate analysis  on these results where we performed 20 tests and recovered 12 positive results. All 12 positive results had a q value of 0.03 or less, and the probability of a false positive among them is 0.36. Because this is less than 1.0, no false positives are expected. The number of amino acid sites on each branch inferred to be under positive selection via Bayes Empirical Bayes is provided Additional file 1: Table S4 (the identity of these sites and more details are provided in Additional file 1: Table S5). These results showed substantial differences in natural selection on lineages leading to taxa with different swimming performance. We found significant evidence of positive selection on the ancestral lineages of all three performance categories (branches b2, b17, b3, and b4 in Fig. 2), lineages leading to moderate-performance swimmers (b6 and b13), and low-performance swimmers (b16, b23, b24, b25, and b11). Conversely, strong purifying (negative) selection was identified on lineages of high-performance swimmers (b7, b8, b9, b10, b14, and b18). Fourteen sites that appeared to be positively selected occur on multiple branches; however, they are not associated with branches leading to a particular performance group (see details in Additional file 1: Table S5).
The above tests assess the strength of the evidence for positive selection rather than the strength of selection itself. An indicator of the strength of positive selection on a branch can be obtained from the product of the proportion of sites having ω > 1 and the estimated ω value for those sites. In Fig. 2, the branch width is shown proportional to the strength of positive selection based on this measure (except for b12, b23 and b24 which had extremely high ω and a cap of ω = 50 is used). Among those branches harboring sites under positive selection, branches ancestral to more than one swimming category (b2, b3, b4, b12, b17) and branches leading to low-performance fishes (b11 and b23) exhibited greater selection than branches leading to moderate-performance fishes (b6 and b13), and substantially more than any branches associated with high-performance fishes. This is counter to the notion that high performance swimmers will have been most strongly influenced by positive selection for enhanced OXPHOS performance assuming enhanced OXPHOS efficiency is the major contributor to high performance. However, there are alternative contributors to high performance such as increased mitochondrial density per tissue mass, more closely packed inner mitochondrial membrane cristae, increased metabolic enzyme activity, increased expression of genes involved in a number of biological pathways such as glycolysis, protein biosynthesis, and cytoskeletal structure. Under any of the later scenarios, an increase in OXPHOS efficiency might not be necessary for high performance.
The physico-chemical properties analysis as implemented by TreeSAAP  does not provide inferences about particular branches on the phylogeny, but does provide information about changes of particular amino acid sites across the whole tree. A large number of amino acid sites were identified as having substitutions with significant physico-chemical differences; presented in Additional file 1: Figure S2 and Table S4. A mitochondrial genome-wide (only the protein coding genes) sliding window analysis (window size = 20 amino acids) revealed that the vast majority of windows contained substitutions that exhibited between 0 and 15 properties with significant differences (Additional file 1: Figure S2). Windows with substitutions exhibiting between 15 and 30 significant properties occurred at much lower frequency. Although windows with the greatest number of significant property differences were observed in the genes for MT-ND1, MT-CO1, MT-ND4 and MT-CYB, there were no other apparent patterns of variation in the distribution of such sites across the genome (Additional file 1: Figure S2). Positively-selected sites identified by both methods are summarized in Additional file 1: Table S4, S5.
The two lineages, seahorses and flatfishes, on which positive selection was most pronounced, have also experienced extraordinary morphological evolution [61, 62]. Thus, it is possible that the positive selection signal detected on mitochondrial OXPHOS genes in these lineages is unrelated to energy demands, but is simply a consequence of rapid genomic evolution in these groups. The probability of nucleotide change is indeed higher in these two lineages than for the other major groups as indicated by five nuclear genes that are not directly involved in OXPHOS system (rag1, rhodopsin, tmo4c4, mhy6, and zic1). However, a positive selection signal was found for only two genes, rag1 and rhodopsin, on one branch leading to flatfishes (LRTs = 72.53, ω = 50.88, three sites 199, 425, and 605 under selection with bayes empirical bayes (BEB) probability of 0.684, 0.501, and 0.958, respectively). Because there was no consistent pattern of selection on these non-OXPHOS genes, it is suggested that the positive selection detected on mitochondrial genes does not reflect a genome-wide pattern of divergence but may be related to adaptation of OXPHOS efficiency.
Structural position of positively selected sites
The position of particular amino acids in the tertiary and quaternary structure of a protein may allow inferences about the function of individual residues. In particular, those sites near the catalytic core or other functionally important regions, or those in physical proximity (likely to interact with) other amino acids, would seem most likely to influence protein function.
Complex I performs the first step and is the largest and most complicated enzyme complex in the OXPHOS pathway. It catalyzes the transfer of two electrons from NADH to ubiquinone (Q), coupled to the translocation of four protons across the inner mitochondrial membrane. It is also a major source of reactive oxygen species in mitochondria. Complex I exhibits an L-shaped architecture with a membrane arm and a hydrophilic peripheral arm that protrudes into the mitochondrial matrix (Fig. 3a). The membrane arm consists of 7 core mitochondrial NADH dehydrogenase (MT-ND) gene encoded subunits.
Site 12, one of the sites identified in subunit ND1 as being positively selected (Fig. 3b), is included in the region associated with MELAS (mitochondrial encephalomyopathy, lactic acidosis, and stroke-like episodes)/DEAF enhancer/hypertension  and sudden infant death . The other positively selected sites, 172 and 192, are in the alpha helix close to a critical site, 184, that is associated with adult onset dystonia  (in yellow in Fig. 3b).
Subunit ND5 contains an unusual structural element, the helix HL (see enlarged Fig. 3c, indicated in yellow), that extends nearly the entire length of the membrane domain and coordinates conformational changes. On the opposite side of the membrane domain, a series of β-hairpins (βH element) (indicated in white in Fig. 3c, d) from neighboring subunits contribute to conformational changes and stability of the complex. Positively-selected sites include 515 [ND5], which is adjacent to helix HL; 53 [ND5] and 54 [ND4] are in βH elements; 63 [ND5] is in the seven-residue loop connecting two βH elements; 86 [ND2] (not shown here) and 71 [ND5] are in regions (including residues 83 [ND2] and 88 [ND5]) demonstrated to have significant negative effects on function .
In addition to the sites listed above, there are several sites that appear to be under positive selection, yet their location in the molecular structure does not provide any clear suggestion of functional significance. These residues could be involved in stabilizing the tertiary or quaternary structure of the various multi-subunit complexes or facilitate their assembly. These include residues 199 [ND5], 269 [ND5], 425 [ND5], 515 [ND5], 600 [ND5], 20 [ND4], 86 [ND2], 168 [ND2], 359 [ND2], 360 [ND2], 408 [ND2], 23 [ND4L], 4 [ND6], 8 [ND6], 11 [ND6], 90 [ND6], 138 [ND6], 94 [ND3] and 98 [ND3] (not shown in figures for clarity).
Complex III is an intermediate component of the respiratory chain, which transfers electrons from reduced ubiquinone to cytochrome c, coupled to proton translocation across the mitochondrial membrane . The MT-CYB gene encoded cytochrome b forms the active redox center: a cavity surrounded by the transmembrane helices A (residues 33–54), D (172–204), and E (221–245), and the amphipathic surface helix a (65–72) (indicated in white in Fig. 4). Residues 221 and 194 are also close enough to contact the active site inhibitor . The positively selected sites 194 and 235 (indicated in cyan in Fig. 4) are within these regions, suggesting important functional effects.
The pocket bound by stigmatellin, one OXPHOS inhibitor, is formed by the end of helix C, the helix cd1, the helix ef linker, and the end of helix F. Specific residues of known importance include 271, 275, 125–129, 138–153 (indicated in green in Fig. 4). The positively selected sites 126 and 297 are close to this area.
Complex IV is the terminal component of the respiratory chain, in which electrons received from cytochrome c reduce molecular oxygen to water and protons are pumped into the intermembrane space. Six positively selected sites were identified in complex IV: 178 in subunit COX1, 54 and 187 in subunit COX2, and 55, 155, and 171 in subunit COX3 (Fig. 5a).
Subunit I (encoded by MT-CO1 gene) is largely embedded in the membrane with three redox centers: heme a, heme a3, and CuB. One of the CuB liganded residues is 240. It has been suggested that residue 244 and 240 are close enough for binding interaction, thereby forming an unusual cross-link structure . One positively selected site 178 is very close to 244 and 240 (Fig. 5b). Two proton translocation pathways, D-pathway and K-pathway, are identified in subunit I. Three key residues of D-pathway 91, 98, and 242, and one key residue of K-pathway 319 are well defined from mutagenesis studies  (key residues are indicated in white in Fig. 5b). The positively selected site 178 (indicated in cyan in Fig. 5b) is close to (within 30 Å) the four key residues and could potentially affect proton translocation.
Subunit II (encoded by MT-CO2 gene) has two trans-membrane helices and a hydrophilic beta strand extending to the extra-membrane domain, housing the CuA center. The first residue crucial for electron to enter the oxidase complex cytochrome c is 106 (indicated in magenta in Fig. 5c). The other CuA liganded residues identified through mutagenesis include 196, 200, 198, 161, 204, and 207 (indicated in magenta in Fig. 5c). One positively selected site 187 is very close to the CuA center and the other positively selected site 54 is in a linking strand between two helices of the entry site which includes critical residue 62 (indicated in cyan in Fig. 5c).
Subunit III (encoded by MT-CO3 gene) is fully embedded in the membrane domain. No key residue of this subunit has been identified through mutagenesis. However, it has been shown that this subunit stabilizes the integrity of the binuclear center in subunit I; for example, when this gene is deleted, only a partially assembled complex results .
The sites described above as having potential functional significance were more or less uniformly distributed among branches of the tree and showed no association with particular swimming performance groups. Many sites exhibited evidence of being under positive selection yet did not appear in structural locations that would suggest special functional importance. However, the interaction between amino acids among different proteins or among sites in the same protein is known to have significant functional effects (e.g., [70, 71]). Such interactions could affect protein assembly, stability or specific function. We used 4 Å as the nominal upper limit for weak interactions between amino acid sites as described in . None of the sites identified here as under positive selection were found to be within 4 Å of any other amino acid sites.
We examined evidence for adaptation of mitochondrial OXPHOS genes in fishes with different swimming performance. Selection was investigated for specific amino acids sites across the whole phylogenetic tree for these species, as well as for amino acid sites on individual branches of the tree. The results show a strong signal of positive selection on branches in the ancestral parts of the tree, and branches leading to low- and moderate-performance swimmers. However, no evidence of positive selection was observed within clades of the high-performance tunas or billfishes. We did not observe a disproportionate effect of selection on any particular gene as all genes exhibited some positively selected sites but this varied across branches of the tree. We also show that many of the sites identified as being under positive selection occur in structural regions where they may have effects on OXPHOS function. Positively selected sites in other regions may also have functional significance, but their potential effects are less clear.
Our results show the strongest signal of positive selection on branches leading to the lowest performance fishes, while purifying selection was identified on the branches of high performance fishes. This result is in strong contrast to previous studies that have exclusively examined the hypothesized association of positive selection on OXPHOS with increased energetic demands required by higher organismal performance (reviewed in ). Our results suggest that a fairly efficient OXPHOS system had evolved under positive selection in the ancestors of the two major groups. Selective modification of OXPHOS on branches leading directly to tunas and billfishes may have further facilitated the evolution of high performance swimming. But once very high performance swimming had evolved in these lineages, purifying selection appears to have predominated on the OXPHOS system as it existed at that time. Conversely, substantial modification due to positive selection occurred in the lower performance lineages of flatfishes and the seahorse + pipefish group. The moderate-high performance jacks and mackerels exhibited moderate-high conservation, with limited positive selection. Thus, the strength of the positive selection signal is inversely proportional to swimming performance within both of these taxonomic groups.
A common expectation is that positive selection should lead to enhanced organismal performance. Moreover, it might be expected that positive selection will lead to increased functional efficiency of OXPHOS in response to the increased ATP demands associated with enhanced performance. In the absence of information on the exact functional significance of individual substitutions, the effects of positive selection to increase or decrease OXPHOS efficiency in these taxa remains unknown. However, our results do not match simple expectations. Positive selection appears to be heterogeneous, fluctuating over phylogenetic time scales, with no simple relationship between the strength of selection and organismal performance. However, one clear pattern is that there is strong functional constraint (negative selection) in these high performance systems and directional selection on these lower performance groups.
Our results appear to contrast with previous studies where positive selection or evolutionary rate variation was found on lineages leading to organisms with high ATP demands. For example, MT-CO1  and MT-CO2  exhibited accelerated d N in the lineage leading to hominids that appears to be associated with increased brain size. Grossman et al.  found accelerated rate variation for MT-CO4, a nuclear gene in OXPHOS complex IV in catarrhine ancestors of hominids in the period between 18 and 40 Mya and then decelerated in the descendant hominid lineages. On the lineage leading to bats, the only mammals capable of powered flight, eight OXPHOS genes were found to have undergone positive selection . Foote et al.  found two positively selected amino acid sites, which could influence overall metabolic performance, in the mitochondrial genes of killer whales (Orcinus orca). In addition, the MT-ND2 and MT-ND5 genes of highly migratory Pacific salmon exhibit evidence of positive selection .
Our results are consistent with previous studies of the high-performance taxa examined here. Dalziel et al.  examined one mitochondrial OXPHOS gene (MT-CO2) on several branches among high-performance fishes, including billfishes and tunas. They found ω was not increased in lineages leading to the tunas but was significantly increased in the lineage preceding the billfish (including several amino acid sites). However, the phylogeny used in  does not include the flatfish or seahorse clades, which are now recognized as close relatives of billfishes and tunas, respectively. Little et al.  examined three mitochondrial OXPHOS genes (MT-CO1, 2, and 3) on the single lineage leading to billfishes and found positive selection along that lineage when flatfishes were excluded but no selection was detected when flatfishes were included. The latter result is consistent with our findings and as suggested by Little et al.  emphasizes the importance of dense phylogenetic sampling for the analysis of positive selection. We extended these two studies by examining 10 additional mitochondrial OXPHOS genes in the high-performance taxa, and also by explicitly investigating selection on related low performance taxa.
With respect to negative selection on high performance fishes, we speculate that once a reasonably efficient OXPHOS system evolved it may have become difficult to change in the high performance groups. The high performance system might be expected to have a much lower tolerance for non-synonymous substitutions as most (even slight) changes would be likely to have negative functional effects. Conversely, lower performance swimmers may have a much broader tolerance for non-synonymous substitutions because the OXPHOS system is under much lower performance demands. We envision the high performance swimmers as occupying a local optimum on a fitness landscape , but their OXPHOS system is so fine-tuned that substitutions that would allow them to cross fitness valleys and reach higher peaks could be strongly deleterious in the short term. On the other hand, lower performance fishes might readily cross such valleys without significant fitness costs because OXPHOS efficiency will be less critical and they may then climb other fitness peaks due to positive selection.
It is also possible that organismal fitness in taxa with low energy demands may be increased by a modified OXPHOS regulatory system or even a reduction in OXPHOS efficiency. In such cases, we would expect to see evidence of positive selection on OXPHOS genes in these taxa. Because there may be trade-offs between OXPHOS rate or efficiency and deleterious effects, reducing OXPHOS efficiency may be adaptive in systems where ATP demand is chronically low. For example, maintenance of a strong chemiosmotic gradient in organisms with low ATP demand may cause increased production of reactive oxygen species (ROS)  leading to oxidative stress . Therefore reduction of the H+/2e ratio (increased slippage) due to altered protein structures could be adaptive in low performance species. OXPHOS regulation is highly complex and involves mechanisms independent of structural OXPHOS proteins. Elevated expression of OXPHOS genes and others involved in aerobic respiration could clearly increase aerobic capacity in the absence of selection on specific OXPHOS variants. However, it is clear that positive selection acts on OXPHOS proteins and the effects of positive selection are frequently associated with the evolution of differences in ATP demand.
There are other issues that could affect the efficiency of ATP generation of OXPHOS in high-performance tunas and billfishes. Tunas and billfishes exhibit high metabolic rates to generate ATP to support their high swimming performance (thus high OXPHOS efficiency). However, they are among the few endothermic species of fishes [43, 79, 80] that generate great amount of heat with high metabolic rates (thus reducing ATP generation). Thus, depending on the extent to which body heat is derived from OXPHOS uncoupling, there might be trade-offs between ATP generation and heat generation in these endothermic high-performance fishes. In addition, it has been shown that the heat through mitochondrial proton leak differs among organisms and even differs in cells and tissues in the same organism . It would be helpful to have empirical data about the P/O ratio (how many ATP molecules made from ADP for each oxygen atom consumed), the H+/O ratio (how many protons do mitochondria pump from matrix into inter-membrane space for each oxygen atom consumed), and the H+/ATP ratio (how many protons flow back to the matrix for each ADP molecule phosphorylated to ATP) in these fish species.
High mitochondrial respiration capacity can be achieved through means other than modification of OXPHOS efficiency via genetic variation, such as increased mitochondrial density per tissue mass, more closely packed inner mitochondrial membrane cristae, increased metabolic enzyme activity, increased expression of genes involved in a number of different but related biological pathways such as glycolysis, the Krebs cycle, and fatty acid metabolism. In addition, swimming performance, as in our qualitatively defined groups, may have arisen due to many different morphological and physiological characteristics beyond OXPHOS function. We do not assert that all variation in swimming performance is due to OXPHOS variation, only that some unknown fraction of the performance variation is positively associated with natural selection. A number of recent studies have presented evidence for direct effects of OXPHOS variation on organismal fitness, many with an identified mechanistic basis [81–86]. While these studies employ tractable model organisms such as Drosophila, they provide strong evidence that OXPHOS variation affects organismal performance and fitness in the wild.
In conclusion, we found that patterns of natural selection on mitochondrial OXPHOS genes are complex among fishes with different swimming performance. The type and direction of selection are heterogeneous through evolutionary time and vary in ways that would not be readily predicted based solely on organismal performance. The most striking result was extensive positive selection on low-performance swimmers. Although examination of the most recent lineages indicates that positive selection is inversely proportional to organismal performance, the broader pattern indicates a complex interplay between organismal adaptations, ATP demand and OXPHOS function through evolutionary time.
Availability of supporting data
The data sets supporting the results of this article are included in its additional files.
Mitochondrially encoded cytochrome c oxidase I
Mitochondrially encoded cytochrome c oxidase II
Mitochondrially encoded cytochrome c oxidase III
Mitochondrially encoded NADH dehydrogenase 1
Mitochondrially encoded NADH dehydrogenase 2
Mitochondrially encoded NADH dehydrogenase 3
Mitochondrially encoded NADH dehydrogenase 4
Mitochondrially encoded NADH dehydrogenase 4 L
Mitochondrially encoded NADH dehydrogenase 5
Mitochondrially encoded NADH dehydrogenase 6
Mitochondrial I: NADH dehydrogenase (ubiquinone) subunits
Mitochondrially encoded cytochrome b
Mitochondrially encoded ATP synthase 6
Mitochondrially encoded ATP synthase 8
Rand DM. Thermal habit, metabolic rate and the evolution of mitochondrial DNA. Trends Ecol Evol. 1994;9(4):125–31.
Ballard JWO, Whitlock MC. The incomplete natural history of mitochondria. Mol Ecol. 2004;13(4):729–44.
da Fonseca RR, Johnson WE, O’Brien SJ, Ramos MJ, Antunes A. The adaptive evolution of the mammalian mitochondrial genome. BMC Genet. 2008;9(1):119.
Ballard J, Melvin R. Linking the mitochondrial genotype to the organismal phenotype. Mol Ecol. 2010;19(8):1523–39.
Hill GE. Cellular respiration: the nexus of stress, condition, and ornamentation. Integr Comp Biol. 2014;54(4):645–57.
Ballard JWO, Rand DM. The population biology of mitochondrial DNA and its phylogenetic implications. Annu Rev Ecol Evol Syst. 2005;621–642.
Wallace DC. Mitochondrial DNA mutations in disease and aging. Environ Mol Mutagen. 2010;51(5):440–50.
Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, et al. Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci. 2003;100(1):171–6.
Rand D, Haney R, Fry A. Cytonuclear coevolution: the genomics of cooperation. Trends Ecol Evol. 2004;19(12):645–53.
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–6.
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.
Bazin E, Glémin S, Galtier N. Population size does not influence mitochondrial genetic diversity in animals. Science. 2006;312(5773):570–2.
Ballard JWO, Melvin RG, Katewa SD, Maas K. Mitochondrial DNA variation is associated with measurable differences in life‐history traits and mitochondrial metabolism in Drosophila simulans. Evolution. 2007;61(7):1735–47.
Meiklejohn CD, Montooth KL, Rand DM. Positive and negative selection on the mitochondrial genome. Trends Genet. 2007;23(6):259–63.
Hassanin A, Ropiquet A, Couloux A, Cruaud C. Evolution of the mitochondrial genome in mammals living at high altitude: new insights from a study of the tribe Caprini (Bovidae, Antilopinae). J Mol Evol. 2009;68(4):293–310.
Garvin MR, Bielawski JP, Sazanov LA, Gharrett AJ. Review and meta‐analysis of natural selection in mitochondrial complex I in metazoans. J Zool Syst Evol Res. 2015;53(1):1–17.
Doan JW, Schmidt TR, Wildman DE, Uddin M, Goldberg A, Hüttemann M, et al. Coadaptive evolution in cytochrome c oxidase: 9 of 13 subunits show accelerated rates of nonsynonymous substitution in anthropoid primates. Mol Phylogenet Evol. 2004;33(3):944.
Shen Y-Y, Liang L, Zhu Z-H, Zhou W-P, Irwin DM, Zhang Y-P. Adaptive evolution of energy metabolism genes and the origin of flight in bats. Proc Natl Acad Sci. 2010;107(19):8666–71.
Welch AJ, Bedoya-Reina OC, Carretero-Paulet L, Miller W, Rode KD, Lindqvist C. Polar bears exhibit genome-wide signatures of bioenergetic adaptation to life in the Arctic environment. Genome Biol Evol. 2014;6(2):433–50.
Beamish F. Swimming capacity. New York: Academic; 1978.
Block BA, Booth D, Carey FG. Direct measurement of swimming speeds and depth of blue marlin. J Exp Biol. 1992;166(1):267–84.
Lauder GV. Fish locomotion: recent advances and new directions. Ann Rev Mar Sci. 2015;7:521–45.
Blake R. Energetics of leaping in dolphins and other aquatic animals. J Mar Biol Assoc U K. 1983;63(01):61–70.
Vogel S. Life in moving fluids: the physical biology of flow. Princeton: Princeton University Press; 1994.
Webb PW. The biology of fish swimming, mechanics and physiology of animal swimming. 1994. p. 4562.
Videler J, Nolet B. Costs of swimming measured at optimum speed: scale effects, differences between swimming styles, taxonomic groups and submerged and surface swimming. Comp Biochem Physiol A Physiol. 1990;97(2):91–9.
Wainwright PC, Bellwood DR, Westneat MW. Ecomorphology of locomotion in labrid fishes. Environ Biol Fish. 2002;65(1):47–62.
Walker JA, Westneat MW. Performance limits of labriform propulsion and correlates with fin shape and motion. J Exp Biol. 2002;205(2):177–87.
Graham J, Koehrn F, Dickson K. Distribution and relative proportions of red muscle in scombrid fishes: consequences of body size and relationships to locomotion and endothermy. Can J Zool. 1983;61(9):2087–96.
Moon T, Altringham J, Johnston I. Energetics and power output of isolated fish fast muscle fibres performing oscillatory work. J Exp Biol. 1991;158(1):261–73.
Tang J, Wardle C. Power output of two sizes of Atlantic salmon (Salmo salar) at their maximum sustained swimming speeds. J Exp Biol. 1992;166(1):33–46.
Shiels H, Galli G, Block B. Cardiac function in an endothermic fish: cellular mechanisms for overcoming acute thermal challenges during diving. Proc R Soc B Biol Sci. 2015;282(1800):20141989.
Brill RW, Bushnell PG. Metabolic and cardiac scope of high energy demand teleosts, the tunas. Can J Zool. 1991;69(7):2002–9.
Farrell AP. Invited perspectives in physiological zoology. From hagfish to tuna: a perspective on cardiac function in fish. Physiol Zool. 1991;64(5):1137–64.
Moyes CD, Mathieu-Costello OA, Brill RW, Hochachka PW. Mitochondrial metabolism of cardiac and skeletal muscles from a fast (Katsuwonus pelamis) and a slow (Cyprinus carpio) fish. Can J Zool. 1992;70(6):1246–53.
Dickson KA. Unique adaptations of the metabolic biochemistry of tunas and billfishes for life in the pelagic environment. Environ Biol Fish. 1995;42(1):65–97.
Evans M, Bernatchez L. Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and nonparallel physiological divergence. J Evol Biol. 2012;25(9):1823–34.
Brand M. The efficiency and plasticity of mitochondrial energy transduction. Biochem Soc Trans. 2005;33(5):897–904.
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–31.
Little A, Lougheed S, Moyes C. Evolution of mitochondrial-encoded cytochrome oxidase subunits in endothermic fish: The importance of taxon-sampling in codon-based models. Mol Phylogenet Evol. 2012;63(3):679–84.
Garvin MR, Bielawski JP, Gharrett AJ. Positive Darwinian selection in the piston that powers proton pumps in complex I of the mitochondria of Pacific salmon. PLoS One. 2011;6(9), e24127.
Teacher AG, Andre C, Merila J, Wheat CW. Whole mitochondrial genome scan for population structure and selection in the Atlantic herring. BMC Evol Biol. 2012;12(1):248.
Little A, Lougheed S, Moyes C. Evolutionary affinity of billfishes (Xiphiidae and Istiophoridae) and flatfishes (Plueronectiformes): independent and trans-subordinal origins of endothermy in teleost fishes. Mol Phylogenet Evol. 2010;56(3):897–904.
Betancur-R R, Broughton RE, Wiley EO, Carpenter K, López JA, Li C, et al. The tree of life and a new classification of bony fishes. PLoS Currents. 2013;5.
Miya M, Friedman M, Satoh TP, Takeshima H, Sado T, Iwasaki W, et al. Evolutionary origin of the scombridae (tunas and mackerels): members of a paleogene adaptive radiation with 14 other pelagic fish families. PLoS One. 2013;8(9), e73535.
Li C, Ortí G, Zhang G, Lu G. A practical approach to phylogenomics: the phylogeny of ray-finned fish (Actinopterygii) as a case study. BMC Evol Biol. 2007;7(1):44.
Streelman J, Karl SA. Reconstructing labroid evolution with single–copy nuclear DNA. Proc R Soc London, Ser B. 1997;264(1384):1011–20.
Kupczok A, Schmidt HA, von Haeseler A. Accuracy of phylogeny reconstruction methods combining overlapping gene data sets. Algorithms Mol Biol. 2010;5(1):1–17.
Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.
Huelsenbeck JP, Ronquist F, Nielsen R, Bollback JP. Bayesian inference of phylogeny and its impact on evolutionary biology. Science. 2001;294(5550):2310–4.
Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19(12):1572–4.
Maddison WP, Maddison D. Mesquite: a modular system for evolutionary analysis. 2001.
Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24(8):1586–91.
Woolley S, Johnson J, Smith MJ, Crandall KA, McClellan DA. TreeSAAP: selection on amino acid properties using phylogenetic trees. Bioinformatics. 2003;19(5):671–2.
Walker JE. The NADH: ubiquinone oxidoreductase (complex I) of respiratory chains. Q Rev Biophys. 1992;25(03):253–324.
Baradaran R, Berrisford JM, Minhas GS, Sazanov LA. Crystal structure of the entire respiratory complex I. Nature. 2013;494(7438):443–8.
Singh M, Gupta A, Lakra W. In silico 3-D structure prediction of cytochrome b protein of sisorid catfish Glyptothorax ngapang. Indian J Biotechnol. 2012;11(2):156–62.
Yoshikawa S, Shinzawa-Itoh K, Nakashima R, Yaono R, Yamashita E, Inoue N, et al. Redox-coupled crystal structural changes in bovine heart cytochrome c oxidase. Science. 1998;280(5370):1723–9.
Tsukihara T, Aoyama H, Yamashita E, Tomizaki T, Yamaguchi H, Shinzawa-Itoh K, et al. The whole structure of the 13-subunit oxidized cytochrome c oxidase at 2.8 Å. SCIENCE. 1996;272:1136–44.
Benjamini Y. Discovering the false discovery rate. J R Stat Soc Ser B (Stat Methodol). 2010;72(4):405–16.
Wilson A, Orr J. The evolutionary origins of Syngnathidae: pipefishes and seahorses. J Fish Biol. 2011;78(6):1603–23.
Yin M, Blaxter J. Morphological changes during growth and starvation of larval cod (Gadus morhua L.) and flounder (Platichthysflesus L.). J Exp Mar Biol Ecol. 1986;104(1):215–28.
Campos Y, Martin MA, Rubio JC, Olmo MC, Cabello A, Arenas J. Bilateral striatal necrosis and MELAS associated with a new T3308C mutation in the mitochondrial ND1 gene. Biochem Biophys Res Commun. 1997;238(2):323–5.
Opdal SH, Egeland T, Musse MA, Rognum TO. Possible role of mtDNA mutations in sudden infant death. Pediatr Neurol. 2002;27(1):23–9.
Simon DK, Friedman J, Breakefield XO, Jankovic J, Brin MF, Provias J, et al. A heteroplasmic mitochondrial complex I gene mutation in adult-onset dystonia. Neurogenetics. 2003;4(4):199–205.
Efremov RG, Sazanov LA. Structure of the membrane domain of respiratory complex I. Nature. 2011;476(7361):414–20.
Scheffler IE. Mitochondrial electron transport and oxidative phosphorylation. Mitochondria. 1999;141–245.
Zhang Z, Huang L, Shulmeister VM, Chi Y-I, Kim KK, Hung L-W, et al. Electron transfer by domain movement in cytochrome bc1. Nature. 1998;392(6677):677–84.
Richter O-M, Ludwig B. Cytochrome c oxidase—structure, function, and physiology of a redox-driven molecular machine, Reviews of physiology, biochemistry and pharmacology. Berlin: Springer; 2003. p. 47–74.
Gershoni M, Fuchs A, Shani N, Fridman Y, Corral-Debrinski M, Aharoni A, et al. Coevolution predicts direct interactions between mtDNA-encoded and nDNA-encoded subunits of oxidative phosphorylation complex I. J Mol Biol. 2010;404(1):158–71.
Willett CS, Burton RS. Environmental influences on epistatic interactions: viabilities of cytochrome c genotypes in interpopulation crosses. Evolution. 2003;57(10):2286–92.
Schmidt TR, Wu W, Goodman M, Grossman LI. Evolution of nuclear-and mitochondrial-encoded subunit interaction in cytochrome c oxidase. Mol Biol Evol. 2001;18(4):563–9.
Andrews TD, Easteal S. Evolutionary rate acceleration of cytochrome c oxidase subunit I in simian primates. J Mol Evol. 2000;50(6):562–8.
Adkins RM, Honeycutt RL. Evolution of the primate cytochrome c oxidase subunit II gene. J Mol Evol. 1994;38(3):215–31.
Grossman LI, Schmidt TR, Wildman DE, Goodman M. Molecular evolution of aerobic energy metabolism in primates. Mol Phylogenet Evol. 2001;18(1):26–36.
Foote AD, Morin PA, Durban JW, Pitman RL, Wade P, Willerslev E, et al. Positive selection on the killer whale mitogenome. Biol Lett. 2011;7(1):116–8.
Gavrilets S. Evolution and speciation on holey adaptive landscapes. Trends Ecol Evol. 1997;12(8):307–12.
Schieber M, Chandel NS. ROS function in redox signaling and oxidative stress. Curr Biol. 2014;24(10):R453–62.
Block BA, Finnerty JR, Stewart AF, Kidd J. Evolution of endothermy in fish: mapping physiological traits on a molecular phylogeny. Science. 1993;260:210.
Dickson KA, Graham JB. Evolution and consequences of endothermy in fishes. Physiol Biochem Zool. 2004;77(6):998–1018.
Montooth KL, Abt DN, Hofmann JW, Rand DM. Comparative genomics of Drosophila mtDNA: novel features of conservation and change across functional domains and lineages. J Mol Evol. 2009;69(1):94–114.
Pichaud N, Ballard JWO, Tanguay RM, Blier PU. Mitochondrial haplotype divergences affect specific temperature sensitivity of mitochondrial respiration. J Bioenerg Biomembr. 2013;45(1–2):25–35.
Meiklejohn CD, Holmbeck MA, Siddiq MA, Abt DN, Rand DM, Montooth KL. An incompatibility between a mitochondrial tRNA and its nuclear-encoded tRNA synthetase compromises development and fitness in Drosophila. PLoS Genet. 2013;9(1), e1003238.
Pichaud N, Ballard JWO, Tanguay RM, Blier PU. Naturally occurring mitochondrial DNA haplotypes exhibit metabolic differences: insight into functional properties of mitochondria. Evolution. 2012;66(10):3189–97.
Cooper BS, Burrus C, Ji C, Hahn MW, Montooth K. Similar efficacies of selection shape mitochondrial and nuclear genes in Drosophila melanogaster and Homo sapiens. BioRxiv. 2014. doi:10.1101/010355.
Greenlee KJ, Montooth KL, Helm BR. Predicting performance and plasticity in the development of respiratory structures and metabolic systems. Integr Comp Biol. 2014;54(2):307–22.
This work was supported by NSF award DEB-0732988 (to REB). We thank the Kansas University Biodiversity Institute & Natural History Museum Ichthyology Collections for some of the fish samples. We thank Z. Hao and T.E. Dowling for helpful discussions, and J.C. Cureton II, L.J. Weider, J.P. Masly, R. Knapp, and C.M. Lewis for comments on the manuscript.
Both authors declare that they have no competing interests.
FZ and RB designed the study, FZ collected and analyzed the data, FZ and RB wrote the manuscript. Both authors read and approved the final manuscript.
Additional file 1: Table S1.
Species list with mitochondrial genome accession number from NCBI. Table S2. PCR primers used for nuclear markers used in this study. Table S3. GenBank accession number of nuclear gene (rag1, rag2, tmo4c4, zic1, myh6, and btbd7) sequences of studied species. Figure S1. Phylogeny of six fish groups based on Bayesian and maximum likelihood methods. Figure S2. Summary of selection on all mitochondrial genes. Table S4. Number of amino acid sites showing positive selection signals in each mitochondrial gene on relevant branches. Table S5. Positively selected amino acid sites for each gene on various branches. (DOCX 351 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhang, F., Broughton, R.E. Heterogeneous natural selection on oxidative phosphorylation genes among fishes with extreme high and low aerobic performance. BMC Evol Biol 15, 173 (2015). https://doi.org/10.1186/s12862-015-0453-7
- Swimming Performance
- Amino Acid Site
- OXPHOS Complex
- Slide Window Size
- OXPHOS Gene