Genome-scale identification of nucleosome positions in S. cerevisiae

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

Global nucleosome positioning regulates salicylic acid mediated transcription in Arabidopsis thaliana

BACKGROUND

The nucleosome positioning regulates the gene expression and many other DNA-related processes in eukaryotes. Genome-wide mapping of nucleosome positions and correlation of genome-wide nucleosomal remodeling with the changes in the gene expression can help us understanding gene regulation on genome level.

RESULTS

In the present study, we correlate the gene expression and the genomic nucleosomal remodeling in response to salicylic acid (SA) treatment in A. thaliana. We have mapped genome-wide nucleosomes by performing tiling microarray using 146 bp mononucleosomal template DNA. The average nucleosomal coverage is approximately 346 bp per nucleosome both under the control and the SA-treated conditions. The nucleosomal coverage is more in the coding region than in the 5' regulatory regions. We observe approximately 50% nucleosomal remodeling on SA treatment where significant nucleosomal depletion and nucleosomal enrichment around the transcription start site (TSS) occur in SA induced genes and SA repressed genes respectively in response to SA treatment. Especially in the case of the SA-induced group, the nucleosomal remodeling over the minimal promoter in response to SA is especially significant in the Non-expresser of PR1 (NPR1)-dependent genes. A detailed investigation of npr1-1 mutant confirms a distinct role of NPR1 in the nucleosome remodeling over the core promoter. We have also identified several motifs for various hormonal responses; including ABRE elements in the remodeled nucleosomal regions around the promoter region in the SA regulated genes. We have further identified that the W-box and TGACG/C motif, reported to play an important role in SA-mediated induction, are enriched in nucleosome free regions (NFRs) of the promoter region of the SA induced genes.

CONCLUSIONS

This is the first study reporting genome-wide effects of SA treatment on the chromatin architecture of A. thaliana. It also reports significant role of NPR1 in genome-wide nucleosomal remodeling in response to SA.

Role of DNA sequence in chromatin remodeling and the formation of nucleosome-free regions

AT-rich DNA is concentrated in the nucleosome-free regions (NFRs) associated with transcription start sites of most genes. Lorch et al. find that the AT-rich sequences present in many NFRs have little effect on the stability of nucleosomes. These sequences instead facilitate the removal of nucleosomes by the RSC chromatin remodeling complex. RSC activity is stimulated by AT-rich sequences in nucleosomes and inhibited by competition with AT-rich DNA.

AT-rich DNA is concentrated in the nucleosome-free regions (NFRs) associated with transcription start sites of most genes. We tested the hypothesis that AT-rich DNA engenders NFR formation by virtue of its rigidity and consequent exclusion of nucleosomes. We found that the AT-rich sequences present in many NFRs have little effect on the stability of nucleosomes. Rather, these sequences facilitate the removal of nucleosomes by the RSC chromatin remodeling complex. RSC activity is stimulated by AT-rich sequences in nucleosomes and inhibited by competition with AT-rich DNA. RSC may remove NFR nucleosomes without effect on adjacent ORF nucleosomes. Our findings suggest that many NFRs are formed and maintained by an active mechanism involving the ATP-dependent removal of nucleosomes rather than a passive mechanism due to the intrinsic instability of nucleosomes on AT-rich DNA sequences.

A high-throughput ChIP-Seq for large-scale chromatin studies

We present a modified approach of chromatin immuno-precipitation followed by sequencing (ChIP-Seq), which relies on the direct ligation of molecular barcodes to chromatin fragments, thereby permitting experimental scale-up. With Bar-ChIP now enabling the concurrent profiling of multiple DNA–protein interactions, we report the simultaneous generation of 90 ChIP-Seq datasets without any robotic instrumentation. We demonstrate that application of Bar-ChIP to a panel of Saccharomyces cerevisiae chromatin-associated mutants provides a rapid and accurate genome-wide overview of their chromatin status. Additionally, we validate the utility of this technology to derive novel biological insights by identifying a role for the Rpd3S complex in maintaining H3K14 hypo-acetylation in gene bodies. We also report an association between the presence of intragenic H3K4 tri-methylation and the emergence of cryptic transcription in a Set2 mutant. Finally, we uncover a crosstalk between H3K14 acetylation and H3K4 methylation in this mutant. These results show that Bar-ChIP enables biological discovery through rapid chromatin profiling at single-nucleosome resolution for various conditions and protein modifications at once.

ChIP-Seq to Analyze the Binding of Replication Proteins to Chromatin

Chromatin immunoprecipitation (ChIP) is a widely used method to study interactions between proteins and discrete chromosomal loci in vivo. ChIP was originally developed for in vivo analysis of protein associations with candidate DNA sequences known or suspected to bind the protein of interest. The advent of DNA microarrays enabled the unbiased, genome-scale identification of all DNA sequences enriched by ChIP, providing a genomic map of a protein's chromatin binding. This method, termed ChIP-chip, is broadly applicable and has been particularly valuable in DNA replication studies to map potential replication origins in Saccharomyces cerevisiae and other organisms based on the specific association of certain replication proteins with these chromosomal elements, which are distributed throughout the genome. More recently, high-throughput sequencing (HTS) technologies have replaced microarrays as the preferred method for genomic analysis of ChIP experiments, and this combination is termed ChIP-Seq. We present a detailed ChIP-Seq protocol for S. cerevisiae that can be adapted for different HTS platforms and for different organisms. We also outline general schemes for data analysis; however, HTS data analyses usually must be tailored specifically for individual studies, depending on the experimental design, data characteristics, and the genome being analyzed.

Nucleosome Positioning with Set of Key Positions and Nucleosome Affinity

The formation and precise positioning of nucleosome in chromatin occupies a very important role in studying life process. Today, there are many researchers who discovered that the positioning where the location of a DNA sequence fragment wraps around a histone octamer in genome is not random but regular. However, the positioning is closely relevant to the concrete sequence of core DNA. So in this paper, we analyzed the relation between the affinity and sequence structure of core DNA, and extracted the set of key positions. In these positions, the nucleotide sequences probably occupy mainly action in the binding. First, we simplified and formatted the experimental data with the affinity. Then, to find the key positions in the wrapping, we used neural network to analyze the positive and negative effects of nucleosome generation for each position in core DNA sequences. However, we reached a class of weights with every position to describe this effect. Finally, based on the positions with high weights, we analyzed the reason why the chosen positions are key positions, and used these positions to construct a model for nucleosome positioning prediction. Experimental results show the effectiveness of our method.

Nucleosome positioning in yeasts: methods, maps, and mechanisms

Eukaryotic nuclear DNA is packaged into nucleosomes. During the past decade, genome-wide nucleosome mapping across species revealed the high degree of order in nucleosome positioning. There is a conserved stereotypical nucleosome organization around transcription start sites (TSSs) with a nucleosome-depleted region (NDR) upstream of the TSS and a TSS-aligned regular array of evenly spaced nucleosomes downstream over the gene body. As nucleosomes largely impede access to DNA and thereby provide an important level of genome regulation, it is of general interest to understand the mechanisms generating nucleosome positioning and especially the stereotypical NDR-array pattern. We focus here on the most advanced models, unicellular yeasts, and review the progress in mapping nucleosomes and which nucleosome positioning mechanisms are discussed. There are four mechanistic aspects: How are NDRs generated? How are individual nucleosomes positioned, especially those flanking the NDRs? How are nucleosomes evenly spaced leading to regular arrays? How are regular arrays aligned at TSSs? The main candidates for nucleosome positioning determinants are intrinsic DNA binding preferences of the histone octamer, specific DNA binding factors, nucleosome remodeling enzymes, transcription, and statistical positioning. We summarize the state of the art in an integrative model where nucleosomes are positioned by a combination of all these candidate determinants. We highlight the predominance of active mechanisms involving nucleosome remodeling enzymes which may be recruited by DNA binding factors and the transcription machinery. While this mechanistic framework emerged clearly during recent years, the involved factors and their mechanisms are still poorly understood and require future efforts combining in vivo and in vitro approaches.

Replication-guided nucleosome packing and nucleosome breathing expedite the formation of dense arrays

The first level of genome packaging in eukaryotic cells involves the formation of dense nucleosome arrays, with DNA coverage near 90% in yeasts. How cells achieve such high coverage within a short time, e.g. after DNA replication, remains poorly understood. It is known that random sequential adsorption of impenetrable particles on a line reaches high density extremely slowly, due to a jamming phenomenon. The nucleosome-shifting action of remodeling enzymes has been proposed as a mechanism to resolve such jams. Here, we suggest two biophysical mechanisms which assist rapid filling of DNA with nucleosomes, and we quantitatively characterize these mechanisms within mathematical models. First, we show that the 'softness' of nucleosomes, due to nucleosome breathing and stepwise nucleosome assembly, significantly alters the filling behavior, speeding up the process relative to 'hard' particles with fixed, mutually exclusive DNA footprints. Second, we explore model scenarios in which the progression of the replication fork could eliminate nucleosome jamming, either by rapid filling in its wake or via memory of the parental nucleosome positions. Taken together, our results suggest that biophysical effects promote rapid nucleosome filling, making the reassembly of densely packed nucleosomes after DNA replication a simpler task for cells than was previously thought.

Unbiased chromatin accessibility profiling by RED-seq uncovers unique features of nucleosome variants in vivo

BACKGROUND

Differential accessibility of DNA to nuclear proteins underlies the regulation of numerous cellular processes. Although DNA accessibility is primarily determined by the presence or absence of nucleosomes, differences in nucleosome composition or dynamics may also regulate accessibility. Methods for mapping nucleosome positions and occupancies genome-wide (MNase-seq) have uncovered the nucleosome landscapes of many different cell types and organisms. Conversely, methods specialized for the detection of large nucleosome-free regions of chromatin (DNase-seq, FAIRE-seq) have uncovered numerous gene regulatory elements. However, these methods are less successful in measuring the accessibility of DNA sequences within nucelosome arrays.

RESULTS

Here we probe the genome-wide accessibility of multiple cell types in an unbiased manner using restriction endonuclease digestion of chromatin coupled to deep sequencing (RED-seq). Using this method, we identified differences in chromatin accessibility between populations of cells, not only in nucleosome-depleted regions of the genome (e.g., enhancers and promoters), but also within the majority of the genome that is packaged into nucleosome arrays. Furthermore, we identified both large differences in chromatin accessibility in distinct cell lineages and subtle but significant changes during differentiation of mouse embryonic stem cells (ESCs). Most significantly, using RED-seq, we identified differences in accessibility among nucleosomes harboring well-studied histone variants, and show that these differences depend on factors required for their deposition.

CONCLUSIONS

Using an unbiased method to probe chromatin accessibility genome-wide, we uncover unique features of chromatin structure that are not observed using more widely-utilized methods. We demonstrate that different types of nucleosomes within mammalian cells exhibit different degrees of accessibility. These findings provide significant insight into the regulation of DNA accessibility.

cisMEP: an integrated repository of genomic epigenetic profiles and cis-regulatory modules in Drosophila

BACKGROUND

Cis-regulatory modules (CRMs), or the DNA sequences required for regulating gene expression, play the central role in biological researches on transcriptional regulation in metazoan species. Nowadays, the systematic understanding of CRMs still mainly resorts to computational methods due to the time-consuming and small-scale nature of experimental methods. But the accuracy and reliability of different CRM prediction tools are still unclear. Without comparative cross-analysis of the results and combinatorial consideration with extra experimental information, there is no easy way to assess the confidence of the predicted CRMs. This limits the genome-wide understanding of CRMs.

DESCRIPTION

It is known that transcription factor binding and epigenetic profiles tend to determine functions of CRMs in gene transcriptional regulation. Thus integration of the genome-wide epigenetic profiles with systematically predicted CRMs can greatly help researchers evaluate and decipher the prediction confidence and possible transcriptional regulatory functions of these potential CRMs. However, these data are still fragmentary in the literatures. Here we performed the computational genome-wide screening for potential CRMs using different prediction tools and constructed the pioneer database, cisMEP (cis-regulatory module epigenetic profile database), to integrate these computationally identified CRMs with genomic epigenetic profile data. cisMEP collects the literature-curated TFBS location data and nine genres of epigenetic data for assessing the confidence of these potential CRMs and deciphering the possible CRM functionality.

CONCLUSIONS

cisMEP aims to provide a user-friendly interface for researchers to assess the confidence of different potential CRMs and to understand the functions of CRMs through experimentally-identified epigenetic profiles. The deposited potential CRMs and experimental epigenetic profiles for confidence assessment provide experimentally testable hypotheses for the molecular mechanisms of metazoan gene regulation. We believe that the information deposited in cisMEP will greatly facilitate the comparative usage of different CRM prediction tools and will help biologists to study the modular regulatory mechanisms between different TFs and their target genes.

Long non-coding RNA-mediated transcriptional interference of a permease gene confers drug tolerance in fission yeast

Most long non-coding RNAs (lncRNAs) encoded by eukaryotic genomes remain uncharacterized. Here we focus on a set of intergenic lncRNAs in fission yeast. Deleting one of these lncRNAs exhibited a clear phenotype: drug sensitivity. Detailed analyses of the affected locus revealed that transcription of the nc-tgp1 lncRNA regulates drug tolerance by repressing the adjacent phosphate-responsive permease gene transporter for glycerophosphodiester 1 (tgp1⁺). We demonstrate that the act of transcribing nc-tgp1 over the tgp1⁺ promoter increases nucleosome density, prevents transcription factor access and thus represses tgp1⁺ without the need for RNA interference or heterochromatin components. We therefore conclude that tgp1⁺ is regulated by transcriptional interference. Accordingly, decreased nc-tgp1 transcription permits tgp1⁺ expression upon phosphate starvation. Furthermore, nc-tgp1 loss induces tgp1⁺ even in repressive conditions. Notably, drug sensitivity results directly from tgp1⁺ expression in the absence of the nc-tgp1 RNA. Thus, transcription of an lncRNA governs drug tolerance in fission yeast.

The presence of long non-coding RNAs (lncRNAs) is pervasive across genomes, yet few lncRNAs have clearly established mechanisms of action. Here the authors demonstrate that the fission yeast lncRNA nc-tgp1 regulates expression of the drug tolerance gene tgp1₊ via⁺ transcriptional interference.

Chromatin accessibility: a window into the genome

Transcriptional activation throughout the eukaryotic lineage has been tightly linked with disruption of nucleosome organization at promoters, enhancers, silencers, insulators and locus control regions due to transcription factor binding. Regulatory DNA thus coincides with open or accessible genomic sites of remodeled chromatin. Current chromatin accessibility assays are used to separate the genome by enzymatic or chemical means and isolate either the accessible or protected locations. The isolated DNA is then quantified using a next-generation sequencing platform. Wide application of these assays has recently focused on the identification of the instrumental epigenetic changes responsible for differential gene expression, cell proliferation, functional diversification and disease development. Here we discuss the limitations and advantages of current genome-wide chromatin accessibility assays with especial attention on experimental precautions and sequence data analysis. We conclude with our perspective on future improvements necessary for moving the field of chromatin profiling forward.

Sequence analysis of origins of replication in the Saccharomyces cerevisiae genomes

DNA replication is a highly precise process that is initiated from origins of replication (ORIs) and is regulated by a set of regulatory proteins. The mining of DNA sequence information will be not only beneficial for understanding the regulatory mechanism of replication initiation but also for accurately identifying ORIs. In this study, the GC profile and GC skew were calculated to analyze the compositional bias in the Saccharomyces cerevisiae genome. We found that the GC profile in the region of ORIs is significantly lower than that in the flanking regions. By calculating the information redundancy, an estimation of the correlation of nucleotides, we found that the intensity of adjoining correlation in ORIs is dramatically higher than that in flanking regions. Furthermore, the relationships between ORIs and nucleosomes as well as transcription start sites were investigated. Results showed that ORIs are usually not occupied by nucleosomes. Finally, we calculated the distribution of ORIs in yeast chromosomes and found that most ORIs are in transcription terminal regions. We hope that these results will contribute to the identification of ORIs and the study of DNA replication mechanisms.

Predicting nucleosome positioning based on geometrically transformed Tsallis entropy

As the fundamental unit of eukaryotic chromatin structure, nucleosome plays critical roles in gene expression and regulation by controlling physical access to transcription factors. In this paper, based on the geometrically transformed Tsallis entropy and two index-vectors, a valid nucleosome positioning information model is developed to describe the distribution of A/T-riched and G/C-riched dimeric and trimeric motifs along the DNA duplex. When applied to train the support vector machine, the model achieves high AUCs across five organisms, which have significantly outperformed the previous studies. Besides, we adopt the concept of relative distance to describe the probability of arbitrary DNA sequence covered by nucleosome. Thus, the average nucleosome occupancy profile over the S.cerevisiae genome is calculated. With our peak detection model, the isolated nucleosomes along genome sequence are located. When compared with some published results, it shows that our model is effective for nucleosome positioning. The index-vector component is identified to be an important influencing factor of nucleosome organizations.

Chd1 co-localizes with early transcription elongation factors independently of H3K36 methylation and releases stalled RNA polymerase II at introns

Background

Chromatin consists of ordered nucleosomal arrays that are controlled by highly conserved adenosine triphosphate (ATP)-dependent chromatin remodeling complexes. One such remodeler, chromodomain helicase DNA binding protein 1 (Chd1), is believed to play an integral role in nucleosomal organization, as the loss of Chd1 is known to disrupt chromatin. However, the specificity and basis for the functional and physical localization of Chd1 on chromatin remains largely unknown.

Results

Using genome-wide approaches, we found that the loss of Chd1 significantly disrupted nucleosome arrays within the gene bodies of highly transcribed genes. We also found that Chd1 is physically recruited to gene bodies, and that its occupancy specifically corresponds to that of the early elongating form of RNA polymerase, RNAPII Ser 5-P. Conversely, RNAPII Ser 5-P occupancy was affected by the loss of Chd1, suggesting that Chd1 is associated with early transcription elongation. Surprisingly, the occupancy of RNAPII Ser 5-P was affected by the loss of Chd1 specifically at intron-containing genes. Nucleosome turnover was also affected at these sites in the absence of Chd1. We also found that deletion of the histone methyltransferase for H3K36 (SET2) did not affect either Chd1 occupancy or nucleosome organization genome-wide.

Conclusions

Chd1 is specifically recruited onto the gene bodies of highly transcribed genes in an elongation-dependent but H3K36me3-independent manner. Chd1 co-localizes with the early elongating form of RNA polymerase, and affects the occupancy of RNAPII only at genes containing introns, suggesting a role in relieving splicing-related pausing of RNAPII.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-8935-7-32) contains supplementary material, which is available to authorized users.

An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding

We present a computational model of transcription factor motion that explains both the observed rapid target finding of transcription factors, and how this motion influences protein and genome structure. Using the Smoldyn software, we modelled transcription factor motion arising from a combination of unrestricted 3D diffusion in the nucleoplasm, sliding along the DNA filament, and transferring directly between filament sections by intersegmental transfer. This presents a fine-grain picture of the way in which transcription factors find their targets two orders of magnitude faster than 3D diffusion alone allows. Eukaryotic genomes contain sections of nucleosome free regions (NFRs) around the promoters; our model shows that the presence and size of these NFRs can be explained as their acting as antennas on which transcription factors slide to reach their targets. Additionally, our model shows that intersegmental transfer may have shaped the quaternary structure of transcription factors: sequence specific DNA binding proteins are unusually enriched in dimers and tetramers, perhaps because these allow intersegmental transfer, which accelerates target site finding. Finally, our model shows that a ‘hopping’ motion can emerge from 3D diffusion on small scales. This explains the apparently long sliding lengths that have been observed for some DNA binding proteins observed in vitro. Together, these results suggest that transcription factor diffusion dynamics help drive the evolution of protein and genome structure.

High-resolution digital profiling of the epigenome

The widespread adoption of short-read DNA sequencing as a digital epigenomic readout platform has motivated the development of genome-wide tools that achieve base-pair resolution. New methods for footprinting and affinity purification of nucleosomes, RNA polymerases, chromatin remodellers and transcription factors have increased the resolution of epigenomic profiling by two orders of magnitude, leading to new insights into how the chromatin landscape affects gene regulation. These digital epigenomic tools have also been applied to directly profile both turnover kinetics and transcription in situ. In this Review, we describe how these new genome-wide tools allow interrogation of diverse aspects of the epigenome.

Genome-Wide Analysis of Nucleosome Positions, Occupancy, and Accessibility in Yeast: Nucleosome Mapping, High-Resolution Histone ChIP, and NCAM

Because histones bind DNA very tightly, the location on DNA and the level of occupancy of a given DNA sequence by nucleosomes can profoundly affect accessibility of non-histone proteins to chromatin, affecting virtually all DNA-dependent processes, such as transcription, DNA repair, DNA replication and recombination. Therefore, it is often necessary to determine positions and occupancy of nucleosomes to understand how DNA-dependent processes are regulated. Recent technological advances made such analyses feasible on a genome-wide scale at high resolution. In addition, we have recently developed a method to measure nuclease accessibility of nucleosomes on a global scale. This unit describes methods to map nucleosome positions, to determine nucleosome density, and to determine nuclease accessibility of nucleosomes using deep sequencing.

Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy

BACKGROUND

An organism's DNA sequence is one of the key factors guiding the positioning of nucleosomes within a cell's nucleus. Sequence-dependent bending anisotropy dictates how DNA is wrapped around a histone octamer. One of the best established sequence patterns consistent with this anisotropy is the periodic occurrence of AT-containing dinucleotides (WW) and GC-containing dinucleotides (SS) in the nucleosomal locations where DNA is bent in the minor and major grooves, respectively. Although this simple pattern has been observed in nucleosomes across eukaryotic genomes, its use for prediction of nucleosome positioning was not systematically tested.

RESULTS

We present a simple computational model, termed the W/S scheme, implementing this pattern, without using any training data. This model accurately predicts the rotational positioning of nucleosomes both in vitro and in vivo, in yeast and human genomes. About 65 - 75% of the experimentally observed nucleosome positions are predicted with the precision of one to two base pairs. The program is freely available at http://people.rit.edu/fxcsbi/WS_scheme/. We also introduce a simple and efficient way to compare the performance of different models predicting the rotational positioning of nucleosomes.

CONCLUSIONS

This paper presents the W/S scheme to achieve accurate prediction of rotational positioning of nucleosomes, solely based on the sequence-dependent anisotropic bending of nucleosomal DNA. This method successfully captures DNA features critical for the rotational positioning of nucleosomes, and can be further improved by incorporating additional terms related to the translational positioning of nucleosomes in a species-specific manner.

Active RNAP pre-initiation sites are highly mutated by cytidine deaminases in yeast, with AID targeting small RNA genes

Cytidine deaminases are single stranded DNA mutators diversifying antibodies and restricting viral infection. Improper access to the genome leads to translocations and mutations in B cells and contributes to the mutation landscape in cancer, such as kataegis. It remains unclear how deaminases access double stranded genomes and whether off-target mutations favor certain loci, although transcription and opportunistic access during DNA repair are thought to play a role. In yeast, AID and the catalytic domain of APOBEC3G preferentially mutate transcriptionally active genes within narrow regions, 110 base pairs in width, fixed at RNA polymerase initiation sites. Unlike APOBEC3G, AID shows enhanced mutational preference for small RNA genes (tRNAs, snoRNAs and snRNAs) suggesting a putative role for RNA in its recruitment. We uncover the high affinity of the deaminases for the single stranded DNA exposed by initiating RNA polymerases (a DNA configuration reproduced at stalled polymerases) without a requirement for specific cofactors.

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.

Large replication skew domains delimit GC-poor gene deserts in human

Besides their large-scale organization in isochores, mammalian genomes display megabase-sized regions, spanning both genes and intergenes, where the strand nucleotide composition asymmetry decreases linearly, possibly due to replication activity. These so-called skew-N domains cover about a third of the human genome and are bordered by two skew upward jumps that were hypothesized to compose a subset of "master" replication origins active in the germline. Skew-N domains were shown to exhibit a particular gene organization. Genes with CpG-rich promoters likely expressed in the germline are over represented near the master replication origins, with large genes being co-oriented with replication fork progression, which suggests some coordination of replication and transcription. In this study, we describe another skew structure that covers ∼13% of the human genome and that is bordered by putative master replication origins similar to the ones flanking skew-N domains. These skew-split-N domains have a shape reminiscent of a N, but split in half, leaving in the center a region of null skew whose length increases with domain size. These central regions (median size ∼860 kb) have a homogeneous composition, i.e. both a null and constant skew and a constant and low GC content. They correspond to heterochromatin gene deserts found in low-GC isochores with an average gene density of 0.81 promoters/Mb as compared to 7.73 promoters/Mb genome wide. The analysis of epigenetic marks and replication timing data confirms that, in these late replicating heterochomatic regions, the initiation of replication is likely to be random. This contrasts with the transcriptionally active euchromatin state found around the bordering well positioned master replication origins. Altogether skew-N domains and skew-split-N domains cover about 50% of the human genome.

Nucleosomal occupancy changes locally over key regulatory regions during cell differentiation and reprogramming

Chromatin structure determines DNA accessibility. We compare nucleosome occupancy in mouse and human embryonic stem cells (ESCs), induced-pluripotent stem cells (iPSCs) and differentiated cell types using MNase-seq. To address variability inherent in this technique, we developed a bioinformatic approach to identify regions of difference (RoD) in nucleosome occupancy between pluripotent and somatic cells. Surprisingly, most chromatin remains unchanged; a majority of rearrangements appear to affect a single nucleosome. RoDs are enriched at genes and regulatory elements, including enhancers associated with pluripotency and differentiation. RoDs co-localize with binding sites of key developmental regulators, including the reprogramming factors Klf4, Oct4/Sox2 and c-Myc. Nucleosomal landscapes in ESC enhancers are extensively altered, exhibiting lower nucleosome occupancy in pluripotent cells than in somatic cells. Most changes are reset during reprogramming. We conclude that changes in nucleosome occupancy are a hallmark of cell differentiation and reprogramming and likely identify regulatory regions essential for these processes.

Changes in chromatin structure impact gene expression programs by modulating accessibility to the transcription machinery. Here, West et al. explore differences in nucleosome occupancy between mammalian pluripotent and somatic cells and uncover regulatory regions likely to play key roles in determining cell identity.

Drosophila Brahma complex remodels nucleosome organizations in multiple aspects

ATP-dependent chromatin remodeling complexes regulate nucleosome organizations. In Drosophila, gene Brm encodes the core Brahma complex, the ATPase subunit of SWI/SNF class of chromatin remodelers. Its role in modulating the nucleosome landscape in vivo is unclear. In this study, we knocked down Brm in Drosophila third instar larvae to explore the changes in nucleosome profiles and global gene transcription. The results show that Brm knockdown leads to nucleosome occupancy changes throughout the entire genome with a bias in occupancy decrease. In contrast, the knockdown has limited impacts on nucleosome position shift. The knockdown also alters another important physical property of nucleosome positioning, fuzziness. Nucleosome position shift, gain or loss and fuzziness changes are all enriched in promoter regions. Nucleosome arrays around the 5' ends of genes are reorganized in five patterns as a result of Brm knockdown. Intriguingly, the concomitant changes in the genes adjacent to the Brahma-dependent remodeling regions have important roles in development and morphogenesis. Further analyses reveal abundance of AT-rich motifs for transcription factors in the remodeling regions.

Computational discovery of feature patterns in nucleosomal DNA sequences

The identification of important factors that affect nucleosome formation is critical to clarify nucleosome-forming mechanisms and the role of the nucleosome in gene regulation. Various features reported in the literature led to our hypothesis that multiple features can together contribute to nucleosome formation. Therefore, we compiled 779 features and developed a pattern discovery and scoring algorithm FFN (Finding Features for Nucleosomes) to identify feature patterns that are differentially enriched in nucleosome-forming sequences and nucleosome-depletion sequences. Applying FFN to genome-wide nucleosome occupancy data in yeast and human, we identified statistically significant feature patterns that may influence nucleosome formation, many of which are common to the two species. We found that both sequence and structural features are important in nucleosome occupancy prediction. We discovered that, even for the same feature combinations, variations in feature values may lead to differences in predictive power. We demonstrated that the identified feature patterns could be used to assist nucleosomal sequence prediction.

RSC-dependent constructive and destructive interference between opposing arrays of phased nucleosomes in yeast

RSC and SWI/SNF are related ATP-dependent chromatin remodeling machines that move nucleosomes, regulating access to DNA. We addressed their roles in nucleosome phasing relative to transcription start sites in yeast. SWI/SNF has no effect on phasing at the global level. In contrast, RSC depletion results in global nucleosome repositioning: Both upstream and downstream nucleosomal arrays shift toward the nucleosome-depleted region (NDR), with no change in spacing, resulting in a narrower and partly filled NDR. The global picture of RSC-depleted chromatin represents the average of a range of chromatin structures, with most genes showing a shift of the +1 or the -1 nucleosome into the NDR. Using RSC ChIP data reported by others, we show that RSC occupancy is highest on the coding regions of heavily transcribed genes, though not at their NDRs. We propose that RSC has a role in restoring chromatin structure after transcription. Analysis of gene pairs in different orientations demonstrates that phasing patterns reflect competition between phasing signals emanating from neighboring NDRs. These signals may be in phase, resulting in constructive interference and a regular array, or out of phase, resulting in destructive interference and fuzzy positioning. We propose a modified barrier model, in which a stable complex located at the NDR acts as a bidirectional phasing barrier. In RSC-depleted cells, this barrier has a smaller footprint, resulting in narrower NDRs. Thus, RSC plays a critical role in organizing yeast chromatin.

The intriguing complexities of mammalian gene regulation: how to link enhancers to regulated genes. Are we there yet?

The information encoded in genomes supports the differentiation and function of the more than 200 unique cell types, which exist in various mammalian species. The major mechanism driving cellular differentiation and specification is differential gene expression regulation. Cis-acting enhancers and silencers appear to have key roles in regulating the expression of mammalian genes. However, these cis-acting elements are often located very far away from the regulated gene. Therefore, it is hard to find all of them and link them to the regulated gene. An intriguing and unresolved issue of the field is to identify all of the enhancers of a particular gene and link these short regulatory sequences to the genes they regulate and thus, reliably identify gene regulatory enhancer networks. Recent advances in molecular biological methods coupled with Next-Generation Sequencing (NGS) technologies have opened up new possibilities in this area of genomics. In this review we summarize the technological advances, bioinformatics challenges and the potential molecular mechanisms allowing the construction of enhancer networks operating in specific cell types and/or activated by various signals.

Terminate and make a loop: regulation of transcriptional directionality

•

Transcriptional directionality is controlled by premature transcription termination.

•

Transcriptional directionality is enforced by gene looping.

•

mRNA-specific termination signals and factors are required for gene looping.

Bidirectional promoters are a common feature of many eukaryotic organisms from yeast to humans. RNA Polymerase II that is recruited to this type of promoter can start transcribing in either direction using alternative DNA strands as the template. Such promiscuous transcription can lead to the synthesis of unwanted transcripts that may have negative effects on gene expression. Recent studies have identified transcription termination and gene looping as critical players in the enforcement of promoter directionality. Interestingly, both mechanisms share key components. Here, we focus on recent findings relating to the transcriptional output of bidirectional promoters.

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.

Histone hypoacetylation-activated genes are repressed by acetyl-CoA- and chromatin-mediated mechanism

Transcriptional activation is typically associated with increased acetylation of promoter histones. However, this paradigm does not apply to transcriptional activation of all genes. In this study we have characterized a group of genes that are repressed by histone acetylation. These histone hypoacetylation-activated genes (HHAAG) are normally repressed during exponential growth, when the cellular level of acetyl-CoA is high and global histone acetylation is also high. The HHAAG are induced during diauxic shift, when the levels of acetyl-CoA and global histone acetylation decrease. The histone hypoacetylation-induced activation of HHAAG is independent of Msn2/Msn4. The repression of HSP12, one of the HHAAG, is associated with well-defined nucleosomal structure in the promoter region, while histone hypoacetylation-induced activation correlates with delocalization of positioned nucleosomes or with reduced nucleosome occupancy. Correspondingly, unlike the majority of yeast genes, HHAAG are transcriptionally upregulated when expression of histone genes is reduced. Taken together, these results suggest a model in which histone acetylation is required for proper positioning of promoter nucleosomes and repression of HHAAG.

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.

A step-by-step protocol for formaldehyde-assisted isolation of regulatory elements from Arabidopsis thaliana

The control of gene expression by transcriptional regulators and other types of functionally relevant DNA transactions such as chromatin remodeling and replication underlie a vast spectrum of biological processes in all organisms. DNA transactions require the controlled interaction of proteins with DNA sequence motifs which are often located in nucleosome-depleted regions (NDRs) of the chromatin. Formaldehyde-assisted isolation of regulatory elements (FAIRE) has been established as an easy-to-implement method for the isolation of NDRs from a number of eukaryotic organisms, and it has been successfully employed for the discovery of new regulatory segments in genomic DNA from, for example, yeast, Drosophila, and humans. Until today, however, FAIRE has only rarely been employed in plant research and currently no detailed FAIRE protocol for plants has been published. Here, we provide a step-by-step FAIRE protocol for NDR discovery in Arabidopsis thaliana. We demonstrate that NDRs isolated from plant chromatin are readily amenable to quantitative polymerase chain reaction and next-generation sequencing. Only minor modification of the FAIRE protocol will be needed to adapt it to other plants, thus facilitating the global inventory of regulatory regions across species.

Genome-wide nucleosome positioning is orchestrated by genomic regions associated with DNase I hypersensitivity in rice

Nucleosome positioning dictates the DNA accessibility for regulatory proteins, and thus is critical for gene expression and regulation. It has been well documented that only a subset of nucleosomes are reproducibly positioned in eukaryotic genomes. The most prominent example of phased nucleosomes is the context of genes, where phased nucleosomes flank the transcriptional starts sites (TSSs). It is unclear, however, what factors determine nucleosome positioning in regions that are not close to genes. We mapped both nucleosome positioning and DNase I hypersensitive site (DHS) datasets across the rice genome. We discovered that DHSs located in a variety of contexts, both genic and intergenic, were flanked by strongly phased nucleosome arrays. Phased nucleosomes were also found to flank DHSs in the human genome. Our results suggest the barrier model may represent a general feature of nucleosome organization in eukaryote genomes. Specifically, regions bound with regulatory proteins, including intergenic regions, can serve as barriers that organize phased nucleosome arrays on both sides. Our results also suggest that rice DHSs often span a single, phased nucleosome, similar to the H2A.Z-containing nucleosomes observed in DHSs in the human genome.

The fundamental unit of chromatin is the nucleosome, which consists of 147 bp of DNA wrapped around a histone octamer containing four core histones (H3, H4, H2A, and H2B). Nucleosome positioning in the genome affects the DNA accessibility for regulatory proteins, and thus is critical for gene expression and regulation. Genomic regions associated with regulatory proteins are associated with a pronounced sensitivity to DNase I digestion, and are thus called DNase I hypersensitive sites (DHSs). It is well known that only a subset of nucleosomes are reproducibly positioned in eukaryotic genomes. However, it is less clear what factors determine genome-wide nucleosome positioning, especially in intergenic regions. We mapped both nucleosome positioning and DHS datasets across the rice genome. We discovered that DHSs located in a variety of contexts, both genic and intergenic, were flanked by strongly phased nucleosome arrays. We confirmed the same association of DHSs with phased nucleosomes in the human genome. We conclude that genomic loci associated with a diverse set of regulatory proteins are major determinants of nucleosome phasing, and this is true in both genic and intergenic regions.

DNA-encoded nucleosome occupancy is associated with transcription levels in the human malaria parasite Plasmodium falciparum

Background

In eukaryotic organisms, packaging of DNA into nucleosomes controls gene expression by regulating access of the promoter to transcription factors. The human malaria parasite Plasmodium falciparum encodes relatively few transcription factors, while extensive nucleosome remodeling occurs during its replicative cycle in red blood cells. These observations point towards an important role of the nucleosome landscape in regulating gene expression. However, the relation between nucleosome positioning and transcriptional activity has thus far not been explored in detail in the parasite.

Results

Here, we analyzed nucleosome positioning in the asexual and sexual stages of the parasite’s erythrocytic cycle using chromatin immunoprecipitation of MNase-digested chromatin, followed by next-generation sequencing. We observed a relatively open chromatin structure at the trophozoite and gametocyte stages, consistent with high levels of transcriptional activity in these stages. Nucleosome occupancy of genes and promoter regions were subsequently compared to steady-state mRNA expression levels. Transcript abundance showed a strong inverse correlation with nucleosome occupancy levels in promoter regions. In addition, AT-repeat sequences were strongly unfavorable for nucleosome binding in P. falciparum, and were overrepresented in promoters of highly expressed genes.

Conclusions

The connection between chromatin structure and gene expression in P. falciparum shares similarities with other eukaryotes. However, the remarkable nucleosome dynamics during the erythrocytic stages and the absence of a large variety of transcription factors may indicate that nucleosome binding and remodeling are critical regulators of transcript levels. Moreover, the strong dependency between chromatin structure and DNA sequence suggests that the P. falciparum genome may have been shaped by nucleosome binding preferences. Nucleosome remodeling mechanisms in this deadly parasite could thus provide potent novel anti-malarial targets.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-347) contains supplementary material, which is available to authorized users.

A genome-wide analysis of genetic diversity in Trypanosoma cruzi intergenic regions

BACKGROUND

Trypanosoma cruzi is the causal agent of Chagas Disease. Recently, the genomes of representative strains from two major evolutionary lineages were sequenced, allowing the construction of a detailed genetic diversity map for this important parasite. However this map is focused on coding regions of the genome, leaving a vast space of regulatory regions uncharacterized in terms of their evolutionary conservation and/or divergence.

METHODOLOGY

Using data from the hybrid CL Brener and Sylvio X10 genomes (from the TcVI and TcI Discrete Typing Units, respectively), we identified intergenic regions that share a common evolutionary ancestry, and are present in both CL Brener haplotypes (TcII-like and TcIII-like) and in the TcI genome; as well as intergenic regions that were conserved in only two of the three genomes/haplotypes analyzed. The genetic diversity in these regions was characterized in terms of the accumulation of indels and nucleotide changes.

PRINCIPAL FINDINGS

Based on this analysis we have identified i) a core of highly conserved intergenic regions, which remained essentially unchanged in independently evolving lineages; ii) intergenic regions that show high diversity in spite of still retaining their corresponding upstream and downstream coding sequences; iii) a number of defined sequence motifs that are shared by a number of unrelated intergenic regions. A fraction of indels explains the diversification of some intergenic regions by the expansion/contraction of microsatellite-like repeats.

Applications of alignment-free methods in epigenomics

Epigenetic mechanisms play an important role in the regulation of cell type-specific gene activities, yet how epigenetic patterns are established and maintained remains poorly understood. Recent studies have supported a role of DNA sequences in recruitment of epigenetic regulators. Alignment-free methods have been applied to identify distinct sequence features that are associated with epigenetic patterns and to predict epigenomic profiles. Here, we review recent advances in such applications, including the methods to map DNA sequence to feature space, sequence comparison and prediction models. Computational studies using these methods have provided important insights into the epigenetic regulatory mechanisms.

Training-free atomistic prediction of nucleosome occupancy

Nucleosomes alter gene expression by preventing transcription factors from occupying binding sites along DNA. DNA methylation can affect nucleosome positioning and so alter gene expression epigenetically (without changing DNA sequence). Conventional methods to predict nucleosome occupancy are trained on observed DNA sequence patterns or known DNA oligonucleotide structures. They are statistical and lack the physics needed to predict subtle epigenetic changes due to DNA methylation. The training-free method presented here uses physical principles and state-of-the-art all-atom force fields to predict both nucleosome occupancy along genomic sequences as well as binding to known positioning sequences. Our method calculates the energy of both nucleosomal and linear DNA of the given sequence. Based on the DNA deformation energy, we accurately predict the in vitro occupancy profile observed experimentally for a 20,000-bp genomic region as well as the experimental locations of nucleosomes along 13 well-established positioning sequence elements. DNA with all C bases methylated at the 5 position shows less variation of nucleosome binding: Strong binding is weakened and weak binding is strengthened compared with normal DNA. Methylation also alters the preference of nucleosomes for some positioning sequences but not others.

Nucleosomes are stably evicted from enhancers but not promoters upon induction of certain pro-inflammatory genes in mouse macrophages

Chromatin is thought to act as a barrier for binding of cis-regulatory transcription factors (TFs) to their sites on DNA and recruitment of the transcriptional machinery. Here we have analyzed changes in nucleosome occupancy at the enhancers as well as at the promoters of three pro-inflammatory genes when they are induced by bacterial lipopolysaccharides (LPS) in primary mouse macrophages. We find that nucleosomes are removed from the distal enhancers of IL12B and IL1A, as well as from the distal and proximal enhancers of IFNB1, and that clearance of enhancers correlates with binding of various cis-regulatory TFs. We further show that for IFNB1 the degree of nucleosome removal correlates well with the level of induction of the gene under different conditions. Surprisingly, we find that nucleosome occupancy at the promoters of IL12B and IL1A does not change significantly when the genes are induced, and that a considerably fraction of the cells is occupied by nucleosomes at any given time. We hypothesize that competing nucleosomes at the promoters of IL12B and IL1A may play a role in limiting the size of transcriptional bursts in individual cells, which may be important for controlling cytokine production in a population of immune cells.

Distinct role of core promoter architecture in regulation of light-mediated responses in plant genes

In the present study, we selected four distinct classes of light-regulated promoters. The light-regulated promoters can be distinctly grouped into either TATA-box-containing or TATA-less (initiator-containing) promoters. Further, using either native promoters or their swapped versions of core promoter elements, we established that TATA-box and Inr (Initiator) elements have distinct mechanisms which are involved in light-mediated regulation, and these elements are not swappable. We identified that mutations in either functional TATA-box or Inr elements lead to the formation of nucleosomal structure. The nucleotide diversity in either the TATA-box or Inr element in Arabidopsis ecotypes proposes that the nucleotide variation in core promoters can alter the gene expression. We show that motif overrepresentation in light-activated promoters encompasses different specific regulatory motifs present downstream of TSS (transcription start site), and this might serve as a key factor in regulating light promoters which are parallel with these elements. Finally, we conclude that the TATA-box or Inr element does not act in isolation, but our results clearly suggests the probable involvement of other distinct core promoter elements in concurrence with the TATA-box or Inr element to impart selectivity to light-mediated transcription.

Dynamic remodeling of histone modifications in response to osmotic stress in Saccharomyces cerevisiae

Background

Specific histone modifications play important roles in chromatin functions; i.e., activation or repression of gene transcription. This participation must occur as a dynamic process. Nevertheless, most of the histone modification maps reported to date provide only static pictures that link certain modifications with active or silenced states. This study, however, focuses on the global histone modification variation that occurs in response to the transcriptional reprogramming produced by a physiological perturbation in yeast.

Results

We did a genome-wide chromatin immunoprecipitation analysis for eight specific histone modifications before and after saline stress. The most striking change was rapid acetylation loss in lysines 9 and 14 of H3 and in lysine 8 of H4, associated with gene repression. The genes activated by saline stress increased the acetylation levels at these same sites, but this acetylation process was quantitatively minor if compared to that of the deacetylation of repressed genes. The changes in the tri-methylation of lysines 4, 36 and 79 of H3 and the di-methylation of lysine 79 of H3 were slighter than those of acetylation. Furthermore, we produced new genome-wide maps for seven histone modifications, and we analyzed, for the first time in S. cerevisiae, the genome-wide profile of acetylation of lysine 8 of H4.

Conclusions

This research reveals that the short-term changes observed in the post-stress methylation of histones are much more moderate than those of acetylation, and that the dynamics of the acetylation state of histones during activation or repression of transcription is a much quicker process than methylation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-247) contains supplementary material, which is available to authorized users.

High-density nucleosome occupancy map of human chromosome 9p21-22 reveals chromatin organization of the type I interferon gene cluster

Genome-wide investigations have dramatically increased our understanding of nucleosome positioning and the role of chromatin in gene regulation, yet some genomic regions have been poorly represented in human nucleosome maps. One such region is represented by human chromosome 9p21-22, which contains the type I interferon gene cluster that includes 16 interferon alpha genes and the single interferon beta, interferon epsilon, and interferon omega genes. A high-density nucleosome mapping strategy was used to generate locus-wide maps of the nucleosome organization of this biomedically important locus at a steady state and during a time course of infection with Sendai virus, an inducer of interferon gene expression. Detailed statistical and computational analysis illustrates that nucleosomes in this locus exhibit preferences for particular dinucleotide and oligomer DNA sequence motifs in vivo, which are similar to those reported for lower eukaryotic nucleosome-DNA interactions. These data were used to visualize the region's chromatin architecture and reveal features that are common to the organization of all the type I interferon genes, indicating a common nucleosome-mediated gene regulatory paradigm. Additionally, this study clarifies aspects of the dynamic changes that occur with the nucleosome occupying the transcriptional start site of the interferon beta gene after virus infection.

Mapping yeast transcriptional networks

The term “transcriptional network” refers to the mechanism(s) that underlies coordinated expression of genes, typically involving transcription factors (TFs) binding to the promoters of multiple genes, and individual genes controlled by multiple TFs. A multitude of studies in the last two decades have aimed to map and characterize transcriptional networks in the yeast Saccharomyces cerevisiae. We review the methodologies and accomplishments of these studies, as well as challenges we now face. For most yeast TFs, data have been collected on their sequence preferences, in vivo promoter occupancy, and gene expression profiles in deletion mutants. These systematic studies have led to the identification of new regulators of numerous cellular functions and shed light on the overall organization of yeast gene regulation. However, many yeast TFs appear to be inactive under standard laboratory growth conditions, and many of the available data were collected using techniques that have since been improved. Perhaps as a consequence, comprehensive and accurate mapping among TF sequence preferences, promoter binding, and gene expression remains an open challenge. We propose that the time is ripe for renewed systematic efforts toward a complete mapping of yeast transcriptional regulatory mechanisms.

The Mi-2 homolog Mit1 actively positions nucleosomes within heterochromatin to suppress transcription

Mit1 is the putative chromatin remodeling subunit of the fission yeast Snf2/histone deacetylase (HDAC) repressor complex (SHREC) and is known to repress transcription at regions of heterochromatin. However, how Mit1 modifies chromatin to silence transcription is largely unknown. Here we report that Mit1 mobilizes histone octamers in vitro and requires ATP hydrolysis and conserved chromatin tethering domains, including a previously unrecognized chromodomain, to remodel nucleosomes and silence transcription. Loss of Mit1 remodeling activity results in nucleosome depletion at specific DNA sequences that display low intrinsic affinity for the histone octamer, but its contribution to antagonizing RNA polymerase II (Pol II) access and transcription is not restricted to these sites. Genetic epistasis analyses demonstrate that SHREC subunits and the transcription-coupled Set2 histone methyltransferase, which is involved in suppression of cryptic transcription at actively transcribed regions, cooperate to silence heterochromatic transcripts. In addition, we have demonstrated that Mit1's remodeling activity contributes to SHREC function independently of Clr3's histone deacetylase activity on histone H3 K14. We propose that Mit1 is a chromatin remodeling factor that cooperates with the Clr3 histone deacetylase of SHREC and other chromatin modifiers to stabilize heterochromatin structure and to prevent access to the transcriptional machinery.

Transcriptional enhancers: from properties to genome-wide predictions

Cellular development, morphology and function are governed by precise patterns of gene expression. These are established by the coordinated action of genomic regulatory elements known as enhancers or cis-regulatory modules. More than 30 years after the initial discovery of enhancers, many of their properties have been elucidated; however, despite major efforts, we only have an incomplete picture of enhancers in animal genomes. In this Review, we discuss how properties of enhancer sequences and chromatin are used to predict enhancers in genome-wide studies. We also cover recently developed high-throughput methods that allow the direct testing and identification of enhancers on the basis of their activity. Finally, we discuss recent technological advances and current challenges in the field of regulatory genomics.

GC-rich DNA elements enable replication origin activity in the methylotrophic yeast Pichia pastoris

The well-studied DNA replication origins of the model budding and fission yeasts are A/T-rich elements. However, unlike their yeast counterparts, both plant and metazoan origins are G/C-rich and are associated with transcription start sites. Here we show that an industrially important methylotrophic budding yeast, Pichia pastoris, simultaneously employs at least two types of replication origins--a G/C-rich type associated with transcription start sites and an A/T-rich type more reminiscent of typical budding and fission yeast origins. We used a suite of massively parallel sequencing tools to map and dissect P. pastoris origins comprehensively, to measure their replication dynamics, and to assay the global positioning of nucleosomes across the genome. Our results suggest that some functional overlap exists between promoter sequences and G/C-rich replication origins in P. pastoris and imply an evolutionary bifurcation of the modes of replication initiation.

A dynamic interplay of nucleosome and Msn2 binding regulates kinetics of gene activation and repression following stress

The transcription factor Msn2 mediates a significant proportion of the environmental stress response, in which a common cohort of genes changes expression in a stereotypic fashion upon exposure to any of a wide variety of stresses. We have applied genome-wide chromatin immunoprecipitation and nucleosome profiling to determine where Msn2 binds under stressful conditions and how that binding affects, and is affected by, nucleosome positioning. We concurrently determined the effect of Msn2 activity on gene expression following stress and demonstrated that Msn2 stimulates both activation and repression. We found that some genes responded to both intermittent and continuous Msn2 nuclear occupancy while others responded only to continuous occupancy. Finally, these studies document a dynamic interplay between nucleosomes and Msn2 such that nucleosomes can restrict access of Msn2 to its canonical binding sites while Msn2 can promote reposition, expulsion and recruitment of nucleosomes to alter gene expression. This interplay may allow the cell to discriminate between different types of stress signaling.