|Home | About | Journals | Submit | Contact Us | Français|
The TOR kinase is a major regulator of growth in eukaryotes. Many components of the TOR pathway are implicated in cancer and metabolic diseases in humans. Analysis of the evolution of TOR and its pathway may provide fundamental insight into the evolution of growth regulation in eukaryotes and provide a practical framework on which experimental evidence can be compared between species. Here we performed phylogenetic analyses on the components of the TOR pathway and determined their point of invention. We find that the two TOR complexes and a large part of the TOR pathway originated before the Last Eukaryotic Common Ancestor and form a core to which new inputs have been added during animal evolution. In addition, we provide insight into how duplications and sub-functionalization of the S6K, RSK, SGK and PKB kinases shaped the complexity of the TOR pathway. In yeast we identify novel AGC kinases that are orthologous to the S6 kinase. These results demonstrate how a vital signaling pathway can be both highly conserved and flexible in eukaryotes.
The online version of this article (doi:10.1007/s00239-011-9469-9) contains supplementary material, which is available to authorized users.
The target of rapamycin (TOR) is a major regulator of growth in eukaryotes (Thomas and Hall 1997). It integrates both intracellular signals that depend on nutrient availability, and extracellular signals, such as growth factors (Goberdhan et al. 2009; Oldham and Hafen 2003; Wullschleger et al. 2006). Therefore, dysfunction of TOR or other proteins in the TOR pathway is involved in organismal and cancer development in mammals (Menon and Manning 2008). The TOR protein is a kinase that participates in two distinct protein complexes. The TOR Complex 1 (TORC1) promotes translation by phosphorylating the S6 subunit of the ribosomal complex via ribosomal S6 kinase (S6K) and by phosphorylating the eukaryotic translation initiation factor 4E binding protein (4E-BP), which then dissociates from eIF4E allowing translation (Mahfouz et al. 2006; Tee and Blenis 2005). TORC1 is activated by the Ras homolog enriched in brain (Rheb) G-protein, which in turn is regulated by the tuberous sclerosis complex (TSC1/2 complex) (Inoki et al. 2003; Zhang et al. 2003). The best characterized activation route of TORC1 in animals is the insulin signaling pathway, which includes well characterized oncogenes and tumor suppressors such as phosphoinositide-3-kinase (PI3K), phosphoinositide-dependent kinase 1 (PDK1), phosphatase and tensin homolog (PTEN) and protein kinase B (PKB) (also known as AKT) (Engelman et al. 2006; Ma and Blenis 2009; Menon and Manning 2008; Shaw and Cantley 2006). It is not fully understood how TORC1 is regulated by nutrient levels (Avruch et al. 2009), although it is known that Rag GTPases play an important role in amino acid regulated TORC1 activity (Kim et al. 2008; Sancak et al. 2008). In addition S6K lowers PKB activity (Harrington et al. 2004; Kockel et al. 2010), creating a negative feedback by which TORC1 inhibits its own activation.
The TOR Complex 2 (TORC2) regulates cytoskeleton rearrangement in response to growth (Jacinto et al. 2004; Loewith et al. 2002). The regulation of TORC2 is mostly independent from TORC1 but involves PI3K signaling, activation by the ribosome (Zinzalla et al. 2011) and Rab GTPases (Tatebe and Shiozaki 2010). TORC2 positively regulates TORC1 via PKB (Cybulski and Hall 2009).
TOR is a conserved kinase and has been functionally characterized in animals (Hall 2008), fungi (yeast) (Lorberg and Hall 2004) and plants (Deprost et al. 2007). However, not all of the TOR complex subunits or TOR pathway components are equally conserved. For instance, the yeast Saccharomyces cerevisiae has retained the Rheb G-protein (Urano 2000), but lacks the Rheb regulatory genes Tuberous sclerosis 1 (TSC1) and Tuberous sclerosis 2 (TSC2), whereas Schizosaccharomyces pombe has retained all three genes (Matsumoto et al. 2002). Plants appear to lack Rheb and the TSC1/2 complex altogether (Díaz-Troya et al. 2008; Vernoud et al. 2003) despite the presence of TORC1 (Deprost et al. 2007; Menand et al. 2002).
To obtain novel insights into the evolution of the TOR pathway we performed phylogenetic analyses on the components of the TOR pathway, starting from the mammalian TOR pathway. We are specifically interested in the distribution of each subunit over all eukaryotes and the point of invention of single genes or even modules in evolution. To that end we searched in 64 diverse eukaryotic genomes for homologous and orthologous genes for components of the TOR pathway. The recent publication of genomes of highly divergent eukaryotic organisms such as Naegleria gruberi (Fritz-Laylin et al. 2010) now permits us to make detailed phylogenetic analyses of genes that play a key role in the TOR pathway and determine the earliest possible point of invention in evolution for each subunit. For a full overview of the species included in our analysis and how they are related to each other see Figure S4. We analyzed each evolutionary or phylogenetic profile in light of distinct components of the pathway (e.g., the TOR complexes, Rheb regulation) and the pathway as a whole.
Importantly, we found that TOR, all subunits of TORC1 and 2 and a large part of the TOR pathway components form an evolutionary core. New regulatory inputs, such as insulin and TNFα signaling, have been added to the core TOR pathway during animal evolution. We show that TORC1 and TORC2 appear to behave as independent evolutionary modules, even though the majority of the subunits are shared between the two complexes. We infer the presence of a large common evolutionary core, including Rheb and TSC2, in the Last Eukaryotic Common Ancestor (LECA), the last ancestral eukaryote that gave rise to all current eukaryotic species. We provide an in-depth phylogenetic analysis of TOR, Rheb, TSC2 and translationally controlled tumor protein (TCTP) as well as other components of the TOR signaling pathway. In addition, we reveal the remarkable role of the duplications of an ancestral AGC kinase that gave rise to the S6K, ribosomal protein S6 kinase (RSK), PKB and serum/glucocorticoid regulated kinase (SGK) kinases, in the increasing complexity of the TOR pathway. We conclude that a vital signaling pathway can be both highly conserved and flexible in eukaryotes.
TOR functions as part of two distinct protein complexes: the TORC1 and TORC2 complexes (Loewith et al. 2002). In mammalian cells TORC1 contains mTOR, Lethal with Sec Thirteen (LST8) and Raptor, while TORC2 contains mTOR, LST8, Raptor independent companion of TOR (Rictor) and stress-activated MAP kinase interacting protein (SIN1) (Fig. 1). TOR and LST8 are both present in genomes in all major eukaryotic lineages and therefore form the evolutionary core of the TOR complexes (Fig. 1). In addition, we also observe that TOR and LST8 co-occur with either Raptor (TORC1) or Rictor (TORC2) or both, indicating that both TOR complexes are old and were likely present in LECA, the common ancestor of all current eukaryotic species.
We find TORC1 together with TORC2 in all major lineages, except plants, which possess only TORC1 (Fig. 1). Interestingly we detect TORC2, but not TORC1 in the ciliates Tetrahymena thermophila and Paramecium tetraurelia. It therefore appears that the two distinct TOR complexes are decoupled in evolution as either one can be lost while the other is maintained.
We do not detect any of the TOR complex subunits in the microsporidium Encephalitozoon cuniculi and the apicomplexa Plasmodium falciparum, Cryptosporidium parvum and Theileria parva indicating at least two independent loss events for TOR signaling in eukaryotes. All four species are intracellular parasites and have reduced genomes and cellular structure. Host–parasite interactions might have replaced the function of TOR in these organisms as growth of the parasite is directly linked to conditions in the host cell.
Distinct protein complexes that share subunits (i.e., hyperlinks) can provide a selective reason for maintaining duplicate copies of these shared subunits (Shevchenko et al. 2008). Therefore, we could expect to find duplications of TOR in some species, resulting in a dedicated TOR for each of the two TOR complexes. Indeed we find duplications of TOR in Saccharomyces cerevisiae, Candida glabrata, Schizosaccharomyces pombe, Batrachochytrium dendrobatidis, Populus trichocarpa, Emiliania huxleyi, Trypanosoma brucei, Leishmania major, Phytophthora infestans and Phytophthora sojae (see Figure S1). The duplication of TOR in S. cerevisiae and C. glabrata likely originates from the Whole Genome Duplication event. In S. cerevisiae both TOR1 and TOR2 can be part of TORC1, while TOR2 is specific for the TORC2 (Loewith et al. 2002), which indicates that the two TOR duplicates are not completely functionally diverged. The TOR duplication in S. pombe and B. dendrobatidis are lineage-specific duplications and occurred independently from each other and from the duplications in S. cerevisiae and C. glabrata. While we have no functional descriptions for B. dendrobatidis, it has been shown for S. pombe, that TOR1 and TOR2 function as part of TORC1 (Hartmuth and Petersen 2009), while TOR1 is specific for TORC2 (Otsubo and Yamamato 2008). Note, that the naming of TOR1 and TOR2 in S. cerevisiae and S. pombe can cause confusion as the genes resulted from independent duplication events, the naming of TOR1 and TOR2 in both yeasts does not reflect one-to-one orthologous relationships but is based on order of discovery (Soulard et al. 2009). Surprisingly, LST8 has not been duplicated in any of the species examined. This raises an interesting question: why do the two hyperlinks TOR and LST8 behave differently in evolution? We hypothesize that duplication and subsequent functional divergence of LST8 may have implications for the structural integrity of the two TOR complexes, while minor modifications to the TOR duplicate genes increased functional divergence without compromising complex stability.
The Rheb G-protein is one of the major regulators of TOR activity in animals and directly regulates the activity of TORC1 but not TORC2 (Cybulski and Hall 2009). Rheb is a Ras-like small GTPase and the sequences of small GTPases are highly conserved (Colicelli 2004; Wennerberg et al. 2005). We previously reconstructed the phylogeny of the Ras-like small GTPases (van Dam et al. 2011). From this phylogeny we identified Rheb orthologs (see “Methods” section; Figure S2) and derived a phylogenetic profile of orthologs (Fig. 1). We identified Rheb orthologs in all animals and fungi (except in C. glabrata, Eremothecium gossypii and E. cuniculi). Additionally we identify orthologs in distantly related organisms such as diatoms, oomycetes, the amoeba Dictyostelium discoideum, the heterolobosida Naegleria gruberi and the red algae Cyanidioschyzon merolae. This is the first time, to our knowledge, that a G-protein belonging to the Ras-like subfamily of small GTPases has been identified in the Archaeplastida (i.e., plants and red and green algae).
The identification of Rheb orthologs in distantly related species strongly suggests that Rheb originated in or before LECA. We observe Rheb orthologs in species that also contain TORC1, which indicates that the regulation of TOR by Rheb is strongly conserved. However, we do not observe the opposite, e.g., species that have TORC1 do not necessarily have a Rheb ortholog. The most notable of these species are the green algae and plants, but also the yeasts C. glabrata and E. gossypii. The presence of a Rheb ortholog in C. merolae indicates a loss of Rheb in the ancestor of the green algae and plants.
Next we investigated the phylogeny of the only known regulator of Rheb, the GTPase activating protein (GAP) TSC2 (data on TSC2 is also presented in (van Dam et al. 2011) in light of Ras signaling and other RapGAP domain containing proteins). TSC2 integrates many inputs such as MAPK/Ras signaling via RSK1 (Ma et al. 2005; Roux et al. 2004), Wnt signaling via GSK3β (Inoki et al. 2006) and insulin signaling via PI3K and PKB (Zoncu et al. 2011) in animals. The TSC2 GAP domain occurs in all major eukaryotic super groups, except the excavates, indicating it much older than previously suggested (Serfontein et al. 2010) and likely originated in or before the LECA. We find the TSC2 GAP domain orthologs in species that also have a Rheb ortholog, including the red algae C. merolae mentioned above. The only exception is the ciliate T. thermophila. Therefore, we predict that the TSC2 GAP orthologs will regulate the Rheb orthologs in D. discoideum, C. merolae and the Phytophthora species.
Furthermore, we find that while TSC1 orthologs are always observed together with TSC2 orthologs in the same genomes, TSC2 can be found on its own in additional eukaryotic species (Fig. 1). Interestingly in some of these species (D. discoideum, C. merolae, P. infestans, P. sojae, Phaeodactylum tricornutum) we were able to identify the GAP domain but not the Tuberin domain that is necessary to dimerize with TSC1. Strikingly, we find TSC1 orthologs in animals and fungi, which is the same phylogenetic distribution as the tuberin domain of TSC2. Therefore, it is likely that TSC1 itself and the ability of TSC2 to dimerize with TSC1 via the tuberin domain are inventions in the Opisthokont ancestor (i.e., in the animal and fungal ancestor).
The absence of the TSC1/2 complex in C. elegans and S. cerevisiae (Fig. 1) suggests that caution should be taken when comparing regulatory mechanisms of TOR between these two and other species. For instance, regulatory mechanisms for Rheb and TOR discovered in animals, such as in D. melanogaster do not necessarily hold for C. elegans and vice versa.
TCTP (also known as TPT1) has been reported to be the Guanine Exchange Factor (GEF) for Rheb (Hsu et al. 2007) but this function has been debated by us and others (Rehmann et al. 2008; Wang et al. 2008). Hence it is interesting to study the phylogenetic profiles of TCTP and TSC2 and compare them to Rheb (Fig. 1). We constructed the phylogenetic profile of TCTP and found that TCTP is present in nearly all of the eukaryotic genomes we studied including plants, chromalveolata and excavata. Interestingly we observe eukaryotic species which have a TCTP ortholog but do not have the Rheb G-protein (i.e., green plants and algae, apicomplexa, ciliates, C. glabrata, E. gossypii, see Fig. 1) and vice versa (i.e., the choanoflagellate M. brevicollis and the diatom P. tricornutum). In addition we observe in the apicomplexa (i.e., a group of unicellular intracellular parasites including Plasmodium falciparum), that even though they contain a TCTP ortholog, they do not only lack Rheb but also TOR.
In a previous study on Ras Guanine Exchange Factors (RasGEFs) we observed a strong evolutionary link between the RasGEFs and their respective Ras-like GTPases (van Dam et al. 2009). The RasGEF functional domain, the CDC25 homology domain, and the Ras, Rap and Ral GTPases were found to be present or absent together in a diverse set of eukaryotic genomes. In contrast we observe no such strong presence/absence pattern for TCTP and Rheb. It therefore seems that there is no evolutionary linkage between TCTP and Rheb.
Further doubts about TCTP’s GEF activity for Rheb arise from experimental evidence of TCTP function in Arabidopsisthaliana when put in an evolutionary framework. Berkowitz et al. (2008) studied the function of the ortholog of TCTP in A. thaliana. They found that TCTP acts as an important regulator of growth, and implied that TCTP functions in TOR activity, which resembles the situation in animals (Hsu et al. 2007). The similarity between plant TCTP and animal TCTP is interesting because Arabidopsis does not have a Rheb ortholog (Vernoud et al. 2003 and this study). Berkowitz et al. postulate that Arabidopsis TCTP regulates another GTPase (either a Rhop- or Rab-like G-protein) which might function in an equivalent way to Rheb. We are of the opinion that the results of Berkowitz et al., that TCTP functions as a positive regulator on TOR activity in A. thaliana in the absence of Rheb, indicate that TCTP might in fact not be a RhebGEF, but instead suggest that TCTP regulates TOR via an alternative route.
We extended our phylogenetic study to include upstream and downstream components of the mammalian TOR pathway and thereby put Rheb and TORC in a wider biological context. We focus on the mammalian TOR pathway because TOR signaling in animals and particularly mammals is the most extensively studied intact TOR pathway (e.g., not lacking key components such as the TSC1/2 complex in S. cerevisiae or TORC2 in A. thaliana). The TOR pathway was assembled from literature to reflect current consensus. For each protein we constructed the phylogenetic profile and determined the point of invention (e.g., the age of a protein. See Table S2 for all phylogenetic profiles). We depict the age of each protein along the metazoan natural history towards LECA in the representation of the TOR pathway (Fig. 2).
Recently, Serfontein et al. (2010) published an evolutionary survey on the components of the TOR pathway in a representative selection of eukaryotic genomes. Our results concerning the evolution of the TSC1/TSC2-TOR pathway underline some observations made by Serfontein and coworkers but differ considerably in others. While Serfontein and coworkers find that the evolutionary “core” of the pathway that was present in LECA consisted out of the TORC1 complex (TOR, LST8, Raptor), AMP-activated protein kinase (AMPK), PI3K and PTEN and S6K, we find that TORC2 (TOR, LST8, Rictor, SIN1), Rheb, TSC2, PDK1 and the remaining AGC kinases PKB, RSK and SGK are also part of this evolutionary core. We show that Rheb and TSC2 have not been “bolted on” in evolution but are in fact part of the evolutionary core that originated in or before LECA. We do, however, observe other regulatory elements of the TOR pathway that have been added at a later evolutionary stage. We made two observations concerning the evolution of new TOR signaling regulation.
The first observation is that the regulation of TOR activity by insulin is an animal-specific addition to the pathway. The second observation is that the more recently invented TSC1 introduces novel regulatory input onto the Rheb-TOR cascade (i.e., FIP200 and IKKβ/TNFα). These two observations suggest that although TOR signaling is highly conserved within eukaryotes, it is also flexible enough to accept new inputs and can be adapted to suit new environments (e.g., multicellular tissues).
The animal-specific addition of insulin signaling to the TOR pathway could be an adaptation to account for cellular growth in a multicellular environment. The growth and divisions of individual cells need to be regulated system wide and insulin might very well play the role of system-wide growth control via TOR. Stimulation by insulin has been a fruitful method to investigate TOR signaling in animal models, but as a model of TOR activation and regulation it is distinct from our knowledge of TOR in non-animal model species. It would be interesting to investigate if or in what way TOR is regulated by system-wide growth control (other than nutrient availability) in plants, a group of organisms that have separately acquired multicellularity.
Ma et al. (Ma and Blenis 2009) suggested that PRAS40 might represent a conserved PKB regulation route to TORC1 while TSC1/2 represented a newer additional pathway in higher eukaryotes because TSC1 and TSC2 were previously not found in lower eukaryotes like S. cerevisiae. However, due to the availability of many newly sequenced eukaryotic genomes and with more extensive phylogenetic profiling as done here, it becomes evident that PRAS40 represents an animal-specific additional route of PKB-TORC1 activation and that the Rheb/TSC1/TSC2 route is in fact much older.
Serfontein and coworkers (Serfontein et al. 2010) have detected PKB orthologs only in animals, amoebozoa and excavate species, but not in plants and chromalveolates. However, they have included only those sequences that included both the kinase domain and a Pleckstrin Homology (PH) domain that is characteristic of the animal PKB. However, the PH domain is known to be very promiscuous in eukaryotes (Basu et al. 2008) and is therefore not suitable to use as a restriction criterion for orthology.
We find that the AGC kinases PKB and SGK are paralogs that have arisen from a duplication event in the animal and choanoflagellar ancestor, but we also find PKB/SGK orthologs in chromalveolates and excavates. PKB and SGK share common ancestry with the AGC kinases S6K and RSK, which are also involved in TOR signaling. This suggests that the AGC kinases and duplication events play a striking and complex evolutionary role in the TOR pathway. We therefore focused on the evolution of the AGC kinases.
In the mTOR pathway, the AGC family kinases S6K, RSK1, PKB and SGK1 are located both upstream and downstream of TOR. S6K and RSK1 arose from a duplication event in the ancestor of animals and fungi (Opisthokont ancestor) while PKB and SGK1 arose from a duplication event in the ancestor of animals and the closely related choanoflagellate Monosiga brevicollis (Shalchian-Tabrizi et al. 2008) (see Fig. 3a; Figure S3). The S6K-RSK and PKB-SGK ancestral genes themselves have arisen from a gene duplication in or before LECA. The evolutionary relation between S6K, RSK1, PKB and SGK1 make it uniquely possible to reconstruct the evolution of their regulatory interactions within the TOR pathway. In Fig. 3 we have reconstructed the TOR pathway at several points in evolution based on events in the evolution of the AGC kinases and experimental characterization of orthologous genes in H. sapiens, S. cerevisiae, S. pombe and A. thaliana.
The PKB and SGK genes duplicated from a single ancestral gene in the filozoan ancestor, i.e., the ancestor of animals and choanoflagellates (see Fig. 3a; Figure S3). However, additional new components have been invented specifically in animals, such as the insulin and TNFα signaling pathways. Therefore, we reconstructed the TOR pathway in the metazoan ancestor. PKB and SGK are both activated by PDK1 (Alessi et al. 1997; Burgering and Coffer 1995; Kobayashi and Cohen 1999; Park et al. 1999) and the TORC2 complex (Hresko and Mueckler 2005; Jones et al. 2009; Sarbassov et al. 2005; Soukas et al. 2009). In addition, PKB is activated by translocation to the plasma membrane upon PI3K activation (Stokoe 1997) (Fig. 3b). It is therefore very likely that the PKB/SGK ancestral protein was also activated by PDK1 and TORC2 ancestral proteins. PI3K-dependent translocation of PKB is likely a new function acquired by PKB.
PKB inhibits the TSC1/2 complex by phosphorylating TSC2 (Dan et al. 2002; Potter et al. 2002) and inhibits FOXO transcription factors by directly phosphorylating them (Burgering and Kops 2002). SGK also inhibits FOXO by phosphorylation (Brunet et al. 2001) but has not been reported to phosphorylate TSC2. We can partly derive the ancestral functions by comparing PKB and SGK functions to the gene functions of the co-orthologous genes in S. pombe and S. cerevisiae.
In S. pombe and S. cerevisiae there are three co-orthologous genes to both PKB and SGK (Sck1, Sck2 and Gad8 in S. pombe, YPK1, YPK2 and SCH9 in S. cerevisiae, see Figure S3). The S. pombe and S. cerevisiae PKB/SGK1 orthologs Sck1, Sck2, Gad8 and SCH9 are implied to have function in oxidative stress responses and aging (Chen and Runge 2009; Ikeda et al. 2008), similar to PKB and SGK1 in animals. Therefore, the role of PKB and SGK in regulating longevity is conserved and likely an ancestral function in the Opisthokont ancestor (Fig. 3c).
In contrast to stress response and aging, the origin of TSC2 phosphorylation by PKB is not immediately apparent. In S. cerevisiae TSC2 has been lost and we have been unable to find any references that implicate the S. pombe PKB/SGK orthologs Sck1, Sck2 in growth regulation via TOR (Sck1/2) or that Gad8 has been associated with the TSC1/2 complex. This makes it difficult to determine if TSC2 phosphorylation is an ancestral function or whether it has been specifically acquired by PKB. Nevertheless, there are similarities in function of PKB with RSK, and we can therefore reconstruct the ancestral function of PKB and SGK by comparing their functions to their paralogs RSK and S6K.
Similar to PKB, RSK also inhibits the TSC1/2 complex by phosphorylating TSC2 in mammals (Roux et al. 2004). Because the GAP domain of TSC2 is conserved throughout the eukaryotic lineage, the most plausible scenario is that PKB and RSK inhibition of TSC2 is an ancestral function of the PKB-SGK-S6K-RSK ancestral gene (henceforth we will refer to this ancestor as the ancestral AGC kinase for brevity). In this scenario the TSC2 regulation is an ancestral function maintained by PKB and lost by SGK (Fig. 3b).
Interestingly, similar to the filozoan PKB-SGK duplication event, fungi seem to have undergone a similar duplication event of the ancestral PKB-SGK kinase. YPK1 and YPK2 function in S. cerevisiae can be rescued by rat SGK, but not mouse PKB or rat S6K (Casamayor et al. 1999). Therefore, although the PKB-SGK and YPK-SCH9 duplication events in animals and fungi occurred independently, the resulting animal and fungal paralogs appear to have evolved in a functionally similar way.
The S6K and RSK genes duplicated from a single ancestral gene in the fungal and animal ancestor (Fig. 3c). Like PKB and SGK, both are regulated by PDK1 (Jensen 1999; Pullen et al. 1998). However, S6K and RSK are not regulated by the TORC2 complex (Jacinto and Lorberg 2008). S6K is regulated by the TORC1 complex instead (Burnett et al. 1998), while RSK is regulated via MAPK signaling (Carrière et al. 2008). While it is very likely that the S6K-RSK ancestral protein was regulated by the PDK1 ancestral protein, the regulation of S6K by TORC1 and RSK by MAPK initially obscures whether the ancestral kinase was either activated by TORC1 or MAPK or both.
We can infer whether the ancestral protein was activated by TORC1 and/or MAPK by inferring the ancestral function from experimental evidence for TOR signaling in the plant A. thaliana. RSK1 and S6K are co-orthologs to the A. thaliana S6K kinase. The S6K kinase of A. thaliana is regulated by the Arabidopsis TOR complex (Mahfouz et al. 2006), which has the same composition as TORC1 in other organisms, and we can therefore infer that the S6K-RSK ancestral protein was likely activated by TORC1. Thus S6K maintained the ancestral regulation by TORC1, i.e., RSK has lost the regulation by TORC1 (the dashed edge in Fig. 3c from TORC1 to RSK). Activation by MAPK of the ancestral protein cannot be inferred from the Arabidopsis S6K as there is no published link between Arabidopsis MAPK and Arabidopsis TOR signaling but we cannot exclude loss of this function in plants (the dashed edge from MAPK in Fig. 3c). Previously we deduced that RSK regulation of TORC1 activity via TSC2 is likely an ancestral function from the symmetry with PKB. Therefore, TSC2 regulation by RSK is an ancestral function that has been lost in S6K (the dashed edge in Fig. 3c from S6K to TORC1).
Intriguingly, both S. pombe and S. cerevisiae have a one-to-one ortholog to S6K that has not been fully characterized yet, or only superficially (see Fig. 3a; Figure S3). The S. cerevisiae locus YBR028C codes for a kinase, but has not been reported in the literature. The S6K ortholog in S. pombe, psk1, has been reported to be involved in phenylarsine oxide resistance and disruption of the psk1 gene did not result in growth defects (Mukai 1995). Instead, the S6K-like cellular function in S. cerevisiae has been ascribed to SCH9 (Urban et al. 2007), which according to our analysis is an ortholog of mammalian PKB and SGK. SCH9 shows that the AGC kinases are capable of performing cellular functions that have been ascribed to their paralogs, possibly increasing the complexity of the roles the AGC kinases play in TOR signaling. Given that the poorly characterized YBR028C gene in S. cerevisiae is clearly orthologous to S6K, we suggest that there might be a substantial role for this gene to be uncovered in TOR signaling.
Above we have described the ancestral states of the S6K-RSK and PKB-SGK ancestral genes. We observe symmetric functions between the two ancestral genes and therefore we are able to (partly) reconstruct the functions of the ancestral AGC kinase (e.g., the ancestral gene of S6K, RSK, PKB and SGK). All four kinases are regulated by PDK1 and we can therefore infer that the ancestral AGC kinase was also regulated by the PDK1 ancestral gene (Fig. 3d). Above, we also deduced that the S6K-RSK and PKB-SGK ancestral genes possibly regulated TORC1 activity via the TSC2 ancestral gene as it is a shared function of both RSK and PKB (Fig. 3d).
Phosphorylation of the S6 ribosomal subunit by the S6K-RSK ancestral kinase cannot be reconstructed beyond the LECA, because the PKB and SGK kinases do not share this function and we are therefore unable to determine if the S6 activation was a function of the ancestral AGC kinase that has been lost by the PKB-SGK ancestral gene or that it is an acquired function of the S6K-RSK ancestral gene (the dashed edge from S6K/RSK in Fig. 3d).
Our phylogenetic reconstruction suggests the existence of a TORC1 and a TORC2 complex in LECA that functions in conjunction with Rheb/TSC2 to activate at least two distinct AGC kinases (Fig. 3d). These AGC kinases arose from duplication and thus required (partially) independent regulation by TOR for their sub-functionalization (compare Fig. 3e and d). TORC1 and TORC2-specific proteins like Raptor and Rictor most likely contributed situation-specific activation of AGC kinases, i.e., determine context-relevant downstream outputs.
The positive feedback loop that emerges in the reconstruction of the ancestral TOR pathway in LECA (Fig. 3d, e) imposes a network structure that is undesirable as the feedback loop could easily result in a constitutively activated TOR and therefore uncontrolled growth. In animals negative feedback from S6K and TOR to PKB or more upstream elements in the insulin pathway have been documented (Findlay et al. 2005; Kockel et al. 2010; Manning 2004). It is very likely, given the importance of proper regulation of the TOR pathway that negative feedback routes were also in place in LECA. However, we are unable to reconstruct these negative feedbacks back to LECA. An explanation for this could be that the negative feedbacks in the TOR pathway have been subject to change in evolution. Subsequent duplications of the ancestral AGC kinase and sub-functionalization, however, might have provided the opportunity to increase the possibility for more precise or additional regulation of TOR activity.
The TOR pathway is a universal regulator of cell growth in eukaryote species. TOR is the master regulator and integrates many inputs such as growth signals and nutrient availability in a cell. We show that the two TOR complexes are highly conserved and were already present in the last common ancestor of all eukaryotes. We analyzed the origin and evolution of each subunit of the TOR complexes separately as well as other components of the TOR signaling pathway. We show that TORC1 and TORC2 behave as separate evolutionary modules that can be individually lost [i.e., loss of either Raptor (TORC1) or Rictor (TORC2) or both complexes as a whole]. We find that the TOR pathway has a conserved “core” to which new inputs have been added early in animal evolution, such as insulin and TNFα signaling. We also find evidence that the core itself has been extensively modified in evolution by duplications of ancestral AGC kinases that gave rise to S6K, RSK, SGK and PKB. The evolution of TOR and conversely the whole TOR pathway demonstrates that a vital signaling pathway can be both highly conserved and flexible in eukaryotes and can be adapted to fulfill changing requirements of growth regulation by eukaryotic organisms.
We acquired protein sequences (proteomes) of 64 divergent eukaryotic species from EnsEMBL (Hubbard et al. 2009), JGI, the Broad institute or their respective genome project sites. Proteomes filtered for “best model” were downloaded when available and longest sequences were selected in case of multiple transcripts per gene. We have selected a wide range of animal and fungal genomes as most research on TOR signaling is being done in either animal or fungal model organisms. We also included a wide range of genomes belonging to other major super groups, such as the Archeaplastida, Chromalveolates and Excavates, to be able to accurately time the origin of each TOR pathway component. For a full overview of genomes, source and version information see Table S1. For a species tree of the used genomes, see Figure S4.
Orthology was determined automatically by applying MCL (Van Dongen 2008) on InParanoid (Remm et al. 2001) species–species comparisons. We performed an all versus all BLAST (Altschul et al. 1990) on the whole set of 64 genomes with the options –p blastp –m 8 –v 1000000 –b 1000000. The BLAST results were then split into separate data files as required for InParanoid. InParanoid was applied to all possible combinations of the 64 genomes using the default settings (score_cutoff = 50, outgroup_cutoff = 50, seq_overlap_cutoff = 0.5, conf_cutoff = 0.05, group_overlap_cutoff = 0.5, grey_zone = 0). A matrix dataset for MCL was constructed from the InParanoid analyses by constructing edges between genes for each InParanoid cluster taking the lowest confidence value of both genes. Edges were only drawn between genes of different species. The MCL analysis was run with parameters –abc –I 1.5 –write-graph.
Phylogenetic profiles were constructed from the InParanoid-MCL clusters by determining which species were (not) represented in each cluster. The phylogenetic profiles of genes of interest were manually verified.
The sequences of the selected genomes were searched using the Pfam (Finn et al. 2008) HMM profile for the Ras family (Pfam accession PF00071.12, Pfam version 23) and hmmsearch of the HMMER package (Eddy 1998) version 2.3.2. All sequences with a bitscore larger than 0 were selected. Due to the high sequence similarity of Ras to other small GTPases many other small GTPases are included in this set. An alignment of all sequences was made using the MAFFT program (Katoh et al. 2002) with the –globalpair option. A neighbor joining tree was constructed using the Quicktree program (Howe et al. 2002). A sub-tree was selected which contained all Ras-like subfamily members but no other small GTPase. The sequences were gathered from the initial alignment as manual inspection of the alignment produced from the subset showed it was suboptimal to the initial alignment. Subsequently a phylogenetic tree was constructed over all Ras-like subfamily members using RAxML (Stamatakis et al. 2005) (-T 4 -x 488761235 -f a -N 1000 -m PROTGAMMAIWAG). From the resulting tree a sub-tree was selected that contained Rheb and as many genes from the searched species as possible without including other known small GTPases. All phylogenetic trees were visualized using Dendroscope (Huson et al. 2007). The data on Rheb phylogeny has also been presented in (van Dam et al. 2011).
Protein sequences belonging to the cluster of orthologous groups that contains S6K, RSK, PKB and SGK, were aligned using MAFFT (Katoh et al. 2002) with the –globalpair option. The resulting alignment was analyzed and a segment that showed high sequence similarity between all sequences was used for further analysis (positions 893–2050, corresponding to the kinase domain and PKC terminal domain). A phylogenetic tree, including bootstrap analysis, was constructed using RAxML (Stamatakis et al. 2005) (-T 8 -x 78382369 -f a -N 1000 -m PROTGAMMAIWAG).
TSC2 orthologs were initially identified using the automated orthology determination as described above. Closer examination revealed that the sequences only had the GAP domain sequence in common. Since the GAP domain of TSC2 belongs to a larger family of GAP domains, the Ran/RapGAP domain family, we used phylogenetic methods to faithfully determine true orthology based on the RapGAP domain sequences. We gathered RapGAP domain sequences from the sequence set by using a custom-made HMM model and hmmsearch of the HMMER package (Eddy 1998) version 2.3.2. The domain sequences were aligned using MAFFT (Katoh et al. 2002) –globalpair. A phylogenetic tree was built using RAxML (Stamatakis et al. 2005) (-T 6 -x 23421421 -f a -N 1000 -m PROTGAMMAIWAG). The data on TSC2 phylogeny has also been presented in (van Dam et al. 2011).
We want to thank Prof. Michael Hall and Thomas Sturgill for their suggestions and critical appraisal, Jos Boekhorst, Gabino Sanchez-Perez, Like Fokkens and Michael Seidl for their help in performing the analyses and support. The sequence data of selected genomes were produced by the US Department of Energy Joint Genome Institute http://www.jgi.doe.gov/, in collaboration with the user community, or the Fungal Genome Initiative of the Broad Institute. For a full overview of the genomes and references see Table S1. This work was supported by the BioRange program of the Netherlands Bioinformatics Centre (NBIC), which is supported by a BSIK grant through the Netherlands Genomics Initiative (NGI).
Open Access This article is distributed under the terms of the Creative Commons Attribution Noncommercial License which permits any noncommercial use, distribution, and reproduction in any medium, provided the original author(s) and source are credited.
Teunis J. P. van Dam, Phone: +31-243619774, Fax: +31-243619395, Email: moc.liamg@maDnav.P.J.sinueT.
Fried J. T. Zwartkruis, Email: ln.thcertucmu@siurktrawZ.T.J.G.
Johannes L. Bos, Email: ln.thcertucmu@soB.L.J.
Berend Snel, Email: ln.uu@lenS.B.