Genome-wide mapping of nucleosomes in yeast using paired-end sequencing

Sequence and chromatin features guide DNA double-strand break resection initiation

DNA double-strand breaks (DSBs) are cytotoxic genome lesions that must be accurately and efficiently repaired to ensure genome integrity. In yeast, the Mre11-Rad50-Xrs2 (MRX) complex nicks 5'-terminated DSB ends to initiate nucleolytic processing of DSBs for repair by homologous recombination. How MRX-DNA interactions support 5' strand-specific nicking and how nicking is influenced by the chromatin context have remained elusive. Using a deep sequencing-based assay, we mapped MRX nicks at single-nucleotide resolution next to multiple DSBs in the yeast genome. We observed that the DNA end-binding Ku70-Ku80 complex directed DSB-proximal nicks and that repetitive MRX cleavage extended the length of resection tracts. We identified a sequence motif and a DNA meltability profile that is preferentially nicked by MRX. Furthermore, we found that nucleosomes as well as transcription impeded MRX incisions. Our findings suggest that local DNA sequence and chromatin features shape the activity of this central DSB repair complex.

Identification and analysis of putative tRNA genes in baculovirus genomes

Transfer RNAs (tRNAs) genes are both coded for and arranged along some viral genomes representing the entire virosphere and seem to play different biological functions during infection, other than transferring the correct amino acid to a growing peptide chain. Baculovirus genome description and annotation has focused mostly on protein-coding genes, microRNA, and homologous regions. Here we carried out a large-scale in silico search for putative tRNA genes in baculovirus genomes. Ninety-six of 257 baculovirus genomes analyzed was found to contain at least one putative tRNA gene. We found great diversity in primary and secondary structure, in location within the genome, in intron presence and size, and in anti-codon identity. In some cases, genes of tRNA-containing genomes were found to have a bias for the codons specified by the tRNAs present in such genomes. Moreover, analysis revealed that most of the putative tRNA genes possessed conserved motifs for tRNA type 2 promoters, including the A-box and B-box motifs with few mismatches from the eukaryotic canonical motifs. From publicly available small RNA deep sequencing datasets of baculovirus-infected insect cells, we found evidence that a putative Autographa californica multiple nucleopolyhedrovirus Gln-tRNA gene was transcribed and modified with the addition of the non-templated 3'-CCA tail found at the end of all tRNAs. Further research is needed to determine the expression and functionality of these viral tRNAs.

Trypanosoma brucei histones are heavily modified with combinatorial post-translational modifications and mark Pol II transcription start regions with hyperacetylated H2A

Trypanosomes diverged from the main eukaryotic lineage about 600 million years ago, and display some unusual genomic and epigenetic properties that provide valuable insight into the early processes employed by eukaryotic ancestors to regulate chromatin-mediated functions. We analysed Trypanosoma brucei core histones by high mass accuracy middle-down mass spectrometry to map core histone post-translational modifications (PTMs) and elucidate cis-histone combinatorial PTMs (cPTMs). T. brucei histones are heavily modified and display intricate cPTMs patterns, with numerous hypermodified cPTMs that could contribute to the formation of non-repressive euchromatic states. The Trypanosoma brucei H2A C-terminal tail is hyperacetylated, containing up to five acetylated lysine residues. MNase-ChIP-seq revealed a striking enrichment of hyperacetylated H2A at Pol II transcription start regions, and showed that H2A histones that are hyperacetylated in different combinations localised to different genomic regions, suggesting distinct epigenetic functions. Our genomics and proteomics data provide insight into the complex epigenetic mechanisms used by this parasite to regulate a genome that lacks the transcriptional control mechanisms found in later-branched eukaryotes. The findings further demonstrate the complexity of epigenetic mechanisms that were probably shared with the last eukaryotic common ancestor.

Distinct structural groups of histone H3 and H4 residues have divergent effects on chronological lifespan in Saccharomyces cerevisiae

We have performed a comprehensive analysis of the involvement of histone H3 and H4 residues in the regulation of chronological lifespan in yeast and identify four structural groups in the nucleosome that influence lifespan. We also identify residues where substitution with an epigenetic mimic extends lifespan, providing evidence that a simple epigenetic switch, without possible additional background modifications, causes longevity. Residues where substitution result in the most pronounced lifespan extension are all on the exposed face of the nucleosome, with the exception of H3E50, which is present on the lateral surface, between two DNA gyres. Other residues that have a more modest effect on lifespan extension are concentrated at the extremities of the H3-H4 dimer, suggesting a role in stabilizing the dimer in its nucleosome frame. Residues that reduce lifespan are buried in the histone handshake motif, suggesting that these mutations destabilize the octamer structure. All residues exposed on the nucleosome disk face and that cause lifespan extension are known to interact with Sir3. We find that substitution of H4K16 and H4H18 cause Sir3 to redistribute from telomeres and silent mating loci to secondary positions, often enriched for Rap1, Abf1 or Reb1 binding sites, whereas H3E50 does not. The redistribution of Sir3 in the genome can be reproduced by an equilibrium model based on primary and secondary binding sites with different affinities for Sir3. The redistributed Sir3 cause transcriptional repression at most of the new loci, including of genes where null mutants were previously shown to extend chronological lifespan. The transcriptomic profiles of H4K16 and H4H18 mutant strains are very similar, and compatible with a DNA replication stress response. This is distinct from the transcriptomic profile of H3E50, which matches strong induction of oxidative phosphorylation. We propose that the different groups of residues are involved in binding to heterochromatin proteins, in destabilizing the association of the nucleosome DNA, disrupting binding of the H3-H4 dimer in the nucleosome, or disrupting the structural stability of the octamer, each category impacting on chronological lifespan by a different mechanism.

A systematic genome-wide account of binding sites for the model transcription factor Gcn4

Sequence-specific DNA-binding transcription factors are central to gene regulation. They are often associated with consensus binding sites that predict far more genomic sites than are bound in vivo. One explanation is that most sites are blocked by nucleosomes, such that only sites in nucleosome-depleted regulatory regions are bound. We compared the binding of the yeast transcription factor Gcn4 in vivo using published ChIP-seq data (546 sites) and in vitro, using a modified SELEX method (“G-SELEX”), which utilizes short genomic DNA fragments to quantify binding at all sites. We confirm that Gcn4 binds strongly to an AP-1-like sequence (TGACTCA) and weakly to half-sites. However, Gcn4 binds only some of the 1078 exact matches to this sequence, even in vitro. We show that there are only 166 copies of the high-affinity RTGACTCAY site (exact match) in the yeast genome, all occupied in vivo, largely independently of whether they are located in nucleosome-depleted or nucleosomal regions. Generally, RTGACTCAR/YTGACTCAY sites are bound much more weakly and YTGACTCAR sites are unbound, with biological implications for determining induction levels. We conclude that, to a first approximation, Gcn4 binding can be predicted using the high-affinity site, without reference to chromatin structure. We propose that transcription factor binding sites should be defined more precisely using quantitative data, allowing more accurate genome-wide prediction of binding sites and greater insight into gene regulation.

Measuring the buffering capacity of gene silencing in Saccharomyces cerevisiae

Gene silencing, once established, is stably maintained for several generations. Despite the high fidelity of the inheritance of the silent state, individual components of silenced chromatin are in constant flux. Models suggest that silent loci can tolerate fluctuations in Sir proteins and histone acetylation levels, but the level of tolerance is unknown. To understand the quantitative relationships between H4K16 acetylation, Sir proteins, and silencing, we developed assays to quantitatively alter a H4K16 acetylation mimic allele and Sir protein levels and measure the effects of these changes on silencing. Our data suggest that a two- to threefold change in levels of histone marks and specific Sir proteins affects the stability of the silent state of a large chromatin domain.

Gene silencing in budding yeast is mediated by Sir protein binding to unacetylated nucleosomes to form a chromatin structure that inhibits transcription. Transcriptional silencing is characterized by the high-fidelity transmission of the silent state. Despite its relative stability, the constituent parts of the silent state are in constant flux, giving rise to a model that silent loci can tolerate such fluctuations without functional consequences. However, the level of tolerance is unknown, and we developed methods to measure the threshold of histone acetylation that causes the silent chromatin state to switch to the active state as well as to measure the levels of the enzymes and structural proteins necessary for silencing. We show that loss of silencing required 50 to 75% acetyl-mimic histones, though the precise levels were influenced by silencer strength and upstream activating sequence (UAS) enhancer/promoter strength. Measurements of repressor protein levels necessary for silencing showed that reducing SIR4 gene dosage two- to threefold significantly weakened silencing, though reducing the gene copy numbers for Sir2 or Sir3 to the same extent did not significantly affect silencing suggesting that Sir4 was a limiting component in gene silencing. Calculations suggest that a mere twofold reduction in the ability of acetyltransferases to acetylate nucleosomes across a large array of nucleosomes may be sufficient to generate a transcriptionally silent domain.

Rapid and inexpensive preparation of genome-wide nucleosome footprints from model and non-model organisms

MNase-seq (micrococcal nuclease sequencing) is used to map nucleosome positions in eukaryotic genomes to study the relationship between chromatin structure and DNA-dependent processes. Current protocols require at least two days to isolate nucleosome-protected DNA fragments. We have developed a streamlined protocol for S. cerevisiae and other fungi which takes only three hours. Modified protocols were developed for wild fungi and mammalian cells. This method for rapidly producing sequencing-ready nucleosome footprints from several organisms makes MNase-seq faster and easier, with less chemical waste.

•

A fast way to prepare micrococcal nuclease nucleosome footprints for MNase-seq

•

Eliminates use of phenol and chloroform and reduces the amount of cells required

•

Adaptable for a variety of organisms

tRNAs as a Driving Force of Genome Evolution in Yeast

Transfer RNAs (tRNAs) are widely known for their roles in the decoding of the linear mRNA information into amino acid sequences of proteins. They are also multifunctional platforms in the translation process and have other roles beyond translation, including sensing amino acid abundance, interacting with the general stress response machinery, and modulating cellular adaptation, survival, and death. In this mini-review, we focus on the emerging role of tRNA genes in the organization and modification of the genomic architecture of yeast and the role of tRNA misexpression and decoding infidelity in genome stability, evolution, and adaption. We discuss published work showing how quickly tRNA genes can mutate to meet novel translational demands, how tRNAs speed up genome evolution, and how tRNA genes can be sites of genomic instability. We highlight recent works showing that loss of tRNA decoding fidelity and small alterations in tRNA expression have unexpected and profound impacts on genome stability. By dissecting these recent evidence, we hope to lay the groundwork that prompts future investigations on the mechanistic interplay between tRNAs and genome modification that likely triggers genome evolution.

Global histone protein surface accessibility in yeast indicates a uniformly loosely packed genome with canonical nucleosomes

Background

The vast majority of methods available to characterize genome-wide chromatin structure exploit differences in DNA accessibility to nucleases or chemical crosslinking. We developed a novel method to gauge genome-wide accessibility of histone protein surfaces within nucleosomes by assessing reactivity of engineered cysteine residues with a thiol-specific reagent, biotin-maleimide (BM).

Results

Yeast nuclei were obtained from cells expressing the histone mutant H2B S116C, in which a cysteine resides near the center of the external flat protein surface of the nucleosome. BM modification revealed that nucleosomes are generally equivalently accessible throughout the S. cerevisiae genome, including heterochromatic regions, suggesting limited, higher-order chromatin structures in which this surface is obstructed by tight nucleosome packing. However, we find that nucleosomes within 500 bp of transcription start sites exhibit the greatest range of accessibility, which correlates with the density of chromatin remodelers. Interestingly, accessibility is not well correlated with RNA polymerase density and thus the level of gene expression. We also investigated the accessibility of cysteine mutations designed to detect exposure of histone surfaces internal to the nucleosome thought to be accessible in actively transcribed genes: H3 102, is at the H2A–H2B dimer/H3–H4 tetramer interface, and H3 A110C, resides at the H3–H3 interface. However, in contrast to the external surface site, we find that neither of these internal sites were found to be appreciably exposed.

Conclusions

Overall, our finding that nucleosomes surfaces within S. cerevisiae chromatin are equivalently accessible genome-wide is consistent with a globally uncompacted chromatin structure lacking substantial higher-order organization. However, we find modest differences in accessibility that correlate with chromatin remodelers but not transcription, suggesting chromatin poised for transcription is more accessible than actively transcribed or intergenic regions. In contrast, we find that two internal sites remain inaccessible, suggesting that such non-canonical nucleosome species generated during transcription are rapidly and efficiently converted to canonical nucleosome structure and thus not widely present in native chromatin.

Interplay among ATP-Dependent Chromatin Remodelers Determines Chromatin Organisation in Yeast

Cellular DNA is packaged into chromatin, which is composed of regularly-spaced nucleosomes with occasional gaps corresponding to active regulatory elements, such as promoters and enhancers, called nucleosome-depleted regions (NDRs). This chromatin organisation is primarily determined by the activities of a set of ATP-dependent remodeling enzymes that are capable of moving nucleosomes along DNA, or of evicting nucleosomes altogether. In yeast, the nucleosome-spacing enzymes are ISW1 (Imitation SWitch protein 1), Chromodomain-Helicase-DNA-binding (CHD)1, ISW2 (Imitation SWitch protein 2) and INOsitol-requiring 80 (INO80); the nucleosome eviction enzymes are the SWItching/Sucrose Non-Fermenting (SWI/SNF) family, the Remodeling the Structure of Chromatin (RSC) complexes and INO80. We discuss the contributions of each set of enzymes to chromatin organisation. ISW1 and CHD1 are the major spacing enzymes; loss of both enzymes results in major chromatin disruption, partly due to the appearance of close-packed di-nucleosomes. ISW1 and CHD1 compete to set nucleosome spacing on most genes. ISW1 is dominant, setting wild type spacing, whereas CHD1 sets short spacing and may dominate on highly-transcribed genes. We propose that the competing remodelers regulate spacing, which in turn controls the binding of linker histone (H1) and therefore the degree of chromatin folding. Thus, genes with long spacing bind more H1, resulting in increased chromatin compaction. RSC, SWI/SNF and INO80 are involved in NDR formation, either directly by nucleosome eviction or repositioning, or indirectly by affecting the size of the complex that resides in the NDR. The nature of this complex is controversial: some suggest that it is a RSC-bound “fragile nucleosome”, whereas we propose that it is a non-histone transcription complex. In either case, this complex appears to serve as a barrier to nucleosome formation, resulting in the formation of phased nucleosomal arrays on both sides.

Q-Nuc: a bioinformatics pipeline for the quantitative analysis of nucleosomal profiles

Nucleosomal profiling is an effective method to determine the positioning and occupancy of nucleosomes, which is essential to understand their roles in genomic processes. However, the positional randomness across the genome and its relationship with nucleosome occupancy remains poorly understood. Here we present a computational method that segments the profile into nucleosomal domains and quantifies their randomness and relative occupancy level. Applying this method to published data, we find on average ~ 3-fold differences in the degree of positional randomness between regions typically considered "well-ordered", as well as an unexpected predominance of only two types of domains of positional randomness in yeast cells. Further, we find that occupancy levels between domains actually differ maximally by ~ 2-3-fold in both cells, which has not been described before. We also developed a procedure by which one can estimate the sequencing depth that is required to identify nucleosomal positions even when regional positional randomness is high. Overall, we have developed a pipeline to quantitatively characterize domain-level features of nucleosome randomness and occupancy genome-wide, enabling the identification of otherwise unknown features in nucleosomal organization.

A method for assessing histone surface accessibility genome-wide

The assembly of DNA into nucleosomes and higher order chromatin structures serves not only as a means of compaction but also organizes the genome to facilitate crucial processes such as cell division and regulation of gene expression. Chromatin structure generally limits access to DNA, with the accessibility of DNA in chromatin being regulated through post translational modification of the histone proteins as well as the activity of chromatin remodeling proteins and architectural chromatin factors. There is great interest in assessing chromatin accessibility genome-wide to identify functional elements associated with enhancers, promoters, and other discontinuities in the compacted chromatin structure associated with gene expression. As the vast majority of techniques rely upon assessment of the exposure of the underlying DNA, we describe here a general method that can be used to assess exposure of internal and external histone protein surfaces. We demonstrate the feasibility of this method, in the organism S. cerevisiae. Our method relies on substitution of residues residing on selected histone protein surfaces with cysteine, and assessment of exposure by reaction with a thiol specific reagent, biotin-maleimide. We demonstrate that modified nucleosomes can be efficiently excised from nuclei treated with the reagent via a one-step purification process. After library preparation and deep sequencing, selected nucleosomes are typically ~25-fold enriched over background signals and exhibit phasing with respect to transcription start sites in yeast that is identical to an unselected population.

tRNA Genes Affect Chromosome Structure and Function via Local Effects

The genome is packaged and organized in an ordered, nonrandom manner, and specific chromatin segments contact nuclear substructures to mediate this organization. tRNA genes (tDNAs) are binding sites for transcription factors and architectural proteins and are thought to play an important role in the organization of the genome. In this study, we investigate the roles of tDNAs in genomic organization and chromosome function by editing a chromosome so that it lacked any tDNAs. Surprisingly our analyses of this tDNA-less chromosome show that loss of tDNAs does not grossly affect chromatin architecture or chromosome tethering and mobility. However, loss of tDNAs affects local nucleosome positioning and the binding of SMC proteins at these loci. The absence of tDNAs also leads to changes in centromere clustering and a reduction in the frequency of long-range HML-HMR heterochromatin clustering with concomitant effects on gene silencing. We propose that the tDNAs primarily affect local chromatin structure, which results in effects on long-range chromosome architecture.

Contrasting roles of the RSC and ISW1/CHD1 chromatin remodelers in RNA polymerase II elongation and termination

Most yeast genes have a nucleosome-depleted region (NDR) at the promoter and an array of regularly spaced nucleosomes phased relative to the transcription start site. We have examined the interplay between RSC (a conserved essential SWI/SNF-type complex that determines NDR size) and the ISW1, CHD1, and ISW2 nucleosome spacing enzymes in chromatin organization and transcription, using isogenic strains lacking all combinations of these enzymes. The contributions of these remodelers to chromatin organization are largely combinatorial, distinct, and nonredundant, supporting a model in which the +1 nucleosome is positioned by RSC and then used as a reference nucleosome by the spacing enzymes. Defective chromatin organization correlates with altered RNA polymerase II (Pol II) distribution. RSC-depleted cells exhibit low levels of elongating Pol II and high levels of terminating Pol II, consistent with defects in both termination and initiation, suggesting that RSC facilitates both. Cells lacking both ISW1 and CHD1 show the opposite Pol II distribution, suggesting elongation and termination defects. These cells have extremely disrupted chromatin, with high levels of closely packed dinucleosomes involving the second (+2) nucleosome. We propose that ISW1 and CHD1 facilitate Pol II elongation by separating closely packed nucleosomes.

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.

Single-Molecule Analysis Reveals Linked Cycles of RSC Chromatin Remodeling and Ace1p Transcription Factor Binding in Yeast

It is unknown how the dynamic binding of transcription factors (TFs) is molecularly linked to chromatin remodeling and transcription. Using single-molecule tracking (SMT), we show that the chromatin remodeler RSC speeds up the search process of the TF Ace1p for its response elements (REs) at the CUP1 promoter. We quantified smFISH mRNA data using a gene bursting model and demonstrated that RSC regulates transcription bursts of CUP1 only by modulating TF occupancy but does not affect initiation and elongation rates. We show by SMT that RSC binds to activated promoters transiently, and based on MNase-seq data, that RSC does not affect the nucleosomal occupancy at CUP1. Therefore, transient binding of Ace1p and rapid bursts of transcription at CUP1 may be dependent on short repetitive cycles of nucleosome mobilization. This type of regulation reduces the transcriptional noise and ensures a homogeneous response of the cell population to heavy metal stress.

GT-rich promoters can drive RNA pol II transcription and deposition of H2A.Z in African trypanosomes

Genome‐wide transcription studies are revealing an increasing number of “dispersed promoters” that, unlike “focused promoters”, lack well‐conserved sequence motifs and tight regulation. Dispersed promoters are nevertheless marked by well‐defined chromatin structures, suggesting that specific sequence elements must exist in these unregulated promoters. Here, we have analyzed regions of transcription initiation in the eukaryotic parasite Trypanosoma brucei, in which RNA polymerase II transcription initiation occurs over broad regions without distinct promoter motifs and lacks regulation. Using a combination of site‐specific and genome‐wide assays, we identified GT‐rich promoters that can drive transcription and promote the targeted deposition of the histone variant H2A.Z in a genomic context‐dependent manner. In addition, upon mapping nucleosome occupancy at high resolution, we find nucleosome positioning to correlate with RNA pol II enrichment and gene expression, pointing to a role in RNA maturation. Nucleosome positioning may thus represent a previously unrecognized layer of gene regulation in trypanosomes. Our findings show that even highly dispersed, unregulated promoters contain specific DNA elements that are able to induce transcription and changes in chromatin structure.

The proto-chromatosome: A fundamental subunit of chromatin?

Eukaryotic DNA is packaged into regularly spaced nucleosomes, resembling beads on a string. Each bead contains ∼147 bp wrapped around a core histone octamer. Linker histone (H1) binds to the linker DNA to drive chromatin folding. Micrococcal nuclease (MNase) digestion studies reveal 2 mono-nucleosomal intermediates: the core particle (∼147 bp) and the chromatosome (∼160 bp; a core particle with additional DNA protected by H1). We have recently developed an improved method for mapping nucleosomes, using exonuclease III to remove residual linker (MNase-Exo-seq). (1) We discovered 2 new intermediate particles corresponding to core particles with ∼7 bp of linker protruding from one side (∼154 bp) or both sides (∼161 bp), which are formed in the absence of H1. We propose that these "proto-chromatosomes" are stabilized by core histone-DNA contacts in the linker, ∼7 bp from the nucleosome boundaries. These contacts may determine the topography of the H1 binding site.

Nucleosomes influence multiple steps during replication initiation

Eukaryotic replication origin licensing, activation and timing are influenced by chromatin but a mechanistic understanding is lacking. Using reconstituted nucleosomal DNA replication assays, we assessed the impact of nucleosomes on replication initiation. To generate distinct nucleosomal landscapes, different chromatin-remodeling enzymes (CREs) were used to remodel nucleosomes on origin-DNA templates. Nucleosomal organization influenced two steps of replication initiation: origin licensing and helicase activation. Origin licensing assays showed that local nucleosome positioning enhanced origin specificity and modulated helicase loading by influencing ORC DNA binding. Interestingly, SWI/SNF- and RSC-remodeled nucleosomes were permissive for origin licensing but showed reduced helicase activation. Specific CREs rescued replication of these templates if added prior to helicase activation, indicating a permissive chromatin state must be established during origin licensing to allow efficient origin activation. Our studies show nucleosomes directly modulate origin licensing and activation through distinct mechanisms and provide insights into the regulation of replication initiation by chromatin.

Well-positioned nucleosomes punctuate polycistronic pol II transcription units and flank silent VSG gene arrays in Trypanosoma brucei

Background

The compaction of DNA in chromatin in eukaryotes allowed the expansion of genome size and coincided with significant evolutionary diversification. However, chromatin generally represses DNA function, and mechanisms coevolved to regulate chromatin structure and its impact on DNA. This included the selection of specific nucleosome positions to modulate accessibility to the DNA molecule. Trypanosoma brucei, a member of the Excavates supergroup, falls in an ancient evolutionary branch of eukaryotes and provides valuable insight into the organization of chromatin in early genomes.

Results

We have mapped nucleosome positions in T. brucei and identified important differences compared to other eukaryotes: The RNA polymerase II initiation regions in T. brucei do not exhibit pronounced nucleosome depletion, and show little evidence for defined −1 and +1 nucleosomes. In contrast, a well-positioned nucleosome is present directly on the splice acceptor sites within the polycistronic transcription units. The RNA polyadenylation sites were depleted of nucleosomes, with a single well-positioned nucleosome present immediately downstream of the predicted sites. The regions flanking the silent variant surface glycoprotein (VSG) gene cassettes showed extensive arrays of well-positioned nucleosomes, which may repress cryptic transcription initiation. The silent VSG genes themselves exhibited a less regular nucleosomal pattern in both bloodstream and procyclic form trypanosomes. The DNA replication origins, when present within silent VSG gene cassettes, displayed a defined nucleosomal organization compared with replication origins in other chromosomal core regions.

Conclusions

Our results indicate that some organizational features of chromatin are evolutionarily ancient, and may already have been present in the last eukaryotic common ancestor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0121-9) contains supplementary material, which is available to authorized users.

Genome-wide mapping of nucleosome positions in Saccharomyces cerevisiae in response to different nitrogen conditions

Well-organized chromatin is involved in a number of various transcriptional regulation and gene expression. We used genome-wide mapping of nucleosomes in response to different nitrogen conditions to determine both nucleosome profiles and gene expression events in Saccharomyces cerevisiae. Nitrogen conditions influence general nucleosome profiles and the expression of nitrogen catabolite repression (NCR) sensitive genes. The nucleosome occupancy of TATA-containing genes was higher compared to TATA-less genes. TATA-less genes in high or low nucleosome occupancy, showed a significant change in gene coding regions when shifting cells from glutamine to proline as the sole nitrogen resource. Furthermore, a correlation between the expression of nucleosome occupancy induced NCR sensitive genes or TATA containing genes in NCR sensitive genes, and nucleosome prediction were found when cells were cultured in proline or shifting from glutamine to proline as the sole nitrogen source compared to glutamine. These results also showed that variation of nucleosome occupancy accompany with chromatin-dependent transcription factor could influence the expression of a series of genes involved in the specific regulation of nitrogen utilization.

Genome-wide profiling of nucleosome sensitivity and chromatin accessibility in Drosophila melanogaster

Nucleosomal DNA is thought to be generally inaccessible to DNA-binding factors, such as micrococcal nuclease (MNase). Here, we digest Drosophila chromatin with high and low concentrations of MNase to reveal two distinct nucleosome types: MNase-sensitive and MNase-resistant. MNase-resistant nucleosomes assemble on sequences depleted of A/T and enriched in G/C-containing dinucleotides, whereas MNase-sensitive nucleosomes form on A/T-rich sequences found at transcription start and termination sites, enhancers and DNase I hypersensitive sites. Estimates of nucleosome formation energies indicate that MNase-sensitive nucleosomes tend to be less stable than MNase-resistant ones. Strikingly, a decrease in cell growth temperature of about 10°C makes MNase-sensitive nucleosomes less accessible, suggesting that observed variations in MNase sensitivity are related to either thermal fluctuations of chromatin fibers or the activity of enzymatic machinery. In the vicinity of active genes and DNase I hypersensitive sites nucleosomes are organized into periodic arrays, likely due to 'phasing' off potential barriers formed by DNA-bound factors or by nucleosomes anchored to their positions through external interactions. The latter idea is substantiated by our biophysical model of nucleosome positioning and energetics, which predicts that nucleosomes immediately downstream of transcription start sites are anchored and recapitulates nucleosome phasing at active genes significantly better than sequence-dependent models.

The ISW1 and CHD1 ATP-dependent chromatin remodelers compete to set nucleosome spacing in vivo

Adenosine triphosphate-dependent chromatin remodeling machines play a central role in gene regulation by manipulating chromatin structure. Most genes have a nucleosome-depleted region at the promoter and an array of regularly spaced nucleosomes phased relative to the transcription start site. In vitro, the three known yeast nucleosome spacing enzymes (CHD1, ISW1 and ISW2) form arrays with different spacing. We used genome-wide nucleosome sequencing to determine whether these enzymes space nucleosomes differently in vivo We find that CHD1 and ISW1 compete to set the spacing on most genes, such that CHD1 dominates genes with shorter spacing and ISW1 dominates genes with longer spacing. In contrast, ISW2 plays a minor role, limited to transcriptionally inactive genes. Heavily transcribed genes show weak phasing and extreme spacing, either very short or very long, and are depleted of linker histone (H1). Genes with longer spacing are enriched in H1, which directs chromatin folding. We propose that CHD1 directs short spacing, resulting in eviction of H1 and chromatin unfolding, whereas ISW1 directs longer spacing, allowing H1 to bind and condense the chromatin. Thus, competition between the two remodelers to set the spacing on each gene may result in a highly dynamic chromatin structure.

Genome-wide cooperation by HAT Gcn5, remodeler SWI/SNF, and chaperone Ydj1 in promoter nucleosome eviction and transcriptional activation

Chaperones, nucleosome remodeling complexes, and histone acetyltransferases have been implicated in nucleosome disassembly at promoters of particular yeast genes, but whether these cofactors function ubiquitously, as well as the impact of nucleosome eviction on transcription genome-wide, is poorly understood. We used chromatin immunoprecipitation of histone H3 and RNA polymerase II (Pol II) in mutants lacking single or multiple cofactors to address these issues for about 200 genes belonging to the Gcn4 transcriptome, of which about 70 exhibit marked reductions in H3 promoter occupancy on induction by amino acid starvation. Examining four target genes in a panel of mutants indicated that SWI/SNF, Gcn5, the Hsp70 cochaperone Ydj1, and chromatin-associated factor Yta7 are required downstream from Gcn4 binding, whereas Asf1/Rtt109, Nap1, RSC, and H2AZ are dispensable for robust H3 eviction in otherwise wild-type cells. Using ChIP-seq to interrogate all 70 exemplar genes in single, double, and triple mutants implicated Gcn5, Snf2, and Ydj1 in H3 eviction at most, but not all, Gcn4 target promoters, with Gcn5 generally playing the greatest role and Ydj1 the least. Remarkably, these three cofactors cooperate similarly in H3 eviction at virtually all yeast promoters. Defective H3 eviction in cofactor mutants was coupled with reduced Pol II occupancies for the Gcn4 transcriptome and the most highly expressed uninduced genes, but the relative Pol II levels at most genes were unaffected or even elevated. These findings indicate that nucleosome eviction is crucial for robust transcription of highly expressed genes but that other steps in gene activation are more rate-limiting for most other yeast genes.

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.

Novel nucleosomal particles containing core histones and linker DNA but no histone H1

Eukaryotic chromosomal DNA is assembled into regularly spaced nucleosomes, which play a central role in gene regulation by determining accessibility of control regions. The nucleosome contains ∼147 bp of DNA wrapped ∼1.7 times around a central core histone octamer. The linker histone, H1, binds both to the nucleosome, sealing the DNA coils, and to the linker DNA between nucleosomes, directing chromatin folding. Micrococcal nuclease (MNase) digests the linker to yield the chromatosome, containing H1 and ∼160 bp, and then converts it to a core particle, containing ∼147 bp and no H1. Sequencing of nucleosomal DNA obtained after MNase digestion (MNase-seq) generates genome-wide nucleosome maps that are important for understanding gene regulation. We present an improved MNase-seq method involving simultaneous digestion with exonuclease III, which removes linker DNA. Remarkably, we discovered two novel intermediate particles containing 154 or 161 bp, corresponding to 7 bp protruding from one or both sides of the nucleosome core. These particles are detected in yeast lacking H1 and in H1-depleted mouse chromatin. They can be reconstituted in vitro using purified core histones and DNA. We propose that these ‘proto-chromatosomes’ are fundamental chromatin subunits, which include the H1 binding site and influence nucleosome spacing independently of H1.

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.

Heavy transcription of yeast genes correlates with differential loss of histone H2B relative to H4 and queued RNA polymerases

Eukaryotic chromatin is composed of nucleosomes, which contain nearly two coils of DNA wrapped around a central histone octamer. The octamer contains an H3-H4 tetramer and two H2A-H2B dimers. Gene activation is associated with chromatin disruption: a wider nucleosome-depleted region (NDR) at the promoter and reduced nucleosome occupancy over the coding region. Here, we examine the nature of disrupted chromatin after induction, using MNase-seq to map nucleosomes and subnucleosomes, and a refined high-resolution ChIP-seq method to map H4, H2B and RNA polymerase II (Pol II) genome-wide. Over coding regions, induced genes show a differential loss of H2B relative to H4, which correlates with Pol II density and the appearance of subnucleosomes. After induction, Pol II is surprisingly low at the promoter, but accumulates on the gene and downstream of the termination site, implying that dissociation is very slow. Thus, induction-dependent chromatin disruption reflects both eviction of H2A-H2B dimers and the presence of queued Pol II elongation complexes. We propose that slow Pol II dissociation after transcription is a major factor in chromatin disruption and that it may be of critical importance in gene regulation.

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.