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

A unified model for yeast transcript definition

Identifying genes in the genomic context is central to a cell's ability to interpret the genome. Yet, in general, the signals used to define eukaryotic genes are poorly described. Here, we derived simple classifiers that identify where transcription will initiate and terminate using nucleic acid sequence features detectable by the yeast cell, which we integrate into a Unified Model (UM) that models transcription as a whole. The cis-elements that denote where transcription initiates function primarily through nucleosome depletion, and, using a synthetic promoter system, we show that most of these elements are sufficient to initiate transcription in vivo. Hrp1 binding sites are the major characteristic of terminators; these binding sites are often clustered in terminator regions and can terminate transcription bidirectionally. The UM predicts global transcript structure by modeling transcription of the genome using a hidden Markov model whose emissions are the outputs of the initiation and termination classifiers. We validated the novel predictions of the UM with available RNA-seq data and tested it further by directly comparing the transcript structure predicted by the model to the transcription generated by the cell for synthetic DNA segments of random design. We show that the UM identifies transcription start sites more accurately than the initiation classifier alone, indicating that the relative arrangement of promoter and terminator elements influences their function. Our model presents a concrete description of how the cell defines transcript units, explains the existence of nongenic transcripts, and provides insight into genome evolution.

Mechanisms and functions of ATP-dependent chromatin-remodeling enzymes

Chromatin provides both a means to accommodate a large amount of genetic material in a small space and a means to package the same genetic material in different chromatin states. Transitions between chromatin states are enabled by chromatin-remodeling ATPases, which catalyze a diverse range of structural transformations. Biochemical evidence over the last two decades suggests that chromatin-remodeling activities may have emerged by adaptation of ancient DNA translocases to respond to specific features of chromatin. Here, we discuss such evidence and also relate mechanistic insights to our understanding of how chromatin-remodeling enzymes enable different in vivo processes.

Functional analyses of NSF1 in wine yeast using interconnected correlation clustering and molecular analyses

Analyzing time-course expression data captured in microarray datasets is a complex undertaking as the vast and complex data space is represented by a relatively low number of samples as compared to thousands of available genes. Here, we developed the Interdependent Correlation Clustering (ICC) method to analyze relationships that exist among genes conditioned on the expression of a specific target gene in microarray data. Based on Correlation Clustering, the ICC method analyzes a large set of correlation values related to gene expression profiles extracted from given microarray datasets. ICC can be applied to any microarray dataset and any target gene. We applied this method to microarray data generated from wine fermentations and selected NSF1, which encodes a C₂H₂ zinc finger-type transcription factor, as the target gene. The validity of the method was verified by accurate identifications of the previously known functional roles of NSF1. In addition, we identified and verified potential new functions for this gene; specifically, NSF1 is a negative regulator for the expression of sulfur metabolism genes, the nuclear localization of Nsf1 protein (Nsf1p) is controlled in a sulfur-dependent manner, and the transcription of NSF1 is regulated by Met4p, an important transcriptional activator of sulfur metabolism genes. The inter-disciplinary approach adopted here highlighted the accuracy and relevancy of the ICC method in mining for novel gene functions using complex microarray datasets with a limited number of samples.

Compositional and structural analysis of selected chromosomal domains from Saccharomyces cerevisiae

Chromatin is the template for replication and transcription in the eukaryotic nucleus, which needs to be defined in composition and structure before these processes can be fully understood. We report an isolation protocol for the targeted purification of specific genomic regions in their native chromatin context from Saccharomyces cerevisiae. Subdomains of the multicopy ribosomal DNA locus containing transcription units of RNA polymerases I, II or III or an autonomous replication sequence were independently purified in sufficient amounts and purity to analyze protein composition and histone modifications by mass spectrometry. We present and discuss the proteomic data sets obtained for chromatin in different functional states. The native chromatin was further amenable to electron microscopy analysis yielding information about nucleosome occupancy and positioning at the single-molecule level. We also provide evidence that chromatin from virtually every single copy genomic locus of interest can be purified and analyzed by this technique.

The BAH domain of Rsc2 is a histone H3 binding domain

Bromo-adjacent homology (BAH) domains are commonly found in chromatin-associated proteins and fall into two classes; Remodels the Structure of Chromatin (RSC)-like or Sir3-like. Although Sir3-like BAH domains bind nucleosomes, the binding partners of RSC-like BAH domains are currently unknown. The Rsc2 subunit of the RSC chromatin remodeling complex contains an RSC-like BAH domain and, like the Sir3-like BAH domains, we find Rsc2 BAH also interacts with nucleosomes. However, unlike Sir3-like BAH domains, we find that Rsc2 BAH can bind to recombinant purified H3 in vitro, suggesting that the mechanism of nucleosome binding is not conserved. To gain insight into the Rsc2 BAH domain, we determined its crystal structure at 2.4 Å resolution. We find that it differs substantially from Sir3-like BAH domains and lacks the motifs in these domains known to be critical for making contacts with histones. We then go on to identify a novel motif in Rsc2 BAH that is critical for efficient H3 binding in vitro and show that mutation of this motif results in defective Rsc2 function in vivo. Moreover, we find this interaction is conserved across Rsc2-related proteins. These data uncover a binding target of the Rsc2 family of BAH domains and identify a novel motif that mediates this interaction.

Gene promoters dictate histone occupancy within genes

Spt6 is a transcriptional elongation factor and histone chaperone that reassembles transcribed chromatin. Genome-wide H3 mapping showed that Spt6 preferentially maintains nucleosomes within the first 500 bases of genes and helps define nucleosome-depleted regions in 5' and 3' flanking sequences. In Spt6-depleted cells, H3 loss at 5' ends correlates with reduced pol II density suggesting enhanced transcription elongation. Consistent with its 'Suppressor of Ty' (Spt) phenotype, Spt6 inactivation caused localized H3 eviction over 1-2 nucleosomes at 5' ends of Ty elements. H3 displacement differed between genes driven by promoters with 'open'/DPN and 'closed'/OPN chromatin conformations with similar pol II densities. More eviction occurred on genes with 'closed' promoters, associated with 'noisy' transcription. Moreover, swapping of 'open' and 'closed' promoters showed that they can specify distinct downstream patterns of histone eviction/deposition. These observations suggest a novel function for promoters in dictating histone dynamics within genes possibly through effects on transcriptional bursting or elongation rate.

Post-zygotic and inter-individual structural genetic variation in a presumptive enhancer element of the locus between the IL10Rβ and IFNAR1 genes

Although historically considered as junk-DNA, tandemly repeated sequence motifs can affect human phenotype. For example, variable number tandem repeats (VNTR) with embedded enhancers have been shown to regulate gene transcription. The post-zygotic variation is the presence of genetically distinct populations of cells in an individual derived from a single zygote, and this is an understudied aspect of genome biology. We report somatically variable VNTR with sequence properties of an enhancer, located upstream of IFNAR1. Initially, SNP genotyping of 63 monozygotic twin pairs and multiple tissues from 21 breast cancer patients suggested a frequent post-zygotic mosaicism. The VNTR displayed a repeated 32 bp core motif in the center of the repeat, which was flanked by similar variable motifs. A total of 14 alleles were characterized based on combinations of segments, which showed post-zygotic and inter-individual variation, with up to 6 alleles in a single subject. Somatic variation occurred in ∼24% of cases. In this hypervariable region, we found a clustering of transcription factor binding sites with strongest sequence similarity to mouse Foxg1 transcription factor binding motif. This study describes a VNTR with sequence properties of an enhancer that displays post-zygotic and inter-individual genetic variation. This element is within a locus containing four related cytokine receptors: IFNAR2, IL10Rβ, IFNAR1 and IFNGR2, and we hypothesize that it might function in transcriptional regulation of several genes in this cluster. Our findings add another level of complexity to the variation among VNTR-based enhancers. Further work may unveil the normal function of this VNTR in transcriptional control and its possible involvement in diseases connected with these receptors, such as autoimmune conditions and cancer.

Dissociation of the H3K36 demethylase Rph1 from chromatin mediates derepression of environmental stress-response genes under genotoxic stress in Saccharomyces cerevisiae

The H3K36 demethylase Rph1 is a transcriptional repressor for stress-responsive genes in yeast. Rph1-mediated transcriptional repression is relieved by phosphorylation of Rph1, reduced Rph1 level, and dissociation of Rph1 from chromatin with genotoxic stress. Rph1 may function as a regulatory node in different stress-signaling pathways.

Cells respond to environmental signals by altering gene expression through transcription factors. Rph1 is a histone demethylase containing a Jumonji C (JmjC) domain and belongs to the C₂H₂ zinc-finger protein family. Here we investigate the regulatory network of Rph1 in yeast by expression microarray analysis. More than 75% of Rph1-regulated genes showed increased expression in the rph1-deletion mutant, suggesting that Rph1 is mainly a transcriptional repressor. The binding motif 5′-CCCCTWA-3′, which resembles the stress response element, is overrepresented in the promoters of Rph1-repressed genes. A significant proportion of Rph1-regulated genes respond to DNA damage and environmental stress. Rph1 is a labile protein, and Rad53 negatively modulates Rph1 protein level. We find that the JmjN domain is important in maintaining protein stability and the repressive effect of Rph1. Rph1 is directly associated with the promoter region of targeted genes and dissociated from chromatin before transcriptional derepression on DNA damage and oxidative stress. Of interest, the master stress-activated regulator Msn2 also regulates a subset of Rph1-repressed genes under oxidative stress. Our findings confirm the regulatory role of Rph1 as a transcriptional repressor and reveal that Rph1 might be a regulatory node connecting different signaling pathways responding to environmental stresses.

Structure of the transcriptional network controlling white-opaque switching in Candida albicans

The human fungal pathogen Candida albicans can switch between two phenotypic cell types, termed 'white' and 'opaque'. Both cell types are heritable for many generations, and the switch between the two types occurs epigenetically, that is, without a change in the primary DNA sequence of the genome. Previous work identified six key transcriptional regulators important for white-opaque switching: Wor1, Wor2, Wor3, Czf1, Efg1, and Ahr1. In this work, we describe the structure of the transcriptional network that specifies the white and opaque cell types and governs the ability to switch between them. In particular, we use a combination of genome-wide chromatin immunoprecipitation, gene expression profiling, and microfluidics-based DNA binding experiments to determine the direct and indirect regulatory interactions that form the switch network. The six regulators are arranged together in a complex, interlocking network with many seemingly redundant and overlapping connections. We propose that the structure (or topology) of this network is responsible for the epigenetic maintenance of the white and opaque states, the switching between them, and the specialized properties of each state.

Nucleosome free regions in yeast promoters result from competitive binding of transcription factors that interact with chromatin modifiers

Because DNA packaging in nucleosomes modulates its accessibility to transcription factors (TFs), unraveling the causal determinants of nucleosome positioning is of great importance to understanding gene regulation. Although there is evidence that intrinsic sequence specificity contributes to nucleosome positioning, the extent to which other factors contribute to nucleosome positioning is currently highly debated. Here we obtained both in vivo and in vitro reference maps of positions that are either consistently covered or free of nucleosomes across multiple experimental data-sets in Saccharomyces cerevisiae. We then systematically quantified the contribution of TF binding to nucleosome positiong using a rigorous statistical mechanics model in which TFs compete with nucleosomes for binding DNA. Our results reconcile previous seemingly conflicting results on the determinants of nucleosome positioning and provide a quantitative explanation for the difference between in vivo and in vitro positioning. On a genome-wide scale, nucleosome positioning is dominated by the phasing of nucleosome arrays over gene bodies, and their positioning is mainly determined by the intrinsic sequence preferences of nucleosomes. In contrast, larger nucleosome free regions in promoters, which likely have a much more significant impact on gene expression, are determined mainly by TF binding. Interestingly, of the 158 yeast TFs included in our modeling, we find that only 10–20 significantly contribute to inducing nucleosome-free regions, and these TFs are highly enriched for having direct interations with chromatin remodelers. Together our results imply that nucleosome free regions in yeast promoters results from the binding of a specific class of TFs that recruit chromatin remodelers.

The DNA of all eukaryotic organisms is packaged into nucleosomes, which cover roughly of the genome. As nucleosome positioning profoundly affects DNA accessibility to other DNA binding proteins such as transcription factors (TFs), it plays an important role in transcription regulation. However, to what extent nucleosome positioning is guided by intrinsic DNA sequence preferences of nucleosomes, and to what extent other DNA binding factors play a role, is currently highly debated. Here we use a rigorous biophysical model to systematically study the relative contributions of intrinsic sequence preferences and competitive binding of TFs to nucleosome positioning in yeast. We find that, on the one hand, the phasing of the many small spacers within dense nucleosome arrays that cover gene bodies are mainly determined by intrinsic sequence preferences. On the other hand, larger nucleosome free regions (NFRs) in promoters are explained predominantly by TF binding. Strikingly, we find that only 10–20 TFs make a significant contribution to explaining NFRs, and these TFs are highly enriched for directly interacting with chromatin modifiers. Thus, the picture that emerges is that binding by a specific class of TFs recruits chromatin modifiers which mediate local nucleosome expulsion.

Inferring gene expression from ribosomal promoter sequences, a crowdsourcing approach

The Gene Promoter Expression Prediction challenge consisted of predicting gene expression from promoter sequences in a previously unknown experimentally generated data set. The challenge was presented to the community in the framework of the sixth Dialogue for Reverse Engineering Assessments and Methods (DREAM6), a community effort to evaluate the status of systems biology modeling methodologies. Nucleotide-specific promoter activity was obtained by measuring fluorescence from promoter sequences fused upstream of a gene for yellow fluorescence protein and inserted in the same genomic site of yeast Saccharomyces cerevisiae. Twenty-one teams submitted results predicting the expression levels of 53 different promoters from yeast ribosomal protein genes. Analysis of participant predictions shows that accurate values for low-expressed and mutated promoters were difficult to obtain, although in the latter case, only when the mutation induced a large change in promoter activity compared to the wild-type sequence. As in previous DREAM challenges, we found that aggregation of participant predictions provided robust results, but did not fare better than the three best algorithms. Finally, this study not only provides a benchmark for the assessment of methods predicting activity of a specific set of promoters from their sequence, but it also shows that the top performing algorithm, which used machine-learning approaches, can be improved by the addition of biological features such as transcription factor binding sites.

A highly redundant gene network controls assembly of the outer spore wall in S. cerevisiae

The spore wall of Saccharomyces cerevisiae is a multilaminar extracellular structure that is formed de novo in the course of sporulation. The outer layers of the spore wall provide spores with resistance to a wide variety of environmental stresses. The major components of the outer spore wall are the polysaccharide chitosan and a polymer formed from the di-amino acid dityrosine. Though the synthesis and export pathways for dityrosine have been described, genes directly involved in dityrosine polymerization and incorporation into the spore wall have not been identified. A synthetic gene array approach to identify new genes involved in outer spore wall synthesis revealed an interconnected network influencing dityrosine assembly. This network is highly redundant both for genes of different activities that compensate for the loss of each other and for related genes of overlapping activity. Several of the genes in this network have paralogs in the yeast genome and deletion of entire paralog sets is sufficient to severely reduce dityrosine fluorescence. Solid-state NMR analysis of partially purified outer spore walls identifies a novel component in spore walls from wild type that is absent in some of the paralog set mutants. Localization of gene products identified in the screen reveals an unexpected role for lipid droplets in outer spore wall formation.

The cell wall of fungi is a complex extracellular matrix and an important target for antifungal drugs. Assembly of the wall during spore formation in baker's yeast is a useful model for fungal wall development. The outermost layers of the spore wall are composed of a polymer of dityrosine connected to an underlying polysaccharide layer. The assembly pathway of this dityrosine polymer is not known. Using a genetic approach we reveal a network of genes that function redundantly to control dityrosine layer synthesis. Solid state NMR analysis of spore walls from wild-type and mutant cells reveals a novel constituent of the spore wall that may link the dityrosine to the underlying polysaccharides and a role for lipid droplets in the incorporation of this new component into the spore wall.

Interplays between ATM/Tel1 and ATR/Mec1 in sensing and signaling DNA double-strand breaks

DNA double-strand breaks (DSBs) are highly hazardous for genome integrity because they have the potential to cause mutations, chromosomal rearrangements and genomic instability. The cellular response to DSBs is orchestrated by signal transduction pathways, known as DNA damage checkpoints, which are conserved from yeasts to humans. These pathways can sense DNA damage and transduce this information to specific cellular targets, which in turn regulate cell cycle transitions and DNA repair. The mammalian protein kinases ATM and ATR, as well as their budding yeast corresponding orthologs Tel1 and Mec1, act as master regulators of the checkpoint response to DSBs. Here, we review the early steps of DSB processing and the role of DNA-end structures in activating ATM/Tel1 and ATR/Mec1 in an orderly and reciprocal manner.

ATAF1 transcription factor directly regulates abscisic acid biosynthetic gene NCED3 in Arabidopsis thaliana

ATAF1, an Arabidopsis thaliana NAC transcription factor, plays important roles in plant adaptation to environmental stress and development. To search for ATAF1 target genes, we used protein binding microarrays and chromatin-immunoprecipitation (ChIP). This identified T[A,C,G]CGT[A,G] and TT[A,C,G]CGT as ATAF1 consensus binding sequences. Co-expression analysis across publicly available microarray experiments identified 25 genes co-expressed with ATAF1. The promoter regions of ATAF1 co-expressors were significantly enriched for ATAF1 binding sites, and TTGCGTA was identified in the promoter of the key abscisic acid (ABA) phytohormone biosynthetic gene NCED3. ChIP-qPCR and expression analysis showed that ATAF1 binding to the NCED3 promoter correlated with increased NCED3 expression and ABA hormone levels. These results indicate that ATAF1 regulates ABA biosynthesis.

The yeast Snt2 protein coordinates the transcriptional response to hydrogen peroxide-mediated oxidative stress

Regulation of gene expression is a vital part of the cellular stress response, yet the full set of proteins that orchestrate this regulation remains unknown. Snt2 is a Saccharomyces cerevisiae protein whose function has not been well characterized that was recently shown to associate with Ecm5 and the Rpd3 deacetylase. Here, we confirm that Snt2, Ecm5, and Rpd3 physically associate. We then demonstrate that cells lacking Rpd3 or Snt2 are resistant to hydrogen peroxide (H2O2)-mediated oxidative stress and use chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) to show that Snt2 and Ecm5 recruit Rpd3 to a small number of promoters and in response to H2O2, colocalize independently of Rpd3 to the promoters of stress response genes. By integrating ChIP-seq and expression analyses, we identify target genes that require Snt2 for proper expression after H2O2. Finally, we show that cells lacking Snt2 are also resistant to nutrient stress imparted by the TOR (target of rapamycin) pathway inhibitor rapamycin and identify a common set of genes targeted by Snt2 and Ecm5 in response to both H2O2 and rapamycin. Our results establish a function for Snt2 in regulating transcription in response to oxidative stress and suggest Snt2 may also function in multiple stress pathways.

Nucleosome maps of the human cytomegalovirus genome reveal a temporal switch in chromatin organization linked to a major IE protein

Human CMV (hCMV) establishes lifelong infections in most of us, causing developmental defects in human embryos and life-threatening disease in immunocompromised individuals. During productive infection, the viral >230,000-bp dsDNA genome is expressed widely and in a temporal cascade. The hCMV genome does not carry histones when encapsidated but has been proposed to form nucleosomes after release into the host cell nucleus. Here, we present hCMV genome-wide nucleosome occupancy and nascent transcript maps during infection of permissive human primary cells. We show that nucleosomes occupy nuclear viral DNA in a nonrandom and highly predictable fashion. At early times of infection, nucleosomes associate with the hCMV genome largely according to their intrinsic DNA sequence preferences, indicating that initial nucleosome formation is genetically encoded in the virus. However, as infection proceeds to the late phase, nucleosomes redistribute extensively to establish patterns mostly determined by nongenetic factors. We propose that these factors include key regulators of viral gene expression encoded at the hCMV major immediate-early (IE) locus. Indeed, mutant virus genomes deficient for IE1 expression exhibit globally increased nucleosome loads and reduced nucleosome dynamics compared with WT genomes. The temporal nucleosome occupancy differences between IE1-deficient and WT viruses correlate inversely with changes in the pattern of viral nascent and total transcript accumulation. These results provide a framework of spatial and temporal nucleosome organization across the genome of a major human pathogen and suggest that an hCMV major IE protein governs overall viral chromatin structure and function.

DNA-binding specificity changes in the evolution of forkhead transcription factors

The evolution of transcriptional regulatory networks entails the expansion and diversification of transcription factor (TF) families. The forkhead family of TFs, defined by a highly conserved winged helix DNA-binding domain (DBD), has diverged into dozens of subfamilies in animals, fungi, and related protists. We have used a combination of maximum-likelihood phylogenetic inference and independent, comprehensive functional assays of DNA-binding capacity to explore the evolution of DNA-binding specificity within the forkhead family. We present converging evidence that similar alternative sequence preferences have arisen repeatedly and independently in the course of forkhead evolution. The vast majority of DNA-binding specificity changes we observed are not explained by alterations in the known DNA-contacting amino acid residues conferring specificity for canonical forkhead binding sites. Intriguingly, we have found forkhead DBDs that retain the ability to bind very specifically to two completely distinct DNA sequence motifs. We propose an alternate specificity-determining mechanism whereby conformational rearrangements of the DBD broaden the spectrum of sequence motifs that a TF can recognize. DNA-binding bispecificity suggests a previously undescribed source of modularity and flexibility in gene regulation and may play an important role in the evolution of transcriptional regulatory networks.

Cell cycle-dependent transcription factors control the expression of yeast telomerase RNA

Telomerase is a specialized ribonucleoprotein that adds repeated DNA sequences to the ends of eukaryotic chromosomes to preserve genome integrity. Some secondary structure features of the telomerase RNA are very well conserved, and it serves as a central scaffold for the binding of associated proteins. The Saccharomyces cerevisiae telomerase RNA, TLC1, is found in very low copy number in the cell and is the limiting component of the known telomerase holoenzyme constituents. The reasons for this low abundance are unclear, but given that the RNA is very stable, transcriptional control mechanisms must be extremely important. Here we define the sequences forming the TLC1 promoter and identify the elements required for its low expression level, including enhancer and repressor elements. Within an enhancer element, we found consensus sites for Mbp1/Swi4 association, and chromatin immunoprecipitation (ChIP) assays confirmed the binding of Mbp1 and Swi4 to these sites of the TLC1 promoter. Furthermore, the enhancer element conferred cell cycle-dependent regulation to a reporter gene, and mutations in the Mbp1/Swi4 binding sites affected the levels of telomerase RNA and telomere length. Finally, ChIP experiments using a TLC1 RNA-binding protein as target showed cell cycle-dependent transcription of the TLC1 gene. These results indicate that the budding yeast TLC1 RNA is transcribed in a cell cycle-dependent fashion late in G1 and may be part of the S phase-regulated group of genes involved in DNA replication.

Functional improvement of Saccharomyces cerevisiae to reduce volatile acidity in wine

Control of volatile acidity (VA) is a major issue for wine quality. In this study, we investigated the production of VA by a deletion mutant of the fermentation stress response gene AAF1 in the budding yeast Saccharomyces cerevisiae. Fermentations were carried out in commercial Chardonnay grape must to mimic industrial wine-making conditions. We demonstrated that a wine yeast strain deleted for AAF1 reduced acetic acid levels in wine by up to 39.2% without increasing the acetaldehyde levels, revealing a potential for industrial application. Deletion of the cytosolic aldehyde dehydrogenase gene ALD6 also reduced acetic acid levels dramatically, but increased the acetaldehyde levels by 41.4%, which is not desired by the wine industry. By comparison, ALD4 and the AAF1 paralog RSF2 had no effects on acetic acid production in wine. Deletion of AAF1 was detrimental to the growth of ald6Δ and ald4Δald6Δ mutants, but had no effect on acetic acid production. Overexpression of AAF1 dramatically increased acetic acid levels in wine in an Ald6p-dependent manner, indicating that Aaf1p regulates acetic acid production mainly via Ald6p. Overexpression of AAF1 in an ald4Δald6Δ strain produced significantly more acetic acid in wine than the ald4Δald6Δ mutant, suggesting that Aaf1p may also regulate acetic acid synthesis independently of Ald4p and Ald6p.

UpSETing chromatin during non-coding RNA production

The packaging of eukaryotic DNA into nucleosomal arrays permits cells to tightly regulate and fine-tune gene expression. The ordered disassembly and reassembly of these nucleosomes allows RNA polymerase II (RNAPII) conditional access to the underlying DNA sequences. Disruption of nucleosome reassembly following RNAPII passage results in spurious transcription initiation events, leading to the production of non-coding RNA (ncRNA). We review the molecular mechanisms involved in the suppression of these cryptic initiation events and discuss the role played by ncRNAs in regulating gene expression.

Mitotic exit and separation of mother and daughter cells

Productive cell proliferation involves efficient and accurate splitting of the dividing cell into two separate entities. This orderly process reflects coordination of diverse cytological events by regulatory systems that drive the cell from mitosis into G1. In the budding yeast Saccharomyces cerevisiae, separation of mother and daughter cells involves coordinated actomyosin ring contraction and septum synthesis, followed by septum destruction. These events occur in precise and rapid sequence once chromosomes are segregated and are linked with spindle organization and mitotic progress by intricate cell cycle control machinery. Additionally, critical paarts of the mother/daughter separation process are asymmetric, reflecting a form of fate specification that occurs in every cell division. This chapter describes central events of budding yeast cell separation, as well as the control pathways that integrate them and link them with the cell cycle.

Dynamic nucleosome organization at hox promoters during zebrafish embryogenesis

Nucleosome organization at promoter regions plays an important role in regulating gene activity. Genome-wide studies in yeast, flies, worms, mammalian embryonic stem cells and transformed cell lines have found well-positioned nucleosomes flanking a nucleosome depleted region (NDR) at transcription start sites. This nucleosome arrangement depends on DNA sequence (cis-elements) as well as DNA binding factors and ATP-dependent chromatin modifiers (trans-factors). However, little is understood about how the nascent embryonic genome positions nucleosomes during development. This is particularly intriguing since the embryonic genome must undergo a broad reprogramming event upon fusion of sperm and oocyte. Using four stages of early embryonic zebrafish development, we map nucleosome positions at the promoter region of 37 zebrafish hox genes. We find that nucleosome arrangement at the hox promoters is a progressive process that takes place over several stages. At stages immediately after fertilization, nucleosomes appear to be largely disordered at hox promoter regions. At stages after activation of the embryonic genome, nucleosomes are detectable at hox promoters, with positions becoming more uniform and more highly occupied. Since the genomic sequence is invariant during embryogenesis, this progressive change in nucleosome arrangement suggests that trans-factors play an important role in organizing nucleosomes during embryogenesis. Separating hox genes into expressed and non-expressed groups shows that expressed promoters have better positioned and occupied nucleosomes, as well as distinct NDRs, than non-expressed promoters. Finally, by blocking the retinoic acid-signaling pathway, we disrupt early hox gene transcription, but observe no effect on nucleosome positions, suggesting that active hox transcription is not a driving force behind the arrangement of nucleosomes at the promoters of hox genes during early development.

Determinants of nucleosome positioning

Nucleosome positioning is critical for gene expression and most DNA-related processes. Here we review the dominant patterns of nucleosome positioning that have been observed and summarize the current understanding of their underlying determinants. The genome-wide pattern of nucleosome positioning is determined by the combination of DNA sequence, ATP-dependent nucleosome remodeling enzymes and transcription factors that include activators, components of the preinitiation complex and elongating RNA polymerase II. These determinants influence each other such that the resulting nucleosome positioning patterns are likely to differ among genes and among cells in a population, with consequent effects on gene expression.

A compendium of nucleosome and transcript profiles reveals determinants of chromatin architecture and transcription

Nucleosomes in all eukaryotes examined to date adopt a characteristic architecture within genes and play fundamental roles in regulating transcription, yet the identity and precise roles of many of the trans-acting factors responsible for the establishment and maintenance of this organization remain to be identified. We profiled a compendium of 50 yeast strains carrying conditional alleles or complete deletions of genes involved in transcriptional regulation, histone biology, and chromatin remodeling, as well as compounds that target transcription and histone deacetylases, to assess their respective roles in nucleosome positioning and transcription. We find that nucleosome patterning in genes is affected by many factors, including the CAF-1 complex, Spt10, and Spt21, in addition to previously reported remodeler ATPases and histone chaperones. Disruption of these factors or reductions in histone levels led genic nucleosomes to assume positions more consistent with their intrinsic sequence preferences, with pronounced and specific shifts of the +1 nucleosome relative to the transcription start site. These shifts of +1 nucleosomes appear to have functional consequences, as several affected genes in Ino80 mutants exhibited altered expression responses. Our parallel expression profiling compendium revealed extensive transcription changes in intergenic and antisense regions, most of which occur in regions with altered nucleosome occupancy and positioning. We show that the nucleosome-excluding transcription factors Reb1, Abf1, Tbf1, and Rsc3 suppress cryptic transcripts at their target promoters, while a combined analysis of nucleosome and expression profiles identified 36 novel transcripts that are normally repressed by Tup1/Cyc8. Our data confirm and extend the roles of chromatin remodelers and chaperones as major determinants of genic nucleosome positioning, and these data provide a valuable resource for future studies.

The genome in eukaryotic cells is packaged into nucleosomes, which play critical roles in regulating where and when different genes are expressed. For example, nucleosomes can physically block access of transcription factor to sites on DNA or direct regulatory proteins to DNA. Consistent with these roles, nucleosomes assume a stereotypical pattern around genes: they are depleted at the promoter region that marks the start of genes and assume a regularly spaced array within genes. To identify factors involved in this organization, we generated high-resolution nucleosome and transcriptome maps for 50 loss-of-function mutants with known or suspected roles in nucleosome biology in budding yeast. We show that nucleosome organization is determined by the combined effects of many factors that often exert opposing forces on nucleosomes. We further demonstrate that specific nucleosomes can be positioned independently within genes and that repositioning of nucleosomes at the start of genes may affect expression of these genes in response to environmental stimuli. Data mining of this extensive resource allowed us to show that general transcription factors act as insulators at diverging promoters to prevent the formation of cryptic transcripts, and also revealed 36 novel transcripts regulated by the Tup1/Cyc8 complex.

A role for the nucleoporin Nup170p in chromatin structure and gene silencing

Embedded in the nuclear envelope, nuclear pore complexes (NPCs) not only regulate nuclear transport but also interface with transcriptionally active euchromatin, largely silenced heterochromatin, as well as the boundaries between these regions. It is unclear what functional role NPCs play in establishing or maintaining these distinct chromatin domains. We report that the yeast NPC protein Nup170p interacts with regions of the genome that contain ribosomal protein and subtelomeric genes, where it functions in nucleosome positioning and as a repressor of transcription. We show that the role of Nup170p in subtelomeric gene silencing is linked to its association with the RSC chromatin-remodeling complex and the silencing factor Sir4p, and that the binding of Nup170p and Sir4p to subtelomeric chromatin is cooperative and necessary for the association of telomeres with the nuclear envelope. Our results establish the NPC as an active participant in silencing and the formation of peripheral heterochromatin.

Bayesian hierarchical model of protein-binding microarray k-mer data reduces noise and identifies transcription factor subclasses and preferred k-mers

MOTIVATION

Sequence-specific transcription factors (TFs) regulate the expression of their target genes through interactions with specific DNA-binding sites in the genome. Data on TF-DNA binding specificities are essential for understanding how regulatory specificity is achieved.

RESULTS

Numerous studies have used universal protein-binding microarray (PBM) technology to determine the in vitro binding specificities of hundreds of TFs for all possible 8 bp sequences (8mers). We have developed a Bayesian analysis of variance (ANOVA) model that decomposes these 8mer data into background noise, TF familywise effects and effects due to the particular TF. Adjusting for background noise improves PBM data quality and concordance with in vivo TF binding data. Moreover, our model provides simultaneous identification of TF subclasses and their shared sequence preferences, and also of 8mers bound preferentially by individual members of TF subclasses. Such results may aid in deciphering cis-regulatory codes and determinants of protein-DNA binding specificity.

AVAILABILITY AND IMPLEMENTATION

Source code, compiled code and R and Python scripts are available from http://thebrain.bwh.harvard.edu/hierarchicalANOVA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Differential binding of the related transcription factors Pho4 and Cbf1 can tune the sensitivity of promoters to different levels of an induction signal

Transcription factors that belong to the same family typically have similar, but not identical, binding specificities. As such, they can be expected to compete differentially for binding to different variants of their binding sites. Pho4 is a yeast factor whose nuclear concentration is up-regulated in low phosphate, while the related factor, Cbf1, is constitutively expressed. We constructed 16 GFP-reporter genes containing all palindromic variants of the motif NNCACGTGNN, and determined their activities at a range of phosphate concentrations. Pho4 affinity did not explain expression data well except under fully induced conditions. However, reporter activity was quantitatively well explained under all conditions by a model in which Cbf1 itself has modest activating activity, and Pho4 and Cbf1 compete with one another. Chromatin immunoprecipitation and computational analyses of natural Pho4 target genes, along with the activities of the reporter constructs, indicates that genes differ in their sensitivity to intermediate induction signals in part because of differences in their affinity for Cbf1. The induction sensitivity of both natural Pho4 target genes and reporter genes was well explained only by a model that assumes a role for Cbf1 in remodeling chromatin. Our analyses highlight the importance of taking into account the activities of related transcription factors in explaining system-wide gene expression data.

DNA looping facilitates targeting of a chromatin remodeling enzyme

ATP-dependent chromatin remodeling enzymes are highly abundant and play pivotal roles regulating DNA-dependent processes. The mechanisms by which they are targeted to specific loci have not been well understood on a genome-wide scale. Here, we present evidence that a major targeting mechanism for the Isw2 chromatin remodeling enzyme to specific genomic loci is through sequence-specific transcription factor (TF)-dependent recruitment. Unexpectedly, Isw2 is recruited in a TF-dependent fashion to a large number of loci without TF binding sites. Using the 3C assay, we show that Isw2 can be targeted by Ume6- and TFIIB-dependent DNA looping. These results identify DNA looping as a mechanism for the recruitment of a chromatin remodeling enzyme and define a function for DNA looping. We also present evidence suggesting that Ume6-dependent DNA looping is involved in chromatin remodeling and transcriptional repression, revealing a mechanism by which the three-dimensional folding of chromatin affects DNA-dependent processes.

Global analysis of Drosophila Cys₂-His₂ zinc finger proteins reveals a multitude of novel recognition motifs and binding determinants

Cys₂-His₂ zinc finger proteins (ZFPs) are the largest group of transcription factors in higher metazoans. A complete characterization of these ZFPs and their associated target sequences is pivotal to fully annotate transcriptional regulatory networks in metazoan genomes. As a first step in this process, we have characterized the DNA-binding specificities of 129 zinc finger sets from Drosophila using a bacterial one-hybrid system. This data set contains the DNA-binding specificities for at least one encoded ZFP from 70 unique genes and 23 alternate splice isoforms representing the largest set of characterized ZFPs from any organism described to date. These recognition motifs can be used to predict genomic binding sites for these factors within the fruit fly genome. Subsets of fingers from these ZFPs were characterized to define their orientation and register on their recognition sequences, thereby allowing us to define the recognition diversity within this finger set. We find that the characterized fingers can specify 47 of the 64 possible DNA triplets. To confirm the utility of our finger recognition models, we employed subsets of Drosophila fingers in combination with an existing archive of artificial zinc finger modules to create ZFPs with novel DNA-binding specificity. These hybrids of natural and artificial fingers can be used to create functional zinc finger nucleases for editing vertebrate genomes.

Genetic landscape of open chromatin in yeast

Chromatin regulation underlies a variety of DNA metabolism processes, including transcription, recombination, repair, and replication. To perform a quantitative genetic analysis of chromatin accessibility, we obtained open chromatin profiles across 96 genetically different yeast strains by FAIRE (formaldehyde-assisted isolation of regulatory elements) assay followed by sequencing. While 5∼10% of open chromatin region (OCRs) were significantly affected by variations in their underlying DNA sequences, subtelomeric areas as well as gene-rich and gene-poor regions displayed high levels of sequence-independent variation. We performed quantitative trait loci (QTL) mapping using the FAIRE signal for each OCR as a quantitative trait. While individual OCRs were associated with a handful of specific genetic markers, gene expression levels were associated with many regulatory loci. We found multi-target trans-loci responsible for a very large number of OCRs, which seemed to reflect the widespread influence of certain chromatin regulators. Such regulatory hotspots were enriched for known regulatory functions, such as recombinational DNA repair, telomere replication, and general transcription control. The OCRs associated with these multi-target trans-loci coincided with recombination hotspots, telomeres, and gene-rich regions according to the function of the associated regulators. Our findings provide a global quantitative picture of the genetic architecture of chromatin regulation.

Evaluation of methods for modeling transcription factor sequence specificity

Genomic analyses often involve scanning for potential transcription factor (TF) binding sites using models of the sequence specificity of DNA binding proteins. Many approaches have been developed to model and learn a protein's DNA-binding specificity, but these methods have not been systematically compared. Here we applied 26 such approaches to in vitro protein binding microarray data for 66 mouse TFs belonging to various families. For nine TFs, we also scored the resulting motif models on in vivo data, and found that the best in vitro-derived motifs performed similarly to motifs derived from the in vivo data. Our results indicate that simple models based on mononucleotide position weight matrices trained by the best methods perform similarly to more complex models for most TFs examined, but fall short in specific cases (<10% of the TFs examined here). In addition, the best-performing motifs typically have relatively low information content, consistent with widespread degeneracy in eukaryotic TF sequence preferences.

A genome-wide cis-regulatory element discovery method based on promoter sequences and gene co-expression networks

Background

Deciphering cis-regulatory networks has become an attractive yet challenging task. This paper presents a simple method for cis-regulatory network discovery which aims to avoid some of the common problems of previous approaches.

Results

Using promoter sequences and gene expression profiles as input, rather than clustering the genes by the expression data, our method utilizes co-expression neighborhood information for each individual gene, thereby overcoming the disadvantages of current clustering based models which may miss specific information for individual genes. In addition, rather than using a motif database as an input, it implements a simple motif count table for each enumerated k-mer for each gene promoter sequence. Thus, it can be used for species where previous knowledge of cis-regulatory motifs is unknown and has the potential to discover new transcription factor binding sites. Applications on Saccharomyces cerevisiae and Arabidopsis have shown that our method has a good prediction accuracy and outperforms a phylogenetic footprinting approach. Furthermore, the top ranked gene-motif regulatory clusters are evidently functionally co-regulated, and the regulatory relationships between the motifs and the enriched biological functions can often be confirmed by literature.

Conclusions

Since this method is simple and gene-specific, it can be readily utilized for insufficiently studied species or flexibly used as an additional step or data source for previous transcription regulatory networks discovery models.

Chromatin dynamics at the replication fork: there's more to life than histones

Before each division, eukaryotic cells face the daunting task of completely and accurately replicating a heterogeneous, chromatinized genome and repackaging both resulting daughters. Because replication requires strand separation, interactions between the DNA and its many associated proteins--including histones--must be transiently broken to allow the passage of the replication fork. Here, we will discuss the disruption and re-establishment of chromatin structure during replication, and the consequences of these processes for epigenetic inheritance.

Functional Implications of Local DNA Structures in Regulatory Motifs

The three-dimensional structure of DNA has been proposed to be a major determinant for functional transcription factors (TFs) and DNA interaction. Here, we use hydroxyl radical cleavage pattern as a measure of local DNA structure. We compared the conservation between DNA sequence and structure in terms of information content and attempted to assess the functional implications of DNA structures in regulatory motifs. We used statistical methods to evaluate the structural divergence of substituting a single position within a binding site and applied them to a collection of putative regulatory motifs. The following are our major observations: (i) we observed more information in structural alignment than in the corresponding sequence alignment for most of the transcriptional factors; (ii) for each TF, majority of positions have more information in the structural alignment as compared to the sequence alignment; (iii) we further defined a DNA structural divergence score (SD score) for each wild-type and mutant pair that is distinguished by single-base mutation. The SD score for benign mutations is significantly lower than that of switch mutations. This indicates structural conservation is also important for TFBS to be functional and DNA structures will provide previously unappreciated information for TF to realize the binding specificity.

Fission yeast hotspot sequence motifs are also active in budding yeast

In most organisms, including humans, meiotic recombination occurs preferentially at a limited number of sites in the genome known as hotspots. There has been substantial progress recently in elucidating the factors determining the location of meiotic recombination hotspots, and it is becoming clear that simple sequence motifs play a significant role. In S. pombe, there are at least five unique sequence motifs that have been shown to produce hotspots of recombination, and it is likely that there are more. In S. cerevisiae, simple sequence motifs have also been shown to produce hotspots or show significant correlations with hotspots. Some of the hotspot motifs in both yeasts are known or suspected to bind transcription factors (TFs), which are required for the activity of those hotspots. Here we show that four of the five hotspot motifs identified in S. pombe also create hotspots in the distantly related budding yeast S. cerevisiae. For one of these hotspots, M26 (also called CRE), we identify TFs, Cst6 and Sko1, that activate and inhibit the hotspot, respectively. In addition, two of the hotspot motifs show significant correlations with naturally occurring hotspots. The conservation of these hotspots between the distantly related fission and budding yeasts suggests that these sequence motifs, and others yet to be discovered, may function widely as hotspots in many diverse organisms.

Genomic insights of protein arginine methyltransferase Hmt1 binding reveals novel regulatory functions

Background

Protein arginine methylation is a post-translational modification involved in important biological processes such as transcription and RNA processing. This modification is catalyzed by both type I and II protein arginine methyltransferases (PRMTs). One of the most conserved type I PRMTs is PRMT1, the homolog of which is Hmt1 in Saccharomyces cerevisiae. Hmt1 has been shown to play a role in various gene expression steps, such as promoting the dynamics of messenger ribonucleoprotein particle (mRNP) biogenesis, pre-mRNA splicing, and silencing of chromatin. To determine the full extent of Hmt1’s involvement during gene expression, we carried out a genome-wide location analysis for Hmt1.

Results

A comprehensive genome-wide binding profile for Hmt1 was obtained by ChIP-chip using NimbleGen high-resolution tiling microarrays. Of the approximately 1000 Hmt1-binding sites found, the majority fall within or proximal to an ORF. Different occupancy patterns of Hmt1 across genes with different transcriptional rates were found. Interestingly, Hmt1 occupancy is found at a number of other genomic features such as tRNA and snoRNA genes, thereby implicating a regulatory role in the biogenesis of these non-coding RNAs. RNA hybridization analysis shows that Hmt1 loss-of-function mutants display higher steady-state tRNA abundance relative to the wild-type. Co-immunoprecipitation studies demonstrate that Hmt1 interacts with the TFIIIB component Bdp1, suggesting a mechanism for Hmt1 in modulating RNA Pol III transcription to regulate tRNA production.

Conclusions

The genome-wide binding profile of Hmt1 reveals multiple potential new roles for Hmt1 in the control of eukaryotic gene expression, especially in the realm of non-coding RNAs. The data obtained here will provide an important blueprint for future mechanistic studies on the described occupancy relationship for genomic features bound by Hmt1.

Tbf1 and Vid22 promote resection and non-homologous end joining of DNA double-strand break ends

The repair of DNA double-strand breaks (DSBs) is crucial for maintaining genome stability. The Saccharomyces cerevisiae protein Tbf1, which is characterized by a Myb domain and is related to mammalian TRF1 and TRF2, has been proposed to act as a transcriptional activator. Here, we show that Tbf1 and its interacting protein Vid22 are new players in the response to DSBs. Inactivation of either TBF1 or VID22 causes hypersensitivity to DSB-inducing agents and shows strong negative interactions with mutations affecting homologous recombination. Furthermore, Tbf1 and Vid22 are recruited to an HO-induced DSB, where they promote both resection of DNA ends and repair by non-homologous end joining. Finally, inactivation of either Tbf1 or Vid22 impairs nucleosome eviction around the DSB, suggesting that these proteins promote efficient repair of the break by influencing chromatin identity in its surroundings.

Natural variation in the yeast glucose-signaling network reveals a new role for the Mig3p transcription factor

The Crabtree effect, in which fermentative metabolism is preferred at the expense of respiration, is a hallmark of budding yeast's glucose response and a model for the Warburg effect in human tumors. While the glucose-responsive transcriptional repressors Mig1p and Mig2p play well-characterized roles in the Crabtree effect, little function for the related Mig3p transcription factor has been uncovered, despite numerous investigations of laboratory yeast strains. Here we studied a wild isolate of Saccharomyces cerevisiae to uncover a critical role for Mig3p that has been lost in S288c-derived laboratory strains. We found that Mig3p affects the expression of hundreds of glucose-responsive genes in the oak strain YPS163, both during growth under standard conditions and upon ethanol treatment. Our results suggest that Mig3p may act as a multifunctional activator/repressor that plays separate roles under standard vs. stress conditions and that this function has been largely lost in the lab strains. Population analysis suggests that the lab strain and several wild strains harbor mutations that diminish Mig3p function. Thus, by expanding our attention to multiple genetic backgrounds, we have uncovered an important missing link in a key metabolic response.

Stabilization of the promoter nucleosomes in nucleosome-free regions by the yeast Cyc8-Tup1 corepressor

The yeast Cyc8 (also known as Ssn6)–Tup1 complex regulates gene expression through a variety of mechanisms, including positioning of nucleosomes over promoters of some target genes to limit accessibility to the transcription machinery. To further define the functions of Cyc8–Tup1 in gene regulation and chromatin remodeling, we performed genome-wide profiling of changes in nucleosome organization and gene expression that occur upon loss of CYC8 or TUP1 and observed extensive nucleosome alterations in both promoters and gene bodies of derepressed genes. Our improved nucleosome profiling and analysis approaches revealed low-occupancy promoter nucleosomes (P nucleosomes) at locations previously defined as nucleosome-free regions. In the absence of CYC8 or TUP1, this P nucleosome is frequently lost, whereas nucleosomes are gained at −1 and +1 positions, accompanying up-regulation of downstream genes. Our analysis of public ChIP-seq data revealed that Cyc8 and Tup1 preferentially bind TATA-containing promoters, which are also enriched in genes derepressed upon loss of CYC8 or TUP1. These results suggest that stabilization of the P nucleosome on TATA-containing promoters may be a central feature of the repressive chromatin architecture created by the Cyc8–Tup1 corepressor, and that releasing the P nucleosome contributes to gene activation.

Timing of transcriptional quiescence during gametogenesis is controlled by global histone H3K4 demethylation

Gametes are among the most highly specialized cells produced during development. Although gametogenesis culminates in transcriptional quiescence in plants and animals, regulatory mechanisms controlling this are unknown. Here, we confirm that gamete differentiation in the single-celled yeast Saccharomyces cerevisiae is accompanied by global transcriptional shutoff following the completion of meiosis. We show that Jhd2, a highly conserved JARID1-family histone H3K4 demethylase, activates protein-coding gene transcription in opposition to this programmed transcriptional shutoff, sustaining the period of productive transcription during spore differentiation. Moreover, using genome-wide nucleosome, H3K4me, and transcript mapping experiments, we demonstrate that JHD2 globally represses intergenic noncoding transcription during this period. The widespread transcriptional defects of JHD2 mutants are associated with precocious differentiation and the production of stress-sensitive spores, demonstrating that Jhd2 regulation of the global postmeiotic transcriptional program is critical for the production of healthy meiotic progeny.

CHD1 remodelers regulate nucleosome spacing in vitro and align nucleosomal arrays over gene coding regions in S. pombe

Nucleosome positioning governs access to eukaryotic genomes. Many genes show a stereotypic organisation at their 5'end: a nucleosome free region just upstream of the transcription start site (TSS) followed by a regular nucleosomal array over the coding region. The determinants for this pattern are unclear, but nucleosome remodelers are likely critical. Here we study the role of remodelers in global nucleosome positioning in S. pombe and the corresponding changes in expression. We find a striking evolutionary shift in remodeler usage between budding and fission yeast. The S. pombe RSC complex does not seem to be involved in nucleosome positioning, despite its prominent role in S. cerevisiae. While S. pombe lacks ISWI-type remodelers, it has two CHD1-type ATPases, Hrp1 and Hrp3. We demonstrate nucleosome spacing activity for Hrp1 and Hrp3 in vitro, and that together they are essential for linking regular genic arrays to most TSSs in vivo. Impaired arrays in the absence of either or both remodelers may lead to increased cryptic antisense transcription, but overall gene expression levels are only mildly affected.

Deciphering the transcriptional cis-regulatory code

Information about developmental gene expression resides in defined regulatory elements, called enhancers, in the non-coding part of the genome. Although cells reliably utilize enhancers to orchestrate gene expression, a cis-regulatory code that would allow their interpretation has remained one of the greatest challenges of modern biology. In this review, we summarize studies from the past three decades that describe progress towards revealing the properties of enhancers and discuss how recent approaches are providing unprecedented insights into regulatory elements in animal genomes. Over the next years, we believe that the functional characterization of regulatory sequences in entire genomes, combined with recent computational methods, will provide a comprehensive view of genomic regulatory elements and their building blocks and will enable researchers to begin to understand the sequence basis of the cis-regulatory code.

DNA sequence preferences of transcriptional activators correlate more strongly than repressors with nucleosomes

Transcription factors (TFs) and histone octamers are two abundant classes of DNA binding proteins that coordinate the transcriptional program in cells. Detailed studies of individual TFs have shown that TFs bind to nucleosome-occluded DNA sequences and induce nucleosome disruption/repositioning, while recent global studies suggest this is not the only mechanism used by all TFs. We have analyzed to what extent the intrinsic DNA binding preferences of TFs and histones play a role in determining nucleosome occupancy, in addition to nonintrinsic factors such as the enzymatic activity of chromatin remodelers. The majority of TFs in budding yeast have an intrinsic sequence preference overlapping with nucleosomal histones. TFs with intrinsic DNA binding properties highly correlated with those of histones tend to be associated with gene activation and might compete with histones to bind to genomic DNA. Consistent with this, we show that activators induce more nucleosome disruption upon transcriptional activation than repressors.

Basic leucine zipper transcription factor Hac1 binds DNA in two distinct modes as revealed by microfluidic analyses

A quantitative understanding of how transcription factors interact with genomic target sites is crucial for reconstructing transcriptional networks in vivo. Here, we use Hac1, a well-characterized basic leucine zipper (bZIP) transcription factor involved in the unfolded protein response (UPR) as a model to investigate interactions between bZIP transcription factors and their target sites. During the UPR, the accumulation of unfolded proteins leads to unconventional splicing and subsequent translation of HAC1 mRNA, followed by transcription of UPR target genes. Initial candidate-based approaches identified a canonical cis-acting unfolded protein response element (UPRE-1) within target gene promoters; however, subsequent studies identified a large set of Hac1 target genes lacking this UPRE-1 and containing a different motif (UPRE-2). Using a combination of unbiased and directed microfluidic DNA binding assays, we established that Hac1 binds in two distinct modes: (i) to short (6-7 bp) UPRE-2-like motifs and (ii) to significantly longer (11-13 bp) extended UPRE-1-like motifs. Using a genetic screen, we demonstrate that a region of extended homology N-terminal to the basic DNA binding domain is required for this dual site recognition. These results establish Hac1 as the first bZIP transcription factor known to adopt more than one binding mode and unify previously conflicting and discrepant observations of Hac1 function into a cohesive model of UPR target gene activation. Our results also suggest that even structurally simple transcription factors can recognize multiple divergent target sites of very different lengths, potentially enriching their downstream target repertoire.

Massively parallel measurements of molecular interaction kinetics on a microfluidic platform

Quantitative biology requires quantitative data. No high-throughput technologies exist capable of obtaining several hundred independent kinetic binding measurements in a single experiment. We present an integrated microfluidic device (k-MITOMI) for the simultaneous kinetic characterization of 768 biomolecular interactions. We applied k-MITOMI to the kinetic analysis of transcription factor (TF)-DNA interactions, measuring the detailed kinetic landscapes of the mouse TF Zif268, and the yeast TFs Tye7p, Yox1p, and Tbf1p. We demonstrated the integrated nature of k-MITOMI by expressing, purifying, and characterizing 27 additional yeast transcription factors in parallel on a single device. Overall, we obtained 2,388 association and dissociation curves of 223 unique molecular interactions with equilibrium dissociation constants ranging from 2 × 10(-6) M to 2 × 10(-9) M, and dissociation rate constants of approximately 6 s(-1) to 8.5 × 10(-3) s(-1). Association rate constants were uniform across 3 TF families, ranging from 3.7 × 10(6) M(-1) s(-1) to 9.6 × 10(7) M(-1) s(-1), and are well below the diffusion limit. We expect that k-MITOMI will contribute to our quantitative understanding of biological systems and accelerate the development and characterization of engineered systems.

P-value-based regulatory motif discovery using positional weight matrices

To analyze gene regulatory networks, the sequence-dependent DNA/RNA binding affinities of proteins and noncoding RNAs are crucial. Often, these are deduced from sets of sequences enriched in factor binding sites. Two classes of computational approaches exist. The first describe binding motifs by sequence patterns and search the patterns with highest statistical significance for enrichment. The second class uses the more powerful position weight matrices (PWMs). Instead of maximizing the statistical significance of enrichment, they maximize a likelihood. Here we present XXmotif (eXhaustive evaluation of matriX motifs), the first PWM-based motif discovery method that can optimize PWMs by directly minimizing their P-values of enrichment. Optimization requires computing millions of enrichment P-values for thousands of PWMs. For a given PWM, the enrichment P-value is calculated efficiently from the match P-values of all possible motif placements in the input sequences using order statistics. The approach can naturally combine P-values for motif enrichment, conservation, and localization. On ChIP-chip/seq, miRNA knock-down, and coexpression data sets from yeast and metazoans, XXmotif outperformed state-of-the-art tools, both in numbers of correctly identified motifs and in the quality of PWMs. In segmentation modules of D. melanogaster, we detect the known key regulators and several new motifs. In human core promoters, XXmotif reports most previously described and eight novel motifs sharply peaked around the transcription start site, among them an Initiator motif similar to the fly and yeast versions. XXmotif's sensitivity, reliability, and usability will help to leverage the quickly accumulating wealth of functional genomics data.

Differences in local genomic context of bound and unbound motifs

Understanding gene regulation is a major objective in molecular biology research. Frequently, transcription is driven by transcription factors (TFs) that bind to specific DNA sequences. These motifs are usually short and degenerate, rendering the likelihood of multiple copies occurring throughout the genome due to random chance as high. Despite this, TFs only bind to a small subset of sites, thus prompting our investigation into the differences between motifs that are bound by TFs and those that remain unbound. Here we constructed vectors representing various chromatin- and sequence-based features for a published set of bound and unbound motifs representing nine TFs in the budding yeast Saccharomyces cerevisiae. Using a machine learning approach, we identified a set of features that can be used to discriminate between bound and unbound motifs. We also discovered that some TFs bind most or all of their strong motifs in intergenic regions. Our data demonstrate that local sequence context can be strikingly different around motifs that are bound compared to motifs that are unbound. We concluded that there are multiple combinations of genomic features that characterize bound or unbound motifs.

Genome-wide nucleosome specificity and directionality of chromatin remodelers

How chromatin remodelers cooperate to organize nucleosomes around the start and end of genes is not known. We determined the genome-wide binding of remodeler complexes SWI/SNF, RSC, ISW1a, ISW1b, ISW2, and INO80 to individual nucleosomes in Saccharomyces, and determined their functional contributions to nucleosome positioning through deletion analysis. We applied ultra-high-resolution ChIP-exo mapping to Isw2 to determine its subnucleosomal orientation and organization on a genomic scale. Remodelers interacted with selected nucleosome positions relative to the start and end of genes and produced net directionality in moving nucleosomes either away or toward nucleosome-free regions at the 5' and 3' ends of genes. Isw2 possessed a subnucleosomal organization in accord with biochemical and crystallographic-based models that place its linker binding region within promoters and abutted against Reb1-bound locations. Together, these findings reveal a coordinated position-specific approach taken by remodelers to organize genic nucleosomes into arrays.

Nutritional control of growth and development in yeast

Availability of key nutrients, such as sugars, amino acids, and nitrogen compounds, dictates the developmental programs and the growth rates of yeast cells. A number of overlapping signaling networks--those centered on Ras/protein kinase A, AMP-activated kinase, and target of rapamycin complex I, for instance--inform cells on nutrient availability and influence the cells' transcriptional, translational, posttranslational, and metabolic profiles as well as their developmental decisions. Here I review our current understanding of the structures of the networks responsible for assessing the quantity and quality of carbon and nitrogen sources. I review how these signaling pathways impinge on transcriptional, metabolic, and developmental programs to optimize survival of cells under different environmental conditions. I highlight the profound knowledge we have gained on the structure of these signaling networks but also emphasize the limits of our current understanding of the dynamics of these signaling networks. Moreover, the conservation of these pathways has allowed us to extrapolate our finding with yeast to address issues of lifespan, cancer metabolism, and growth control in more complex organisms.