- Open Access
Origin and evolution of the peroxisomal proteome
Biology Direct volume 1, Article number: 8 (2006)
Peroxisomes are ubiquitous eukaryotic organelles involved in various oxidative reactions. Their enzymatic content varies between species, but the presence of common protein import and organelle biogenesis systems support a single evolutionary origin. The precise scenario for this origin remains however to be established. The ability of peroxisomes to divide and import proteins post-translationally, just like mitochondria and chloroplasts, supports an endosymbiotic origin. However, this view has been challenged by recent discoveries that mutant, peroxisome-less cells restore peroxisomes upon introduction of the wild-type gene, and that peroxisomes are formed from the Endoplasmic Reticulum. The lack of a peroxisomal genome precludes the use of classical analyses, as those performed with mitochondria or chloroplasts, to settle the debate. We therefore conducted large-scale phylogenetic analyses of the yeast and rat peroxisomal proteomes.
Our results show that most peroxisomal proteins (39–58%) are of eukaryotic origin, comprising all proteins involved in organelle biogenesis or maintenance. A significant fraction (13–18%), consisting mainly of enzymes, has an alpha-proteobacterial origin and appears to be the result of the recruitment of proteins originally targeted to mitochondria. Consistent with the findings that peroxisomes are formed in the Endoplasmic Reticulum, we find that the most universally conserved Peroxisome biogenesis and maintenance proteins are homologous to proteins from the Endoplasmic Reticulum Assisted Decay pathway.
Altogether our results indicate that the peroxisome does not have an endosymbiotic origin and that its proteins were recruited from pools existing within the primitive eukaryote. Moreover the reconstruction of primitive peroxisomal proteomes suggests that ontogenetically as well as phylogenetically, peroxisomes stem from the Endoplasmic Reticulum.
This article was reviewed by Arcady Mushegian, Gáspár Jékely and John Logsdon
Open peer review
Reviewed by Arcady Mushegian, Gáspar Jékely and John Logsdon. For the full reviews, please go to the Reviewers' comments section.
Peroxisomes were first isolated from liver and biochemically characterized by the group of de Duve . Later it became clear that these organelles can differ substantially between species with respect to their enzyme content. The conversion of fatty acids into carbohydrates through the glyoxylate cycle is the hallmark of glyoxysomes present in plants, protozoa and yeasts. Part of the glycolysis is compartmentalized in the glycosomes of Trypanosomatids. Photorespiration is typical for plant peroxisomes while peroxisomes of various yeasts can oxidize alkanes or methanol. Despite this diversity all these organelles belong to the same microbody family. This became clear with the discovery that they share targeting codes (PTS1 and PTS2) for the import of proteins and with the identification of similar sets of proteins responsible for their biogenesis and maintenance . Although the unity within the microbody family has thus firmly been established, their evolutionary origin remains a matter of debate . Strong arguments support the view of peroxisomes as autonomous organelles with an endosymbiotic origin: i) matrix enzymes are synthesized on free polyribosomes and post-translationally imported into the organelles, ii) peroxisomes have their own protein import machinery, like mitochondria and chloroplasts, and iii) peroxisomes have been shown to divide .
Recent discoveries, however, have challenged this view. First, after several generations the lacking of peroxisomes in some mutants is reversible upon the introduction of the wild-type gene . Second, it has been observed that new peroxisomes originate from the ER . These observations are at odds with the autonomy of peroxisomes and therefore weaken the case for an endosymbiotic origin. Here we address the issue of peroxisomal evolution by phylogenetic analysis of peroxisomal proteins. To this end we collected an exhaustive set of proteins with an experimentally determined peroxisomal location in the yeast Saccharomyces cerevisiae and the rodent Rattus norvegicus, and performed phylogenetic analyses to investigate whether the peroxisomal proteome contains a significant evolutionary signal just as has been shown for mitochondria [7, 8].
Results and discussion
From databases and experimental literature we collected 62 yeast and 51 rat proteins with a peroxisomal location or function (Table 1). Since our lists include proteins from various large-scale proteomics analyses [9–11], as well as from individual studies under various conditions, we consider them to be representative samples of peroxisomal proteomes. Phylogenies (see materials and methods) of peroxisomal proteins were reconstructed to determine their origin. We consider a protein to be of eukaryotic origin when it has no homologs in prokaryotes, or when the prokaryotic branches within the tree are mono-phyletic as in Figure 1a. In the latter case the protein is classified as of "ancient origin" in Table 1, even though one could argue that in the case of PEX1 the protein resulted from a gene duplication at the origin of the eukaryotes. Although in the case of PEX1 the relatively short branch length of CDC48 suggests that CDC48 is the "ancestral protein" and PEX1 is derived, in general such a distinction is not easy to make and in this analysis we did not distinguish between genes that are, or are not duplicated at the origin of the eukaryotes. A protein is considered of bacterial or archaeal origin when it clusters "within" a prokaryotic branch, implying horizontal transfer between the taxa (Figure 1b). Unresolved cases imply the existence of homology to prokaryotic sequences without a tree that specifically supports a bacterial or archaeal origin. For the families with resolved phylogenies we observed a clear dichotomy in terms of evolutionary origin and functional roles: all proteins with a specific bacterial origin have enzymatic functions while most proteins (90%) with eukaryotic origin are functioning in peroxisome organization and biogenesis. Like in the proteins with bacterial ancestry also among the proteins with bacterial homologs for which we cannot establish bacterial ancestry (the unresolved cases) a clear preponderance (85%) of enzymes can be observed (Table 1).
Peroxisomal proteins of eukaryotic origin and an evolutionary link with the E.R
The largest fraction of peroxisomal proteins is of eukaryotic origin: 58.1% of the yeast proteome, 39.2% of the rat proteome (Figure 1c). These include the so-called Pex proteins that are involved in peroxisomal biogenesis and maintenance that are most consistently present in all microbodies, underlining their essential role. Interestingly, five of the six most ancient Pex proteins (see below) show homology with the ERAD (Endoplasmic Reticulum Associated Decay) system, which pulls proteins from the ER membrane and ubiquitinylates them in preparation for degradation in the proteasome  (Figure 2). Pex1 and Pex6, AAA cassette containing proteins, have evolved from Cdc48/p97  (Figure 1a), a protein central to the ERAD pathway which is also involved in Golgi vesicle fusion and spindle body disassembly after mitosis; Pex2 and Pex10, ubiquitin ligase domain (RING domain) containing proteins, contain homology to the ERAD ubiquitin ligase Hrd1; the TPR repeats of Pex5 are homologous to the SEL1 repeats of the Hrd1 interacting protein Hrd3; Pex4 contains an E2 ubiquitin conjugating enzyme domain and is homologous to the ERAD ubiquitin conjugating enzymes Ubc1, Ubc6 and Ubc7. In the cases of PEX2/10, PEX5 and PEX4 the levels of sequence identity between the shared domains and the short regions of homology preclude the reconstruction of reliable phylogenies to argue that these proteins have descended from a protein involved in ERAD, as it is the case for Cdc48/p97-PEX1/6. Here it is the number of homologous relations between ERAD and the most ancient PEX proteins that hint at an evolutionary relation. Although there are some systems known that use a TPR repeat protein together with a protein containing an E2 ubiquitin conjugating enzyme and a protein with a RING domain, like the Anaphase Promoting Complex/Cyclosome , to our knowledge there is no system other than ERAD that uses those domains together with an AAA+ ATPase. Nevertheless, we cannot exclude that PEX1, PEX2/10, PEX5 and PEX4 do not originate from a single molecular system like ERAD, specifically because the TPR repeat in HRD3 is classified in a different class of TPR repeats than the TPR repeat of PEX5 (Figure 2).
The similarities in amino acid sequence between ERAD and the most ancient PEX proteins extend into similarities in function and sub-cellular location (Figure 2). Pex1 and Pex6 (both AAA containing proteins) are needed to extract the cycling PTS1 receptor Pex5 from the peroxisomal membrane to facilitate a new cycle of Pex5-mediated protein import . Ubiquitinylation of Pex5 is part of this process. In both cases, the ERAD and the peroxisomal AAA proteins operate in the cytoplasm and are recruited to the membrane by organelle-specific anchor proteins: Cdc48/p97 to the ER membrane by VIMP, Pex1 and Pex6 to the peroxisomal membrane by Pex15 (in yeast) and Pex26 (in mammals) . This resemblance in ancient proteins with similar functions and the link with the universal endomembrane compartment of the eukaryotic cell suggest that the peroxisome is an invention that took place within the eukaryotic lineage itself. Also Erdmann and Schliebs  have recently linked the homology between AAA+ domains of ERAD and PEX1, and the presence of E2 and E3 domains involved in ubiquitinylation in the PEX proteins, to a mechanism of protein import in the Peroxisomal matrix that would be similar to ERAD, without proposing a direct evolutionary descent of Peroxisomal import from ERAD however.
For the other PEX proteins we did not find indications that they were also recruited from pre-existing cellular systems. Their distribution and phylogenies do suggest that they originate from separate events post-dating the origin of the five of the six core PEX proteins from ERAD.
We have visualized the retargeting during evolution of peroxisomal proteins from various cellular locations in a cartoon (Figure 3). The group of proteins of eukaryotic origin also contains certain household proteins with dual or plural functions with respect to organelles. The ER located or associated proteins Erg1, Erg6, Emp24, Rho1 and the multipurpose dynamin Vps1 have also been implicated in peroxisomal functions [19, 20].
Recruitment to the peroxisome of proteins of alpha-proteobacterial origin
Remarkably, the second largest fraction of proteins, 17–18%, has an alpha-proteobacterial origin (Figure 1c). This is similar to what has been found for mitochondria [7, 8], and, at first sight appears to be at odds with a eukaryotic origin of the peroxisome. There are strong indications that these proteins have been retargeted from the mitochondria (Figure 3, scenario I), rather than having evolved directly from an independent endosymbiont, an observation that is consistent with the high degree of retargeting observed for proteins derived from the proto-mitochondrion in general . Six of the eight S. cerevisiae peroxisomal proteins of alpha-proteobacterial origin are closely related to mitochondrial proteins. Thioesterase (Tes1p) is located in both the peroxisome and mitochondrion of S. cerevisiae . In other cases the orthologs or paralogs of a peroxisomal protein are mitochondrial: i), the peroxisomal glycerole-P dehydrogenase Gpd1p has a paralog in yeast (Gpd2p) with a cytoplasmic and mitochondrial localization ; ii) the peroxisomal Fat2p is paralogous to the mitochondrial long-chain fatty acid CoA ligases iii), the orthologous group consisting of Eci1p, Dci1p and 3,2-transenoyl CoA isomerase is peroxisomal in yeast and human, has a mitochondrial paralog in mammals ; and iv), the nudix phosphatase family (Npy1p) of which the yeast, human and plant orthologs are peroxisomal has a paralogous group in metazoa that is mitochondrial according to GFP-fusion studies in mouse  and to Mitoprot  (p = 0.97). The phylogenetic tree (Figure 1b) indicates a single origin from the alpha-proteobacteria of both mitochondrial and peroxisomal proteins of this family. The two remaining cases of yeast peroxisomal proteins of alpha-proteobacterial origin are Fox2p and Pcd1p. For these no homologs with experimental evidence of mitochondrial location were found, although Pcd1p does have a bona-fide mitochondrial targeting signal (P = 0.97 in Mitoprot).
With respect to the rat peroxisome, there are two proteins of alpha-proteobacterial descent that do not have orthologs in the yeast peroxisome. One of these presents cases of dual targeting: some isoforms of peroxisomal bile acid thioestherase BAAT have been detected in mitochondria and the cytoplasm in human liver .
Recruitment to the peroxisome of mitochondrial proteins of other origins
There are also peroxisomal proteins with homologs in the mitochondrion that do not have a (detectable) alpha-proteobacterial origin: Idp3p, Cta1p, Faa1p, Cit2p, Fis1p and Faa2p [21, 27, 28] (Figure 3, Scenario II). In contrast to proteins of alpha-proteobacterial origin, here one cannot simply argue that the mitochondrial localization preceded the peroxisomal one. At least for one of these proteins, Cit2p, a peroxisomal protein from the citrate synthase family, a phylogenetic analysis reveals its ancestral location. The other two members of this family in S. cerevisiae, Cit1p and Cit3p, are mitochondrial and so are their homologs in Homo sapiens, Arabidopsis thaliana and Caenorhabditis elegans. The phylogeny of this family in fungi indicates that Cit1p and Cit2p originate from a recent gene duplication, after which Cit2p lost its mitochondrial targeting signal (Figure 4), indicating that the peroxisomal location is secondary. That the retargeting of proteins between mitochondria and peroxisomes frequently happens during evolution is also indicated by the case of alanine:glyoxylate aminotransferase (AGT), whose peroxisomal or mitochondrial location is species-dependent and related to diet in mammals . In humans, where AGT is peroxisomal, a single point mutation miss-localizes the protein to the mitochondrion, leading to the hereditary kidney stone disease: primary hyperoxaluria type 1 (PH1).
There are peroxisomal rat proteins, like dihydroxyacetone phosphate acyl transferase(DAPT) and alkyl-dihydroxyacetonephosphate synthase (gi-12002203) whose phylogenetic trees suggest an ancestry from within the actinomycetales while the Alanine-Glyoxylate aminotransferase (AGT) appears derived from the cyanobacteria. We do not have an obvious evolutionary scenario for the origin for such proteins with a bacterial but not alpha-proteobacterial ancestry. In any case, the finding of peroxisomal proteins with such diverse origins underscores the ease at which the peroxisomal proteome can recruit new proteins.
Reconstruction of ancestral states of the peroxisomal proteome
To investigate the order of protein recruitment to the peroxisome we reconstructed the evolution of the peroxisomal proteome based on the absence/presence of genes among sequenced genomes and assuming a parsimonious scenario (Figure 5). First we reconstructed the minimal peroxisome of the opisthokont, the common ancestor of metazoa and fungi, by including proteins present in both yeast and rat peroxisomal proteomes or proteins that are present in only one of the two proteomes but whose orthologs in plants have a (putative) peroxisomal location in the Araperox database . In addition, we reconstructed the protein content of the common ancestor of all known peroxisomes, glycosomes and glyoxysomes from proteins that, besides being present in the opisthokont peroxisome, are present in genomes from plants and kinetoplastida (Trypanosoma brucei and Leishmania major). This core-set comprises six PEX proteins (Pex1p, Pex2p, Pex4p, Pex5p, Pex10p, Pex14p) and four proteins involved in fatty acid metabolism and transport (Pxa1p/Pxa2p, Fox2p, Faa2p). We also included the peroxisomal hallmark protein catalase (Cta1p), even though it is absent from most glycosomes and kinetoplastidial genomes because it is found in the glycosomes of the non-pathogenic trypanosomatid Crithidia . Similarly Fox1p, which catalyzes the first step of long-chain fatty acid beta-oxidation, was included despite its absence from kinetoplastida, because the concomitant loss from peroxisomes of Fox1 (the enzyme generating H2O2) and catalase (the enzyme detoxifying H2O2) has been observed in species such as Neurospora crassa .
Although the specific functional role of many PEX proteins remains to be established, and it is therewith hard to asses whether e.g. the reconstructed ancestral opisthokont PEX proteins are functionally coherent and complete, at least the sub-set present in the ancestral eukaryotic peroxisome appears functionally coherent. All of the six universal PEX proteins are specifically involved in the PEX5 pathway for the import of proteins into the peroxisome.
The earliest tractable function of peroxisomes appears herewith to be the beta-oxidation of fatty acids. This pathway already contains at least one protein of alpha-proteobacterial descent (Fox2p), indicating that the presence of long-chain fatty acid beta-oxidation in the peroxisome followed the endosymbiosis of mitochondria. The proteins with detectable origin in the ancestral peroxisome that are not involved in beta-oxidation are all of eukaryotic origin. Most of the present-day species variability is found in the enzymes housed in peroxisomes, a significant fraction of which has an alpha-proteobacterial origin and has entered the primitive eukaryote with the mitochondrial ancestor . Note that the recruitment of proteins with an endosymbiotic origin to peroxisomes is not an exceptional event. Nine proteins in the glycosomes of the kinetoplastida T. brucei and Leishmania mexicana are derived from chloroplasts from which they can be traced back to the cyanobacteria . Somehow it seems rather easy to (re)locate proteins to microbodies which may be related to the simplicity of the PTS1 targeting code. This 'grab what you can get' principle may have contributed to the observed versatility and species variability.
The phylogenetic analysis of the rat and yeast peroxisomal proteomes reveals that the largest fraction of peroxisomal proteins originated within the eukaryotic lineage and that the significant fraction of peroxisomal proteins which stems from the alpha-proteobacteria is likely the result of a secondary retargeting from the mitochondrion. The most widespread and ancient set of peroxisomal proteins is mainly composed of eukaryotic proteins involved in peroxisome biogenesis and organization. Most of these core proteins are evolutionarily related to the Endoplasmic Reticulum Assisted Decay pathway, suggesting an evolutionary origin of the peroxisomes from the endomembrane system. While this manuscript was under review a common evolutionary origin of the Peroxisome and the ER was also proposed by Schluter and coworkers  based on the homologies in Figure 2 between ERAD and Pex proteins (Figure 2) while this homology has also been observed by Erdman and Schliebs . In the former analysis, full length homologs with Bacterial proteins were not included and the authors could not exclude that such proteins were indeed donated by an early symbiont. The retargeting of enzymes documented in this paper solves the paradox of a eukaryotic organelle with bacterial enzymes. Recent experimental work indicates that some peroxisomal proteins first enter the ER thereby capturing part of the ER membrane for subsequent formation of the organelle . These observations are consistent with our findings that the oldest PEX proteins are homologous to proteins of the ERAD pathway, suggesting that evolutionarily as well as ontogenetically peroxisomes are in fact offshoots from the ER.
Manually curated sets of 62 S. cerevisiae and 50 R. norvegicus proteins with experimental evidence of peroxisomal location were compiled from the literature [9–11, 19] and from the Saccharomyces Genome  and Swiss-Prot  databases. For the purpose of this paper we consider a protein to be peroxisomal when it permanently resides in the peroxisomal matrix or membrane, or when it is a cytoplasmic protein but has a dedicated function in peroxisomal protein import and/or biogenesis.
Protein sequences encoded by 144 publicly available complete genomes were obtained from Swissprot , except for Plasmodium falciparum, Schizosaccharomyce pombe, Candida albicans, Encephalitozoon cuniculi (Genbank, ), Homo sapiens, Rattus norvegicus and Mus musculus (EBI, ).
For every yeast and rat peroxisomal protein, homologous sequences (E < 0.01) were retrieved using Smith-Waterman comparisons against the aforementioned 144 complete proteomes. Only sequences that aligned with at least one third of the query sequence were selected. Sequences were aligned using MUSCLE . Neighbour Joining (NJ) trees were made using Kimura distances as implemented in ClustalW . Positions with gaps were excluded and 1000 bootstrap iterations were performed. Maximum Likelihood (ML) trees were derived using PhyML v2.1b1 , with a four rate gamma-distribution model, before and after excluding from the alignment positions with gaps in 10% or more of the sequences. In all cases NJ and ML trees were manually examined to search for consistent patterns indicating the origin of the peroxisomal proteins. Trees in which eukaryotic proteins clustered together, within the Bacteria or the Archaea and with a specific prokaryotic out-group were classified as having that phylogenetic origin (e.g. Figure 1b). Trees were only regarded as resolved when both the NJ tree and the ML tree agreed to the level of resolution required, e.g. a specific bacterial group as a sister clade of the peroxisomal group of proteins, or when at least the ML tree had the level of resolution required while the NJ tree did not point to another origin of a protein.
Reconstruction of yeast, rat peroxisomal metabolisms and their ancestral states
Annotated biochemical and cellular functions of the yeast and rat peroxisomal proteins were mapped onto metabolic KEGG maps  and are represented in Figure 5, indicating their phylogenetic origin by a color-code. Proteins known or predicted to be membrane-associated are depicted close to the membrane. The minimal ancestral opisthokont peroxisome was reconstructed by combining proteins that are present in both yeast and rat peroxisomal proteomes or that are present in only one of the two proteomes but have orthologs in plants with a peroxisomal location or are described as putative peroxisomal proteins in Araperox database . The minimal ancestral eukaryotic proteome is formed by those proteins of the ancestral opisthokont proteome that are also found in the genomes of plants, Typanosoma brucei and Leishmania major. Catalase and Fox1 that are absent from glycosomes were nevertheless included for the reasons explained in the results and discussion section.
Reviewer's report 1
Arcady Mushegian. Bioinformatics center. Stowers institute for medical research. Kansas City. Missouri. USA.
We now indicate in the table those peroxisomal proteins that show homology with components of the ERAD pathway. In order to provide more detail to figure 2.A we have included the role of ubiquitine in both the processes of ERAD and peroxisomal import. Moreover we have re-arranged the proteins and arrows so that their mechanism of action is clearer. Nevertheless the exact functioning of either ERAD or the PEX5 pathway for protein import has not completely been resolved.
2. What is the identity of actinomycete-like, cyanobacteria-like, and "unresolved" components of peroxisome in Figure 1? What is the explanation for the existence of the first two groups? What can be said about functions of "resolved" vs "unresolved" groups – any trends there?
The identity of the actinomycete-like and cyanobacteria-like proteins is indicated now in the text, they can also be seen in Figure 3 and the table. We do not have a plausible explanation for the origin of these proteins and that is now mentioned in the text. We discuss as well the observed functional dichotomy observed in the proteins with prokaryotic or eukaryotic ancestry for the resolved cases as well as the preponderance of enzymes in the unresolved cases.
3. On the ancestral reconstruction: which parsimony was used – unweighted or weighted? Are opisthocont and eukaryotic sets of PEX genes functionally coherent, or are there missing components?
We used a simple parsimony approach in deciding where certain proteins appeared in evolution: a protein is supposed to have been present at the root of the smallest partition containing all genomes that have that gene. The functional coherence of the PEX subsets is difficult to assess, since many PEX proteins have no specific function assigned. However, at least for the ancestral eukaryotic peroxisome the subset of PEX proteins recovered are all involved in peroxisomal protein import, as indicated in the text.
Reviewer's report 2
Gáspár Jékely. European Molecular Biology Laboratory. Heidelberg, Germany.
This paper makes a compelling argument for the autogenous evolutionary origin of the peroxisome. Although this was not a surprise given recent cell biological findings showing that peroxisomes grow from the endoplasmic reticulum, the autogenous origin of the organelle is now clearly backed by the systematic bioinformatic analysis of its proteome. Most interestingly Gabaldón et al. found that some components of the peroxisomal proteome (the Pex5 pathway) are evolutionary derivatives of the endoplasmic reticulum assisted decay (ERAD) pathway.
The paper is technically sound and well written, I only have a few comments.
1) I have a problem about how the authors define that a protein has eukaryotic origin. For example the Cdc48/Sec18/Pex6 family seems to have descended from archaebacterial AAA ATPases. What the tree shown in Fig. 1A shows is rather that the multiplication of this ancestral ATPase leading to several paralogs was an eukaryotic event. So the protein family clearly has prokaryotic origin, it is the formation of distinct paralogs that occurred during eukaryote evolution. This should be explained better in the text and this group should be referred to differently, like 'originated by eukaryote-specific duplication'.
We specifically want to make a distinction between horizontally transferred genes and "ancient genes" that were already present at the evolutionary split between the lineage leading to the Archaea and the one leading to the eukaryotes. Although in the case of CDC48 and Pex1 a case can indeed be made that CDC48 represents the ancestral function, given its level of sequence identity with the Archaeal sequences, and that PEX1 resulted from a gene duplication, such a clear scenario is rarely present. We have put more emphasis on the distinction between horizontally transferred genes and genes already present in ancient eukaryotes in the text, and mention the CDC48 duplication explicitly now.
2) The reconstruction of the ancestral state of the peroxisomal proteome hinges on the accepted topology of the eukaryotic tree. If Kinetoplastids are not early branching but the root lies between animals and plants, then one would probably get a different picture. This alternative reconstruction should also be presented and/or the effect of tree topology on the results should be discussed.
The consequences of using an alternative topology in the reconstruction of the ancestral proteome are now mentioned in the figure legend. They indeed lead to a larger set of ancestral Peroxisomal proteins.
3) Several of the eukaryote-specific Pex proteins are not discussed in the text. One is left wondering what could have been the evolutionary origin of these proteins. If it is not clear for most of them, this should be mentioned briefly.
We tried hard to ascertain the origin of all Pex proteins, unfortunately for the Pex proteins not discussed in the text we could not find homologies with other proteins of known function or these were too weak to be considered reliable. We now explicitly mention this fact.
Reviewer's report 3
John M. Logsdon, Jr., Department of Biological Sciences, Roy J. Carver Center for Comparative Genomics, University of Iowa, Iowa City, IA 52242 USA
This paper reports the "phylogenomic" analysis of peroxisomal proteins with an aim to distinguish between an endosymbiotic vs. endogenous origin of this organelle in eukaryotic cells. This has been a long-standing question in the evolution of eukaryotic cells and these authors have provided a compelling analysis that rejects the hypothesis that the peroxisome is of endosymbiotic origin. Instead, the data indicate an endogenous origin of peroxisomes from the endoplasmic reticulum.
The authors compiled a curated set of peroxisomal proteins from two major model systems in which global proteomic studies of the peroxisome have been done: yeast and rat. These protein sets, thus, represent a large fraction of the peroxisomal proteome. The authors then used a systematic and rigorous analysis procedure to identify all of the homologs of these proteins from among available complete genomes (prokaryotic and eukaryotic). For all peroxisomal proteins and their homologs, phylogenetic trees were reconstructed and the topologies were evaluated to determine the evolutionary history of each peroxisomal gene. The analysis methodology used is appropriately robust to the questions asked.
1) Although I wholly recommend the publication of this work in Biology Direct, it should be noted that, during the process of review, another paper reporting the same conclusions has appeared as an "Advance Access" publication at Molecular Biology and Evolution:
A. Schlüter, R. Ripp, S. Fourcade, J. L. Mandel, O. Poch, A. Pujol, "The Evolutionary Origin of Peroxisomes: An ER-Peroxisome Connection". I am satisfied that the approaches taken here are sufficiently different than those used by Schlüter et al. and thus merit separate publication. However, it would be helpful for the authors here to refer to the Schlüter et al. paper in their revision and to compare and contrast their approaches and results if at all possible. In addition, it is suggested that the authors consider changing their title so as to not so closely resemble the Schlüter et al. paper.
The Schlüter paper addresses the origin of Peroxisomal proteins without bacterial homologs. As can be seen from our analysis there is actually a conflict in the conclusion one can draw from on the one hand the presence of proteins with alpha-proteobacterial ancestry and on the other hand from the presence of proteins with ER ancestry. One can only resolve this by addressing the retargeting of proteins with alpha-proteobacterial ancestry as we have done in our analysis. We explicitly refer to the Schlüter paper and its observation of the link with the ER in the conclusion, and have changed the title of our paper.
2) Figure 3, and the verbiage associated with it (last paragraph of "Peroxisomal proteins..." section), is confusing and should either be clarified (expanded) or deleted. The figure seems too abstract to be useful. What do the dashed arrows mean?
Figure 3 depicts the moving of DNA and protein localization in evolution, which is rather central to the manuscript. We have rephrased the legend, including an explanation of the meaning of the dashed arrows.
3) Figure 4 is unnecessary to the main message of the paper and could instead be included as a supplement. In fact, it would seem that the phylogeny of this gene family would be a more relevant figure, given the verbiage in the manuscript.
We have left Figure 4 in. It includes the most relevant part of the phylogeny of the citrate synthase genes and does illustrate how the retargeting of proteins has continued in recent evolution and is even visible in the sequences.
4) The phylogenetic trees that are provided as supplementary data are supplied as a single webpage http://www.cmbi.ru.nl/~jagabald/peroxisome_trees.html with the trees given in newick format. Although providing these data in a supplemental format is perfectly acceptable, the authors should provide graphic versions of each tree. Indeed, the abbreviations used for the sequences/taxa in these trees are apparently not defined anywhere in the manuscript. Thus, a key to the taxa is at a minimum required, but even better would be a clear labeling of all of the taxon names on all of trees.
The taxonomic names of the species and the trees are being included.
de Duve D: The peroxisome: a new cytoplasmic organelle. Proc R Soc Lond B Biol Sci 1969,173(30):71-83.
Subramani S: Components involved in peroxisome import, biogenesis, proliferation, turnover, and movement. Physiol Rev 1998,78(1):171-188.
Latruffe N, Vamecq J: Evolutionary aspects of peroxisomes as cell organelles, and of genes encoding peroxisomal proteins. Biol Cell 2000,92(6):389-395. 10.1016/S0248-4900(00)01083-2
Lazarow PB, Fujiki Y: Biogenesis of peroxisomes. Annu Rev Cell Biol 1985, 1: 489-530. 10.1146/annurev.cb.01.110185.002421
Erdmann R, Kunau WH: A genetic approach to the biogenesis of peroxisomes in the yeast Saccharomyces cerevisiae. Cell Biochem Funct 1992,10(3):167-174. 10.1002/cbf.290100306
Hoepfner D, Schildknegt D, Braakman I, Philippsen P, Tabak HF: Contribution of the endoplasmic reticulum to peroxisome formation. Cell 2005,122(1):85-95. 10.1016/j.cell.2005.04.025
Kurland CG, Andersson SG: Origin and evolution of the mitochondrial proteome. Microbiol Mol Biol Rev 2000,64(4):786-820. 10.1128/MMBR.64.4.786-820.2000
Gabaldón T, Huynen MA: Reconstruction of the proto-mitochondrial metabolism. Science 2003,301(5633):609. 10.1126/science.1085463
Schafer H, Nau K, Sickmann A, Erdmann R, Meyer HE: Identification of peroxisomal membrane proteins of Saccharomyces cerevisiae by mass spectrometry. Electrophoresis 2001,22(14):2955-2968. 10.1002/1522-2683(200108)22:14<2955::AID-ELPS2955>3.0.CO;2-U
Yi EC, Marelli M, Lee H, Purvine SO, Aebersold R, Aitchison JD, Goodlett DR: Approaching complete peroxisome characterization by gas-phase fractionation. Electrophoresis 2002,23(18):3205-3216. 10.1002/1522-2683(200209)23:18<3205::AID-ELPS3205>3.0.CO;2-Y
Kikuchi M, Hatano N, Yokota S, Shimozawa N, Imanaka T, Taniguchi H: Proteomic analysis of rat liver peroxisome: presence of peroxisome-specific isozyme of Lon protease. J Biol Chem 2004,279(1):421-428. 10.1074/jbc.M305623200
Jarosch E, Lenk U, Sommer T: Endoplasmic reticulum-associated protein degradation. Int Rev Cytol 2003, 223: 39-81.
Rabinovich E, Kerem A, Frohlich KU, Diamant N, Bar-Nun S: AAA-ATPase p97/Cdc48p, a cytosolic chaperone required for endoplasmic reticulum-associated protein degradation. Mol Cell Biol 2002,22(2):626-634. 10.1128/MCB.22.2.626-634.2002
Harper JW, Burton JL, Solomon MJ: The anaphase-promoting complex: it's not just for mitosis any more. Genes Dev 2002,16(17):2179-2206. 10.1101/gad.1013102
Platta HW, Grunau S, Rosenkranz K, Girzalsky W, Erdmann R: Functional role of the AAA peroxins in dislocation of the cycling PTS1 receptor back to the cytosol. Nat Cell Biol 2005,7(8):817-822. 10.1038/ncb1281
Ye Y, Shibata Y, Yun C, Ron D, Rapoport TA: A membrane protein complex mediates retro-translocation from the ER lumen into the cytosol. Nature 2004,429(6994):841-847. 10.1038/nature02656
Birschmann I, Stroobants AK, van den Berg M, Schafer A, Rosenkranz K, Kunau WH, Tabak HF: Pex15p of Saccharomyces cerevisiae provides a molecular basis for recruitment of the AAA peroxin Pex6p to peroxisomal membranes. Mol Biol Cell 2003,14(6):2226-2236. 10.1091/mbc.E02-11-0752
Erdmann R, Schliebs W: Peroxisomal matrix protein import: the transient pore model. Nat Rev Mol Cell Biol 2005,6(9):738-742.
Marelli M, Smith JJ, Jung S, Yi E, Nesvizhskii AI, Christmas RH, Saleem RA, Tam YY, Fagarasanu A, Goodlett DR, Aebersold R, Rachubinski RA, Aitchison JD: Quantitative mass spectrometry reveals a role for the GTPase Rho1p in actin organization on the peroxisome membrane. J Cell Biol 2004,167(6):1099-1112. 10.1083/jcb.200404119
Hoepfner D, van den Berg M, Philippsen P, Tabak HF, Hettema EH: A role for Vps1p, actin, and the Myo2p motor in peroxisome abundance and inheritance in Saccharomyces cerevisiae. J Cell Biol 2001,155(6):979-990. 10.1083/jcb.200107028
Sickmann A, Reinders J, Wagner Y, Joppich C, Zahedi R, Meyer HE, Schonfisch B, Perschil I, Chacinska A, Guiard B, Rehling P, Pfanner N, Meisinger C: The proteome of Saccharomyces cerevisiae mitochondria. Proc Natl Acad Sci U S A 2003,100(23):13207-13212. 10.1073/pnas.2135385100
Valadi A, Granath K, Gustafsson L, Adler L: Distinct intracellular localization of Gpd1p and Gpd2p, the two yeast isoforms of NAD+-dependent glycerol-3-phosphate dehydrogenase, explains their different contributions to redox-driven glycerol production. J Biol Chem 2004,279(38):39677-39685. 10.1074/jbc.M403310200
Janssen U, Fink T, Lichter P, Stoffel W: Human mitochondrial 3,2-trans-enoyl-CoA isomerase (DCI): gene structure and localization to chromosome 16p13.3. Genomics 1994,23(1):223-228. 10.1006/geno.1994.1480
Abdelraheim SR, Spiller DG, McLennan AG: Mammalian NADH diphosphatases of the Nudix family: cloning and characterization of the human peroxisomal NUDT12 protein. Biochem J 2003,374(Pt 2):329-335. 10.1042/BJ20030441
Scharfe C, Zaccaria P, Hoertnagel K, Jaksch M, Klopstock T, Dembowski M, Lill R, Prokisch H, Gerbitz KD, Neupert W, Mewes HW, Meitinger T: MITOP, the mitochondrial proteome database: 2000 update. Nucleic Acids Res 2000,28(1):155-158. 10.1093/nar/28.1.155
Solaas K, Ulvestad A, Soreide O, Kase BF: Subcellular organization of bile acid amidation in human liver: a key issue in regulating the biosynthesis of bile salts. J Lipid Res 2000,41(7):1154-1162.
Elgersma Y, van Roermund CW, Wanders RJ, Tabak HF: Peroxisomal and mitochondrial carnitine acetyltransferases of Saccharomyces cerevisiae are encoded by a single gene. Embo J 1995,14(14):3472-3479.
Koch A, Yoon Y, Bonekamp NA, McNiven MA, Schrader M: A Role for Fis1 in Both Mitochondrial and Peroxisomal Fission in Mammalian Cells. Mol Biol Cell 2005.
Birdsey GM, Lewin J, Cunningham AA, Bruford MW, Danpure CJ: Differential enzyme targeting as an evolutionary adaptation to herbivory in carnivora. Mol Biol Evol 2004,21(4):632-646. 10.1093/molbev/msh054
Danpure CJ, Lumb MJ, Birdsey GM, Zhang X: Alanine:glyoxylate aminotransferase peroxisome-to-mitochondrion mistargeting in human hereditary kidney stone disease. Biochim Biophys Acta 2003,1647(1-2):70-75.
Reumann S, Ma C, Lemke S, Babujee L: AraPerox. A database of putative Arabidopsis proteins from plant peroxisomes. Plant Physiol 2004,136(1):2587-2608. 10.1104/pp.104.043695
Soares MJ, De Souza W: Cytoplasmic organelles of trypanosomatids: a cytochemical and stereological study. J Submicrosc Cytol Pathol 1988,20(2):349-361.
Thieringer R, Kunau WH: The beta-oxidation system in catalase-free microbodies of the filamentous fungus Neurospora crassa. Purification of a multifunctional protein possessing 2-enoyl-CoA hydratase, L-3-hydroxyacyl-CoA dehydrogenase, and 3-hydroxyacyl-CoA epimerase activities. J Biol Chem 1991,266(20):13110-13117.
Hannaert V, Saavedra E, Duffieux F, Szikora JP, Rigden DJ, Michels PA, Opperdoes FR: Plant-like traits associated with metabolism of Trypanosoma parasites. Proc Natl Acad Sci U S A 2003,100(3):1067-1071. 10.1073/pnas.0335769100
Schluter A, Ripp R, Fourcade S, Mandel JL, Poch O, Pujol A: The Evolutionary Origin of Peroxisomes: An ER-Peroxisome Connection. Mol Biol Evol 2006.
Saccharomyces Genome Database
IPI: IPI database.[http://www.ebi.ac.uk/IPI]
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
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994,22(22):4673-4680.
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
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 2000,28(1):27-30. 10.1093/nar/28.1.27
material S: Supplementary material.Supplemantary material. [http://www.cmbi.ru.nl/~jagabald/peroxisome_trees.html]
Letunic I, Copley RR, Schmidt S, Ciccarelli FD, Doerks T, Schultz J, Ponting CP, Bork P: SMART 4.0: towards genomic data integration. Nucleic Acids Res 2004, 32 Database issue: D142-4. 10.1093/nar/gkh088
Soding J: Protein homology detection by HMM-HMM comparison. Bioinformatics 2005,21(7):951-960. 10.1093/bioinformatics/bti125
Altschul SF, Koonin EV: Iterated profile searches with PSI-BLAST--a tool for discovery in protein databases. Trends Biochem Sci 1998,23(11):444-447. 10.1016/S0968-0004(98)01298-5
Shiozawa K, Maita N, Tomii K, Seto A, Goda N, Akiyama Y, Shimizu T, Shirakawa M, Hiroaki H: Structure of the N-terminal domain of PEX1 AAA-ATPase. Characterization of a putative adaptor-binding domain. J Biol Chem 2004,279(48):50060-50068. 10.1074/jbc.M407837200
Baldauf SL, Roger AJ, Wenk-Siefert I, Doolittle WF: A kingdom-level phylogeny of eukaryotes based on combined protein data. Science 2000,290(5493):972-977. 10.1126/science.290.5493.972
We thank the referees for their many helpful suggestions ideas and corrections. We also thank Eugene Koonin for his comments on the manuscript. This work was supported by a grant from the Netherlands organization for scientific research (NWO). TG is recipient of an EMBO post-doctoral fellowship, ALTF 402-2005.
TG and MAH conceived the study, participated in its design and coordination. TG, MAH, BS and FvZ participated in the computational analysis of peroxisomal proteins. TG, MAH, BS and HT drafted the manuscript. WH provided important insights in the study and the results. All authors read and approved the final manuscript.