In Silico Biology 4, 0050 (2004); ©2004, Bioinformation Systems e.V.  


In silico reconstruction of nutrient-sensing signal transduction pathways in Aspergillus nidulans

Vignesh Muthuvijayan and Mark R. Marten*




University of Maryland, Baltimore County (UMBC)
Department of Chemical and Biochemical Engineering
1000 Hilltop Circle
Baltimore, MD 21250
USA



*  Corresponding author; Phone: +1-410-455 3439; Fax: +1-410-455 1049; Email: marten@umbc.edu





Edited by E. Wingender; received August 05, 2004; revised October 29, 2004; accepted October 31, 2004; published November 12, 2004



Abstract

We report here probable nutrient-sensing signal transduction pathways in Aspergillus nidulans, a model filamentous fungus, based on sequence homology studies with known Saccharomyces cerevisiae and Schizosaccharomyces pombe proteins. Specifically, we identified A. nidulans homologs for yeast proteins involved in (1) filamentation-invasion, (2) cAMP-PKA, (3) pheromone response, (4) cell integrity and (5) TOR signaling pathways. We have also studied autophagy, one of the most important cellular responses regulated by TOR signaling. The Basic Local Alignment Search Tool program "blastp" was used to assess the homology of proteins. We note that by using a highly conservative approach, 70% of the S. cerevisiae signal transduction proteins (107 proteins out of 153 proteins studied) have significant homologs in A. nidulans. Using a slightly less conservative approach, we are able to identify homologs for as high as 91% of the S. cerevisiae signal transduction proteins (139 proteins out of 153 proteins studied). The filamentation-invasion, cell integrity and TOR signaling pathways showed greatest similarity with S. cerevisiae, while the cAMP-PKA and pheromone response pathways showed greater similarity with S. pombe. Based on these results, probable pathways in A. nidulans were constructed using well-established S. cerevisiae and S. pombe models.

Key words: BLAST, signal transduction, bioinformatics, nutrient sensing, Aspergillus nidulans, comparative genomics



Introduction

Filamentous fungi are a diverse group of heterotrophic microorganisms that are agriculturally and medically important and are also widely used in the production of food and beverages [1], organic acids [2], enzymes [1, 3], antibiotics [4] and other pharmaceuticals [5]; as agents of biological control of pest insects and weeds, and in biomass conversion [6]. Hence filamentous fungi are currently responsible for about half the world’s pharmaceutical and biotechnology market [3]. Also, filamentous fungi have a significant medical importance as they are known to be serious human pathogens [7], especially to immuno-compromised patients [8] and have been reported to account for up to 40% of the deaths from hospital-acquired infections [9].

Aspergillus is a ubiquitous filamentous fungus found in nature [10]. It is commonly isolated from soil, plant debris, and indoor air environment [10]. The genus Aspergillus includes over 185 species, out of which approximately 20 have been identified as opportunistic pathogens [10]. The fungus also causes allergic diseases in asthmatics and patients suffering from cystic fibrosis [11]. Hence study of Aspergillus is of critical significance in the medical field. Aspergillus nidulans is a commonly used model filamentous fungus [1, 12]. It is closely related to a large number of other Aspergillus species of industrial and medical significance - e. g., A. niger, A. oryzae, A. flavus, and A. fumigatus [12]. Hence a better understanding of A. nidulans physiology will yield important information on critical mechanisms in other industrially and medically significant filamentous fungi.

Signal transduction is the cascade of processes by which an extra-cellular signal interacts with a receptor at the cell surface, causing a change in the level of a second messenger and ultimately effects a change in the cell’s function [13]. Recent studies show signal transduction pathways regulate fungal virulence in both human (Candida albicans and Cryptococcus neoformans [14]) and plant pathogens (Ustilago maydis, Magnaporthe grisea and Cryphonectria parasitica [14]). Also, many of these pathways have been identified as potential targets for new antifungal drugs [14]. Hence understanding fungal signal transduction pathways may aid in identifying targets for novel antifungal compounds.

The signal transduction pathways of model yeasts, Saccharomyces cerevisiae and Schizosaccharomyces pombe, have been studied extensively [15]. In contrast, information on signal transduction pathways in filamentous fungi is fragmentary [15]. For example, in S. cerevisiae about 2% of the proteins have been identified as kinases, classified under more than 30 families [16]. Whereas, in A. nidulans, a model filamentous fungus, only 18 kinases have been annotated. As the A. nidulans genome is approximately twice the size of S. cerevisiae genome, we predict the filamentous fungus contains at least as many, probably more, protein kinases as S. cerevisiae. This implies that many enzymes involved in the signal transduction pathways of the filamentous fungi have not yet been identified.

Earlier, Han and Prade [17] reported the reconstruction of the A. nidulans salt stress-controlling pathway based on homology analysis with known S. cerevisiae genes. To accomplish this they identified sequence homologs of the S. cerevisiae High Osmolarity Glycerol (HOG) pathway genes in the A. nidulans EST database [17]. The results they obtained by computational techniques were validated by gene deletion studies [17]. But this study was limited by the EST information available. Now, with the availability of the entire A. nidulans genome [12], sequence homology studies can be conducted more effectively. A similar approach has been conducted with Aspergillus fumigatus to identify the existence of a sexual cycle [18] and sequence homologs of the S. cerevisiae pheromone response pathway were identified from an incomplete genome sequence of A. fumigatus [18].

Here we have studied various nutrient sensing signal transduction pathways in A. nidulans, using sequence homology studies to identify homologs for S. cerevisiae and S. pombe proteins involved in these pathways. Based on these observations, we have reconstructed probable A. nidulans signal transduction pathways.



Methods

We have used bioinformatics techniques for identifying new putative signal transduction proteins in the model fungus A. nidulans. As the genomes of A. nidulans [12], S. cerevisiae [19] and S. pombe [20] are now available, we performed sequence homology studies using the stand-alone BLAST program (Basic Local Alignment Search Tools) [21, 22]. We downloaded the A. nidulans genome [12] and re-formatted the FASTA database to a BLAST database using the "formatdb" program [23]. We compared the sequence homology of the signal transduction proteins in S. cerevisiae and S. pombe with the A. nidulans database created using "blastp" program [21, 22]. This method is very effective as it is fast and the reliability of the results is scored based on statistical techniques [21, 24].


Significance threshold

Percentage sequence identity and statistical score, such as expectation value (E-value) of BLAST or FASTA, are widely used measures for sequence comparison [25]. As it has been well-established that use of a statistical score such as E-value of BLAST is superior to percentage sequence identity in detecting remote homology [26], we have used the expectation value from the BLAST output as the measure for identifying the significance of the homology. Expectation value is defined as the number of different alignments with scores equivalent to or better than raw alignment score that are expected to occur in a database search by chance [23]. The statistical significance threshold for reporting matches against database sequences was E = 0.001. We have chosen stringent E-value thresholds for assigning homology based on the previous literature [25, 27, 28]. For example, Poggeler [18] had identified a probable pheromone response pathway in A. fumigatus using sequence homology with the S. cerevisiae pheromone response pathway. They had used an E-value of 10-6 as the threshold to assign homologs for all types of proteins. Other studies have used thresholds of E = 10-10 [29] and E = 10-15 [27]. We have taken a more conservative approach and we have used E = 10-15 as the threshold for assigning homology for only general proteins. Tian et al. used enzyme function conservation to derive the threshold of sequence identity necessary to transfer function from an enzyme of known function to an unknown protein to be E = 10-50. Hence, our significance threshold for enzymes in general was chosen to be E = 10-50 [25]. For identification of kinase homologs Wang et al. [28] showed an even lower E-value of 10-80 was required. Thus for kinase identifications we have used E = 10-80. We chose the threshold of E = 10-6 used by Poggeler [18] as a less conservative approach.



Results and discussion


Validation

For validation, we constructed a BLAST database using the S. cerevisiae FASTA genome. Using this database, we identified the S. cerevisiae homologs for the 18 annotated (4 experimentally and 14 electronically) A. nidulans kinases [30]. The results obtained are shown in Table 1. As can be seen, S. cerevisiae homologs of the A. nidulans kinases either have the same function or belong to the same family depending on the expectation value. For example, AnkA and its homolog Swe1p belong to the same family; KapR and its homolog Sra1p have the same function. We also observe that the homolog for NimA is Kin3p, a S. cerevisiae protein similar to NimA [31]. The only contradictory result is the high homology between Pho85p and PhoA(M1)/PhoA(M47). Even though they show high sequence homology, Pho85p is involved in phosphate and glycogen metabolism [32] while PhoA(M1)/PhoA(M47) are involved in phosphate limited growth [33]. Also, our results indicate that S. cerevisiae Hog1p shows significant homology with two A. nidulans kinases, namely, HogA and MpkC. This is probably because Hog1p in S. cerevisiae is involved in osmoregulation [34] as HogA of A. nidulans [17], and is induced under stress conditions [32]. As the function of MpkC is not yet clearly understood, we predict that MpkC may be involved in a stress response mechanism.

Based on the BLAST significance threshold, we have identified the proteins that are likely involved in A. nidulans signal transduction pathways. The results are shown in the Tables 2 through 9. Using this information, we have used the well-established models of S. cerevisiae [35] and S. pombe [14] to construct the A. nidulans signal transduction.

Table 1: Validation based on comparison of annotated kinases in A. nidulans with their S. cerevisiae homolog
(Click for viewing complete Table 1)



Filamentation-invasion pathway

Under specific culture conditions, diploid S. cerevisiae cells undergo a dimorphic switch and differentiate to form pseudohyphae, growing as extended and connected cells. Starvation of nitrogen has been identified to induce pseudohyphal differentiation [48]. In S. cerevisiae, the switch to pseudohyphal differentiation probably involves two signal transduction pathways [14, 35]. One of them is the filamentation-invasion pathway [14, 35] and the other is cAMP-Protein Kinase A pathway [14].

Filamentation-invasion pathway (Fig. 1) contains a mitogen-activated protein kinase (MAPK) cascade. These cascades generally contain three protein kinases that act in series, namely mitogen-activated protein kinase kinase kinase (MAPKKK or MEKK), mitogen-activated protein kinase kinase (MAPKK or MEK) and mitogen-activated protein kinase (MAPK). When the cascade is activated, MEKK phosphorylates the MEK, which in turn phosphorylates the MAPK [35]. These MAPK cascades regulate transcription factors by MAPK-mediated phosphorylation [14, 15, 35]. The components of the MAPK cascade of the S. cerevisiae filamentation-invasion pathway are Ste11p (MEKK), Ste7p (MEK) and Kss1p (MAPK). Ste20p acts as the upstream kinase. The MAPK cascade mediates signal transduction from two small GTP binding proteins, Ras2p and Cdc42p [35].

Our results show that 11 of the 20 S. cerevisiae proteins involved in filamentation-invasion pathway show significant A. nidulans homologs (Table 2, Fig. 1). Furthermore, an additional 6 S. cerevisiae filamentation-invasion proteins have homologs in A. nidulans under less conservative thresholds (Table 2, Fig. 1). These results imply A. nidulans is likely to have a pathway similar to S. cerevisiae filamentation-invasion pathway.

As A. nidulans is a filamentous fungus, it seems likely that the pathway homologous to the S. cerevisiae filamentation-invasion pathway would trigger a different cellular response. We predict that the A. nidulans pathway may be involved in the regulation of asexual sporulation. In support of this prediction we note a number of consistent observations from the literature. Asexual sporulation in A. nidulans is known to be regulated by two different mechanisms similar to the pseudohyphal differentiation in S. cerevisiae [14]. Initiation of conidiation is contingent upon cells having undergone a defined period of growth to reach developmental competence [49]. This is an intrinsically programmed process that is independent of external changes in growth medium [50]. Although conidiation largely occurs in a medium-independent manner, it can be induced in hyphae grown in submerged culture by nutrient (carbon or nitrogen)-limiting condition [51]. Conidiation in confluent cultures of A. nidulans depends upon inoculum density; and higher initial cell densities promote sporulation [33]. Induction of morphogenesis involves temporal and spatial regulation of several hundred genes, many of which function in conidiophore assembly or conidiospore differentiation [14]. The sequential expression of bristle (brlA), abacus (abaA), and wet-white conidia (wetA) establishes a central regulatory pathway (brlA --> abaA --> wetA) required for the transition from vegetative growth to asexual reproduction [52, 53]. We observe that the S. cerevisiae transcription factor Tec1p, which is involved in the regulation of filamentation-invasion proteins is homologous to AbaA in A. nidulans (Table 2). AbaA controls the temporal and spatial specificity of A. nidulans [34]. Expression of abaA leads to activation of brlA and wetA and cessation of vegetative growth [34]. AbaA also regulates the expression of numerous sporulation specific genes [34]. This implies that the signal transduction pathway homologous to filamentation-invasion pathway may regulate the A. nidulans asexual sporulation.

Table 2: S. cerevisiae filamentation-invasion pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 2)



Figure 1: In silico reconstruction of the S. cerevisiae filamentation-invasion pathway [35] in A. nidulans. A. nidulans sequence homologs of S. cerevisiae proteins are listed in square brackets. The proteins listed in bold are homologous to more than one S. cerevisiae protein. Homology of the underlined proteins is less than the significance threshold. Proteins with just the annotation numbers are hypothetical proteins. --> activation, --| inhibition, ? No homologs identified.



cAMP-PKA Pathway

The S. cerevisiae cAMP-PKA pathway plays a vital role in regulating carbon metabolism and cell cycle progression and also functions in parallel with the filamentation-invasion pathway to regulate pseudohyphal growth [14]. This nutrient-sensing pathway involves a novel G protein-coupled receptor (GPCR), G-proteins, Ga protein Gpa2p, adenylyl cyclase, cAMP, and cAMP-dependent protein kinase A (PKA) [54, 55]. In S. cerevisiae, Gpa2p regulates cAMP levels [56, 57], which then modulate the activity of cAMP-dependent protein kinase (PKA); Tpk2p, one of the three catalytic subunits of PKA, specifically controls the filamentation response [58, 59].

Except for the S. cerevisiae GPCR, all other cAMP-PKA pathway proteins had homologs in A. nidulans (Table 3, Fig. 2). The absence of an A. nidulans homolog for the GPCR is not surprising as this gene shows no sequence homology with a functional homolog from S. pombe. Hence, A. nidulans may also have a GPCR functional homolog that shows little sequence homology to Gpr1.

It has been established that in A. nidulans, a pathway involving a heterotrimeric G protein a subunit encoded by fadA regulates mycelium proliferation and antagonizes conidiophore development [60]. Our results indicate that Gba1 encoded by fadA is homologous to Gpa2 of S. cerevisiae which is involved in the cAMP-mediated signal transduction mechanism [54, 55]. This implies that the cAMP-PKA pathway in A. nidulans functions in parallel with the filamentation-invasion pathway homolog to regulate conidiation. Also, protein kinase A regulatory subunit has been isolated in A. nidulans [30]. The possible involvement of a cAMP-dependent protein kinase, PkaC, in sporulation has been studied in the related fungus Aspergillus niger [61]. The cAMP signaling pathway in the filamentous fungus Neurospora crassa is known to regulate the conidiation, morphogenesis, mating and stress tolerance [14]. Hence the A. nidulans cAMP signaling pathway may regulate similar cellular responses. Together, these results imply the cAMP-PKA pathway in A. nidulans is involved in regulation of conidiation.

Table 3: S. cerevisiae cAMP-PKA pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 3)


Figure 2: In silico reconstruction of the S. cerevisiae cAMP-PKA pathway [37, 54, 62] in A. nidulans. A. nidulans sequence homologs of S. cerevisiae proteins are listed in square brackets. | | | | Binding. Notations as in Fig. 1.


Since we observed the lack of sequence homolog for GPCR in A. nidulans which is similar to the S. pombe, we decided to identify the A. nidulans homologs for S. pombe cAMP-PKA pathway proteins. The results are shown in Table 4. As we can see from the results, the central part of the signaling mechanism is very well conserved between A. nidulans and S. pombe. This comparison between S. pombe and A. nidulans genome also doesn’t show a sequence homolog for the glucose receptor protein. This indicates that the glucose receptor protein molecule in A. nidulans doesn’t show a sequence homology with glucose receptor proteins in both S. pombe and S. cerevisiae. It has been established that the cAMP-PKA pathways in S. pombe and S. cerevisiae are very similar [14]. The lack of sequence homology of the glucose receptor proteins can mainly be attributed to asparagine-rich third intracellular loop found in S. cerevisiae Gpr1 and its absence in S. pombe Git3 . Hence, we can predict that the cAMP-PKA pathway in A. nidulans is also similar to the cAMP-PKA pathways in S. cerevisiae and S. pombe but has a glucose receptor protein which doesn’t show significant sequence homology. From our results in Table 3, we find that the all the three cAMP-dependent protein kinase subunits show only one sequence homolog in A. nidulans. The same protein from A. nidulans shows homology with the cAMP-dependent protein kinase catalytic subunit of S. pombe. This probably means that A. nidulans has only one cAMP-dependent protein kinase catalytic subunit similar to the S. pombe genome and not three cAMP-dependent protein kinase catalytic subunits as in the case of S. cerevisiae. Due to the probable high similarity of A. nidulans cAMP-PKA pathway to the S. pombe cAMP-PKA pathway we have reconstructed the S. pombe cAMP-PKA pathway in A. nidulans. This is shown in Fig. 3.

Table 4: S. pombe cAMP-PKA pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 4)


Figure 3: In silico reconstruction of S. pombe cAMP-PKA pathway [14, 63] in A. nidulans. A. nidulans sequence homologs of S. pombe proteins are listed in square brackets. Notations as in Fig. 1.



Pheromone response pathway

In S. cerevisiae, the pheromone response pathway induces polarized growth towards the mating partner and cell cycle arrest [35]. In addition, expression of proteins needed for cell adhesion, cell fusion and nuclear fusion are increased [35]. While this is not a nutrient sensing pathway, we include it here because of its similarity to the filamentation-invasion pathway.

The S. cerevisiae pheromone response pathway is also a MAPK cascade, which uses components common to the filamentation-invasion pathway [14]. The Ste2p and Ste3p receptors sense the pheromone and activate the MAPK cascade through the upstream protein kinase Ste20p [35]. Ste20p activates Ste11p (MEKK), which phosphorylates Ste7p (MEK), which in turn phosphorylates Fus3p (MAPK). S. cerevisiae Fus3p regulates cell cycle arrest and the transcription of pheromone response genes [35]. In addition the pathway is controlled by a scaffold protein, Ste5p [35].

Our results show that 12 of the 23 S. cerevisiae proteins involved in the pheromone response pathway have significant A. nidulans homologs (Table 5, Fig. 4). In addition, eight other S. cerevisiae pheromone response pathway proteins had less conservative homologs. One of these was the scaffold protein Ste5p. This may imply that the A. nidulans pheromone response pathway is more closely related to that in S. pombe, where the MAPK pathway doesn’t have any scaffold protein [15]. The pheromone response pathway in A. nidulans likely regulates the transcription of genes involved in mating.

Table 5: S. cerevisiae pheromone response pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 5)


Figure 4: In silico reconstruction of the S. cerevisiae pheromone response pathway [35] in A. nidulans. A. nidulans sequence homologs of S. cerevisiae proteins are listed in square brackets. Notation as in Fig. 1.


Since A. nidulans didn’t show a significant sequence homolog for the scaffold protein of S. cerevisiae, we wanted to compare the similarity between the pheromone response pathway in A. nidulans and S. pombe. The results obtained by comparing the pheromone response proteins in S. pombe with the A. nidulans genome is shown in Table 6. The results show that the S. pombe proteins which are involved in the MAPK cascade of pheromone response pathway [15] have significant homologs in A. nidulans. This appears to confirm our prediction that the MAPK cascade of the A. nidulans pheromone pathway may be more closely related to S. pombe than to S. cerevisiae. Our results in Table 6 show that almost all the proteins that are involved in S. pombe pheromone response pathway have significant homologs in A. nidulans. In contrast, most of the S. pombe proteins involved in mating and meiosis do not have significant homologs in A. nidulans. The probable A. nidulans pheromone response pathway is shown in Fig. 5.

Table 6: S. pombe pheromone response pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 6)


Figure 5: In silico reconstruction of S. pombe pheromone response pathway [14, 15] in A. nidulans. A. nidulans sequence homologs of S. pombe proteins are listed in square brackets. Notation as in Fig. 1.



Cell integrity pathway

In S. cerevisiae, the cell integrity pathway regulates the expression of cell wall genes [35] and is stimulated by a variety of factors which include cell cycle regulation, heat stress, hypotonic stress and nutrient availability [35]. This pathway is also called as the Protein Kinase C pathway as it is mediated by members of the family of phospholipid-dependent, serine/threonine-specific protein kinases called protein kinase C. The cell integrity pathway in S. cerevisiae is highly conserved and it is very similar to the mammalian cell integrity pathway [35].

A number of membrane proteins, namely, Wsc1p, Wsc2p, Wsc3p, Wsc4p and Mid2p, have been reported to provide input signals to the cell integrity pathway [35, 64]. These signals activate the MAPK cascade which is mediated by the upstream protein kinase protein kinase C [35]. This kinase activates Bck1p (MEKK), which phosphorylates the functionally redundant threonine/tyrosine kinases Mkk1p and Mkk2p (MEK) [64, 65], which in turn phosphorylates Slt2p (MAPK). Slt2p regulates the transcription of cell wall genes [35]. In addition PKC may be activated through the Tor2p-Rom2p-Rho1p route [64].

Our results show that 15 of the 22 S. cerevisiae cell integrity pathway proteins have significant A. nidulans homologs (Table 7, Fig. 6). In addition, five other S. cerevisiae cell integrity pathway proteins show less conservative A. nidulans homologs.

Homolog of Rho1p in A. nidulans, RhoA has been isolated [66]. Pkc1 homolog Kpc1 has been cloned and characterized from A. niger [67] and it has been identified in A. nidulans based on similarity [34]. Our results show that AN4189.2 shows homology with both Mkk1p and Mkk2p, which act as the mitogen activated protein kinase kinases. This could be due to two reasons. First, since Mkk1p and Mkk2p have redundant functions [65] A. nidulans may use only one gene for the same purpose. Another possible explanation is that both Mkk1p and Mkk2p have high sequence homology and result in the same homolog in A. nidulans although there are two proteins present. Cell cycle regulation of cell integrity pathway is mediated by the activity of the cyclin-dependent kinase Cdc28p (CDK in Fig. 6, A. nidulans homolog is CDC2_EMENI) in complex with G1 cyclins. The downstream substrates of Slt2p are Rlm1p and SBF which regulate the transcription of cell wall related genes [35].

Table 7: S. cerevisiae cell integrity pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 7)


Figure 6: In silico reconstruction of the S. cerevisiae cell integrity pathway [35] in A. nidulans. A. nidulans sequence homologs of S. cerevisiae proteins are listed in square brackets. CDK – cyclin dependent kinase. Notation as in Fig. 1.



Target of rapamycin pathway

The target of rapamycin (TOR) pathway is one of the central nutrient-sensing signal transduction pathways in eukaryotic cells [68] and is conserved from yeasts to humans [69]. In S. cerevisiae, TOR responds to nutrient limitation and mediates temporal control of cell growth by positively regulating anabolic processes such as translation, transcription and ribosome biosynthesis, and by negative regulation of catabolic processes such as protein and RNA degradation [70]. S. cerevisiae TOR also mediates spatial control of cell growth by regulating polarization of the actin cytoskeleton [71-73]. Furthermore TOR regulates the expression of proteins required for adaptation to starvation of carbon or nitrogen and utilization of poor nitrogen sources [70, 74]. Recent studies have shown that TOR signaling in S. cerevisiae is required for proper pseudohyphal development [69, 74]. However, it is not clear how TOR signaling affects the filamentation [74]. In S. cerevisiae, TOR also plays a key role in regulating autophagy [75-77] based on the nutrient availability.

In S. cerevisiae, TOR signaling consists of two TOR protein kinases namely Tor1p and Tor2p [68]. S. cerevisiae TOR kinases are activated by nutrient availability and are inactivated by either starvation or the presence of rapamycin [74]. Rapamycin inhibits TOR function by forming a complex with the prolyl isomerase FKBP12 and interacting with the TOR kinases at a conserved FRB (FKBP12-rapamycin binding) domain distinct from the kinase domain [74]. Other members of the TOR signaling cascade include type 2A protein phosphatases (PP2A), a type 2A related phosphatase Sit4p and Tap42p, a direct target of the S. cerevisiae TOR proteins [78]. The complex formed between Tap42p and the phosphatases play a central role in TOR signaling cascade [78]. TOR regulates the formation of these complexes by reversible phosphorylation and dephosphorylation of Tap42p [78].

We identified A. nidulans homologs for all the S. cerevisiae TOR signaling pathway proteins (Table 8, Fig. 7). We found only one homolog for both Tor1p and Tor2p of S. cerevisiae in A. nidulans genome. This is probably because A. nidulans is higher on the evolution ladder and it has been observed that in general higher organisms have a single TOR kinase [79-81].

Table 8: S. cerevisiae TOR pathway proteins and their A. nidulans homologs
(Click for viewing complete Table 8)


Figure 7: In silico reconstruction of the S. cerevisiae TOR signaling pathway [78] in A. nidulans. A. nidulans sequence homologs of S. cerevisiae proteins are listed in square brackets. Reversible arrows represent reversible phosphorylation/dephosphorylation. Notation as in Fig. 1.



Autophagy

Here we concentrate on one of the most important cellular responses to TOR signaling, namely the regulation of autophagy. Autophagy is a cellular response to nutrient starvation and is essential for cell survival under nutrient depleted conditions [82]. During autophagy, the bulk of the cytoplasmic components are non-selectively enclosed within a double-membrane structure called autophagosome, and this is transported in to the vacuole/lysosome to be degraded [83-85]. Autophagy is also characterized as type II programmed cell death and plays a critical role during development [77]. Autophagic dysfunction is also associated with various diseases [77]. For example, investigators have found a correlation between defects in autophagy and carcinogenesis [86]. In addition, defective autophagy is also correlated with numerous neurodegenerative diseases such as Huntington’s and Parkinson’s diseases [75].

The genes essential for autophagy, apg and aut genes, have been isolated [87, 88]. It has been shown that the apg mutant cells lose viability under nutrient starvation [87]. There are two types of autophagy: microautophagy and macroautophagy [77]. Microautophagy involves uptake of cytoplasm or whole organelles directly at the vacuole membrane [77]. While macroautophagy involves a sequestration event that is initiated separate from the vacuole membrane [77]. Macroautophagy can be broken down into a few basic steps: signaling, sequestration of cytoplasm, completion of vesicle formation, targeting of the completed vesicle to the lysosome/vacuole followed by docking and fusion, and breakdown (Fig. 8). The proteins involved in each of these steps have been identified in S. cerevisiae (Fig. 8). We note that homologs of the S. cerevisiae autophagy genes have been discovered in a wide range of cells including human, Drosophila, Dictyostelium, Arabidopsis, and Caenorhabditis elegans [89].

Our results show that 35 out of the 51 S. cerevisiae autophagy proteins have significant A. nidulans homologs (Table 9, Fig. 8). In addition, we observe that 9 S. cerevisiae autophagy proteins have less conservative homologs in A. nidulans. These results imply that A. nidulans autophagy is similar to S. cerevisiae autophagy. We also observe that with the exception of Apg7, none of the proteins involved in Apg12-Apg5 conjugation had a significant A. nidulans homolog. Whereas all the proteins involved in Aut7 lipid conjugation had a significant A. nidulans homolog.

Table 9: S. cerevisiae autophagy proteins and their A. nidulans homologs
(Click for viewing complete Table 9)


Figure 8: Schematic of the four major steps in the macroautophagy pathway in S. cerevisiae, adapted from [82]. Listed proteins are from S. cerevisiae [82], and correspond to each step in the pathway. Using these S. cerevisiae proteins, we have performed a BLAST search and identified homologous proteins in our model fungal strain (A. nidulans A4). Notation as in Fig. 1.



TOR-mediated regulation of autophagy

In S. cerevisiae, the rapamycin treatment mimics nutrient starvation [90] and induces autophagy showing that the TOR signaling cascade is involved in the regulation of autophagy. The mechanism by which the TOR pathway regulates autophagy is well-understood in S. cerevisiae [77] and appears to act by at least by two mechanisms [77].

A putative nutrient sensor on the cell surface appears to transduce the TOR pathway (Fig. 9). Under starvation conditions, Tor2 protein is inhibited, resulting in the dephosphorylation of Tap42 and its subsequent displacement from the protein phosphatase 2A catalytic subunits (i.e. Sit4) and their resulting activation [91]. These activated phosphatases control the activity of transcriptional activators (Fig. 9). Another mechanism of TOR regulation involves protein kinase Apg1 and phosphoprotein Apg13 [92], which are both involved in vesicle biogenesis. Under nutrient-rich conditions, Apg13 is hyper-phosphorylated and interacts weakly with Apg1. Inhibition of TOR under nutrient starvation leads to dephosphorylation of Apg13, thereby increasing its affinity to Apg1. Interaction with Apg13 is thought to regulate Apg1 kinase activity (Fig. 9).

Significant homologs for Sit4 (E-value 10-108) and TOR (E-value 0.00) are seen in A. nidulans genome (Table 8). Our results show the presence of a significant homolog for Apg1 kinase, but the homolog for Apg13 shows less homology (Table 9). This is not surprising because homologs of proteins that interact with Apg1, including Apg13, are absent in higher eukaryotes [77]. Hence, we assume this is the reason for the lack of a significant A. nidulans homolog for S. cerevisiae Apg13.


Figure 9: Schematic of Tor mediated regulation of autophagy, adapted from [76]. The sequence homologs of the S. cerevisiae protein in A. nidulans genome are listed in the square brackets. Notation as in Fig. 1.




Conclusions

Our results show that most of the signal transduction proteins in S. cerevisiae and S. pombe have significant homologs in A. nidulans indicating that the signal transduction pathways are well conserved. Our work has helped in identifying the probable signal transduction proteins in A. nidulans and in reconstructing the possible pathways using in silico techniques. These identifications can be used for further analysis of the signal transduction pathways. Based on the results obtained from our studies, further studies on deletion mutants can be performed to confirm function and physiological responses in A. nidulans under activation of the different signal transduction pathways. As we have seen from earlier studies [17, 18], the results obtained from computational techniques are highly reliable and have been experimentally validated. But our results had a few discrepancies. The most important was that some of the A. nidulans proteins showed significant homology with more than one yeast protein. For example, the mitogen-activated protein kinase in S. cerevisiae filamentation-invasion pathway, pheromone response pathway and cell integrity pathway have the same A. nidulans homolog. This is likely because the MAPKs have high sequence homology amongst themselves. This implies that some other A. nidulans proteins may act as the functional homolog of the S. cerevisiae MAPKs. On the whole, this study has provided insight on the signal transduction pathways involved in A. nidulans, a model filamentous fungus and will aid in understanding similar pathways in other industrially and medically significant Aspergilli.



Acknowledgements

We would like to thank Broad Institute of MIT and Harvard for allowing us to use the A. nidulans genome for performing in silico analysis of the genome sequence.




References


  1. MacCabe, A. P., Orejas, M., Tamayo, E. N., Villanueva, A. and Ramon, D. (2002). Improving extracellular production of food-use enzymes from Aspergillus nidulans. J. Biotechnol. 96, 43-54.

  2. Zidwick, M. J. (1991). Organic Acids. In: Biotechnology of Filamentous Fungi: Technology and Products, Finkelstein DB, B. C. (eds.), Butterworth-Heinemann, Boston, pp. 303-334.

  3. Arora, D. K., et al. (1992). Introduction. In: Handbook of Applied Mycology, D. K. Arora, R. P. E., K. G. Mukerji (eds.), Marcel Dekker, New York, pp. 3-4.

  4. Rambosek, J. A. (1991). Molecular Biology and Biochemistry of b-Lactam Antibiotics. In: Biotechnology of Filamentous Fungi: Technology and Products, Finkelstein DB, B. C. (eds.), Butterworth-Heinemann, Boston, pp. 221-240.

  5. Masurekar, P. S. (1991). Therapeutic Metabolites. In: Biotechnology of Filamentous Fungi: Technology and Products, Finkelstein DB, B. C. (eds.), Butterworth-Heinemann, Boston, pp. 241-302.

  6. O'Donnell, K. and Peterson, S. W. (1991). Isolation, Preservation, and Taxonomy. In: Biotechnology of Filamentous Fungi: Technology and Products, Finkelstein DB, B. C. (eds.), Butterworth-Heinemann, Boston, pp. 7-40.

  7. White, S., McIntyre, M., Berry, D. R. and McNeil, B. (2002). The autolysis of industrial filamentous fungi. Crit. Rev. Biotechnol. 22, 1-14.

  8. Garber, G. (2001). An overview of fungal infections. Drugs 61, Suppl. 1, 1-12.

  9. Sternberg, S. (1994). The emerging fungal threat. Science 266, 1632-1634.

  10. Toxic Mold, Black Mold, Mold Problem and Mold Damage Information Center. http://www.mold.ph/aspergillus.htm.

  11. Mold Help. http://www.mold-help.org/aspergillus.htm.

  12. Aspergillus nidulans Database. http://www.genome.wi.mit.edu/annotation/fungi/aspergillus/.

  13. Medical Dictionary Database. http://www.books.md/index.html.

  14. Lengeler, K. B., Davidson, R. C., D'souza, C., Harashima, T., Shen, W. C., Wang, P., Pan, X., Waugh, M. and Heitman, J. (2000). Signal transduction cascades regulating fungal development and virulence. Microbiol. Mol. Biol. Rev. 64, 746-785.

  15. Banuett, F. (1998). Signalling in the yeasts: an informational cascade with links to the filamentous fungi. Microbiol. Mol. Biol. Rev. 62, 249-274.

  16. Hunter, T. and Plowman, G. D. (1997). The protein kinases of budding yeast: six score and more. Trends Biochem. Sci. 22, 18-22.

  17. Han, K. H. and Prade, R. A. (2002). Osmotic stress-coupled maintenance of polar growth in Aspergillus nidulans. Mol. Microbiol. 43, 1065-1078.

  18. Poggeler, S. (2002). Genomic evidence for mating abilities in the asexual pathogen Aspergillus fumigatus. Curr. Genet. 42, 153-160.

  19. The S. cerevisae genome. ftp://ftp.ncbi.nlm.nih.gov/blast/db/FASTA/.

  20. The S. pombe Genome Project. http://www.sanger.ac.uk/Projects/S_pombe/.

  21. Altschul, S. F. and Gish, W. (1996). Local alignment statistics. Methods Enzymol. 266, 460-480.

  22. Altschul, S. F., Gish, W., Miller, W., Myers, E. W. and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403-410.

  23. Basic Local Alignment Search Tools. http://www.ncbi.nlm.nih.gov/blast/.

  24. Altschul, S. F., Boguski, M. S., Gish, W. and Wootton, J. C. (1994). Issues in searching molecular sequence databases. Nat. Genet. 6, 119-129.

  25. Tian, W. and Skolnick, J. (2003). How well is enzyme function conserved as a function of pairwise sequence identity? J. Mol. Biol. 333, 863-882.

  26. Brenner, S. E., Chothia, C. and Hubbard, T. J. (1998). Assessing sequence comparison method with reliable structurally identified distant evolutionary relationships. Proc. Natl. Acad. Sci. USA 95, 6073-6078.

  27. Sorrells, M. E., La Rota, M., Bermudez-Kandianis, C. E., Greene, R. A., Kantety, R., Munkvold, J. D., Miftahudin, Mahmoud, A., Ma, X., Gustafson, P. J., Qi, L. L., Echalier, B., Gill, B. S., Matthews, D. E., Lazo, G. R., Chao, S., Anderson, O. D., Edwards, H., Linkiewicz, A. M., Dubcovsky, J., Akhunov, E. D., Dvorak, J., Zhang, D., Nguyen, H. T., Peng, J., Lapitan, N. L., Gonzalez-Hernandez, J. L., Anderson, J. A., Hossain, K., Kalavacharla, V., Kianian, S. F., Choi, D. W., Close, T. J., Dilbirligi, M., Gill, K. S., Steber, C., Walker-Simmons, M. K., McGuire, P. E. and Qualset, C. O. (2003). Comparative DNA sequence analysis of wheat and rice genomes. Genome Res. 13, 1818-1827.

  28. Wang, D., Harper, J. F. and Gribskov, M. (2003). Systematic trans-genomic comparison of protein kinases between Arabidopsis and Saccharomyces cerevisiae. Plant Physiol. 132, 2152-2165.

  29. Gardiner, J. and Marc, J. (2003). Putative microtubule-associated proteins from Arabidopsis genome. Protoplasma 222, 61-74.

  30. Smith, C. M., Shindyalov, I. N., Veretnik, S., Gribskov, M., Taylor, S. S., Ten Eyck, L. F. and Bourne, P. E. (1997). The protein kinase resource. Trends Biochem. Sci. 22, 444-446.

  31. Schweitzer, B. and Philippsen, P. (1992). NPK1, a nonessential protein kinase gene in Saccharomyces cerevisiae with similarity to Aspergillus nidulans nimA. Mol. Gen. Genet. 234, 164-167.

  32. Saccharomyces Genome Database. http://www.yeastgenome.org/.

  33. Bussink, H. J. and Osmani, S. A. (1998). A cyclin-dependent kinase family member (PHOA) is required to link developmental fate to environmental conditions in Aspergillus nidulans. EMBO J. 17, 3990-4003.

  34. Swiss-Prot protein knowledgebase. http://us.expasy.org/sprot.

  35. Gustin, M. C., Albertyn, J., Alexander, M. and Davenport, K. (1998). MAP Kinase Pathways in the Yeast Saccharomyces cerevisiae. Microbiol. Mol. Biol. Rev. 62, 1264-1300.

  36. Hofmann, A. F. and Harris, S. D. (2000). The Aspergillus nidulans uvsB gene encodes an ATM-related kinase required for multiple facets of the DNA damage response. Genetics 154, 1577-1586.

  37. Kato, R. and Ogawa, H. (1994). An essential gene, ESR1, is required for mitotic cell growth, DNA repair and meiotic recombination in Saccharomyces cerevisiae. Nucleic Acids Res. 22, 3104-3112.

  38. Shimizu, K. and Keller, N. P. (2001). Genetic involvement of a cAMP-dependent protein kinase in a G protein signaling pathway regulating morphological and chemical transitions in Aspergillus nidulans. Genetics 157, 591-600.

  39. Osmani, A. H., van Peij, N., Mischke, M., O'Connell, M. J. and Osmani, S. A. (1994). A single p34cdc2 protein kinase (encoded by nimXcdc2) is required at G1 and G2 in Aspergillus nidulans. J. Cell Sci. 107, 1519-1528.

  40. Efimov, V. P. and Morris, N. R. (1998). A screen for dynein synthetic lethals in Aspergillus nidulans identifies spindle assembly checkpoint genes and other genes involved in mitosis. Genetics 149, 101-116.

  41. Roberts, B. T., Farr, K. A. and Hoyt, M. A. (1994). The Saccharomyces cerevisiae checkpoint gene BUB1 encodes a novel protein kinase. Mol. Cell. Biol. 14, 8282-8291.

  42. Bussink, H. J. and Osmani, S. A. (1999). A mitogen-activated protein kinase (MPKA) is involved in polarized growth in the filamentous fungus, Aspergillus nidulans. FEMS Microbiol. Lett. 173, 117-125.

  43. Mazzoni, C., Zarov, P., Rambourg, A. and Mann, C. (1993). The SLT2 (MPK1) MAP kinase homolog is involved in polarized cell growth in Saccharomyces cerevisiae. J. Cell Biol. 123, 1821-1833.

  44. Joseph, J. D. and Means, A. R. (2000). Identification and characterization of two Ca2+/CaM-dependent protein kinases required for normal nuclear division in Aspergillus nidulans. J. Biol. Chem. 275, 38230-38238.

  45. Kraus, P. R. and Harris, S. D. (2001). The Aspergillus nidulans snt genes are required for the regulation of septum formation and cell cycle checkpoints. Genetics 159, 557-569.

  46. Booher, R. N., Deshaies, R. J. and Kirschner, M. W. (1993). Properties of Saccharomyces cerevisiae wee1 and its differential regulation of p34CDC28 in response to G1 and G2 cyclins. EMBO J. 12, 3417-3426.

  47. Bruno, K. S., Morrell, J. L., Hamer, J. E. and Staiger, C. J. (2001). SEPH, a Cdc7p orthologue from Aspergillus nidulans, functions upstream of actin ring formation during cytokinesis. Mol. Microbiol. 42, 3-12.

  48. Gimeno, C. J., Ljungdahl, P. O., Styles, C. A. and Fink, G. R. (1992). Unipolar cell divisions in the yeast S. cerevisiae lead to filamentous growth: regulation by starvation and RAS. Cell 68, 1077-1090.

  49. Axelrod, D. E., Gealt, M. and Pastushok, M. (1973). Gene control of developmental competence in Aspergillus nidulans. Dev. Biol. 34, 9-15.

  50. Pastushok, M. and Axelrod, D. E. (1976). Effect of glucose, ammonium and media maintenance on the time of conidiophore initiation by surface colonies of Aspergillus nidulans. J. Gen. Microbiol. 94, 221-224.

  51. Skromne, I., Sanchez, O. and Aguirre, J. (1995). Starvation stress modulates the expression of the Aspergillus nidulans brlA regulatory gene. Microbiology 141, 21-28.

  52. Adams, T. H., Boylan, M. T. and Timberlake, W. E. (1988). brlA is necessary and sufficient to direct conidiophore development in Aspergillus nidulans. Cell 51, 353-362.

  53. Marshall, M. A. and Timberlake, W. E. (1991). Aspergillus nidulans wetA activates spore-specific gene expression. Mol. Cell. Biol. 11, 55-62.

  54. Thevelein, J. M. and de Winde, J. H. (1999). Novel sensing mechanisms and targets for the cAMP-protein kinase A pathway in the yeast Saccharomyces cerevisiae. Mol. Microbiol. 33, 904-918.

  55. Müller, D., Exler, S., Aguilera-Vázquez, L., Guerrero-Martín, E. and Reuss, M. (2003). Cyclic AMP mediates the cell cycle dynamics of energy metabolism in Saccharomyces cerevisiae. Yeast 20, 351-367.

  56. Lorenz, M. C. and Heitman, J. (1997). Yeast pseudohyphal growth is regulated by GPA2, a G protein a homolog. EMBO J. 16, 7008-7018.

  57. Lorenz, M. C., Pan, X., Harashima, T., Cardenas, M. E., Xue, Y., Hirsch, J. P. and Heitman, J. (2000). The G protein-coupled receptor Gpr1 is a nutrient sensor that regulates pseudohyphal differentiation in Saccharomyces cerevisiae. Genetics 154, 609-622.

  58. Pan, X. and Heitman, J. (1999). Cyclic AMP-Dependent protein kinase regulates pseudohyphal differentiation in Saccharomyces cerevisiae. Mol. Cell. Biol. 19, 4874-4887.

  59. Robertson, L. S. and Fink, G. R. (1998). The three yeast A kinases have specific signaling functions in pseudohyphal growth. Proc. Natl. Acad. Sci. USA 95, 13783-13787.

  60. Yu, J. H., Wieser, J. and Adams, T. H. (1996). The Aspergillus FlbA RGS domain protein antagonizes G protein signaling to block proliferation and allow development. EMBO J. 15, 5184-5190.

  61. Bencina, M., Panneman, H., Ruijter, G. J., Legisa, M. and Visser, J. (1997). Characterization and overexpression of the Aspergillus niger gene encoding the cAMP-dependent protein kinase catalytic subunit. Microbiology 143, 1211-1220.

  62. Kawasaki, L., Sanchez, O., Shiozaki, K. and Aguirre, J. (2002). SakA MAP kinase is involved in stress signal transduction, sexual development and spore viability in Aspergillus nidulans. Mol. Microbiol. 45, 1153-1163.

  63. D'Souza, C. A. and Heitman, J. (2001). Conserved cAMP signaling cascades regulate fungal development and virulence. FEMS Microbiol. Rev. 25, 349-364.

  64. Heinisch, J. J., Lorberg, A., Schmitz, H. P. and Jacoby, J. J. (1999). The protein kinase C-mediated MAP kinase pathway involved in maintenance of cellular integrity in Saccharomyces cerevisiae. Mol. Microbiol. 32, 671-680.

  65. Irie, K., Takase, M., Lee, K. S., Levin, D. E., Araki, H., Matsumoto, K. and Oshima, Y. (1993). MKK1 and MKK2, which encode Saccharomyces cerevisiae mitogen-activated protein kinase-kinase homologs, function in the pathway mediated by protein kinase C. Mol. Cell. Biol. 13, 3076-3083.

  66. Guest, G. M., Lin, X. and Momany, M. (2004). Aspergillus nidulans RhoA is involved in polar growth, branching, and cell wall synthesis. Fungal Genet. Biol. 41, 13-22.

  67. Morawetz, R., Lendenfeld, T., Mischak, H., Mühlbauer, M., Gruber, F., Goodnight, J., de Graaff, L. H., Visser, J., Mushinski, J. F. and Kubicek, C. P. (1996). Cloning and characterisation of genes (pkc1 and pkcA) encoding protein kinase C homologues from Trichoderma reesei and Aspergillus niger. Mol. Gen. Genet. 250, 17-28.

  68. Lorberg, A. and Hall, M. N. (2004). TOR: the first 10 years. Curr. Top. Microbiol. Immunol. 279, 1-18.

  69. Cutler, N. S., Pan, X., Heitman, J. and Cardenas, M. E. (2001). The TOR signal transduction cascade controls cellular differentiation in response to nutrients. Mol. Biol. Cell 12, 4103-4113.

  70. Cardenas, M. E., Cutler, N. S., Lorenz, M. C., Di Como, C. J. and Heitman, J. (1999). The TOR signaling cascade regulates gene expression in response to nutrients. Genes Dev. 13, 3271-3279.

  71. Bickle, M., Delley, P. A., Scmidt, A. and Hall, M. N. (1998). Cell wall integrity modulates RHO1 activity via the exchange factor ROM2. EMBO J. 17, 2235-2245.

  72. Schmelzle, T., Helliwell, S. B. and Hall, M. N. (2002). Yeast protein kinases and the RHO1 exchange factor TUS1 are novel components of the cell integrity pathway in yeast. Mol. Cell. Biol. 22, 1329-1339.

  73. Schmidt, A., Bickle, M., Beck, T. and Hall, M. N. (1997). The yeast phosphotidylinositol kinase homolog TOR2 activates RHO1 and RHO2 via the exchange factor ROM2. Cell 88, 531-542.

  74. Rohde, J. R. and Cardenas, M. E. (2004). Nutrient signaling through TOR kinases controls gene expression and cellular differentiation in fungi. Curr. Top. Microbiol. Immunol. 279, 53-72.

  75. Klionsky, D. J. and Emr, S. D. (2000). Autophagy as a regulated pathway of cellular degradation. Science 290, 1717-1721.

  76. Kamada, Y., Sekito, T. and Ohsumi, Y. (2004). Autophagy in yeast: a TOR-mediated response to nutrient starvation. Curr. Top. Microbiol. Immunol. 279, 73-84.

  77. Wang, C. W. and Klionsky, D. J. (2003). The molecular mechanism of autophagy. Mol. Med. 9, 65-76.

  78. Duvel, K. and Broach, J. R. (2004). The role of phosphatases in TOR signaling in yeast. Curr. Top. Microbiol. Immunol. 279, 19-38.

  79. Menand, B., Meyer, C. and Robaglia, C. (2004). Plant growth and the TOR pathway. Curr. Top. Microbiol. Immunol. 279, 97-113.

  80. Long, X., Muller, F. and Avruch, J. (2004). TOR action in mammalian cells and in Caenorhabditis elegans. Curr Top Microbiol. Immunol. 279, 115-138.

  81. Neufeld, T. P. (2004). Genetic analysis of TOR signaling in Drosophila. Curr. Top. Microbiol. Immunol. 279, 139-152.

  82. Abeliovich, H. and Klionsky, D. (2001). Autophagy in Yeast: Mechanistic Insights and Physiological Function. Microbiol. Mol. Biol. Rev. 65, 463-479.

  83. Takeshige, K., Baba, M., Tsuboi, S., Noda, T. and Ohsumi, Y. (1992). Autophagy in yeast demonstrated with proteinase-deficient mutants and conditions for its induction. J. Cell Biol. 119, 301-311.

  84. Baba, M., Takashige, K., Baba, N. and Ohsumi, Y. (1994). Ultrastructural analysis of the autophagic process in yeast: detection of autophagosmes and their characterization. J. Cell Biol. 124, 903-913.

  85. Klionsky, D. J. and Ohsumi, Y. (1999). Vacuolar import of proteins and organelles from cytoplasm. Annu. Rev. Cell. Dev. Biol. 15, 1-32.

  86. Liang, X. H., Jackson, S., Seaman, M., Brown, K., Kempkes, B., Hibshoosh, H. and Levine, B. (1999). Induction of autophagy and inhibition of tumorigenesis by beclin 1. Nature 402, 672-676.

  87. Tsukada, M. and Ohsumi, Y. (1993). Isolation and characterization of autophagy-defective mutants of Saccharomyces cerevisiae. FEBS Lett. 333, 169-174.

  88. Thumm, M., Egner, R., Koch, B., Schlumpberger, M., Straub, M., Veenhuis, M. and Wolf, D. H. (1994). Isolation of autophagocytosis mutants of Saccharomyces cerevesiae. FEBS Lett. 349, 275-280.

  89. Reggiori, F. and Klionsky, D. J. (2002). Autophagy in the eukaryotic cell. Eukaryot. Cell 1, 11-21.

  90. Barbet, N. C., Schneider, U., Helliwell, S. B., Stansfield, I., Tuite, M. F. and Hall, M. N. (1996). TOR controls translation initiation and early G1 progression in yeast. Mol. Biol. Cell 7, 25-42.

  91. Beck, T. and Hall, M. N. (1999). The TOR signaling pathway controls the nuclear localization of nutrient regulated transcription factors. Nature 402, 689-692.

  92. Noda, T. and Ohsumi, Y. (1998). Tor, a phosphatidylinositol kinase homologue, controls autophagy in yeast. J. Biol. Chem. 273, 3963-3966.