Active PHO5 chromatin encompasses variable numbers of nucleosomes at individual promoters

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here, we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Antisense non-coding transcription represses the PHO5 model gene at the level of promoter chromatin structure

Pervasive transcription of eukaryotic genomes generates non-coding transcripts with regulatory potential. We examined the effects of non-coding antisense transcription on the regulation of expression of the yeast PHO5 gene, a paradigmatic case for gene regulation through promoter chromatin remodeling. A negative role for antisense transcription at the PHO5 gene locus was demonstrated by leveraging the level of overlapping antisense transcription through specific mutant backgrounds, expression from a strong promoter in cis, and use of the CRISPRi system. Furthermore, we showed that enhanced elongation of PHO5 antisense leads to a more repressive chromatin conformation at the PHO5 gene promoter, which is more slowly remodeled upon gene induction. The negative effect of antisense transcription on PHO5 gene transcription is mitigated upon inactivation of the histone deacetylase Rpd3, showing that PHO5 antisense RNA acts via histone deacetylation. This regulatory pathway leads to Rpd3-dependent decreased recruitment of the RSC chromatin remodeling complex to the PHO5 gene promoter upon induction of antisense transcription. Overall, the data in this work reveal an additional level in the complex regulatory mechanism of PHO5 gene expression by showing antisense transcription-mediated repression at the level of promoter chromatin structure remodeling.

Partitioned usage of chromatin remodelers by nucleosome-displacing factors

Nucleosome-displacing-factors (NDFs) in yeast, similar to pioneer factors in higher eukaryotes, can open closed chromatin and generate nucleosome-depleted regions (NDRs). NDRs in yeast are also affected by ATP-dependent chromatin remodelers (CRs). However, how NDFs and CRs coordinate in nucleosome invasion and NDR formation is still unclear. Here, we design a high-throughput method to systematically study the interplay between NDFs and CRs. By combining an integrated synthetic oligonucleotide library with DNA methyltransferase-based, single-molecule nucleosome mapping, we measure the impact of CRs on NDRs generated by individual NDFs. We find that CRs are dispensable for nucleosome invasion by NDFs, and they function downstream of NDF binding to modulate the NDR length. A few CRs show high specificity toward certain NDFs; however, in most cases, CRs are recruited in a factor-nonspecific and NDR length-dependent manner. Overall, our study provides a framework to investigate how NDFs and CRs cooperate to regulate chromatin opening.

Chromatin accessibility in yeast is regulated by nucleosome-displacing-factors (NDFs) and chromatin remodelers (CRs). Chen et al. show that NDFs first invade into nucleosomes and then recruit CRs to modulate the NDR length. NDF-specific and NDR length-dependent recruitment of CRs allow partitioned usage of CRs by NDFs.

Nucleosomal proofreading of activator-promoter interactions

Specificity in transcriptional regulation is imparted by transcriptional activators that bind to specific DNA sequences from which they stimulate transcription. Specificity may be increased by slowing down the kinetics of regulation: by increasing the energy for dissociation of the activator-DNA complex or decreasing activator concentration. In general, higher dissociation energies imply longer DNA dwell times of the activator; the activator-bound gene may not readily turn off again. Lower activator concentrations entail longer pauses between binding events; the activator-unbound gene is not easily turned on again and activated transcription occurs in stochastic bursts. We show that kinetic proofreading of activator-DNA recognition-insertion of an energy-dissipating delay step into the activation pathway for transcription-reconciles high specificity of transcriptional regulation with fast regulatory kinetics. We show that kinetic proofreading results from the stochastic removal and reformation of promoter nucleosomes, at a distance from equilibrium.

Nucleosome positioning and spacing: from genome-wide maps to single arrays

The positioning of nucleosomes relative to DNA and their neighboring nucleosomes represents a fundamental layer of chromatin organization. Changes in nucleosome positioning and spacing affect the accessibility of DNA to regulatory factors and the formation of higher order chromatin structures. Sequencing of mononucleosomal fragments allowed mapping nucleosome positions on a genome-wide level in many organisms. This revealed that successions of evenly spaced and well-positioned nucleosomes-so called phased nucleosome arrays-occur at the 5' end of many active genes and in the vicinity of transcription factor and other protein binding sites. Phased arrays arise from the interplay of barrier elements on the DNA, which position adjacent nucleosomes, and the nucleosome spacing activity of ATP-dependent chromatin remodelers. A shortcoming of classic mononucleosomal mapping experiments is that they only reveal nucleosome spacing and array regularity at select sites in the genome with well-positioned nucleosomes. However, new technological approaches elucidate nucleosome array structure throughout the genome and with single-cell resolution. In the future, it will be interesting to see whether changes in nucleosome array regularity and spacing contribute to the formation of higher order chromatin structures and the spatial organization of the genome in vivo.

Chromatin structure analysis of single gene molecules by psoralen cross-linking and electron microscopy

Nucleosomes occupy a central role in regulating eukaryotic gene expression by blocking access of transcription factors to their target sites on chromosomal DNA. Analysis of chromatin structure and function has mostly been performed by probing DNA accessibility with endonucleases. Such experiments average over large numbers of molecules of the same gene, and more recently, over entire genomes. However, both digestion and averaging erase the structural variation between molecules indicative of dynamic behavior, which must be reconstructed for any theory of regulation. Solution of this problem requires the structural analysis of single gene molecules. In this chapter, we describe a method by which single gene molecules are purified from the yeast Saccharomyces cerevisiae and cross-linked with psoralen, allowing the determination of nucleosome configurations by transmission electron microscopy. We also provide custom analysis software that semi-automates the analysis of micrograph data. This single-gene technique enables detailed examination of chromatin structure at any genomic locus in yeast.

Integrated DNA methylation and chromatin structural analysis at single-molecule resolution

Chromatin limits the accessibility of DNA to trans-acting factors in transcription, replication, and repair. Although transcriptional variation between cells in a population may contribute to survival and disease, most assays of chromatin structure recover only population averages. We have developed DNA methyltransferases (MTases) as probing agents of DNA accessibility in chromatin, either expressed in vivo in budding yeast or as recombinant enzymatic probes of nuclei isolated from mammalian cells. In this chapter, we focus on the use of recombinant MTase (M) M.CviPI to probe chromatin accessibility in nuclei isolated from mammalian cell lines and animal tissue. This technique, named methylation accessibility protocol for individual templates (MAPit), reports protein-DNA interactions at single-molecule resolution. The single-molecule readout allows identification of chromatin subpopulations and rare epigenetic variants within a cell population. Furthermore, the use of M.CviPI in mammalian systems gives a comprehensive view of both chromatin structure and endogenous DNA methylation in a single assay.

Nucleosomal promoter variation generates gene expression noise

Gene product molecule numbers fluctuate over time and between cells, confounding deterministic expectations. The molecular origins of this noise of gene expression remain unknown. Recent EM analysis of single PHO5 gene molecules of yeast indicated that promoter molecules stochastically assume alternative nucleosome configurations at steady state, including the fully nucleosomal and nucleosome-free configuration. Given that distinct configurations are unequally conducive to transcription, the nucleosomal variation of promoter molecules may constitute a source of gene expression noise. This notion, however, implies an untested conjecture, namely that the nucleosomal variation arises de novo or intrinsically (i.e., that it cannot be explained as the result of the promoter's deterministic response to variation in its molecular surroundings). Here, we show--by microscopically analyzing the nucleosome configurations of two juxtaposed physically linked PHO5 promoter copies--that the configurational variation, indeed, is intrinsically stochastic and thus, a cause of gene expression noise rather than its effect.

Epigenetics reloaded: the single-cell revolution

Mechanistically, how epigenetic states are inherited through cellular divisions remains an important open question in the chromatin field and beyond. Defining the heritability of epigenetic states and the underlying chromatin-based mechanisms within a population of cells is complicated due to cell heterogeneity combined with varying levels of stability of these states; thus, efforts must be focused toward single-cell analyses. The approaches presented here constitute the forefront of epigenetics research at the single-cell level using classic and innovative methods to dissect epigenetics mechanisms from the limited material available in a single cell. This review further outlines exciting future avenues of research to address the significance of epigenetic heterogeneity and the contributions of microfluidics technologies to single-cell isolation and analysis.

The yeast PHO5 promoter: from single locus to systems biology of a paradigm for gene regulation through chromatin

Chromatin dynamics crucially contributes to gene regulation. Studies of the yeast PHO5 promoter were key to establish this nowadays accepted view and continuously provide mechanistic insight in chromatin remodeling and promoter regulation, both on single locus as well as on systems level. The PHO5 promoter is a context independent chromatin switch module where in the repressed state positioned nucleosomes occlude transcription factor sites such that nucleosome remodeling is a prerequisite for and not consequence of induced gene transcription. This massive chromatin transition from positioned nucleosomes to an extensive hypersensitive site, together with respective transitions at the co-regulated PHO8 and PHO84 promoters, became a prime model for dissecting how remodelers, histone modifiers and chaperones co-operate in nucleosome remodeling upon gene induction. This revealed a surprisingly complex cofactor network at the PHO5 promoter, including five remodeler ATPases (SWI/SNF, RSC, INO80, Isw1, Chd1), and demonstrated for the first time histone eviction in trans as remodeling mode in vivo. Recently, the PHO5 promoter and the whole PHO regulon were harnessed for quantitative analyses and computational modeling of remodeling, transcription factor binding and promoter input-output relations such that this rewarding single-locus model becomes a paradigm also for theoretical and systems approaches to gene regulatory networks.

Reconfiguration of nucleosome-depleted regions at distal regulatory elements accompanies DNA methylation of enhancers and insulators in cancer

It is well established that cancer-associated epigenetic repression occurs concomitant with CpG island hypermethylation and loss of nucleosomes at promoters, but the role of nucleosome occupancy and epigenetic reprogramming at distal regulatory elements in cancer is still poorly understood. Here, we evaluate the scope of global epigenetic alterations at enhancers and insulator elements in prostate and breast cancer cells using simultaneous genome-wide mapping of DNA methylation and nucleosome occupancy (NOMe-seq). We find that the genomic location of nucleosome-depleted regions (NDRs) is mostly cell type specific and preferentially found at enhancers in normal cells. In cancer cells, however, we observe a global reconfiguration of NDRs at distal regulatory elements coupled with a substantial reorganization of the cancer methylome. Aberrant acquisition of nucleosomes at enhancer-associated NDRs is associated with hypermethylation and epigenetic silencing marks, and conversely, loss of nucleosomes with demethylation and epigenetic activation. Remarkably, we show that nucleosomes remain strongly organized and phased at many facultative distal regulatory elements, even in the absence of a NDR as an anchor. Finally, we find that key transcription factor (TF) binding sites also show extensive peripheral nucleosome phasing, suggesting the potential for TFs to organize NDRs genome-wide and contribute to deregulation of cancer epigenomes. Together, our findings suggest that “decommissioning” of NDRs and TFs at distal regulatory elements in cancer cells is accompanied by DNA hypermethylation susceptibility of enhancers and insulator elements, which in turn may contribute to an altered genome-wide architecture and epigenetic deregulation in malignancy.

Single-cell nucleosome mapping reveals the molecular basis of gene expression heterogeneity

Nucleosomes, the basic unit of chromatin, have a critical role in the control of gene expression. Nucleosome positions have generally been determined by examining bulk populations of cells and then correlated with overall gene expression. Here, we describe a technique to determine nucleosome positioning in single cells by virtue of the ability of the nucleosome to protect DNA from GpC methylation. In the acid phosphatase inducible PHO5 gene, we find that there is significant cell-to-cell variation in nucleosome positions and shifts in nucleosome positioning correlate with changes in gene expression. However, nucleosome positioning is not absolute, and even with major shifts in gene expression, some cells fail to change nucleosome configuration. Mutations of the PHO5 promoter that introduce a poly(dA:dT) tract-stimulated gene expression under nonpermissive conditions led to shifts of positioned nucleosomes similar to induction of PHO5. By contrast, mutations that altered AA/TT/AT periodicity reduced gene expression upon PHO5 induction and stabilized nucleosomes in most cells, suggesting that enhanced nucleosome affinity for DNA antagonizes chromatin remodelers. Finally, we determined nucleosome positioning in two regions described as "fuzzy" or nucleosome-free when examined in a bulk assay. These regions consisted of distinct nucleosomes with a larger footprint for potential location and an increase population of cells lacking a nucleosome altogether. These data indicate an underlying complexity of nucleosome positioning that may contribute to the flexibility and heterogeneity of gene expression.

Linking stochastic fluctuations in chromatin structure and gene expression

Electron microscopy of single gene molecules and mathematical modeling shows that a promoter stochastically transitions between transcriptionally favorable and unfavorable nucleosome configurations, providing a mechanism for transcriptional bursting.

The number of mRNA and protein molecules expressed from a single gene molecule fluctuates over time. These fluctuations have been attributed, in part, to the random transitioning of promoters between transcriptionally active and inactive states, causing transcription to occur in bursts. However, the molecular basis of transcriptional bursting remains poorly understood. By electron microscopy of single PHO5 gene molecules from yeast, we show that the “activated” promoter assumes alternative nucleosome configurations at steady state, including the maximally repressive, fully nucleosomal, and the maximally non-repressive, nucleosome-free, configuration. We demonstrate that the observed probabilities of promoter nucleosome configurations are obtained from a simple, intrinsically stochastic process of nucleosome assembly, disassembly, and position-specific sliding; and we show that gene expression and promoter nucleosome configuration can be mechanistically coupled, relating promoter nucleosome dynamics and gene expression fluctuations. Together, our findings suggest a structural basis for transcriptional bursting, and offer new insights into the mechanism of transcriptional regulation and the kinetics of promoter nucleosome transitions.

In eukaryotes, such as plants, fungi, and animals, the DNA is wrapped around basic protein cores called nucleosomes at more or less regular intervals. This wrapping discourages transcription, the first step in gene expression. By isolating PHO5 gene molecules from yeast cells and analyzing their structure by electron microscopy, we provide evidence that the “nucleosomes” completely unwrap and then re-wrap in an intrinsically stochastic manner. Only nucleosomes that wrap the regulatory sequences of the gene (promoter) were observed to unspool; no such unspooling was found across the body of the gene. Random unwrapping and re-wrapping generates an ensemble of alternative promoter nucleosome configurations, some conducive to transcription, others not. Mounting evidence suggests that transcription occurs in bursts, where transcripts are released in close succession, interrupted by intervals of transcriptional inactivity; this may lead to significant stochastic fluctuations in gene expression. Although the mechanism of this behavior is not understood, our findings now provide a structural basis for it, suggesting that spooling and unspooling of promoter DNA from the nucleosomes determines the fundamental frequency of transcriptional bursting.

The effect of an intervening promoter nucleosome on gene expression

Nucleosomes, which are the basic packaging units of chromatin, are stably positioned in promoters upstream of most stress-inducible genes. These promoter nucleosomes are generally thought to repress gene expression due to exclusion; they prevent transcription factors from accessing their target sites on the DNA. However, the role of promoter nucleosomes that do not directly occlude transcription factor binding sites is not obvious. Here, we varied the stability of a non-occluding nucleosome positioned between a transcription factor binding site and the TATA box region in an inducible yeast promoter and measured downstream gene expression level. We found that gene expression level depends on the occupancy of the non-occluding nucleosome in a non-monotonic manner. We postulated that a non-occluding nucleosome can serve both as a vehicle of and a barrier to chromatin remodeling activity and built a quantitative, nonequilibrium model to explain the observed nontrivial effect of the intervening nucleosome. Our work sheds light on the dual role of nucleosome as a repressor and an activator and expands the standard model of gene expression to include irreversible promoter chromatin transitions.

Epigenetic diversity of Kaposi's sarcoma-associated herpesvirus

Spontaneous lytic reactivation of Kaposi’s sarcoma–associated herpesvirus (KSHV) occurs at a low rate in latently infected cells in disease and culture. This suggests imperfect epigenetic maintenance of viral transcription programs, perhaps due to variability in chromatin structure at specific loci across the population of KSHV episomal genomes. To characterize this locus-specific chromatin structural diversity, we used MAPit single-molecule footprinting, which simultaneously maps endogenous CG methylation and accessibility to M.CviPI at GC sites. Diverse chromatin structures were detected at the LANA, RTA and vIL6 promoters. At each locus, chromatin ranged from fully closed to fully open across the population. This diversity has not previously been reported in a virus. Phorbol ester and RTA transgene induction were used to identify chromatin conformations associated with reactivation of lytic transcription, which only a fraction of episomes had. Moreover, certain chromatin conformations correlated with CG methylation patterns at the RTA and vIL6 promoters. This indicated that some of the diverse chromatin conformations at these loci were epigenetically distinct. Finally, by comparing chromatin structures from a cell line infected with constitutively latent virus, we identified products of lytic replication. Our findings show that epigenetic drift can restrict viral propagation by chromatin compaction at latent and lytic promoters.

Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules

DNA methylation and nucleosome positioning work together to generate chromatin structures that regulate gene expression. Nucleosomes are typically mapped using nuclease digestion requiring significant amounts of material and varying enzyme concentrations. We have developed a method (NOMe-seq) that uses a GpC methyltransferase (M.CviPI) and next generation sequencing to generate a high resolution footprint of nucleosome positioning genome-wide using less than 1 million cells while retaining endogenous DNA methylation information from the same DNA strand. Using a novel bioinformatics pipeline, we show a striking anti-correlation between nucleosome occupancy and DNA methylation at CTCF regions that is not present at promoters. We further show that the extent of nucleosome depletion at promoters is directly correlated to expression level and can accommodate multiple nucleosomes and provide genome-wide evidence that expressed non-CpG island promoters are nucleosome-depleted. Importantly, NOMe-seq obtains DNA methylation and nucleosome positioning information from the same DNA molecule, giving the first genome-wide DNA methylation and nucleosome positioning correlation at the single molecule, and thus, single cell level, that can be used to monitor disease progression and response to therapy.

Chromatin and transcription in yeast

Understanding the mechanisms by which chromatin structure controls eukaryotic transcription has been an intense area of investigation for the past 25 years. Many of the key discoveries that created the foundation for this field came from studies of Saccharomyces cerevisiae, including the discovery of the role of chromatin in transcriptional silencing, as well as the discovery of chromatin-remodeling factors and histone modification activities. Since that time, studies in yeast have continued to contribute in leading ways. This review article summarizes the large body of yeast studies in this field.

Simultaneous single-molecule detection of endogenous C-5 DNA methylation and chromatin accessibility using MAPit

Bisulfite genomic sequencing provides a single-molecule view of cytosine methylation states. After deamination, each cloned molecule contains a record of methylation within its sequence. The full power of this technique is harnessed by treating nuclei with an exogenous DNMT prior to DNA extraction. This exogenous methylation marks regions of accessibility and footprints nucleosomes, as well as other DNA-binding proteins. Thus, each cloned molecule records not only the endogenous methylation present (at CG sites, in mammals), but also the exogenous (GC, when using the Chlorella virus protein M.CviPI). We term this technique MAPit, methylation accessibility protocol for individual templates.

DNA methyltransferase accessibility protocol for individual templates by deep sequencing

A single-molecule probe of chromatin structure can uncover dynamic chromatin states and rare epigenetic variants of biological importance that bulk measures of chromatin structure miss. In bisulfite genomic sequencing, each sequenced clone records the methylation status of multiple sites on an individual molecule of DNA. An exogenous DNA methyltransferase can thus be used to image nucleosomes and other protein-DNA complexes. In this chapter, we describe the adaptation of this technique, termed Methylation Accessibility Protocol for individual templates, to modern high-throughput sequencing, which both simplifies the workflow and extends its utility.

Signal-dependent dynamics of transcription factor translocation controls gene expression

Information about environmental stimuli is often transmitted using common signaling molecules, but the mechanisms that ensure signaling specificity are not entirely known. Here we show that the identities and intensities of different stresses are transmitted by modulation of the amplitude, duration or frequency of nuclear translocation of the Saccharomyces cerevisiae general stress response transcription factor Msn2. Through artificial control of the dynamics of Msn2 translocation, we reveal how distinct dynamical schemes differentially affect reporter gene expression. Using a simple model, we predict stress-induced reporter gene expression from single-cell translocation dynamics. We then demonstrate that the response of natural target genes to dynamical modulation of Msn2 translocation is influenced by differences in the kinetics of promoter transitions and transcription factor binding properties. Thus, multiple environmental signals can trigger qualitatively different dynamics of a single transcription factor and influence gene expression patterns.

WIF1 is a frequent target for epigenetic silencing in squamous cell carcinoma of the cervix

Aberrant activation of the Wnt/β-catenin signaling axis is a prominent oncogenic mechanism in numerous cancers including cervical cancer. Wnt inhibitory factor-1 (WIF1) is a secreted protein that binds Wnt and antagonizes Wnt activity. While the WIF1 gene is characterized as a target for epigenetic silencing in some tumor types, WIF1 expression has not been examined in human cervical tissue and cervical cancer. Here, we show that WIF1 is unmethylated and its gene product is expressed in normal cervical epithelium and some cultured cervical tumor lines. In contrast, several cervical cancer lines contained dense CpG methylation within the WIF1 gene, and expression of both WIF1 transcript and protein was restored by culturing cells in the presence of the global DNA demethylating agent 5-aza-2'-deoxycytidine. Using single-molecule MAPit methylation footprinting, we observed differences in chromatin structure within the WIF1 promoter region between cell lines that express and those that do not express WIF1, consistent with transcriptional activity and repression, respectively. The WIF1 promoter was aberrantly methylated in ∼60% (10 of 17) high-grade highly undifferentiated squamous cell cervical tumors examined, whereas paired normal tissue showed significantly lower levels of CpG methylation. WIF1 protein was not detectable by immunohistochemistry in tumors with quantitatively high levels of WIF1 methylation. Of note, WIF1 protein was not detectable in two of the seven unmethylated cervical tumors examined, suggesting other mechanisms may contribute WIF1 repression. Our findings establish the WIF1 gene as a frequent target for epigenetic silencing in squamous cell carcinoma of the cervix.

Simultaneous single-molecule mapping of protein-DNA interactions and DNA methylation by MAPit

Sites of protein binding to DNA are inferred from footprints or spans of protection against a probing reagent. In most protocols, sites of accessibility to a probe are detected by mapping breaks in DNA strands. As discussed in this unit, such methods obscure molecular heterogeneity by averaging cuts at a given site over all DNA strands in a sample population. The DNA methyltransferase accessibility protocol for individual templates (MAPit), an alternative method described in this unit, localizes protein-DNA interactions by probing with cytosine-modifying DNA methyltransferases followed by bisulfite sequencing. Sequencing individual DNA products after amplification of bisulfite-converted sequences permits assignment of the methylation status of every enzyme target site along a single DNA strand. Use of the GC-methylating enzyme M.CviPI allows simultaneous mapping of chromatin accessibility and endogenous CpG methylation. MAPit is therefore the only footprinting method that can detect subpopulations of molecules with distinct patterns of protein binding or chromatin architecture and correlate them directly with the occurrence of endogenous methylation. Additional advantages of MAPit methylation footprinting as well as considerations for experimental design and potential sources of error are discussed.

Determinants and dynamics of genome accessibility

In eukaryotes, all DNA-templated reactions occur in the context of chromatin. Nucleosome packaging inherently restricts DNA accessibility for regulatory proteins but also provides an opportunity to regulate DNA-based processes through modulating nucleosome positions and local chromatin structure. Recent advances in genome-scale methods are yielding increasingly detailed profiles of the genomic distribution of nucleosomes, their modifications and their modifiers. The picture now emerging is one in which the dynamic control of genome accessibility is governed by contributions from DNA sequence, ATP-dependent chromatin remodelling and nucleosome modifications. Here we discuss the interplay of these processes by reviewing our current understanding of how chromatin access contributes to the regulation of transcription, replication and repair.

Quantitative analysis of the transcription control mechanism

Activated PHO5 promoter chromatin at steady state represents a statistical ensemble of distinct structures.

The extent of promoter nucleosome loss depends on the strength of the transcriptional activator of PHO5, indicative of continuous disassembly and reassembly of nucleosomes at the induced promoter.

PHO5 promoter nucleosome loss and expression are exponentially related, pointing at two or more steps of the expression process that are activator controlled.

The intrinsic noise profile of PHO5 expression permits quantitative distinction between alternative regulatory architectures.

The assumption of two activator-controlled steps, promoter nucleosome removal and assembly of the transcription machinery, is necessary and sufficient to account for the quantitative relationship between PHO5 expression, intrinsic noise, and promoter nucleosome loss.

The search for factors that interact with transcriptional activators has provided clues to the possible mechanisms of transcriptional regulation. However, the promiscuity of activator interactions has supported alternative regulatory hypotheses whose relative likelihood could not be assessed. A possible solution to the problem is suggested by recent theoretical work showing that the regulatory architecture of the gene expression process is reflected in the magnitude and frequency spectrum of steady-state fluctuations in gene expression (intrinsic noise). In this study, we analyzed promoter chromatin structure and the magnitude of intrinsic noise over a wide range of expression values for the PHO5 gene of yeast. We show that the relationship between gene expression and promoter nucleosome loss is exponential rather than linear, suggesting activator control of at least two steps of the expression process. Interpretation of our data with regard to a stochastic model of chromatin remodeling and gene expression allowed us to quantitatively distinguish between competing regulatory hypotheses, demonstrating the value of noise measurements for analysis of the gene regulatory mechanism. We show that the assumption of two activator-controlled steps, nucleosome removal and assembly of the transcription machinery, is necessary and sufficient to account for the quantitative relationship between our experimental observables.

Gene transcription requires a sequence of promoter state transitions, including chromatin remodeling, assembly of the transcription machinery, and clearance of the promoter by RNA polymerase. The rate-limiting steps in this sequence are regulated by transcriptional activators that bind at specific promoter elements. As the transition kinetics of individual promoters cannot be observed, the identity of the activator-controlled steps has remained a matter of speculation. In this study, we investigated promoter chromatin structure, and the intrinsic noise of expression over a wide range of expression values for the PHO5 gene of yeast. Interpretation of our results with regard to a stochastic model of promoter chromatin remodeling and gene expression suggests that the regulatory architecture of the gene expression process is measurably reflected in its intrinsic noise profile. Our chromatin structure and noise analyses indicate that the activator of PHO5 transcription stimulates the rates of promoter nucleosome disassembly, and assembly of the transcription machinery after nucleosome removal, but no other rates of the expression process.

MethylViewer: computational analysis and editing for bisulfite sequencing and methyltransferase accessibility protocol for individual templates (MAPit) projects

Bisulfite sequencing is a widely-used technique for examining cytosine DNA methylation at nucleotide resolution along single DNA strands. Probing with cytosine DNA methyltransferases followed by bisulfite sequencing (MAPit) is an effective technique for mapping protein-DNA interactions. Here, MAPit methylation footprinting with M.CviPI, a GC methyltransferase we previously cloned and characterized, was used to probe hMLH1 chromatin in HCT116 and RKO colorectal cancer cells. Because M.CviPI-probed samples contain both CG and GC methylation, we developed a versatile, visually-intuitive program, called MethylViewer, for evaluating the bisulfite sequencing results. Uniquely, MethylViewer can simultaneously query cytosine methylation status in bisulfite-converted sequences at as many as four different user-defined motifs, e.g. CG, GC, etc., including motifs with degenerate bases. Data can also be exported for statistical analysis and as publication-quality images. Analysis of hMLH1 MAPit data with MethylViewer showed that endogenous CG methylation and accessible GC sites were both mapped on single molecules at high resolution. Disruption of positioned nucleosomes on single molecules of the PHO5 promoter was detected in budding yeast using M.CviPII, increasing the number of enzymes available for probing protein-DNA interactions. MethylViewer provides an integrated solution for primer design and rapid, accurate and detailed analysis of bisulfite sequencing or MAPit datasets from virtually any biological or biochemical system.

Bisulfite sequencing of DNA

Exact positions of 5-methylcytosine (m(5)C) on a single strand of DNA can be determined by bisulfite genomic sequencing (BGS). Treatment with bisulfite ion preferentially deaminates unmethylated cytosines, which are then converted to uracil upon desulfonation. Amplifying regions of interest from deaminated DNA and sequencing products cloned from amplicons permits determination of methylation at single-nucleotide resolution along single DNA molecules, which is not possible with other methylation analysis techniques. This unit describes a BGS technique suitable for most DNA sources, including formaldehyde-fixed tissue. Considerations for experimental design and common sources of error are discussed.

Gene regulation by nucleosome positioning

To achieve high compaction, most genomic DNA in eukaryotes is incorporated into nucleosomes; however, regulatory factors and transcriptional machinery must gain access to chromatin to extract genetic information. This conflict is partially resolved by a particular arrangement of nucleosome locations on the genome. Across all eukaryotic species, promoters and other regulatory sequences are more nucleosome-depleted, whereas transcribed regions tend to be occupied with well-positioned, high-density nucleosomal arrays. This nucleosome positioning pattern, as well as its dynamic regulation, facilitates the access of transcription factors to their target sites and plays a crucial role in determining the transcription level, cell-to-cell variation and activation or repression dynamics.

Nucleosome interactions and stability in an ordered nucleosome array model system

Although it is well established that the majority of eukaryotic DNA is sequestered as nucleosomes, the higher-order structure resulting from nucleosome interactions as well as the dynamics of nucleosome stability are not as well understood. To characterize the structural and functional contribution of individual nucleosomal sites, we have developed a chromatin model system containing up to four nucleosomes, where the array composition, saturation, and length can be varied via the ordered ligation of distinct mononucleosomes. Using this system we find that the ligated tetranucleosomal arrays undergo intra-array compaction. However, this compaction is less extensive than for longer arrays and is histone H4 tail-independent, suggesting that well ordered stretches of four or fewer nucleosomes do not fully compact to the 30-nm fiber. Like longer arrays, the tetranucleosomal arrays exhibit cooperative self-association to form species composed of many copies of the array. This propensity for self-association decreases when the fraction of nucleosomes lacking H4 tails is systematically increased. However, even tetranucleosomal arrays with only two octamers possessing H4 tails recapitulate most of the inter-array self-association. Varying array length shows that systems as short as dinucleosomes demonstrate significant self-association, confirming that relatively few determinants are required for inter-array interactions and suggesting that in vivo multiple interactions of short runs of nucleosomes might contribute to complex fiber-fiber interactions. Additionally, we find that the stability of nucleosomes toward octamer loss increases with array length and saturation, suggesting that in vivo stretches of ordered, saturated nucleosomes could serve to protect these regions from histone ejection.

In vitro reconstitution of PHO5 promoter chromatin remodeling points to a role for activator-nucleosome competition in vivo

The yeast PHO5 promoter is a classical model for studying the role of chromatin in gene regulation. To enable biochemical dissection of the mechanism leading to PHO5 activation, we reconstituted the process in vitro. Positioned nucleosomes corresponding to the repressed PHO5 promoter state were assembled using a yeast extract-based in vitro system. Addition of the transactivator Pho4 yielded an extensive DNase I-hypersensitive site resembling induced PHO5 promoter chromatin. Importantly, this remodeling was energy dependent. In contrast, little or no chromatin remodeling was detected at the PHO8 or PHO84 promoter in this in vitro system. Only the PHO5 promoter harbors a high-affinity intranucleosomal Pho4 binding site (UASp) where Pho4 binding can compete with nucleosome formation, prompting us to test the importance of such competition for chromatin remodeling by analysis of UASp mutants in vivo. Indeed, the intranucleosomal location of the UASp element was critical, but not essential, for complete remodeling at the PHO5 promoter in vivo. Further, binding of just the Gal4 DNA binding domain to an intranucleosomal site could increase PHO5 promoter opening. These data establish an auxiliary role for DNA binding competition between Pho4 and histones in PHO5 promoter chromatin remodeling in vivo.

SWI/SNF has intrinsic nucleosome disassembly activity that is dependent on adjacent nucleosomes

The ATP-dependent chromatin remodeling complex SWI/SNF regulates transcription and has been implicated in promoter nucleosome eviction. Efficient nucleosome disassembly by SWI/SNF alone in biochemical assays, however, has not been directly observed. Employing a model system of dinucleosomes rather than mononucleosomes, we demonstrate that remodeling leads to ordered and efficient disassembly of one of the two nucleosomes. An H2A/H2B dimer is first rapidly displaced, and then, in a slower reaction, an entire histone octamer is lost. Nucleosome disassembly by SWI/SNF did not require additional factors such as chaperones or acceptors of histones. Observations in single molecules as well as bulk measurement suggest that a key intermediate in this process is one in which a nucleosome is moved toward the adjacent nucleosome. SWI/SNF recruited by the transcriptional activator Gal4-VP16 preferentially mobilizes the proximal nucleosome and destabilizes the adjacent nucleosome.

Nucleosome-depleted regions in cell-cycle-regulated promoters ensure reliable gene expression in every cell cycle

Many promoters in eukaryotes have nucleosome-depleted regions (NDRs) containing transcription factor binding sites. However, the functional significance of NDRs is not well understood. Here, we examine NDR function in two cell cycle-regulated promoters, CLN2pr and HOpr, by varying nucleosomal coverage of the binding sites of their activator, Swi4/Swi6 cell-cycle box (SCB)-binding factor (SBF), and probing the corresponding transcriptional activity in individual cells with time-lapse microscopy. Nucleosome-embedded SCBs do not significantly alter peak expression levels. Instead, they induce bimodal, "on/off" activation in individual cell cycles, which displays short-term memory, or epigenetic inheritance, from the mother cycle. In striking contrast, the same SCBs localized in NDR lead to highly reliable activation, once in every cell cycle. We further demonstrate that the high variability in Cln2p expression induced by the nucleosomal SCBs reduces cell fitness. Therefore, we propose that the NDR function in limiting stochasticity in gene expression promotes the ubiquity and conservation of promoter NDR. PAPERCLIP:

SWI/SNF and Asf1p cooperate to displace histones during induction of the saccharomyces cerevisiae HO promoter

Regulation of the Saccharomyces cerevisiae HO promoter has been shown to require the recruitment of chromatin-modifying and -remodeling enzymes. Despite this, relatively little is known about what changes to chromatin structure occur during the course of regulation at HO. Here, we used indirect end labeling in synchronized cultures to show that the chromatin structure is disrupted in a region that spans bp -600 to -1800 relative to the transcriptional start site. Across this region, there is a loss of canonical nucleosomes and a reduction in histone DNA cross-linking, as monitored by chromatin immunoprecipitation. The ATPase Snf2 is required for these alterations, but the histone acetyltransferase Gcn5 is not. This suggests that the SWI/SNF complex is directly involved in nucleosome removal at HO. We also present evidence indicating that the histone chaperone Asf1 assists in this. These observations suggest that SWI/SNF-related complexes in concert with histone chaperones act to remove histone octamers from DNA during the course of gene regulation.

Single-molecule analysis of chromatin: changing the view of genomes one molecule at a time

Wrapping DNA into chromatin provides a wealth of regulatory mechanisms that ensure normal growth and development in eukaryotes. Our understanding of chromatin structure, including nucleosomes and non-histone protein-DNA interactions, has benefited immensely from nuclease and chemical digestion techniques. DNA-bound proteins, such as histones or site-specific factors, protect DNA against nuclease cleavage and generate large nucleosomal or small regulatory factor footprints. Chromatin subject to distinct modes of regulation often coincides with sites of nuclease hypersensitivity or nucleosome positioning. An inherent limitation of cleavage-based analyses has been the inability to reliably analyze regions of interest when levels of digestion depart from single-hit kinetics. Moreover, cleavage-based techniques provide views that are averaged over all the molecules in a sample population. Therefore, in cases of occupancy of multiple regulatory elements by factors, one cannot define whether the factors are bound to the same or different molecules in the population. The recent development of DNA methyltransferase-based, single-molecule MAP-IT technology overcomes limitations of ensemble approaches and has opened numerous new avenues in chromatin research. Here, we review the strengths, limitations, applications and future prospects of MAP-IT ranging from structural issues to mechanistic questions in eukaryotic chromatin regulation.

DNA methyltransferase probing of chromatin structure within populations and on single molecules

Non-invasive methods for mapping chromatin structure are necessary for creating an accurate view of genome function and dynamics in vivo. Ectopic induction of cytosine-5 DNA methyltransferases (C5 MTases) in Saccharomyces cerevisiae is a powerful technique for probing chromatin structure with minimal disruption to yeast physiology. Accessibility of MTases to their cognate sites is impaired based on the strength and span of the protein-DNA interaction to be probed. Methylated cytosines that resist chemical deamination are detected positively by the PCR-based technique of bisulfite genomic sequencing. PCR amplicons can be sequenced directly yielding an average m(5)C frequency or accessibility of each target site within the population, a technique termed methyltransferase accessibility protocol (MAP). More recently, the sequencing of cloned molecules in MAP for individual templates (MAPit) enables assignment of the methylation status of each target site along a continuous DNA strand from a single cell. The unique capability to score methylation at multiple sites in single molecules permits detection of inherent structural variability in chromatin. Here, MAPit analysis of the repressed and induced PHO5 promoter of budding yeast, using a C5 MTase with dinucleotide recognition specificity, reveals considerable cell-to-cell heterogeneity in chromatin structure. Substantial variation is observed in the extent to which the MTase gains entry to each of the nucleosomes positioned at PHO5, suggesting differences in their intrinsic thermodynamic stability in vivo. MAPit should be readily adaptable to the analysis of chromatin structure and non-histone protein-DNA interactions in a variety of model systems.

A quantitative model of transcription factor-activated gene expression

A challenge facing biology is to develop quantitative, predictive models of gene regulation. Eukaryotic promoters contain transcription factor binding sites of differing affinity and accessibility, but we understand little about how these variables combine to generate a fine-tuned, quantitative transcriptional response. Here we used the PHO5 promoter in budding yeast to quantify the relationship between transcription factor input and gene expression output, termed the gene-regulation function (GRF). A model that captures variable interactions between transcription factors, nucleosomes and the promoter faithfully reproduced the observed quantitative changes in the GRF that occur upon altering the affinity of transcription factor binding sites, and implicates nucleosome-modulated accessibility of transcription factor binding sites in increasing the diversity of gene expression profiles. This work establishes a quantitative framework that can be applied to predict GRFs of other eukaryotic genes.

Nucleosome retention and the stochastic nature of promoter chromatin remodeling for transcription

The rate-limiting step of transcriptional activation in eukaryotes, and thus the critical point for gene regulation, is unknown. Combining biochemical analyses of the chromatin transition at the transcriptionally induced PHO5 promoter in yeast with modeling based on a small number of simple assumptions, we demonstrate that random removal and reformation of promoter nucleosomes can account for stochastic and kinetic properties of PHO5 expression. Our analysis suggests that the disassembly of promoter nucleosomes is rate limiting for PHO5 expression, and supports a model for the underlying mechanism of promoter chromatin remodeling, which appears to conserve a single nucleosome on the promoter at all times.

Mediator requirement downstream of chromatin remodeling during transcriptional activation of CHA1 in yeast

Mediator complex is essential for transcription by RNA polymerase II in eukaryotes. Although chromatin remodeling is an integral part of transcriptional activation at many promoters, whether Mediator is required for this function has not been determined. Here we have used the yeast CHA1 gene to study the role of Mediator in chromatin remodeling and recruitment of the transcription machinery. We show by chromatin immunoprecipitation that Mediator subunits are recruited to the induced CHA1 promoter. Inactivation of Mediator at 37 degrees C in yeast harboring the srb4-138 (med17) ts mutation severely reduces CHA1 activation and prevents recruitment to the induced CHA1 promoter of Med18/Srb5, from the head module of Mediator, and Med14/Rgr1, which bridges the middle and tail modules. In contrast, recruitment of Med15/Gal11 from the tail module is unaffected in med17 ts yeast at 37 degrees C. Recruitment of TATA-binding protein (TBP) is severely compromised in the absence of functional Mediator, whereas Kin28 and polymerase II recruitment are reduced but to a lesser extent. Induced levels of histone H3K4me3 at the CHA1 promoter are not diminished by inactivation of Mediator, whereas recruitment of Paf1 and of Ser2- and Ser5-phosphorylated forms of Rbp1 are reduced but not eliminated. Loss of histone H3 from the induced CHA1 promoter is seen in wild type yeast but is greatly reduced by loss of intact Mediator. In contrast, Swi/Snf recruitment and nucleosome remodeling are unaffected by loss of Mediator function. Thus, Mediator is required for recruitment of the transcription machinery subsequent to chromatin remodeling during CHA1 induction.

Engineered apoptotic nucleases for chromatin research

We have created new genomics tools for chromatin research by genetically engineering the human and mouse major apoptotic nucleases that are responsible for internucleosomal DNA cleavage, DNA fragmentation factor (DFF). Normally, in its inactive form, DFF is a heterodimer composed of a 45-kDa chaperone inhibitor subunit (DFF45 or ICAD), and a 40-kDa latent endonuclease subunit (DFF40 or CAD). Upon caspase-3 cleavage of DFF45, DFF40 forms active endonuclease homo-oligomers. Although Saccharomyces cerevisiae lacks DFF, expression of caspase-3 is lethal in this organism, but expression of the highly sequence-specific tobacco etch virus protease (TEVP) is harmless. Therefore, we inserted TEVP cleavage sites immediately downstream of the two caspase-3 cleavage sites within DFF45, generating a novel form of DFF (DFF-T) whose nuclease activity proved to be exclusively under the control of TEVP. We demonstrate that co-expression of TEVP and DFF-T under galactose control results in nucleosomal DNA laddering and cell death in S. cerevisiae. We also created synthetic DFF genes with optimized codons for high-level expression in Eschericia coli or S. cerevisiae. We further demonstrate the excellence of the synthetic gene products for in vitro mapping of the nucleosome positions and hypersensitive sites in specific genes such as the yeast PHO5.

Chromatin disassembly from the PHO5 promoter is essential for the recruitment of the general transcription machinery and coactivators

The disassembly of promoter nucleosomes appears to be a general property of highly transcribed eukaryotic genes. We have previously shown that the disassembly of chromatin from the promoters of the Saccharomyces cerevisiae PHO5 and PHO8 genes, mediated by the histone chaperone anti-silencing function 1 (Asf1), is essential for transcriptional activation upon phosphate depletion. This mechanism of transcriptional regulation is shared with the ADY2 and ADH2 genes upon glucose removal. Promoter chromatin disassembly by Asf1 is required for recruitment of TBP and RNA polymerase II, but not the Pho4 and Pho2 activators. Furthermore, accumulation of SWI/SNF and SAGA at the PHO5 promoter requires promoter chromatin disassembly. By contrast, the requirement for SWI/SNF and SAGA to facilitate Pho4 activator recruitment to the nucleosome-buried binding site in the PHO5 promoter occurs prior to chromatin disassembly and is distinct from the stable recruitment of SWI/SNF and SAGA that occurs after chromatin disassembly.

Continuous histone H2B and transcription-dependent histone H3 exchange in yeast cells outside of replication

We investigated the dynamics of histone-DNA interactions in yeast by using inducible forms of epitope-tagged histones H2B and H3. Chromatin assembly of newly synthesized histones was assessed by chromatin immunoprecipitation in G1-arrested cells to prevent replication-coupled histone incorporation. We find that while histone deposition within a subtelomeric region is strictly linked to DNA replication, histone H2B is continuously incorporated at the promoter and coding regions of both transcriptionally active and inactive loci. In contrast, incorporation of histone H3 occurs only at active genes, being predominant at the promoter and showing a dynamics along the gene that inversely correlates with the average nucleosomal density. Similar results were obtained with N-terminally truncated H2B and H3 variants. We infer that replication-independent incorporation of H2B and H3 are distinct events, each occurring independently of the histone tail, and that nucleosome loss at active promoters reflects a dynamic equilibrium between histone deposition and dissociation.

Single-molecule and population probing of chromatin structure using DNA methyltransferases

Probing chromatin structure with DNA methyltransferases offers advantages over more commonly used nuclease-based and chromatin immunoprecipitation methods for detection of nucleosomes and non-histone protein-DNA interactions. Here, we describe two related methods in which the readout of MTase accessibility is obtained by assaying 5-methylcytosine in DNA through the PCR-based technique of bisulfite genomic sequencing. The methyltransferase accessibility protocol (MAP) determines the relative frequency at which the enzyme accesses each of its target sites over an entire population of PCR amplified product. While MAP yields much quantitative information about relative accessibility of a region of chromatin, a complementary single-molecule view of methyltransferase accessibility, termed MAP for individual templates (MAP-IT), is provided by analysis of cloned PCR products. Absolute rather than relative methylation frequencies in a region are obtained by summing the methylation status at each site over a cohort of clones. Moreover, as the integrity of individual molecules is maintained in MAP-IT, unique information about the distribution of multiple footprints along continuous regions is gleaned. In principle, the population MAP and single-molecule MAP-IT strategies can be used to analyze chromatin structure in a variety of model systems. Here, we describe the application of MAP in living Saccharomyces cerevisiae cells and MAP-IT in the analysis of a mammalian tumor suppressor gene in nuclei. This application of MAP-IT provides the first means to simultaneously determine CpG methylation of mammalian genes and their overlying chromatin structure in the same single DNA molecule.

Constitutive nucleosome depletion and ordered factor assembly at the GRP78 promoter revealed by single molecule footprinting

Chromatin organization and transcriptional regulation are interrelated processes. A shortcoming of current experimental approaches to these complex events is the lack of methods that can capture the activation process on single promoters. We have recently described a method that combines methyltransferase M.SssI treatment of intact nuclei and bisulfite sequencing allowing the representation of replicas of single promoters in terms of protected and unprotected footprint modules. Here we combine this method with computational analysis to study single molecule dynamics of transcriptional activation in the stress inducible GRP78 promoter. We show that a 350–base pair region upstream of the transcription initiation site is constitutively depleted of nucleosomes, regardless of the induction state of the promoter, providing one of the first examples for such a promoter in mammals. The 350–base pair nucleosome-free region can be dissected into modules, identifying transcription factor binding sites and their combinatorial organization during endoplasmic reticulum stress. The interaction of the transcriptional machinery with the GRP78 core promoter is highly organized, represented by six major combinatorial states. We show that the TATA box is frequently occupied in the noninduced state, that stress induction results in sequential loading of the endoplasmic reticulum stress response elements, and that a substantial portion of these elements is no longer occupied following recruitment of factors to the transcription initiation site. Studying the positioning of nucleosomes and transcription factors at the single promoter level provides a powerful tool to gain novel insights into the transcriptional process in eukaryotes.

Control of gene expression and transcription are complex and well-coordinated processes. Most current experimental approaches to understanding the underlying mechanisms, which include binding of transcription factors to regulatory regions of genes, and changes in the structure and composition of chromatin, rely on studies of populations of cells and cannot capture the transcription activation process on single promoters. The authors describe the use of a footprinting method which enables analysis of chromatin structure and binding of factors on single DNA molecules. This is applied to study the activation process of GRP78, a protein which is important for the induction of a response to endoplasmic reticulum stress. By combining the footprinting method and computational analyses, the authors define functional modules on the GRP78 promoter and show that it exists in few major combinatorial states, reflecting its high level of organization. These results provide novel insights into the activation of GRP78 which could not be gleaned using conventional methods. They also demonstrate the use of the method as a unique and powerful tool to study the transcriptional process in eukaryotes, which remains a major source of interest and challenge for the scientific community.