- Research article
- Open Access
Conservation of the TRAPPII-specific subunits of a Ypt/Rab exchanger complex
BMC Evolutionary Biology volume 7, Article number: 12 (2007)
Ypt/Rab GTPases and their GEF activators regulate intra-cellular trafficking in all eukaryotic cells. In S. cerivisiae, the modular TRAPP complex acts as a GEF for the Golgi gatekeepers: Ypt1 and the functional pair Ypt31/32. While TRAPPI, which acts in early Golgi, is conserved from fungi to animals, not much is known about TRAPPII, which acts in late Golgi and consists of TRAPPI plus three additional subunits.
Here, we show a phylogenetic analysis of the three TRAPPII-specific subunits. One copy of each of the two essential subunits, Trs120 and Trs130, is present in almost every fully sequenced eukaryotic genome. Moreover, the primary, as well as the predicted secondary, structure of the Trs120- and Trs130-related sequences are conserved from fungi to animals. The mammalian orthologs of Trs120 and Trs130, NIBP and TMEM1, respectively, are candidates for human disorders. Currently, NIBP is implicated in signaling, and TMEM1 is suggested to have trans-membrane domains (TMDs) and to function as a membrane channel. However, we show here that the yeast Trs130 does not function as a trans-membrane protein, and the human TMEM1 does not contain putative TMDs. The non-essential subunit, Trs65, is conserved only among many fungi and some unicellular eukaryotes. Multiple alignment analysis of each TRAPPII-specific subunit revealed conserved domains that include highly conserved amino acids.
We suggest that the function of both NIBP and TMEM1 in the regulation of intra-cellular trafficking is conserved from yeast to man. The conserved domains and amino acids discovered here can be used for functional analysis that should help to resolve the differences in the assigned functions of these proteins in fungi and animals.
In all eukaryotic cells, intra-cellular trafficking connects the cell with its environment by the orderly transport of membranes and proteins via the exocytic and endocytic pathways. In the exocytic pathway, proteins destined to be secreted or presented on the plasma membrane (PM) are transported from the endoplasmic reticulum (ER), through the Golgi apparatus, to the PM. In the endocytic pathway, proteins from the environment or the PM are shuttled via a set of endosomes to lysosomes. The machinery and the mechanisms of intra-cellular trafficking are highly conserved among all eukaryotes (orthologs), and some protein components are also conserved between the various steps of the pathways (paralogs) .
Ypt/Rab GTPases are key regulators of this protein trafficking. They are conserved both within a single genome between transport steps and across large phylogenetic distances [2–4]. Saccharomycescerivisiae cells contain 11 Ypts, whereas human cells have ~70 Rabs [5–7]. In S. cerivisiae, Ypt1 and the functional pair Ypt31/32 regulate entry into and exit from the Golgi, respectively [8, 9]. The mammalian Rab1A and Rab1B share ~70% identity with S. cerivisiae Ypt1, and Rab1A can replace Ypt1 in yeast knockout cells . The mammalian Rab11A, Rab11B and Rab25 share ~60% identity with S. cerivisiae Ypt31/32 and regulate the same transport steps: exit from the Golgi and endosome-to-Golgi transport [11, 12].
Ypt/Rabs are activated by specific nucleotide exchangers, called guanine nucleotide exchange factors (GEFs). GEFs for different Ypt/Rab paralogs do not share sequence similarity and are therefore harder to identify. In S. cerivisiae, the multi-subunit complex TRAPP was identified as the GEF for both Ypt1 [13, 14], and Ypt31/32 . TRAPP is a modular complex that exists in two forms: TRAPPI and TRAPPII . We recently showed that TRAPPI acts as a Ypt1 GEF, whereas TRAPPII functions as a Ypt31/32 GEF .
The S. cerivisiae TRAPPI was shown to function in ER-to-Golgi transport. It contains seven subunits that co-precipitate as a ~300 kDa complex from yeast cell lysates . There is considerable evidence for the structural and functional conservation of the TRAPPI complex (excluding the Trs85 subunit) from fungi to animals  (Table 1).
TRAPPII is a ~1000 kDa complex that contains three subunits in addition to those of TRAPPI: two large, essential subunits (>1000 amino acids), Trs120 and Trs130, and one small, non-essential subunit, Trs65 . In yeast, Trs130 was implicated in late Golgi transport, whereas Trs120 was suggested to function in endosome-to-Golgi transport [15, 18]. The conservation of the TRAPPII complex is less clear than that of TRAPPI . Blast analyses show that the closest mammalian homologues of Trs120 and Trs130, NIBP and TMEM1, respectively, share a ~20% identity with the S. cerivisiae proteins over about a third of the protein (Table 1) as compared to 23–34% identity for five of the seven TRAPPI-specific proteins and 56% for Bet3 (Trs85 is not conserved, ). Only TMEM1 was shown to co-precipitate with the human TRAPP complex , and there is no functional evidence for a role for either of these proteins in protein trafficking. Instead, NIBP and TMEM1 were implicated in very different functions, NF-kappaB signaling and as a membrane channel, respectively [21–24].
There is no obvious mammalian ortholog for Trs65. However, current evidence supports a role for Trs65 in the S. cerivisiae TRAPPII complex. First, deletion of TRS65 is synthetic lethal with deletion of the non-essential TRAPPI subunit, TRS33. This synthetic lethality can be rescued by over-expression of Ypt31 [25, 26]. In addition, Trs65/Kre11 was shown to function in cell-wall biogenesis , a process that is dependent on protein transport.
Because of the low sequence conservation by Blast analysis of the TRAPPII-specific subunits (Table 1), and the discrepancy between their suggested functions in fungal and mammalian cells, it is important to establish any evolutionary conservation more precisely. Here, we show a phylogenetic analysis of the TRAPPII-specific subunits. We found that nearly every fully sequenced eukaryotic proteome present in the NCBI non-redundant protein database (as for July 06) contains one Trs130 and one Trs120 homolog. We found three exceptions, which are noted in Results. In contrast, Trs65 is conserved among many, but not all, unicellular eukaryotes, and not in any fully sequenced multi-cellular proteome. The conservation of Trs120 and Trs130 is supported by conserved predicted secondary structure of whole proteins, and conserved domains that include highly conserved amino acids, which should also help in future functional analysis. In addition, our analysis suggests that neither Trs130 nor TMEM1 is a trans-membrane protein, arguing against the suggestion that TMEM1 functions as a membrane channel. Understanding the function of the human homologues is important, because they are candidates for several human disorders (Table 1).
Phylogenetic trees of TRAPPII-specific subunits
We searched the non-redundant protein NCBI database on 07/2006 for sequences related to the three S. cerivisiae TRAPPII-specific subunits: Trs120, Trs130, and Trs65 (Table S1-A in Additional File 1). We identified 35 Trs120-related protein sequences and 31 TRS130-related protein sequences. We identified no more than one ortholog for each of these proteins in all examined genomes from fungi to vertebrates. In the 14 fully and 23 largely sequenced genomes, we always identified exactly one copy. We found sixteen Trs65-related sequences, also no more than one per genome, but only in fungi. Further, even some fully sequenced fungi (e.g., Schizosaccharomyces pombe) lack any recognizable ortholog of this protein (Figure 1 and Table S1-A in Additional File 1).
In addition to the sequences found in the NCBI database at the time of our search, we also looked at a 24 other largely completed genomes (including 19 additional fungi and 5 protists) that were at varying stages of assembly and annotation as of December 06. In spite of the incomplete nature of these projects, we still found Trs120 and Trs130 in 21 genomes, and Trs65 in 17 genomes. We show these results in Table S1-B (in Additional File 1).
There are a few interesting exceptions to our observations of the universality of Trs120- and Trs130-related sequences in eukaryotes. The three eukaryotic genomes that did not contain Trs120 and Trs130 are: the fungus Encephalitozoon ecuniculi, the Amoebazoan Entamoeba histolytica and the Archaeplastid Cyanidioschyzon merolae. It seems reasonable to suggest that the absence of Trs120 and Trs130 from the proteome of the first organism can be explained by the fact that it is an obligate intracellular parasite that might lack some essential cellular machinery. In agreement, three essential TRAPPI subunits are also missing in its genomes (Table S1-C in Additional File 1). E, histolytica CMerolae, however, are well-sequenced genomes, the latter is a red alga with a very small genome. Interestingly, the two genomes are also missing one essential TRAPPI subunit each (Table S1-C in Additional File 1). If the absence of TRAPP subunits in these organisms is confirmed, it might shed light on the evolution of the intra-cellular trafficking machinery.
The full-length protein sequences from the NCBI non-redundant (nr) database were used to construct un-rooted phylogenetic trees for each TRAPPII-specific subunit employing programs described in Methods. We used two different methods to construct the trees: PAUP's distance method, and PHYML's maximum likelihood (ML) method. Both methods yielded essentially identical topologies, except for branches near the center of the tree representing far-flung subunit-family members. Only branch lengths differed significantly between the two methods, with the ML method giving consistently ~4 fold longer branch lengths. Because we could estimate branch retention probabilities using PAUP, we show just the PAUP method results in Figure 1. The reliability of each branch position was assessed by distance-based bootstrap analysis and indicated by red (weak) to green (strong) colors. In almost all cases branches were maintained in almost all bootstrap trials, excepting only those close to the center of the tree, which reflect the details of ancient divergence between phylogenetic groups. The colored pie sections delineate groups of proteins in accord with tree branching and taxonomic relationships. There are seven basic groups: animals (A), plants (P), slime molds (M), and four fungal (F) groups Saccharomycetes (FS), Pezizimycotina (FP), Basidiomycota (FB), and Schizosaccharomycetes (FZ). Grouping in the tree corresponds in almost all cases to recognized taxonomic groupings . The sole exception was Yarrowia lipolytica (FS_Yl), which is a Saccharomycete, but falls (weakly) with the Pezizomycotinas on the Trs130 tree.
Multiple alignment analysis of TRAPPII-specific subunits
Multiple alignments and the domain structure of the TRAPPII-specific subunits are shown in Figures 2, 3, 4 and Figures S1-S3 (see Additional File 1). Domains were defined by inspection of the raw alignment and Figures S1-S3. Highly conserved (HC) amino acids in the alignment (see Methods) are shown at the bottom of the diagram as short vertical lines, and are detailed in the corresponding supplementary tables (Tables S2-S4, respectively, see Additional File 1). There is a number of particularly interesting HC amino acids, including those that are invariant and rare (e.g., prolines) or invariant and potentially catalytic (e.g., aspartic or glutamic acids). These HC amino acids are indicated with asterisks in the multiple alignment and in the supplemental tables (Tables S1-S3, see Additional File 1). Finally, conserved domains, as indicated by alignment strength and HC amino acids are boxed (Figures S1-S3) and the alignment of these conserved domains from representative organisms is shown in Figures 2, 3, 4.
Multiple alignment analysis of Trs120-related proteins is shown in Figures 2 and S1. This figure shows five well-conserved domains across the length of Trs120, each characterized by high-quality sequence similarities and by the presence of HC amino acids in the alignment. These domains vary in size from 20 to 165 amino acids. Larger-scale sequence conservation and HC amino acids agree with one another: HC amino acids (detailed in Table S2) are found only in the otherwise conserved regions (except for a lone conserved glycine in a region of weak sequence conservation, C-terminal to domain TRS120-V.)
Trs120 is essential for S. cerivisiae cell viability. However, truncation of the C-terminal third of Trs120 in S. cerivisiae (amino acids 808 through 1289, see arrows in Figures 2B and S1), which includes part of conserved domain IV and all of domain V, results in a temperature-sensitive growth phenotype. Deletions further N-terminal of this point are lethal . Therefore, we predict that domains I, II, and III and part of IV are important for the essential interactions of Trs120.
Multiple alignment analysis of Trs130-related proteins is shown in Figures 3 and S2. In this figure, the worms Caenorhabditiselegans and Caenorhabditisbriggsiae are slightly different in that they share a C-terminal deletion. Alignment and HC amino acids show three Trs130 domains ranging from 33–314 amino acids. We found no significant HC amino acids (detailed in Table S3) outside of these domains. Domains I and II contain many highly unusual HC amino acids. Consequently, we predict that these domains are critical for essential interactions or functions of Trs130.
Domain III is unlikely to be strongly important. It is shorter and contains no significant HC amino acids. Truncation C-terminal to the arrow shown in this domain in S. cerivisiae results in a temperature-sensitive growth phenotype [15, 26]. While it is true that this domain might be important for the stability of the Trs130 protein, based on the low level of the protein in mutant cells , it might not be true for all organisms, since the two worm species share a C-terminal truncation of Trs130 (confirmed as described in Methods) that includes part of Trs130-II and all of Trs130-III.
Multiple alignments of Trs65-related proteins are shown in Figures 4 and S3. The protein is common for Saccharomyceta and Pezizomycotina subphyla (and at least one Basidiomycota), but not all fungi. Table S1 (see Additional File 1) shows that Trs65 also occurs in at least two other non-fungal unicellular eukaryotes. Since Trs65 is only found in some fungi, we used all three fungal groups shown in Figure 1 that contain Trs65 sequences. HC amino acids (detailed in Table S4) and alignment qualities show that these organisms share four conserved domains varying in size 40–140 amino acids. Although no mutagenesis data exists for these proteins, there are some highly unusual amino acids in domains II, III, and IV. Therefore we suggest that they are important for Trs65 function.
The predicted secondary structure of Trs120 and Trs130 is conserved
We suspected that protein secondary structure was conserved in addition to its primary sequence. Since no crystal data is currently known for the three TRAPPII-specific subunits, we predicted their secondary structures using Prof . This program employs multiple sequence alignments to predict the helical, beta-sheet, or coiled-coil nature of each position in the primary sequence of a group of proteins. We examined these predictions for all the Trs120, Trs130, and Trs65 proteins, using the clustal alignment employed to construct the trees and alignments in Figures 1, 2, 3, 4 for the input to Prof. Color-coded bars, representing the predicted secondary structures (helix-red; coil-blue; and beta sheet-green) are shown in Figures 5, S4 and S5, along with the conserved domains and HC amino acids described in Figures 2, 3, 4. To be sure that differences among evolutionary groups were in agreement with these results, we also chose three groups of species that contain a large number of sequenced proteins (6–13) and that are evolutionary far from each other: animals (A), and two fungi groups, Saccharomycetes (FS), and Pezizomycotina (FP), and show these results in Figures S4 and S5.
Surprisingly, the predicted secondary structure of Trs120 and Trs130 is highly conserved not only in domains that are conserved for primary structure, but also along the whole sequence of these proteins, even where there is low sequence identity between the groups (Figures 5 and S4). Therefore, the predicted secondary structure of these two proteins is conserved better than their primary structure. Similar analysis of Trs65 among the Saccharomycetes and Pezizomycotina fungal groups shows conservation of predicted secondary structure among members of the same group. However, these predictions do not agree well between the two subphyla (Figure S5). This suggests that unlike Trs120 and Trs130, the interactions, and perhaps the specific functions, of Trs65 are not well conserved even among all fungi.
For the three protein families, we searched all studied sequences (from all organisms) for potential structural motifs using SMART [30, 31]. We found some isolated SMART hits, but none were consistent across a whole group (i.e., all animals, all of a fungal subphylum, etc). Notably, none of the isolated hits overlapped the conserved domains shown in Figs 2, 3, 4. Portions of the conserved domains of TRAPPII-specific subunits seem good prospects for inclusion in future releases of such motif-finding engines.
Trs130 is a membrane-associated protein
Based on sequence analysis, the S. cerevisiae Trs130 and its human homolog were suggested in 1997 to contain trans-membrane domains (TMDs). In fact, the human homologue was named TMEM1, for Trans-Membrane Epilepsy Myoclonus . However, the yeast TRAPP complex is only membrane associated, because it can be removed from membranes by salt but not by detergent [32, 33], even though this was not shown specifically for Trs130.
Because the S. cereviciae Trs130 contains two putative trans-membrane domains , we wished to determine whether this protein is a trans-membrane or a membrane-associated protein. Trs130 isolated from yeast cell lysates is found in the P100 (pellet of 100,000 × g) particulate fraction of cell lysates (Figure 6, panel A), which floats with membranes on an OptiPrep gradient (Figure 6, panel B). If Trs130 were a trans-membrane protein, it could be extracted from the membrane by detergent, but not by salt. However, treatment of the P100 fraction with salt, and not with detergent, yields soluble Trs130 (Figure 6, panel C). Therefore, the yeast Trs130 behaves as a membrane-associated protein, and not as a trans-membrane protein. These results imply that the observed hydrophobic domains do not function as TMDs, but rather are important for the proper folding or interactions of Trs130.
This idea that TMEM1 and Trs130 are trans-membrane proteins was based on sequence analysis, so we revisited this sequence analysis for the human TMEM1. Previously, sequence analysis of the human TMEM1 suggested that the protein contained two (using SOSUI, ) or multiple  trans-membrane domains (TMDs). We tested the human TMEM1 sequence for putative TMDs using SOSUI  and SMART  on 6/06, and no TMDs were detected by these programs. Note that SMART uses TMHMM2.0 to predict TMDs, a method which has been shown to be among the most accurate at such tasks . We suspect that the TMD-detecting algorithms have become more discriminating in the intervening years and suggest that, like the S. cerivisiae Trs130, the human protein is a membrane-associated protein, not a trans-membrane protein as suggested by its original name.
We show here that the TRAPPII-specific subunits, Trs120 and Trs130, are conserved from fungi to animals. This conclusion is based on sequence and predicted secondary structure analyses. Multiple alignments analysis reveals conserved domains with highly conserved amino acids in each protein family. Predicted secondary structure analysis reveals conservation in the clustal/HC-derived domains as well as in the intervening sequences. The conservation of the predicted secondary structure and the fact that the secondary structure of these proteins is conserved better than their primary structure supports the analysis presented here. This suggests that TRAPPII, like TRAPPI, is a conserved complex, and that the function of both TRAPP complexes, in the regulation of trafficking through the Golgi apparatus, is also conserved.
In contrast, the non-essential subunit, Trs65, is conserved only between some fungi. The Trs65-related proteins also share conserved domains with highly conserved amino acids. However, their secondary structures as predicted by Prof are not highly conserved between two evolutionary distant fungal groups, Saccharomycetes (FS) and Pezizomycotina (FP). The fact that Trs65 is not essential for viability in S. cerivisiae, and is only weakly conserved among some unicellular eukaryotes, suggests that Trs65 has a fungi-specific function. Indeed, it was implicated in S. cerivisiae in cell-wall biogenesis . Alternatively, other organisms may have a functionally similar subunit that does not share sequence similarity with Trs65.
The S. cerivisiae TRAPP complexes act as activators, GEFs, for the Golgi Ypts [13, 16, 32]. In general, Ypt/Rab GEFs are large protein complexes , and there is no similarity between GEFs for different Ypt/Rab paralogs. This suggests the recruitment of Ypt/Rab GEFs from divergent families. The reason for this might stem from the divergence of factors that regulate these GEFs, about which, very little is currently known. Thus, the diversity of the Ypt/Rab GEFs suggests that GEFs that function in different cellular compartments are regulated by diverse upstream factors. Alternatively, Ypt/Rab GEFs might have diverged beyond recognition, which is hard to believe in light of the conservation of the whole trafficking machinery.
Unlike the Ypt/Rabs, for which the number between fungi and mammalian cells increases by about seven fold, the number of TRAPPI and TRAPPII subunits per genome remains constant through evolution (except for two mammalian paralogs of Trs33, A and B, see Table 1). This ratio of Rab/GEF suggests that in mammalian cells TRAPP acts as a GEF for more Rabs than in S. cerivisiae, perhaps for whole groups of Ypt1/Rab1- and Ypt31/Rab11-related proteins . Alternatively, more animal GEFs evolved during evolution that do not share sequence similarity with the fungal GEFs.
Domains and functions
The conservation of the primary and the secondary structures of Trs120 and Trs130 suggests that they cooperatively execute an important cellular function. We propose that conservation of secondary structure of the well-conserved domains (defined by sequence conservation and the presence of highly conserved amino acids) probably reflects conservation of catalytic functions and/or protein-protein interactions of these domains. The latter are known to place constraints on the divergence of residues in contact surfaces. In contrast, conservation of the secondary structure of the intervening sequences might be required for the three-dimensional orientation of the well-conserved domains.
We have previously suggested a role for the TRAPP complex in the coordination of entry into and exit from the Golgi [14, 16], and for the TRAPPII-specific subunits, in the specificity switch of the GEF activity of TRAPP from Ypt1 to Ypt 31/32 (, and Figure 7). Based on our functional analysis of the S. cerivisiae Trs120 and Trs130, we propose a number of protein-protein interactions for the two essential TRAPPII-specific subunits. We predict at least two protein-protein interactions for Trs120, with Trs130 and with TRAPPI . The four Trs120 conserved domains (I-III, and part of IV), which are essential for S. cerevisiae, are candidates for these interactions. For Trs130, we predict at least two protein-protein interactions and a possible catalytic role in the Ypt31 GEF activity of TRAPPII . The two conserved Trs130 domains that are essential for S. cerevisiae viability, I and II, are candidates for these functions. The identification of Trs120- and Trs130-conserved domains and HC amino acids will allow future functional analysis of the separate domains of these large proteins.
Proposed roles for the mammalian Trs120- and Trs130-related proteins
We suggest that the role of the two TRAPP complexes as Ypt/Rab GEFs and their coordination function are conserved from fungi to animals. However, the mammalian orthologs of Trs120 and Trs130 were implicated in other cellular processes. The mammalian Trs120-related protein, NIBP was implicated in NF-kappaB signaling . While the suggested functions of the fungal and mammalian orthologs of Trs120 do not overlap, they do not contradict each other, and it is possible that trafficking through the Golgi in mammalian cells is regulated by NF-kappaB signaling.
In contrast, the mammalian Trs130-related protein, TMEM1 (for Trans Membrane Epilepsy Myoclonus), was implicated in a function distinct from protein trafficking. While the S. cerevisiae TRAPP, including Trs130, is a membrane associated complex, TMEM1 was suggested to be a trans-membrane protein and to function as a membrane channel [23, 24]. Our analyses suggest that both the S. cereviciae Trs130 (by membrane-extraction analysis) and the human TMEM1 (by sequence analysis) are not trans-membrane proteins. Based on our alignment and predicted secondary structure analyses, we propose that the function of Trs130 as a specificity switch for TRAPP's GEF activity is conserved. These ideas can now be tested.
Trs120 and Trs130 and human disorders
Regulation of critical steps of the exocytic pathway is bound to be crucial in understanding the basis of every disease that is connected with secretion of a substance or presentation of a receptor on the plasma membrane. Indeed, several TRAPPI subunits were implicated in human diseases. For example, a missense mutation in human Bet5/TRAPPC1 results in the expression of antigenic peptides in melanoma ; and the SEDL/Trs20 gene is responsible for SEDT, an Xlinked skeletal disorder . The two mammalian orthologs of Trs120 and Trs130 were also implicated in human disorders: NIBP, the Trs120 ortholog, was implicated in neurogenerative disorders, based on its connection to NF-kappaB signaling . The Trs130 ortholog, TMEM1, is a candidate for several human disorders, including certain types of epilepsy, autoimmune, and holoprosencephaly disorders, based on genetic linkage studies mapping genes responsible for these disorders to a chromosome region that includes TMEM1 [23, 24, 40]. Therefore, resolving the probably conserved function of the mammalian NIBP and TMEM1 proteins should help elucidate the basis for these disorders.
TRAPP is a modular protein complex that regulates entry into and exit from the Golgi, which is a cellular compartment central to multiple trafficking pathways. TRAPPI, the complex that functions in the entry to the Golgi, is highly conserved. We have recently shown that the two yeast essential TRAPPII-specific subunits are required for changing the activity of TRAPPI to that of TRAPPII . Here, we explored the conservation of the TRAPPII-specific subunits. We show that the primary and secondary structures of the two essential subunits, Trs120 and Trs130, are conserved from yeast to man, whereas the non-essential Trs65 is conserved only among some unicellular eukaryotes. We suggest that the role of the human orthologs of Trs120 and Trs130, NIBP and TMEM1, respectively, in the regulation of intra-cellular trafficking is also conserved. Conserved domains and highly conserved amino acids in these domains should help future functional studies of these proteins in higher eukaryotes.
Identification of TRAPPII-specific subunit sequences
We searched the non-redundant protein database from NCBI on 7/06 . Trs120-containing proteins were identified by iteratively performing BLAST analysis of this database with the S. cerivisiae Trs120 sequence with an expectation value of 10-5. Hits from the first and subsequent rounds of BLAST analysis were used to perform additional analyses until no more sequences were discovered. An identical procedure was used for the identification of Trs130- and Trs65-containing proteins, again using the S. cerivisiae sequences as seeds.
Table S1-A shows all identified sequences including abbreviated names used in all figures, protein sequence length, accession numbers, and annotations from the NCBI database. Full-length versions of each sequence were identified for all further work, discarding truncations or deletions. We used these Trs120, Trs130, and Trs65 proteins to search 24 additional, largely completed genomes in various stages of completion and annotation. We searched in both called proteins sequences and in the raw genome using tblastn. The sequence source and the presence or absence of Trs120, Trs130, and Trs65 in these genomes is shown in Table S1-B.
For Trs130, the two worm sequences, C. briggsae and C. elegans, seemed to be truncated relative to the remaining Trs130-containing proteins. We searched both completed genomes in all six reading frames using the universal translation table for sequences corresponding to the C-terminal remainder of any of the other animal proteins. We found no significant blast hits.
Construction of the phylogenetic trees
The identified Trs120, Trs130, and Trs65 proteins were aligned by clustal , using the default settings for slow/accurate alignments (gap penalty of 10, gap extension cost of 0.2, 30% delay for divergent sequences, 4 space gap separation distance, without end gap separation, with residue-specific penalties, and using the Gonnet series protein weight matrix). The aligned sequences were manually trimmed on the N- and C-terminal ends to remove weak or ambiguous alignments. Phylogenetic analysis was performed using PAUP 4.0b10  and PHYML . Distance trees were created using the heuristic distance search for optimal trees by PAUP and the maximum likelihood method by PHYML. The PAUP tree was started with neighbor joining, and branch swapping used the TBR algorithm. We used default settings for the ML method (JTT substitution model, 1 substitution rate category, optimization on for branch lengths and topology). To estimate the reliability of the PAUP tree, bootstrap analysis was performed with 1000 replicates of full heuristic searches, using the same weighting parameters employed in the initial analysis. The trees were drawn by TreeView version 1.6.6  and manually modified in a general-purpose graphics editor. Branches were color coded from red to green to reflect their persistence in bootstrap analysis in 50–100% of sampled trees. We constructed unrooted trees for the TRAPPII-specific subunits because we could not find any sequences that are related to them, even at an extremely lenient blast cutoff (e<10-2), and are unambiguously diverged from all eukaryotes, e.g., from archaeal or bacterial proteomes.
Domain structure and identification of Highly Conserved (HC) amino acids
Domains were defined by inspection of the raw clustal files and of the graphical multiple alignments shown in Figures S1-S3. Domains were consistently present in all sequences in all groups of organisms examined.
We examined the clustal-generated alignments for highly-conserved (HC) amino acids, using popular amino acid groupings . Alignment positions were considered HC if more than 90% of the residues belonged to an amino acid group. Since some groups (MILV, MILF, and SAT) are quite common, we discarded these positions if they were more than 20 amino acids away from any other HC position.
Prediction of secondary structures
We predicted sequence secondary structures using Prof , a secondary structure predictor that uses multiple alignments to inform its prediction about the helical, beta-sheet, or coiled nature of each position in the primary sequence. We used the same clustal-generated alignment used for phylogenetic trees and annotated alignments for Prof analysis. To be sure that its results were consistent, we gave Prof either all the aligned sequences (ALL) just those from animals (A), just those from Saccharomycetes (FS), or just those from Pezizomycotina (FP). We wrote a custom perl program to draw these predicted secondary structures as red (helix), green (beta-sheet) or blue (coiled) bars alongside HC and domain annotations.
The following yeast strain was used in this study: wild type, NSY991 (VSY459; MATa leu2-3,112 his3-200 trp1-901 lys2-801 suc2-9 ura3-52 TRS130-HA:HIS3MX6) . Antibodies used in this study: Mouse monoclonal Anti-HA (clone12CA5, Roche); affinity purified rabbit anti-Ypt31 ; rabbit anti-EMP47 (gift from H. Riezman, ); and horseraddish Peroxidase linked Anti-rabbit and anti-mouse IgG (Amersham Biosciences). All chemical reagents were purchased from Sigma (St. Louis, MO), unless otherwise noted. Iodixanol density gradients, Optiprep™, were purchased from Axis-Shield PoC AS (Oslo, Norway).
Preparation of cell lysates and protein analyses
Yeast cells were grown in rich (YPD) medium . Yeast cell extracts were prepared as previously described . Cell breakage buffers were supplemented with an EDTA-free protease inhibitor cocktail (Roche Diagnostics, Indianapolis, IN). Protein concentrations were determined by BioRad protein assay (BioRad). 10 μg of yeast whole-cell lysates were loaded on 7.5–10% SDS-PAGE. Gels were run, and proteins were transferred to PVDF membranes and subjected to immuno-blot analysis. Quantification of protein bands was done using the AlphaEase FC and Alpha-Imager (Alpha Innotech Corporation).
Membrane attachment analysis
Cell fractionation: Yeast cell lysates were prepared as described above with the following alterations: Frozen cell pellets (22–25 OD600 units) were resuspended and broken in 100 μl spheroplast buffer  by vortexing 3 times for 2 minutes each at 4°C. 100 μl of spheroplast buffer was added and the supernatant was separated from cell debris as previously described . 100 μl of the supernatant was saved to be used as total cell lysates and the other 100 μl was taken for further centrifugation at 100,000 × g for 30 minutes at 4°C. The supernatant (S100) was separated from the pellet (P100). P100 was resuspended in spheroplast buffer to the same volume as S100. 75 μg total cell lysates and an equivalent of 1.5 OD600 units for P100 and S100 were subjected to immuno-blot analysis.
Iodixanol Density Gradient: Resuspended P100 fractions were analyzed on 30% iodixanol density gradients (Optiprep™, Axis-Shield PoC AS) according to the protocol provided by the company. Fractions were analyzed by immuno-blot analysis.
Membrane Extraction: P100 was resuspended in 100 μl B88 alone, or supplemented with 1% Triton X-100, or with 0.5 M NaCl as previously described . Samples were centrifuged at 100,000 × g and subjected to immuno-blot analysis.
Guanine nucleotide Exchange Factor
Trans Membrane Domain.
Bock JB, Matern HT, Peden AA, Scheller RH: A genomic perspective on membrane compartment organization. Nature. 2001, 409 (6822): 839-841. 10.1038/35057024.
Pfeffer SR: Structural clues to Rab GTPase functional diversity. J Biol Chem. 2005, 280 (16): 15485-15488. 10.1074/jbc.R500003200.
Segev N: Ypt/rab gtpases: regulators of protein trafficking. Sci STKE. 2001, 2001 (100): RE11-
Zerial M, McBride H: Rab proteins as membrane organizers. Nat Rev Mol Cell Biol. 2001, 2 (2): 107-117. 10.1038/35052055.
Buvelot Frei S, Rahl PB, Nussbaum M, Briggs BJ, Calero M, Janeczko S, Regan AD, Chen CZ, Barral Y, Whittaker GR, Collins RN: Bioinformatic and comparative localization of rab proteins reveals functional insights into the uncharacterized GTPases ypt10p and ypt11p. Mol Cell Biol. 2006, 26 (19): 7299-7317. 10.1128/MCB.02405-05.
Colicelli J: Human RAS superfamily proteins and related GTPases. Sci STKE. 2004, 2004 (250): RE13-10.1126/stke.2502004re13.
Pereira-Leal JB, Seabra MC: The mammalian Rab family of small GTPases: definition of family and subfamily sequence motifs suggests a mechanism for functional specificity in the Ras superfamily. J Mol Biol. 2000, 301 (4): 1077-1087. 10.1006/jmbi.2000.4010.
Jedd G, Richardson C, Litt R, Segev N: The Ypt1 GTPase is essential for the first two steps of the yeast secretory pathway. J Cell Biol. 1995, 131 (3): 583-590. 10.1083/jcb.131.3.583.
Jedd G, Mulholland J, Segev N: Two new Ypt GTPases are required for exit from the yeast trans-Golgi compartment. J Cell Biol. 1997, 137 (3): 563-580. 10.1083/jcb.137.3.563.
Haubruck H, Prange R, Vorgias C, Gallwitz D: The ras-related mouse ypt1 protein can functionally replace the YPT1 gene product in yeast. Embo J. 1989, 8 (5): 1427-1432.
Wilcke M, Johannes L, Galli T, Mayau V, Goud B, Salamero J: Rab11 regulates the compartmentalization of early endosomes required for efficient transport from early endosomes to the trans-golgi network. J Cell Biol. 2000, 151 (6): 1207-1220. 10.1083/jcb.151.6.1207.
Chen W, Feng Y, Chen D, Wandinger-Ness A: Rab11 is required for trans-golgi network-to-plasma membrane transport and a preferential target for GDP dissociation inhibitor. Mol Biol Cell. 1998, 9 (11): 3241-3257.
Wang W, Sacher M, Ferro-Novick S: TRAPP stimulates guanine nucleotide exchange on Ypt1p. J Cell Biol. 2000, 151 (2): 289-296. 10.1083/jcb.151.2.289.
Jones S, Newman C, Liu F, Segev N: The TRAPP complex is a nucleotide exchanger for Ypt1 and Ypt31/32. Mol Biol Cell. 2000, 11 (12): 4403-4411.
Sacher M, Barrowman J, Wang W, Horecka J, Zhang Y, Pypaert M, Ferro-Novick S: TRAPP I implicated in the specificity of tethering in ER-to-Golgi transport. Mol Cell. 2001, 7 (2): 433-442. 10.1016/S1097-2765(01)00190-3.
Morozova N, Liang Y, Tokarev AA, Chen SH, Cox R, Andrejic J, Lipatova Z, Sciorra VA, Emr SD, Segev N: TRAPPII subunits are required for the specificity switch of a Ypt-Rab GEF. Nat Cell Biol. 2006, 8 (11): 1263-1269. 10.1038/ncb1489.
Kim YG, Raunser S, Munger C, Wagner J, Song YL, Cygler M, Walz T, Oh BH, Sacher M: The architecture of the multisubunit TRAPP I complex suggests a model for vesicle tethering. Cell. 2006, 127 (4): 817-830. 10.1016/j.cell.2006.09.029.
Cai H, Zhang Y, Pypaert M, Walker L, Ferro-Novick S: Mutants in trs120 disrupt traffic from the early endosome to the late Golgi. J Cell Biol. 2005, 171 (5): 823-833. 10.1083/jcb.200505145.
Loh E, Peter F, Subramaniam VN, Hong W: Mammalian Bet3 functions as a cytosolic factor participating in transport from the ER to the Golgi apparatus. J Cell Sci. 2005, 118 (Pt 6): 1209-1222. 10.1242/jcs.01723.
Gavin AC, Bosche M, Krause R, Grandi P, Marzioch M, Bauer A, Schultz J, Rick JM, Michon AM, Cruciat CM, Remor M, Hofert C, Schelder M, Brajenovic M, Ruffner H, Merino A, Klein K, Hudak M, Dickson D, Rudi T, Gnau V, Bauch A, Bastuck S, Huhse B, Leutwein C, Heurtier MA, Copley RR, Edelmann A, Querfurth E, Rybin V, Drewes G, Raida M, Bouwmeester T, Bork P, Seraphin B, Kuster B, Neubauer G, Superti-Furga G: Functional organization of the yeast proteome by systematic analysis of protein complexes. Nature. 2002, 415 (6868): 141-147. 10.1038/415141a.
Hu WH, Pendergast JS, Mo XM, Brambilla R, Bracchi-Ricard V, Li F, Walters WM, Blits B, He L, Schaal SM, Bethea JR: NIBP, a novel NIK and IKK(beta)-binding protein that enhances NF-(kappa)B activation. J Biol Chem. 2005, 280 (32): 29233-29241. 10.1074/jbc.M501670200.
Lafreniere RG, Kibar Z, Rochefort DL, Han FY, Fon EA, Dube MP, Kang X, Baird S, Korneluk RG, Rommens JM, Rouleau GA: Genomic structure of the human GT334 (EHOC-1) gene mapping to 21q22.3. Gene. 1997, 198 (1-2): 313-321. 10.1016/S0378-1119(97)00333-8.
Nagamine K, Kudoh J, Kawasaki K, Minoshima S, Asakawa S, Ito F, Shimizu N: Genomic organization and complete nucleotide sequence of the TMEM1 gene on human chromosome 21q22.3. Biochem Biophys Res Commun. 1997, 235 (1): 185-190. 10.1006/bbrc.1997.6758.
Yamakawa K, Mitchell S, Hubert R, Chen XN, Colbern S, Huo YK, Gadomski C, Kim UJ, Korenberg JR: Isolation and characterization of a candidate gene for progressive myoclonus epilepsy on 21q22.3. Hum Mol Genet. 1995, 4 (4): 709-716. 10.1093/hmg/4.4.709.
Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, Chen Y, Cheng X, Chua G, Friesen H, Goldberg DS, Haynes J, Humphries C, He G, Hussein S, Ke L, Krogan N, Li Z, Levinson JN, Lu H, Menard P, Munyana C, Parsons AB, Ryan O, Tonikian R, Roberts T, Sdicu AM, Shapiro J, Sheikh B, Suter B, Wong SL, Zhang LV, Zhu H, Burd CG, Munro S, Sander C, Rine J, Greenblatt J, Peter M, Bretscher A, Bell G, Roth FP, Brown GW, Andrews B, Bussey H, Boone C: Global mapping of the yeast genetic interaction network. Science. 2004, 303 (5659): 808-813. 10.1126/science.1091317.
Sciorra VA, Audhya A, Parsons AB, Segev N, Boone C, Emr SD: Synthetic genetic array analysis of the PtdIns 4-kinase Pik1p identifies components in a Golgi-specific Ypt31/rab-GTPase signaling pathway. Mol Biol Cell. 2005, 16 (2): 776-793. 10.1091/mbc.E04-08-0700.
Brown JL, Kossaczka Z, Jiang B, Bussey H: A mutational analysis of killer toxin resistance in Saccharomyces cerevisiae identifies new genes involved in cell wall (1-->6)-beta-glucan synthesis. Genetics. 1993, 133 (4): 837-849.
Adl SM, Simpson AG, Farmer MA, Andersen RA, Anderson OR, Barta JR, Bowser SS, Brugerolle G, Fensome RA, Fredericq S, James TY, Karpov S, Kugrens P, Krug J, Lane CE, Lewis LA, Lodge J, Lynn DH, Mann DG, McCourt RM, Mendoza L, Moestrup O, Mozley-Standridge SE, Nerad TA, Shearer CA, Smirnov AV, Spiegel FW, Taylor MF: The new higher level classification of eukaryotes with emphasis on the taxonomy of protists. J Eukaryot Microbiol. 2005, 52 (5): 399-451. 10.1111/j.1550-7408.2005.00053.x.
Soding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21 (7): 951-960. 10.1093/bioinformatics/bti125.
Letunic I, Goodstadt L, Dickens NJ, Doerks T, Schultz J, Mott R, Ciccarelli F, Copley RR, Ponting CP, Bork P: Recent improvements to the SMART domain-based sequence annotation resource. Nucleic Acids Res. 2002, 30 (1): 242-244. 10.1093/nar/30.1.242.
Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci U S A. 1998, 95 (11): 5857-5864. 10.1073/pnas.95.11.5857.
Jones S, Richardson CJ, Litt RJ, Segev N: Identification of regulators for Ypt1 GTPase nucleotide cycling. Mol Biol Cell. 1998, 9 (10): 2819-2837.
Sacher M, Barrowman J, Schieltz D, Yates JR, Ferro-Novick S: Identification and characterization of five new subunits of TRAPP. Eur J Cell Biol. 2000, 79 (2): 71-80. 10.1078/S0171-9335(04)70009-6.
Dolinski K, Balakrishnan R, Christie KR, Costanzo MC, Dwight SS, Engel SR, Fisk DG, Hirschman JE, Hong EL, Issel-Tarver L, Sethuraman A, Theesfeld CL, Binkley G, Lane C, Schroeder M, Dong S, Weng S, Andrada R, Botstein D, Cherry JM: Saccharomyces Genome Database. [http://www.yeastgenome.org/]
Mitaku S, Hirokawa T, Tsuji T: Amphiphilicity index of polar amino acids as an aid in the characterization of amino acid preference at membrane-water interfaces. Bioinformatics. 2002, 18 (4): 608-616. 10.1093/bioinformatics/18.4.608.
Ikeda M, Arai M, Lao DM, Shimizu T: Transmembrane topology prediction methods: a re-assessment and improvement by a consensus method using a dataset of experimentally-characterized transmembrane topologies. In Silico Biol. 2002, 2 (1): 19-33.
Segev N: Ypt and Rab GTPases: insight into functions through novel interactions. Curr Opin Cell Biol. 2001, 13 (4): 500-511. 10.1016/S0955-0674(00)00242-8.
Proteome BioKnowledge Library. [http://www.incyte.com/sequence/proteome/index.shtml]
Gedeon AK, Colley A, Jamieson R, Thompson EM, Rogers J, Sillence D, Tiller GE, Mulley JC, Gecz J: Identification of the gene (SEDL) causing X-linked spondyloepiphyseal dysplasia tarda. Nat Genet. 1999, 22 (4): 400-404. 10.1038/11976.
Muenke M, Bone LJ, Mitchell HF, Hart I, Walton K, Hall-Johnson K, Ippel EF, Dietz-Band J, Kvaloy K, Fan CM: Physical mapping of the holoprosencephaly critical region in 21q22.3, exclusion of SIM2 as a candidate gene for holoprosencephaly, and mapping of SIM2 to a region of chromosome 21 important for Down syndrome. Am J Hum Genet. 1995, 57 (5): 1074-1079.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2003, 31 (1): 23-27. 10.1093/nar/gkg057.
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. 10.1093/nar/22.22.4673.
Swofford DL: Phylogenetic Analysis Using Parsimony (*and Other Methods). 2003, Sunderland, Massachusetts , Sinauer Associates, 4.0
Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33 (Web Server issue): W557-9. 10.1093/nar/gki352.
Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12 (4): 357-358.
Schroder S, Schimmoller F, Singer-Kruger B, Riezman H: The Golgi-localization of yeast Emp47p depends on its di-lysine motif but is not affected by the ret1-1 mutation in alpha-COP. J Cell Biol. 1995, 131 (4): 895-912. 10.1083/jcb.131.4.895.
Rose M, Winston F, Heiter P: Methods in Yeast Genetics. 1988, Cold Spring Harbor, NY , Cold Spring Harbor Laboratory
Chen SH, Chen S, Tokarev AA, Liu F, Jedd G, Segev N: Ypt31/32 GTPases and their novel F-box effector protein Rcy1 regulate protein recycling. Mol Biol Cell. 2005, 16 (1): 178-192. 10.1091/mbc.E04-03-0258.
Nothwehr SF, Hindes AE: The yeast VPS5/GRD2 gene encodes a sorting nexin-1-like protein required for localizing membrane proteins to the late Golgi. J Cell Sci. 1997, 110 ( Pt 9): 1063-1072.
Cox R, Chen SH, Yoo E, Segev N: Conservation of the TRAPP II-specific subunits of a Ypt/Rab exchanger complex. [http://www.uic.edu/~nava/papers/TRAPPII]
Barton GJ: ALSCRIPT: a tool to format multiple sequence alignments. Protein Eng. 1993, 6 (1): 37-40. 10.1093/protein/6.1.37.
Kummel D, Muller JJ, Roske Y, Misselwitz R, Bussow K, Heinemann U: The structure of the TRAPP subunit TPC6 suggests a model for a TRAPP subcomplex. EMBO Rep. 2005, 6 (8): 787-793. 10.1038/sj.embor.7400463.
We thank H. Riezman for anti-Emp47 antibody, and B. Kay for advice, regarding the presentation of the results. This research was supported by grant GM-45444 from National Institutes of Health to N. Segev.
RC performed all sequence analyses, SHC and EY carried the protein analyses, NS with the help of the other authors, was involved in project planning, data analysis, and writing. All authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Supplementary Methods, Figures and Tables. Supplementary Figures S1-S5 show whole-sequence multiple alignments of Trs120, Trs130, and Trs65 proteins as well as predicted secondary structures of these subunits by taxonomic group. Supplementary Tables S1-S4 detail TRAPP II-specific subunits found in each genome and the highly conserved (HC) amino acids found in each subunit. Supplementary Methods applying to these supplementary data are described. (PDF 974 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Cox, R., Chen, S.H., Yoo, E. et al. Conservation of the TRAPPII-specific subunits of a Ypt/Rab exchanger complex. BMC Evol Biol 7, 12 (2007) doi:10.1186/1471-2148-7-12
- Unicellular Eukaryote
- Mammalian Orthologs
- Multiple Alignment Analysis
- Yeast Cell Lysate
- Future Functional Analysis