- Research article
- Open Access
Conservation, duplication, and loss of the Tor signaling pathway in the fungal kingdom
BMC Genomicsvolume 11, Article number: 510 (2010)
The nutrient-sensing Tor pathway governs cell growth and is conserved in nearly all eukaryotic organisms from unicellular yeasts to multicellular organisms, including humans. Tor is the target of the immunosuppressive drug rapamycin, which in complex with the prolyl isomerase FKBP12 inhibits Tor functions. Rapamycin is a gold standard drug for organ transplant recipients that was approved by the FDA in 1999 and is finding additional clinical indications as a chemotherapeutic and antiproliferative agent. Capitalizing on the plethora of recently sequenced genomes we have conducted comparative genomic studies to annotate the Tor pathway throughout the fungal kingdom and related unicellular opisthokonts, including Monosiga brevicollis, Salpingoeca rosetta, and Capsaspora owczarzaki.
Interestingly, the Tor signaling cascade is absent in three microsporidian species with available genome sequences, the only known instance of a eukaryotic group lacking this conserved pathway. The microsporidia are obligate intracellular pathogens with highly reduced genomes, and we hypothesize that they lost the Tor pathway as they adapted and streamlined their genomes for intracellular growth in a nutrient-rich environment. Two TOR paralogs are present in several fungal species as a result of either a whole genome duplication or independent gene/segmental duplication events. One such event was identified in the amphibian pathogen Batrachochytrium dendrobatidis, a chytrid responsible for worldwide global amphibian declines and extinctions.
The repeated independent duplications of the TOR gene in the fungal kingdom might reflect selective pressure acting upon this kinase that populates two proteinaceous complexes with different cellular roles. These comparative genomic analyses illustrate the evolutionary trajectory of a central nutrient-sensing cascade that enables diverse eukaryotic organisms to respond to their natural environments.
The nutrient-sensing t arget o f r apamycin (Tor) pathway is highly conserved among eukaryotes and governs several essential cellular processes including protein synthesis, ribosome biogenesis, autophagy, and cytoskeletal organization [1–3]. In the fungal kingdom, the Tor pathway has been best studied in the budding yeast Saccharomyces cerevisiae[2–4], the fission yeast Schizosaccharomyces pombe[5, 6], and the human pathogen Candida albicans[5–8]. However, little is known about Tor in basal fungal lineages, including the Zygomycota and Chytridiomycota.
In S. cerevisiae and S. pombe, two Tor paralogs form distinct complexes known as Tor Complex 1 (TORC1) and Tor Complex 2 (TORC2) [9–12], while in most other species, including humans, a single Tor protein can populate both complexes [11–14]. Interestingly, S. cerevisiae Tor2 can complement the loss of Tor1, but Tor1 cannot complement the loss of Tor2 [15, 16]. Two Tor paralogs have also been identified in a metazoan, the silkworm Bombyx mori and three Tor paralogs were identified in the trypanosomatid parasites Leishmania major and the related species Trypanosoma brucei, the first reported Tor triumvirates.
The ScTORC1 is sensitive to rapamycin and controls protein synthesis, mRNA synthesis and degradation, ribosome biogenesis, and autophagy. TORC2 is insensitive to rapamycin and is involved in the control of actin polarization and cell wall integrity [9, 16]. TORC1 consists of Tor1 or Tor2, Kog1 , Tco89, and Lst8 , while TORC2 contains Tor2, Lst8, Avo1, Avo2, and Avo3  (Figure 1A). Recently, the EGO-GTPase complex and its orthologs were shown to convey amino acid signals for TORC1 activation in yeast, insects, and mammals [23–26]. In S. cerevisiae, immediate effectors of TORC1 include the PP2A-like phosphatase Sit4  and the AGC kinase Sch9 . In S. pombe, Tor participates in other cellular functions including nutrient signaling , cell growth and differentiation , mitotic commitment , and sexual development  (Figure 1B).
The structure of the Tor proteins is remarkably conserved (Figure 2) and features several domains for protein-protein interactions including N-terminal HEAT (H untingtin, E longation factor 3, PP2A A subunit, T or) repeats , a FAT (F RAP, A TM, and T RRAP) domain [33, 34], and an FRB (F K506-R apamycin B inding)  domain. The kinase and the FATC (FAT domain at the C-terminus) [32, 36] domains are present in the C-terminal region.
The Tor inhibitor rapamycin blocks cell proliferation and is currently used as an immunosuppressive drug for organ and tissue transplant recipients and a chemotherapy agent against a variety of solid cancers [37–40]. Rapamycin binds to the prolyl isomerase FKBP12 to form a protein-drug complex that then interacts with the Tor FRB domain in a ternary complex [41, 42]. S. cerevisiae cells treated with rapamycin display phenotypes associated with nutrient depletion including G1 cell cycle arrest, cellular volume expansion, protein synthesis inhibition, glycogen accumulation, and autophagy [35, 43].
Genetic analysis in S. cerevisiae characterizing rapamycin-resistant mutants led to the identification of FKBP12 as the intracellular receptor for rapamycin and defined Tor1 and Tor2 as the targets of the FKBP12-rapamycin complex . Subsequent studies resulted in the identification and characterization of the TOR1 and TOR2 gene products [15, 44] as well as the elucidation of the Tor signaling cascade [8, 45–47]. These pioneering studies aided in the identification of the mammalian Tor ortholog and characterization of this highly conserved signaling cascade [8, 48, 49]. Remarkably, expression of the human FKBP12 ortholog in yeast fpr1 (the gene that encodes FKBP12) deletion mutants complements to restore rapamycin sensitivity, and hybrid Tor proteins consisting of the yeast N-terminal domain fused to the mammalian Tor kinase domain are also functional in yeast . Thus, the Tor pathway has been functionally and structurally conserved from yeasts to humans over the billion years of evolution separating the two species from their last common ancestor.
We now know that the fungal and metazoan kingdoms are both within the opisthokont lineage of eukaryotes and are thus more closely related to each other than the vast majority of eukaryotic organisms [51, 52]. Moreover, these two highly successful kingdoms shared a last common ancestor as recently as one billion years ago, much more recently than most eukaryotes. Recent initiatives, particularly the UNICORN project, have facilitated the sequencing of several opisthokonts related to the last common ancestor of the metazoans and fungi, providing an interesting window into the evolution of both fungi and animals [53–56]. Thus, studies on the evolutionary trajectory of the fungal kingdom in general, and of the Tor signaling cascade in particular, promise to reveal insights about how orthologous pathways function in the more complex milieu of multicellular metazoan organisms.
In this study, we have capitalized upon the wealth of available genomic information by annotating the Tor pathway in several fungal organisms in which this pathway has not been described. Our study included selected species in the major groups of the fungal kingdom with genome sequences available: the basidiomycete Pleurotus ostreatus; the Mucorales zygomycetes Mucor circinelloides, Rhizopus oryzae, and Phycomyces blakesleeanus; the chytridiomycetes Spizellomyces punctatus and Batrachochytrium dendrobatidis; the microsporidian species Encephalitozoon cuniculi, Enterocytozoon bieneusi, and Nosema ceranae; and the related non-fungal opisthokonts Capsaspora owczarzaki, Salpingoeca rosetta, and Monosiga brevicollis. Whereas the Tor pathway is conserved throughout the eukaryotes, strikingly, microsporidian species with their highly reduced and compacted genomes lack the entire Tor pathway. We have also investigated gene and genome duplication events that resulted in two Tor homologs in S. cerevisiae, S. pombe, P. ostreatus, and B. dendrobatidis, and the loss of a second Tor homolog following a whole genome duplication event in R. oryzae.
Results and Discussion
Conservation of the Tor signaling pathway molecular components
The Tor pathway is well conserved among nearly all eukaryotic species examined to date. Tor is essential for life and is the target of the potent drug rapamycin in fungi, humans, and other eukaryotic organisms [2, 12, 19, 41]. Tor originated early during the eukaryotic radiation, as it is present in the basal eukaryotes Giardia lamblia, L. major and T. brucei, and also in plants . Here we focused on the Tor pathway in the fungal kingdom and other representative unicellular opisthokont species outside of the fungi. Several genes encoding Tor complex components, upstream regulators, and downstream effectors were identified in all major groups throughout the fungal kingdom (Figure 1, Table 1). Because we used BLASTp reciprocal best hits (RBH) to identify orthologs, in some cases the absence of an ortholog in our results may represent only a failure to detect it with this method. Further, functional homologs that are not sufficiently similar enough in sequence to be identified may be present. Overall, a high degree of pathway conservation is observed in the Tor signaling cascade throughout the fungal kingdom, with the exception of microsporidia in which all Tor pathway components are absent. In addition, several pathway components are conserved in the related unicellular opisthokonts M. brevicollis, S. rosetta, and C. owczarzaki (Table 2).
In S. cerevisiae, Tor1 and Tor2 interact with several proteins to form TORC1 and TORC2. TORC1 contains Tor1 or Tor2, Kog1, Tco89, and Lst8. Kog1 functions as a substrate-recruiting subunit in mammalian TORC1 (mTORC1) [13, 59–61]. A gene encoding a putative Kog1 homolog was identified in all species included in this study, with the exception of the microsporidia and S. rosetta (Tables 1 and 2). Mutation of TCO89 results in hypersensitivity to rapamycin and affects cellular integrity in S. cerevisiae. Tco89 is only found in S. cerevisiae and could have resulted from a specialization of the pathway exclusive to Saccharomyces species (Table 1). However, while a BLASTp search did not identify Tco89 homologs in other species, S. pombe has been reported to contain a functional homolog [GenBank:NP_588232] , and this suggests that functional homologs may exist in other fungal species as well. However, no homologs were identified when using S. pombe Tco89 as a query sequence. Lst8 binds to the Tor kinase domain in S. cerevisiae to stimulate catalytic activity  and a putative Lst8 homolog is conserved in most species analyzed except microsporidia and S. rosetta (Tables 1 and 2).
Additional BLASTp searches identified a single Tor homolog in the following species: Lachancea thermotolerans CBS6340 [GenBank:XP_002552336], Pichia pastoris GS115 [GenBank:XP_002491471], Pichia stipitis CBS6054 [GenBank:XP_001385651], Debaryomyces hansenii CBS767 [GenBank:XP_002770885], Yarrowia lipolytica CLIB122 [GenBank:XP_505106], Podospora anserina strain S mat+ [GenBank:XP_001903968], Chaetomium globosum CBS148.51 [GenBank:XP_001226647], Magnaporthe oryzae 70-15 [GenBank:XP_001414541], Gibberella zeae PH-1 [GenBank:XP_388309], Aspergillus fumigatus Af293 [GenBank:XP_755360], Aspergillus flavus NRRL3357 [GenBank:XP_002377897], Neosartorya fischeri NRRL181 [GenBank:XP_001260509], Aspergillus terreus NIH2624 [GenBank:XP_001213640], Aspergillus nidulans FGSCA4 [GenBank:XP_663586], Aspergillus clavatus NRRL1 [GenBank:XP_001275326], Aspergillus oryzae RIB40 [GenBank: XP_001826216], and Aspergillus niger CBS513.88 [GenBank:XP_001397781]. These ascomycetes were not included in further studies.
In yeast and mammals, TORC2 controls spatial aspects of growth. TORC2 includes Tor2, Avo1, Avo2, Avo3, and Bit61 [9, 22]. Avo1 plays an essential role in actin cytoskeleton polarization  and is conserved throughout the organisms that were the focus of this study (Table 1) except in S. rosetta and C. owczarzaki (Table 2). However, due to high divergence of the AVO1 sequence amongst the known homologs, it is possible that the S. rosetta and C. owczarzaki homologs cannot be identified using the BLASTp algorithm. Avo2 is a nonessential substrate adaptor for TORC2  and potential homologs were identified in most species studied, with the exceptions of S. pombe, P. ostreatus, the three microsporidian species (Table 1), and M. brevicollis (Table 2). In S. cerevisiae, Avo3 (also know as Tsc11) controls cytoskeletal dynamics [9, 22], and homologs are conserved throughout the species examined (Tables 1 and 2). Bit61, a nonessential protein that associates with TORC2 , was only identified in S. cerevisiae (Tables 1 and 2).
Upstream regulators of Tor include Tsc1 and Tsc2 (Tuberous Sclerosis 1 and 2) and the GTPase Rheb (R as h omolog e nhanced in b rain). In S. pombe and other eukaryotes, Tsc1 and Tsc2 form a GTPase-activating complex that negatively regulates the action of Rheb to activate TORC1 [26, 62, 63]. Tsc1 and Tsc2 putative homologs were identified in most species in the study, with the exception of S. cerevisiae (Table 1); however, Tsc1 homologs were not identified in B. dendrobatidis, S. rosetta, or C. owczarzaki (Tables 1 and 2). In mammals and insects, four Rag GTPases (RagA-D) bind raptor (Kog1) to mediate TORC1 signaling in response to amino acids [24, 26]. RagA and RagB are orthologs of S. cerevisiae Gtr1 whereas RagC and RagD are orthologs of S. cerevisiae Gtr2 . Gtr1 and Gtr2 form a complex along with Ego1 and Ego3 known as EGOC/GSE , which is regulated by the GTP exchange factor Vam6 . Interestingly we found presumptive Gtr1 and Gtr2 orthologs in species in this study with the exception of microsporidia (Table 1), and Gtr1 was not identified in M. brevicollis (Table 2). However, Ego1 and Ego3, which anchor Gtr1 and Gtr2 to endosomal and vacuolar membranes, are unique to S. cerevisiae and were not identified in the other species (Table 1). This suggests that the specific EGO complex architecture might be restricted to Saccharomyces, though it is possible that functional homologs may be present in other species.
In fungi, several downstream effectors are targets of Tor signaling. The PP2A phosphatase Sit4 and its regulatory subunit Tap42 regulate the expression of several TORC1-controlled genes and Gcn2-regulated translation [3, 64, 65]. A Sit4 homolog was identified in R. oryzae (Table 1), though a homolog of the regulatory subunit Tap42 was not (Table 1). In S. rosetta, a Sit4 homolog was not identified, though Tap42 was (Table 2). Both of these situations could indicate rewiring of the pathway, or alternatively these findings could be a result of insufficient similarity to detect with a RBH search method. Sch9, a member of the AGC kinase family and the functional homolog of the p70 S6 kinase, is a direct substrate of S. cerevisiae TORC1 involved in ribosome biogenesis [28, 66]. Sch9 was not identified in RBH BLASTp searches with S. cerevisiae Sch9 in M. circinelloides, P. blakesleeanus, the microsporidia (Table 1), or S. rosetta (Table 2), but was identified in all other species. Putative homologs of Ypk2 kinase, a regulator of ceramide synthesis controlled by TORC2 signaling , were only identified in S. cerevisiae, S. pombe, (Table 1), M. brevicollis, S. rosetta, and C. owczarzaki (Table 2). Notably, the failure to detect many homologs in M. brevicollis, S. rosetta, and C. owczarzaki could be due to their evolutionary distance from the other organisms in the study and their gene divergence.
Evolutionary conservation of Tor kinase protein structure
The Tor kinases are phosphatidylinositol kinase-related kinases (PIKKs)  and aside from their known protein kinase function also serve as evolutionary constrained scaffolds with several protein-protein interaction domains that mediate multi-protein complex formation . Tor protein architecture was annotated using SMART [70, 71] and we found that it has been maintained throughout the species of this study (Figure 2). N-terminal HEAT repeats are present in varying numbers and serve as scaffolding structures for protein-protein interactions [32, 72, 73]. Similarly, the FAT/toxic domain and the accompanying C-terminal FATC domain also participate in protein scaffolding and are present and conserved in all species subject to analysis. Over-expression of the FAT domain has a dominant negative effect on cell growth in yeast and this effect can be suppressed by over-expression of phospholipase C . The kinase domain has similarities to both lipid and protein kinases, but it is a bona fide protein kinase. The FRB (F KBP12-r apamycin b inding) domain is a highly conserved 100 amino acid region of Tor with several residues that are known to be required for binding the FKBP12-rapamycin drug complex .
In addition, the FRB domain contains several residues necessary for both phosphatidic acid binding and rapamycin interaction [74, 75]. In mammals, phosphatidic acid binding to the FRB region may promote the assembly of mTOR complexes. It has been proposed that the potent effect of the Tor inhibitor rapamycin may be due to direct competition for overlapping phosphatidic acid binding sites within the FRB domain . Mutation of the conserved amino acid residues S1975, W2041, or F2048 in the FRB domain confers rapamycin resistance in S. cerevisiae[4, 42, 77]. Strikingly, the corresponding amino acid residues are conserved in all species examined (Figure 3). In mTOR, phosphatidic acid binding is disrupted by mutations of L2031, F2309, and Y2105 , corresponding to S. cerevisiae Tor2 L1971, F1979, and Y2045, respectively, all of which are conserved in the studied fungal species (Figure 3).
Duplication of Tor in the fungal kingdom
Approximately 100 million years ago, the S. cerevisiae ancestor underwent a whole genome duplication (WGD) event, and species that descend from this evolutionary event retained duplicated subsets of genes. S. cerevisiae retained approximately 8% of duplicated genes [78, 79], while other fungi in this lineage have maintained varying numbers of gene duplicates, such as Candida glabrata in which a smaller number (~2%) of these pairs are maintained . Interestingly, several of these species, including S. cerevisiae, S. paradoxus, S. mikatae, S. kudriavzevii, S. bayanus, S. castellii, and C. glabrata have retained two TOR paralogs, whereas species outside of the WGD clade, including Kluyveromyces lactis and Ashbya gossypii, have only a single TOR gene.
The fates of duplicated gene paralogs can be explained through the duplication, degeneration, and complementation (DDC) model of gene duplication . In species with a single Tor, this protein populates both TORC1 and TORC2, each of which has distinct functions. Following duplication, each paralog likely sub-functionalized to carry out some, but not all, of its previous functions so that between the two paralogs, each process is performed. In S. cerevisiae, for example, Tor1 exclusively functions within TORC1, while Tor2 preferentially populates TORC2 but can also function in TORC1. S. pombe has two Tor paralogs as a result of an independent segmental gene duplication event, and each paralog has distinct roles not necessarily equivalent to those of S. cerevisiae Tor1 and Tor2. SpTor1 is not essential, whereas SpTor2 is essential for growth [11, 29, 82].
We hypothesize that similar models could apply in species encoding two Tor homologs resulting from independent gene duplication events such as P. ostreatus and B. dendrobatidis. Based on their level of identity with S. cerevisiae Tor homologs, the B. dendrobatidis genes BDEG_08293 and BDEG_05727 have been designated TOR1 and TOR2, respectively. However, these gene names are not necessarily based on functional similarity with ScTor1 and ScTor2 but rather indicate only that there are two TOR homologs in B. dendrobatidis. A similar naming challenge is evident in S. pombe, where the functions of SpTor1 and SpTor2 are not equivalent with those of ScTor1 and ScTor2.
While several species, such as those in the Saccharomyces clade, have two Tor paralogs resulting from ancient WGD events, other species have acquired a second paralog through independent segmental gene duplication events, including B. dendrobatidis and several fission yeast species closely related to S. pombe (Figure 4). Segmental gene duplications have played a major role in the evolution of species and can result in the duplication of single genes or large blocks of genes . In S. cerevisiae, the duplication blocks with TOR include three paralogous gene pairs that have been retained from the ancestral WGD event: PTK1/PTK2, TOR1/TOR2, and MNN4/YJR061W  (Figure 5A). Further, C. glabrata, an ascomycete closely related to S. cerevisiae that diverged following the WGD, has also maintained a second TOR paralog and synteny is apparent in the region flanking the TOR1 and TOR2 genes, though one must extend the analysis to >40 kb from the TOR genes to detect it (See additional file 1: Supplemental figure 1). However, because C. glabrata and S. cerevisiae have retained different subsets of genes from the ancestral WGD, the duplicated genes surrounding TOR in these species are not homologs of one another. Further, C. glabrata has maintained fewer of the synteny blocks found in S. cerevisiae, and the block containing the TOR genes (Block 42) is not conserved in C. glabrata.
In contrast, genomic regions surrounding the Tor paralogs in S. pombe (Figure 5B) and B. dendrobatidis (Figure 5C) are not syntenically conserved. That is, the duplicated genes are not flanked by other duplicated gene pairs, and the two TOR paralogs likely resulted from independent segmental gene duplication events. Further, there are no indications of WGDs in these species. In the Schizosaccharomycotina, the TOR gene duplication occurred prior to speciation as evidenced by the fact that S. pombe, S. octosporus, S. japonicus, and S. cryophilus have each retained a second Tor paralog; the Tor orthologs across these four species are more closely related to each other than to those from other species (Figure 4). No syntenic conservation between the genomic regions containing TOR paralogs was detected in any species of the Schizosaccharomyces group, supporting the occurrence of an independent segmental gene duplication in their common ancestor (see additional file 2: supplemental Figure 2).
The Taphrinomycotina is a monophyletic taxon of the Ascomycota containing the Schizosaccharomycetes, Taphrinomycetes, Neolectomycetes, and Pneumocystidomycetes. Genomes are only available for species within the Schizosaccharomycetes (S. pombe, S. octosporus, S. japonicus, and S. cryophilus) and the Taphrinomycetes (Pneumocystis carinii). While the studied Schizosaccharomyces spp. contain two Tor homologs, a single homolog was identified in P. carinii with tBLASTx using the highly conserved S. cerevisiae FRB domain of Tor1 based on sequence data produced by the Pneumocystis Genome Project (funded by the NIH NIAID) that was obtained from http://pgp.cchmc.org at the time of publication. Based on the phylogenetic placement of these two groups  and the available genomic information, we hypothesize that the segmental TOR duplication is unique to the Schizosaccharomyces spp. and occurred following their divergence from the common ancestor of the Taphrinomycota, as P. carinii, the other species within this taxon with available genome data, appears to contain a single Tor homolog.
Interestingly, the zygomycete R. oryzae underwent a whole genome duplication event and despite retaining 12% of the resulting duplicated genes  compared to 8% retained genes in a separate WGD event in S. cerevisiae[78, 79], two Tor paralogs have not been retained; only a single Tor homolog was identified. Currently, the Zygomycota is the only group studied without an identified species containing two Tor homologs; however, there are not enough available genomes to make a definitive conclusion in this group. Other species with independent TOR duplications likely remain to be identified in all of the major fungal lineages, and as more genome sequences become available this will be addressed.
During this study, we identified three additional species that possibly have two TOR paralogs: M. brevicollis, Thecamonas trahens, and Allomyces macrogynus, species that are included in the Origins of Multicellularity Project . While analysis is limited because of the state of the M. brevicollis (choanoflagellate) genome, it appears that an additional bonafide Tor protein is encoded with the conserved Tor protein architecture. Thecamonas trahens is a unicellular apusomonad that also contains two TOR paralogs. Interestingly, A. macrogynus is a fungus that may contain a second paralog; however, it appears that the FRB domain overlaps with the kinase domain. This may represent a case in which TOR was duplicated but is in the process of losing one of the copies. These cases are a springboard for future studies involving the duplication of TOR within and beyond the fungal kingdom.
Loss of the Tor signaling pathway in microsporidia
While components of the Tor signaling cascade are conserved in eukaryotes from yeasts to humans, the microsporidia, a group of highly specialized obligate intracellular pathogens completely lack all of the pathway components investigated in this study. E. cuniculi has previously been noted to lack the Tor kinase , and we found that all pathway components including Tor were missing not only in E. cuniculi but also in its microsporidian relatives E. bieneusi and N. ceranae (Table 2). Microsporidian genomes are highly compacted and show a marked reduced gene content. The genome of the human pathogen E. cuniculi has 1997 genes distributed throughout 11 chromosomes in a 2.9 Mb genome , less than half the number of protein-coding genes found in Escherichia coli. For comparison, the genome of the human pathogen E. bieneusi genome is ~6 Mb with 3,804 genes  and the honeybee pathogen N. ceranae genome is ~7.9 Mb with 2,614 predicted genes . Genome reduction has been well studied in intracellular bacteria: during the specialization process from a free-living to intracellular lifestyle, massive gene loss can occur [91–95]. We hypothesize that as a consequence of their specialization to an obligate intracellular lifestyle, the need for nutrient sensing was relaxed and these obligate intracellular microsporidia acquired nutrients from within the host cytoplasm. Thus, they no longer required the Tor pathway, which is dedicated to sensing changes in the nutrient composite of the environmental milieu. However, E. cuniculi has retained some of the genes encoding signaling pathways involved in nutrient sensing in other species, including homologs of S. cerevisiae Ras1 and Ras2 (data not shown).
Interestingly, P. carinii is an obligate pathogen with a reduced genome estimated to be ~8 Mb , which unlike microsporidia is not obligately intracellular. A single Tor homolog was identified [Pneumocystis Genome Project:cap3_it1_grp346_contig490] using the S. cerevisiae Tor1 FRB domain in tBLASTx, supporting our hypothesis that the microsporidia have streamlined their genomes as they adopted an obligate intracellular life cycle. P. carinii likely requires Tor because it must survive in a metabolically active form outside of host cells, undergoing cell division in the extracellular milieu of the infected host lung.
Many microsporidian genomes have undergone extreme genome reduction and are among the eukaryotes with the smallest genomes. In addition to fewer genes than other species (1997 in E. cuniculi compared to 6607 in S. cerevisiae), E. cuniculi has shorter genes, small intergenic distances, and very few introns [88, 97]. This species (and to an even greater extent, E. bieneusi) lacks several genes that are necessary for biosynthetic pathways and the tricarboxylic acid cycle in free-living organisms [88, 89, 98]. Remarkably, E. cuniculi must acquire ATP from its host through a series of ATP transporters anciently acquired from Rickettsia- or Chlamydia-like bacteria , and this is an example of how E. cuniculi, and possibly other microsporidia, can hijack mechanisms of other intracellular bacteria, fungi, or the host to survive in and adapt to their highly specialized intracellular lifestyle.
The Tor pathway is highly conserved with some exceptions among fungi, including several pathogens. Similarly, the molecular structural organization of the Tor kinases has been remarkably conserved as well as the presumed ability to bind the FKBP12-rapamycin complex. This feature, combined with the essential nature of the rapamycin-sensitive TORC1 pathway for cell growth, should enable the development of rapamycin-based strategies for antifungal therapies. Duplications of the Tor protein occurred in most fungal groups examined, resulting from either independent segmental gene duplication events or a WGD event. The maintenance of two Tor homologs can be explained through the DDC model of gene duplication, in which the paralogs subfunctionalize and together support all of the functions of the pre-duplicated protein. In addition, species outside of the fungal kingdom with two or more TOR homologs include the metazoan B. mori, the apusomonad protozoan T. trahens, the trypanosome protozoans L. major and T. brucei, and possibly the choanoflagellate M. brevicollis. This suggests that the Tor gene has been independently duplicated multiple times since the last common ancestor of the Opisthokonta and Excavata lineages. The Tor pathway, as well as the Tor protein itself, is highly conserved in eukaryotes, so the observation that three microsporidian species are missing the entire pathway is striking. It would appear that the specialized obligate intracellular life cycle of microsporidia allows for the acquisition of nutrients from the host cell, obviating a requirement for this central nutrient-sensing pathway that is essential in all known cases of other eukaryotes for axenic growth.
Identification of pathway component homologs
Tor homologs were identified using the highly conserved 100 amino acid FRB domain of S. cerevisiae Tor1 in BLASTp reciprocal best-hit searches between the specific species database and the Saccharomyces Genome Database (SGD). All queries for other Tor pathway components were retrieved from SGD and homologs were identified by BLASTp reciprocal best-hit searches, with the exception of Tsc1 and Tsc2. Tsc1 and Tsc2 protein query sequences were retrieved from the S. pombe GeneDB, and homologs were identified using BLASTp searches between the S. pombe genome database and the species database. Species identification numbers are listed in Supplemental Table 1. Protein architecture of each homolog was elucidated with Simple Modular Architecture Research Tool (SMART) analysis [70, 71].
Accession numbers are provided for putative Tor pathway components (see Additional file 3, supplemental Table 1 for the accession numbers of putative Tor pathway components) from the following public genomic databases that were used to carry out BLAST: S. cerevisiae, SGD, http://www.yeastgenome.org; S. pombe, GeneDB, http://www.genedb.org/genedb/pombe/index.jsp; P. ostreatus, JGI, http://genome.jgi-psf.org/PleosPC15_1/PleosPC15_1.home.html; M. circinelloides, JGI and the Mucor Genome Project, http://genome.jgi-psf.org/Mucci1/Mucci1.home.html; R. oryzae, Broad Institute, http://www.broadinstitute.org/annotation/genome/rhizopus_oryzae/MultiHome.html; P. blakesleeanus, JGI, http://genome.jgi-psf.org/Phybl1/Phybl1.home.html; B. dendrobatidis, Broad Institute, http://www.broadinstitute.org/annotation/genome/batrachochytrium_dendrobatidis/MultiHome.html; S. punctatus, Broad Institute and the UNICORN initiative, http://www.broadinstitute.org/annotation/genome/multicellularity_project/MultiHome.html; M. brevicollis, JGI, http://genome.jgi-psf.org/Monbr1/Monbr1.home.html; S. rosetta, Broad Institute and the UNICORN initiative, http://www.broadinstitute.org/annotation/genome/multicellularity_project/MultiHome.html; C. owczarzaki, Broad Institute and the UNICORN initiative, http://www.broadinstitute.org/annotation/genome/multicellularity_project/MultiHome.html; E. cuniculi, E. bieneusi, and N. ceranae, NCBI, http://www.ncbi.nlm.nih.gov/sutils/genom_table.cgi?organism=fungi.
Phylogenetic analysis of Tor
The ProtTest program  was used to select an appropriate model to construct the phylogenetic relationship of the Tor proteins of the studied species. For this study, the best model for phylogenetic construction of the Tor protein using maximum likelihood (ML) was JTT+I+G+F . The ML tree was generated using PhyML with the JTT+I+G+F mode with 500 bootstrap replicates [102, 103]. The Tor amino acid sequences from all species were aligned with CLUSTAL W .
Synteny analysis of TOR paralogs
Tor homologs in fungal species were identified with BLAST  using TOR1 and TOR2 of S. cerevisiae as query sequences. The TOR sequence and flanking regions (up to 100 kb) were extracted from corresponding contigs or chromosomal sequences. To analyze the synteny of TOR and flanking regions, we constructed a database of one TOR and flanking region sequence and used the other TOR homolog and flanking region sequence to perform tBLASTx searches in the database (BLOSUM62 matrix; E-value = 1e-3; gap cost: existence 11, extension 1; and low complexity regions filter). The tBLASTx results were parsed and formatted using an in-house Perl script (available upon request). The parsed and formatted tBLASTx results were then imported into Artemis software  for synteny analysis. We employed a cutoff value of 100 bp to filter short non-significant matches between query and database sequences.
De Virgilio C, Loewith R: The TOR signalling network from yeast to man. Int J Biochem Cell Biol. 2006, 38: 1476-1481. 10.1016/j.biocel.2006.02.013.
Bastidas RJ, Reedy JL, Morales-Johansson H, Heitman J, Cardenas ME: Signaling cascades as drug targets in model and pathogenic fungi. Curr Opin Investig Drugs. 2008, 9: 856-864.
Rohde JR, Bastidas R, Puria R, Cardenas ME: Nutritional control via Tor signaling in Saccharomyces cerevisiae. Curr Opin Microbiol. 2008, 11: 153-160. 10.1016/j.mib.2008.02.013.
Heitman J, Movva NR, Hall MN: Targets for cell cycle arrest by the immunosuppressant rapamycin in yeast. Science. 1991, 253: 905-909. 10.1126/science.1715094.
Otsubo Y, Yamamato M: TOR signaling in fission yeast. Crit Rev Biochem Mol Biol. 2008, 43: 277-283. 10.1080/10409230802254911.
Petersen J: TOR signalling regulates mitotic commitment through stress-activated MAPK and Polo kinase in response to nutrient stress. Biochem Soc Trans. 2009, 37: 273-277. 10.1042/BST0370273.
Bastidas RJ, Heitman J, Cardenas ME: The protein kinase Tor1 regulates adhesin gene expression in Candida albicans. PLoS Pathog. 2009, 5: e1000294-10.1371/journal.ppat.1000294.
Rohde JR, Cardenas ME: Nutrient signaling through TOR kinases controls gene expression and cellular differentiation in fungi. Curr Top Microbiol Immunol. 2004, 279: 53-72.
Loewith R, Jacinto E, Wullschleger S, Lorberg A, Crespo JL, Bonenfant D, Oppliger W, Jenoe P, Hall MN: Two TOR complexes, only one of which is rapamycin sensitive, have distinct roles in cell growth control. Mol Cell. 2002, 10: 457-468. 10.1016/S1097-2765(02)00636-6.
Wedaman KP, Reinke A, Anderson S, Yates J, McCaffery JM, Powers T: Tor kinases are in distinct membrane-associated protein complexes in Saccharomyces cerevisiae. Mol Biol Cell. 2003, 14: 1204-1220. 10.1091/mbc.E02-09-0609.
Matsuo T, Otsubo Y, Urano J, Tamanoi F, Yamamoto M: Loss of the TOR kinase Tor2 mimics nitrogen starvation and activates the sexual development pathway in fission yeast. Mol Cell Biol. 2007, 27: 3154-3164. 10.1128/MCB.01039-06.
Hayashi T, Hatanaka M, Nagao K, Nakaseko Y, Kanoh J, Kokubu A, Ebe M, Yanagida M: Rapamycin sensitivity of the Schizosaccharomyces pombe tor2 mutant and organization of two highly phosphorylated TOR complexes by specific and common subunits. Genes Cells. 2007, 12: 1357-1370. 10.1111/j.1365-2443.2007.01141.x.
Kim DH, Sarbassov DD, Ali SM, King JE, Latek RR, Erdjument-Bromage H, Tempst P, Sabatini DM: mTOR interacts with raptor to form a nutrient-sensitive complex that signals to the cell growth machinery. Cell. 2002, 110: 163-175. 10.1016/S0092-8674(02)00808-5.
Sarbassov DD, Ali SM, Kim DH, Guertin DA, Latek RR, Erdjument-Bromage H, Tempst P, Sabatini DM: Rictor, a novel binding partner of mTOR, defines a rapamycin-insensitive and raptor-independent pathway that regulates the cytoskeleton. Curr Biol. 2004, 14: 1296-1302. 10.1016/j.cub.2004.06.054.
Helliwell SB, Wagner P, Kunz J, Deuter-Reinhard M, Henriquez R, Hall MN: TOR1 and TOR2 are structurally and functionally similar but not identical phosphatidylinositol kinase homologues in yeast. Mol Biol Cell. 1994, 5: 105-118.
Zheng XF, Florentino D, Chen J, Crabtree GR, Schreiber SL: TOR kinase domains are required for two distinct functions, only one of which is inhibited by rapamycin. Cell. 1995, 82: 121-130. 10.1016/0092-8674(95)90058-6.
Zhou S, Zhou Q, Liu Y, Wang S, Wen D, He Q, Wang W, Bendena WG, Li S: Two Tor genes in the silkworm Bombyx mori. Insect Mol Biol. 2010
Madeira da Silva L, Beverley SM: Expansion of the target of rapamycin (TOR) kinase family and function in Leishmania shows that TOR3 is required for acidocalcisome biogenesis and animal infectivity. Proc Natl Acad Sci USA. 2010, 107: 11965-11970. 10.1073/pnas.1004599107.
Barquilla A, Crespo JL, Navarro M: Rapamycin inhibits trypanosome cell growth by preventing TOR complex 2 formation. Proc Natl Acad Sci USA. 2008, 105: 14579-14584. 10.1073/pnas.0802668105.
Adami A, Garcia-Alvarez B, Arias-Palomo E, Barford D, Llorca O: Structure of TOR and its complex with KOG1. Mol Cell. 2007, 27: 509-516. 10.1016/j.molcel.2007.05.040.
Reinke A, Anderson S, McCaffery JM, Yates J, Aronova S, Chu S, Fairclough S, Iverson C, Wedaman KP, Powers T: TOR complex 1 includes a novel component, Tco89p (YPL180w), and cooperates with Ssd1p to maintain cellular integrity in Saccharomyces cerevisiae. J Biol Chem. 2004, 279: 14752-14762. 10.1074/jbc.M313062200.
Wullschleger S, Loewith R, Oppliger W, Hall MN: Molecular organization of target of rapamycin complex 2. J Biol Chem. 2005, 280: 30697-30704. 10.1074/jbc.M505553200.
Binda M, Peli-Gulli MP, Bonfils G, Panchaud N, Urban J, Sturgill TW, Loewith R, De Virgilio C: The Vam6 GEF controls TORC1 by activating the EGO complex. Mol Cell. 2009, 35: 563-573. 10.1016/j.molcel.2009.06.033.
Kim E, Goraksha-Hicks P, Li L, Neufeld TP, Guan KL: Regulation of TORC1 by Rag GTPases in nutrient response. Nat Cell Biol. 2008, 10: 935-945. 10.1038/ncb1753.
Gao M, Kaiser CA: A conserved GTPase-containing complex is required for intracellular sorting of the general amino-acid permease in yeast. Nat Cell Biol. 2006, 8: 657-667. 10.1038/ncb1419.
Sancak Y, Peterson TR, Shaul YD, Lindquist RA, Thoreen CC, Bar-Peled L, Sabatini DM: The Rag GTPases bind raptor and mediate amino acid signaling to mTORC1. Science. 2008, 320: 1496-1501. 10.1126/science.1157535.
Rohde JR, Campbell S, Zurita-Martinez SA, Cutler NS, Ashe M, Cardenas ME: TOR controls transcriptional and translational programs via Sap-Sit4 protein phosphatase signaling effectors. Mol Cell Biol. 2004, 24: 8332-8341. 10.1128/MCB.24.19.8332-8341.2004.
Urban J, Soulard A, Huber A, Lippman S, Mukhopadhyay D, Deloche O, Wanke V, Anrather D, Ammerer G, Riezman H: Sch9 is a major target of TORC1 in Saccharomyces cerevisiae. Mol Cell. 2007, 26: 663-674. 10.1016/j.molcel.2007.04.020.
Weisman R, Choder M: The fission yeast TOR homolog, tor1+, is required for the response to starvation and other stresses via a conserved serine. J Biol Chem. 2001, 276: 7027-7032. 10.1074/jbc.M010446200.
Alvarez B, Moreno S: Fission yeast Tor2 promotes cell growth and represses cell differentiation. J Cell Sci. 2006, 119: 4475-4485. 10.1242/jcs.03241.
Petersen J, Nurse P: TOR signalling regulates mitotic commitment through the stress MAP kinase pathway and the Polo and Cdc2 kinases. Nat Cell Biol. 2007, 9: 1263-1272. 10.1038/ncb1646.
Kobe B, Gleichmann T, Horne J, Jennings IG, Scotney PD, Teh T: Turn up the HEAT. Structure. 1999, 7: R91-97. 10.1016/S0969-2126(99)80060-4.
Alarcon CM, Heitman J, Cardenas ME: Protein kinase activity and identification of a toxic effector domain of the target of rapamycin TOR proteins in yeast. Mol Biol Cell. 1999, 10: 2531-2546.
Bosotti R, Isacchi A, Sonnhammer EL: FAT: a novel domain in PIK-related kinases. Trends Biochem Sci. 2000, 25: 225-227. 10.1016/S0968-0004(00)01563-2.
Barbet NC, Schneider U, Helliwell SB, Stansfield I, Tuite MF, Hall MN: TOR controls translation initiation and early G1 progression in yeast. Mol Biol Cell. 1996, 7: 25-42.
Dames SA, Mulet JM, Rathgeb-Szabo K, Hall MN, Grzesiek S: The solution structure of the FATC domain of the protein kinase target of rapamycin suggests a role for redox-dependent structural and cellular stability. J Biol Chem. 2005, 280: 20558-20564. 10.1074/jbc.M501116200.
Miller RE, Larkin JM: Combination Systemic Therapy for Advanced Renal Cell Carcinoma. Oncologist. 2009, 14: 1218-1224. 10.1634/theoncologist.2009-0105.
Tsang CK, Qi H, Liu LF, Zheng XF: Targeting mammalian target of rapamycin (mTOR) for health and diseases. Drug Discov Today. 2007, 12: 112-124. 10.1016/j.drudis.2006.12.008.
Kruck S, Bedke J, Hennenlotter J, Ohneseit PA, Kuehs U, Senger E, Sievert KD, Stenzl A: Activation of mTOR in renal cell carcinoma is due to increased phosphorylation rather than protein overexpression. Oncol Rep. 2010, 23: 159-163.
Lieberthal W, Levine JS: The role of the mammalian target of rapamycin (mTOR) in renal disease. J Am Soc Nephrol. 2009, 20: 2493-2502. 10.1681/ASN.2008111186.
Choi J, Chen J, Schreiber SL, Clardy J: Structure of the FKBP12-rapamycin complex interacting with the binding domain of human FRAP. Science. 1996, 273: 239-242. 10.1126/science.273.5272.239.
Lorenz MC, Heitman J: TOR mutations confer rapamycin resistance by preventing interaction with FKBP12-rapamycin. J Biol Chem. 1995, 270: 27531-27537. 10.1074/jbc.270.10.5483.
Kamada Y, Funakoshi T, Shintani T, Nagano K, Ohsumi M, Ohsumi Y: Tor-mediated induction of autophagy via an Apg1 protein kinase complex. J Cell Biol. 2000, 150: 1507-1513. 10.1083/jcb.150.6.1507.
Kunz J, Henriquez R, Schneider U, Deuter-Reinhard M, Movva NR, Hall MN: Target of rapamycin in yeast, TOR2, is an essential phosphatidylinositol kinase homolog required for G1 progression. Cell. 1993, 73: 585-596. 10.1016/0092-8674(93)90144-F.
Cutler NS, Heitman J, Cardenas ME: TOR kinase homologs function in a signal transduction pathway that is conserved from yeast to mammals. Mol Cell Endocrinol. 1999, 155: 135-142. 10.1016/S0303-7207(99)00121-5.
Duvel K, Broach JR: The role of phosphatases in TOR signaling in yeast. Curr Top Microbiol Immunol. 2004, 279: 19-38.
Dann SG, Thomas G: The amino acid sensitive TOR pathway from yeast to mammals. FEBS Lett. 2006, 580: 2821-2829. 10.1016/j.febslet.2006.04.068.
Sabatini DM, Erdjument-Bromage H, Lui M, Tempst P, Snyder SH: RAFT1: a mammalian protein that binds to FKBP12 in a rapamycin-dependent fashion and is homologous to yeast TORs. Cell. 1994, 78: 35-43. 10.1016/0092-8674(94)90570-3.
Brown EJ, Albers MW, Shin TB, Ichikawa K, Keith CT, Lane WS, Schreiber SL: A mammalian protein targeted by G1-arresting rapamycin-receptor complex. Nature. 1994, 369: 756-758. 10.1038/369756a0.
Alarcon CM, Cardenas ME, Heitman J: Mammalian RAFT1 kinase domain provides rapamycin-sensitive TOR function in yeast. Genes Dev. 1996, 10: 279-288. 10.1101/gad.10.3.279.
Baldauf SL, Palmer JD: Animals and fungi are each other's closest relatives: congruent evidence from multiple proteins. Proc Natl Acad Sci USA. 1993, 90: 11558-11562. 10.1073/pnas.90.24.11558.
Wainright PO, Hinkle G, Sogin ML, Stickel SK: Monophyletic origins of the metazoa: an evolutionary link with fungi. Science. 1993, 260: 340-342. 10.1126/science.8469985.
King N, Westbrook MJ, Young SL, Kuo A, Abedin M, Chapman J, Fairclough S, Hellsten U, Isogai Y, Letunic I: The genome of the choanoflagellate Monosiga brevicollis and the origin of metazoans. Nature. 2008, 451: 783-788. 10.1038/nature06617.
King N: The unicellular ancestry of animal development. Dev Cell. 2004, 7: 313-325. 10.1016/j.devcel.2004.08.010.
Ruiz-Trillo I, Burger G, Holland PW, King N, Lang BF, Roger AJ, Gray MW: The origins of multicellularity: a multi-taxon genome initiative. Trends Genet. 2007, 23: 113-118. 10.1016/j.tig.2007.01.005.
Ruiz-Trillo I, Roger AJ, Burger G, Gray MW, Lang BF: A phylogenomic investigation into the origin of metazoa. Mol Biol Evol. 2008, 25: 664-672. 10.1093/molbev/msn006.
Morrison HG, Zamora G, Campbell RK, Sogin ML: Inferring protein function from genomic sequence: Giardia lamblia expresses a phosphatidylinositol kinase-related kinase similar to yeast and mammalian TOR. Comp Biochem Physiol B Biochem Mol Biol. 2002, 133: 477-491. 10.1016/S1096-4959(02)00218-X.
Menand B, Meyer C, Robaglia C: Plant growth and the TOR pathway. Curr Top Microbiol Immunol. 2004, 279: 97-113.
Nojima H, Tokunaga C, Eguchi S, Oshiro N, Hidayat S, Yoshino K, Hara K, Tanaka N, Avruch J, Yonezawa K: The mammalian target of rapamycin (mTOR) partner, raptor, binds the mTOR substrates p70 S6 kinase and 4E-BP1 through their TOR signaling (TOS) motif. J Biol Chem. 2003, 278: 15461-15464. 10.1074/jbc.C200665200.
Yonezawa K, Tokunaga C, Oshiro N, Yoshino K: Raptor, a binding partner of target of rapamycin. Biochem Biophys Res Commun. 2004, 313: 437-441. 10.1016/j.bbrc.2003.07.018.
Schalm SS, Fingar DC, Sabatini DM, Blenis J: TOS motif-mediated raptor binding regulates 4E-BP1 multisite phosphorylation and function. Curr Biol. 2003, 13: 797-806. 10.1016/S0960-9822(03)00329-4.
Urano J, Comiso MJ, Guo L, Aspuria PJ, Deniskin R, Tabancay AP, Kato-Stankiewicz J, Tamanoi F: Identification of novel single amino acid changes that result in hyperactivation of the unique GTPase, Rheb, in fission yeast. Mol Microbiol. 2005, 58: 1074-1086. 10.1111/j.1365-2958.2005.04877.x.
Aspuria PJ, Sato T, Tamanoi F: The TSC/Rheb/TOR signaling pathway in fission yeast and mammalian cells: temperature sensitive and constitutive active mutants of TOR. Cell Cycle. 2007, 6: 1692-1695. 10.4161/cc.6.14.4478.
Jiang Y, Broach JR: Tor proteins and protein phosphatase 2A reciprocally regulate Tap42 in controlling cell growth in yeast. EMBO J. 1999, 18: 2782-2792. 10.1093/emboj/18.10.2782.
Cherkasova VA, Hinnebusch AG: Translational control by TOR and TAP42 through dephosphorylation of eIF2alpha kinase GCN2. Genes Dev. 2003, 17: 859-872. 10.1101/gad.1069003.
Jorgensen P, Nishikawa JL, Breitkreutz BJ, Tyers M: Systematic identification of pathways that couple cell growth and division in yeast. Science. 2002, 297: 395-400. 10.1126/science.1070850.
Aronova S, Wedaman K, Aronov PA, Fontes K, Ramos K, Hammock BD, Powers T: Regulation of ceramide biosynthesis by TOR complex 2. Cell Metab. 2008, 7: 148-158. 10.1016/j.cmet.2007.11.015.
Keith CT, Schreiber SL: PIK-related kinases: DNA repair, recombination, and cell cycle checkpoints. Science. 1995, 270: 50-51. 10.1126/science.270.5233.50.
Alvarez-Ponce D, Aguade M, Rozas J: Network-level molecular evolutionary analysis of the insulin/TOR signal transduction pathway across 12 Drosophila genomes. Genome Res. 2009, 19: 234-242. 10.1101/gr.084038.108.
Schultz J, Milpetz F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domains. Proc Natl Acad Sci USA. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.
Letunic I, Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic Acids Res. 2009, 37: D229-232. 10.1093/nar/gkn808.
Andrade MA, Bork P: HEAT repeats in the Huntington's disease protein. Nat Genet. 1995, 11: 115-116. 10.1038/ng1095-115.
Andrade MA, Petosa C, O'Donoghue SI, Muller CW, Bork P: Comparison of ARM and HEAT protein repeats. J Mol Biol. 2001, 309: 1-18. 10.1006/jmbi.2001.4624.
Liang J, Choi J, Clardy J: Refined structure of the FKBP12-rapamycin-FRB ternary complex at 2.2 A resolution. Acta Crystallogr D Biol Crystallogr. 1999, 55: 736-744. 10.1107/S0907444998014747.
Veverka V, Crabbe T, Bird I, Lennie G, Muskett FW, Taylor RJ, Carr MD: Structural characterization of the interaction of mTOR with phosphatidic acid and a novel class of inhibitor: compelling evidence for a central role of the FRB domain in small molecule-mediated regulation of mTOR. Oncogene. 2008, 27: 585-595. 10.1038/sj.onc.1210693.
Toschi A, Lee E, Xu L, Garcia A, Gadir N, Foster DA: Regulation of mTORC1 and mTORC2 complex assembly by phosphatidic acid: competition with rapamycin. Mol Cell Biol. 2009, 29: 1411-1420. 10.1128/MCB.00782-08.
Chen J, Zheng XF, Brown EJ, Schreiber SL: Identification of an 11-kDa FKBP12-rapamycin-binding domain within the 289-kDa FKBP12-rapamycin-associated protein and characterization of a critical serine residue. Proc Natl Acad Sci USA. 1995, 92: 4947-4951. 10.1073/pnas.92.11.4947.
Wolfe KH, Shields DC: Molecular evidence for an ancient duplication of the entire yeast genome. Nature. 1997, 387: 708-713. 10.1038/42711.
Seoighe C, Wolfe KH: Extent of genomic rearrangement after genome duplication in yeast. Proc Natl Acad Sci USA. 1998, 95: 4447-4452. 10.1073/pnas.95.8.4447.
Dujon B, Sherman D, Fischer G, Durrens P, Casaregola S, Lafontaine I, De Montigny J, Marck C, Neuveglise C, Talla E: Genome evolution in yeasts. Nature. 2004, 430: 35-44. 10.1038/nature02579.
van Hoof A: Conserved functions of yeast genes support the duplication, degeneration and complementation model for gene duplication. Genetics. 2005, 171: 1455-1461. 10.1534/genetics.105.044057.
Weisman R, Roitburg I, Schonbrun M, Harari R, Kupiec M: Opposite effects of tor1 and tor2 on nitrogen starvation responses in fission yeast. Genetics. 2007, 175: 1153-1162. 10.1534/genetics.106.064170.
Koszul R, Caburet S, Dujon B, Fischer G: Eucaryotic genome evolution through the spontaneous duplication of large chromosomal segments. EMBO J. 2004, 23: 234-243. 10.1038/sj.emboj.7600024.
Liu Y, Leigh JW, Brinkmann H, Cushion MT, Rodriguez-Ezpeleta N, Philippe H, Lang BF: Phylogenomic analyses support the monophyly of Taphrinomycotina, including Schizosaccharomyces fission yeasts. Mol Biol Evol. 2009, 26: 27-34. 10.1093/molbev/msn221.
Ma LJ, Ibrahim AS, Skory C, Grabherr MG, Burger G, Butler M, Elias M, Idnurm A, Lang BF, Sone T: Genomic analysis of the basal lineage fungus Rhizopus oryzae reveals a whole-genome duplication. PLoS Genet. 2009, 5: e1000549-10.1371/journal.pgen.1000549.
Ruiz-Trillo I, Burger G, Holland PWH, King N, Lang BF, Roger AJ, Gray MW: The origins of multicellularity: a multi-taxon genome initiative. Trends Genet. 2007, 23: 113-118. 10.1016/j.tig.2007.01.005.
Miranda-Saavedra D, Stark MJ, Packer JC, Vivares CP, Doerig C, Barton GJ: The complement of protein kinases of the microsporidium Encephalitozoon cuniculi in relation to those of Saccharomyces cerevisiae and Schizosaccharomyces pombe. BMC Genomics. 2007, 8: 309-10.1186/1471-2164-8-309.
Katinka MD, Duprat S, Cornillot E, Metenier G, Thomarat F, Prensier G, Barbe V, Peyretaillade E, Brottier P, Wincker P: Genome sequence and gene compaction of the eukaryote parasite Encephalitozoon cuniculi. Nature. 2001, 414: 450-453. 10.1038/35106579.
Akiyoshi DE, Morrison HG, Lei S, Feng X, Zhang Q, Corradi N, Mayanja H, Tumwine JK, Keeling PJ, Weiss LM, Tzipori S: Genomic survey of the non-cultivatable opportunistic human pathogen, Enterocytozoon bieneusi. PLoS Pathog. 2009, 5: e1000261-10.1371/journal.ppat.1000261.
Cornman RS, Chen YP, Schatz MC, Street C, Zhao Y, Desany B, Egholm M, Hutchison S, Pettis JS, Lipkin WI, Evans JD: Genomic analyses of the microsporidian Nosema ceranae, an emergent pathogen of honey bees. PLoS Pathog. 2009, 5: e1000466-10.1371/journal.ppat.1000466.
Gomez-Valero L, Rocha EP, Latorre A, Silva FJ: Reconstructing the ancestor of Mycobacterium leprae: the dynamics of gene loss and genome reduction. Genome Res. 2007, 17: 1178-1185. 10.1101/gr.6360207.
Boussau B, Karlberg EO, Frank AC, Legault BA, Andersson SG: Computational inference of scenarios for alpha-proteobacterial genome evolution. Proc Natl Acad Sci USA. 2004, 101: 9722-9727. 10.1073/pnas.0400975101.
Darby AC, Cho NH, Fuxelius HH, Westberg J, Andersson SG: Intracellular pathogens go extreme: genome evolution in the Rickettsiales. Trends Genet. 2007, 23: 511-520. 10.1016/j.tig.2007.08.002.
Merhej V, Royer-Carenzi M, Pontarotti P, Raoult D: Massive comparative genomic analysis reveals convergent evolution of specialized bacteria. Biol Direct. 2009, 4: 13-10.1186/1745-6150-4-13.
Yus E, Maier T, Michalodimitrakis K, van Noort V, Yamada T, Chen WH, Wodke JA, Guell M, Martinez S, Bourgeois R: Impact of genome reduction on bacterial metabolism and its regulation. Science. 2009, 326: 1263-1268. 10.1126/science.1177263.
Cushion MT, Kaselis M, Stringer SL, Stringer JR: Genetic stability and diversity of Pneumocystis carinii infecting rat colonies. Infect Immun. 1993, 61: 4801-4813.
Cushion MT: Comparative genomics of Pneumocystis carinii with other protists: implications for life style. J Eukaryot Microbiol. 2004, 51: 30-37. 10.1111/j.1550-7408.2004.tb00158.x.
Keeling PJ, Corradi N, Morrison HG, Haag KL, Ebert D, Weiss LM, Akiyoshi DE, Tzipori S: The Reduced Genome of the Parasitic Microsporidian Enterocytozoon bieneusi Lacks Genes for Core Carbon Metabolism. Genome Biol Evol. 2010, 2: 304-309. 10.1093/gbe/evq022.
Tsaousis AD, Kunji ER, Goldberg AV, Lucocq JM, Hirt RP, Embley TM: A novel route for ATP acquisition by the remnant mitochondria of Encephalitozoon cuniculi. Nature. 2008, 453: 553-556. 10.1038/nature06903.
Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.
Jones DT, Taylor WR, Thornton JM: The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992, 8: 275-282.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.
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: W557-559. 10.1093/nar/gki352.
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: 4673-4680. 10.1093/nar/22.22.4673.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Carver T, Berriman M, Tivey A, Patel C, Bohme U, Barrell BG, Parkhill J, Rajandream MA: Artemis and ACT: viewing, annotating and comparing sequences stored in a relational database. Bioinformatics. 2008, 24: 2672-2676. 10.1093/bioinformatics/btn529.
We thank Soo Chan Lee, Yonathan Lewit, Anna Averette, and Sheri Frank of the Heitman lab for advice and technical support, Nicolas Corradi and Patrick Keeling of the University of British Columbia for advice and genome access for the microsporidia, and Jason Stajich of the University of California Riverside and Christina Cuomo of the Broad Institute for advice. We are grateful to Antonio Pisabarro, Lucía Ramirez, and others involved in the P. ostreatus genome project, the Mucor Genome Project and JGI for the Mucor genome, the Pneumocystis Genome Project, the Origins of Multicellularity Project, and the Broad Institute and JGI for making publicly accessible several of the genomes used in this study. This work was supported by the National Institute of Allergy and Infectious Disease of the National Institutes of Health (grant AI050438-08 to M.C. and J.H.) and the National Cancer Institute of the National Institutes of Health (grant CA114107 to M.C.).
CAS carried out the homolog searches and phylogenetic analysis and drafted the manuscript. RJB performed homolog searches and participated in the design of the study. WL performed the synteny and phylogenetic analysis and helped draft the manuscript. JH and MEC conceived of the study, participated in its design and helped to draft the manuscript. All authors read and approved the final manuscript.