Functional profiling of the Saccharomyces cerevisiae genome

Glutathione, Altruistic Metabolite in Fungi

Glutathione (GSH; gamma-L-glutamyl-L-cysteinyl-glycine), a non-protein thiol with a very low redox potential (E'0 = 240 mV for thiol-disulfide exchange), is present in high concentration up to 10 mM in yeasts and filamentous fungi. GSH is concerned with basic cellular functions as well as the maintenance of mitochondrial structure, membrane integrity, and in cell differentiation and development. GSH plays key roles in the response to several stress situations in fungi. For example, GSH is an important antioxidant molecule, which reacts non-enzymatically with a series of reactive oxygen species. In addition, the response to oxidative stress also involves GSH biosynthesis enzymes, NADPH-dependent GSH-regenerating reductase, glutathione S-transferase along with peroxide-eliminating glutathione peroxidase and glutaredoxins. Some components of the GSH-dependent antioxidative defence system confer resistance against heat shock and osmotic stress. Formation of protein-SSG mixed disulfides results in protection against desiccation-induced oxidative injuries in lichens. Intracellular GSH and GSH-derived phytochelatins hinder the progression of heavy metal-initiated cell injuries by chelating and sequestering the metal ions themselves and/or by eliminating reactive oxygen species. In fungi, GSH is mobilized to ensure cellular maintenance under sulfur or nitrogen starvation. Moreover, adaptation to carbon deprivation stress results in an increased tolerance to oxidative stress, which involves the induction of GSH-dependent elements of the antioxidant defence system. GSH-dependent detoxification processes concern the elimination of toxic endogenous metabolites, such as excess formaldehyde produced during the growth of the methylotrophic yeasts, by formaldehyde dehydrogenase and methylglyoxal, a by-product of glycolysis, by the glyoxalase pathway. Detoxification of xenobiotics, such as halogenated aromatic and alkylating agents, relies on glutathione S-transferases. In yeast, these enzymes may participate in the elimination of toxic intermediates that accumulate in stationary phase and/or act in a similar fashion as heat shock proteins. GSH S-conjugates may also form in a glutathione S-transferases-independent way, e.g. through chemical reaction between GSH and the antifugal agent Thiram. GSH-dependent detoxification of penicillin side-chain precursors was shown in Penicillium sp. GSH controls aging and autolysis in several fungal species, and possesses an anti-apoptotic feature.

The Mycobacterium tuberculosis complex transcriptome of attenuation

Although the deletion of RD1 is likely correlated to attenuation from virulence for members of the Mycobacterium tuberculosis (MTB) complex, the reasons for this phenotype remain to be fully explained. As genomic variation is responsible for at least a component of variability in gene expression, we looked to the in vitro global expression profile of the RD1 artificial knockout from M. tuberculosis H37Rv (H37Rv:deltaRD1) for clues to elucidate its phenotypic shift towards attenuation. By comparing the transcriptome of H37Rv:deltaRD1 to that of virulent H37Rv, 15 regulated genes located in nine different regions outside of RD1 have been identified, capturing an effect of RD1's deletion on the rest of the genome. To assess whether these regulations are characteristic of attenuated MTB in general, expression profiles of natural RD1 mutants (BCG Russia, BCG Pasteur, and M. microti) as well as the 'avirulent' M. tuberculosis H37Ra, whose RD1 region is genomically intact, were obtained. Results indicate that attenuated strains lack the expression of RD1 genes including cfp10 and esat6, whether through deletion or reduced expression. Furthermore, comparative transcriptomics reveals the concurrent down-regulation of several gene neighborhoods beyond RD1. The potential relevance of these other expression changes towards MTB virulence is discussed.

Functional genomic hypothesis generation and experimentation by a robot scientist

The question of whether it is possible to automate the scientific process is of both great theoretical interest and increasing practical importance because, in many scientific areas, data are being generated much faster than they can be effectively analysed. We describe a physically implemented robotic system that applies techniques from artificial intelligence to carry out cycles of scientific experimentation. The system automatically originates hypotheses to explain observations, devises experiments to test these hypotheses, physically runs the experiments using a laboratory robot, interprets the results to falsify hypotheses inconsistent with the data, and then repeats the cycle. Here we apply the system to the determination of gene function using deletion mutants of yeast (Saccharomyces cerevisiae) and auxotrophic growth experiments. We built and tested a detailed logical model (involving genes, proteins and metabolites) of the aromatic amino acid synthesis pathway. In biological experiments that automatically reconstruct parts of this model, we show that an intelligent experiment selection strategy is competitive with human performance and significantly outperforms, with a cost decrease of 3-fold and 100-fold (respectively), both cheapest and random-experiment selection.

Multidrug efflux pumps and resistance: regulation and evolution

Multidrug transporters have long puzzled researchers because of their ability to extrude multiple structurally dissimilar toxic chemotherapeutic agents. They appear to be essentially ubiquitous with many microorganisms possessing large arsenals of predicted multidrug efflux transporters. Recent reports on the structures of multidrug efflux transporters and their cognate regulators have suggested that the basis of multidrug recognition is the presence of large flexible hydrophobic cavities capable of accommodating different compounds that are bound via hydrogen bonding and/or electrostatic interactions. This structural data provides the context for understanding other questions about the evolution of multidrug transporters and their natural physiological roles.

Mapping eIF5A binding sites for Dys1 and Lia1: in vivo evidence for regulation of eIF5A hypusination

The evolutionarily conserved factor eIF5A is the only protein known to undergo hypusination, a unique posttranslational modification triggered by deoxyhypusine synthase (Dys1). Although eIF5A is essential for cell viability, the function of this putative translation initiation factor is still obscure. To identify eIF5A-binding proteins that could clarify its function, we screened a two-hybrid library and identified two eIF-5A partners in S. cerevisiae: Dys1 and the protein encoded by the gene YJR070C, named Lia1 (Ligand of eIF5A). The interactions were confirmed by GST pulldown. Mapping binding sites for these proteins revealed that both eIF5A domains can bind to Dys1, whereas the C-terminal domain is sufficient to bind Lia1. We demonstrate for the first time in vivo that the N-terminal alpha-helix of Dys1 can modulate enzyme activity by inhibiting eIF5A interaction. We suggest that this inhibition be abrogated in the cell when hypusinated and functional eIF5A is required.

High-resolution yeast phenomics resolves different physiological features in the saline response

We present a methodology for gene functional prediction based on extraction of physiologically relevant growth variables from all viable haploid yeast knockout mutants. This quantitative phenomics approach, here applied to saline cultivation, identified marginal but functionally important phenotypes and allowed the precise determination of time to adapt to an environmental challenge, rate of growth, and efficiency of growth. We identified approximately 500 salt-sensitive gene deletions, the majority of which were previously uncharacterized and displayed salt sensitivity for only one of the three physiological features. We also report a high correlation to protein-protein interaction data; in particular, several salt-sensitive subcellular networks indicating functional modules were revealed. In contrast, no correlation was found between gene dispensability and gene expression. It is proposed that high-resolution phenomics will be instrumental in systemwide descriptions of intragenomic functional networks.

Gene loss, protein sequence divergence, gene dispensability, expression level, and interactivity are correlated in eukaryotic evolution

Lineage-specific gene loss, to a large extent, accounts for the differences in gene repertoires between genomes, particularly among eukaryotes. We derived a parsimonious scenario of gene losses for eukaryotic orthologous groups (KOGs) from seven complete eukaryotic genomes. The scenario involves substantial gene loss in fungi, nematodes, and insects. Based on this evolutionary scenario and estimates of the divergence times between major eukaryotic phyla, we introduce a numerical measure, the propensity for gene loss (PGL). We explore the connection among the propensity of a gene to be lost in evolution (PGL value), protein sequence divergence, the effect of gene knockout on fitness, the number of protein-protein interactions, and expression level for the genes in KOGs. Significant correlations between PGL and each of these variables were detected. Genes that have a lower propensity to be lost in eukaryotic evolution accumulate fewer substitutions in their protein sequences and tend to be essential for the organism viability, tend to be highly expressed, and have many interaction partners. The dependence between PGL and gene dispensability and interactivity is much stronger than that for sequence evolution rate. Thus, propensity of a gene to be lost during evolution seems to be a direct reflection of its biological importance.

Yeast Genes That Enhance the Toxicity of a Mutant Huntingtin Fragment or -Synuclein

Genome-wide screens were performed in yeast to identify genes that enhance the toxicity of a mutant huntingtin fragment or of alpha-synuclein. Of 4850 haploid mutants containing deletions of nonessential genes, 52 were identified that were sensitive to a mutant huntingtin fragment, 86 that were sensitive to alpha-synuclein, and only one mutant that was sensitive to both. Genes that enhanced toxicity of the mutant huntingtin fragment clustered in the functionally related cellular processes of response to stress, protein folding, and ubiquitin-dependent protein catabolism, whereas genes that modified alpha-synuclein toxicity clustered in the processes of lipid metabolism and vesicle-mediated transport. Genes with human orthologs were overrepresented in our screens, suggesting that we may have discovered conserved and nonoverlapping sets of cell-autonomous genes and pathways that are relevant to Huntington's disease and Parkinson's disease.

Applicability of tandem affinity purification MudPIT to pathway proteomics in yeast

A combined multidimensional chromatography-mass spectrometry approach known as "MudPIT" enables rapid identification of proteins that interact with a tagged bait while bypassing some of the problems associated with analysis of polypeptides excised from SDS-polyacrylamide gels. However, the reproducibility, success rate, and applicability of MudPIT to the rapid characterization of dozens of proteins have not been reported. We show here that MudPIT reproducibly identified bona fide partners for budding yeast Gcn5p. Additionally, we successfully applied MudPIT to rapidly screen through a collection of tagged polypeptides to identify new protein interactions. Twenty-five proteins involved in transcription and progression through mitosis were modified with a new tandem affinity purification (TAP) tag. TAP-MudPIT analysis of 22 yeast strains that expressed these tagged proteins uncovered known or likely interacting partners for 21 of the baits, a figure that compares favorably with traditional approaches. The proteins identified here comprised 102 previously known and 279 potential physical interactions. Even for the intensively studied Swi2p/Snf2p, the catalytic subunit of the Swi/Snf chromatin remodeling complex, our analysis uncovered a new interacting protein, Rtt102p. Reciprocal tagging and TAP-MudPIT analysis of Rtt102p revealed subunits of both the Swi/Snf and RSC complexes, identifying Rtt102p as a common interactor with, and possible integral component of, these chromatin remodeling machines. Our experience indicates it is feasible for an investigator working with a single ion trap instrument in a conventional molecular/cellular biology laboratory to carry out proteomic characterization of a pathway, organelle, or process (i.e. "pathway proteomics") by systematic application of TAP-MudPIT.

Toward a functional annotation of the human genome using artificial transcription factors

We have developed a novel, high-throughput approach to collecting randomly perturbed gene-expression profiles from the human genome.A human 293 cell library that stably expresses randomly chosen zinc-finger transcription factors was constructed, and the expression profile of each cell line was obtained using cDNA microarray technology.Gene expression profiles from a total of 132 cell lines were collected and analyzed by (1) a simple clustering method based on expression-profile similarity, and (2) the shortest-path analysis method. These analyses identified a number of gene groups, and further investigation revealed that the genes that were grouped together had close biological relationships. The artificial transcription factor-based random genome perturbation method thus provides a novel functional genomic tool for annotation and classification of genes in the human genome and those of many other organisms.

Principles of transcriptional control in the metabolic network of Saccharomyces cerevisiae

Cellular networks are subject to extensive regulation, which modifies the availability and efficiency of connections between components in response to external conditions. Thus far, studies of large-scale networks have focused on their connectivity, but have not considered how the modulation of this connectivity might also determine network properties. To address this issue, we analyzed how the coordinated expression of enzymes shapes the metabolic network of Saccharomyces cerevisiae. By integrating large-scale expression data with the structural description of the metabolic network, we systematically characterized the transcriptional regulation of metabolic pathways. The analysis revealed recurrent patterns, which may represent design principles of metabolic gene regulation. First, we find that transcription regulation biases metabolic flow toward linearity by coexpressing only distinct branches at metabolic branchpoints. Second, individual isozymes were often separately coregulated with distinct processes, providing a means of reducing crosstalk between pathways using a common reaction. Finally, transcriptional regulation defined a hierarchical organization of metabolic pathways into groups of varying expression coherence. These results emphasize the utility of incorporating regulatory information when analyzing properties of large-scale cellular networks.

Negative regulation of calcineurin signaling by Hrr25p, a yeast homolog of casein kinase I

Calcineurin is a Ca2+/calmodulin-regulated protein phosphatase required for Saccharomyces cerevisiae to respond to a variety of environmental stresses. Calcineurin promotes cell survival during stress by dephosphorylating and activating the Zn-finger transcription factor Crz1p/Tcn1p. Using a high-throughput assay, we screened 119 yeast kinases for their ability to phosphorylate Crz1p in vitro and identified the casein kinase I homolog Hrr25p. Here we show that Hrr25p negatively regulates Crz1p activity and nuclear localization in vivo. Hrr25p binds to and phosphorylates Crz1p in vitro and in vivo. Overexpression of Hrr25p decreases Crz1p-dependent transcription and antagonizes its Ca2+-induced nuclear accumulation. In the absence of Hrr25p, activation of Crz1p by Ca2+/calcineurin is potentiated. These findings represent the first identification of a negative regulator for Crz1p, and establish a novel physiological role for Hrr25p in antagonizing calcineurin signaling.

Rds3p is required for stable U2 snRNP recruitment to the splicing apparatus

Rds3p is a well-conserved 12-kDa protein with five CxxC zinc fingers that has been implicated in the activation of certain drug transport genes and in the pre-mRNA splicing pathway. Here we show that Rds3p resides in the yeast spliceosome and is essential for splicing in vitro. Rds3p purified from yeast stably associates with at least five U2 snRNP proteins, Cus1p, Hsh49p, Hsh155p, Rse1p, and Ist3p/Snu17p, and with the Yra1p RNA export factor. A mutation upstream of the first Rds3p zinc finger causes the conditional release of the putative branchpoint nucleotide binding protein, Ist3p/Snu17p, and weakens Rse1p interaction with the Rds3p complex. The resultant U2 snRNP particle migrates exceptionally slowly in polyacrylamide gels, suggestive of a disorganized structure. U2 snRNPs depleted of Rds3p fail to form stable prespliceosomes, although U2 snRNA stability is not affected. Metabolic depletion of Yra1p blocks cell growth but not splicing, suggesting that Yra1p association with Rds3p relates to Yra1p's role in RNA trafficking. Together these data establish Rds3p as an essential component of the U2 snRNP SF3b complex and suggest a new link between the nuclear processes of pre-mRNA splicing and RNA export.

Pex30p, Pex31p, and Pex32p form a family of peroxisomal integral membrane proteins regulating peroxisome size and number in Saccharomyces cerevisiae

The peroxin Pex23p of the yeast Yarrowia lipolytica exhibits high sequence similarity to the hypothetical proteins Ylr324p, Ygr004p, and Ybr168p encoded by the Saccharomyces cerevisiae genome. Ylr324p, Ygr004p, and Ybr168p are integral to the peroxisomal membrane and act to control peroxisome number and size. Synthesis of Ylr324p and Ybr168p, but not of Ygr004p, is induced during incubation of cells in oleic acid-containing medium, the metabolism of which requires intact peroxisomes. Cells deleted for YLR324w exhibit increased numbers of peroxisomes, whereas cells deleted for YGR004w or YBR168w exhibit enlarged peroxisomes. Ylr324p and Ybr168p cannot functionally substitute for one another or for Ygr004p, whereas Ygr004p shows partial functional redundancy with Ylr324p and Ybr168p. Ylr324p, Ygr004p, and Ybr168p interact within themselves and with Pex28p and Pex29p, which have been shown also to regulate peroxisome size and number. Systematic deletion of genes demonstrated that PEX28 and PEX29 function upstream of YLR324w, YGR004w, and YBR168w in the regulation of peroxisome proliferation. Our data suggest a role for Ylr324p, Ygr004p, and Ybr168p--now designated Pex30p, Pex31p, and Pex32p, respectively--together with Pex28p and Pex29p in controlling peroxisome size and proliferation in Saccharomyces cerevisiae.

Comprehensive transposon mutant library of Pseudomonas aeruginosa

We have developed technologies for creating saturating libraries of sequence-defined transposon insertion mutants in which each strain is maintained. Phenotypic analysis of such libraries should provide a virtually complete identification of nonessential genes required for any process for which a suitable screen can be devised. The approach was applied to Pseudomonas aeruginosa, an opportunistic pathogen with a 6.3-Mbp genome. The library that was generated consists of 30,100 sequence-defined mutants, corresponding to an average of five insertions per gene. About 12% of the predicted genes of this organism lacked insertions; many of these genes are likely to be essential for growth on rich media. Based on statistical analyses and bioinformatic comparison to known essential genes in E. coli, we estimate that the actual number of essential genes is 300-400. Screening the collection for strains defective in two defined multigenic processes (twitching motility and prototrophic growth) identified mutants corresponding to nearly all genes expected from earlier studies. Thus, phenotypic analysis of the collection may produce essentially complete lists of genes required for diverse biological activities. The transposons used to generate the mutant collection have added features that should facilitate downstream studies of gene expression, protein localization, epistasis, and chromosome engineering.

The proteome of Saccharomyces cerevisiae mitochondria

We performed a comprehensive approach to determine the proteome of Saccharomyces cerevisiae mitochondria. The proteins of highly pure yeast mitochondria were separated by several independent methods and analyzed by tandem MS. From >20 million MS spectra, 750 different proteins were identified, indicating an involvement of mitochondria in numerous cellular processes. All known components of the oxidative phosphorylation machinery, the tricarboxylic acid cycle, and the stable mitochondria-encoded proteins were found. Based on the mitochondrial proteins described in the literature so far, we calculate that the identified proteins represent approximately 90% of all mitochondrial proteins. The function of a quarter of the identified proteins is unknown. The mitochondrial proteome will provide an important database for the analysis of new mitochondrial and mitochondria-associated functions and the characterization of mitochondrial diseases.

Novel cytochromes P450 applications arising from the directed-evolution of recombinant micro-organisms

AIMS

Directed (forced) evolution of cytochromes P450 (overall 2700 CYP isoforms in non-recombinant biota) is a method that has been investigated in yeasts (and other micro-organisms) by aerobically growing brewers' yeast Saccharomyces cerevisiae in very high glucose (20%) media.

METHODS AND RESULTS

Mitochondrial repression subverts cytochrome oxidase biosynthesis into manifest cytochromes P450 accumulation in brewers' yeast. A similar phenomenon is observed with the acridine-induced petit mutant. Cytochromes P450 EC 1.14.14.1 (and mimics) display a range of redox iron-mediated bioconversions in food processing, with mixed function oxidase (O2:mono-oxygenase) intervention results. Unfortunately these enzymes generate reactive oxygen species (ROS) through redox electron recycling, whilst isoform CYP 1A1 can activate precarcinogens such as benzo(a)pyrene to the ultimate (proximate) carcinogen that binds to nuclear DNA.

CONCLUSIONS

In conclusion, another 5000 CYP isoforms, for example, might be identified in micro-organisms and many more made to order through recombinant DNA technology and utilized both in vitro and in vivo for aimed bioconversions in industry and in the environment, as part of the impact of greener-approach supporting strategies to minimize global pollution.

Global analysis of protein localization in budding yeast

A fundamental goal of cell biology is to define the functions of proteins in the context of compartments that organize them in the cellular environment. Here we describe the construction and analysis of a collection of yeast strains expressing full-length, chromosomally tagged green fluorescent protein fusion proteins. We classify these proteins, representing 75% of the yeast proteome, into 22 distinct subcellular localization categories, and provide localization information for 70% of previously unlocalized proteins. Analysis of this high-resolution, high-coverage localization data set in the context of transcriptional, genetic, and protein-protein interaction data helps reveal the logic of transcriptional co-regulation, and provides a comprehensive view of interactions within and between organelles in eukaryotic cells.

Comparative genomics, minimal gene-sets and the last universal common ancestor

Comparative genomics, using computational and experimental methods, enables the identification of a minimal set of genes that is necessary and sufficient for sustaining a functional cell. For most essential cellular functions, two or more unrelated or distantly related proteins have evolved; only about 60 proteins, primarily those involved in translation, are common to all cellular life. The reconstruction of ancestral life-forms is based on the principle of evolutionary parsimony, but the size and composition of the reconstructed ancestral gene-repertoires depend on relative rates of gene loss and horizontal gene-transfer. The present estimate suggests a simple last universal common ancestor with only 500-600 genes.

An analysis of determinants of amino acids substitution rates in bacterial proteins

The variation of amino acid substitution rates in proteins depends on several variables. Among these, the protein's expression level, functional category, essentiality, or metabolic costs of its amino acid residues may play an important role. However, the relative importance of each variable has not yet been evaluated in comparative analyses. To this aim, we made regression analyses combining data available on these variables and on evolutionary rates, in two well-documented model bacteria, Escherichia coli and Bacillus subtilis. In both bacteria, the level of expression of the protein in the cell was by far the most important driving force constraining the amino acids substitution rate. Subsequent inclusion in the analysis of the other variables added little further information. Furthermore, when the rates of synonymous substitutions were included in the analysis of the E. coli data, only the variable expression levels remained statistically significant. The rate of nonsynonymous substitution was shown to correlate with expression levels independently of the rate of synonymous substitution. These results suggest an important direct influence of expression levels, or at least codon usage bias for translation optimization, on the rates of nonsynonymous substitutions in bacteria. They also indicate that when a control for this variable is included, essentiality plays no significant role in the rate of protein evolution in bacteria, as is the case in eukaryotes.

DNA helicase gene interaction network defined using synthetic lethality analyzed by microarray

We describe a new synthetic lethality analysis by microarray (SLAM) technique that uses approximately 4,600 Saccharomyces cerevisiae haploid deletion mutants with molecular 'bar codes' (TAGs). We used SGS1 and SRS2, two 3'-->5' DNA helicase genes, as 'queries' to identify their redundant and unique biological functions. We introduced these 'query mutations' into a haploid deletion pool by integrative transformation to disrupt the query gene in every cell, generating a double mutant pool. Optimization of integrative transformation efficiency was essential to the success of SLAM. Synthetic interactions defined a DNA helicase genetic network and predicted a role for SRS2 in processing damaged replication forks but, unlike SGS1, not in rDNA replication, DNA topology or lagging strand synthesis. SGS1 and SRS2 have synthetic defects with MRC1 but not RAD9, suggesting that SGS1 and SRS2 function in a parallel pathway with MRC1 to transduce the DNA replication stress signal to the general DNA damage checkpoint pathway. Both helicase genes have rad51-reversible synthetic defects with 5'-->3' DNA helicase RRM3, suggesting that RRM3 helps prevent formation of toxic recombination intermediates. SLAM detects synthetic lethality efficiently and ranks candidate genetic interactions, making it an especially useful method.

Cardiolipin is not required for Bax-mediated cytochrome c release from yeast mitochondria

Cardiolipin (CL) is an inner mitochondrial membrane phospholipid that contributes to optimal mitochondrial function and is gaining widespread attention in studies of mitochondria-mediated apoptosis. Divergent hypotheses describing the role of CL in cytochrome c release and apoptosis have evolved. We addressed this controversy directly by comparing the spontaneous- and Bax-mediated cytochrome c release from mitochondria isolated from two strains of Saccharomyces cerevisiae: one lacking CL-synthase and therefore CL (DeltaCRD1) and the other, its corresponding wild type (WT). We demonstrated by liquid chromatography-mass spectrometry that the main yeast CL species [(16:1)2(18:1)2] differs in fatty acid composition from mammalian CL [(18:2)4], and we verified the absence of the yeast CL species in the DeltaCRD1 strain. We also demonstrated that the mitochondrial association of Bax and the resulting cytochrome c release is not dependent on the CL content of the yeast mitochondrial membranes. Bax inserted equally into both WT and DeltaCRD1 mitochondrial membranes under conditions that lead to the release of cytochrome c from both strains of yeast mitochondria. Furthermore, using models of synthetic liposomes and isolated yeast mitochondria, we found that cytochrome c was bound more "loosely" to the CL-deficient systems compared with when CL is present. These data challenge recent studies implicating that CL is required for Bax-mediated pore formation leading to the release of proteins from the mitochondrial intermembrane space. In contrast, they support our recently proposed two-step mechanism of cytochrome c release, which suggests that CL is required for binding cytochrome c to the inner mitochondrial membrane.

Integrating ‘omic’ information: a bridge between genomics and systems biology

The availability of genome sequences for several organisms, including humans, and the resulting first-approximation lists of genes, have allowed a transition from molecular biology to 'modular biology'. In modular biology, biological processes of interest, or modules, are studied as complex systems of functionally interacting macromolecules. Functional genomic and proteomic ('omic') approaches can be helpful to accelerate the identification of the genes and gene products involved in particular modules, and to describe the functional relationships between them. However, the data emerging from individual omic approaches should be viewed with caution because of the occurrence of false-negative and false-positive results and because single annotations are not sufficient for an understanding of gene function. To increase the reliability of gene function annotation, multiple independent datasets need to be integrated. Here, we review the recent development of strategies for such integration and we argue that these will be important for a systems approach to modular biology.

Overexpression analysis of plant transcription factors

Transcription factors (TFs) play important roles in plant development and its response to the environment. A variety of reverse genetics tools have been developed to study TF function, the two most commonly used ones being knockout and overexpression. Because of the unique characteristics and modes of action of TFs, the overexpression strategy has been particularly effective in revealing TF function. Thus, a number of overexpression-based methodologies - constitutive expression, tissue-specific expression, chemically inducible expression and overexpression of modified TFs - have been developed and are used in the analysis of TF function.

New tools for functional mammalian cancer genetics

Knowledge of the function of individual genes that encode components of cell-signalling pathways is crucial to our understanding of normal growth control and its deregulation in cancer, but we have functional information for only 15% of human genes at present. Several new technologies have recently become available to identify gene function in mammalian cells using high-throughput genetic screens. These new tools will make it possible to identify new and innovative classes of anticancer drugs, including those that show synthetic lethal interactions with cancer-specific mutations.

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Functional genomics of cell morphology using RNA interference: pick your style, broad or deep

Several new studies have used RNA interference to screen for protein functions affecting cell shape, mitosis and cytokinesis of Drosophila cells in culture. One broad survey of nearly 1,000 proteins and three studies focused on cytoskeletal and motor proteins have identified key proteins essential for these processes in animal cells.

Rrp47p is an exosome-associated protein required for the 3' processing of stable RNAs

Related exosome complexes of 3'-->5' exonucleases are present in the nucleus and the cytoplasm. Purification of exosome complexes from whole-cell lysates identified a Mg(2+)-labile factor present in substoichiometric amounts. This protein was identified as the nuclear protein Yhr081p, the homologue of human C1D, which we have designated Rrp47p (for rRNA processing). Immunoprecipitation of epitope-tagged Rrp47p confirmed its interaction with the exosome and revealed its association with Rrp6p, a 3'-->5' exonuclease specific to the nuclear exosome fraction. Northern analyses demonstrated that Rrp47p is required for the exosome-dependent processing of rRNA and small nucleolar RNA (snoRNA) precursors. Rrp47p also participates in the 3' processing of U4 and U5 small nuclear RNAs (snRNAs). The defects in the processing of stable RNAs seen in rrp47-Delta strains closely resemble those of strains lacking Rrp6p. In contrast, Rrp47p is not required for the Rrp6p-dependent degradation of 3'-extended nuclear pre-mRNAs or the cytoplasmic 3'-->5' mRNA decay pathway. We propose that Rrp47p functions as a substrate-specific nuclear cofactor for exosome activity in the processing of stable RNAs.

Regulation of ENA1 Na(+)-ATPase gene expression by the Ppz1 protein phosphatase is mediated by the calcineurin pathway

Saccharomyces cerevisiae strains lacking the Ppz1 protein phosphatase are salt tolerant and display increased expression of the ENA1 Na(+)-ATPase gene, a major determinant for sodium extrusion, while cells devoid of the similar Ppz2 protein do not show these phenotypes. However, a ppz1 ppz2 mutant displays higher levels of ENA1 expression than the ppz1 strain. We show here that the increased activity of the ENA1 promoter in a ppz1 ppz2 mutant maps to two regions: one region located at -751 to -667, containing a calcineurin-dependent response element (CDRE), and one downstream region (-573 to -490) whose activity responds to intracellular alkalinization. In contrast, the increased ENA1 expression in a ppz1 mutant is mediated solely by an intact calcineurin/Crz1 signaling pathway, on the basis that (i) this effect maps to a single region that contains the CDRE and (ii) it is blocked by the calcineurin inhibitor FK506, as well as by deletion of the CNB1 or CRZ1 gene. The calcineurin dependence of the increased ENA1 expression of a ppz1 mutant would suggest that Ppz1 could negatively regulate calcineurin activity. In agreement with this notion, a ppz1 strain is calcium sensitive, and this mutation does not result in a decrease in the calcium hypertolerance of a cnb1 mutant. It has been shown that ENA1 can be induced by alkalinization of the medium and that a ppz1 ppz2 strain has a higher intracellular pH. However, we present several lines of evidence that show that the gene expression profile of a ppz1 mutant does not involve an alkalinization effect. In conclusion, we have identified a novel role for calcineurin, but not alkalinization, in the control of ENA1 expression in ppz1 mutants.

Widespread cytoplasmic mRNA transport in yeast: identification of 22 bud-localized transcripts using DNA microarray analysis

Cytoplasmic mRNA localization provides a means of generating cell asymmetry and segregating protein activity. Previous studies have identified two mRNAs that localize to the bud tips of the yeast Saccharomyces cerevisiae. To identify additional localized mRNAs, we immunoprecipitated the RNA transport components She2p, She3p, and Myo4p and performed DNA microarray analysis of their associated RNAs. A secondary screen, using a GFP-tagged RNA reporter assay, identified 22 mRNAs that are localized to bud tips. These messages encode a wide variety of proteins, including several involved in stress responses and cell wall maintenance. Many of these proteins are asymmetrically localized to buds. However, asymmetric localization also occurs in the absence of RNA transport, suggesting the existence of redundant protein localization mechanisms. In contrast to findings in metazoans, the untranslated regions are dispensable for mRNA localization in yeast. This study reveals an unanticipated widespread use of RNA transport in budding yeast.

A genomewide screen in Saccharomyces cerevisiae for genes that suppress the accumulation of mutations

A genomewide screen of a collection of 4,847 yeast gene deletion mutants was carried out to identify the genes required for suppressing mutations in the CAN1 forward-mutation assay. The primary screens and subsequent analysis allowed (i) identification of 18 known mutator mutants, providing a solid means for checking the efficiency of the screen, and (ii) identification of a number of genes not known previously to be involved in suppressing mutations. Among the previously uncharacterized mutation-suppressing genes were six genes of unknown function including four (CSM2, SHU2, SHU1, and YLR376c) encoding proteins that interact with each other and promote resistance to killing by methyl methanesulfonate, one gene (EGL1) previously identified as suppressing Ty1 mobility and recombination between repeated sequences, and one gene (YLR154c) that was not associated with any known processes. In addition, five genes (TSA1, SOD1, LYS7, SKN7, and YAP1) implicated in the oxidative-stress responses were found to play a significant role in mutation suppression. Furthermore, TSA1, which encodes thioredoxin peroxidase, was found to strongly suppress gross chromosomal rearrangements. These results provide a global view of the nonessential genes involved in preventing mutagenesis. Study of such genes should provide useful clues in identification of human genes potentially involved in cancer predisposition and in understanding their mechanisms of action.

Equilibrium distribution of mutators in the single fitness peak model

This Letter develops an analytically tractable model for determining the equilibrium distribution of mismatch repair deficient strains in unicellular populations. The approach is based on the single fitness peak model, which has been used in Eigen's quasispecies equations in order to understand various aspects of evolutionary dynamics. As with the quasispecies model, our model for mutator-nonmutator equilibrium undergoes a phase transition in the limit of infinite sequence length. This "repair catas-trophe" occurs at a critical repair error probability of epsilon(r)=L(via)/L, where L(via) denotes the length of the genome controlling viability, while L denotes the overall length of the genome. The repair catastrophe therefore occurs when the repair error probability exceeds the fraction of deleterious mutations. Our model also gives a quantitative estimate for the equilibrium fraction of mutators in Escherichia coli.

Tackling an essential problem in functional proteomics of Saccharomyces cerevisiae

Gene inactivation is the cornerstone of functional genetic analysis, but the analysis of essential genes requires conditional inactivation of the gene product. A new study has adapted a simple method for creating conditional alleles to allow large-scale analysis of essential genes in Saccharomyces cerevisiae and has identified a role in DNA replication for a newly identified protein complex.

Gene inactivation is the cornerstone of functional genetic analysis, but the analysis of essential genes requires conditional inactivation of the gene product. A new study has adapted a simple method for creating conditional alleles to allow large-scale analysis of essential genes in Saccharomyces cerevisiae and has identified a role in DNA replication for a newly identified protein complex.

Molecular characterization of Ypi1, a novel Saccharomyces cerevisiae type 1 protein phosphatase inhibitor

The Saccharomyces cerevisiae open reading frame YFR003c encodes a small (155-amino acid) hydrophilic protein that we identified as a novel, heat-stable inhibitor of type 1 protein phosphatase (Ypi1). Ypi1 interacts physically in vitro with both Glc7 and Ppz1 phosphatase catalytic subunits, as shown by pull-down assays. Ypi1 inhibits Glc7 but appears to be less effective toward Ppz1 phosphatase activity under the conditions tested. Ypi1 contains a 48RHNVRW53 sequence, which resembles the characteristic consensus PP1 phosphatase binding motif. A W53A mutation within this motif abolishes both binding to and inhibition of Glc7 and Ppz1 phosphatases. Deletion of YPI1 is lethal, suggesting a relevant role of the inhibitor in yeast physiology. Cells overexpressing Ypi1 display a number of phenotypes consistent with an inhibitory role of this protein on Glc7, such as decreased glycogen content and an increased growth defect in a slt2/mpk1 mitogen-activated protein kinase-deficient background. Taking together, these results define Ypi1 as the first inhibitory subunit of Glc7 identified in budding yeast.

Methylthioadenosine phosphorylase regulates ornithine decarboxylase by production of downstream metabolites

The gene encoding methylthioadenosine phosphorylase (MTAP), the initial enzyme in the methionine salvage pathway, is deleted in a variety of human tumors and acts as a tumor suppressor gene in cell culture (Christopher, S. A., Diegelman, P., Porter, C. W., and Kruger, W. D. (2002) Cancer Res. 62, 6639-6644). Overexpression of the polyamine biosynthetic enzyme ornithine decarboxylase (ODC) is frequently observed in tumors and has been shown to be tumorigenic in vitro and in vivo. In this paper, we demonstrate a novel regulatory pathway in which the methionine salvage pathway products inhibit ODC activity. We show that in Saccharomyces cerevisiae the MEU1 gene encodes MTAP and that Meu1delta cells have an 8-fold increase in ODC activity, resulting in large elevations in polyamine pools. Mutations in putative salvage pathway genes downstream of MTAP also cause elevated ODC activity and elevated polyamines. The addition of the penultimate salvage pathway compound 4-methylthio-2-oxobutanoic acid represses ODC levels in both MTAP-deleted yeast and human tumor cell lines, indicating that 4-methylthio-2-oxobutanoic acid acts as a negative regulator of polyamine biosynthesis. Expression of MTAP in MTAP-deleted MCF-7 breast adenocarcinoma cells results in a significant reduction of ODC activity and reduction in polyamine levels. Taken together, our results show that products of the methionine salvage pathway regulate polyamine biosynthesis and suggest that MTAP deletion may lead to ODC activation in human tumors.

Probing the yeast proteome for RNA-processing factors

A method has been developed to identify proteins required for the biogenesis of non-coding RNA in yeast, using a microarray to screen for aberrant patterns of RNA processing in mutant strains, and new proteins involved in the processing of ribosomal and non-coding RNAs have been found.

A method has been developed to identify proteins required for the biogenesis of non-coding RNA in yeast, using a microarray to screen for aberrant patterns of RNA processing in mutant strains, and new proteins involved in the processing of ribosomal and non-coding RNAs have been found.

A novel eukaryotic factor for cytosolic Fe-S cluster assembly

Iron regulatory protein 1 (IRP1) is regulated through the assembly/disassembly of a [4Fe-4S] cluster, which interconverts IRP1 with cytosolic aconitase. A genetic screen to isolate Saccharomyces cerevisiae strains bearing mutations in genes required for the conversion of IRP1 to c-aconitase led to the identification of a previously uncharacterized, essential gene, which we call CFD1 (cytosolic Fe-S cluster deficient). CFD1 encodes a highly conserved, putative P-loop ATPase. A non-lethal mutation of CFD1 (cfd1-1) reduced c-aconitase specific activity in IRP1-transformed yeast by >90%, although IRP1 in these cells could be readily converted to c-aconitase in vitro upon incubation with iron alone. IRP1-transformed cfd1-1 yeast lacked EPR-detectable Fe-S clusters in c-aconitase, pointing to a defect in Fe-S cluster assembly. The specific activity of another cytosolic Fe-S protein, Leu1p, was also inhibited by >90% in cfd1-1 yeast, whereas activity of mitochondrial Fe-S proteins was not inhibited. Consistent with a cytosolic site of activity, Cfd1p was localized in the cytoplasm. To our knowledge, Cfd1p is the first cytoplasmic Fe-S cluster assembly factor described in eukaryotes.

Genome-wide expression profiling of the response to polyene, pyrimidine, azole, and echinocandin antifungal agents in Saccharomyces cerevisiae

Antifungal compounds exert their activity through a variety of mechanisms, some of which are poorly understood. Novel approaches to characterize the mechanism of action of antifungal agents will be of great use in the antifungal drug development process. The aim of the present study was to investigate the changes in the gene expression profile of Saccharomyces cerevisiae following exposure to representatives of the four currently available classes of antifungal agents used in the management of systemic fungal infections. Microarray analysis indicated differential expression of 0.8, 4.1, 3.0, and 2.6% of the genes represented on the Affymetrix S98 yeast gene array in response to ketoconazole, amphotericin B, caspofungin, and 5-fluorocytosine (5-FC), respectively. Quantitative real time reverse transcriptase-PCR was used to confirm the microarray analyses. Genes responsive to ketoconazole, caspofungin, and 5-FC were indicative of the drug-specific effects. Ketoconazole exposure primarily affected genes involved in ergosterol biosynthesis and sterol uptake; caspofungin exposure affected genes involved in cell wall integrity; and 5-FC affected genes involved in DNA and protein synthesis, DNA damage repair, and cell cycle control. In contrast, amphotericin B elicited changes in gene expression reflecting cell stress, membrane reconstruction, transport, phosphate uptake, and cell wall integrity. Genes with the greatest specificity for a particular drug were grouped together as drug-specific genes, whereas genes with a lack of drug specificity were also identified. Taken together, these data shed new light on the mechanisms of action of these classes of antifungal agents and demonstrate the potential utility of gene expression profiling in antifungal drug development.

Chemical genomic profiling of biological networks using graph theory and combinations of small molecule perturbations

Genome-wide measurements of multiple experimental samples yield rich fingerprints for comparison and interpretation. Here, a two-dimensional matrix of the cellular effects of all possible pairwise combinations of 24 small molecules, each with a different structure and bioactivity, was used to profile otherwise isogenic deletion strains of the yeast Saccharomyces cerevisiae. Using principles from graph theory, we derived a discrete model of the data for each strain by encoding the information in the form of a binary adjacency matrix. This matrix was used to construct a graph composed of nodes representing small molecules and edges connecting combinations that inhibited cell cycle progression. Computation of a set of graph theoretic descriptors for each chemical genetic network provided a topological fingerprint that showed genotype-dependent fluctuations. Because the structure of the genetic network determines the structure of the chemical genetic network, multidimensional chemical genomic profiling can be used for the characterization of perturbations in biological networks or the networks themselves. This application of small molecules could be useful for discerning the molecular basis of highly complex biological phenotypes, including those involved in the susceptibility to or etiology of human disease.

Expression deconvolution: a reinterpretation of DNA microarray data reveals dynamic changes in cell populations

Cells grow in dynamically evolving populations, yet this aspect of experiments often goes unmeasured. A method is proposed for measuring the population dynamics of cells on the basis of their mRNA expression patterns. The population's expression pattern is modeled as the linear combination of mRNA expression from pure samples of cells, allowing reconstruction of the relative proportions of pure cell types in the population. Application of the method, termed expression deconvolution, to yeast grown under varying conditions reveals the population dynamics of the cells during the cell cycle, during the arrest of cells induced by DNA damage and the release of arrest in a cell cycle checkpoint mutant, during sporulation, and following environmental stress. Using expression deconvolution, cell cycle defects are detected and temporally ordered in 146 yeast deletion mutants; six of these defects are independently experimentally validated. Expression deconvolution allows a reinterpretation of the cell cycle dynamics underlying all previous microarray experiments and can be more generally applied to study most forms of cell population dynamics.

Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors

Natural genetic variation can cause significant differences in gene expression, but little is known about the polymorphisms that affect gene regulation. We analyzed regulatory variation in a cross between laboratory and wild strains of Saccharomyces cerevisiae. Clustering and linkage analysis defined groups of coregulated genes and the loci involved in their regulation. Most expression differences mapped to trans-acting loci. Positional cloning and functional assays showed that polymorphisms in GPA1 and AMN1 affect expression of genes involved in pheromone response and daughter cell separation, respectively. We also asked whether particular classes of genes were more likely to contain trans-regulatory polymorphisms. Notably, transcription factors showed no enrichment, and trans-regulatory variation seems to be broadly dispersed across classes of genes with different molecular functions.

Yeast functional analysis: identification of two essential genes involved in ER to Golgi trafficking

We screened for genes potentially involved in the secretory and vacuolar pathways a collection of 61 yeast strains, each bearing an essential orphan gene regulated by the tetO7-CYC1 promoter that can be down-regulated by doxycycline. After down-regulating the expression of these genes, we performed systematic Western blot analysis for markers of the secretory and vacuolar pathways that undergo post-translational modifications in their intracellular trafficking. Accumulation of protein precursors, revealed by Western immunoblot analysis, indicates defects in the secretory pathway or in associated biochemical modifications. After screening the whole collection, we identified two genes involved in ER to Golgi trafficking: RER2, a cis-prenyl transferase, and USE1, the function of which was unknown. We demonstrated that repression of USE1 also leads to BiP secretion, and therefore likely affects retrograde, in addition to anterograde, ER to Golgi trafficking. The collection also includes two essential genes involved in intracellular trafficking that were conveniently repressed without resulting growth or trafficking defects.

A chemical genomics approach to understanding drug action

The complete collection of yeast deletion strains represents a unique, living biological computer for understanding gene function. The molecular 'barcodes' present in each of the deletion strains allow a quantitative ranking of the importance of any gene under any experimental condition of choice. In this article, some of the recent results generated from experiments that exploit the yeast deletion collection to understand mechanisms of drug action are discussed.

Large-scale essential gene identification in Candida albicans and applications to antifungal drug discovery

Candida albicans is the primary fungal pathogen of humans. Despite the need for novel drugs to combat fungal infections [Sobel, J.D. (2000) Clin Infectious Dis 30: 652], antifungal drug discovery is currently limited by both the availability of suitable drug targets and assays to screen corresponding targets. A functional genomics approach based on the diploid C. albicans genome sequence, termed GRACETM (gene replacement and conditional expression), was used to assess gene essentiality through a combination of gene replacement and conditional gene expression. In a systematic application of this approach, we identify 567 essential genes in C. albicans. Interestingly, evaluating the conditional phenotype of all identifiable C. albicans homologues of the Saccharomyces cerevisiae essential gene set [Giaever, G., Chu, A.M., Ni, L., Connelly, C., Riles, L., Veronneau, S., et al. (2002) Nature 418: 387-391] by GRACE revealed only 61% to be essential in C. albicans, emphasizing the importance of performing such studies directly within the pathogen. Construction of this conditional mutant strain collection facilitates large-scale examination of terminal phenotypes of essential genes. This information enables preferred drug targets to be selected from the C. albicans essential gene set by phenotypic information derived both in vitro, such as cidal versus static terminal phenotypes, as well as in vivo through virulence studies using conditional strains in an animal model of infection. In addition, the combination of phenotypic and bioinformatic analyses further improves drug target selection from the C. albicans essential gene set, and their respective conditional mutant strains may be directly used as sensitive whole-cell assays for drug screening.

Ribosome assembly in eukaryotes

Ribosome synthesis is a highly complex and coordinated process that occurs not only in the nucleolus but also in the nucleoplasm and the cytoplasm of eukaryotic cells. Based on the protein composition of several ribosomal subunit precursors recently characterized in yeast, a total of more than 170 factors are predicted to participate in ribosome biogenesis and the list is still growing. So far the majority of ribosomal factors have been implicated in RNA maturation (nucleotide modification and processing). Recent advances gave insight into the process of ribosome export and assembly. Proteomic approaches have provided the first indications for a ribosome assembly pathway in eukaryotes and confirmed the dynamic character of the whole process.

Functional distinctions between IMP dehydrogenase genes in providing mycophenolate resistance and guanine prototrophy to yeast

IMP dehydrogenase (IMPDH) catalyzes the rate-limiting step in the de novo synthesis of GTP. Yeast with mutations in the transcription elongation machinery are sensitive to inhibitors of this enzyme such as 6-azauracil and mycophenolic acid, at least partly because of their inability to transcriptionally induce IMPDH. To understand the molecular basis of this drug-sensitive phenotype, we have dissected the expression and function of a four-gene family in yeast called IMD1 through IMD4. We show here that these family members are distinct, despite a high degree of amino acid identity between the proteins they encode. Extrachromosomal copies of IMD1, IMD3, or IMD4 could not rescue the drug-sensitive phenotype of IMD2 deletants. When overexpressed, IMD3 or IMD4 weakly compensated for deletion of IMD2. IMD1 is transcriptionally silent and bears critical amino acid substitutions compared with IMD2 that destroy its function, offering strong evidence that it is a pseudogene. The simultaneous deletion of all four IMD genes was lethal unless growth media were supplemented with guanine. This suggests that there are no other essential functions of the IMPDH homologs aside from IMP dehydrogenase activity. Although neither IMD3 nor IMD4 could confer drug resistance to cells lacking IMD2, either alone was sufficient to confer guanine prototrophy. The special function of IMD2 was provided by its ability to be transcriptionally induced and the probable intrinsic drug resistance of its enzymatic activity.

Early phosphorylation events in biotic stress

Several kinases are activated during the defence response following microbial elicitation. While studying the regulation of these kinases in greater detail, it has become clear that the means by which phosphorylation events transmit specific information to the cell is highly complex. To gain a better understanding of the molecular events leading to a response, it will be increasingly important to identify not only the protein targets of phosphorylation but also the specific sites of phosphorylation. New developments in peptide-based phosphoproteome analysis appear to hold the promise of achieving these goals at the whole-cell level.