Nucleosome depletion at yeast terminators is not intrinsic and can occur by a transcriptional mechanism linked to 3'-end formation

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate

The RNA polymerase II (Pol II) elongation rate influences poly(A) site selection, with slow and fast Pol II derivatives causing upstream and downstream shifts, respectively, in poly(A) site utilization. In yeast, depletion of either of the histone chaperones FACT or Spt6 causes an upstream shift of poly(A) site use that strongly resembles the poly(A) profiles of slow Pol II mutant strains. Like slow Pol II mutant strains, FACT- and Spt6-depleted cells exhibit Pol II processivity defects, indicating that both Spt6 and FACT stimulate the Pol II elongation rate. Poly(A) profiles of some genes show atypical downstream shifts; this subset of genes overlaps well for FACT- or Spt6-depleted strains but is different from the atypical genes in Pol II speed mutant strains. In contrast, depletion of histone H3 or H4 causes a downstream shift of poly(A) sites for most genes, indicating that nucleosomes inhibit the Pol II elongation rate in vivo. Thus, chromatin-based control of the Pol II elongation rate is a potential mechanism, distinct from direct effects on the cleavage/polyadenylation machinery, to regulate alternative polyadenylation in response to genetic or environmental changes.

How is polyadenylation restricted to 3'-untranslated regions?

Polyadenylation occurs at numerous sites within 3'-untranslated regions (3'-UTRs) but rarely within coding regions. How does Pol II travel through long coding regions without generating poly(A) sites, yet then permits promiscuous polyadenylation once it reaches the 3'-UTR? The cleavage/polyadenylation (CpA) machinery preferentially associates with 3'-UTRs, but it is unknown how its recruitment is restricted to 3'-UTRs during Pol II elongation. Unlike coding regions, 3'-UTRs have long AT-rich stretches of DNA that may be important for restricting polyadenylation to 3'-UTRs. Recognition of the 3'-UTR could occur at the DNA (AT-rich), RNA (AU-rich), or RNA:DNA hybrid (rU:dA- and/or rA:dT-rich) level. Based on the nucleic acid critical for 3'-UTR recognition, there are three classes of models, not mutually exclusive, for how the CpA machinery is selectively recruited to 3'-UTRs, thereby restricting where polyadenylation occurs: (1) RNA-based models suggest that the CpA complex directly (or indirectly through one or more intermediary proteins) binds long AU-rich stretches that are exposed after Pol II passes through these regions. (2) DNA-based models suggest that the AT-rich sequence affects nucleosome depletion or the elongating Pol II machinery, resulting in dissociation of some elongation factors and subsequent recruitment of the CpA machinery. (3) RNA:DNA hybrid models suggest that preferential destabilization of the Pol II elongation complex at rU:dA- and/or rA:dT-rich duplexes bridging the nucleotide addition and RNA exit sites permits preferential association of the CpA machinery with 3'-UTRs. Experiments to provide evidence for one or more of these models are suggested.

Condition-specific 3' mRNA isoform half-lives and stability elements in yeast

Alternative polyadenylation generates numerous 3' mRNA isoforms that can differ in their stability, structure, and function. These isoforms can be used to map mRNA stabilizing and destabilizing elements within 3' untranslated regions (3'UTRs). Here, we examine how environmental conditions affect 3' mRNA isoform turnover and structure in yeast cells on a transcriptome scale. Isoform stability broadly increases when cells grow more slowly, with relative half-lives of most isoforms being well correlated across multiple conditions. Surprisingly, dimethyl sulfate probing reveals that individual 3' isoforms have similar structures across different conditions, in contrast to the extensive structural differences that can exist between closely related isoforms in an individual condition. Unexpectedly, most mRNA stabilizing and destabilizing elements function only in a single growth condition. The genes associated with some classes of condition-specific stability elements are enriched for different functional categories, suggesting that regulated mRNA stability might contribute to adaptation to different growth environments. Condition-specific stability elements do not result in corresponding condition-specific changes in steady-state mRNA isoform levels. This observation is consistent with a compensatory mechanism between polyadenylation and stability, and it suggests that condition-specific mRNA stability elements might largely reflect condition-specific regulation of mRNA 3' end formation.

Sumoylation is Largely Dispensable for Normal Growth but Facilitates Heat Tolerance in Yeast

Numerous proteins are sumoylated in normally growing yeast and SUMO conjugation levels rise upon exposure to several stress conditions. We observe high levels of sumoylation also during early exponential growth and when nutrient-rich medium is used. However, we find that reduced sumoylation (∼75% less than normal) is remarkably well-tolerated, with no apparent growth defects under nonstress conditions or under osmotic, oxidative, or ethanol stresses. In contrast, strains with reduced activity of Ubc9, the sole SUMO conjugase, are temperature-sensitive, implicating sumoylation in the heat stress response, specifically. Aligned with this, a mild heat shock triggers increased sumoylation which requires functional levels of Ubc9, but likely also depends on decreased desumoylation, since heat shock reduces protein levels of Ulp1, the major SUMO protease. Furthermore, we find that a ubc9 mutant strain with only ∼5% of normal sumoylation levels shows a modest growth defect, has abnormal genomic distribution of RNA polymerase II (RNAPII), and displays a greatly expanded redistribution of RNAPII after heat shock. Together, our data implies that SUMO conjugations are largely dispensable under normal conditions, but a threshold level of Ubc9 activity is needed to maintain transcriptional control and to modulate the redistribution of RNAPII and promote survival when temperatures rise.

Involvement of the SAGA and TFIID coactivator complexes in transcriptional dysregulation caused by the separation of core and tail Mediator modules

Regulation of RNA polymerase II transcription requires the concerted efforts of several multisubunit coactivator complexes, which interact with the RNA polymerase II preinitiation complex to stimulate transcription. We previously showed that separation of the Mediator core from Mediator's tail module results in modest overactivation of genes annotated as highly dependent on TFIID for expression. However, it is unclear if other coactivators are involved in this phenomenon. Here, we show that the overactivation of certain genes by Mediator core/tail separation is blunted by disruption of the Spt-Ada-Gcn5-Acetyl transferase complex through the removal of its structural Spt20 subunit, though this downregulation does not appear to completely depend on reduced Spt-Ada-Gcn5-Acetyl transferase association with the genome. Consistent with the enrichment of TFIID-dependent genes among genes overactivated by Mediator core/tail separation, depletion of the essential TFIID subunit Taf13 suppressed the overactivation of these genes when Med16 was simultaneously removed. As with Spt-Ada-Gcn5-Acetyl transferase, this effect did not appear to be fully dependent on the reduced genomic association of TFIID. Given that the observed changes in gene expression could not be clearly linked to alterations in Spt-Ada-Gcn5-Acetyl transferase or TFIID occupancy, our data may suggest that the Mediator core/tail connection is important for the modulation of Spt-Ada-Gcn5-Acetyl transferase and/or TFIID conformation and/or function at target genes.

GC content, but not nucleosome positioning, directly contributes to intron splicing efficiency in Paramecium

Eukaryotic genes are interrupted by introns that must be accurately spliced from mRNA precursors. With an average length of 25 nt, the more than 90,000 introns of Paramecium tetraurelia stand among the shortest introns reported in eukaryotes. The mechanisms specifying the correct recognition of these tiny introns remain poorly understood. Splicing can occur cotranscriptionally, and it has been proposed that chromatin structure might influence splice site recognition. To investigate the roles of nucleosome positioning in intron recognition, we determined the nucleosome occupancy along the P. tetraurelia genome. We show that P. tetraurelia displays a regular nucleosome array with a nucleosome repeat length of ∼151 bp, among the smallest periodicities reported. Our analysis has revealed that introns are frequently associated with inter-nucleosomal DNA, pointing to an evolutionary constraint favoring introns at the AT-rich nucleosome edge sequences. Using accurate splicing efficiency data from cells depleted for nonsense-mediated decay effectors, we show that introns located at the edge of nucleosomes display higher splicing efficiency than those at the center. However, multiple regression analysis indicates that the low GC content of introns, rather than nucleosome positioning, is associated with high splicing efficiency. Our data reveal a complex link between GC content, nucleosome positioning, and intron evolution in Paramecium.

Landscape of transcription termination in Arabidopsis revealed by single-molecule nascent RNA sequencing

BACKGROUND

The dynamic process of transcription termination produces transient RNA intermediates that are difficult to distinguish from each other via short-read sequencing methods.

RESULTS

Here, we use single-molecule nascent RNA sequencing to characterize the various forms of transient RNAs during termination at genome-wide scale in wildtype Arabidopsis and in atxrn3, fpa, and met1 mutants. Our data reveal a wide range of termination windows among genes, ranging from ~ 50 nt to over 1000 nt. We also observe efficient termination before downstream tRNA genes, suggesting that chromatin structure around the promoter region of tRNA genes may block pol II elongation. 5' Cleaved readthrough transcription in atxrn3 with delayed termination can run into downstream genes to produce normally spliced and polyadenylated mRNAs in the absence of their own transcription initiation. Consistent with previous reports, we also observe long chimeric transcripts with cryptic splicing in fpa mutant; but loss of CG DNA methylation has no obvious impact on termination in the met1 mutant.

CONCLUSIONS

Our method is applicable to establish a comprehensive termination landscape in a broad range of species.

Connection of core and tail Mediator modules restrains transcription from TFIID-dependent promoters

The Mediator coactivator complex is divided into four modules: head, middle, tail, and kinase. Deletion of the architectural subunit Med16 separates core Mediator (cMed), comprising the head, middle, and scaffold (Med14), from the tail. However, the direct global effects of tail/cMed disconnection are unclear. We find that rapid depletion of Med16 downregulates genes that require the SAGA complex for full expression, consistent with their reported tail dependence, but also moderately overactivates TFIID-dependent genes in a manner partly dependent on the separated tail, which remains associated with upstream activating sequences. Suppression of TBP dynamics via removal of the Mot1 ATPase partially restores normal transcriptional activity to Med16-depleted cells, suggesting that cMed/tail separation results in an imbalance in the levels of PIC formation at SAGA-requiring and TFIID-dependent genes. We propose that the preferential regulation of SAGA-requiring genes by tailed Mediator helps maintain a proper balance of transcription between these genes and those more dependent on TFIID.

Composed of over two dozen subunits, the Mediator complex plays several roles in RNA polymerase II (RNAPII) transcription in eukaryotes. In yeast, deletion of Med16, which splits Mediator into two stable subcomplexes, both increases and decreases transcript levels, suggesting that Med16 might play a repressive role. However, the direct effects of Med16 removal on RNAPII transcription have not been assessed, owing to the use of deletion mutants and measurement of steady-state RNA levels in prior studies. Here, using a combination of inducible protein depletion and analysis of nascent RNA, we find that Med16 removal 1) downregulates a small group of genes reported to be highly dependent on the SAGA complex and 2) upregulates a larger set of genes reported to be more dependent on the TFIID complex in a manner dependent on another component of Mediator. We find that artificially altering the balance of transcription pre-initiation complex (PIC) formation toward SAGA-requiring promoters and away from TFIID-dependent promoters partially restores normal transcription, indicating a contribution of altered PIC formation to the transcriptional alterations observed with Med16 loss. Taken together, our results indicate that the structural integrity of Mediator is important for maintaining balanced transcription between different gene classes.

The Active Mechanism of Nucleosome Depletion by Poly(dA:dT) Tracts In Vivo

Poly(dA:dT) tracts cause nucleosome depletion in many species, e.g., at promoters and replication origins. Their intrinsic biophysical sequence properties make them stiff and unfavorable for nucleosome assembly, as probed by in vitro nucleosome reconstitution. The mere correlation between nucleosome depletion over poly(dA:dT) tracts in in vitro reconstituted and in in vivo chromatin inspired an intrinsic nucleosome exclusion mechanism in vivo that is based only on DNA and histone properties. However, we compile here published and new evidence that this correlation does not reflect mechanistic causation. (1) Nucleosome depletion over poly(dA:dT) in vivo is not universal, e.g., very weak in S. pombe. (2) The energy penalty for incorporating poly(dA:dT) tracts into nucleosomes is modest (<10%) relative to ATP hydrolysis energy abundantly invested by chromatin remodelers. (3) Nucleosome depletion over poly(dA:dT) is much stronger in vivo than in vitro if monitored without MNase and (4) actively maintained in vivo. (5) S. cerevisiae promoters evolved a strand-biased poly(dA) versus poly(dT) distribution. (6) Nucleosome depletion over poly(dA) is directional in vivo. (7) The ATP dependent chromatin remodeler RSC preferentially and directionally displaces nucleosomes towards 5′ of poly(dA). Especially distribution strand bias and displacement directionality would not be expected for an intrinsic mechanism. Together, this argues for an in vivo mechanism where active and species-specific read out of intrinsic sequence properties, e.g., by remodelers, shapes nucleosome organization.

Regulation of alternative polyadenylation in the yeast Saccharomyces cerevisiae by histone H3K4 and H3K36 methyltransferases

Adjusting DNA structure via epigenetic modifications, and altering polyadenylation (pA) sites at which precursor mRNA is cleaved and polyadenylated, allows cells to quickly respond to environmental stress. Since polyadenylation occurs co-transcriptionally, and specific patterns of nucleosome positioning and chromatin modifications correlate with pA site usage, epigenetic factors potentially affect alternative polyadenylation (APA). We report that the histone H3K4 methyltransferase Set1, and the histone H3K36 methyltransferase Set2, control choice of pA site in Saccharomyces cerevisiae, a powerful model for studying evolutionarily conserved eukaryotic processes. Deletion of SET1 or SET2 causes an increase in serine-2 phosphorylation within the C-terminal domain of RNA polymerase II (RNAP II) and in the recruitment of the cleavage/polyadenylation complex, both of which could cause the observed switch in pA site usage. Chemical inhibition of TOR signaling, which causes nutritional stress, results in Set1- and Set2-dependent APA. In addition, Set1 and Set2 decrease efficiency of using single pA sites, and control nucleosome occupancy around pA sites. Overall, our study suggests that the methyltransferases Set1 and Set2 regulate APA induced by nutritional stress, affect the RNAP II C-terminal domain phosphorylation at Ser2, and control recruitment of the 3′ end processing machinery to the vicinity of pA sites.

Active transcription and Orc1 drive chromatin association of the AAA+ ATPase Pch2 during meiotic G2/prophase

Pch2 is an AAA+ protein that controls DNA break formation, recombination and checkpoint signaling during meiotic G2/prophase. Chromosomal association of Pch2 is linked to these processes, and several factors influence the association of Pch2 to euchromatin and the specialized chromatin of the ribosomal (r)DNA array of budding yeast. Here, we describe a comprehensive mapping of Pch2 localization across the budding yeast genome during meiotic G2/prophase. Within non-rDNA chromatin, Pch2 associates with a subset of actively RNA Polymerase II (RNAPII)-dependent transcribed genes. Chromatin immunoprecipitation (ChIP)- and microscopy-based analysis reveals that active transcription is required for chromosomal recruitment of Pch2. Similar to what was previously established for association of Pch2 with rDNA chromatin, we find that Orc1, a component of the Origin Recognition Complex (ORC), is required for the association of Pch2 to these euchromatic, transcribed regions, revealing a broad connection between chromosomal association of Pch2 and Orc1/ORC function. Ectopic mitotic expression is insufficient to drive recruitment of Pch2, despite the presence of active transcription and Orc1/ORC in mitotic cells. This suggests meiosis-specific ‘licensing’ of Pch2 recruitment to sites of transcription, and accordingly, we find that the synaptonemal complex (SC) component Zip1 is required for the recruitment of Pch2 to transcription-associated binding regions. Interestingly, Pch2 binding patterns are distinct from meiotic axis enrichment sites (as defined by Red1, Hop1, and Rec8). Inactivating RNAPII-dependent transcription/Orc1 does not lead to effects on the chromosomal abundance of Hop1, a known chromosomal client of Pch2, suggesting a complex relationship between SC formation, Pch2 recruitment and Hop1 chromosomal association. We thus report characteristics and dependencies for Pch2 recruitment to meiotic chromosomes, and reveal an unexpected link between Pch2, SC formation, chromatin and active transcription.

Meiosis is a specialized cellular division program that is required to produce haploid reproductive cells, also known as gametes. To allow meiosis to occur faithfully, several processes centred around DNA breakage and recombination are needed. Pch2, an AAA+ ATPase enzyme is important to coordinate several of these processes. Here, we analyze the genome-wide association of Pch2 to budding yeast meiotic chromosomes. Our results show that Pch2 is recruited to a subset of actively transcribed genes, and we find that active RNAPII transcription contributes to Pch2 chromosomal association. In addition, we reveal a general contribution of Orc1, a subunit of the ORC assembly, to Pch2 chromosomal recruitment. These findings thus reveal a connection between Pch2, Orc1 and RNAPII activity during meiosis.

Requirements for RNA polymerase II preinitiation complex formation in vivo

Transcription by RNA polymerase II requires assembly of a preinitiation complex (PIC) composed of general transcription factors (GTFs) bound at the promoter. In vitro, some GTFs are essential for transcription, whereas others are not required under certain conditions. PICs are stable in the absence of nucleotide triphosphates, and subsets of GTFs can form partial PICs. By depleting individual GTFs in yeast cells, we show that all GTFs are essential for TBP binding and transcription, suggesting that partial PICs do not exist at appreciable levels in vivo. Depletion of FACT, a histone chaperone that travels with elongating Pol II, strongly reduces PIC formation and transcription. In contrast, TBP-associated factors (TAFs) contribute to transcription of most genes, but TAF-independent transcription occurs at substantial levels, preferentially at promoters containing TATA elements. PICs are absent in cells deprived of uracil, and presumably UTP, suggesting that transcriptionally inactive PICs are removed from promoters in vivo.

A strong structural correlation between short inverted repeat sequences and the polyadenylation signal in yeast and nucleosome exclusion by these inverted repeats

DNA sequences that read the same from 5′ to 3′ in either strand are called inverted repeat sequences or simply IRs. They are found throughout a wide variety of genomes, from prokaryotes to eukaryotes. Despite extensive research, their in vivo functions, if any, remain unclear. Using Saccharomyces cerevisiae, we performed genome-wide analyses for the distribution, occurrence frequency, sequence characteristics and relevance to chromatin structure, for the IRs that reportedly have a cruciform-forming potential. Here, we provide the first comprehensive map of these IRs in the S. cerevisiae genome. The statistically significant enrichment of the IRs was found in the close vicinity of the DNA positions corresponding to polyadenylation [poly(A)] sites and ~ 30 to ~ 60 bp downstream of start codon-coding sites (referred to as ‘start codons’). In the former, ApT- or TpA-rich IRs and A-tract- or T-tract-rich IRs are enriched, while in the latter, different IRs are enriched. Furthermore, we found a strong structural correlation between the former IRs and the poly(A) signal. In the chromatin formed on the gene end regions, the majority of the IRs causes low nucleosome occupancy. The IRs in the region ~ 30 to ~ 60 bp downstream of start codons are located in the + 1 nucleosomes. In contrast, fewer IRs are present in the adjacent region downstream of start codons. The current study suggests that the IRs play similar roles in Escherichia coli and S. cerevisiae to regulate or complete transcription at the RNA level.

A unified computational framework for modeling genome-wide nucleosome landscape

Nucleosomes form the fundamental building blocks of eukaryotic chromatin, and previous attempts to understand the principles governing their genome-wide distribution have spurred much interest and debate in biology. In particular, the precise role of DNA sequence in shaping local chromatin structure has been controversial. This paper rigorously quantifies the contribution of hitherto-debated sequence features-including G+C content, 10.5 bp periodicity, and poly(dA:dT) tracts-to three distinct aspects of genome-wide nucleosome landscape: occupancy, translational positioning and rotational positioning. Our computational framework simultaneously learns nucleosome number and nucleosome-positioning energy from genome-wide nucleosome maps. In contrast to other previous studies, our model can predict both in vitro and in vivo nucleosome maps in Saccharomyces cerevisiae. We find that although G+C content is the primary determinant of MNase-derived nucleosome occupancy, MNase digestion biases may substantially influence this GC dependence. By contrast, poly(dA:dT) tracts are seen to deter nucleosome formation, regardless of the experimental method used. We further show that the 10.5 bp nucleotide periodicity facilitates rotational but not translational positioning. Applying our method to in vivo nucleosome maps demonstrates that, for a subset of genes, the regularly-spaced nucleosome arrays observed around transcription start sites can be partially recapitulated by DNA sequence alone. Finally, in vivo nucleosome occupancy derived from MNase-seq experiments around transcription termination sites can be mostly explained by the genomic sequence. Implications of these results and potential extensions of the proposed computational framework are discussed.

DMS-Seq for In Vivo Genome-wide Mapping of Protein-DNA Interactions and Nucleosome Centers

Protein-DNA interactions provide the basis for chromatin structure and gene regulation. Comprehensive identification of protein-occupied sites is thus vital to an in-depth understanding of genome function. Dimethyl sulfate (DMS) is a chemical probe that has long been used to detect footprints of DNA-bound proteins in vitro and in vivo. Here, we describe a genomic footprinting method, dimethyl sulfate sequencing (DMS-seq), which exploits the cell-permeable nature of DMS to obviate the need for nuclear isolation. This feature makes DMS-seq simple in practice and removes the potential risk of protein re-localization during nuclear isolation. DMS-seq successfully detects transcription factors bound to cis-regulatory elements and non-canonical chromatin particles in nucleosome-free regions. Furthermore, an unexpected preference of DMS confers on DMS-seq a unique potential to directly detect nucleosome centers without using genetic manipulation. We expect that DMS-seq will serve as a characteristic method for genome-wide interrogation of in vivo protein-DNA interactions.

Chemically induced proximity in biology and medicine

Proximity, or the physical closeness of molecules, is a pervasive regulatory mechanism in biology. For example, most posttranslational modifications such as phosphorylation, methylation, and acetylation promote proximity of molecules to play deterministic roles in cellular processes. To understand the role of proximity in biologic mechanisms, chemical inducers of proximity (CIPs) were developed to synthetically model biologically regulated recruitment. Chemically induced proximity allows for precise temporal control of transcription, signaling cascades, chromatin regulation, protein folding, localization, and degradation, as well as a host of other biologic processes. A systematic analysis of CIPs in basic research, coupled with recent technological advances utilizing CRISPR, distinguishes roles of causality from coincidence and allows for mathematical modeling in synthetic biology. Recently, induced proximity has provided new avenues of gene therapy and emerging advances in cancer treatment.

Mechanism of Apoptosis Induction by Mycoplasmal Nuclease MGA_0676 in Chicken Embryo Fibroblasts

MGA_0676 has been characterized as a Mycoplasma gallisepticum nuclease that can induce apoptosis of chicken cells. However, the mechanism by which MGA_0676 induces apoptosis has remained unclear. In this study, we evaluated MGA_0676-induced apoptosis and internalization in immortalized chicken embryo fibroblasts (DF-1) and cancer cell lines. The internalization of MGA_0676 was proven through caveolin-mediated endocytosis by blocking the endocytosis with specific inhibitors or with siRNA. We identified the Thif domain of NEDD8-activating enzyme E1 regulatory subunit (NAE) in DF-1 as the target region interacting with the SNC domain of MGA_0676. The interaction between the Thif and SNC domains was observed co-located in the perinuclear and nuclear of DF-1. We found that the interaction between NAE and MGA_0676 increased the ability of apoptosis and accelerated the process of cullin neddylation in DF-1 cells, in turn activating NF-κB. This resulted in the observed aggregation of NF-κB in the nuclei of DF-1 cells. Moreover, the apoptosis induced by MGA_0676 decreased significantly when NF-κB was inhibited by siRNA or BAY 11-7082 or when NAE was silenced by siRNA. Overall, our results demonstrate that MGA_0676 is internalized through caveolin-mediated endocytosis, interacts with SNC-dependent Thif to accelerate the process of cullin neddylation and activates NF-κB in DF-1 cells, ultimately playing a key role in apoptosis in chicken cells. Our results indicate MGA_0676 constitutes a critical etiological virulence factor of the respiratory disease caused by M. gallisepticum. This study also opens a venue to investigate MGA_0676 as a potential candidate as pro-apoptotic drug in cancer studies.

Identification of Two Distinct Classes of the Human INO80 Complex Genome-Wide

Chromatin remodeling and histone modifying enzymes play a critical role in shaping the regulatory output of a cell. Although much is known about these classes of proteins, identifying the mechanisms by which they coordinate gene expression programs remains an exciting topic of investigation. One factor that may contribute to the targeting and activity of chromatin regulators is local chromatin landscape. We leveraged genomic approaches and publically-available datasets to characterize the chromatin landscape at targets of the human INO80 chromatin remodeling complex (INO80-C). Our data revealed two classes of INO80-C targets with distinct chromatin signatures. The predominant INO80-C class was enriched for open chromatin, H3K27ac, and representative subunits from each of the three INO80-C modules (RUVBL1, RUVBL2, MCRS1, YY1). We named this class Canonical INO80. Notably, we identified an unexpected class of INO80-C targets that contained only the INO80 ATPase and harbored a repressive chromatin signature characterized by inaccessible chromatin, H3K27me3, and the methyltransferase EZH2. We named this class Non-Canonical INO80 (NC-INO80). Biochemical approaches indicated that INO80-C and the H3K27 acetyltransferase P300 physically interact, suggesting INO80-C and P300 may jointly coordinate chromatin accessibility at Canonical INO80 sites. No interaction was detected between INO80-C and EZH2, indicating INO80-C and EZH2 may engage in a separate form of regulatory crosstalk at NC-INO80 targets. Our data indicate that INO80-C is more compositionally heterogenous at its genomic targets than anticipated. Moreover, our data suggest there is an important link between INO80-C and histone modifying enzymes that may have consequences in developmental and pathological contexts.

Precise genome-wide mapping of single nucleosomes and linkers in vivo

We developed a chemical cleavage method that releases single nucleosome dyad-containing fragments, allowing us to precisely map both single nucleosomes and linkers with high accuracy genome-wide in yeast. Our single nucleosome positioning data reveal that nucleosomes occupy preferred positions that differ by integral multiples of the DNA helical repeat. By comparing nucleosome dyad positioning maps to existing genomic and transcriptomic data, we evaluated the contributions of sequence, transcription, and histones H1 and H2A.Z in defining the chromatin landscape. We present a biophysical model that neglects DNA sequence and shows that steric occlusion suffices to explain the salient features of nucleosome positioning.

Subtracting the sequence bias from partially digested MNase-seq data reveals a general contribution of TFIIS to nucleosome positioning

BACKGROUND

TFIIS stimulates RNA cleavage by RNA polymerase II and promotes the resolution of backtracking events. TFIIS acts in the chromatin context, but its contribution to the chromatin landscape has not yet been investigated. Co-transcriptional chromatin alterations include subtle changes in nucleosome positioning, like those expected to be elicited by TFIIS, which are elusive to detect. The most popular method to map nucleosomes involves intensive chromatin digestion by micrococcal nuclease (MNase). Maps based on these exhaustively digested samples miss any MNase-sensitive nucleosomes caused by transcription. In contrast, partial digestion approaches preserve such nucleosomes, but introduce noise due to MNase sequence preferences. A systematic way of correcting this bias for massively parallel sequencing experiments is still missing.

RESULTS

To investigate the contribution of TFIIS to the chromatin landscape, we developed a refined nucleosome-mapping method in Saccharomyces cerevisiae. Based on partial MNase digestion and a sequence-bias correction derived from naked DNA cleavage, the refined method efficiently mapped nucleosomes in promoter regions rich in MNase-sensitive structures. The naked DNA correction was also important for mapping gene body nucleosomes, particularly in those genes whose core promoters contain a canonical TATA element. With this improved method, we analyzed the global nucleosomal changes caused by lack of TFIIS. We detected a general increase in nucleosomal fuzziness and more restricted changes in nucleosome occupancy, which concentrated in some gene categories. The TATA-containing genes were preferentially associated with decreased occupancy in gene bodies, whereas the TATA-like genes did so with increased fuzziness. The detected chromatin alterations correlated with functional defects in nascent transcription, as revealed by genomic run-on experiments.

CONCLUSIONS

The combination of partial MNase digestion and naked DNA correction of the sequence bias is a precise nucleosomal mapping method that does not exclude MNase-sensitive nucleosomes. This method is useful for detecting subtle alterations in nucleosome positioning produced by lack of TFIIS. Their analysis revealed that TFIIS generally contributed to nucleosome positioning in both gene promoters and bodies. The independent effect of lack of TFIIS on nucleosome occupancy and fuzziness supports the existence of alternative chromatin dynamics during transcription elongation.

Yeast Terminator Function Can Be Modulated and Designed on the Basis of Predictions of Nucleosome Occupancy

The design of improved synthetic parts is a major goal of synthetic biology. Mechanistically, nucleosome occupancy in the 3' terminator region of a gene has been found to correlate with transcriptional expression. Here, we seek to establish a predictive relationship between terminator function and predicted nucleosome positioning to design synthetic terminators in the yeast Saccharomyces cerevisiae. In doing so, terminators improved net protein output from these expression cassettes nearly 4-fold over their original sequence with observed increases in termination efficiency to 96%. The resulting terminators were indeed depleted of nucleosomes on the basis of mapping experiments. This approach was successfully applied to synthetic, de novo, and native terminators. The mode of action of these modifications was mainly through increased termination efficiency, rather than half-life increases, perhaps suggesting a role in improved mRNA maturation. Collectively, these results suggest that predicted nucleosome depletion can be used as a heuristic approach for improving terminator function, though the underlying mechanism remains to be shown.

MNase-Sensitive Complexes in Yeast: Nucleosomes and Non-histone Barriers

Micrococcal nuclease (MNase) is commonly used to map nucleosomes genome-wide, but nucleosome maps are affected by the degree of digestion. It has been proposed that many yeast promoters are not nucleosome-free but instead occupied by easily digested, unstable, "fragile" nucleosomes. We analyzed the histone content of all MNase-sensitive complexes by MNase-ChIP-seq and sonication-ChIP-seq. We find that yeast promoters are predominantly bound by non-histone protein complexes, with little evidence for fragile nucleosomes. We do detect MNase-sensitive nucleosomes elsewhere in the genome, including at transcription termination sites. However, they have high A/T content, suggesting that MNase sensitivity does not indicate instability, but rather the preference of MNase for A/T-rich DNA, such that A/T-rich nucleosomes are digested faster than G/C-rich nucleosomes. We confirm our observations by analyzing ChIP-exo, chemical mapping, and ATAC-seq data from other laboratories. Thus, histone ChIP-seq experiments are essential to distinguish nucleosomes from other DNA-binding proteins that protect against MNase.

The Ground State and Evolution of Promoter Region Directionality

Eukaryotic promoter regions are frequently divergently transcribed in vivo, but it is unknown whether the resultant antisense RNAs are a mechanistic by-product of RNA polymerase II (Pol II) transcription or biologically meaningful. Here, we use a functional evolutionary approach that involves nascent transcript mapping in S. cerevisiae strains containing foreign yeast DNA. Promoter regions in foreign environments lose the directionality they have in their native species. Strikingly, fortuitous promoter regions arising in foreign DNA produce equal transcription in both directions, indicating that divergent transcription is a mechanistic feature that does not imply a function for these transcripts. Fortuitous promoter regions arising during evolution promote bidirectional transcription and over time are purged through mutation or retained to enable new functionality. Similarly, human transcription is more bidirectional at newly evolved enhancers and promoter regions. Thus, promoter regions are intrinsically bidirectional and are shaped by evolution to bias transcription toward coding versus non-coding RNAs.

Establishing nucleosome architecture and stability at promoters: Roles of pioneer transcription factors and the RSC chromatin remodeler

Improvements in deep sequencing, together with methods to rapidly deplete essential transcription factors (TFs) and chromatin remodelers, have recently led to a more detailed picture of promoter nucleosome architecture in yeast and its relationship to transcriptional regulation. These studies revealed that ∼40% of all budding yeast protein-coding genes possess a unique promoter structure, where we propose that an unusually unstable nucleosome forms immediately upstream of the transcription start site (TSS). This "fragile" nucleosome (FN) promoter architecture relies on the combined action of the essential RSC (Remodels Structure of Chromatin) nucleosome remodeler and pioneer transcription factors (PTFs). FNs are associated with genes whose expression is high, coupled to cell growth, and characterized by low cell-to-cell variability (noise), suggesting that they may promote these features. Recent studies in metazoans suggest that the presence of dynamic nucleosomes upstream of the TSS at highly expressed genes may be conserved throughout evolution.

Transcriptionally Driven DNA Replication Program of the Human Parasite Leishmania major

Faithful inheritance of eukaryotic genomes requires the orchestrated activation of multiple DNA replication origins (ORIs). Although origin firing is mechanistically conserved, how origins are specified and selected for activation varies across different model systems. Here, we provide a complete analysis of the nucleosomal landscape and replication program of the human parasite Leishmania major, building on a better evolutionary understanding of replication organization in Eukarya. We found that active transcription is a driving force for the nucleosomal organization of the L. major genome and that both the spatial and the temporal program of DNA replication can be explained as associated to RNA polymerase kinetics. This simple scenario likely provides flexibility and robustness to deal with the environmental changes that impose alterations in the genetic programs during parasitic life cycle stages. Our findings also suggest that coupling replication initiation to transcription elongation could be an ancient solution used by eukaryotic cells for origin maintenance.

Exploration of nucleosome positioning patterns in transcription factor function

The binding of transcription factors (TFs) triggers activation of specific chromatin regions through the recruitment and activation of RNA polymerase. Unique nucleosome positioning (NP) occurs during gene expression and has been suggested to be involved in various other chromatin functions. However, the diversity of NP that can occur for each function has not been clarified. Here we used MNase-Seq data to evaluate NP around 258 cis-regulatory elements in the mouse genome. Principal component analysis of the 258 elements revealed that NP consisted of five major patterns. Furthermore, the five NP patterns had predictive power for the level of gene expression. We also demonstrated that selective NP patterns appeared around TF binding sites. These results suggest that the NP patterns are correlated to specific functions on chromatin.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions

Transcription factors canonically bind nucleosome-free DNA, making the positioning of nucleosomes within regulatory regions crucial to the regulation of gene expression. Using the assay of transposase accessible chromatin (ATAC-seq), we observe a highly structured pattern of DNA fragment lengths and positions around nucleosomes in Saccharomyces cerevisiae, and use this distinctive two-dimensional nucleosomal “fingerprint” as the basis for a new nucleosome-positioning algorithm called NucleoATAC. We show that NucleoATAC can identify the rotational and translational positions of nucleosomes with up to base-pair resolution and provide quantitative measures of nucleosome occupancy in S. cerevisiae, Schizosaccharomyces pombe, and human cells. We demonstrate the application of NucleoATAC to a number of outstanding problems in chromatin biology, including analysis of sequence features underlying nucleosome positioning, promoter chromatin architecture across species, identification of transient changes in nucleosome occupancy and positioning during a dynamic cellular response, and integrated analysis of nucleosome occupancy and transcription factor binding.

Annotation of genomics data using bidirectional hidden Markov models unveils variations in Pol II transcription cycle

DNA replication, transcription and repair involve the recruitment of protein complexes that change their composition as they progress along the genome in a directed or strand-specific manner. Chromatin immunoprecipitation in conjunction with hidden Markov models (HMMs) has been instrumental in understanding these processes, as they segment the genome into discrete states that can be related to DNA-associated protein complexes. However, current HMM-based approaches are not able to assign forward or reverse direction to states or properly integrate strand-specific (e.g., RNA expression) with non-strand-specific (e.g., ChIP) data, which is indispensable to accurately characterize directed processes. To overcome these limitations, we introduce bidirectional HMMs which infer directed genomic states from occupancy profiles de novo. Application to RNA polymerase II-associated factors in yeast and chromatin modifications in human T cells recovers the majority of transcribed loci, reveals gene-specific variations in the yeast transcription cycle and indicates the existence of directed chromatin state patterns at transcribed, but not at repressed, regions in the human genome. In yeast, we identify 32 new transcribed loci, a regulated initiation–elongation transition, the absence of elongation factors Ctk1 and Paf1 from a class of genes, a distinct transcription mechanism for highly expressed genes and novel DNA sequence motifs associated with transcription termination. We anticipate bidirectional HMMs to significantly improve the analyses of genome-associated directed processes.

Topoisomerases, chromatin and transcription termination

In eukaryotes transcription is complicated by the DNA being packed in nucleosomes and by supercoils induced by opening of the DNA double helix during elongation. Here we discuss our recent genome-wide work regarding topoisomerases and their role in chromatin remodeling during the transcription cycle and we report a novel function for topoisomerases in transcription termination.

Genome-wide mapping of yeast RNA polymerase II termination

Yeast RNA polymerase II (Pol II) terminates transcription of coding transcripts through the polyadenylation (pA) pathway and non-coding transcripts through the non-polyadenylation (non-pA) pathway. We have used PAR-CLIP to map the position of Pol II genome-wide in living yeast cells after depletion of components of either the pA or non-pA termination complexes. We show here that Ysh1, responsible for cleavage at the pA site, is required for efficient removal of Pol II from the template. Depletion of Ysh1 from the nucleus does not, however, lead to readthrough transcription. In contrast, depletion of the termination factor Nrd1 leads to widespread runaway elongation of non-pA transcripts. Depletion of Sen1 also leads to readthrough at non-pA terminators, but in contrast to Nrd1, this readthrough is less processive, or more susceptible to pausing. The data presented here provide delineation of in vivo Pol II termination regions and highlight differences in the sequences that signal termination of different classes of non-pA transcripts.

Transcription termination is an important regulatory event for both non-coding and coding transcripts. Using high-throughput sequencing, we have mapped RNA Polymerase II's position in the genome after depletion of termination factors from the nucleus. We found that depletion of Ysh1 and Sen1 cause build up of polymerase directly downstream of coding and non-coding genes, respectively. Depletion of Nrd1 causes an increase in polymerase that is distributed up to 1,000 bases downstream of non-coding genes. The depletion of Nrd1 helped us to identify more than 250 unique termination regions for non-coding RNAs. Within this set of newly identified non-coding termination regions, we are further able to classify them based on sequence motif similarities, suggesting a functional role for different terminator motifs. The role of these factors in transcriptional termination of coding and/or non-coding transcripts can be inferred from the effect of polymerase's position downstream of given termination sites. This method of depletion and sequencing can be used to further elucidate other factors whose importance to transcription has yet to be determined.

Identifying and mitigating bias in next-generation sequencing methods for chromatin biology

Next-generation sequencing (NGS) technologies have been used in diverse ways to investigate various aspects of chromatin biology by identifying genomic loci that are bound by transcription factors, occupied by nucleosomes or accessible to nuclease cleavage, or loci that physically interact with remote genomic loci. However, reaching sound biological conclusions from such NGS enrichment profiles requires many potential biases to be taken into account. In this Review, we discuss common ways in which biases may be introduced into NGS chromatin profiling data, approaches to diagnose these biases and analytical techniques to mitigate their effect.

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.

Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast

We measured half-lives of 21,248 mRNA 3' isoforms in yeast by rapidly depleting RNA polymerase II from the nucleus and performing direct RNA sequencing throughout the decay process. Interestingly, half-lives of mRNA isoforms from the same gene, including nearly identical isoforms, often vary widely. Based on clusters of isoforms with different half-lives, we identify hundreds of sequences conferring stabilization or destabilization upon mRNAs terminating downstream. One class of stabilizing element is a polyU sequence that can interact with poly(A) tails, inhibit the association of poly(A)-binding protein, and confer increased stability upon introduction into ectopic transcripts. More generally, destabilizing and stabilizing elements are linked to the propensity of the poly(A) tail to engage in double-stranded structures. Isoforms engineered to fold into 3' stem-loop structures not involving the poly(A) tail exhibit even longer half-lives. We suggest that double-stranded structures at 3' ends are a major determinant of mRNA stability.

Simultaneous mapping of transcript ends at single-nucleotide resolution and identification of widespread promoter-associated non-coding RNA governed by TATA elements

Understanding the relationships between regulatory factor binding, chromatin structure, cis-regulatory elements and RNA-regulation mechanisms relies on accurate information about transcription start sites (TSS) and polyadenylation sites (PAS). Although several approaches have identified transcript ends in yeast, limitations of resolution and coverage have remained, and definitive identification of TSS and PAS with single-nucleotide resolution has not yet been achieved. We developed SMORE-seq (simultaneous mapping of RNA ends by sequencing) and used it to simultaneously identify the strongest TSS for 5207 (90%) genes and PAS for 5277 (91%) genes. The new transcript annotations identified by SMORE-seq showed improved distance relationships with TATA-like regulatory elements, nucleosome positions and active RNA polymerase. We found 150 genes whose TSS were downstream of the annotated start codon, and additional analysis of evolutionary conservation and ribosome footprinting suggests that these protein-coding sequences are likely to be mis-annotated. SMORE-seq detected short non-coding RNAs transcribed divergently from more than a thousand promoters in wild-type cells under normal conditions. These divergent non-coding RNAs were less evident at promoters containing canonical TATA boxes, suggesting a model where transcription initiation at promoters by RNAPII is bidirectional, with TATA elements serving to constrain the directionality of initiation.

Independent RNA polymerase II preinitiation complex dynamics and nucleosome turnover at promoter sites in vivo

Transcription by all three eukaryotic RNA polymerases involves the assembly of a large preinitiation complex (PIC) at gene promoters. The PIC comprises several general transcription factors (GTFs), including TBP, and the respective RNA polymerase. It has been suggested that some GTFs remain stably bound at active promoters to facilitate multiple transcription events. Here we used two complementary approaches to show that, in G1-arrested yeast cells, TBP exchanges very rapidly even at the most highly active RNA Pol II promoters. A similar situation is observed at RNA Pol III promoters. In contrast, TBP remains stably bound at RNA Pol I promoters. We also provide evidence that, unexpectedly, PIC dynamics are neither the cause nor the consequence of nucleosome exchange at most of the RNA Pol II promoters we analyzed. These results point to a stable reinitiation complex at RNA Pol I promoters and suggest independent PIC and nucleosome turnover at many RNA Pol II promoters.

Nucleosome-free region dominates histone acetylation in targeting SWR1 to promoters for H2A.Z replacement

The histone variant H2A.Z is a genome-wide signature of nucleosomes proximal to eukaryotic regulatory DNA. Whereas the multisubunit chromatin remodeler SWR1 is known to catalyze ATP-dependent deposition of H2A.Z, the mechanism of SWR1 recruitment to S. cerevisiae promoters has been unclear. A sensitive assay for competitive binding of dinucleosome substrates revealed that SWR1 preferentially binds long nucleosome-free DNA and the adjoining nucleosome core particle, allowing discrimination of gene promoters over gene bodies. Analysis of mutants indicates that the conserved Swc2/YL1 subunit and the adenosine triphosphatase domain of Swr1 are mainly responsible for binding to substrate. SWR1 binding is enhanced on nucleosomes acetylated by the NuA4 histone acetyltransferase, but recognition of nucleosome-free and nucleosomal DNA is dominant over interaction with acetylated histones. Such hierarchical cooperation between DNA and histone signals expands the dynamic range of genetic switches, unifying classical gene regulation by DNA-binding factors with ATP-dependent nucleosome remodeling and posttranslational histone modifications.

A unified model for yeast transcript definition

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

Gene promoters dictate histone occupancy within genes

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

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

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

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

Determinants of nucleosome positioning

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

ATR-like kinase Mec1 facilitates both chromatin accessibility at DNA replication forks and replication fork progression during replication stress

Faithful DNA replication is essential for normal cell division and differentiation. In eukaryotic cells, DNA replication takes place on chromatin. This poses the critical question as to how DNA replication can progress through chromatin, which is inhibitory to all DNA-dependent processes. Here, we developed a novel genome-wide method to measure chromatin accessibility to micrococcal nuclease (MNase) that is normalized for nucleosome density, the NCAM (normalized chromatin accessibility to MNase) assay. This method enabled us to discover that chromatin accessibility increases specifically at and ahead of DNA replication forks in normal S phase and during replication stress. We further found that Mec1, a key regulatory ATR-like kinase in the S-phase checkpoint, is required for both normal chromatin accessibility around replication forks and replication fork rate during replication stress, revealing novel functions for the kinase in replication stress response. These results suggest a possibility that Mec1 may facilitate DNA replication fork progression during replication stress by increasing chromatin accessibility around replication forks.

Histone H3R2 symmetric dimethylation and histone H3K4 trimethylation are tightly correlated in eukaryotic genomes

The preferential in vitro interaction of the PHD finger of RAG2, a subunit of the V(D)J recombinase, with histone H3 tails simultaneously trimethylated at lysine 4 and symmetrically dimethylated at arginine 2 (H3R2me2sK4me3) predicted the existence of the previously unknown histone modification H3R2me2s. Here, we report the in vivo identification of H3R2me2s . Consistent with the binding specificity of the RAG2 PHD finger, high levels of H3R2me2sK4me3 are found at antigen receptor gene segments ready for rearrangement. However, this double modification is much more general; it is conserved throughout eukaryotic evolution. In mouse, H3R2me2s is tightly correlated with H3K4me3 at active promoters throughout the genome. Mutational analysis in S. cerevisiae reveals that deposition of H3R2me2s requires the same Set1 complex that deposits H3K4me3. Our work suggests that H3R2me2sK4me3, not simply H3K4me3 alone, is the mark of active promoters and that factors that recognize H3K4me3 will have their binding modulated by their preference for H3R2me2s.

Genome-wide chromatin remodeling identified at GC-rich long nucleosome-free regions

To gain deeper insights into principles of cell biology, it is essential to understand how cells reorganize their genomes by chromatin remodeling. We analyzed chromatin remodeling on next generation sequencing data from resting and activated T cells to determine a whole-genome chromatin remodeling landscape. We consider chromatin remodeling in terms of nucleosome repositioning which can be observed most robustly in long nucleosome-free regions (LNFRs) that are occupied by nucleosomes in another cell state. We found that LNFR sequences are either AT-rich or GC-rich, where nucleosome repositioning was observed much more prominently in GC-rich LNFRs - a considerable proportion of them outside promoter regions. Using support vector machines with string kernels, we identified a GC-rich DNA sequence pattern indicating loci of nucleosome repositioning in resting T cells. This pattern appears to be also typical for CpG islands. We found out that nucleosome repositioning in GC-rich LNFRs is indeed associated with CpG islands and with binding sites of the CpG-island-binding ZF-CXXC proteins KDM2A and CFP1. That this association occurs prominently inside and also prominently outside of promoter regions hints at a mechanism governing nucleosome repositioning that acts on a whole-genome scale.

Genome-wide nucleosome positioning during embryonic stem cell development

We determined genome-wide nucleosome occupancies in mouse embryonic stem cells and their neural progenitor and embryonic fibroblast counterparts to assess features associated with nucleosome positioning during lineage commitment. Cell-type- and protein-specific binding preferences of transcription factors to sites with either low (Myc, Klf4 and Zfx) or high (Nanog, Oct4 and Sox2) nucleosome occupancy as well as complex patterns for CTCF were identified. Nucleosome-depleted regions around transcription start and transcription termination sites were broad and more pronounced for active genes, with distinct patterns for promoters classified according to CpG content or histone methylation marks. Throughout the genome, nucleosome occupancy was correlated with certain histone methylation or acetylation modifications. In addition, the average nucleosome repeat length increased during differentiation by 5-7 base pairs, with local variations for specific regions. Our results reveal regulatory mechanisms of cell differentiation that involve nucleosome repositioning.

Transcriptional activation of yeast genes disrupts intragenic nucleosome phasing

Nucleosomes often undergo extensive rearrangement when genes are activated for transcription. We have shown previously, using paired-end sequencing of yeast nucleosomes, that major changes in chromatin structure occur when genes are activated by 3-aminotriazole (3AT), an inducer of the transcriptional activator Gcn4. Here, we provide a global analysis of these data. At the genomic level, nucleosomes are regularly phased relative to the transcription start site. However, for a subset of 234 strongly induced genes, this phasing is much more irregular after induction, consistent with the loss of some nucleosomes and the re-positioning of the remaining nucleosomes. To address the nature of this rearrangement, we developed the inter-nucleosome distance auto-correlation (DAC) function. At long range, DAC analysis indicates that nucleosomes have an average spacing of 162 bp, consistent with the reported repeat length. At short range, DAC reveals a 10.25-bp periodicity, implying that nucleosomes in overlapping positions are rotationally related. DAC analysis of the 3AT-induced genes suggests that transcription activation coincides with rearrangement of nucleosomes into irregular arrays with longer spacing. Sequence analysis of the +1 nucleosomes belonging to the 45 most strongly activated genes reveals a distinctive periodic oscillation in the A/T-dinucleotide occurrence that is present throughout the nucleosome and extends into the linker. This unusual pattern suggests that the +1 nucleosomes might be prone to sliding, thereby facilitating transcription.

A map of nucleosome positions in yeast at base-pair resolution

The exact positions of nucleosomes along genomic DNA can influence many aspects of chromosome function, yet existing methods for mapping nucleosomes do not provide the necessary single base pair accuracy to determine these positions. Here we develop and apply a new approach for direct mapping of nucleosome centers based on chemical modification of engineered histones. The resulting map locates nucleosome positions genome-wide in unprecedented detail and accuracy. It reveals novel aspects of the in vivo nucleosome organization that are linked to transcription factor binding, RNA polymerase pausing, and the higher order structure of the chromatin fiber.

A functional evolutionary approach to identify determinants of nucleosome positioning: a unifying model for establishing the genome-wide pattern

Although the genomic pattern of nucleosome positioning is broadly conserved, quantitative aspects vary over evolutionary timescales. We identify the cis and trans determinants of nucleosome positioning using a functional evolutionary approach involving S. cerevisiae strains containing large genomic regions from other yeast species. In a foreign species, nucleosome depletion at promoters is maintained over poly(dA:dT) tracts, whereas internucleosome spacing and all other aspects of nucleosome positioning tested are not. Interestingly, the locations of the +1 nucleosome and RNA start sites shift in concert. Strikingly, in a foreign species, nucleosome-depleted regions occur fortuitously in coding regions, and they often act as promoters that are associated with a positioned nucleosome array linked to the length of the transcription unit. We suggest a three-step model in which nucleosome remodelers, general transcription factors, and the transcriptional elongation machinery are primarily involved in generating the nucleosome positioning pattern in vivo.