Research article | Open | Published:
Phylogenomics of the oxidative phosphorylation in fungi reveals extensive gene duplication followed by functional divergence
BMC Evolutionary Biologyvolume 9, Article number: 295 (2009)
Oxidative phosphorylation is central to the energy metabolism of the cell. Due to adaptation to different life-styles and environments, fungal species have shaped their respiratory pathways in the course of evolution. To identify the main mechanisms behind the evolution of respiratory pathways, we conducted a phylogenomics survey of oxidative phosphorylation components in the genomes of sixty fungal species.
Besides clarifying orthology and paralogy relationships among respiratory proteins, our results reveal three parallel losses of the entire complex I, two of which are coupled to duplications in alternative dehydrogenases. Duplications in respiratory proteins have been common, affecting 76% of the protein families surveyed. We detect several instances of paralogs of genes coding for subunits of respiratory complexes that have been recruited to other multi-protein complexes inside and outside the mitochondrion, emphasizing the role of evolutionary tinkering.
Processes of gene loss and gene duplication followed by functional divergence have been rampant in the evolution of fungal respiration. Overall, the core proteins of the respiratory pathways are conserved in most lineages, with major changes affecting the lineages of microsporidia, Schizosaccaromyces and Saccharomyces/Kluyveromyces due to adaptation to anaerobic life-styles. We did not observe specific adaptations of the respiratory metabolism common to all pathogenic species.
Oxidative phosphorylation (OXPHOS) is the primary energy-producing pathway in aerobic organisms . It functions by coupling the energy obtained from the oxidation of certain metabolic substrates to the phosphorylation of adenosine biphosphate (ADP) to produce ATP. This is achieved by a process of electronic transference through an intricate assembly of more than 20 discrete carriers. These carriers are mainly grouped into four membrane-embedded protein complexes, named Complex I through Complex IV, which form the electron transport chain (ETC). Some of the complexes in this chain are able to use the energy liberated by the electron transfer to the pumping of protons across the membrane, thereby generating a proton gradient. Finally, the energy obtained from the dissipation of this gradient is used by a fifth protein complex, ATP-synthase or Complex V, to synthesize ATP.
In eukaryotes, the oxidative phosphorylation machinery resides in the inner membrane of the mitochondrion. Molecular phylogenies of eukaryotic OXPHOS components indicate that the core subunits of the complexes were inherited from the alpha-proteobacterial ancestor of mitochondria [2, 3]. In contrast, other subunits might have different origins and show complex phylogenetic distributions . Besides providing important information on how complex systems evolve, knowledge about lineage-specific variations may serve to identify novel components or interactions. For instance, the evolutionary analysis of Complex I across a set of eighteen eukaryotes, lead to the prediction that the so-far uncharacterised human protein B17L was involved in Complex I function . This protein was later found to be participating as a chaperone in Complex I assembly and a mutation in this gene was identified in patients showing severe encephalopathy .
Fungi is the group of eukaryotic organisms that is best sampled in terms of fully sequenced genomes [6, 7]. The adaptation of this kingdom to a diversity of environments is reflected in a high metabolic variability that also affects the respiratory pathway [3, 8]. Indeed, the adaptation to oxygen-limited conditions or to high levels of oxidative stress during certain phases of their life cycle may have been crucial in the emergence of fermentative or pathogenic lifestyles. A recent comparative genomics study  has provided a comprehensive view of the patterns of presence and absence of OXPHOS components in 27 fungal species. Here we extend the analyses to 60 fully-sequenced fungal genomes and use a phylogenetics approach that enables us not only to obtain reliable orthology relationships but also to trace the history of duplications of OXPHOS components and related pathways during fungal evolution. In particular, we wanted to assess the role that gene duplication and functional divergence has played in the evolution of this pathway. A prediction of the gene-balance hypothesis is that independent duplications of protein complexes are likely to have deleterious effects , thereby constraining this mode of evolution in a pathway that is mostly composed of large complexes. Moreover, we wanted to test whether some loss or duplications of OXPHOS components could be associated to specific phenotypes such as virulence or adaptation to anaerobic environments. Altogether, our results show a relatively high rate of duplication events that affect 76% of the protein families surveyed. Interestingly, some of these duplications have been directly followed by processes of functional divergence, sometimes involving the recruitment of one of the duplicates to other multi-protein complexes.
Results and Discussion
Phylogenomic profiling of the OXPHOS pathway
Sequences of fungal proteins annotated as OXPHOS components were retrieved from the KEGG database  and used as queries for blastp searches against the proteins encoded in 60 fully-sequenced fungal genomes (see figure 1 and Methods section). A phylogenetic analysis was performed on each set of homologous proteins to derive a phylogenetic tree. This tree was used to establish orthology and paralogy relationships using a species-overlap algorithm that has been described earlier . This phylogeny-based approach to orthology detection, approaches more closely the original definition of orthology and reflects more appropriately the complex evolutionary relationships within protein families [13, 14]. The presence of the different components of the respiratory pathway in the species surveyed is summarized in figures 2, 3 and 4. Overall, our results agree with those reported by Lavin et. al in the 27 species that both surveys have in common . That the two approaches render so similar results, indicates that, despite using different approaches both methods have a similar stringency in the detection of OXPHOS components in this taxonomic range. In addition, our study extends the information on the distribution of components of the respiratory pathway to 33 additional species. The main advantage of our approach, however, is more qualitative than quantitative. By performing phylogenetic analyses on every protein family, we can readily obtain information on duplication events affecting components of the respiratory pathway, an important evolutionary process that was ignored in the previous study. Recognizing gene duplications is important, since this process is considered one of the main processes that drive functional innovation . Our study reveals that duplication events have affected the OXPHOS pathway extensively. Overall, we detect duplications in 76% of the families surveyed. These results were similar (duplications in 75% of the families surveyed), when more stringent cut-offs for homology detection were applied (see figures S1, S2 and S3 in the additional file 1). Such high proportion of duplications is not the result of errors in the annotation or assembly of the genomes. We controlled for this by inspecting manually every duplication case to discard dubious cases. Moreover, even when species-specific duplications in which duplicates had more than 95% identity at the nucleotide level or all duplications from the recently assembled genomes Postia placenta and Puccina graminis were not taken into account, the fraction of OXPHOS families with a duplication event remained high (71% and 74%, respectively). Zygomycota, in particular, present the highest proportion of duplicated proteins in the OXPHOS pathway. For instance, we found duplications in 60% of the genes involved in Rhizopus oryzae OXPHOS pathway. A large percentage of these duplications (82%) can be mapped specifically to the R. oryzae lineage or to the lineage preceding the separation of R. oryzae and Phycomyces blakesleeanus and thus are specific of Zygomycota species. This large amount of lineage specific duplications seems to be general in R. oryzae and P. blakesleeanus (unpublished observation from our group). An interesting possibility is that the ancestors of these organisms underwent a Whole Genome Duplication (WGD) event, similar to that described for Saccharomyces . This possibility has recently been confirmed for R. oryzae , in a comprehensive study that catalogues duplicated regions where the gene order is conserved. Consistently with our results, a duplication of nearly all subunits of the protein complexes associated with respiratory electron transport chains is detected, although our phylogeny-based approach detects additional, more ancestral, duplications that are not associated to the WGD event.
Complete loss of the OXPHOS pathway in microsporidia and two additional independent losses of Complex I coupled with alternative dehydrogenase expansions in Schizosaccharomyces and Saccharomycetales
Our results confirm earlier findings of a complete loss of the OXPHOS pathway in microsporidia  and the absence of most components of Complex I in Schizosaccharomyces and Saccharomycetales . We are able to find most of the subunits of complex I in the Taphrinomycotina species Pneumocystis carinii, suggesting that the event of gene loss occurred after the diversification of Pneumocystis and Schizosaccharomyces lineages. The apparent multiple absences of Complex I subunits, and those of other complexes, in P. carinii is probably related to a low coverage of the genome sequence for this organism. Similarly, the presence of a complete repertoire of Complex I subunits in all species in the Candida cluster and the lack of this complex in all surveyed species from the Saccharomyces/Kluyveromyces clade, situates the loss of Complex I in the latter lineage. Remarkably, the two independent losses of Complex I in the Taphrinomycotina and Saccharomyces/Kluyveromyces clades are concomitant with independent expansions of their alternative NADH dehydrogenases repertoire by virtue of gene duplications. Alternative NADH dehydrogenases bypass Complex I electron transport, oxidizing NADH without pumping of protons. The duplication of alternative NADH dehydrogenases (Figure 5) might have provided a selective advantage for yeast species using predominantly fermentative metabolism, due to adaptation to anaerobic environments. Excess of NADH causes a problem under fermentative anaerobic growth, since it prevents further oxidation of substrates due to a lack of a sufficient NAD+ pool to accept electrons. Thus, the diversification of pathways to further oxidize NADH would have been beneficial in such conditions. The loss of Complex I in the same evolutionary periods might also be related to adaptation to fermentative growth. It is unclear which of the processes preceded the other or whether both processes were concomitant. A higher taxon sampling within the Saccharomycotina and Taphrinomycotina might help to solve this issue in the future. Also coupled with Complex I loss, and in line with adaptations to anaerobic environments in the abovementioned lineages, we observe the loss of alternative oxidases.
Duplications in Alternative Oxidases are not necessarily coupled with a pathogenic life-style
Alternative oxidases catalyze the cyanide-resistant alternative pathway of mitochondrial respiration in some fungi, plants and several protists. This pathway directly transfers electrons from the ubiquinone pool to oxygen, thereby bypassing complex III and cytochrome c oxidase . Alternative oxidases are common in yeasts but limited almost exclusively to non-fermentative and crabtree-negative yeasts. Alternative oxidases participate in energy production but also in antioxidant defense of cells. It has been shown that alternative oxidases represent an important factor for the survival of pathogenic fungi inside macrophages . Considering this, it could be postulated that the duplication of these enzymes might have played a role in the emergence of pathogenesis in several mammal fungal pathogens. In our survey we detect several copies of alternative oxidases in 13 species. Some of these duplications seem to have occurred quite recently in their respective lineages, such as the duplication that lead to AOX1 and AOX2 (orf19.4774 and orf19.4773 in C. albicans) involving some Candida species, which can be mapped before the speciation of C. tropicalis, C. dubliniensis and C. albicans. Although many of these duplications do affect pathogenic genera such as Candida and Aspergillus, there are notable exceptions such as the intra-specific duplications found in the generally non-pathogenic species Yarrowia lypolytica or Coprinus cinereus. Conversely, we find pathogenic species such as Histoplasma capsulatum or Cryptococcus neoformans that have been shown to survive in macrophages  and nevertheless present a single alternative oxidase. Taken all together, our results suggest that a single copy of alternative oxidase gene is sufficient to protect fungal pathogens against macrophages and rather points to alternative selective advantages for the duplication of this gene. Conversely, alternative adaptations might be behind the emergence of the ability to survive inside macrophages in certain lineages. For instance, the presence of a polysaccharide capsule in Cryptococcus has been shown to confer resistance to oxidative stress .
Extensive duplication followed by functional divergence in the fungal OXPHOS pathway
According to the gene balance hypothesis , the duplication of genes that encode for subunits of multi-protein complexes should have a higher chance of being deleterious due to dosage effects. As a result, one would expect to find few duplication events in the OXPHOS system, as this is mainly formed by intricate complexes. Contrary to that expectation, we find numerous cases of duplications in OXPHOS proteins, which overall affect 66 (76%) of the proteins surveyed. These have occurred at different moments in fungal evolution. At least for the genes duplicated during the Whole Genome Duplication event (WGD) occurred in the yeast lineage about 80 Myr ago, the rate of gene loss of duplicated OXPHOS genes is not higher than the overall rate for Saccharomyces cerevisiae. Indeed, our study finds six yeast proteins (all complex II subunits, Qcr6p and Cox5p), whose duplication is mapped to the WGD event. These represent 7% of the OXPHOS proteins, meaning that for 93% of the nuclear OXPHOS proteins supposedly duplicated in WGD were subsequently lost, a rate of gene loss that is roughly similar to the 88% estimated for the whole S. cerevisiae genome . It must be noted that one of the duplications affected the whole Complex II, meaning that the duplication was conservative in terms of stoichiometry of the different subunits. However, one of the duplicated subunits has been subsequently lost in the Saccharomyces sensu stricto species, suggesting the four duplicates do not form an alternative complex II. A possible reconciliation between the extensive rate of gene duplications and the gene balance hypothesis is that functional divergence directly followed the duplication event, thereby facilitating the retention of both duplicates . Differences in the expression patterns of some of the WGD duplicates, point to a functional specialization of each duplicate. For instance, Cox5 (YNL052W) is expressed during aerobic growth whereas its paralog (YIL111W) is expressed under anaerobic growth . Similarly, the duplicate of the SDHA complex II subunit (YJL045W) is specifically expressed during the diauxic shift. Several other observations suggest that functional divergence processes have been common after duplication of OXPHOS protein families (see below).
Evolutionary cross-talk between the OXPHOS complexes and other multi-protein complexes
Several instances of paralogy relationships between complex I subunits and other mitochondrial multi-protein complexes have been previously reported . This is the case for the NDUFA11 subunit, which is paralogous to the Tim17/22 family as well as that of NI8M (NDUFA2) and NUZM, which are paralogous to L43 and L2 subunits of the mitochondrial ribosome. It has been suggested that OXPHOS proteins with paralogs in other complexes would play a structural role rather than being involved in proton or electron transport, since ribosomes and the import machinery do not display those functions . Similarly, we find several instances of paralogs of OXPHOS subunits that play a role in other complexes. Interestingly, another evolutionary connection between OXPHOS and the mitochondrial import machinery (MIM) is evidenced by the fact that the MIM subunit TIM18 (YOR297C) is a paralog of the Complex II subunit SDHD (YDR178W). Yet another paralog of the same Complex II subunit, which originated from a more recent duplication in the Saccharomycotina lineage (YLR164W) encodes for a mitochondrial inner membrane protein of yet unknown function. Paralogies to the protein import system in the mitochondrion extend to the two subunits of the Mitochondrial processing peptidase (MPP), an essential processing enzyme that cleaves the N-terminal targeting sequences from mitochondrially imported proteins . Indeed the large and small subunits MAS1 (YLR163C) and MAS2 (YHR024C) are homologous to QCR1 (YBL045C) and QCR2 (YPR191W) subunits of Complex III. Paralogy relationships to other multi-protein complexes extend beyond mitochondria. Indeed, several paralogs of Complex V subunits have been described as components of complexes from other cell compartments. For instance, the alpha and beta subunits of the F1 sector of the mitochondrial ATP synthase (YBL099W, YJR121W) are paralogous to the A and B subunits of the vacuolar ATP synthase (YDL185W, YBR127C). Vacuolar ATP synthases are found in the membranes of a large number of organelles which include endosomes, lysosomes and secretory vesicles. This duplication, however is not specific to fungi, since both paralogous groups have representatives in Arabidopsis thaliana and Homo sapiens (see phylogenetic trees additional file 1: figure S4), which indicates that the duplication preceded the diversification of plants and opisthokonts.
Yeast ACPM is possibly not a complex I remnant but a Saccharomycotina-specific paralog of complex I Acyl-carrier protein
Although previously identified as a remnant Complex I subunit in the Complex I devoid organism S. cerevisiae , our current phylogenomic analysis suggest that this protein might actually be a paralog originated from an ancient duplication that is specifically conserved in Saccharomycotina (Figure 6). This paralogy relationship is supported by approximate Likelihood Ratio Tests (aLRT) analyses of the duplication node (0.79, shown in the figure) as well as by bootstrap analyses (74% bootstrap support, not shown in the figure). This finding clarifies apparent inconsistencies in the function of Complex I acyl-carrier protein and the isolated ACPM protein (YKL192C). Indeed S. cerevisiae ACPM protein has been found to participate in the synthesis of octanoic acid, a precursor of lipoic acid , whereas the acyl carrier Complex I subunit in Neurospora crassa seems not to participate in this process . The fact that a mammalian homolog of this protein family has also been identified as a Complex I subunit in bovine mitochondria , suggests that association with Complex I is an ancestral feature of the family. Taken together, these results indicate that a process of functional divergence might have occurred after the duplication event diverting the new duplicate for specializing in the synthesis of octanoid acid. This specialization is presumably present in all Saccharomycotina species including the Candida group, which additionally possesses the true acyl carrier Complex I subunit.
Altogether our results shed light on how processes of gene loss, duplication and functional divergence have shaped the core of the respiratory pathway in fungi. Although most fungal organisms present a similar overall composition in terms of respiratory complexes, extensive differences in what particular units have been lost or duplicated in each complex, might help explaining differences found at the physiological level. This continuous evolution of OXPHOS components seems to be common in other groups of organisms [4, 29, 30], emphasizing the plasticity of this central energetic pathway.
Proteins encoded in 60 fully-sequenced fungal genomes were downloaded from several databases (figure 1). For consistency, we used in our analysis the species names as provided by the database source. Some of these species have been renamed and the corresponding new names and synonyms are listed in the additional file 1 (Additional table S1). Additionally, genomes from Homo sapiens and Arabidopsis thaliana were downloaded from ensembl http://www.ensembl.org. The final database comprises 626,834 unique protein sequences.
Reconstruction of the presence/absence matrix
Fungal proteins annotated as being part of the OXPHOS pathway were downloaded from the KEGG database (map 00190) . In addition, 6 proteins that were identified in the literature as belonging to complex I but were not present in the KEGG database were downloaded from UniProt and included in the analyses (NI9M, NURM, NUWM, NUXM and NUZM). The resulting 85 proteins were used to perform a blast search against a database of fungal proteins encoded in 60 fungal genomes (figure 1). Low complexity filters were used in the blast search. To detect homology, we used the same parameters that have been used previously in the same taxonomic range . In brief, only significant hits (E-val < 10-3) that aligned with a continuous region covering more than one third of the query sequence were selected. Note that the use of low complexity filters in the blast can reduce significantly the length of continuous regions of homology. Sets of homologous sequences were aligned and used to reconstruct a Maximum Likelihood tree from which orthology relationships were inferred (see below). These orthology relationships were used to build a presence/absence matrix (figures 2, 3 and 4) in which for each OXPHOS component the species with a corresponding ortholog are indicated. Putative absences in the matrix were double-checked by tBlastN  searches against the corresponding genome sequence and Blast searches from family members of more related species. These hits were checked manually and whenever they were considered orthologous to the already identified members they were added to the list.
We used a similar pipeline to that described in . Sets of homologous proteins were aligned using MUSCLE 3.6  with default parameters. Positions in the alignment with gaps in more than 10% of the sequences were trimmed with trimAl . Finally, PhyML aLRT version [36, 37] was used to derive Maximum Likelihood (ML) trees. Four different evolutionary models were used for each seed sequence (JTT, WAG, Blosum62 and VT). In all cases, a discrete gamma-distribution model with four rate categories plus invariant positions was used, estimating the gamma parameter and the fraction of invariant positions from the data. The evolutionary model best fitting the data was determined by comparing the likelihood of the used models according to the AIC criterion . Orthology and paralogy relationships among members of a family were inferred from the analysis of their corresponding phylogenetic trees, using a previously described algorithm that has been described before and has been shown to be accurate [7, 12]. Phylogeny-based methods are considered to better reflect the actual complexity of orthology relationships than pair-wise methods such as best-bidirectional hits . All phylogenetic trees are provided in the additional file 1: figure S4 as well as a list of proteins used as a seed in our analyses (Additional file 1: table S2). All duplications where manually checked to discard possible cases of spurious duplications. This was done by manually inspecting the alignments and the nucleotide sequences of the relevant duplicates. Moreover, the corresponding genome browsers or assembly data were searched to analyze the sequence context of the duplicates. Highly similar sequences in which one is only partially sequenced or in a small contig can be taken as possible source of errors. These cases were discarded. Total counts of duplications were also computed discarding the fraction of duplications that is expected to be more sensitive to error annotation: lineage-specific duplications with highly similar duplicates and duplications found in recently assembled genomes (see main text).
Saraste M: Oxidative phosphorylation at the fin de siecle. Science. 1999, 283 (5407): 1488-1493. 10.1126/science.283.5407.1488.
Gabaldón T, Huynen MA: Reconstruction of the proto-mitochondrial metabolism. Science. 2003, 301 (5633): 609-10.1126/science.1085463.
Gabaldón T, Huynen MA: Shaping the mitochondrial proteome. Biochim Biophys Acta. 2004, 1659 (2-3): 212-220. 10.1016/j.bbabio.2004.07.011.
Gabaldón T, Rainey D, Huynen MA: Tracing the Evolution of a Large Protein Complex in the Eukaryotes, NADH:Ubiquinone Oxidoreductase (Complex I). J Mol Biol. 2005, 348 (4): 857-870. 10.1016/j.jmb.2005.02.067.
Ogilvie I, Kennaway NG, Shoubridge EA: A molecular chaperone for mitochondrial complex I assembly is mutated in a progressive encephalopathy. J Clin Invest. 2005, 115 (10): 2784-2792. 10.1172/JCI26020.
Galagan JE, Henn MR, Ma LJ, Cuomo CA, Birren B: Genomics of the fungal kingdom: insights into eukaryotic biology. Genome Res. 2005, 15 (12): 1620-1631. 10.1101/gr.3767105.
Marcet-Houben M, Gabaldón T: The tree versus the forest: the fungal tree of life and the topological diversity within the yeast phylome. PLoS ONE. 2009, 4 (2): e4357-10.1371/journal.pone.0004357.
Bullerwell CE, Lang BF: Fungal evolution: the case of the vanishing mitochondrion. Curr Opin Microbiol. 2005, 8 (4): 362-369. 10.1016/j.mib.2005.06.009.
Lavin JL, Oguiza JA, Ramirez L, Pisabarro AG: Comparative genomics of the oxidative phosphorylation system in fungi. Fungal Genet Biol. 2008, 45 (9): 1248-1256. 10.1016/j.fgb.2008.06.005.
Papp B, Pal C, Hurst LD: Dosage sensitivity and the evolution of gene families in yeast. Nature. 2003, 424 (6945): 194-197. 10.1038/nature01771.
Okuda S, Yamada T, Hamajima M, Itoh M, Katayama T, Bork P, Goto S, Kanehisa M: KEGG Atlas mapping for global analysis of metabolic pathways. Nucleic Acids Res. 2008, W423-426. 10.1093/nar/gkn282. 36 Web Server
Huerta-Cepas J, Dopazo H, Dopazo J, Gabaldón T: The human phylome. Genome Biol. 2007, 8 (6): R109-
Fitch WM: Distinguishing homologous from analogous proteins. Syst Zool. 1970, 19 (2): 99-113. 10.2307/2412448.
Gabaldón T: Large-scale assignment of orthology: back to phylogenetics?. Genome Biol. 2008, 9 (10): 235-10.1186/gb-2008-9-10-235.
Ohno S: Evolution by gene duplication. 1970, London: Allen & Unwin
Kellis M, Birren BW, Lander ES: Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature. 2004, 428 (6983): 617-624. 10.1038/nature02424.
Ma LJ, Ibrahim AS, Skory C, Grabherr MG, Burger G, Butler M, Elias M, Idnurm A, Lang BF, Sone T, et al: Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication. PLoS Genet. 2009, 5 (7): e1000549-10.1371/journal.pgen.1000549.
Katinka MD, Duprat S, Cornillot E, Metenier G, Thomarat F, Prensier G, Barbe V, Peyretaillade E, Brottier P, Wincker P, et al: Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature. 2001, 414 (6862): 450-453. 10.1038/35106579.
Veiga A, Arrabaca JD, Loureiro-Dias MC: Cyanide-resistant respiration, a very frequent metabolic pathway in yeasts. FEMS Yeast Res. 2003, 3 (3): 239-245. 10.1016/S1567-1356(03)00036-9.
Magnani T, Soriani FM, Martins VP, Nascimento AM, Tudella VG, Curti C, Uyemura SA: Cloning and functional expression of the mitochondrial alternative oxidase of Aspergillus fumigatus and its induction by oxidative stress. FEMS Microbiol Lett. 2007, 271 (2): 230-238. 10.1111/j.1574-6968.2007.00716.x.
Johnson CH, Prigge JT, Warren AD, McEwen JE: Characterization of an alternative oxidase activity of Histoplasma capsulatum. Yeast. 2003, 20 (5): 381-388. 10.1002/yea.968.
Zaragoza O, Chrisman CJ, Castelli MV, Frases S, Cuenca-Estrella M, Rodriguez-Tudela JL, Casadevall A: Capsule enlargement in Cryptococcus neoformans confers resistance to oxidative stress suggesting a mechanism for intracellular survival. Cell Microbiol. 2008, 10 (10): 2043-2057. 10.1111/j.1462-5822.2008.01186.x.
Lynch M, Katju V: The altered evolutionary trajectories of gene duplicates. Trends Genet. 2004, 20 (11): 544-549. 10.1016/j.tig.2004.09.001.
Hodge MR, Kim G, Singh K, Cumsky MG: Inverse regulation of the yeast COX5 genes by oxygen and heme. Mol Cell Biol. 1989, 9 (5): 1958-1964.
Gakh O, Cavadini P, Isaya G: Mitochondrial processing peptidases. Biochim Biophys Acta. 2002, 1592 (1): 63-77. 10.1016/S0167-4889(02)00265-3.
Brody S, Oh C, Hoja U, Schweizer E: Mitochondrial acyl carrier protein is involved in lipoic acid synthesis in Saccharomyces cerevisiae. FEBS Lett. 1997, 408 (2): 217-220. 10.1016/S0014-5793(97)00428-6.
Schneider R, Massow M, Lisowsky T, Weiss H: Different respiratory-defective phenotypes of Neurospora crassa and Saccharomyces cerevisiae after inactivation of the gene encoding the mitochondrial acyl carrier protein. Curr Genet. 1995, 29 (1): 10-17. 10.1007/BF00313188.
Runswick MJ, Fearnley IM, Skehel JM, Walker JE: Presence of an acyl carrier protein in NADH:ubiquinone oxidoreductase from bovine heart mitochondria. FEBS Lett. 1991, 286 (1-2): 121-124. 10.1016/0014-5793(91)80955-3.
De Grassi A, Lanave C, Saccone C: Genome duplication and gene-family evolution: the case of three OXPHOS gene families. Gene. 2008, 421 (1-2): 1-6. 10.1016/j.gene.2008.05.011.
Saccone C, Lanave C, De Grassi A: Metazoan OXPHOS gene families: evolutionary forces at the level of mitochondrial and nuclear genomes. Biochim Biophys Acta. 2006, 1757 (9-10): 1171-1178. 10.1016/j.bbabio.2006.04.021.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, et al: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, D480-484. 36 Database
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Huerta-Cepas J, Bueno A, Dopazo J, Gabaldón T: PhylomeDB: a database for genome-wide collections of gene phylogenies. Nucleic Acids Res. 2008, D491-496. 36 Database
Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004, 5 (1): 113-10.1186/1471-2105-5-113.
Capella-Gutíerrez S, Silla-Martínez JM, Gabaldón T: trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009, 25 (15): 1972-3. 10.1093/bioinformatics/btp348.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55 (4): 539-552. 10.1080/10635150600755453.
Akaike H: Information theory and extension of the maximum likelihood principle. Proceedings of the 2nd international symposium on information theory: 1973; Budapest, Hungary. 1973, 267-281.
Christie KR, Weng S, Balakrishnan R, Costanzo MC, Dolinski K, Dwight SS, Engel SR, Feierbach B, Fisk DG, Hirschman JE, et al: Saccharomyces Genome Database (SGD) provides tools to identify and analyze sequences from Saccharomyces cerevisiae and related sequences from other organisms. Nucleic Acids Res. 2004, D311-314. 10.1093/nar/gkh033. 32 Database
MMH and TG are supported by grants from the ERA-NET pathogenomics network funded by the Spanish Ministry of Science (GEN06-27784).
TG conceived of the study and coordinated the analyses. MMH and GM performed the analyses. TG and MMH wrote the manuscript. All authors read and approved the manuscript.