Microbes tailor their growth rate to nutrient availability. Here, we measured, using liquid chromatography-mass spectrometry, >100 intracellular metabolites in steady-state cultures of Saccharomyces cerevisiae growing at five different rates and in each of five different limiting nutrients. In contrast to gene transcripts, where approximately 25% correlated with growth rate irrespective of the nature of the limiting nutrient, metabolite concentrations were highly sensitive to the limiting nutrient's identity. Nitrogen (ammonium) and carbon (glucose) limitation were characterized by low intracellular amino acid and high nucleotide levels, whereas phosphorus (phosphate) limitation resulted in the converse. Low adenylate energy charge was found selectively in phosphorus limitation, suggesting the energy charge may actually measure phosphorus availability. Particularly strong concentration responses occurred in metabolites closely linked to the limiting nutrient, e.g., glutamine in nitrogen limitation, ATP in phosphorus limitation, and pyruvate in carbon limitation. A simple but physically realistic model involving the availability of these metabolites was adequate to account for cellular growth rate. The complete data can be accessed at the interactive website http://growthrate.princeton.edu/metabolome.

Maintaining balanced growth in a changing environment is a fundamental systems-level challenge for cellular physiology, particularly in microorganisms. While the complete set of regulatory and functional pathways supporting growth and cellular proliferation are not yet known, portions of them are well understood. In particular, cellular proliferation is governed by mechanisms that are highly conserved from unicellular to multicellular organisms, and the disruption of these processes in metazoans is a major factor in the development of cancer. In this paper, we develop statistical methodology to identify quantitative aspects of the regulatory mechanisms underlying cellular proliferation in Saccharomyces cerevisiae. We find that the expression levels of a small set of genes can be exploited to predict the instantaneous growth rate of any cellular culture with high accuracy. The predictions obtained in this fashion are robust to changing biological conditions, experimental methods, and technological platforms. The proposed model is also effective in predicting growth rates for the related yeast Saccharomyces bayanus and the highly diverged yeast Schizosaccharomyces pombe, suggesting that the underlying regulatory signature is conserved across a wide range of unicellular evolution. We investigate the biological significance of the gene expression signature that the predictions are based upon from multiple perspectives: by perturbing the regulatory network through the Ras/PKA pathway, observing strong upregulation of growth rate even in the absence of appropriate nutrients, and discovering putative transcription factor binding sites, observing enrichment in growth-correlated genes. More broadly, the proposed methodology enables biological insights about growth at an instantaneous time scale, inaccessible by direct experimental methods. Data and tools enabling others to apply our methods are available at http://function.princeton.edu/growthrate.

In yeast, the pathways of sulfur assimilation are combinatorially controlled by five transcriptional regulators (three DNA-binding proteins [Met31p, Met32p, and Cbf1p], an activator [Met4p], and a cofactor [Met28p]) and a ubiquitin ligase subunit (Met30p). This regulatory system exerts combinatorial control not only over sulfur assimilation and methionine biosynthesis, but also on many other physiological functions in the cell. Recently we characterized a gene induction system that, upon the addition of an inducer, results in near-immediate transcription of a gene of interest under physiological conditions. We used this to perturb levels of single transcription factors during steady-state growth in chemostats, which facilitated distinction of direct from indirect effects of individual factors dynamically through quantification of the subsequent changes in genome-wide patterns of gene expression. We were able to show directly that Cbf1p acts sometimes as a repressor and sometimes as an activator. We also found circumstances in which Met31p/Met32p function as repressors, as well as those in which they function as activators. We elucidated and numerically modeled feedback relationships among the regulators, notably feedforward regulation of Met32p (but not Met31p) by Met4p that generates dynamic differences in abundance that can account for the differences in function of these two proteins despite their identical binding sites.

High concentrations of ammonium at physiological concentrations of potassium are toxic for the standard laboratory strain of Saccharomyces cerevisiae In the original description of this metabolic phenotype, we focused on the standard laboratory strains of Saccharomyces In this study, we screened a large collection of S. cerevisiae natural isolates and identified one strain that is resistant to high concentrations of ammonium. This strain, K12, was isolated in sake breweries. When the K12 strain was crossed to the standard laboratory strain (FY4), the resulting tetrads displayed 2:2 segregation of the resistance phenotype, suggesting a single gene trait. Using a bulk segregant analysis strategy, we mapped this trait to a 150-kb region on chromosome X containing the TRK1 gene. This gene encodes a transporter required for high-affinity potassium transport in S. cerevisiae Data from reciprocal hemizygosity experiments with TRK1 deletion strains in K12 and BY backgrounds, as well as analysis of the deletion of this gene in the K12 strain, demonstrate that the K12 allele of TRK1 is responsible for ammonium toxicity resistance. Furthermore, we determined the minimal amount of potassium required for both the K12 and laboratory strain needed for growth. These results demonstrate that the gene encoded by the K12 allele of TRK1 has a greater affinity for potassium than the standard allele of TRK1 found in Saccharomyces strains. We hypothesize that this greater-affinity allele of the potassium transporter reduces the flux of ammonium into the yeast cells under conditions of ammonium toxicity. These findings further refine our understanding of ammonium toxicity in yeast and provide an example of using natural variation to understand cellular processes.

Sequencing and annotation of the entire Saccharomyces cerevisiae genome has made it possible to gain a genome-wide perspective on yeast genes and gene products. To make this information available on an ongoing basis, the Saccharomyces Genome Database (SGD) (http://www.yeastgenome.org/) has created the Genome Snapshot (http://db.yeastgenome.org/cgi-bin/genomeSnapShot.pl). The Genome Snapshot summarizes the current state of knowledge about the genes and chromosomal features of S.cerevisiae. The information is organized into two categories: (i) number of each type of chromosomal feature annotated in the genome and (ii) number and distribution of genes annotated to Gene Ontology terms. Detailed lists are accessible through SGD's Advanced Search tool (http://db.yeastgenome.org/cgi-bin/search/featureSearch), and all the data presented on this page are available from the SGD ftp site (ftp://ftp.yeastgenome.org/yeast/).

Genome rearrangements, especially amplifications and deletions, have regularly been observed as responses to sustained application of the same strong selective pressure in microbial populations growing in continuous culture. We studied eight strains of budding yeast (Saccharomyces cerevisiae) isolated after 100-500 generations of growth in glucose-limited chemostats. Changes in DNA copy number were assessed at single-gene resolution by using DNA microarray-based comparative genomic hybridization. Six of these evolved strains were aneuploid as the result of gross chromosomal rearrangements. Most of the aneuploid regions were the result of translocations, including three instances of a shared breakpoint on chromosome 14 immediately adjacent to CIT1, which encodes the citrate synthase that performs a key regulated step in the tricarboxylic acid cycle. Three strains had amplifications in a region of chromosome 4 that includes the high-affinity hexose transporters; one of these also had the aforementioned chromosome 14 break. Three strains had extensive overlapping deletions of the right arm of chromosome 15. Further analysis showed that each of these genome rearrangements was bounded by transposon-related sequences at the breakpoints. The observation of repeated, independent, but nevertheless very similar, chromosomal rearrangements in response to persistent selection of growing cells parallels the genome rearrangements that characteristically accompany tumor progression.

A yeast strain lacking Met4p, the primary transcriptional regulator of the sulfur assimilation pathway, cannot synthesize methionine. This apparently simple auxotroph did not grow well in rich media containing excess methionine, forming small colonies on yeast extract/peptone/dextrose plates. Faster-growing large colonies were abundant when overnight cultures were plated, suggesting that spontaneous suppressors of the growth defect arise with high frequency. To identify the suppressor mutations, we used genome-wide single-nucleotide polymorphism and standard genetic analyses. The most common suppressors were loss-of-function mutations in OPI1, encoding a transcriptional repressor of phospholipid metabolism. Using a new system that allows rapid and specific degradation of Met4p, we could study the dynamic expression of all genes following loss of Met4p. Experiments using this system with and without Opi1p showed that Met4 activates and Opi1p represses genes that maintain levels of S-adenosylmethionine (SAM), the substrate for most methyltransferase reactions. Cells lacking Met4p grow normally when either SAM is added to the media or one of the SAM synthetase genes is overexpressed. SAM is used as a methyl donor in three Opi1p-regulated reactions to create the abundant membrane phospholipid, phosphatidylcholine. Our results show that rapidly growing cells require significant methylation, likely for the biosynthesis of phospholipids.

The genome of budding yeast (Saccharomyces cerevisiae) contains approximately 5800 protein-encoding genes, the majority of which are associated with some known biological function. Yet the extent of amino acid sequence conservation of these genes over all phyla has only been partially examined. Here we provide a more comprehensive overview and visualization of the conservation of yeast genes and a means for browsing and exploring the data in detail, down to the individual yeast gene, at http://yeast-phylogroups.princeton.edu. We used data from the OrthoMCL database, which has defined orthologs from approximately 150 completely sequenced genomes, including diverse representatives of the archeal, bacterial, and eukaryotic domains. By clustering genes based on similar patterns of conservation, we organized and visualized all the protein-encoding genes in yeast as a single heat map. Most genes fall into one of eight major clusters, called "phylogroups." Gene ontology analysis of the phylogroups revealed that they were associated with specific, distinct trends in gene function, generalizations likely to be of interest to a wide range of biologists.

A general method for the dynamic control of single gene expression in eukaryotes, with no off-target effects, is a long-sought tool for molecular and systems biologists. We engineered two artificial transcription factors (ATFs) that contain Cys(2)His(2) zinc-finger DNA-binding domains of either the mouse transcription factor Zif268 (9 bp of specificity) or a rationally designed array of four zinc fingers (12 bp of specificity). These domains were expressed as fusions to the human estrogen receptor and VP16 activation domain. The ATFs can rapidly induce a single gene driven by a synthetic promoter in response to introduction of an otherwise inert hormone with no detectable off-target effects. In the absence of inducer, the synthetic promoter is inactive and the regulated gene product is not detected. Following addition of inducer, transcripts are induced >50-fold within 15 min. We present a quantitative characterization of these ATFs and provide constructs for making their implementation straightforward. These new tools allow for the elucidation of regulatory network elements dynamically, which we demonstrate with a major metabolic regulator, Gcn4p.

The experimental evolution of laboratory populations of microbes provides an opportunity to observe the evolutionary dynamics of adaptation in real time. Until very recently, however, such studies have been limited by our inability to systematically find mutations in evolved organisms. We overcome this limitation by using a variety of DNA microarray-based techniques to characterize genetic changes -- including point mutations, structural changes, and insertion variation -- that resulted from the experimental adaptation of 24 haploid and diploid cultures of Saccharomyces cerevisiae to growth in either glucose, sulfate, or phosphate-limited chemostats for approximately 200 generations. We identified frequent genomic amplifications and rearrangements as well as novel retrotransposition events associated with adaptation. Global nucleotide variation detection in ten clonal isolates identified 32 point mutations. On the basis of mutation frequencies, we infer that these mutations and the subsequent dynamics of adaptation are determined by the batch phase of growth prior to initiation of the continuous phase in the chemostat. We relate these genotypic changes to phenotypic outcomes, namely global patterns of gene expression, and to increases in fitness by 5-50%. We found that the spectrum of available mutations in glucose- or phosphate-limited environments combined with the batch phase population dynamics early in our experiments allowed several distinct genotypic and phenotypic evolutionary pathways in response to these nutrient limitations. By contrast, sulfate-limited populations were much more constrained in both genotypic and phenotypic outcomes. Thus, the reproducibility of evolution varies with specific selective pressures, reflecting the constraints inherent in the system-level organization of metabolic processes in the cell. We were able to relate some of the observed adaptive mutations (e.g., transporter gene amplifications) to known features of the relevant metabolic pathways, but many of the mutations pointed to genes not previously associated with the relevant physiology. Thus, in addition to answering basic mechanistic questions about evolutionary mechanisms, our work suggests that experimental evolution can also shed light on the function and regulation of individual metabolic pathways.

Obtaining reliable and reproducible two-color microarray gene expression data is critically important for understanding the biological significance of perturbations made on a cellular system. Microarray design, RNA preparation and labeling, hybridization conditions and data acquisition and analysis are variables difficult to simultaneously control. A useful tool for monitoring and controlling intra- and inter-experimental variation is Universal Reference RNA (URR), developed with the goal of providing hybridization signal at each microarray probe location (spot). Measuring signal at each spot as the ratio of experimental RNA to reference RNA targets, rather than relying on absolute signal intensity, decreases variability by normalizing signal output in any two-color hybridization experiment.

Loss-of-function alterations of INK4A are commonly observed in lymphoid malignancies, but are consistently absent in pre-B cell leukemias induced by the chimeric oncoprotein E2a-Pbx1 created by t(1;19) chromosomal translocations. We report here that experimental induction of E2a-Pbx1 enhances expression of BMI-1, a lymphoid oncogene whose product functions as a transcriptional repressor of the INK4A-ARF tumor suppressor locus. Bmi-1-deficient hematopoietic progenitors are resistant to transformation by E2a-Pbx1; however, the requirement for Bmi-1 is alleviated in cells deficient for both Bmi-1 and INK4A-ARF. Furthermore, the adverse effects of E2a-Pbx1 on pre-B cell survival and differentiation are partially bypassed by forced expression of p16(Ink4a). These results link E2a-Pbx1 with Bmi-1 on an oncogenic pathway that is likely to play a role in the pathogenesis of human lymphoid leukemias through downregulation of the INK4A-ARF gene.

DNA microarray analysis of gene expression in steady-state chemostat cultures limited for potassium revealed a surprising connection between potassium and ammonium: potassium limits growth only when ammonium is the nitrogen source. Under potassium limitation, ammonium appears to be toxic for Saccharomyces cerevisiae. This ammonium toxicity, which appears to occur by leakage of ammonium through potassium channels, is recapitulated under high-potassium conditions by over-expression of ammonium transporters. Although ammonium toxicity is well established in metazoans, it has never been reported for yeast. To characterize the response to ammonium toxicity, we examined the filtrates of these cultures for compounds whose excretion might serve to detoxify the ammonium (such as urea in mammals). Using liquid chromatography-tandem mass spectrometry to assay for a wide array of metabolites, we detected excreted amino acids. The amounts of amino acids excreted increased in relation to the severity of growth impairment by ammonium, suggesting that amino acid excretion is used by yeast for ammonium detoxification.

Cancer invasion and metastasis have been likened to wound healing gone awry. Despite parallels in cellular behavior between cancer progression and wound healing, the molecular relationships between these two processes and their prognostic implications are unclear. In this study, based on gene expression profiles of fibroblasts from ten anatomic sites, we identify a stereotyped gene expression program in response to serum exposure that appears to reflect the multifaceted role of fibroblasts in wound healing. The genes comprising this fibroblast common serum response are coordinately regulated in many human tumors, allowing us to identify tumors with gene expression signatures suggestive of active wounds. Genes induced in the fibroblast serum-response program are expressed in tumors by the tumor cells themselves, by tumor-associated fibroblasts, or both. The molecular features that define this wound-like phenotype are evident at an early clinical stage, persist during treatment, and predict increased risk of metastasis and death in breast, lung, and gastric carcinomas. Thus, the transcriptional signature of the response of fibroblasts to serum provides a possible link between cancer progression and wound healing, as well as a powerful predictor of the clinical course in several common carcinomas.

We studied the steady-state responses to changes in growth rate of yeast when ethanol is the sole source of carbon and energy. Analysis of these data, together with data from studies where glucose was the carbon source, allowed us to distinguish a "universal" growth rate response (GRR) common to all media studied from a GRR specific to the carbon source. Genes with positive universal GRR include ribosomal, translation, and mitochondrial genes, and those with negative GRR include autophagy, vacuolar, and stress response genes. The carbon source-specific GRR genes control mitochondrial function, peroxisomes, and synthesis of vitamins and cofactors, suggesting this response may reflect the intensity of oxidative metabolism. All genes with universal GRR, which comprise 25% of the genome, are expressed periodically in the yeast metabolic cycle (YMC). We propose that the universal GRR may be accounted for by changes in the relative durations of the YMC phases. This idea is supported by oxygen consumption data from metabolically synchronized cultures with doubling times ranging from 5 to 14 h. We found that the high oxygen consumption phase of the YMC can coincide exactly with the S phase of the cell division cycle, suggesting that oxidative metabolism and DNA replication are not incompatible.

Cytoprotective functions of a 20S proteasome activator were investigated. Saccharomyces cerevisiae Blm10 and human 20S proteasome activator 200 (PA200) are homologs. Comparative genome-wide analyses of untreated diploid cells lacking Blm10 and growing at steady state at defined growth rates revealed downregulation of numerous genes required for accurate chromosome structure, assembly and repair, and upregulation of a specific subset of genes encoding protein-folding chaperones. Blm10 loss or truncation of the Ubp3/Blm3 deubiquitinating enzyme caused massive chromosomal damage and cell death in homozygous diploids after phleomycin treatments, indicating that Blm10 and Ubp3/Blm3 function to stabilize the genome and protect against cell death. Diploids lacking Blm10 also were sensitized to doxorubicin, hydroxyurea, 5-fluorouracil, rapamycin, hydrogen peroxide, methyl methanesulfonate, and calcofluor. Fluorescently tagged Blm10 localized in nuclei, with enhanced fluorescence after DNA replication. After DNA damage that caused a classic G2/M arrest, fluorescence remained diffuse, with evidence of nuclear fragmentation in some cells. Protective functions of Blm10 did not require the carboxyl-terminal region that makes close contact with 20S proteasomes, indicating that protection does not require this contact or the truncated Blm10 can interact with the proteasome apart from this region. Without its carboxyl-terminus, Blm10((-339aa)) localized to nuclei in untreated, nonproliferating (G(0)) cells, but not during G(1) S, G(2), and M. The results indicate Blm10 functions in protective mechanisms that include the machinery that assures proper assembly of chromosomes. These essential guardian functions have implications for ubiquitin-independent targeting in anticancer therapy. Targeting Blm10/PA200 together with one or more of the upregulated chaperones or a conventional treatment could be efficacious.

The subcellular localization of proteins is critical to their biological roles. Moreover, whether a protein is membrane-bound, secreted, or intracellular affects the usefulness of, and the strategies for, using a protein as a diagnostic marker or a target for therapy. We employed a rapid and efficient experimental approach to classify thousands of human gene products as either "membrane-associated/secreted" (MS) or "cytosolic/nuclear" (CN). Using subcellular fractionation methods, we separated mRNAs associated with membranes from those associated with the soluble cytosolic fraction and analyzed these two pools by comparative hybridization to DNA microarrays. Analysis of 11 different human cell lines, representing lymphoid, myeloid, breast, ovarian, hepatic, colon, and prostate tissues, identified more than 5,000 previously uncharacterized MS and more than 6,400 putative CN genes at high confidence levels. The experimentally determined localizations correlated well with in silico predictions of signal peptides and transmembrane domains, but also significantly increased the number of human genes that could be cataloged as encoding either MS or CN proteins. Using gene expression data from a variety of primary human malignancies and normal tissues, we rationally identified hundreds of MS gene products that are significantly overexpressed in tumors compared to normal tissues and thus represent candidates for serum diagnostic tests or monoclonal antibody-based therapies. Finally, we used the catalog of CN gene products to generate sets of candidate markers of organ-specific tissue injury. The large-scale annotation of subcellular localization reported here will serve as a reference database and will aid in the rational design of diagnostic tests and molecular therapies for diverse diseases.

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org/) collects and organizes biological information about the chromosomal features and gene products of the budding yeast Saccharomyces cerevisiae. Although published data from traditional experimental methods are the primary sources of evidence supporting Gene Ontology (GO) annotations for a gene product, high-throughput experiments and computational predictions can also provide valuable insights in the absence of an extensive body of literature. Therefore, GO annotations available at SGD now include high-throughput data as well as computational predictions provided by the GO Annotation Project (GOA UniProt; http://www.ebi.ac.uk/GOA/). Because the annotation method used to assign GO annotations varies by data source, GO resources at SGD have been modified to distinguish data sources and annotation methods. In addition to providing information for genes that have not been experimentally characterized, GO annotations from independent sources can be compared to those made by SGD to help keep the literature-based GO annotations current.

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

Signal transduction pathways guided by cellular receptors commonly exhibit low-level constitutive signaling in a continuous, ligand-independent manner. The dynamic equilibrium of positive and negative regulators establishes such a tonic signal. Ligand-independent signaling by the precursors of mature antigen receptors regulates development of B and T lymphocytes. Here we describe a basal signal that controls gene expression profiles in the Jurkat T cell line and mouse thymocytes. Using DNA microarrays and Northern blots to analyze unstimulated cells, we demonstrate that expression of a cluster of genes, including RAG-1 and RAG-2, is repressed by constitutive signals requiring the adapter molecules LAT and SLP-76. This TCR-like pathway results in constitutive low-level activity of Erk and Abl kinases. Inhibition of Abl by the drug STI-571 or inhibition of signaling events upstream of Erk increases RAG-1 expression. Our data suggest that physiologic gene expression programs depend upon tonic activity of signaling pathways independent of receptor ligation.