The Ume6 regulon coordinates metabolic and meiotic gene expression in yeast

A de novo originated gene depresses budding yeast mating pathway and is repressed by the protein encoded by its antisense strand

Recent transcription profiling studies have revealed an unexpectedly large proportion of antisense transcripts in eukaryotic genomes. These antisense genes seem to regulate gene expression by interacting with sense genes. Previous studies have focused on the non-coding antisense genes, but the possible regulatory role of the antisense protein is poorly understood. In this study, we found that a protein encoded by the antisense gene ADF1 acts as a transcription suppressor, regulating the expression of sense gene MDF1 in Saccharomyces cerevisiae. Based on the evolutionary, genetic, cytological and biochemical evidence, we show that the protein-coding sense gene MDF1 most likely originated de novo from a previously non-coding sequence and can significantly suppress the mating efficiency of baker's yeast in rich medium by binding MATalpha2 and thus promote vegetative growth. These results shed new light on several important issues, including a new sense-antisense interaction mechanism, the de novo origination of a functional gene, and the regulation of yeast mating pathway.

Ume6 is required for the MATa/MATalpha cellular identity and transcriptional silencing in Kluyveromyces lactis

To explore the similarities and differences of regulatory circuits among budding yeasts, we characterized the role of the unscheduled meiotic gene expression 6 (UME6) gene in Kluyveromyces lactis. We found that Ume6 was required for transcriptional silencing of the cryptic mating-type loci HMLalpha and HMRa. Chromatin immunoprecipitation (ChIP) suggested that Ume6 acted directly by binding the cis-regulatory silencers of these loci. Unexpectedly, a MATa ume6 strain was mating proficient, whereas a MATalpha ume6 strain was sterile. This observation was explained by the fact that ume6 derepressed HMLalpha2 only weakly, but derepressed HMRa1 strongly. Consistently, two a/alpha-repressed genes (MTS1 and STE4) were repressed in the MATalpha ume6 strain, but were expressed in the MATa ume6 strain. Surprisingly, ume6 partially suppressed the mating defect of a MATa sir2 strain. MTS1 and STE4 were repressed in the MATa sir2 ume6 double-mutant strain, indicating that the suppression acted downstream of the a1/alpha2-repressor. We show that both STE12 and the MATa2/HMRa2 genes were overexpressed in the MATa sir2 ume6 strain. Consistent with the idea that this deregulation suppressed the mating defect, ectopic overexpression of Ste12 and a2 in a MATa sir2 strain resulted in efficient mating. In addition, Ume6 served as a block to polyploidy, since ume6/ume6 diploids mated as pseudo a-strains. Finally, Ume6 was required for repression of three meiotic genes, independently of the Rpd3 and Sin3 corepressors.

Transcriptional responses of candida albicans to epithelial and endothelial cells

Candida albicans interacts with oral epithelial cells during oropharyngeal candidiasis and with vascular endothelial cells when it disseminates hematogenously. We set out to identify C. albicans genes that govern interactions with these host cells in vitro. The transcriptional response of C. albicans to the FaDu oral epithelial cell line and primary endothelial cells was determined by microarray analysis. Contact with epithelial cells caused a decrease in transcript levels of genes related to protein synthesis and adhesion, whereas contact with endothelial cells did not significantly influence any specific functional category of genes. Many genes whose transcripts were increased in response to either host cell had not been previously characterized. We constructed mutants with homozygous insertions in 22 of these uncharacterized genes to investigate their function during host-pathogen interaction. By this approach, we found that YCK2, VPS51, and UEC1 are required for C. albicans to cause normal damage to epithelial cells and resist antimicrobial peptides. YCK2 is also necessary for maintenance of cell polarity. VPS51 is necessary for normal vacuole formation, resistance to multiple stressors, and induction of maximal endothelial cell damage. UEC1 encodes a unique protein that is required for resistance to cell membrane stress. Therefore, some C. albicans genes whose transcripts are increased upon contact with epithelial or endothelial cells are required for the organism to damage these cells and withstand the stresses that it likely encounters during growth in the oropharynx and bloodstream.

A genome-wide screen for regulators of TORC1 in response to amino acid starvation reveals a conserved Npr2/3 complex

TORC1 is a central regulator of cell growth in response to amino acid availability, yet little is known about how it is regulated. Here, we performed a reverse genetic screen in yeast for genes necessary to inactivate TORC1. The screen consisted of monitoring the expression of a TORC1 sensitive GFP-based transcriptional reporter in all yeast deletion strains using flow cytometry. We find that in response to amino acid starvation, but not to carbon starvation or rapamycin treatment, cells lacking NPR2 and NPR3 fail to fully (1) activate transcription factors Gln3/Gat1, (2) dephosphorylate TORC1 effector Npr1, and (3) repress ribosomal protein gene expression. Both mutants show proliferation defects only in media containing a low quality nitrogen source, such as proline or ammonia, whereas no defects are evident when cells are grown in the presence of glutamine or peptone mixture. Proliferation defects in npr2Δ and npr3Δ cells can be completely rescued by artificially inhibiting TORC1 by rapamycin, demonstrating that overactive TORC1 in both strains prevents their ability to adapt to an environment containing a low quality nitrogen source. A biochemical purification of each demonstrates that Npr2 and Npr3 form a heterodimer, and this interaction is evolutionarily conserved since the human homologs of NPR2 and NPR3 (NPRL2 and NPRL3, respectively) also co-immunoprecipitate. We conclude that, in yeast, the Npr2/3 complex mediates an amino acid starvation signal to TORC1.

Before a eukaryotic cell commits to cell division it must be large enough so that both daughter cells would be of viable size. The control of cell size is largely mediated by nutritional input signals via an evolutionarily conserved protein complex termed TORC1. In particular, TORC1 has been shown to sense the level of amino acids and its activity correlates with the level of amino acids present in the media. Yet, it is largely unknown how TORC1 senses amino acids. Here we demonstrate that the evolutionarily conserved Npr2/3 complex mediates the amino acid scarcity signal to TORC1. Cells lacking NPR2 and NPR3 genes fail to inactivate TORC1 when amino acids are scarce. Overactive TORC1 prevents these cells from adapting to an amino acid scarce environment, and, as a result, these cells are unable to proliferate in media that is not rich in amino acids. Artificially inhibiting TORC1 with rapamycin can completely rescue these defects. These results provide insight into how cells sense amino acid deficiency. Moreover, as deletions of NPR2 have been implicated in tumor growth, these results offer a fertile ground to study the role overactive TORC1 might play in those cancers.

Untargeted tail acetylation of histones in chromatin: lessons from yeast

Dynamic acetylation of lysine residues in the amino-terminal tails of the core histones is functionally important for the regulation of diverse DNA-dependent processes in the nucleus, including replication, transcription, and DNA repair. The targeted and untargeted activities of histone lysine acetylases (KATs) and deacetylases (HDACs) both contribute to the dynamics of chromatin acetylation. While the mechanisms and functional consequences of targeted on histone acetylation are well understood, relatively little is known about untargeted histone acetylation. Here, we review the current understanding of the mechanisms by which untargeted KAT and HDAC activities modulate the acetylation state of nucleosomal histones, focusing on results obtained for H3 and H4 in budding yeast. We also highlight unresolved problems in this area, including the question of how a particular steady-state level of untargeted acetylation is set in the absence of cis-dependent mechanisms that instruct the activity of KATs and HDACs.

The Rpd3/HDAC complex is present at the URS1 cis-element with hyperacetylated histone H3

In eukaryotes, the hypoacetylated state of histone N-terminal lysines at many gene-promoters, which is created by histone deacetylases (HDACs), is changed to the hyperacetylated state by the function of histone acetyltransferases (HATs) upon transcription activation. Although much insight has been obtained to date as to how modification of the histone tail regulates gene expression, little is known about how the transition between the unmodified and modified states takes place. In Saccharomyces cerevisiae, the HDAC complex containing Rpd3 (Rpd3L) represses the transcription of several sets of genes through the URS1 cis-element. We found that the histone H3 acetylation level at the URS1 of seven genes (INO1, CAT2, ACS1, YAT1, RIM4, CRC1, and SIP4) was elevated in the presence of Rpd3/HDAC in growth in acetate-containing medium (YPA), suggesting that a mechanism that regulates HDAC activity is present in this organism. The biological significance of this phenomenon is discussed below.

Superoxide anions regulate TORC1 and its ability to bind Fpr1:rapamycin complex

The small natural product rapamycin, when bound to FKBP12, is a potent inhibitor of an evolutionarily conserved Target of Rapamycin Complex 1 (TORC1), which plays a central role in mediating cellular response to nutrient availability. Given the prominent role of TORC1 in cell growth and proliferation, clinical trials have explored the possibility of using rapamycin as an anticancer agent. Unfortunately, the percentage of patients responding favorably has been low, intensifying the need to find biomarkers able to predict rapamycin sensitivity or resistance. In this study, we elucidate the molecular mechanism underlying partial rapamycin resistance in yeast. Using the yeast deletion collection, we identified 15 deletion strains leading to partial rapamycin resistance. Among these were Cu/Zn-superoxide dismutase Sod1, copper transporter Ctr1, and copper chaperone Lys7, suggesting a role for oxidative stress in rapamycin resistance. Further analysis revealed that all 15 strains exhibit elevated levels of superoxide anions, and we show that elevated levels of reactive oxygen species specifically modify TORC1 such that it is no longer able to fully bind FKBP12:rapamycin. Therefore, elevated oxidative stress modifies TORC1 and prevents its binding to the FKBP12:rapamycin complex, ultimately leading to rapamycin resistance. These results warrant an examination into whether similar reasons explain rapamycin resistance observed in various clinical samples.

Effect of starvation on transcriptomes of brain and liver in adult female zebrafish (Danio rerio)

We used microarray and quantitative real-time PCR (qRT-PCR) analyses in adult female zebrafish (Danio rerio) to identify metabolic pathways regulated by starvation in the liver and brain. The transcriptome of whole zebrafish brain showed little response to 21 days of starvation. Only agouti-related protein 1 (agrp1) significantly responded, with increased expression in brains of starved fish. In contrast, a 21-day period of starvation significantly downregulated 466 and upregulated 108 transcripts in the liver, indicating an overall decrease in metabolic activity, reduced lipid metabolism, protein biosynthesis, proteolysis, and cellular respiration, and increased gluconeogenesis. Starvation also regulated expression of many components of the unfolded protein response, the first such report in a species other than yeast (Saccharomyces cerevisiae) and mice (Mus musculus). The response of the zebrafish hepatic transcriptome to starvation was strikingly similar to that of rainbow trout (Oncorhynchus mykiss) and less similar to mouse, while the response of common carp (Cyprinus carpio) differed considerably from the other three species.

Genome-wide expression profiling, in vivo DNA binding analysis, and probabilistic motif prediction reveal novel Abf1 target genes during fermentation, respiration, and sporulation in yeast

The autonomously replicating sequence binding factor 1 (Abf1) was initially identified as an essential DNA replication factor and later shown to be a component of the regulatory network controlling mitotic and meiotic cell cycle progression in budding yeast. The protein is thought to exert its functions via specific interaction with its target site as part of distinct protein complexes, but its roles during mitotic growth and meiotic development are only partially understood. Here, we report a comprehensive approach aiming at the identification of direct Abf1-target genes expressed during fermentation, respiration, and sporulation. Computational prediction of the protein's target sites was integrated with a genome-wide DNA binding assay in growing and sporulating cells. The resulting data were combined with the output of expression profiling studies using wild-type versus temperature-sensitive alleles. This work identified 434 protein-coding loci as being transcriptionally dependent on Abf1. More than 60% of their putative promoter regions contained a computationally predicted Abf1 binding site and/or were bound by Abf1 in vivo, identifying them as direct targets. The present study revealed numerous loci previously unknown to be under Abf1 control, and it yielded evidence for the protein's variable DNA binding pattern during mitotic growth and meiotic development.

Effect of carbon source perturbations on transcriptional regulation of metabolic fluxes in Saccharomyces cerevisiae

BACKGROUND

Control effective flux (CEF) of a reaction is the weighted sum of all fluxes through that reaction, derived from elementary flux modes (EFM) of a metabolic network. Change in CEFs under different environmental conditions has earlier been proven to be correlated with the corresponding changes in the transcriptome. Here we use this to investigate the degree of transcriptional regulation of fluxes in the metabolism of Saccharomyces cerevisiae. We do this by quantifying correlations between changes in CEFs and changes in transcript levels for shifts in carbon source, i.e. between the fermentative carbon source glucose and nonfermentative carbon sources like ethanol, acetate, and lactate. The CEF analysis is based on a simple stoichiometric model that includes reactions of the central carbon metabolism and the amino acid metabolism.

RESULTS

The effect of the carbon shift on the metabolic fluxes was investigated for both batch and chemostat cultures. For growth on glucose in batch (respiro-fermentative) cultures, EFMs with no by-product formation were removed from the analysis of the CEFs, whereas those including any by-products (ethanol, glycerol, acetate, succinate) were omitted in the analysis of growth on glucose in chemostat (respiratory) cultures. This resulted in improved correlations between CEF changes and transcript levels. A regression correlation coefficient of 0.60 was obtained between CEF changes and gene expression changes in the central carbon metabolism for the analysis of 5 different perturbations. Out of 45 data points there were no more than 6 data points deviating from the correlation. Additionally, up- or down-regulation of at least 75% of the genes were in qualitative agreement with the CEF changes for all perturbations studied.

CONCLUSION

The analysis indicates that changes in carbon source are associated with a high degree of hierarchical regulation of metabolic fluxes in the central carbon metabolism as the change in fluxes are correlating directly with the change in transcript levels of genes encoding their corresponding enzymes. For amino acid biosynthesis there was, however, not found to exist a similar correlation, and this may point to either post-transcriptional and/or metabolic regulation, or be due to the absence of a direct perturbation on the amino acid pathways in these experiments.

The unfolded protein response

The unfolded protein response (UPR) is a signal transduction network activated by inhibition of protein folding in the endoplasmic reticulum (ER). The UPR coordinates adaptive responses to this stress situation, including induction of ER resident molecular chaperone and protein foldase expression to increase the protein folding capacity of the ER, induction of phospholipid synthesis, attenuation of general translation, and upregulation of ER-associated degradation to decrease the unfolded protein load of the ER, and an antioxidant response. Upon severe or prolonged ER stress the UPR induces apoptosis to eliminate unhealthy cells from an organism or a population. In this review, I will summarize our current knowledge about signal transduction pathways involved in transducing the unfolded protein signal from the ER to the nucleus or the cytosol.

Ashbya Genome Database 3.0: a cross-species genome and transcriptome browser for yeast biologists

BACKGROUND

The Ashbya Genome Database (AGD) 3.0 is an innovative cross-species genome and transcriptome browser based on release 40 of the Ensembl developer environment.

DESCRIPTION

AGD 3.0 provides information on 4726 protein-encoding loci and 293 non-coding RNA genes present in the genome of the filamentous fungus Ashbya gossypii. A synteny viewer depicts the chromosomal location and orientation of orthologous genes in the budding yeast Saccharomyces cerevisiae. Genome-wide expression profiling data obtained with high-density oligonucleotide microarrays (GeneChips) are available for nearly all currently annotated protein-coding loci in A. gossypii and S. cerevisiae.

CONCLUSION

AGD 3.0 hence provides yeast- and genome biologists with comprehensive report pages including reliable DNA annotation, Gene Ontology terms associated with S. cerevisiae orthologues and RNA expression data as well as numerous links to external sources of information. The database is accessible at http://agd.vital-it.ch/.

Perturbation of the activity of replication origin by meiosis-specific transcription

We have determined the activity of all ARSs on the Saccharomyces cerevisiae chromosome VI as chromosomal replication origins in premeiotic S-phase by neutral/neutral two-dimensional gel electrophoresis. The comparison of origin activity of each origin in mitotic and premeiotic S-phase showed that one of the most efficient origins in mitotic S-phase, ARS605, was completely inhibited in premeiotic S-phase. ARS605 is located within the open reading frame of MSH4 gene that is transcribed specifically during an early stage of meiosis. Systematic analysis of relationships between MSH4 transcription and ARS605 origin activity revealed that transcription of MSH4 inhibited the ARS605 origin activity by removing origin recognition complex from ARS605. Deletion of UME6, a transcription factor responsible for repressing MSH4 during mitotic S-phase, resulted in inactivation of ARS605 in mitosis. Our finding is the first demonstration that the transcriptional regulation on the replication origin activity is related to changes in cell physiology. These results may provide insights into changes in replication origin activity in embryonic cell cycle during early developmental stages.

The GermOnline cross-species systems browser provides comprehensive information on genes and gene products relevant for sexual reproduction

We report a novel release of the GermOnline knowledgebase covering genes relevant for the cell cycle, gametogenesis and fertility. GermOnline was extended into a cross-species systems browser including information on DNA sequence annotation, gene expression and the function of gene products. The database covers eight model organisms and Homo sapiens, for which complete genome annotation data are available. The database is now built around a sophisticated genome browser (Ensembl), our own microarray information management and annotation system (MIMAS) used to extensively describe experimental data obtained with high-density oligonucleotide microarrays (GeneChips) and a comprehensive system for online editing of database entries (MediaWiki). The RNA data include results from classical microarrays as well as tiling arrays that yield information on RNA expression levels, transcript start sites and lengths as well as exon composition. Members of the research community are solicited to help GermOnline curators keep database entries on genes and gene products complete and accurate. The database is accessible at .

A fungal family of transcriptional regulators: the zinc cluster proteins

The trace element zinc is required for proper functioning of a large number of proteins, including various enzymes. However, most zinc-containing proteins are transcription factors capable of binding DNA and are named zinc finger proteins. They form one of the largest families of transcriptional regulators and are categorized into various classes according to zinc-binding motifs. This review focuses on one class of zinc finger proteins called zinc cluster (or binuclear) proteins. Members of this family are exclusively fungal and possess the well-conserved motif CysX(2)CysX(6)CysX(5-12)CysX(2)CysX(6-8)Cys. The cysteine residues bind to two zinc atoms, which coordinate folding of the domain involved in DNA recognition. The first- and best-studied zinc cluster protein is Gal4p, a transcriptional activator of genes involved in the catabolism of galactose in the budding yeast Saccharomyces cerevisiae. Since the discovery of Gal4p, many other zinc cluster proteins have been characterized; they function in a wide range of processes, including primary and secondary metabolism and meiosis. Other roles include regulation of genes involved in the stress response as well as pleiotropic drug resistance, as demonstrated in budding yeast and in human fungal pathogens. With the number of characterized zinc cluster proteins growing rapidly, it is becoming more and more apparent that they are important regulators of fungal physiology.

A scalable method for integration and functional analysis of multiple microarray datasets

MOTIVATION

The diverse microarray datasets that have become available over the past several years represent a rich opportunity and challenge for biological data mining. Many supervised and unsupervised methods have been developed for the analysis of individual microarray datasets. However, integrated analysis of multiple datasets can provide a broader insight into genetic regulation of specific biological pathways under a variety of conditions.

RESULTS

To aid in the analysis of such large compendia of microarray experiments, we present Microarray Experiment Functional Integration Technology (MEFIT), a scalable Bayesian framework for predicting functional relationships from integrated microarray datasets. Furthermore, MEFIT predicts these functional relationships within the context of specific biological processes. All results are provided in the context of one or more specific biological functions, which can be provided by a biologist or drawn automatically from catalogs such as the Gene Ontology (GO). Using MEFIT, we integrated 40 Saccharomyces cerevisiae microarray datasets spanning 712 unique conditions. In tests based on 110 biological functions drawn from the GO biological process ontology, MEFIT provided a 5% or greater performance increase for 54 functions, with a 5% or more decrease in performance in only two functions.

When transcriptome meets metabolome: fast cellular responses of yeast to sudden relief of glucose limitation

Within the first 5 min after a sudden relief from glucose limitation, Saccharomyces cerevisiae exhibited fast changes of intracellular metabolite levels and a major transcriptional reprogramming. Integration of transcriptome and metabolome data revealed tight relationships between the changes at these two levels. Transcriptome as well as metabolite changes reflected a major investment in two processes: adaptation from fully respiratory to respiro-fermentative metabolism and preparation for growth acceleration. At the metabolite level, a severe drop of the AXP pools directly after glucose addition was not accompanied by any of the other three NXP. To counterbalance this loss, purine biosynthesis and salvage pathways were transcriptionally upregulated in a concerted manner, reflecting a sudden increase of the purine demand. The short-term dynamics of the transcriptome revealed a remarkably fast decrease in the average half-life of downregulated genes. This acceleration of mRNA decay can be interpreted both as an additional nucleotide salvage pathway and an additional level of glucose-induced regulation of gene expression.

Natural isolates of Saccharomyces cerevisiae display complex genetic variation in sporulation efficiency

Sporulation is a well-studied process executed with varying efficiency by diverse yeast strains. We developed a high-throughput method to quantify yeast sporulation efficiency and used this technique to analyze a line cross between a high-efficiency oak tree isolate and a low-efficiency wine strain. We find that natural variation in sporulation efficiency mirrors natural variation in higher eukaryotes: it shows divergence between isolated populations, arises from loci of major effect, and exhibits epistasis. We show that the lower sporulation efficiency of the wine strain results from a failure to initiate sporulation, rather than from slower kinetics of meiosis and spore formation. The two strains differentially regulate many genes involved in aerobic respiration, an essential pathway for sporulation, such that the oak tree strain appears better poised to generate energy from this pathway. We also report that a polymorphism in RME1 that affects sporulation efficiency in laboratory strains also cosegregates with significant phenotypic differences in our cross of natural isolates. These results lay the groundwork for the study of variation in sporulation efficiency among natural isolates of yeast.

Extensive low-affinity transcriptional interactions in the yeast genome

Major experimental and computational efforts are targeted at the characterization of transcriptional networks on a genomic scale. The ultimate goal of many of these studies is to construct networks associating transcription factors with genes via well-defined binding sites. Weaker regulatory interactions other than those occurring at high-affinity binding sites are largely ignored and are not well understood. Here I show that low-affinity interactions are abundant in vivo and quantifiable from current high-throughput ChIP experiments. I develop algorithms that predict DNA-binding energies from sequences and ChIP data across a wide dynamic range of affinities and use them to reveal widespread functionality of low-affinity transcription factor binding. Evolutionary analysis suggests that binding energies of many transcription factors are conserved even in promoters lacking classical binding sites. Gene expression analysis shows that such promoters can generate significant expression. I estimate that while only a small percentage of the genome is strongly regulated by a typical transcription factor, up to an order of magnitude more may be involved in weaker interactions. Low-affinity transcription factor-DNA interaction may therefore be important both evolutionarily and functionally.

The mammalian unfolded protein response

In the endoplasmic reticulum (ER), secretory and transmembrane proteins fold into their native conformation and undergo posttranslational modifications important for their activity and structure. When protein folding in the ER is inhibited, signal transduction pathways, which increase the biosynthetic capacity and decrease the biosynthetic burden of the ER to maintain the homeostasis of this organelle, are activated. These pathways are called the unfolded protein response (UPR). In this review, we briefly summarize principles of protein folding and molecular chaperone function important for a mechanistic understanding of UPR-signaling events. We then discuss mechanisms of signal transduction employed by the UPR in mammals and our current understanding of the remodeling of cellular processes by the UPR. Finally, we summarize data that demonstrate that UPR signaling feeds into decision making in other processes previously thought to be unrelated to ER function, e.g., eukaryotic starvation responses and differentiation programs.

Novel response to microtubule perturbation in meiosis

During the mitotic cell cycle, microtubule depolymerization leads to a cell cycle arrest in metaphase, due to activation of the spindle checkpoint. Here, we show that under microtubule-destabilizing conditions, such as low temperature or the presence of the spindle-depolymerizing drug benomyl, meiotic budding yeast cells arrest in G(1) or G(2), instead of metaphase. Cells arrest in G(1) if microtubule perturbation occurs as they enter the meiotic cell cycle and in G(2) if cells are already undergoing premeiotic S phase. Concomitantly, cells down-regulate genes required for cell cycle progression, meiotic differentiation, and spore formation in a highly coordinated manner. Decreased expression of these genes is likely to be responsible for halting both cell cycle progression and meiotic development. Our results point towards the existence of a novel surveillance mechanism of microtubule integrity that may be particularly important during specialized cell cycles when coordination of cell cycle progression with a developmental program is necessary.

Evidence for abundant transcription of non-coding regions in the Saccharomyces cerevisiae genome

Background

Recent studies in a growing number of organisms have yielded accumulating evidence that a significant portion of the non-coding region in the genome is transcribed. We address this issue in the yeast Saccharomyces cerevisiae.

Results

Taking into account the absence of a significantly large yeast EST database, we use microarray expression data collected for genomic regions erroneously believed to be coding to study the expression pattern of non-coding regions in the Saccharomyces cerevisiae genome. We find that at least 164 out of 589 (28%) such regions are expressed under specific biological conditions. In particular, looking at the probes that are located opposing other known genes at the same genomic locus, we find that 88 out of 341 (26%) of these genes support antisense transcription. The expression patterns of these antisense genes are positively correlated. We validate these results using RT-PCR on a sample of 6 non-coding transcripts.

Conclusion

1. The yeast genome is transcribed on a scale larger than previously assumed. 2. Correlated transcription of antisense genes is abundant in the yeast genome. 3. Antisense genes in yeast are non-coding.

Identification of functional transcription factor binding sites using closely related Saccharomyces species

Comparative genomics provides a rapid means of identifying functional DNA elements by their sequence conservation between species. Transcription factor binding sites (TFBSs) may constitute a significant fraction of these conserved sequences, but the annotation of specific TFBSs is complicated by the fact that these short, degenerate sequences may frequently be conserved by chance rather than functional constraint. To identify intergenic sequences that function as TFBSs, we calculated the probability of binding site conservation between Saccharomyces cerevisiae and its two closest relatives under a neutral model of evolution. We found that this probability is <5% for 134 of 163 transcription factor binding motifs, implying that we can reliably annotate binding sites for the majority of these transcription factors by conservation alone. Although our annotation relies on a number of assumptions, mutations in five of five conserved Ume6 binding sites and three of four conserved Ndt80 binding sites show Ume6- and Ndt80-dependent effects on gene expression. We also found that three of five unconserved Ndt80 binding sites show Ndt80-dependent effects on gene expression. Together these data imply that although sequence conservation can be reliably used to predict functional TFBSs, unconserved sequences might also make a significant contribution to a species' biology.

Cyclophilin A is localized to the nucleus and controls meiosis in Saccharomyces cerevisiae

Cyclophilin A is conserved from yeast to humans and mediates the ability of cyclosporine to perturb signal transduction cascades via inhibition of calcineurin. Cyclophilin A also catalyzes cis-trans peptidyl-prolyl isomerization during protein folding or conformational changes; however, cyclophilin A is not essential in yeast or human cells, and the true biological functions of this highly conserved enzyme have remained enigmatic. In Saccharomyces cerevisiae, cyclophilin A becomes essential in cells compromised for the nuclear prolyl-isomerase Ess1, and cyclophilin A physically interacts with two nuclear histone deacetylase complexes, Sin3-Rpd3 and Set3C, which both control meiosis. Here we show that cyclophilin A is localized to the nucleus in yeast cells and governs the meiotic gene program to promote efficient sporulation. The prolyl-isomerase activity of cyclophilin A is required for this meiotic function. We document that cyclophilin A physically associates with the Set3C histone deacetylase and analyze in detail the structure of this protein-protein complex. Genetic studies support a model in which cyclophilin A controls meiosis via Set3C and an additional target. Our findings reveal a novel nuclear role for cyclophilin A in governing the transcriptional program required for the vegetative to meiotic developmental switch in budding yeast.

Identification, mutational analysis, and coactivator requirements of two distinct transcriptional activation domains of the Saccharomyces cerevisiae Hap4 protein

The Hap4 protein of the budding yeast Saccharomyces cerevisiae activates the transcription of genes that are required for growth on nonfermentable carbon sources. Previous reports suggested the presence of a transcriptional activation domain within the carboxyl-terminal half of Hap4 that can function in the absence of Gcn5, a transcriptional coactivator protein and histone acetyltransferase. The boundaries of this activation domain were further defined to a region encompassing amino acids 359 to 476. Within this region, several clusters of hydrophobic amino acids are critical for transcriptional activity. This activity does not require GCN5 or two other components of the SAGA coactivator complex, SPT3 and SPT8, but it does require SPT7 and SPT20. Contrary to previous reports, a Hap4 fragment comprising amino acids 1 to 330 can support the growth of yeast on lactate medium, and when tethered to lexA, can activate a reporter gene with upstream lexA binding sites, demonstrating the presence of a second transcriptional activation domain. In contrast to the C-terminal activation domain, the transcriptional activity of this N-terminal region depends on GCN5. We conclude that the yeast Hap4 protein has at least two transcriptional activation domains with strikingly different levels of dependence on specific transcriptional coactivator proteins.

Rank Difference Analysis of Microarrays (RDAM), a novel approach to statistical analysis of microarray expression profiling data

Background

A key step in the analysis of microarray expression profiling data is the identification of genes that display statistically significant changes in expression signals between two biological conditions.

Results

We describe a new method, Rank Difference Analysis of Microarrays (RDAM), which estimates the total number of truly varying genes and assigns a p-value to each signal variation. Information on a group of differentially expressed genes includes the sensitivity and the false discovery rate. We demonstrate the feasibility and efficiency of our approach by applying it to a large synthetic expression data set and to a biological data set obtained by comparing vegetatively-growing wild type and tor2-mutant yeast strains. In both cases we observed a significant improvement of the power of analysis when our method is compared to another popular nonparametric method.

Conclusions

This study provided a valuable new statistical method to analyze microarray data. We conclude that the good quality of the results obtained by RDAM is mainly due to the quasi-perfect equalization of variation distribution, which is related to the standardization procedure used and to the measurement of variation by rank difference.

The in vivo activity of Ime1, the key transcriptional activator of meiosis-specific genes in Saccharomyces cerevisiae, is inhibited by the cyclic AMP/protein kinase A signal pathway through the glycogen synthase kinase 3-beta homolog Rim11

Phosphorylation is the main mode by which signals are transmitted to key regulators of developmental pathways. The glycogen synthase kinase 3 family plays pivotal roles in the development and well-being of all eukaryotic organisms. Similarly, the budding yeast homolog Rim11 is essential for the exit of diploid cells from the cell cycle and for entry into the meiotic developmental pathway. In this report we show that in vivo, in cells grown in a medium promoting vegetative growth with acetate as the sole carbon source (SA medium), Rim11 phosphorylates Ime1, the master transcriptional activator required for entry into the meiotic cycle and for the transcription of early meiosis-specific genes. We demonstrate that in the presence of glucose, the kinase activity of Rim11 is inhibited. This inhibition could be due to phosphorylation on Ser-5, Ser-8, and/or Ser-12 because in the rim11S5AS8AS12A mutant, Ime1 is incorrectly phosphorylated in the presence of glucose and cells undergo sporulation. We further show that this nutrient signal is transmitted to Rim11 and consequently to Ime1 by the cyclic AMP/protein kinase A signal transduction pathway. Ime1 is phosphorylated in SA medium on at least two residues, Tyr-359 and Ser-302 and/or Ser-306. Ser-302 and Ser-306 are part of a consensus site for the mammalian homolog of Rim11, glycogen synthase kinase 3-beta. Phosphorylation on Tyr-359 but not Ser-302 or Ser-306 is essential for the transcription of early meiosis-specific genes and sporulation. We show that Tyr-359 is phosphorylated by Rim11.

The Rpd3-Sin3 histone deacetylase regulates replication timing and enables intra-S origin control in Saccharomyces cerevisiae

The replication of eukaryotic genomes follows a temporally staged program, in which late origin firing often occurs within domains of altered chromatin structure(s) and silenced genes. Histone deacetylation functions in gene silencing in some late-replicating regions, prompting an investigation of the role of histone deacetylation in replication timing control in Saccharomyces cerevisiae. Deletion of the histone deacetylase Rpd3 or its interacting partner Sin3 caused early activation of late origins at internal chromosomal loci but did not alter the initiation timing of early origins or a late-firing, telomere-proximal origin. By delaying initiation relative to the earliest origins, Rpd3 enables regulation of late origins by the intra-S replication checkpoint. RPD3 deletion suppresses the slow S phase of clb5Delta cells by enabling late origins to fire earlier, suggesting that Rpd3 modulates the initiation timing of many origins throughout the genome. Examination of factors such as Ume6 that function together with Rpd3 in transcriptional repression indicates that Rpd3 regulates origin initiation timing independently of its role in transcriptional repression. This supports growing evidence that for much of the S. cerevisiae genome transcription and replication timing are not linked.

The unfolded protein response represses differentiation through the RPD3-SIN3 histone deacetylase

In Saccharomyces cerevisiae, splicing of HAC1 mRNA is initiated in response to the accumulation of unfolded proteins in the endoplasmic reticulum by the transmembrane kinase-endoribonuclease Ire1p. Spliced Hac1p (Hac1ip) is a negative regulator of differentiation responses to nitrogen starvation, pseudohyphal growth, and meiosis. Here we show that the RPD3-SIN3 histone deacetylase complex (HDAC), its catalytic activity, recruitment of the HDAC to the promoters of early meiotic genes (EMGs) by Ume6p, and the Ume6p DNA-binding site URS1 in the promoters of EMGs are required for nitrogen-mediated negative regulation of EMGs and meiosis by Hac1ip. Co-immunoprecipitation experiments demonstrated that Hac1ip can interact with the HDAC in vivo. Systematic analysis of double deletion strains revealed that HAC1 is a peripheral component of the HDAC. In summary, nitrogen-induced synthesis of Hac1ip and association of Hac1ip with the HDAC are physiological events in the regulation of EMGs by nutrients. These data also define for the first time a gene class that is under negative control by the UPR, and provide the framework for a novel mechanism through which bZIP proteins repress transcription.

Lessons from the genome sequence of Neurospora crassa: tracing the path from genomic blueprint to multicellular organism

We present an analysis of over 1,100 of the approximately 10,000 predicted proteins encoded by the genome sequence of the filamentous fungus Neurospora crassa. Seven major areas of Neurospora genomics and biology are covered. First, the basic features of the genome, including the automated assembly, gene calls, and global gene analyses are summarized. The second section covers components of the centromere and kinetochore complexes, chromatin assembly and modification, and transcription and translation initiation factors. The third area discusses genome defense mechanisms, including repeat induced point mutation, quelling and meiotic silencing, and DNA repair and recombination. In the fourth section, topics relevant to metabolism and transport include extracellular digestion; membrane transporters; aspects of carbon, sulfur, nitrogen, and lipid metabolism; the mitochondrion and energy metabolism; the proteasome; and protein glycosylation, secretion, and endocytosis. Environmental sensing is the focus of the fifth section with a treatment of two-component systems; GTP-binding proteins; mitogen-activated protein, p21-activated, and germinal center kinases; calcium signaling; protein phosphatases; photobiology; circadian rhythms; and heat shock and stress responses. The sixth area of analysis is growth and development; it encompasses cell wall synthesis, proteins important for hyphal polarity, cytoskeletal components, the cyclin/cyclin-dependent kinase machinery, macroconidiation, meiosis, and the sexual cycle. The seventh section covers topics relevant to animal and plant pathogenesis and human disease. The results demonstrate that a large proportion of Neurospora genes do not have homologues in the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe. The group of unshared genes includes potential new targets for antifungals as well as loci implicated in human and plant physiology and disease.

ISWI complexes in Saccharomyces cerevisiae

The imitation switch (ISWI) class of chromatin remodeling ATPase is ubiquitous in eukaryotes. It is becoming clear that these enzymes exist as part of larger complexes and the nature of the associated proteins dictate the function associated with a complex both in biochemical assays and in the cell. Much progress has been made in understanding these relationships in the budding yeast Saccharomyces cerevisiae, containing two ATPases, Isw1p and Isw2p. This has been aided by the ease of genetic manipulation, by a number of systematic screens designed to specifically detect ISWI function and by the plethora of data generated from a number of global screens for function. At present, many functions for yeast Isw1p and Isw2p are related to effects on RNA levels and are associated with the controlled repression of gene expression that crudely fall into three types: displacement of the basal transcription machinery to repress or silence transcription of genes (Isw2 complex and Isw1/Ioc3 complex); control of the activation of expression leading to coordination of transcription elongation; and efficient termination of transcription (Isw1/Ioc4/Ioc2 complex). The latter two functions are regulated by specific phosphorylation of residues within the carboxy terminal domain (CTD) of the largest subunit of RNA polymerase II (RNAPII). Other functions may relate to the ability of ISWI complex to displace transcription factors or enzymes from the template. Other ISWI-containing complexes that have yet to be characterized indicate that much remains to be learnt about yeast ISWI itself and importantly, how the various forms cooperate with different classes of chromatin remodeling ATPase, complexes containing histone acetylases, deacetylases, methylases and both DNA and RNA polymerases.

GermOnline, a cross-species community knowledgebase on germ cell differentiation

GermOnline provides information and microarray expression data for genes involved in mitosis and meiosis, gamete formation and germ line development across species. The database has been developed, and is being curated and updated, by life scientists in cooperation with bioinformaticists. Information is contributed through an online form using free text, images and the controlled vocabulary developed by the GeneOntology Consortium. Authors provide up to three references in support of their contribution. The database is governed by an international board of scientists to ensure a standardized data format and the highest quality of GermOnline's information content. Release 2.0 provides exclusive access to microarray expression data from Saccharomyces cerevisiae and Rattus norvegicus, as well as curated information on approximately 700 genes from various organisms. The locus report pages include links to external databases that contain relevant annotation, microarray expression and proteome data. Conversely, the Saccharomyces Genome Database (SGD), S.cerevisiae GeneDB and Swiss-Prot link to the budding yeast section of GermOnline from their respective locus pages. GermOnline, a fully operational prototype subject-oriented knowledgebase designed for community annotation and array data visualization, is accessible at http://www.germonline.org. The target audience includes researchers who work on mitotic cell division, meiosis, gametogenesis, germ line development, human reproductive health and comparative genomics.

Phylogenetically and spatially conserved word pairs associated with gene-expression changes in yeasts

BACKGROUND

Transcriptional regulation in eukaryotes often involves multiple transcription factors binding to the same transcription control region, and to understand the regulatory content of eukaryotic genomes it is necessary to consider the co-occurrence and spatial relationships of individual binding sites. The determination of conserved sequences (often known as phylogenetic footprinting) has identified individual transcription factor binding sites. We extend this concept of functional conservation to higher-order features of transcription control regions.

RESULTS

We used the genome sequences of four yeast species of the genus Saccharomyces to identify sequences potentially involved in multifactorial control of gene expression. We found 989 potential regulatory 'templates': pairs of hexameric sequences that are jointly conserved in transcription regulatory regions and also exhibit non-random relative spacing. Many of the individual sequences in these templates correspond to known transcription factor binding sites, and the sets of genes containing a particular template in their transcription control regions tend to be differentially expressed in conditions where the corresponding transcription factors are known to be active. The incorporation of word pairs to define sequence features yields more specific predictions of average expression profiles and more informative regression models for genome-wide expression data than considering sequence conservation alone.

CONCLUSIONS

The incorporation of both joint conservation and spacing constraints of sequence pairs predicts groups of target genes that are specific for common patterns of gene expression. Our work suggests that positional information, especially the relative spacing between transcription factor binding sites, may represent a common organizing principle of transcription control regions.

Signal pathway integration in the switch from the mitotic cell cycle to meiosis in yeast

Diploid yeast, like most eukaryotes, can undergo meiotic differentiation to form haploid gametes. Meiotic differentiation and cell growth (proliferation) are mutually exclusive programs, and in yeast the switch between growth and meiosis is controlled by nutritional signals. The signaling pathways that mediate nutritional controls on meiotic initiation fall into three broad classes: those that respond to nutrient starvation, those that respond to non-fermentable carbon sources, and those that respond to glucose. At the onset of meiosis, nutritional signaling pathways converge on transcriptional regulation of two genes: IME1, which encodes a transcription factor; and IME2, which encodes a protein kinase. Transcription of IME1 and IME2 trigger initiation of meiosis, and the expression of these two genes is linked with one other, with expression of later meiotic genes and with early meiotic events such as DNA replication. In addition, the signaling pathways that control IME1 and IME2 expression are themselves integrated through a variety of mechanisms. Thus the signal network that controls the switch from growth to meiotic differentiation provides a signaling code that translates different combinations of extracellular signals into appropriate cellular responses.