Systematic Identification of the Genes Affecting Glycogen Storage in the Yeast Saccharomyces cerevisiae

Implication of the Vacuole as a Determinant of Glycogen Level*,S

Wayne A. Wilson, Zhong Wang and Peter J. Roach{ddagger}

{ddagger} Department of Biochemistry and Molecular Biology and Center for Diabetes Research, Indiana University School of Medicine, Indianapolis, Indiana 46202


    ABSTRACT
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
At the onset of nutrient limitation, the yeast Saccharomyces cerevisiae synthesizes glycogen to serve as a carbon and energy reserve. We undertook a systematic survey for the genes that affect glycogen accumulation by taking advantage of the strain deletion set generated by the Saccharomyces Genome Deletion Project. The strain collection analyzed contained some 4600 diploid homozygous null deletants, representing ~88% of all viable haploid disruptants. We identified 324 strains with low and 242 with elevated glycogen stores, accounting for 12.4% of the genes analyzed. The screen was validated by the identification of many of the genes known already to influence glycogen accumulation. Many of the mutants could be placed into coherent families. For example, 195 or 60% of the hypoaccumulators carry mutations linked to respiratory function, a class of mutants well known to be defective in glycogen storage. The second largest group consists of ~60 genes involved in vesicular trafficking and vacuolar function, including genes encoding 13 of 17 proteins involved in the structure or assembly of the vacuolar ATPase. These data are consistent with our recent findings that the process of autophagy has a significant impact on glycogen storage (Wang, Z., Wilson, W. A., Fujino, M. A., and Roach, P. J. (2001) Antagonistic controls of autophagy and glycogen accumulation by Snf1p, the yeast homolog of AMP-activated protein kinase, and the cyclin-dependent kinase Pho85p. Mol. Cell. Biol. 21, 5742–5752). Autophagy delivers glycogen to the vacuole, and we propose that the impaired vacuolar function associated with ATPase mutants (vma10 or vma22) results in reduced degradation and subsequent hyperaccumulation of glycogen.


One attractive feature of the budding yeast Saccharomyces cerevisiae as an experimental organism has been its ease of genetic manipulation, including the ability to perform genetic screens by which yeast strains with a phenotype of interest can be recognized. Experimentally, the next phase involves identification of the gene(s) responsible, a task that can be time-consuming and sometimes non-trivial. The S. cerevisiae genome contains some 6,200 open reading frames (ORFs),1 and these have been disrupted systematically in the Saccharomyces genome deletion project (1). The availability of the resulting set of deletion strains, each carrying a deletion of one specific ORF, permits a totally different type of screen or survey for genes linked to a particular phenotype. Most importantly, the survey is systematic. The first report of using a partial strain deletion set analyzed growth rates in rich and minimal medium (1). Recently, an analysis of the rapamycin sensitivity of 2,216 haploid disruptants was reported (2), and a systematic study of fluid phase endocytosis using around 700 strains generated by the European Functional Analysis Network has been conducted (3). The first truly genome-wide screen of a defined null mutant collection was reported by Ni and Snyder (4) who analyzed over 4000 strains in a study of polarized growth. Information is also available about viability, based on analysis of the deletion strain set, and is available from the Saccharomyces Genome Database (5) (genome-www.stanford.edu/Saccharomyces) . We report here the use of a specific metabolic end point, the ability to store glycogen, as the basis for a screen of ~4600 homozygous diploid mutants to identify genes that affect glycogen accumulation.

Glycogen serves as a reserve of glucose. Its accumulation is initiated under conditions of nutrient limitation, such as the approach to stationary phase in liquid culture. Limitation for carbon, nitrogen, phosphorous, or sulfur all act as triggers for increased glycogen synthesis (6). Our laboratory has been interested in glycogen as an example, in mammals and yeast, of a compound whose synthesis and utilization is under complex and intricate controls linked to the intracellular energy state, as well as the nutritional status of the environment (see Ref. 7 for a review). Synthesis of glycogen requires the activities of glycogenin, a self-glucosylating initiator protein (encoded by GLG1 and GLG2; see Ref.8), glycogen synthase (GSY1 and GSY2; see Ref.9), which catalyzes bulk synthesis, and the branching enzyme (GLC3; see Ref.10), which introduces the branches characteristic of the mature polysaccharide. Glycogen breakdown requires glycogen phosphorylase (Gph1p; see Ref. 11) and debranching enzyme (Gdb1p; see Refs. 12 and 13) or, under certain conditions, glucoamylase (Sga1p; see Ref. 14). However, the enzymes of glycogen metabolism are under a variety of transcriptional and post-translational controls, and so genes encoding a number of other proteins affect glycogen accumulation (for a review, see Ref. 7). For example, the cyclic AMP pathway controls both gene expression and phosphorylation of key proteins (1517). Starvation, as sensed by the Tor pathway, stimulates glycogen accumulation (18). Signaling through the Snf1p and Pho85p protein kinases has antagonistic effects on glycogen storage (19). The genes mentioned above have been identified by a combination of conventional biochemical and genetic approaches, and there is no guarantee that all relevant genes have been found or any indication as to how many genes affect glycogen storage. The systematic survey described in this work indicated that 566 of ~4600 strains from the homozygous diploid release of the deletion library had glycogen levels that differed from wild type. Of these genes, a surprising fraction, about 10%, had functions related to vesicular trafficking or vacuolar function.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Strains and Media—
The homozygous diploid deletion series (BY4743 strain background) was purchased from Research Genetics. The yeast deletion series comprises a set of mutants where each open reading frame has been disrupted from start to stop codon, and a kanMX marker cassette (conferring resistance to the antibiotic G418) has been inserted. The library of deletions was supplied frozen in 96-well microtiter plates. Each well contained 200 µl of YPD (2% peptone, 2% glucose, 1% yeast extract) medium supplemented with G418 (200 µg/ml) and 15% glycerol. For keying, contamination checks and orientation purposes, two wells per plate contained only medium. One empty well defines the bottom left hand corner, and the other serves as an identifier. Complete details regarding series construction and availability can be found at sequence-www.stanford.edu/group/yeast_deletion_project/deletions3.html and www.resgen.com/products/YEASTD.php3. The deletion set contained 4,639 different strains, of which 82 were flagged by Research Genetics as quality-control failures. Excluding these, the collection is of 4,557 strains that correspond to 88% of all possible viable mutants. Additionally, we constructed strains derived from BY4741 and BY4742, which are also available from Research Genetics and are the MATa and MAT{alpha} parents, respectively, of BY4743. A polymerase chain reaction strategy (20) was used to generate apg1::URA3 in BY4741 and vma10::LEU2 in BY4742. These strains were mated, and tetrad dissection was performed to isolate apg1::URA3 vma10::LEU2 double mutants of both mating types. The genotypes of the strains isolated were MATahis3 leu2 met15 lys2 apg1::URA3 vma10::LEU2 and MAT{alpha} his3 leu2 met15 lys2 apg1::URA3 vma10::LEU2. These MATa and MAT{alpha} strains were then crossed to generate a diploid strain.

Growth Conditions—
For analysis, the master plates were thawed, and the cells were resuspended by pipetting up and down using a 12-channel automatic pipette. A 96-pin microplate replicator was used to transfer an aliquot from each well to a fresh microtiter plate containing 200 µl of appropriate medium per well. To assess glycogen accumulation, SC medium (0.67% yeast nitrogen base, 0.079% Complete Supplement Mix (Bio101 Inc.), 2% glucose) was used. Growth on YPG medium (2% peptone, 1% yeast extract, 3% glycerol) was used to assess respiratory competence. For glycogen accumulation studies, plates were prepared in duplicate. The parental wild-type strain (BY4743) was inoculated into well H1 on each plate to serve as an internal control. Plates were incubated at 30°C for 48 h without shaking.

Genomic Survey for Aberrant Glycogen Accumulation—
The 54 microtiter plates that constitute the deletion series were analyzed for both glycogen accumulation and ability to grow using glycerol as a carbon source. Preliminary studies with strains known to over- or underaccumulate glycogen were used to optimize the growth and staining conditions. Growth of cells in the microplate format was slower than in standard liquid culture in shaking flasks, likely because of a lack of aeration. Thus, to obtain enough cells for analysis of glycogen content by iodine staining, we grew cultures for 48 h. Cells were resuspended by pipetting up and down using a 12-channel automatic pipette. The cell suspensions were then transferred to a 96-well vacuum manifold (Bio-Dot microfiltration unit; Bio-Rad) and filtered onto a nitrocellulose membrane. The membrane was removed and stained for 2 min by exposure to iodine vapor, a time that was found to be suitable for uniform and reproducible results. Each membrane was photographed using a digital camera to create a permanent record, and cultures that stained either more or less intensely than wild type were scored (Fig. 1). Images were downloaded from the camera to an Apple Macintosh G4 microcomputer and processed with Adobe PhotoShop LE. In each case, duplicate plates were grown, harvested, and stained to control for any variability in staining intensity, and only wells that gave the same result on both plates were scored in the final tally. Growth in glycerol was assessed by visual inspection of plates and comparison to the congenic wild-type strain.



View larger version (118K):
[in this window]
[in a new window]
 
FIG. 1. Screen for glycogen accumulation using a 96-well microtiter plate format. Cells were cultured in 96-well microtiter plates and collected by filtration onto nitrocellulose membranes. The nitrocellulose membranes were stained by exposure to iodine vapor to assess glycogen accumulation, and wells staining either more (ovals) or less (squares) intensely than the well containing the wild-type strain (arrowheads) were scored. The empty keying well, which identifies the plate number in the library, is marked with a circle.

 
Mutant Identification—
The deleted gene responsible for the aberrant glycogen phenotype in each case was identified by reference to a spreadsheet compiled from data available at www.resgen.com/products/YEASTD.php3 where a text file detailing the ORF deletion present in each well of each microtiter plate can be found. Information from the Saccharomyces Genome Database (5) (genome-www.stanford.edu/Saccharomyces) and from the YPD database (21, 22) (www.proteome.com/databases/index.html) was used to pair ORF numbers with gene names and functional properties where known.

Enzymatic Determination of Glycogen—
For quantitative determination of glycogen levels, cells were grown for 24 h to early stationary phase (~1 x 108 cells/ml) in 10 ml of SC medium at 30°C. The cell density was checked by counting using a hemocytometer. An aliquot (1 ml) of culture was removed, and the cells were collected by centrifugation (14,000 x g, 1 min, 4°C). The supernatant was aspirated, and the cell pellet was immediately frozen on dry ice and stored at -80°C until use. The culture was maintained for another 24 h, and a second sample (48 h time point) was taken as described. Cell pellets were thawed by addition of 200 µl of 20% (w/v) KOH and boiled in a water bath for 1 h. The solution was neutralized by addition of HCl and processed as described previously (23). The data shown are the mean of at least two independent determinations performed in duplicate.

Determination of Glycogen Synthase Activity—
Glycogen synthase was assayed in extracts prepared from yeast cells by lysis with glass beads using the method of Thomas et al. (24) as described previously (25). Phosphorylation of glycogen synthase converts the enzyme into a less active form that requires the presence of the allosteric activator glucose 6-phosphate to elicit full activity. Thus, the ratio of activity without and with glucose 6-phosphate (-/+ glucose 6-phosphate activity ratio) is an index of the phosphorylation state of the enzyme with high values indicating that dephosphorylated and active glycogen synthase predominates.

Induction of Autophagy, Microscopy, and Image Capture—
To induce autophagy, cells were grown to late logarithmic phase in YPD medium and then starved for nitrogen as described by Takeshige et al. (26) either in the presence or absence of 1 mM phenylmethylsulfonyl fluoride (PMSF). For analysis of vma10, vma22, and apg1 vma10 mutants, cells were grown to logarithmic phase or to saturation in SC medium. Cells were examined under a Nikon Microphot-FXA microscope, equipped with Nomarski optics, using a x60 oil-immersion objective at a magnification of x1200. Images were captured using a Pulnix TM-745 digital camera and the public domain NIH Image 1.62 software (developed at the United States National Institutes of Health and available on the Internet at rsb.info.nih.gov/nih-image/) running on a Macintosh G4 microcomputer. Movies were prepared and edited using QuickTime Pro 5.0.1 software from Apple Computer.


    RESULTS
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
Overview—
The strain deletion set was surveyed, in duplicate, for glycogen accumulation using iodine staining of cells harvested after growth in microtiter plates and filtration through nitrocellulose filters, as described under "Materials and Methods" (see Fig. 1). In addition, the ability to grow on glycerol was monitored. The rationale for including this assay was the well established observation that cells unable to grow using non-fermentable carbon sources cannot store glycogen, usually because of respiratory defects in what are termed petite mutants (27, 28). In this way, we could place mutants in this category, even if the corresponding gene had no known function.

From the 4557 strains examined, we recovered 242 with high glycogen and 324 with low glycogen so that a total of 12.4% of the gene deletions influenced glycogen. We attempted to cluster the mutants into logical families (Table I). In some cases the families are clearly meaningful, with the presence of multiple genes of related function reinforcing the validity of their identification in the survey. In other instances, we grouped proteins according to biochemical function, such as "protein kinases," where common biological functions were not obvious. A large group of 100 mutants, representing almost 18% of the genes identified, corresponded to ORFs of unknown function, about which little or nothing is known. Thirty-five genes involved in general transcription or RNA processing were identified as were 18 that had to do with DNA structure or maintenance of structural integrity. Our thought is that the mutations in these strains are likely to be only very indirectly linked to glycogen metabolism. Glycogen storage is in part a response to stress, which may be an aspect of the phenotype of these groups of mutants. Such could also be true for the genes encoding small ribosomal subunits, all 12 of which were associated with glycogen hyperaccumulation. Because of the large number of mutants identified, it is not possible to give a thorough description of each (refer to the Supplemental Material for a table containing all of the strains identified and all of our annotations).


View this table:
[in this window]
[in a new window]
 
TABLE I Gene families

The mutants that were identified in the screen were assigned to various different families. For details of all family members and annotations, see Supplemental Material.

 
A number of mutants were chosen for additional analyses. The selection was made to represent several of the gene families of Table I, together with some subjective choices (Table II). The further analyses included assessment of glycogen levels as judged by iodine staining of colonies on either YPD or SC plates, enzymatic determination of glycogen in liquid cultures, and measurement of glycogen synthase activity. These secondary measures of glycogen usually reproduced the original screening result, especially noting that yeast accumulate generally much less glycogen when grown in rich as opposed to synthetic medium. Thus, a number of mutants scored as high glycogen on the original survey with synthetic medium were wild type when grown on YPD plates. A notable exception is provided by some of the vacuolar and vesicular trafficking mutants (see below). Of the mutants selected for further study, few exhibited very great changes in total glycogen synthase activity. Glycogen synthase phosphorylation, and activation state, can be monitored via the -/+ glucose 6-phosphate activity ratio, lower values correlating with greater phosphorylation (see "Materials and Methods"). Deletion of genes implicated in direct phosphorylation of glycogen synthase, such as PHO85, leads to elevated activity ratio. Of the genes tested, only one, a WD-40 repeat protein of unknown function, had a significantly elevated activity ratio, making it a candidate to be a constituent of a glycogen synthase kinase.


View this table:
[in this window]
[in a new window]
 
TABLE II Further analysis of selected deletion mutants isolated in the screen

The amount of glycogen and the glycogen synthase activities at 24 and 48 h are normalized to the values obtained with wild-type cells at the same time points. For the qualitative determinations of glycogen on YPD and SC plates, wild-type cells were scored as (+). The activity ratio of glycogen synthase in the wild type strain at the 24-h time point was 0.26 ± 0.02 (n = 6). ND indicates no determination could be made.

 
Genes Implicated Previously in Glycogen Storage—
A number of the genes that were isolated in the screen have been implicated previously in the control of glycogen accumulation (Table III), which we view as a validation of the screening methodology. In an earlier genetic screen for aberrant glycogen accumulation, Cannon et al. (29) had characterized eight so-called glc mutants and identified the corresponding genes. Five GLC genes were found in our screen (Table III). Of the three remaining glc mutants, two were not represented in the deletion series (GLC7 is essential, and GLC5/IRA1 was missing) and one, glc1/ras2, was present but had wild-type levels of glycogen in this genetic background under the growth conditions used.


View this table:
[in this window]
[in a new window]
 
TABLE III Genes isolated by the screen and that had been implicated previously in glycogen accumulation

 
Genes encoding four enzymes involved directly in glycogen metabolism were identified. Strains defective for the degradative enzymes, glycogen phosphorylase, Gph1p, and debranching enzyme, Gdb1p, were hyperaccumulators. The strain lacking GSY2, which encodes the predominant isoform of glycogen synthase, had low glycogen storage as did the strain lacking the glycogen branching enzyme encoded by GLC3. The PFK1 and TPS1 genes encode components of phosphofructokinase and the trehalose synthase complex, respectively. Both pfk1 and tps1 mutations result in elevated intracellular glucose 6-phosphate (30, 31). It is probable that the increased glucose 6-phosphate levels bypass the phosphorylation control of glycogen synthase resulting in deregulated, hyperactive glycogen synthase (see for example Ref. 32).

Four putative or actual targeting/regulatory subunits for the type 1 protein phosphatase Glc7p were identified. Of these, the best characterized is Gac1p, a protein that targets Glc7p to dephosphorylate and activate glycogen synthase (33, 34). The GLC8 gene encodes a protein related in sequence to mammalian I-2, a constituent of a cytosolic form of type I protein phosphatase (29). The putative Glc7p targeting protein Pig2p was isolated first in a two-hybrid screen for proteins that interacted with yeast glycogen synthase (35). Previous studies had not revealed any substantial role for Pig2p in glycogen accumulation, but this may reflect differences in the strain backgrounds used. The fourth protein phosphatase non-catalytic subunit isolated was Reg1p, which Tu and Carlson (36) had shown to physically associate with Glc7p and which had previously been implicated in glucose repression of gene expression (37). Our laboratory had also isolated reg1 mutants as suppressors of the glc7–1 glycogen accumulation defect (38). The role of SNF1 in controlling glycogen synthase is well established and is thought to result both from a negative control of the phosphorylation of glycogen synthase (17) and a positive control of autophagy, which normally preserves glycogen (39). Strains lacking snf1 are therefore defective in glycogen storage. Likewise, the PKA pathway has long been known to play a role in glycogen accumulation, with PKA negatively regulating transcription of GSY2 and other enzymes of glycogen metabolism (1517). Rim15p is proposed to function downstream of PKA, as a regulator of entry into stationary phase, and rim15 mutants had been reported to have low levels of glycogen and the other storage carbohydrate, trehalose (40).

Genes Required for Respiratory Growth—
Of the 324 mutants with reduced glycogen stores, 60% of the disruptions were in genes that were either known to be required for respiratory growth (e.g. those encoding components of the mitochondrion) or were found in the screen to be unable to grow with glycerol as a carbon source. Many mutants that fail to utilize non-fermentable carbon sources actually synthesize glycogen (28), but this glycogen is used to fuel growth upon exhaustion of glucose, because, unlike wild-type cells, respiratory mutants cannot oxidize the ethanol produced in the initial growth phase. It is noteworthy that we found 11 mutants that were unable to grow using glycerol but had increased glycogen. It is possible that these mutants are defective in glycogen breakdown and therefore might be worthy of future study.

Other Protein Kinases and Phosphatases—
Protein kinases linked to glycogen accumulation have the potential to be involved in regulatory pathways, like Snf1p or Rim15p, or the direct phosphorylation of metabolic enzymes. Thirteen conventional protein Ser/Thr kinases were recovered in the screen, six with reduced glycogen (including the snf1 and rim15 mutants discussed above) and seven with increased glycogen. Of those kinase mutants with elevated glycogen, none was associated with a significant increase in the glycogen synthase activity ratio, as would be predicted for a glycogen synthase kinase. An additional protein kinase of the phosphatidylinositol 3-/phosphatidylinositol 4-kinase family was identified, namely Tor1p. The strain with a tor1 mutation hyperaccumulated glycogen, consistent with the known effect of rapamycin treatment, to cause increased glycogen synthesis (18).

In addition to the Glc7p-associated subunits described above, several other protein phosphatase subunits were identified. The deletion mutant of RTS1, which encodes a homolog of the mammalian B' regulatory subunit of type 2A protein phosphatases (41), was found to have greatly increased glycogen. Although type 2A phosphatases had been implicated previously in glycogen storage (4244), this is the first indication of a role for RTS1 in the process. Deletion of YCR079W, which encodes a protein phosphatase of the type 2C family (45), also caused increased glycogen. Ycr079p has been little studied, and its cellular role is unknown, although deletion of YCR079W has been reported to increase the sensitivity of the cell to caffeine, which could indicate a link to the cyclic AMP pathway (46).

Genes Implicated in Vacuole Function or Vesicle Transport—
The second largest family of genes identified has to do with vesicular trafficking or vacuolar function (Table IV). Approximately 10% of the mutants with aberrant glycogen levels fall into this category. This result would have been quite perplexing but for recent work in our laboratory indicating that yeast defective for autophagy are unable to maintain their glycogen stores even though the polysaccharide is synthesized normally (39). Autophagy, in response to starvation signals, is a process whereby cytosol and organelles are engulfed to become autophagosomes that are delivered to the vacuole for the recycling of constituents (reviewed in Refs. 47 and 48). Therefore, a variety of defects in the formation or delivery of autophagosomes or in vacuolar function could be linked with glycogen storage. Correspondingly, genes with a wide variety of specific roles were detected, and it is not simple to rationalize in detail the results for every individual gene.


View this table:
[in this window]
[in a new window]
 
TABLE IV Genes identified in the screen and that have roles in vesicle trafficking or vacuole function

 
A clearly defined subfamily in this category involves components of the vacuolar H+-ATPase that is responsible for acidification of the vacuole (reviewed in Ref. 49). This complex comprises some 13 different gene products plus at least four factors that participate in its assembly. Of these 17 different genes, deletion mutants in 13 were recovered in the screen, 11 with high and two with low glycogen. Deletion mutants in the remaining genes (VMA2, VMA7, VMA8, and VMA11) were present in the deletion series but were scored as wild type. The two mutants that had low glycogen (vma5 and vma16) were the only two of the 17 that failed to grow with glycerol as a carbon source, suggesting a respiratory defect that could explain the low glycogen levels. It has been reported that a characteristic of the vma class of mutants is their inability to grow using non-fermentable carbon sources (50), but presumably there must be strain to strain variation, because most of the vma mutants studied here did grow on glycerol. These ATPase genes represent an obvious functional cluster, making a robust link with glycogen storage.

Although the exact mechanistic link between autophagy and glycogen storage is not completely understood, our previous studies (39) led us to two conclusions. First, autophagy provides a source of intermediates and energy in stationary phase, and its absence may promote premature utilization of cytosolic glycogen. Second, autophagy delivers glycogen to the vacuole where it is actually protected from this cytosolic utilization and is degraded by vacuolar enzymes such as Sga1p. Ohsumi and colleagues (26) developed an assay for autophagy in which yeast growing in rich medium are transferred to a starvation medium, totally lacking nitrogen, and also containing PMSF. The rationale is that starvation will induce autophagy, and the PMSF will inhibit vacuolar proteases thus preventing degradation of autophagosomes delivered to the vacuole. Wild-type cells transferred to starvation medium in the presence of PMSF have full, active vacuoles compared with cells not exposed to PMSF (Fig. 2). This difference is hardly evident in the still shot but is very obvious if the cells are viewed over time when undegraded autophagic bodies can be seen moving rapidly within the confines of the vacuole (see supplemental video for Fig. 2). APG1 encodes a protein kinase that is necessary for autophagy, and analysis of apg1 mutants using the starvation assay reveals empty vacuoles, evident even in the still shot. When monitored over time, there is no motion within the vacuole, but particles can be seen moving outside of the vacuole (see supplemental video for Fig. 2).



View larger version (121K):
[in this window]
[in a new window]
 
FIG. 2. Comparison of the vacuoles of wild-type and apg1 mutant cells subjected to nitrogen starvation. The wild-type strain, BY4743 (panels A and B), and an apg1 mutant (panels C and D) in the BY4743 background were grown to late logarithmic phase in YPD medium. The cells were then transferred to nitrogen starvation medium (26), and either the protease inhibitor PMSF dissolved in isopropanol (panels B and D) or isopropanol (panels A and C) was added. Incubation was continued for 4 h, and the cells were examined under a microscope equipped with Nomarski optics. The scale bar represents 5 µm.

 
We monitored by microscopy two of the mutants defective in the vacuolar ATPase, vma10 and vma22, during growth on SC medium in liquid culture. At the stage of exponential growth, the vma10 and vma22 mutants had more clearly defined vacuoles than wild-type cells (Fig. 3), but typically the level of movement within the vacuole was similar (see video for Fig. 3). A small fraction, 5%, of the mutants had large and full vacuoles, however (not shown). After 24 h, most of the mutants cells, ~60%, had large well defined vacuoles inside of which an accumulation of particles can be seen in rapid motion (see supplemental video for Fig. 3). In fact, these vacuoles resembled the wild-type cells subjected to the autophagy assay in the presence of PMSF (Fig. 2). Autophagy is presumably induced normally upon entry to stationary phase (26, 39), but degradation of autophagic bodies is impaired, because the vacuole cannot acidify because of a defective ATPase (51). A double mutant that lacked both apg1 and vma10 (Fig. 4) resembled the apg1 mutant rather than the vma10 mutant, as has been reported by Nakamura et al. (51). Thus, with autophagy impaired, the inability of the vma mutants to degrade autophagic bodies is not in evidence. Interestingly, analysis of the glycogen content revealed the same epistatic relationship, with glycogen accumulation in the double mutant resembling that of apg1 mutants rather than vma10 mutants (Table V). Therefore, the elevated glycogen storage associated with vma10 and vma22 mutants is most likely explained by hyperaccumulation of vacuolar glycogen.



View larger version (113K):
[in this window]
[in a new window]
 
FIG. 3. Accumulation of autophagic bodies in the vacuole of vma10 and vma22 mutants. The wild-type strain BY4743 (panels A and B) and vma10 (panels C and D) or vma22 (panels E and F) deletion mutants in the BY4743 background were grown in SC medium. Aliquots were removed in the logarithmic phase of growth (panels A, C, and E) and again when the cultures were saturated (panels B, D, and F). The cells were examined under a microscope equipped with Nomarski optics. The scale bar represents 5 µm.

 


View larger version (45K):
[in this window]
[in a new window]
 
FIG. 4. The vacuoles of apg1 vma10 mutants fail to accumulate autophagic bodies. A homozygous diploid strain in which both the APG1 and VMA10 genes were deleted was grown in SC medium. Aliquots were removed during the logarithmic phase of growth (panel A) and at saturation (panel B). The cells were examined under a microscope equipped with Nomarski optics. The scale bar represents 5 µm.

 

View this table:
[in this window]
[in a new window]
 
TABLE V Glycogen accumulation by various yeast strains

The results shown are the mean ± S.E. for three independent determinations.

 

    DISCUSSION
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
The availability of the gene deletion strain collection is symbolic of an era in which not only information but also reagents are available on a genomic scale. In the present case, the entire strain collection occupies only 54 microtiter plates so that, even without investment in robotic instrumentation, many types of phenotype scan can be performed with relative ease (for example, see the work of Ni and Synder (4) or Wiederkehr et al. (3)). It is of interest to compare the approach described in this study with conventional genetic screens. Advantages of the genomic scan are severalfold. One is simplicity. Once a phenotype is scored, gene identification is by reference to a database. Nonetheless, if an important conclusion would be based on a single member of the deletion set, it would obviously be prudent to confirm that the phenotype can be reproduced by an independent gene knock-out. This is necessary as part of ongoing quality control as the deletion set becomes used by the yeast research community. Another advantage is that the screen is systematic and does not depend fundamentally on chance; coverage of the genome with mutations is defined. As has been stated by Ni and Snyder (4), there are several disadvantages to this new mode of screening. First, it is limited to studying loss of function mutations. Second, in the simplest approach, essential genes are obviously not included, and for the homozygous diploid set, nor are genes required for mating. Also, because the duplication of the yeast genome (52) has resulted in some pairs of genes of very high sequence identity, no attempt was made to target both genes of the pair, and so not all are represented. A related point, which applies equally to conventional screens, is that any genes that have redundant functions are unlikely to be detected. Finally, some gene deletions are simply missing from the series at this time. Nonetheless, the collection currently includes ~88% of viable mutants and is a valuable resource for this new approach to genetic screening.

As described under "Results," many genes known already to affect glycogen storage were identified, providing strong validation of the screen. Some genes that were expected did not surface in the screen. One of the most notable was PHO85, which we know from our previous work is involved in the direct regulation of glycogen synthesis (19, 53). In the genetic background of the strain that we use in most of our work, deletion of PHO85 causes a substantial hyperaccumulation of glycogen. However, in the present screen, the pho85 mutant was scored wild type for glycogen. The explanation is the difference in strain backgrounds, because we have found that, in some strains, loss of pho85 causes activation of glycogen synthase but no increase in glycogen accumulation.2 It is likely that PHO85 has other connections to the control of glycogen storage.3 Lee et al. (54) also noted that in some strains PHO85 deletion does not result in glycogen hyperaccumulation.

With our level of knowledge of yeast glycogen metabolism and its control prior to our screen, one would have had difficulty justifying more than 30 genes that affect glycogen storage if respiratory or mitochondrial mutants were excluded. These would have been the genes responsible for glycogen and related metabolism, direct regulators like protein kinases and phosphatases, or genes involved in known regulatory pathways such as those involving cyclic AMP, SNF1 and glucose repression, PHO85, and the Tor pathway. Therefore, the first notable feature of our results is the large number of genes identified in the screen, which is 362 even if the respiratory mutants are excluded. Of these, 100 have unknown function. It is likely that, as biological functions are assigned, these genes will for the most part fit into existing categories, and we may have already defined most of the biologically important families. Given the large number of novel mutations identified here as affecting glycogen storage compared with the small number known previously to influence this process, it is of interest to ask whether the previously identified mutants were those that conferred the strongest glycogen accumulation phenotype. To a certain extent this may be true, because reg1 and pfk1 mutants, which had previously been shown to have a high level of glycogen, were among the mutants that stained most intensely with iodine in the present screen. However, several of the newly identified mutants, such as rts1, vid21, and ado1, stained just as strongly as reg1.

The second striking aspect of the results was that, based on the 35 strains where we followed up with glycogen synthase assays, a large fraction of the changes in glycogen storage did not correlate with the glycogen synthase activity ratio, which indexes its phosphorylation and activation state. Part of the explanation is our current appreciation that not only synthesis but also the ability to retain glycogen is a key determinant of glycogen storage (39). This conclusion, based on recent analysis of SNF1 and its relationship to autophagy, is reinforced by the results of this study.

The most interesting gene category to emerge is that defined by the genes involved in vesicular trafficking and vacuolar function. This is a family where there is a reasonably good connection with the process of autophagy, which we have recently linked to glycogen storage (39). We have proposed that defective autophagy leads to premature depletion of glycogen stores. The mechanism is not fully understood, but at least two points can be made. First, in the absence of autophagy, the cell has lost an important source of intermediary metabolites and may be forced to use other reserves, like glycogen, at an earlier stage than normal. Second, there is a vacuolar glycogen pool (26), and we hypothesized that vacuolar glycogen is protected from cytosolic degradation until late into stationary phase (39). Thus, many of the genes identified with low glycogen in this group likely reflect some defect in autophagy.

Within the vesicular/vacuolar category, the H+-ATPase genes are especially interesting. First, the fact that such a large fraction of the ATPase genes was netted in the screen makes the case very strongly that this vacuolar function is linked to glycogen levels. It has been reported that elimination of any one of the subunits cripples ATPase function, consistent with our identifying genes coding for several ATPase subunits. Without its ATPase, the vacuole cannot acidify and function normally. In particular, the degradative enzymes of the vacuole, including the vacuolar glycosidase Sga1p, will not be fully active (51, 55). Impaired vacuolar function removes a source of intermediates and energy in stationary phase, and one might have predicted a premature utilization of glycogen from the cytosol, much as is seen for those mutants defective in autophagy (39). The fact that defective ATPase correlated with a hyperaccumulation of glycogen is thus best explained by a stabilization of the vacuolar pool. This proposal is also consistent with the observation that the elevated glycogen in a vma10 mutant is eliminated if autophagy is also impaired by mutation of APG1. In conclusion, the study suggests a significant role for the vacuole in the maintenance of yeast glycogen stores.


    ACKNOWLEDGMENTS
 
We thank Drs. Mark Goebl and Ron Wek for many helpful discussions.


    FOOTNOTES
 
Received, October 11, 2001, and in revised form, January 30, 2002.

Published, MCP Papers in Press, February 18, 2002, DOI 10.1074/mcp.M100024-MCP200

1 The abbreviations used are: ORF, open reading frame; PMSF, phenylmethylsulfonyl fluoride. Back

2 M. A. Fujino and P. J. Roach, unpublished observations. Back

3 W. A. Wilson, Z. Wang, M. A. Fujino, and P. J. Roach, unpublished observations. Back

* This work was supported by National Institute of Health Grant DK42576 and the Indiana University Diabetes Research and Training Center (DK20542). Back

S The on-line version of this article (available at http://www.mcponline.org) contains Supplemental Material. Back

{ddagger} To whom correspondence should be addressed: Dept. of Biochemistry and Molecular Biology, Indiana University School of Medicine, 635 Barnhill Dr., Indianapolis, IN 46202-5122. Tel.: 317-274-1582; Fax: 317-274-4686; E-mail: proach{at}iupui.edu.


    REFERENCES
 TOP
 ABSTRACT
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 

  1. Winzeler, E. A., Shoemaker, D. D., Astromoff, A., Liang, H., Anderson, K., Andre, B., Bangham, R., Benito, R., Boeke, J. D., Bussey, H., Chu, A. M., Connelly, C., Davis, K., Dietrich, F., Dow, S. W., El Bakkoury, M., Foury, F., Friend, S. H., Gentalen, E., Giaever, G., Hegemann, J. H., Jones, T., Laub, M., Liao, H., and Davis, R. W. (1999) Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science 285, 901 –906[Abstract/Free Full Text]

  2. Chan, T. F., Carvalho, J., Riles, L., and Zheng, X. F. (2000) A chemical genomics approach toward understanding the global functions of the target of rapamycin protein (TOR). Proc. Natl. Acad. Sci. U. S. A. 97, 13227 –13232[Abstract/Free Full Text]

  3. Wiederkehr, A., Meier, K. D., and Riezman, H. (2001) Identification and characterization of Saccharomyces cerevisiae mutants defective in fluid-phase endocytosis. Yeast 18, 759 –773[CrossRef][Medline]

  4. Ni, L., and Snyder, M. (2001) A genomic study of the bipolar bud site selection pattern in Saccharomyces cerevisiae. Mol. Biol. Cell 12, 2147 –2170[Abstract/Free Full Text]

  5. Cherry, J. M., Adler, C., Ball, C., Chervitz, S. A., Dwight, S. S., Hester, E. T., Jia, Y., Juvik, G., Roe, T., Schroeder, M., Weng, S., and Botstein, D. (1998) SGD: Saccharomyces Genome Database. Nucleic Acids Res. 26, 73 –79[Abstract/Free Full Text]

  6. Lillie, S. H., and Pringle, J. R. (1980) Reserve carbohydrate metabolism in Saccharomyces cerevisiae: responses to nutrient limitation. J. Bacteriol. 143, 1384 –1394[Medline]

  7. Francois, J., and Parrou, J. L. (2001) Reserve carbohydrates metabolism in the yeast Saccharomyces cerevisiae. FEMS Microbiol. Rev. 25, 125 –145[CrossRef][Medline]

  8. Cheng, C., Mu, J., Farkas, I., Huang, D., Goebl, M. G., and Roach, P. J. (1995) Requirement of the self-glucosylating initiator proteins Glg1p and Glg2p for glycogen accumulation in Saccharomyces cerevisiae. Mol. Cell. Biol. 15, 6632 –6640[Abstract]

  9. Farkas, I., Hardy, T. A., Goebl, M. G., and Roach, P. J. (1991) Two glycogen synthase isoforms in Saccharomyces cerevisiae are coded by distinct genes that are differentially controlled. J. Biol. Chem. 266, 15602 –15607[Abstract/Free Full Text]

  10. Rowen, D. W., Meinke, M., and LaPorte, D. C. (1992) GLC3 and GHA1 of Saccharomyces cerevisiae are allelic and encode the glycogen branching enzyme. Mol. Cell. Biol. 12, 22 –29[Abstract]

  11. Hwang, P. K., and Fletterick, R. J. (1986) Convergent and divergent evolution of regulatory sites in eukaryotic phosphorylases. Nature 324, 80 –84[Medline]

  12. Nakayama, A., Yamamoto, K., and Tabata, S. (2000) High expression of glycogen-debranching enzyme in Escherichia coli and its competent purification method. Protein Expr. Purif. 19, 298 –303[CrossRef][Medline]

  13. Teste, M. A., Enjalbert, B., Parrou, J. L., and Francois, J. M. (2000) The Saccharomyces cerevisiae YPR184w gene encodes the glycogen debranching enzyme. FEMS Microbiol. Lett. 193, 105 –110[CrossRef][Medline]

  14. Yamashita, I., and Fukui, S. (1985) Transcriptional control of the sporulation-specific glucoamylase gene in the yeast Saccharomyces cerevisiae. Mol. Cell. Biol. 5, 3069 –3073[Medline]

  15. Smith, A., Ward, M. P., and Garrett, S. (1998) Yeast PKA represses Msn2p/Msn4p-dependent gene expression to regulate growth, stress response and glycogen accumulation. EMBO J. 17, 3556 –3564[Abstract/Free Full Text]

  16. Gorner, W., Durchschlag, E., Martinez-Pastor, M. T., Estruch, F., Ammerer, G., Hamilton, B., Ruis, H., and Schuller, C. (1998) Nuclear localization of the C2H2 zinc finger protein Msn2p is regulated by stress and protein kinase A activity. Genes Dev. 12, 586 –597[Abstract/Free Full Text]

  17. Hardy, T. A., Huang, D., and Roach, P. J. (1994) Interactions between cAMP-dependent and SNF1 protein kinases in the control of glycogen accumulation in Saccharomyces cerevisiae. J. Biol. Chem. 269, 27907 –27913[Abstract/Free Full Text]

  18. 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[Abstract]

  19. Huang, D., Farkas, I., and Roach, P. J. (1996) Pho85p, a cyclin-dependent protein kinase, and the Snf1p protein kinase act antagonistically to control glycogen accumulation in Saccharomyces cerevisiae. Mol. Cell. Biol. 16, 4357 –4365[Abstract]

  20. Wach, A. (1996) PCR-synthesis of marker cassettes with long flanking homology regions for gene disruptions in S. cerevisiae. Yeast 12, 259 –265[CrossRef][Medline]

  21. Costanzo, M. C., Hogan, J. D., Cusick, M. E., Davis, B. P., Fancher, A. M., Hodges, P. E., Kondu, P., Lengieza, C., Lew-Smith, J. E., Lingner, C., Roberg-Perez, K. J., Tillberg, M., Brooks, J. E., and Garrels, J. I. (2000) The yeast proteome database (YPD) and Caenorhabditis elegans proteome database (WormPD): comprehensive resources for the organization and comparison of model organism protein information. Nucleic Acids Res. 28, 73 –76[Abstract/Free Full Text]

  22. Costanzo, M. C., Crawford, M. E., Hirschman, J. E., Kranz, J. E., Olsen, P., Robertson, L. S., Skrzypek, M. S., Braun, B. R., Hopkins, K. L., Kondu, P., Lengieza, C., Lew-Smith, J. E., Tillberg, M., and Garrels, J. I. (2001) YPD, PombePD and WormPD: model organism volumes of the BioKnowledge library, an integrated resource for protein information. Nucleic Acids Res. 29, 75 –79[Abstract/Free Full Text]

  23. Parrou, J. L., and Francois, J. (1997) A simplified procedure for a rapid and reliable assay of both glycogen and trehalose in whole yeast cells. Anal. Biochem. 248, 186 –188[CrossRef][Medline]

  24. Thomas, J. A., Schlender, K. K., and Larner, J. (1968) A rapid filter paper assay for UDPglucose-glycogen glucosyltransferase, including an improved biosynthesis of UDP-14C-glucose. Anal. Biochem. 25, 486 –499[Medline]

  25. Hardy, T. A., and Roach, P. J. (1993) Control of yeast glycogen synthase-2 by COOH-terminal phosphorylation. J. Biol. Chem. 268, 23799 –23805[Abstract/Free Full Text]

  26. 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[Abstract]

  27. Chester, V. E. (1968) Heritable glycogen-storage deficiency in yeast and its induction by ultra-violet light. J. Gen. Microbiol. 51, 49 –56[Medline]

  28. Enjalbert, B., Parrou, J. L., Vincent, O., and Francois, J. (2000) Mitochondrial respiratory mutants of Saccharomyces cerevisiae accumulate glycogen and readily mobilize it in a glucose-depleted medium. Microbiology 146, 2685 –2694[Abstract/Free Full Text]

  29. Cannon, J. F., Pringle, J. R., Fiechter, A., and Khalil, M. (1994) Characterization of glycogen-deficient glc mutants of Saccharomyces cerevisiae. Genetics 136, 485 –503[Abstract/Free Full Text]

  30. Heinisch, J. (1986) Construction and physiological characterization of mutants disrupted in the phosphofructokinase genes of Saccharomyces cerevisiae. Curr. Genet. 11, 227 –234[Medline]

  31. Thevelein, J. M., and Hohmann, S. (1995) Trehalose synthase: guard to the gate of glycolysis in yeast? Trends Biochem. Sci. 20, 3 –10[CrossRef][Medline]

  32. Huang, D., Wilson, W. A., and Roach, P. J. (1997) Glucose-6-P control of glycogen synthase phosphorylation in yeast. J. Biol. Chem. 272, 22495 –22501[Abstract/Free Full Text]

  33. Francois, J. M., Thompson-Jaeger, S., Skroch, J., Zellenka, U., Spevak, W., and Tatchell, K. (1992) GAC1 may encode a regulatory subunit for protein phosphatase type 1 in Saccharomyces cerevisiae. EMBO J. 11, 87 –96[Abstract]

  34. Stuart, J. S., Frederick, D. L., Varner, C. M., and Tatchell, K. (1994) The mutant type 1 protein phosphatase encoded by glc7–1 from Saccharomyces cerevisiae fails to interact productively with the GAC1-encoded regulatory subunit. Mol. Cell. Biol. 14, 896 –905[Abstract]

  35. Cheng, C., Huang, D., and Roach, P. J. (1997) Yeast PIG genes: PIG1 encodes a putative type 1 phosphatase subunit that interacts with the yeast glycogen synthase Gsy2p. Yeast 13, 1 –8[CrossRef][Medline]

  36. Tu, J. L., and Carlson, M. (1995) REG1 binds to protein phosphatase type 1 and regulates glucose repression in Saccharomyces cerevisiae. EMBO J. 14, 5939 –5946[Abstract]

  37. Entian, K. D., and Zimmermann, F. K. (1980) Glycolytic enzymes and intermediates in carbon catabolite repression mutants of Saccharomyces cerevisiae. Mol. Gen. Genet. 177, 345 –350[Medline]

  38. Huang, D., Chun, K. T., Goebl, M. G., and Roach, P. J. (1996) Genetic interactions between REG1/HEX2 and GLC7, the gene encoding the protein phosphatase type 1 catalytic subunit in Saccharomyces cerevisiae. Genetics 143, 119 –127[Abstract/Free Full Text]

  39. Wang, Z., Wilson, W. A., Fujino, M. A., and Roach, P. J. (2001) Antagonistic controls of autophagy and glycogen accumulation by Snf1p, the yeast homolog of AMP-activated protein kinase, and the cyclin-dependent kinase Pho85p. Mol. Cell. Biol. 21, 5742 –5752[Abstract/Free Full Text]

  40. Reinders, A., Burckert, N., Boller, T., Wiemken, A., and DeVirgilio, C. (1998) Saccharomyces cerevisiae cAMP-dependent protein kinase controls entry into stationary phase through the Rim15p protein kinase. Genes Dev. 12, 2943 –2955[Abstract/Free Full Text]

  41. Shu, Y., Yang, H., Hallberg, E., and Hallberg, R. (1997) Molecular genetic analysis of Rts1p, a B' regulatory subunit of Saccharomyces cerevisiae protein phosphatase 2A. Mol. Cell. Biol. 17, 3242 –3253[Abstract]

  42. Peng, Z. Y., Wang, W., Wilson, S. E., Schlender, K. K., Trumbly, R. J., and Reimann, E. M. (1991) Identification of a glycogen synthase phosphatase from yeast Saccharomyces cerevisiae as protein phosphatase 2A. J. Biol. Chem. 266, 10925 –10932[Abstract/Free Full Text]

  43. Posas, F., Clotet, J., Muns, M. T., Corominas, J., Casamayor, A., and Arino, J. (1993) The gene PPG encodes a novel yeast protein phosphatase involved in glycogen accumulation. J. Biol. Chem. 268, 1349 –1354[Abstract/Free Full Text]

  44. Clotet, J., Posas, F., Hu, G. Z., Ronne, H., and Arino, J. (1995) Role of protein phosphatase 2A in the control of glycogen metabolism in yeast. Eur. J. Biochem. 229, 207 –214[Abstract]

  45. Stark, M. J. R. (1996) Yeast protein serine/threonine phosphatases: multiple roles and diverse regulation. Yeast 12, 1647 –1675[CrossRef][Medline]

  46. Sakumoto, N., Mukai, Y., Uchida, K., Kouchi, T., Kuwajima, J., Nakagawa, Y., Sugioka, S., Yamamoto, E., Furuyama, T., Mizubuchi, H., Ohsugi, N., Sakuno, T., Kikuchi, K., Matsuoka, I., Ogawa, N., Kaneko, Y., and Harashima, S. (1999) A series of protein phosphatase gene disruptants in Saccharomyces cerevisiae. Yeast 15, 1669 –1679[CrossRef][Medline]

  47. Kim, J., and Klionsky, D. J. (2000) Autophagy, cytoplasm-to-vacuole targeting pathway, and pexophagy in yeast and mammalian cells. Annu. Rev. Biochem. 69, 303 –342[CrossRef][Medline]

  48. Klionsky, D. J., and Ohsumi, Y. (1999) Vacuolar import of proteins and organelles from the cytoplasm. Annu. Rev. Cell Dev. Biol. 15, 1 –32[CrossRef][Medline]

  49. Graham, L. A., Powell, B., and Stevens, T. H. (2000) Composition and assembly of the yeast vacuolar H+-ATPase complex. J. Exp. Biol. 203, 61 –70[Abstract]

  50. Ohya, Y., Umemoto, N., Tanida, I., Ohta, A., Iida, H., and Anraku, Y. (1991) Calcium-sensitive cls mutants of Saccharomyces cerevisiae showing a Pet-phenotype are ascribable to defects of vacuolar membrane H+-ATPase activity. J. Biol. Chem. 266, 13971 –13977[Abstract/Free Full Text]

  51. Nakamura, N., Matsuura, A., Wada, Y., and Ohsumi, Y. (1997) Acidification of vacuoles is required for autophagic degradation in the yeast, Saccharomyces cerevisiae. J. Biochem. 121, 338 –344[Abstract]

  52. Wolfe, K. H., and Shields, D. C. (1997) Molecular evidence for an ancient duplication of the entire yeast genome. Nature 387, 708 –713[CrossRef][Medline]

  53. Huang, D., Moffat, J., Wilson, W. A., Moore, L., Cheng, C., Roach, P. J., and Andrews, B. (1998) Cyclin partners determine Pho85 protein kinase substrate specificity in vitro and in vivo: control of glycogen biosynthesis by Pcl8 and Pcl10. Mol. Cell. Biol. 18, 3289 –3299[Abstract/Free Full Text]

  54. Lee, M., O’Regan, S., Moreau, J. L., Johnson, A. L., Johnston, L. H., and Goding, C. R. (2000) Regulation of the Pcl7-Pho85 cyclin-cdk complex by Pho81. Mol. Microbiol. 38, 411 –422[CrossRef][Medline]

  55. Pugh, T. A., Shah, J. C., Magee, P. T., and Clancy, M. J. (1989) Characterization and localization of the sporulation glucoamylase of Saccharomyces cerevisiae. Biochim. Biophys. Acta 994, 200 –209[Medline]