Sum1 and Ndt80 proteins compete for binding to middle sporulation element sequences that control meiotic gene expression

Avoiding unscheduled transcription in shared promoters: Saccharomyces cerevisiae Sum1p represses the divergent gene pair SPS18-SPS19 through a midsporulation element (MSE)

The sporulation-specific gene SPS18 shares a common promoter region with the oleic acid-inducible gene SPS19. Both genes are transcribed in sporulating diploid cells, albeit unevenly in favour of SPS18, whereas in haploid cells grown on fatty acids only SPS19 is highly activated. Here, SPS19 oleate-response element (ORE) conferred activation on a basal CYC1-lacZ reporter gene equally in both orientations, but promoter analysis using SPS18-lacZ reporter constructs with deletions identified a repressing fragment containing a midsporulation element (MSE) that could be involved in imposing directionality towards SPS19 in oleic acid-induced cells. In sporulating diploids, MSEs recruit the Ndt80p transcription factor for activation, whereas under vegetative conditions, certain MSEs are targeted by the Sum1p repressor in association with Hst1p and Rfm1p. Quantitative real-time PCR demonstrated that in haploid sum1Δ, hst1Δ, or rfm1Δ cells, oleic acid-dependent expression of SPS18 was higher compared with the situation in wild-type cells, but in the sum1Δ mutant, this effect was diminished in the absence of Oaf1p or Pip2p. We conclude that SPS18 MSE is a functional element repressing the expression of both SPS18 and SPS19, and is a component of a stricture mechanism shielding SPS18 from the dramatic increase in ORE-dependent transcription of SPS19 in oleic acid-grown cells.

Role of Ndt80p in sterol metabolism regulation and azole resistance in Candida albicans

The Ndt80p transcription factor modulates azole tolerance in Candida albicans by controlling the expression of the gene for the drug efflux pump Cdr1p. To date, the contribution of this transcriptional modulator to drug tolerance is not yet well understood. Here, we investigate the role of Ndt80p in mediating fluconazole tolerance by determining its genome-wide occupancy using chromatin immunoprecipitation coupled to high-density tiling arrays. Ndt80p was found to bind a large number of gene promoters with diverse biological functions. Gene ontology analysis of these Ndt80p targets revealed a significant enrichment in gene products related to the cell wall, carbohydrate metabolism, stress responses, hyphal development, multidrug transport, and the cell cycle. Ndt80p was found on the promoters of ergosterol biosynthesis genes, including on the azole target Erg11p. Additionally, expression profiling was used to identify fluconazole-responsive genes that require Ndt80p for their proper expression. We found that Ndt80p is crucial for the expression of numerous fluconazole-responsive genes, especially genes involved in ergosterol metabolism. Therefore, by combining genome-wide location and transcriptional profiling, we have characterized the Ndt80p fluconazole-dependent regulon and demonstrated the key role of this global transcriptional regulator in modulating sterol metabolism and drug resistance in C. albicans.

The Ime2 protein kinase enhances the disassociation of the Sum1 repressor from middle meiotic promoters

Meiotic development in Saccharomyces cerevisiae (sporulation) is controlled by the sequential transcription of temporally distinct sets of meiosis-specific genes. The induction of middle genes controls exit from meiotic prophase, the completion of the nuclear divisions, and spore formation. Middle promoters are controlled through DNA elements termed middle sporulation elements (MSEs) that are bound by the Sum1 repressor during vegetative growth and by the Ndt80 activator during meiosis. It has been proposed that the induction of middle promoters is controlled by competition between Ndt80 and Sum1 for MSE occupancy. Here, we show that the Sum1 repressor can be removed from middle promoters in meiotic cells independent of Ndt80 expression. This process requires the phosphorylation of Sum1 by the meiosis-specific cyclin-dependent kinase-like kinase Ime2. The deletion of HST1, which encodes a Sir2 paralog that interacts with Sum1, bypasses the requirement for this phosphorylation. These findings suggest that in the presence of Ndt80, Sum1 may be displaced from MSEs through a competition-based mechanism but that in the absence of Ndt80, Sum1 is removed from chromatin in a separate pathway requiring the phosphorylation of Sum1 by Ime2 and the inhibition of Hst1.

Environment-specific combinatorial cis-regulation in synthetic promoters

When a cell's environment changes, a large transcriptional response often takes place. The exquisite sensitivity and specificity of these responses are controlled in large part by the combinations of cis-regulatory elements that reside in gene promoters and adjacent control regions. Here, we present a study aimed at accurately modeling the relationship between combinations of cis-regulatory elements and the expression levels they drive in different environments. We constructed four libraries of synthetic promoters in yeast, consisting of combinations of transcription factor binding sites and assayed their expression in four different environments. Thermodynamic models relating promoter sequences to their corresponding four expression levels explained at least 56% of the variation in expression in each library through the different conditions. Analyses of these models suggested that a large fraction of regulated gene expression is explained by changes in the effective concentration of sequence-specific transcription factors, and we show that in most cases, the corresponding transcription factors are expressed in a pattern that is predicted by the thermodynamic models. Our analysis uncovered two binding sites that switch from activators to repressors in different environmental conditions. In both the cases, the switch was not the result of a single transcription factor changing regulatory modes, but most likely due to competition between multiple factors binding to the same site. Our analysis suggests that this mode of regulation allows for large and steep changes in expression in response to changing transcription factor concentrations. Our results demonstrate that many complex changes in gene expression are accurately explained by simple changes in the effective concentrations of transcription factors.

A library of yeast transcription factor motifs reveals a widespread function for Rsc3 in targeting nucleosome exclusion at promoters

The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially approximately 100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.

Control of replication initiation by the Sum1/Rfm1/Hst1 histone deacetylase

BACKGROUND

Replication initiation at origins of replication in the yeast genome takes place on chromatin as a template, raising the question how histone modifications, for instance histone acetylation, influence origin firing. Initiation requires binding of the replication initiator, the Origin Recognition Complex (ORC), to a consensus sequence within origins. In addition, other proteins bind to recognition sites in the vicinity of ORC and support initiation. In previous work, we identified Sum1 as an origin-binding protein that contributes to efficient replication initiation. Sum1 is part of the Sum1/Rfm1/Hst1 complex that represses meiotic genes during vegetative growth via histone deacetylation by the histone deacetylase (HDAC) Hst1.

RESULTS

In this study, we investigated how Sum1 affected replication initiation. We found that it functioned in initiation as a component of the Sum1/Rfm1/Hst1 complex, implying a role for histone deacetylation in origin activity. We identified several origins in the yeast genome whose activity depended on both Sum1 and Hst1. Importantly, sum1Delta or hst1Delta caused a significant increase in histone H4 lysine 5 (H4 K5) acetylation levels, but not other H4 acetylation sites, at those origins. Furthermore, mutation of lysines to glutamines in the H4 tail, which imitates the constantly acetylated state, resulted in a reduction of origin activity comparable to that in the absence of Hst1, showing that deacetylation of H4 was important for full initiation capacity of these origins.

CONCLUSION

Taken together, our results demonstrate a role for histone deacetylation in origin activity and reveal a novel aspect of origin regulation by chromatin. These results suggest recruitment of the Sum1/Rfm1/Hst1 complex to a number of yeast origins, where Hst1 deacetylated H4 K5.

Cdc7-Dbf4 regulates NDT80 transcription as well as reductional segregation during budding yeast meiosis

In budding yeast, as in other eukaryotes, the Cdc7 protein kinase is important for initiation of DNA synthesis in vegetative cells. In addition, Cdc7 has crucial meiotic functions: it facilitates premeiotic DNA replication, and it is essential for the initiation of recombination. This work uses a chemical genetic approach to demonstrate that Cdc7 kinase has additional roles in meiosis. First, Cdc7 allows expression of NDT80, a meiosis-specific transcriptional activator required for the induction of genes involved in exit from pachytene, meiotic progression, and spore formation. Second, Cdc7 is necessary for recruitment of monopolin to sister kinetochores, and it is necessary for the reductional segregation occurring at meiosis I. The use of the same kinase to regulate several distinct meiosis-specific processes may be important for the coordination of these processes during meiosis.

A catalog of neutral and deleterious polymorphism in yeast

The abundance and identity of functional variation segregating in natural populations is paramount to dissecting the molecular basis of quantitative traits as well as human genetic diseases. Genome sequencing of multiple organisms of the same species provides an efficient means of cataloging rearrangements, insertion, or deletion polymorphisms (InDels) and single-nucleotide polymorphisms (SNPs). While inbreeding depression and heterosis imply that a substantial amount of polymorphism is deleterious, distinguishing deleterious from neutral polymorphism remains a significant challenge. To identify deleterious and neutral DNA sequence variation within Saccharomyces cerevisiae, we sequenced the genome of a vineyard and oak tree strain and compared them to a reference genome. Among these three strains, 6% of the genome is variable, mostly attributable to variation in genome content that results from large InDels. Out of the 88,000 polymorphisms identified, 93% are SNPs and a small but significant fraction can be attributed to recent interspecific introgression and ectopic gene conversion. In comparison to the reference genome, there is substantial evidence for functional variation in gene content and structure that results from large InDels, frame-shifts, and polymorphic start and stop codons. Comparison of polymorphism to divergence reveals scant evidence for positive selection but an abundance of evidence for deleterious SNPs. We estimate that 12% of coding and 7% of noncoding SNPs are deleterious. Based on divergence among 11 yeast species, we identified 1,666 nonsynonymous SNPs that disrupt conserved amino acids and 1,863 noncoding SNPs that disrupt conserved noncoding motifs. The deleterious coding SNPs include those known to affect quantitative traits, and a subset of the deleterious noncoding SNPs occurs in the promoters of genes that show allele-specific expression, implying that some cis-regulatory SNPs are deleterious. Our results show that the genome sequences of both closely and distantly related species provide a means of identifying deleterious polymorphisms that disrupt functionally conserved coding and noncoding sequences.

DNA sequence variation makes an important contribution to most traits that vary in natural populations. However, mapping mutations that underlie a trait of interest is a significant challenge. Genome sequencing of multiple organisms provides a complete list of DNA sequence differences responsible for any trait that differs among the organisms. Yet, distinguishing those DNA sequence variants that contribute to a trait from all other variants is not easy. Here, we sequence the genomes of two strains of yeast and, through comparisons with a reference genome, we catalog multiple types of DNA sequence variation among the three strains. Using a variety of comparative genomics methods, we show that a substantial fraction of DNA sequence variations has deleterious effects on fitness. Finally, we show that a subset of deleterious mutations is associated with changes in gene expression levels. Our results imply that comparative genomics methods will be a valuable approach to identifying DNA sequence changes underlying numerous traits of interest.

Interspecies variation reveals a conserved repressor of alpha-specific genes in Saccharomyces yeasts

The mating-type determination circuit in Saccharomyces yeast serves as a classic paradigm for the genetic control of cell type in all eukaryotes. Using comparative genetics, we discovered a central and conserved, yet previously undetected, component of this genetic circuit: active repression of alpha-specific genes in a cells. Upon inactivation of the SUM1 gene in Saccharomyces bayanus, a close relative of Saccharomyces cerevisiae, a cells acquired mating characteristics of alpha cells and displayed autocrine activation of their mating response pathway. Sum1 protein bound to the promoters of alpha-specific genes, repressing their transcription. In contrast to the standard model, alpha1 was important but not required for alpha-specific gene activation and mating of alpha cells in the absence of Sum1. Neither Sum1 protein expression, nor its association with target promoters was mating-type-regulated. Thus, the alpha1/Mcm1 coactivators did not overcome repression by occluding Sum1 binding to DNA. Surprisingly, the mating-type regulatory function of Sum1 was conserved in S. cerevisiae. We suggest that a comprehensive understanding of some genetic pathways may be best attained through the expanded phenotypic space provided by study of those pathways in multiple related organisms.

The Saccharomyces cerevisiae CLB5 promoter contains two middle sporulation elements (MSEs) that are differentially regulated during sporulation

The B-type cyclins Clb5 and Clb6 are essential activators of DNA replication during sporulation in Saccharomyces cerevisiae. The expression of CLB5 is maximally induced during the middle phase of sporulation by the transcription factor Ndt80. We have performed an analysis of the CLB5 promoter and have identified two middle sporulation elements (MSEs) that act as binding sites for Ndt80. Although both MSE sequences bind Ndt80 in vitro, they display differential effectiveness in their ability to function as cis-acting regulatory sequences in vivo. Mutation of both MSE sequences in the CLB5 promoter profoundly reduces the induction of CLB5 transcription during the middle phase of sporulation but results in no obvious defect in progression through meiosis and sporulation, implying that the Ndt80-dependent induction of CLB5 is not required for effective DNA replication or chromosome division.

Functional characterization of variations on regulatory motifs

Transcription factors (TFs) regulate gene expression through specific interactions with short promoter elements. The same regulatory protein may recognize a variety of related sequences. Moreover, once they are detected it is hard to predict whether highly similar sequence motifs will be recognized by the same TF and regulate similar gene expression patterns, or serve as binding sites for distinct regulatory factors. We developed computational measures to assess the functional implications of variations on regulatory motifs and to compare the functions of related sites. We have developed computational means for estimating the functional outcome of substituting a single position within a binding site and applied them to a collection of putative regulatory motifs. We predict the effects of nucleotide variations within motifs on gene expression patterns. In cases where such predictions could be compared to suitable published experimental evidence, we found very good agreement. We further accumulated statistics from multiple substitutions across various binding sites in an attempt to deduce general properties that characterize nucleotide substitutions that are more likely to alter expression. We found that substitutions involving Adenine are more likely to retain the expression pattern and that substitutions involving Guanine are more likely to alter expression compared to the rest of the substitutions. Our results should facilitate the prediction of the expression outcomes of binding site variations. One typical important implication is expected to be the ability to predict the phenotypic effect of variation in regulatory motifs in promoters.

Evolution of new function through a single amino acid change in the yeast repressor Sum1p

The SUM1-1 mutation is an example of a single amino acid change that results in new function. Wild-type Sum1p in Saccharomyces cerevisiae is a DNA-binding repressor that acts locally, whereas mutant Sum1-1p forms an extended repressive chromatin structure. By characterizing a panel of mutations in which various amino acids replaced the critical residue, threonine 988, we found that threonine was required for wild-type function and that in the absence of threonine the association of Sum1p with DNA was reduced. Isoleucine, the amino acid in mutant Sum1-1p, was required for the novel spreading property. Thus, the SUM1-1 mutation results in both a loss and a gain of function. The presence of isoleucine caused Sum1-1p to self-associate, a property that may promote spreading. In addition, isoleucine enabled Sum1-1p to associate with the origin recognition complex (ORC) and accumulate near ORC binding sites. Thus, both threonine and isoleucine at position 988 enable Sum1p to form intermolecular interactions. We propose that interaction domains may be hotspots for gain-of-function mutations because alterations in such domains have the potential to redirect a protein to new sets of binding partners. In addition, self-association of chromatin proteins may promote the formation of extended chromatin structures.

A cis-regulatory logic simulator

Background

A major goal of computational studies of gene regulation is to accurately predict the expression of genes based on the cis-regulatory content of their promoters. The development of computational methods to decode the interactions among cis-regulatory elements has been slow, in part, because it is difficult to know, without extensive experimental validation, whether a particular method identifies the correct cis-regulatory interactions that underlie a given set of expression data. There is an urgent need for test expression data in which the interactions among cis-regulatory sites that produce the data are known. The ability to rapidly generate such data sets would facilitate the development and comparison of computational methods that predict gene expression patterns from promoter sequence.

Results

We developed a gene expression simulator which generates expression data using user-defined interactions between cis-regulatory sites. The simulator can incorporate additive, cooperative, competitive, and synergistic interactions between regulatory elements. Constraints on the spacing, distance, and orientation of regulatory elements and their interactions may also be defined and Gaussian noise can be added to the expression values. The simulator allows for a data transformation that simulates the sigmoid shape of expression levels from real promoters. We found good agreement between sets of simulated promoters and predicted regulatory modules from real expression data. We present several data sets that may be useful for testing new methodologies for predicting gene expression from promoter sequence.

Conclusion

We developed a flexible gene expression simulator that rapidly generates large numbers of simulated promoters and their corresponding transcriptional output based on specified interactions between cis-regulatory sites. When appropriate rule sets are used, the data generated by our simulator faithfully reproduces experimentally derived data sets. We anticipate that using simulated gene expression data sets will facilitate the direct comparison of computational strategies to predict gene expression from promoter sequence. The source code is available online and as additional material. The test sets are available as additional material.

Substitution as a mechanism for genetic robustness: the duplicated deacetylases Hst1p and Sir2p in Saccharomyces cerevisiae

How duplicate genes provide genetic robustness remains an unresolved question. We have examined the duplicated histone deacetylases Sir2p and Hst1p in Saccharomyces cerevisiae and find that these paralogs with non-overlapping functions can provide genetic robustness against null mutations through a substitution mechanism. Hst1p is an NAD⁺-dependent histone deacetylase that acts with Sum1p to repress a subset of midsporulation genes. However, hst1Δ mutants show much weaker derepression of target loci than sum1Δ mutants. We show that this modest derepression of target loci in hst1Δ strains occurs in part because Sir2p substitutes for Hst1p. Sir2p contributes to repression of the midsporulation genes only in the absence of Hst1p and is recruited to target promoters by a physical interaction with the Sum1 complex. Furthermore, when Sir2p associates with the Sum1 complex, the complex continues to repress in a promoter-specific manner and does not spread. Our results imply that after the duplication, SIR2 and HST1 subfunctionalized. The single SIR2/HST1 gene from Kluyveromyces lactis, a closely related species that diverged prior to the duplication, can suppress an hst1Δ mutation in S. cerevisiae as well as interact with Sir4p in S. cerevisiae. In addition, the existence of two distinct protein interaction domains for the Sir and Sum1 complexes was revealed through the analysis of a chimeric Sir2–Hst1 molecule. Therefore, the ability of Sir2p to substitute for Hst1p probably results from a retained but reduced affinity for the Sum1 complex that is a consequence of subfunctionalization via the duplication, degeneration, and complementation mechanism. These results suggest that the evolutionary path of duplicate gene preservation may be an important indicator for the ability of duplicated genes to contribute to genetic robustness.

Gene duplication is an important force in evolution, as it provides a source of new genetic material. However, the mechanisms by which duplicated genes are retained and diverge are understudied at the experimental level. We have examined a pair of duplicated histone deacetylases Hst1p and Sir2p from baker's yeast, which are important for distinct types of gene repression. In this study, we show that before the duplication the ancestral histone deacetylase had both Hst1p- and Sir2p-like functions, and after the duplication Sir2p and Hst1p subfunctionalized, giving rise to two distinct proteins with normally nonoverlapping functions. Despite having partitioned the ancestral functions after the duplication, Sir2p can substitute for Hst1p in its absence by interacting with the normal partner of Hst1p. This study suggests that the evolutionary path of duplicate gene preservation may be an important indicator for the ability of duplicated genes to substitute for one another and hence protect the organism against deleterious mutations.

Frequent gain and loss of functional transcription factor binding sites

Cis-regulatory sequences are not always conserved across species. Divergence within cis-regulatory sequences may result from the evolution of species-specific patterns of gene expression or the flexible nature of the cis-regulatory code. The identification of functional divergence in cis-regulatory sequences is therefore important for both understanding the role of gene regulation in evolution and annotating regulatory elements. We have developed an evolutionary model to detect the loss of constraint on individual transcription factor binding sites (TFBSs). We find that a significant fraction of functionally constrained binding sites have been lost in a lineage-specific manner among three closely related yeast species. Binding site loss has previously been explained by turnover, where the concurrent gain and loss of a binding site maintains gene regulation. We estimate that nearly half of all loss events cannot be explained by binding site turnover. Recreating the mutations that led to binding site loss confirms that these sequence changes affect gene expression in some cases. We also estimate that there is a high rate of binding site gain, as more than half of experimentally identified S. cerevisiae binding sites are not conserved across species. The frequent gain and loss of TFBSs implies that cis-regulatory sequences are labile and, in the absence of turnover, may contribute to species-specific patterns of gene expression.

Research in the field of molecular evolution is focused on understanding the genetic basis of functional differences between species. Protein coding sequences have traditionally been the focus of these studies, as the genetic code enables a detailed study of the strength of selection acting on amino acid sequences. However, from the earliest cross-species sequence comparisons, it was clear that protein sequences among closely related species are too similar to explain the observed phenotypic diversity. This led to the hypothesis that the evolution of gene regulation has played a key role in generating diversity between species. The availability of numerous complete genome sequences has made it possible to begin testing this hypothesis. In this work, the authors use an evolutionary model to identify functional divergence within transcription factor binding sites, the core functional elements involved in gene regulation. Applying this model to the baker's yeast, Saccharomyces cerevisiae, and its three closest relatives, the authors find that a substantial fraction of the ancestral binding sites have been lost in a species-specific manner. In some cases the loss of the binding site creates gene expression differences that may be indicative of species-specific changes in gene regulation. This work provides a useful computational framework that will allow further study of the conservation of cis-regulatory sequences and their role in molecular evolution.

Arg-Pro-X-Ser/Thr is a consensus phosphoacceptor sequence for the meiosis-specific Ime2 protein kinase in Saccharomyces cerevisiae

Ime2 is a meiosis-specific protein kinase in Saccharomyces cerevisiae that is functionally related to cyclin-dependent kinase. Although Ime2 regulates multiple steps in meiosis, only a few of its substrates have been identified. Here we show that Ime2 phosphorylates Sum1, a repressor of meiotic gene transcription, on Thr-306. Ime2 protein kinase assays with Sum1 mutants and synthetic peptides define a consensus Arg-Pro-X-Ser/Thr motif that is required for efficient phosphorylation by Ime2. The carboxyl residue adjacent to the phosphoacceptor (+1 position) also influences the efficiency of Ime2 phosphorylation with alanine being a preferred residue. This information has predictive value in identifying new potential Ime2 targets as shown by the ability of Ime2 to phosphorylate Sgs1 and Gip1 in vitro and could be important in differentiating mitotic and meiotic regulatory pathways.

Protein binding microarrays for the characterization of DNA-protein interactions

A number of important cellular processes, such as transcriptional regulation, recombination, replication, repair, and DNA modification, are performed by DNA binding proteins. Of particular interest are transcription factors (TFs) which, through their sequence-specific interactions with DNA binding sites, modulate gene expression in a manner required for normal cellular growth and differentiation, and also for response to environmental stimuli. Despite their importance, the DNA binding specificities of most DNA binding proteins still remain unknown, since prior technologies aimed at identifying DNA-protein interactions have been laborious, not highly scalable, or have required limiting biological reagents. Recently a new DNA microarray-based technology, termed protein binding microarrays (PBMs), has been developed that allows rapid, high-throughput characterization of the in vitro DNA binding site sequence specificities of TFs, other DNA binding proteins, or synthetic compounds. DNA binding site data from PBMs combined with gene annotation data, comparative sequence analysis, and gene expression profiling, can be used to predict what genes are regulated by a given TF, what the functions are of a given TF and its predicted target genes, and how that TF may fit into the cell's transcriptional regulatory network.

Protein binding microarrays for the characterization of DNA-protein interactions

A number of important cellular processes, such as transcriptional regulation, recombination, replication, repair, and DNA modification, are performed by DNA binding proteins. Of particular interest are transcription factors (TFs) which, through their sequence-specific interactions with DNA binding sites, modulate gene expression in a manner required for normal cellular growth and differentiation, and also for response to environmental stimuli. Despite their importance, the DNA binding specificities of most DNA binding proteins still remain unknown, since prior technologies aimed at identifying DNA-protein interactions have been laborious, not highly scalable, or have required limiting biological reagents. Recently a new DNA microarray-based technology, termed protein binding microarrays (PBMs), has been developed that allows rapid, high-throughput characterization of the in vitro DNA binding site sequence specificities of TFs, other DNA binding proteins, or synthetic compounds. DNA binding site data from PBMs combined with gene annotation data, comparative sequence analysis, and gene expression profiling, can be used to predict what genes are regulated by a given TF, what the functions are of a given TF and its predicted target genes, and how that TF may fit into the cell's transcriptional regulatory network.

Principles of protein-DNA recognition revealed in the structural analysis of Ndt80-MSE DNA complexes

The Saccharomyces cerevisiae transcription factor Ndt80 selectively binds a DNA consensus sequence (the middle sporulation element [MSE]) to activate gene expression after the successful completion of meiotic recombination. Here we report the X-ray crystal structures of Ndt80 bound to ten distinct MSE variants. Comparison of these structures with the structure of Ndt80 bound to a consensus MSE reveals structural principles that determine the DNA binding specificity of this transcription factor. The 5' GC-rich end of the MSE contains distinct 5'-YpG-3' steps that are recognized by arginine side chains through a combination of hydrogen bonding and cation-pi interactions. The 3' AT-rich region is recognized via minor groove contacts that sterically exclude the N2 atom of GC base pairs. The conformation of the AT-rich region is fixed by interactions with the protein that favor recognition of poly(A)-poly(T) versus mixed AT sequences through an avoidance of major groove steric clashes at 5'-ApT-3' steps.

Evidence for ORC-dependent repression of budding yeast genes induced by starvation and other stresses

The highly conserved origin recognition complex (ORC) is required for repressing genes in the silent mating type loci of budding yeast. Here we report that at a non-permissive temperature, the temperature-sensitive orc2-1 mutation induces the expression of more than 500 genes, the majority of which are also induced during starvation of wild-type cells. Many genes induced by starvation or by the orc2-1 mutation are also induced by inactivation of proteins required for chromatin-mediated repression of transcription. Genes induced by the orc2-1 mutation, starvation, or inactivation of repressor proteins, map near ORC-binding loci significantly more frequently compared to all genes. Genes repressed by starvation map near ORC-binding sites less frequently compared to all genes, which suggests they have been evolutionarily excluded from regions of repressive chromatin near ORC-binding sites. Deletion of sequences containing ORC-binding sites near the DAL2 and DAL4 genes in the DAL gene cluster, which are induced by either the orc2-1 mutation or by starvation, constitutively activates these genes and abolishes their activation by the orc2-1 mutation. Our findings suggest a role for ORC in the repression of a large number of budding yeast genes induced by starvation or other aspects of a deleterious environment.

Non-random clustering of stress-related genes during evolution of the S. cerevisiae genome

Background

Coordinately regulated genes often physically cluster in eukaryotic genomes, for reasons that remain unclear.

Results

Here we provide evidence that many S. cerevisiae genes induced by starvation and other stresses reside in non-random clusters, where transcription of these genes is repressed in the absence of stress. Most genes essential for growth or for rapid, post-transcriptional responses to stress in cycling cells map between these gene clusters. Genes that are transcriptionally induced by stresses include a large fraction of rapidly evolving paralogues of duplicated genes that arose during an ancient whole genome duplication event. Many of these rapidly evolving paralogues have acquired new or more specialized functions that are less essential for growth. The slowly evolving paralogues of these genes are less likely to be transcriptionally repressed in the absence of stress, and are frequently essential for growth or for rapid stress responses that may require constitutive expression of these genes in cycling cells.

Conclusion

Our findings suggest that a fundamental organizing principle during evolution of the S. cerevisiae genome has been clustering of starvation and other stress-induced genes in chromosome regions that are transcriptionally repressed in the absence of stress, from which most genes essential for growth or rapid stress responses have been excluded. Chromatin-mediated repression of many stress-induced genes may have evolved since the whole genome duplication in parallel with functions for proteins encoded by these genes that are incompatible with growth. These functions likely provide fitness effects that escape detection in assays of reproductive capacity routinely employed to assess evolutionary fitness, or to identify genes that confer stress-resistance in cycling cells.

Checking Your Breaks: Surveillance Mechanisms of Meiotic Recombination

Numerous DNA double-strand breaks (DSBs) are introduced into the genome in the course of meiotic recombination. This poses a significant hazard to the genomic integrity of the cell. Studies in a number of organisms have unveiled the existence of surveillance mechanisms or checkpoints that couple the formation and repair of DSBs to cell cycle progression. Through these mechanisms, aberrant meiocytes are delayed in their meiotic progression, thereby facilitating repair of meiotic DSBs, or are culled through programmed cell death, thereby protecting the germline from aneuploidies that could lead to spontaneous abortions, birth defects and cancer predisposition in the offspring. Here we summarize recent progress in our understanding of these checkpoints. This review focuses on the surveillance mechanisms of the budding yeast S. cerevisiae, where the molecular details are best understood, but will frequently compare and contrast these mechanisms with observations in other organisms.

Gene expression from random libraries of yeast promoters

Genomewide techniques to assay gene expression and transcription factor binding are in widespread use, but are far from providing predictive rules for the function of regulatory DNA. To investigate more intensively the grammar rules for active regulatory sequence, we made libraries from random ligations of a very restricted set of sequences. Working with the yeast Saccharomyces cerevisiae, we developed a novel screen based on the sensitivity of ascospores lacking dityrosine to treatment with lytic enzymes. We tested two separate libraries built by random ligation of a single type of activator site either for a well-characterized sporulation factor, Ndt80, or for a new sporulation-specific regulatory site that we identified and several neutral spacer elements. This selective system achieved up to 1:10(4) enrichment of the artificial sequences that were active during sporulation, allowing a high-throughput analysis of large libraries of synthetic promoters. This is not practical with methods involving direct screening for expression, such as those based on fluorescent reporters. There were very few false positives, since active promoters always passed the screen when retested. The survival rate of our libraries containing roughly equal numbers of spacers and activators was a few percent that of libraries made from activators alone. The sequences of approximately 100 examples of active and inactive promoters could not be distinguished by simple binary rules; instead, the best model for the data was a linear regression fit of a quantitative measure of gene activity to multiple features of the regulatory sequence.

Analyses of SUM1-1-mediated long-range repression

In Saccharomyces cerevisiae, local repression is promoter specific and localized to a small region on the DNA, while silencing is promoter nonspecific, encompasses large domains of chromatin, and is stably inherited for multiple generations. Sum1p is a local repressor protein that mediates repression of meiosis-specific genes in mitotic cells while the Sir proteins are long-range repressors that stably silence genes at HML, HMR, and telomeres. The SUM1-1 mutation is a dominant neomorphic mutation that enables the mutant protein to be recruited to the HMR locus and repress genes, even in the absence of the Sir proteins. In this study we show that the mutation in Sum1-1p enabled it to spread, and the native HMR barrier blocked it from spreading. Thus, like the Sir proteins, Sum1-1p was a long-range repressor, but unlike the Sir proteins, Sum1-1p-mediated repression was more promoter specific, repressing certain genes better than others. Furthermore, repression mediated by Sum1-1p was not stably maintained or inherited and we therefore propose that Sum1-1p-mediated long-range repression is related but distinct from silencing.

Genome-wide identification of the regulatory targets of a transcription factor using biochemical characterization and computational genomic analysis

Background

A major challenge in computational genomics is the development of methodologies that allow accurate genome-wide prediction of the regulatory targets of a transcription factor. We present a method for target identification that combines experimental characterization of binding requirements with computational genomic analysis.

Results

Our method identified potential target genes of the transcription factor Ndt80, a key transcriptional regulator involved in yeast sporulation, using the combined information of binding affinity, positional distribution, and conservation of the binding sites across multiple species. We have also developed a mathematical approach to compute the false positive rate and the total number of targets in the genome based on the multiple selection criteria.

Conclusion

We have shown that combining biochemical characterization and computational genomic analysis leads to accurate identification of the genome-wide targets of a transcription factor. The method can be extended to other transcription factors and can complement other genomic approaches to transcriptional regulation.

Protein-DNA binding specificity predictions with structural models

Protein-DNA interactions play a central role in transcriptional regulation and other biological processes. Investigating the mechanism of binding affinity and specificity in protein-DNA complexes is thus an important goal. Here we develop a simple physical energy function, which uses electrostatics, solvation, hydrogen bonds and atom-packing terms to model direct readout and sequence-specific DNA conformational energy to model indirect readout of DNA sequence by the bound protein. The predictive capability of the model is tested against another model based only on the knowledge of the consensus sequence and the number of contacts between amino acids and DNA bases. Both models are used to carry out predictions of protein-DNA binding affinities which are then compared with experimental measurements. The nearly additive nature of protein-DNA interaction energies in our model allows us to construct position-specific weight matrices by computing base pair probabilities independently for each position in the binding site. Our approach is less data intensive than knowledge-based models of protein-DNA interactions, and is not limited to any specific family of transcription factors. However, native structures of protein-DNA complexes or their close homologs are required as input to the model. Use of homology modeling can significantly increase the extent of our approach, making it a useful tool for studying regulatory pathways in many organisms and cell types.

Explicit equilibrium modeling of transcription-factor binding and gene regulation

A computational model, GOMER, is presented that predicts transcription-factor binding and incorporates effects of cooperativity and competition.

We have developed a computational model that predicts the probability of transcription factor binding to any site in the genome. GOMER (generalizable occupancy model of expression regulation) calculates binding probabilities on the basis of position weight matrices, and incorporates the effects of cooperativity and competition by explicit calculation of coupled binding equilibria. GOMER can be used to test hypotheses regarding gene regulation that build upon this physically principled prediction of protein-DNA binding.

The Ama1-directed anaphase-promoting complex regulates the Smk1 mitogen-activated protein kinase during meiosis in yeast

Smk1 is a meiosis-specific MAPK homolog in Saccharomyces cerevisiae that regulates the postmeiotic program of spore formation. Similar to other MAPKs, it is activated via phosphorylation of the T-X-Y motif in its regulatory loop, but the signals controlling Smk1 activation have not been defined. Here we show that Ama1, a meiosis-specific activator of the anaphase-promoting complex/cyclosome (APC/C), promotes Smk1 activation during meiosis. A weakened allele of CDC28 suppresses the sporulation defect of an ama1 null strain and increases the activation state of Smk1. The function of Ama1 in regulating Smk1 is independent of the FEAR network, which promotes exit from mitosis and exit from meiosis I through the Cdc14 phosphatase. The data indicate that Cdc28 and Ama1 function in a pathway to trigger Smk1-dependent steps in spore morphogenesis. We propose that this novel mechanism for controlling MAPK activation plays a role in coupling the completion of meiosis II to gamete formation.

Control of replication initiation and heterochromatin formation in Saccharomyces cerevisiae by a regulator of meiotic gene expression

Heterochromatinization at the silent mating-type loci HMR and HML in Saccharomyces cerevisiae is achieved by targeting the Sir complex to these regions via a set of anchor proteins that bind to the silencers. Here, we have identified a novel heterochromatin-targeting factor for HML, the protein Sum1, a repressor of meiotic genes during vegetative growth. Sum1 bound both in vitro and in vivo to HML via a functional element within the HML-E silencer, and sum1Delta caused HML derepression. Significantly, Sum1 was also required for origin activity of HML-E, demonstrating a role of Sum1 in replication initiation. In a genome-wide search for Sum1-regulated origins, we identified a set of autonomous replicative sequences (ARS elements) that bound both the origin recognition complex and Sum1. Full initiation activity of these origins required Sum1, and their origin activity was decreased upon removal of the Sum1-binding site. Thus, Sum1 constitutes a novel global regulator of replication initiation in yeast.

Sum1p, the origin recognition complex, and the spreading of a promoter-specific repressor in Saccharomyces cerevisiae

In Saccharomyces cerevisiae, Sum1p is a promoter-specific repressor. A single amino acid change generates the mutant Sum1-1p, which causes regional silencing at new loci where wild-type Sum1p does not act. Thus, Sum1-1p is a model for understanding how the spreading of repressive chromatin is regulated. When wild-type Sum1p was targeted to a locus where mutant Sum1-1p spreads, wild-type Sum1p did not spread as efficiently as mutant Sum1-1p did, despite being in the same genomic context. Thus, the SUM1-1 mutation altered the ability of the protein to spread. The spreading of Sum1-1p required both an enzymatically active deacetylase, Hst1p, and the N-terminal tail of histone H4, consistent with the spreading of Sum1-1p involving sequential modification of and binding to histone tails, as observed for other silencing proteins. Furthermore, deletion of the N-terminal tail of H4 caused Sum1-1p to return to loci where wild-type Sum1p acts, consistent with the SUM1-1 mutation increasing the affinity of the protein for H4 tails. These results imply that the spreading of repressive chromatin proteins is regulated by their affinities for histone tails. Finally, this study uncovered a functional connection between wild-type Sum1p and the origin recognition complex, and this relationship also contributes to mutant Sum1-1p localization.

Birth of a metabolic gene cluster in yeast by adaptive gene relocation

Although most eukaryotic genomes lack operons, they contain some physical clusters of genes that are related in function despite being unrelated in sequence. How these clusters are formed during evolution is unknown. The DAL cluster is the largest metabolic gene cluster in yeast and consists of six adjacent genes encoding proteins that enable Saccharomyces cerevisiae to use allantoin as a nitrogen source. We show here that the DAL cluster was assembled, quite recently in evolutionary terms, through a set of genomic rearrangements that happened almost simultaneously. Six of the eight genes involved in allantoin degradation, which were previously scattered around the genome, became relocated to a single subtelomeric site in an ancestor of S. cerevisiae and Saccharomyces castellii. These genomic rearrangements coincided with a biochemical reorganization of the purine degradation pathway, which switched to importing allantoin instead of urate. This change eliminated urate oxidase, one of several oxygen-consuming enzymes that were lost by yeasts that can grow vigorously in anaerobic conditions. The DAL cluster is located in a domain of modified chromatin involving both H2A.Z histone exchange and Hst1-Sum1-mediated histone deacetylation, and it may be a coadapted gene complex formed by epistatic selection.

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.

Ime2p and Cdc28p: co-pilots driving meiotic development

Meiosis can be considered an elaboration of the cell division cycle in the sense that meiosis combines cell-cycle processes with programs specific to meiosis. Each phase of the cell division cycle is driven forward by cell-cycle kinases (Cdk) and coordinated with other phases of the cycle through checkpoint functions. Meiotic differentiation is also controlled by these two types of regulation; however, recent study in the budding yeast S. cerevisiae indicates that progression of meiosis is also controlled by a master regulator specific to meiosis, namely the Ime2p kinase. Below, I describe the overlapping roles of Ime2p and Cdk during meiosis in yeast and speculate on how these two kinases cooperate to drive the progression of meiosis.

Inference of combinatorial regulation in yeast transcriptional networks: a case study of sporulation

Decomposing transcriptional regulatory networks into functional modules and determining logical relations between them is the first step toward understanding transcriptional regulation at the system level. Modules based on analysis of genome-scale data can serve as the basis for inferring combinatorial regulation and for building mathematical models to quantitatively describe the behavior of the networks. We present here an algorithm called modem to identify target genes of a transcription factor (TF) from a single expression experiment, based on a joint probabilistic model for promoter sequence and gene expression data. We show how this method can facilitate the discovery of specific instances of combinatorial regulation and illustrate this for a specific case of transcriptional networks that regulate sporulation in the yeast Saccharomyces cerevisiae. Applying this method to analyze two crucial TFs in sporulation, Ndt80p and Sum1p, we were able to delineate their overlapping binding sites. We proposed a mechanistic model for the competitive regulation by the two TFs on a defined subset of sporulation genes. We show that this model accounts for the temporal control of the "middle" sporulation genes and suggest a similar regulatory arrangement can be found in developmental programs in higher organisms.

Analysis of the meiotic role of the mitochondrial ribosomal proteins Mrps17 and Mrpl37 in Saccharomyces cerevisiae

Sporulation in the yeast Saccharomyces cerevisiae is a complex and tightly regulated pathway that involves the induction of a large number of genes. We have identified MRPS17 in a cDNA library enriched for sporulation-specific genes. Homology searches show that the first one-third of Mrps17 has strong sequence similarity to bacterial S17 proteins, suggesting that Mrps17 is a potential mitochondrial ribosomal protein. This is further supported by the fact that mrps17Delta cells are respiratory-deficient and that a Mrps17-GFP fusion localizes to the mitochondria. We have confirmed by Northern blot analysis that both MRPS17 and MRPL37 are strongly induced during the middle stages of sporulation and that this induction is dependent on the presence of a middle sporulation element (MSE) in the promoters of these genes. Interestingly, we found that Mrps17 and Mrpl37, but not other mitochondrial ribosomal proteins, accumulate during the middle stages of sporulation. These results suggest that Mrps17 and Mrpl37 may have additional meiosis-specific roles.

Characterization of critical interactions between Ndt80 and MSE DNA defining a novel family of Ig-fold transcription factors

The Ndt80 protein of the yeast Saccharomyces cerevisiae is the founding member of a new sub-family of proteins in the Ig-fold superfamily of transcription factors. The crystal structure of Ndt80 bound to DNA shows that it makes contacts through several loops on one side of the protein that connect beta-strands which form the beta-sandwich fold common to proteins in this superfamily. However, the DNA-binding domain of Ndt80 is considerably larger than many other members of the Ig-fold superfamily and it appears to make a larger number of contacts with the DNA than these proteins. To determine the contribution of each of these contacts and to examine if the mechanism of Ndt80 DNA binding was similar to other members of the Ig-fold superfamily, amino acid substitutions were introduced at each residue that contacts the DNA and assayed for their effect on Ndt80 activity. Many of the mutations caused significant decreases in DNA-binding affinity and transcriptional activation. Several of these are in residues that are not found in other sub-families of Ig-fold proteins. These additional contacts are likely responsible for Ndt80's ability to bind DNA as a monomer while most other members require additional domains or cofactors to recognize their sites.

Purification and characterization of the DNA binding domain of Saccharomyces cerevisiae meiosis-specific transcription factor Ndt80

Ndt80 is a Saccharomyces cerevisiae meiosis-specific transcription factor responsible for promoting the stage-specific expression of a family of genes referred to as middle sporulation genes. Many members of this gene family are essential for the completion of meiotic chromosome segregation. Thus, Ndt80 is essential for the completion of meiosis. Ndt80 is highly regulated both transcriptionally and post-translationally. To facilitate biochemical analysis of Ndt80, we have expressed the DNA binding domain in Escherichia coli and purified the recombinant protein with an affinity chromatography procedure. In addition we have dissected the amino-terminus of Ndt80 to delimit the functional DNA binding domain. This analysis shows that the amino-terminal 40 amino-acids of Ndt80, although not essential for its DNA binding activity, do have an effect on its ability to bind specifically to its target DNA sequence. In addition, we show that the Ndt80 DNA binding domain can be phosphorylated by the meiosis-specific protein kinase Ime2 in vitro, but contrary to our initial hypothesis this phosphorylation does not significantly affect the affinity of Ndt80 for its target DNA sequence.

The activities of eukaryotic replication origins in chromatin

DNA replication initiates at chromosomal positions called replication origins. This review will focus on the activity, regulation and roles of replication origins in Saccharomyces cerevisiae. All eukaryotic cells, including S. cerevisiae, depend on the initiation (activity) of hundreds of replication origins during a single cell cycle for the duplication of their genomes. However, not all origins are identical. For example, there is a temporal order to origin activation with some origins firing early during the S-phase and some origins firing later. Recent studies provide evidence that posttranslational chromatin modifications, heterochromatin-binding proteins and nucleosome positioning can control the efficiency and/or timing of chromosomal origin activity in yeast. Many more origins exist than are necessary for efficient replication. The availability of excess replication origins leaves individual origins free to evolve distinct forms of regulation and/or roles in chromosomes beyond their fundamental role in DNA synthesis. We propose that some origins have acquired roles in controlling chromatin structure and/or gene expression. These roles are not linked obligatorily to replication origin activity per se, but instead exploit multi-subunit replication proteins with the potential to form context-dependent protein-protein interactions.

The Cdk-activating kinase Cak1p promotes meiotic S phase through Ime2p

CAK1 encodes an essential protein kinase in Saccharomyces cerevisiae that is required for activation of the Cdc28p Cdk. CAK1 also has several CDC28-independent functions that are unique to meiosis. The earliest of these functions is to induce S phase, which is regulated differently in meiosis than in mitosis. In mitosis, Cdc28p controls its own S-phase-promoting activity by signaling the destruction of its inhibitor, Sic1p. In meiosis, Sic1p destruction is signaled by the meiosis-specific Ime2p protein kinase. Our data show that Cak1p is required to activate Ime2p through a mechanism that requires threonine 242 and tyrosine 244 in Ime2p's activation loop. This activation promotes autophosphorylation and accumulation of multiply phosphorylated forms of Ime2p during meiotic development. Consistent with Cak1p's role in activating Ime2p, cells lacking Cak1p are deficient in degrading Sic1p. Deletion of SIC1 or overexpression of IME2 can partially suppress the S-phase defect in cak1 mutant cells, suggesting that Ime2p is a key target of Cak1p regulation. These data show that Cak1p is required for the destruction of Sic1p in meiosis, as in mitosis, but in meiosis, it functions through a sporulation-specific kinase.