# High-dimensional and large-scale phenotyping of yeast mutants

Yoshikazu Ohya, Jun Sese, Masashi Yukawa, Fumi Sano, Yoichiro Nakatani, Taro Saito, Ayaka Saka, Tomoyuki Fukuda, Satoru Ishihara, Satomi Oka, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Yoshikazu Ohya, Jun Sese, Masashi Yukawa, Fumi Sano, Yoichiro Nakatani, et al. High-dimensional and large-scale phenotyping of yeast mutants. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102 (52), pp.19015-19020. 10.1073/pnas.0509436102. hal-02682283

HAL Id: hal-02682283

https://hal.inrae.fr/hal-02682283

Submitted on 1 Jun 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# High-dimensional and large-scale phenotyping of yeast mutants

Yoshikazu Ohya*†‡, Jun Sese§¶, Masashi Yukawa*†¶, Fumi Sano*†, Yoichiro Nakatani†§, Taro L. Saito†‖, Ayaka Saka*†, Tomoyuki Fukuda*, Satoru Ishihara*, Satomi Oka*, Genjiro Suzuki*, Machika Watanabe*, Aiko Hirata*†, Miwaka Ohtani*†§, Hiroshi Sawai*†§, Nicolas Fraysse**, Jean-Paul Latgé**, Jean M. François††, Markus Aebi†‡, Seiji Tanaka§§, Sachiko Muramatsu§§, Hiroyuki Araki§§, Kintake Sonoike*, Satoru Nogami*†, and Shinichi Morishita†§

Departments of *Integrated Biosciences and §Computational Biology, Graduate School of Frontier Sciences, University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8562, Japan; †Institute for Bioinformatics and Research and Development, Japan Science and Technology Corporation, Science Plaza 5-3 Yonbancho, Chiyoda-ku, Tokyo 102-8666, Japan; ‡Department of Computer Science, Graduate School of Information Science and Technology, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan; **Unité des Aspergillus, Institut Pasteur, 25 Rue du Docteur Roux, 75015 Paris, France; ††Centre de Bioingenierie Gilbert Durand, Unité Mixte de Recherche-Centre National de la Recherche Scientifique 5504, Unité Mixte de Recherche-Institut National de la Recherche Agronomique 792, F31077 Toulouse Cedex, France; ‡‡Institute of Microbiology, Swiss Federal Institute of Technology, Eidgenössische Technische Hochschule-Honggerberg, CH-8093 Zurich, Switzerland; and §§Division of Microbial Genetics, National Institute of Genetics, 1111 Yata, Mishima, Shizuoka 411-8540, Japan

Communicated by David Botstein, Princeton University, Princeton, NJ, October 31, 2005 (received for review August 8, 2005)

One of the most powerful techniques for attributing functions to genes in uni- and multicellular organisms is comprehensive analysis of mutant traits. In this study, systematic and quantitative analyses of mutant traits are achieved in the budding yeast *Saccharomyces cerevisiae* by investigating morphological phenotypes. Analysis of fluorescent microscopic images of triple-stained cells makes it possible to treat morphological variations as quantitative traits. Deletion of nearly half of the yeast genes not essential for growth affects these morphological traits. Similar morphological phenotypes are caused by deletions of functionally related genes, enabling a functional assignment of a locus to a specific cellular pathway. The high-dimensional phenotypic analysis of defined yeast mutant strains provides another step toward attributing gene function to all of the genes in the yeast genome.

cell morphology | functional genomics | high-dimensional phenotyping | phenome

One of the ultimate goals of genetics is to reveal relationships between gene function and phenotypic traits. Comprehensive analysis of mutant traits is a very powerful technique for attributing functions to genes in uni- and multicellular organisms. In the budding yeast *Saccharomyces cerevisiae*, a complete set of mutants, each of which carries a precise deletion of one yeast ORF, has been systematically constructed (1). By using these mutant strains combined with microarray and robot technology, genome-wide analyses of various mutant traits, including general growth rate, fitness under a particular condition, and sensitivity to drugs, has been reported (reviewed in ref. 2).

Cell morphology becomes an attractive target for comprehensive analysis, because more powerful methods for fluorescent microscopic imaging analysis in biological research have been emerging after development of high-resolution microscopes and specific fluorescent dyes. Yeast cell morphology reflects various cellular events, including progression through the cell cycle, establishment of cell polarity, and regulation of cell size control. Previous genome-wide studies of yeast morphology were focused on a specific morphology, such as cell size, cell shape, or bud site pattern (3–6), and therefore extracted limited information. Because morphological traits are often judged “by eye,” it has remained difficult to obtain quantitative and reproducible results.

We recently developed an image-processing system that automatically processes digital cell images of each yeast cell (7, 8) to obtain quantitative morphological data of yeast mutant cells. Mannoprotein (as a cell wall component marker), the actin cytoskeleton, and nuclear DNA are specifically stained simultaneously. Cells are then photographed, and fluorescence im-

ages are automatically processed. The obtained images of all yeast mutants and data-mining functions are available at our *Saccharomyces cerevisiae* Morphological Database (SCMD) web site (8, 9).

In this study, we employ high-dimensional and quantitative phenotyping of yeast mutants. This large-scale phenotyping represents the comprehensive detailed analyses of the SCMD morphological data set. Our quantitative and high-dimensional analysis is an improvement of previously published, rather qualitative analyses (10–16). Similar phenotypes are caused by deletions of functionally related genes, enabling a functional assignment of a locus to a specific cellular pathway.

## Materials and Methods

**Strains, Molecular Biology, and Functional Assays.** A set of haploid *S. cerevisiae* *MATA* deletion strains (4,786 strains) was obtained from the European *Saccharomyces cerevisiae* Archive for Functional Analysis (EUROSCARF). For obtaining *fks1*ts data, YOC1085 (*MATA* *fks1*-1144 *fks2*Δ) and its isogenic control strain YOC1001 were used. Each strain was grown in yeast extract/peptone/dextrose medium, and logarithmic-phase cells were fixed. To obtain fluorescent images of the cell-surface manno-protein, actin cytoskeleton, and nuclear DNA, cells were triply stained with fluorescein isothiocyanate-Con A, rhodamine-phalloidin, and 4',6-diamidino-2-phenylindole, respectively. Four deletion strains (*yal016w/tpd3*, *ylr131c/ace2*, *yor290c/snf2*, and *ylr425w/tus1*) could not be processed because of aggregation. Quantitative data were obtained from a total of 4,782 strains. We omitted 62 strains for further analysis because corresponding ORFs have been reported to be merged or deleted according to the updated version of the S288C sequence. Deletion strains of 4,718 ORFs were used. Disruption of the corresponding genes was verified by PCR for 81 randomly sampled strains. Standard molecular biological techniques were used (17).

**Estimation of Mutants.** Statistical tests are not applicable to our data set because we have only one set of data for each deletion strain; thus, we defined morphological normality of a deletion strain as the

<footnote>Conflict of interest statement: No conflicts declared.</footnote>

<footnote>Freely available online through the PNAS open access option.</footnote>

<footnote>Abbreviation: GO, gene ontology.</footnote>

<footnote>¶J.S. and M.Y. contributed equally to this work.</footnote>

<footnote>†To whom correspondence may be addressed. E-mail: ohya@k.u-tokyo.ac.jp or (for bioinformatics questions) moris@cb.k.u-tokyo.ac.jp.</footnote>

<footnote>© 2005 by The National Academy of Sciences of the USA</footnote>
A

<table><thead><tr><th></th><th>Cell shape</th><th>Actin</th><th>Nucleus</th><th>Total</th></tr></thead><tbody><tr><td>Number of parameters</td><td>133</td><td>96</td><td>272</td><td>501</td></tr><tr><td>Unbudded cell<br>with one nucleus</td><td>16</td><td>17</td><td>40</td><td>73</td></tr><tr><td>Budded cell<br>with one nucleus</td><td>55</td><td>32</td><td>66</td><td>153</td></tr><tr><td>Budded cell<br>with two nuclei</td><td>55</td><td>32</td><td>148</td><td>235</td></tr><tr><td>Cell number ratio</td><td>7</td><td>15</td><td>18</td><td>40</td></tr><tr><td>Statistically reliable</td><td>67</td><td>35</td><td>152</td><td>254</td></tr></tbody></table>

<figure><img src="image_3.png" /></figure>

C

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 1.** Comprehensive analysis of morphological phenotypes in yeast. (A) Classification of parameters. A subset of parameters used to characterize individual cells is shown. Among all of the 501 parameters measured, 254 revealed statistically reliable data (see *Supporting Text*). Thresholded at a correlation coefficient value of 0.9, 175 of 254 parameters were considered to be independent (see Fig. 6, which is published as supporting information on the PNAS web site). (B) Data sets were obtained from 126 independent wild-type samples and compared with all data sets from the 4,718 mutant strains. As an example, the distribution of the average of outline length of cells with no bud and a single nucleus is shown. (C) Number of morphological mutants distinct from wild type at the *P* value for both sides indicated. Deletion strains whose morphological normality *P* is less than or equal to the threshold in at least one parameter are counted as morphological mutants</figcaption></figure>

probability that data from wild-type cells would have measurements outside the range of the deletion strain. We estimated the distribution of wild-type data by transformation and discarded 247 parameters by a normality test because the estimated normal distributions for these parameters were not reliable. Assuming that the transformed wild-type distribution follows the estimated normal distribution in the remaining 254 parameters, we estimated the number of mutants by counting deletion strains that show abnormal morphology in at least 1 parameter. The detailed method is described in *Supporting Text*, which is published as supporting information on the PNAS web site.

**Gene Ontology (GO) Analysis.** We downloaded the GO-gene associations and the complete list of yeast ORFs annotated with GO terms from the *Saccharomyces* genome database as of Jan. 27, 2005.

**Functional Prediction.** For the “DNA recombinational repair” group, we selected 14 well defined mutant strains with a deletion in a gene locus affecting DNA repair. For the “cell wall biosynthesis” group, we selected a group of 42 mutant strains characterized by a reduced content of glucan in the cell wall (data not shown). To automatically predict the genes having the query functional module of interest from high-dimensional morphological data, we first selected morphological parameters that characterize the query functional module by using classified clustering (18). The classified clustering allows us to extract the most discriminative combination of parameters between the mutants that are known to be related to the query functional module and the other mutants in terms of interclass variance, a common statistical measure for classifying high-dimensional data. Focusing on the dimensions selected, we then separated the periphery of the mutants relevant to the query functional module from the others by using a support vector machine (19). The mutants in the periphery are predicted to be part of the query functional module.

**URLs.** The *Saccharomyces* genome database is available at www.yeastgenome.org, the GO database at www.geneontology.org, and the *Saccharomyces cerevisiae* Morphological Database (SCMD) at http://scmd.gi.k.u-tokyo.ac.jp.

## Results and Discussion

**Systematic Identification of Yeast Morphological Mutants.** To systematically identify yeast morphological mutants, we first defined parameters that accurately reflect yeast cellular morphology, then we used those parameters to find mutants whose morphology was significantly different from wild-type cells. We defined 501 quantitative parameters representing information about cell shape (visualized by cell wall staining), the actin cytoskeleton, and nuclear morphology of cells at a specific stage of the cell cycle (Fig. 1A; see also Table 1 and Fig. 5, which are published as supporting information on the PNAS web site). Of these 501 parameters, 254 were statistically reliable after power transformation of wild-type data (see *Supporting Text*). Of the 4,718 haploid mutants of the systematically constructed gene deletion collection (1), a total of 2,378 mutant strains exhibit differences from wild-type cells in at least one of the 254 morphological parameters with normal distribution (Fig. 1B and C, *P* < 0.0001 for one side). This estimation indicates that individual deletion of nearly half of the nonessential genes in the genome affects cellular morphology. Two hundred and forty-seven of the 254 parameters identify differences in at least one

<footnote>(see *Supporting Text* for details). (D) Frequency of the number of parameters altered in the strains analyzed. The number of deletion strains in each category with the GO term “biological process unknown” assigned to the corresponding locus is shown in yellow.</footnote>
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. 2.** Phenotypic analysis of mutant strains and correlation to cellular functions. (A) Polarisome mutants are enriched in the mutant class with altered bud morphology; 4,718 mutants were sorted by roundness of bud. Mutant strains that deviate from wild-type phenotype are shown in magenta or green. The positions of polarosome mutant strains and wild type are indicated with yellow and white lines, respectively. A round bud is a common feature for polarosome perturbation. (B) Cellular functions (as defined by GO terms) that are affected in subsets of mutant strains displaying specific morphological phenotypes. Mutant strains significantly different from wild-type cells are shown in green and magenta, respectively. Genes annotated to the GO term but that exhibit wild-type measurements are in black text. The relative positions of mutant strains with deletions in the loci given are indicated by yellow lines.</figcaption></figure>

mutant strain (see Table 2, which is published as supporting information on the PNAS web site). Among the 2,378 morphological mutants identified, 544 carry a deletion in a gene of unknown function (Fig. 1D and Table 3, which is published as supporting information on the PNAS web site), indicating that our high-dimensional analytical approach significantly advances

the phenotyping of the yeast mutant collection significantly. The statistical analysis of all yeast mutant strains analyzed, as well as the Java-based CALMORPH program, will be provided on request.

**Relationship Between Morphology and Gene Function.** One of the ultimate goals of genetics is to reveal relationships between gene
function and phenotypic traits. With the high-dimensional phenotypic traits, we demonstrate the close relationship between cellular morphology of the mutant strains and the functions of the deleted loci. As an example, we describe the analysis of polarisome function. The polarisome (Fig. 2A) is a protein complex that helps to determine cell polarity (20). Accordingly, four of five polarisome mutant strains (as provided by the GO consortium) display a round bud phenotype. The frequency of the round bud phenotype occurs by chance with a probability of $1.51 \times 10^{-8}$ (binomial test, Fig. 2A), indicating that a round bud is a characteristic trait of polarisome perturbation. A similar correlation is observed in 902 cases among a total of 368,808 combinations (1,452 GO groups × 254 parameters, $P < 0.01$, Bonferroni-corrected) (Fig. 2B; see also Tables 4 and 5, which are published as supporting information on the PNAS web site). Particular morphological phenotypes are enriched in 260 of the 1,452 GO groups analyzed. This finding implies that, in the case of the genes belonging to these annotated groups, the phenotype of a mutant strain deficient in such a locus can be predicted.

A specific phenotype can be assigned to a specific function. Fig. 3A shows the example of the collection of strains displaying the phenotype “large bud size before mitosis.” Of 86 strains within this group, 39 mutants carry a deletion of a gene annotated with “DNA metabolism” ($P < 2.2 \times 10^{-16}$ by chance); most of these loci are involved in DNA-damage response (21–26). Among the remaining 48 mutant strains with this phenotype, 16 corresponding loci are annotated to other DNA-related functions (27–29). Therefore, we conclude that mutant strains defective in DNA damage response represent the majority within the phenotypic class “large bud size before mitosis.” Likewise, for 109 of 254 morphological traits, the enrichment of mutant strains with deletions in loci with a specific annotation is observed (binomial test, $P < 0.01$, Bonferroni-corrected) (Fig. 3B; see also Tables 6 and 7, which are published as supporting information on the PNAS web site). Therefore, the function of a given gene can be predicted based on the morphological phenotype of the corresponding mutant strain.

**Function Prediction of the Gene from Morphological Data.** The validity of gene function prediction is improved by applying two or three phenotypic parameters for classification. We have chosen a set of deletion strains as a query group, by automatically selecting multiple phenotypic parameters that discriminate best between the strains of the query group and the rest of the strain collection. These parameters are then used to identify mutant strains in the whole collection that correspond to these phenotypes. We applied this approach to two groups of strains. One group is characterized as deficient in “DNA recombinational repair,” the other as impaired in “cell wall biosynthesis.” For the “DNA recombinational repair” group, we selected 14 well defined mutant strains with a deletion in a gene locus affecting DNA repair. This group of strains is characterized best by the parameters “ratio of cells in nuclear division” and “ratio of large-budded cells to budding cells before completion of nuclear division,” indicating that cells that arrest or delay at mitosis accumulated in these mutant strains. By using these query terms, a set of 19 mutant strains was identified in the collection, and 7 of these strains were included in the group of 14 strains that defined the query parameters (Fig. 4A). Importantly, the group of the 19 identified strains reveals additional common phenotypic characteristics as visualized by the presence of large contiguous patches of color representing similar morphological features (Fig. 4B). Wild-type cells treated with hydroxyurea, depleting the nucleoside triphosphate pool and inhibiting DNA synthesis, show morphological phenotypes similar to those of the group of the 19 mutant strains identified (Fig. 4C). Interestingly, all candidate strains identified show increased sensitivity to the DNA-damaging reagent hydroxyurea or methyl methanesulfo-

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 3.** Enrichment of defined mutant classes within a defined set of morphological phenotypes. (A) Mutant strains affected in DNA metabolism and DNA-related functions are enriched in strains displaying a large bud size before mitosis. Distributions of the mutants with a large bud size before mitosis are shown in magenta. The pie chart on the right represents the ratio of mutants with defects in DNA-related function within the morphological group “large bud size before mitosis.” Genes associated with the GO term “DNA metabolism” and genes involved in DNA-related function are shown in yellow and red, respectively. The pie chart on the left represents the ratio of genes associated with the GO term “DNA metabolism” to total nonessential mutants. (B) Mutants affected in specific cellular functions are enriched in mutant class displaying the defined phenotype. The cellular function is given above each statistical analysis. The enrichment is observed at indicated probability by chance. (Left) The diagram of each analysis displays the frequency of the mutant strains assigned to the given cellular function. (Right) The diagram displays the frequency of such defined strains within strains with the indicated morphological phenotype.</figcaption></figure>

nate (see Table 8, which is published as supporting information on the PNAS web site, and *Supporting Text*). Among these candidates, we further characterized *SAC3* and *DIA2*, whose functions in DNA repair are unknown, suggesting involvement of these genes in DNA repair (see Fig. 7, which is published as supporting information on the PNAS web site, and *Supporting Text*).

The similarity between the inhibitor-induced morphological phenotype and the morphological characteristics of the func-
<figure><img src="image_10.png" /></figure>

I. query: recombinational repair (14 genes)

<table><tr><td>ydr369c</td><td>xrs2*</td><td>double strand break repair</td></tr><tr><td>yer095w</td><td>rad51*</td><td>DNA recombinase</td></tr><tr><td>yg1163c</td><td>rad54*</td><td>DNA repair</td></tr><tr><td>yml032c</td><td>rad52*</td><td>DNA recombinase</td></tr><tr><td>ymr224c</td><td>mre11*</td><td>double strand break repair</td></tr><tr><td>ynl250w</td><td>rad50*</td><td>double strand break repair</td></tr><tr><td>ycr066w</td><td>rad18*</td><td>DNA repair</td></tr><tr><td>ydl013w</td><td>hex3</td><td>DNA recombination</td></tr><tr><td>yhr154w</td><td>rtt107</td><td>DNA transposition</td></tr><tr><td>yhr191c</td><td>ctf8</td><td>DNA replication factor C complex</td></tr><tr><td>yjl115w</td><td>asf1</td><td>DNA damage response</td></tr><tr><td>yll002w</td><td>rtt109</td><td>DNA transposition</td></tr><tr><td>ylr233c</td><td>est1</td><td>telomere length regulation</td></tr><tr><td>ymr078c</td><td>ctf18</td><td>DNA replication factor C complex</td></tr><tr><td>ypr135w</td><td>ctf4</td><td>DNA polymerase a binding protein</td></tr><tr><td>ypr164w</td><td>mms1</td><td>DNA repair</td></tr><tr><td>ydr159w</td><td>sac3</td><td>unknown</td></tr><tr><td>yor080w</td><td>dia2</td><td>unknown</td></tr><tr><td>ynl171c</td><td></td><td>unknown</td></tr></table>

II. query: low glucan content (42 genes)

<table><tr><td>ydl192w</td><td>arf1*</td><td>ER to Golgi transport</td></tr><tr><td>ydr129c</td><td>sac6*</td><td>actin organization, endocytosis</td></tr><tr><td>ydr245w</td><td>mnn10*</td><td>mannosyltransferase</td></tr><tr><td>yg1038c</td><td>och1*</td><td>mannosyltransferase</td></tr><tr><td>yg1223c</td><td>cog1*</td><td>Golgi transport complex</td></tr><tr><td>yor035c</td><td>she4*</td><td>actin organization</td></tr><tr><td>ypl045w</td><td>vps16*</td><td>Golgi to endosome transport</td></tr><tr><td>ygr036c</td><td>cax4</td><td>N-linked glycosylation</td></tr><tr><td>ynl084c</td><td>end3</td><td>actin organization, endocytosis</td></tr><tr><td>ykl007w</td><td>cap1</td><td>barbed-end actin filament capping</td></tr></table>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Fig. 4.** Functional classification of morphological mutants. (A) The descriptive morphological phenotype of a set of 14 recombinational repair mutant strains (I) and 42 mutant strains with low glucan content (II) was defined and used as a query to identify mutant strains that display the appropriate phenotype. Nineteen (for query I) and 10 (for query II) candidate strains were identified. The locus deleted in these strains, the gene name, and the putative function (as given by the GO term) is given. A * indicates mutant strains used as a query. (B) Color-coded representation of the morphological phenotypes displayed by these mutant strains. Morphological difference from wild type is shown for each mutant such that the magnitude is indicated by the intensity of the colors displayed. The brightest colors represent a significant difference (*P* < 0.00001) for one side. Magenta and green indicate that the mutant has a significantly higher or lower value than that of the wild type, respectively. The numerals in the vertical axis represent the following morphological parameters: 1, C12–2.A1B; 2, A7–2.A1B; 3, C117.A1B; 4, C118.A1B; 5, DCV114.A1B; 6, C125.A1B; 7, D118.A1B; 8, CCV115.C; 9, C118.C; 10, DCV112.C; 11, D109.C; 12, D125.C; 13, CCV104.C; 14, DCV151.C; 15, C117.C; 16, DCV106.C; 17, D103.C; 18, A108; 19, DCV145.C; 20, DCV146.C; 21, A107.A1B; 22, C123; 23, D202; 24, D213; 25, D216; 26, A110; 27, A119; 28, D110.A1B; 29, D207; 30, C106.A1B; and 31, D214. A precise parameter description is shown in Fig. 5 and Table 1. (C and D) Scatter plots of indicated mutant phenotypes of each strain analyzed. The candidate strains identified by using the specified group phenotypes are highlighted. The phenotype of wild-type strains treated with hydroxyurea (HU) is shown in C. (D) The phenotype of wild-type and *fks1-1144* mutant cells grown at the temperature indicated is shown.</figcaption></figure>
tionally defined set of mutant strains suggests that morphological
phenotyping can be used to address the target area of a given
chemical. The chemical-induced phenotype can be used to
identify mutant strains with similar phenotypes. Potentially, the
chemical targets the function(s) that are affected by the gene
deletions present in the identified mutant strains.

Next, we selected a group of mutant strains characterized by a reduced content of glucan in the cell wall. This group reveals the identifying parameters “the ratio of cells with large bud after nuclear division” and “roundness of bud.” When these parameters were used, 10 strains were identified, and 7 of these strains belonged to the original query group (Fig. 4A). The phenotype used for the characterization of these strains suggests that cytokinesis or cell separation is incomplete in these mutants. Nine of the 10 strains are sensitive to one of four cell wall damaging drugs (calcofluor white, caffeine, echinocandin B, or SDS) and have altered polysaccharide content in cell wall fractions. Transmission electron microscopic analysis reveals altered cell wall architecture in 8 of 10 strains (see Fig. 8, which is published as supporting information on the PNAS web site, and *Supporting Text*). A similar morphological phenotype is obtained in an *fks1-1144* mutant strain when grown under

semipermissive conditions (Fig. 4D). This mutant allele affects
the essential 1,3-β-glucan synthase activity and results in a
temperature-sensitive phenotype (30). Our results indicate that
cell wall integrity is impaired in the identified mutant strains,
resulting in very similar morphological phenotypes.

In summary, our quantitative, high-dimensional phenotype analysis reveals that deletions of about half of the nonessential genes in yeast affect cell morphology. Therefore, our approach reveals a phenotype for ≈500 yeast mutant strains with a deletion in a locus of unknown function. Importantly, morphological phenotypes strongly correlate with gene function. We propose to use morphological phenotyping of yeast mutant strains to assign potential functions to unknown genes. In addition, chemically induced morphological phenotypes of yeast can be used to rapidly identify cellular targets of drugs.

We thank Tamao Goto, Yuka Kitamura, and Emi Shimoi for technical assistance; Kara Dolinski and Fugaku Aoki for critical reading of the manuscript; and the members of the Ohya and Morishita groups at the University of Tokyo for stimulating discussions. This work was supported by the Institute for Bioinformatics and Research and Development of the Japan Science and Technology Corporation.

1. Winzeler, E. A., Shoemaker, D. D., Astromoff, A., Liang, H., Anderson, K., Andre, B., Bangham, R., Benito, R., Boeke, J. D., Bussey, H., et al. (1999) *Science* **285**, 901–906.

2. Bader, G. D., Heilbut, A., Andrews, B., Tyers, M., Hughes, T. & Boone, C. (2003) *Trends Cell Biol.* **13**, 344–356.

3. Giaever, G., Chu, A. M., Ni, L., Connelly, C., Riles, L., Veronneau, S., Dow, S., Lucau-Danila, A., Anderson, K., Andre, B., et al. (2002) *Nature* **418**, 387–391.

4. Jorgensen, P., Nishikawa, J. L., Breitkreutz, B. J. & Tyers, M. (2002) *Science* **297**, 395–400.

5. Zhang, J., Schneider, C., Ottmers, L., Rodriguez, R., Day, A., Markwardt, J. & Schneider, B. L. (2002) *Curr. Biol.* **12**, 1992–2001.

6. Ni, L. & Snyder, M. (2001) *Mol. Biol. Cell* **12**, 2147–2170.

7. Ohtani, M., Saka, A., Sano, F., Ohya, Y. & Morishita, S. (2004) *J. Bioinform. Comput. Biol.* **1**, 695–709.

8. Saito, T. L., Ohtani, M., Sawai, H., Sano, F., Saka, A., Watanabe, D., Yukawa, M., Ohya, Y. & Morishita, S. (2004) *Nucleic Acids Res.* **32**, D319–D322.

9. Saito, T. L., Sese, J., Nakatani, Y., Sano, F., Yukawa, M., Ohya, Y. & Morishita, S. (2005) *Nucleic Acids Res.* **33**, W753–W757.

10. Christie, K. R., Weng, S., Balakrishnan, R., Costanzo, M. C., Dolinski, K., Dwight, S. S., Engel, S. R., Feierbach, B., Fisk, D. G., Hirschman, J. E., et al. (2004) *Nucleic Acids Res.* **32**, D311–D314.

11. Guldener, U., Munsterkotter, M., Kastenmuller, G., Strack, N., van Helden, J., Lemer, C., Richelles, J., Wodak, S. J., Garcia-Martinez, J., Perez-Ortin, J. E., et al. (2005) *Nucleic Acids Res.* **33**, D364–D368.

12. Spellman, P. T., Sherlock, G., Zhang, M. Q., Iyer, V. R., Anders, K., Eisen, M. B., Brown, P. O., Botstein, D. & Futcher, B. (1998) *Mol. Biol. Cell* **9**, 3273–3297.

13. Huh, W. K., Falvo, J. V., Gerke, L. C., Carroll, A. S., Howson, R. W., Weissman, J. S. & O'Shea, E. K. (2003) *Nature* **425**, 686–691.

14. Uetz, P., Giot, L., Cagney, G., Mansfield, T. A., Judson, R. S., Knight, J. R., Lockshon, D., Narayan, V., Srinivasan, M., Pochart, P., et al. (2000) *Nature* **403**, 623–627.

15. Ito, T., Chiba, T., Ozawa, R., Yoshida, M., Hattori, M. & Sakaki, Y. (2001) *Proc. Natl. Acad. Sci. USA* **98**, 4569–4574.

16. Ho, Y., Gruhler, A., Heilbut, A., Bader, G. D., Moore, L., Adams, S. L., Millar, A., Taylor, P., Bennett, K., Boutilier, K., et al. (2002) *Nature* **415**, 180–183.

17. Burke, D., Dawson, D. & Stearns, T. (2000) *Methods in Yeast Genetics: A Cold Spring Harbor Laboratory Course Manual* (Cold Spring Harbor Lab. Press, Woodbury, New York).

18. Sese, J., Kurokawa, Y., Monden, M., Kato, K. & Morishita, S. (2004) *Bioinformatics* **20**, 3137–3145.

19. Vapnik, V. (1995) *The Nature of the Statistical Learning Theory* (Springer, Berlin).

20. Pruyne, D. & Bretscher, A. (2000) *J. Cell Sci.* **113**, 365–375.

21. Broomfield, S., Hryciw, T. & Xiao, W. (2001) *Mutat. Res.* **486**, 167–184.

22. Symington, L. S. (2002) *Microbiol. Mol. Biol. Rev.* **66**, 630–670.

23. Kouprina, N., Kroll, E., Bannikov, V., Bliskovsky, V., Gizatullin, R., Kirillov, A., Shestopalov, B., Zakharyev, V., Hieter, P., Spencer, F., et al. (1992) *Mol. Cell. Biol.* **12**, 5736–5747.

24. Mullen, J. R., Kaliraman, V., Ibrahim, S. S. & Brill, S. J. (2001) *Genetics* **157**, 103–118.

25. Scholes, D. T., Banerjee, M., Bowen, B. & Curcio, M. J. (2001) *Genetics* **159**, 1449–1465.

26. Hryciw, T., Tang, M., Fontanie, T. & Xiao, W. (2002) *Mol. Genet. Genomics* **266**, 848–857.

27. Ramey, C. J., Howar, S., Adkins, M., Linger, J., Spicer, J. & Tyler, J. K. (2004) *Mol. Cell. Biol.* **24**, 10313–10327.

28. Taggart, A. K. & Zakian, V. A. (2003) *Curr. Opin. Cell Biol.* **15**, 275–280.

29. Wintersberger, U., Kuhne, C. & Karwan, A. (1995) *Yeast* **11**, 929–944.

30. Dijkstra, G. J., Abe, M., Ohya, Y. & Bussey, H. (2002) *Yeast* **19**, 671–690.