A high-resolution atlas of nucleosome occupancy in yeast

Natural antisense transcripts as versatile regulators of gene expression

Long non-coding RNAs (lncRNAs) are emerging as a major class of gene products that have central roles in cell and developmental biology. Natural antisense transcripts (NATs) are an important subset of lncRNAs that are expressed from the opposite strand of protein-coding and non-coding genes and are a genome-wide phenomenon in both eukaryotes and prokaryotes. In eukaryotes, a myriad of NATs participate in regulatory pathways that affect expression of their cognate sense genes. Recent developments in the study of NATs and lncRNAs and large-scale sequencing and bioinformatics projects suggest that whether NATs regulate expression, splicing, stability or translation of the sense transcript is influenced by the pattern and degrees of overlap between the sense-antisense pair. Moreover, epigenetic gene regulatory mechanisms prevail in somatic cells whereas mechanisms dependent on the formation of double-stranded RNA intermediates are prevalent in germ cells. The modulating effects of NATs on sense transcript expression make NATs rational targets for therapeutic interventions.

Exploring the reciprocity between pioneer factors and development

Development is regulated by coordinated changes in gene expression. Control of these changes in expression is largely governed by the binding of transcription factors to specific regulatory elements. However, the packaging of DNA into chromatin prevents the binding of many transcription factors. Pioneer factors overcome this barrier owing to unique properties that enable them to bind closed chromatin, promote accessibility and, in so doing, mediate binding of additional factors that activate gene expression. Because of these properties, pioneer factors act at the top of gene-regulatory networks and drive developmental transitions. Despite the ability to bind target motifs in closed chromatin, pioneer factors have cell type-specific chromatin occupancy and activity. Thus, developmental context clearly shapes pioneer-factor function. Here, we discuss this reciprocal interplay between pioneer factors and development: how pioneer factors control changes in cell fate and how cellular environment influences pioneer-factor binding and activity.

Genome organization across scales: mechanistic insights from in vitro reconstitution studies

Eukaryotic genomes are compacted and organized into distinct three-dimensional (3D) structures, which range from small-scale nucleosome arrays to large-scale chromatin domains. These chromatin structures play an important role in the regulation of transcription and other nuclear processes. The molecular mechanisms that drive the formation of chromatin structures across scales and the relationship between chromatin structure and function remain incompletely understood. Because the processes involved are complex and interconnected, it is often challenging to dissect the underlying principles in the nuclear environment. Therefore, in vitro reconstitution systems provide a valuable approach to gain insight into the molecular mechanisms by which chromatin structures are formed and to determine the cause-consequence relationships between the processes involved. In this review, we give an overview of in vitro approaches that have been used to study chromatin structures across scales and how they have increased our understanding of the formation and function of these structures. We start by discussing in vitro studies that have given insight into the mechanisms of nucleosome positioning. Next, we discuss recent efforts to reconstitute larger-scale chromatin domains and loops and the resulting insights into the principles of genome organization. We conclude with an outlook on potential future applications of chromatin reconstitution systems and how they may contribute to answering open questions concerning chromatin architecture.

Genome-wide regulation of Pol II, FACT, and Spt6 occupancies by RSC in Saccharomyces cerevisiae

RSC (remodels the structure of chromatin) is an essential ATP-dependent chromatin remodeling complex in Saccharomyces cerevisiae. RSC utilizes its ATPase subunit, Sth1, to slide or remove nucleosomes. RSC has been shown to regulate the width of the nucleosome-depleted regions (NDRs) by sliding the flanking nucleosomes away from NDRs. As such, when RSC is depleted, nucleosomes encroach NDRs, leading to transcription initiation defects. In this study, we examined the effects of the catalytic-dead Sth1 on transcription and compared them to those observed during acute and rapid Sth1 depletion by auxin-induced degron strategy. We found that rapid depletion of Sth1 reduces recruitment of TBP and Pol II in highly transcribed genes, as would be expected considering its role in regulating chromatin structure at promoters. In contrast, cells harboring the catalytic-dead Sth1 (sth1-K501R) exhibited a severe reduction in TBP binding, but, surprisingly, also displayed a substantial accumulation in Pol II occupancies within coding regions. The Pol II occupancies further increased upon depleting endogenous Sth1 in the catalytic-dead mutant, suggesting that the inactive Sth1 contributes to Pol II accumulation in coding regions. Notwithstanding the Pol II increase, the ORF occupancies of histone chaperones, FACT and Spt6 were significantly reduced in the mutant. These results suggest a potential role for RSC in recruiting/retaining these chaperones in coding regions. Pol II accumulation despite substantial reductions in TBP, FACT, and Spt6 occupancies in the catalytic-dead mutant could indicate severe transcription elongation and termination defects. Such defects would be consistent with studies showing that RSC is recruited to coding regions in a transcription-dependent manner. Thus, these findings imply a role for RSC in transcription elongation and termination processes, in addition to its established role in transcription initiation.

Quantitative analysis of transcription start site selection reveals control by DNA sequence, RNA polymerase II activity and NTP levels

Transcription start site (TSS) selection is a key step in gene expression and occurs at many promoter positions over a wide range of efficiencies. Here we develop a massively parallel reporter assay to quantitatively dissect contributions of promoter sequence, nucleoside triphosphate substrate levels and RNA polymerase II (Pol II) activity to TSS selection by 'promoter scanning' in Saccharomyces cerevisiae (Pol II MAssively Systematic Transcript End Readout, 'Pol II MASTER'). Using Pol II MASTER, we measure the efficiency of Pol II initiation at 1,000,000 individual TSS sequences in a defined promoter context. Pol II MASTER confirms proposed critical qualities of S. cerevisiae TSS -8, -1 and +1 positions, quantitatively, in a controlled promoter context. Pol II MASTER extends quantitative analysis to surrounding sequences and determines that they tune initiation over a wide range of efficiencies. These results enabled the development of a predictive model for initiation efficiency based on sequence. We show that genetic perturbation of Pol II catalytic activity alters initiation efficiency mostly independently of TSS sequence, but selectively modulates preference for the initiating nucleotide. Intriguingly, we find that Pol II initiation efficiency is directly sensitive to guanosine-5'-triphosphate levels at the first five transcript positions and to cytosine-5'-triphosphate and uridine-5'-triphosphate levels at the second position genome wide. These results suggest individual nucleoside triphosphate levels can have transcript-specific effects on initiation, representing a cryptic layer of potential regulation at the level of Pol II biochemical properties. The results establish Pol II MASTER as a method for quantitative dissection of transcription initiation in eukaryotes.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Predictions of DNA mechanical properties at a genomic scale reveal potentially new functional roles of DNA flexibility

Mechanical properties of DNA have been implied to influence many of its biological functions. Recently, a new high-throughput method, called loop-seq, which allows measuring the intrinsic bendability of DNA fragments, has been developed. Using loop-seq data, we created a deep learning model to explore the biological significance of local DNA flexibility in a range of different species from different kingdoms. Consistently, we observed a characteristic and largely dinucleotide-composition-driven change of local flexibility near transcription start sites. In the presence of a TATA-box, a pronounced peak of high flexibility can be observed. Furthermore, depending on the transcription factor investigated, flanking-sequence-dependent DNA flexibility was identified as a potential factor influencing DNA binding. Compared to randomized genomic sequences, depending on species and taxa, actual genomic sequences were observed both with increased and lowered flexibility. Furthermore, in Arabidopsis thaliana, mutation rates, both de novo and fixed, were found to be associated with relatively rigid sequence regions. Our study presents a range of significant correlations between characteristic DNA mechanical properties and genomic features, the significance of which with regard to detailed molecular relevance awaits further theoretical and experimental exploration.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Improving genome-wide mapping of nucleosomes in Trypanosome cruzi

In Trypanosoma cruzi DNA is packaged into chromatin by octamers of histone proteins that form nucleosomes. Transcription of protein coding genes in trypanosomes is constitutive producing polycistronic units and gene expression is primarily regulated post-transcriptionally. However, chromatin organization influences DNA dependent processes. Hence, determining nucleosome position is of uppermost importance to understand the peculiarities found in trypanosomes. To map nucleosomes genome-wide in several organisms, digestion of chromatin with micrococcal nuclease followed by deep sequencing has been applied. Nonetheless, the special requirements for cell manipulation and the uniqueness of the chromatin organization in trypanosomes entails a customized analytical approach. In this work, we adjusted this broadly used method to the hybrid reference strain, CL Brener. Particularly, we implemented an exhaustive and thorough computational workflow to overcome the difficulties imposed by this complex genome. We tested the performance of two aligners, Bowtie2 and HISAT2, and discuss their advantages and caveats. Specifically, we highlight the relevance of using the whole genome as a reference instead of the commonly used Esmeraldo-like haplotype to avoid spurious alignments. Additionally, we show that using the whole genome refines the average nucleosome representation, but also the quality of mapping for every region represented. Moreover, we show that the average nucleosome organization around trans-splicing acceptor site described before, is not just an average since the same chromatin pattern is detected for most of the represented regions. In addition, we extended the study to a non-hybrid strain applying the experimental and analytical approach to Sylvio-X10 strain. Furthermore, we provide a source code for the construction of 2D plots and heatmaps which are easy to adapt to any T. cruzi strain.

Linking New Alleles at the Oscillator Loci to Flowering and Expansion of Asian Rice

The central oscillator is believed to be the key mechanism by which plants adapt to new environments. However, impacts from hybridization, the natural environment, and human selection have rarely been assessed on the oscillator of a crop. Here, from clearly identified alleles at oscillator loci (OsCCA1/LHY, OsPRR95, OsPRR37, OsPRR59, and OsPRR1) in ten diverse genomes of Oryza sativa, additional accessions, and functional analysis, we show that rice's oscillator was rebuilt primarily by new alleles from recombining parental sequences and subsequent 5' or/and coding mutations. New alleles may exhibit altered transcript levels from that of a parental allele and are transcribed variably among genetic backgrounds and natural environments in RIL lines. Plants carrying more expressed OsCCA1_a and less transcribed OsPRR1_e flower early in the paddy field. 5' mutations are instrumental in varied transcription, as shown by EMSA tests on one deletion at the 5' region of highly transcribed OsPRR1_a. Compared to relatively balanced mutations at oscillator loci of Arabidopsis thaliana, 5' mutations of OsPRR37 (and OsCCA1 to a less degree) were under negative selection while those of OsPRR1 alleles were under strong positive selection. Together, range expansion of Asian rice can be elucidated by human selection on OsPRR1 alleles via local flowering time-yield relationships.

ISW1a modulates cohesin distribution in centromeric and pericentromeric regions

Cohesin is a highly conserved, multiprotein complex whose canonical function is to hold sister chromatids together to ensure accurate chromosome segregation. Cohesin association with chromatin relies on the Scc2-Scc4 cohesin loading complex that enables cohesin ring opening and topological entrapment of sister DNAs. To better understand how sister chromatid cohesion is regulated, we performed a proteomic screen in budding yeast that identified the Isw1 chromatin remodeler as a cohesin binding partner. In addition, we found that Isw1 also interacts with Scc2-Scc4. Lack of Isw1 protein, the Ioc3 subunit of ISW1a or Isw1 chromatin remodeling activity resulted in increased accumulation of cohesin at centromeres and pericentromeres, suggesting that ISW1a may promote efficient translocation of cohesin from the centromeric site of loading to neighboring regions. Consistent with the role of ISW1a in the chromatin organization of centromeric regions, Isw1 was found to be recruited to centromeres. In its absence we observed changes in the nucleosomal landscape at centromeres and pericentromeres. Finally, we discovered that upon loss of RSC functionality, ISW1a activity leads to reduced cohesin binding and cohesion defect. Taken together, our results support the notion of a key role of chromatin remodelers in the regulation of cohesin distribution on chromosomes.

Lola-I is a promoter pioneer factor that establishes de novo Pol II pausing during development

While the accessibility of enhancers is dynamically regulated during development, promoters tend to be constitutively accessible and poised for activation by paused Pol II. By studying Lola-I, a Drosophila zinc finger transcription factor, we show here that the promoter state can also be subject to developmental regulation independently of gene activation. Lola-I is ubiquitously expressed at the end of embryogenesis and causes its target promoters to become accessible and acquire paused Pol II throughout the embryo. This promoter transition is required but not sufficient for tissue-specific target gene activation. Lola-I mediates this function by depleting promoter nucleosomes, similar to the action of pioneer factors at enhancers. These results uncover a level of regulation for promoters that is normally found at enhancers and reveal a mechanism for the de novo establishment of paused Pol II at promoters.

The role of transcript regions and amino acid choice in nucleosome positioning

Eukaryotic DNA is organized and compacted in a string of nucleosomes, DNA-wrapped protein cylinders. The positions of nucleosomes along DNA are not random but show well-known base pair sequence preferences that result from the sequence-dependent elastic and geometric properties of the DNA double helix. Here, we focus on DNA around transcription start sites, which are known to typically attract nucleosomes in multicellular life forms through their high GC content. We aim to understand how these GC signals, as observed in genome-wide averages, are produced and encoded through different genomic regions (mainly 5' UTRs, coding exons, and introns). Our study uses a bioinformatics approach to decompose the genome-wide GC signal into between-region and within-region signals. We find large differences in GC signal contributions between vertebrates and plants and, remarkably, even between closely related species. Introns contribute most to the GC signal in vertebrates, while in plants the exons dominate. Further, we find signal strengths stronger on DNA than on mRNA, suggesting a biological function of GC signals along the DNA itself, as is the case for nucleosome positioning. Finally, we make the surprising discovery that both the choice of synonymous codons and amino acids contribute to the nucleosome positioning signal.

Regulation of IFNβ expression: focusing on the role of its promoter and transcription regulators

IFNβ is a single-copy gene without an intron. Under normal circumstances, it shows low or no expression in cells. It is upregulated only when the body needs it or is stimulated. Stimuli bind to the pattern recognition receptors (PRRs) and pass via various signaling pathways to several basic transcriptional regulators, such as IRFs, NF-кB, and AP-1. Subsequently, the transcriptional regulators enter the nucleus and bind to regulatory elements of the IFNβ promoter. After various modifications, the position of the nucleosome is altered and the complex is assembled to activate the IFNβ expression. However, IFNβ regulation involves a complex network. For the study of immunity and diseases, it is important to understand how transcription factors bind to regulatory elements through specific forms, which elements in cells are involved in regulation, what regulation occurs during the assembly of enhancers and transcription complexes, and the possible regulatory mechanisms after transcription. Thus, this review focuses on the various regulatory mechanisms and elements involved in the activation of IFNβ expression. In addition, we discuss the impact of this regulation in biology.

Functional analysis of a random-sequence chromosome reveals a high level and the molecular nature of transcriptional noise in yeast cells

We measure transcriptional noise in yeast by analyzing chromatin structure and transcription of an 18-kb region of DNA whose sequence was randomly generated. Nucleosomes fully occupy random-sequence DNA, but nucleosome-depleted regions (NDRs) are much less frequent, and there are fewer well-positioned nucleosomes and shorter nucleosome arrays. Steady-state levels of random-sequence RNAs are comparable to yeast mRNAs, although transcription and decay rates are higher. Transcriptional initiation from random-sequence DNA occurs at numerous sites, indicating very low intrinsic specificity of the RNA Pol II machinery. In contrast, poly(A) profiles of random-sequence RNAs are roughly comparable to those of yeast mRNAs, suggesting limited evolutionary restraints on poly(A) site choice. Random-sequence RNAs show higher cell-to-cell variability than yeast mRNAs, suggesting that functional elements limit variability. These observations indicate that transcriptional noise occurs at high levels in yeast, and they provide insight into how chromatin and transcription patterns arise from the evolved yeast genome.

Swc4 protects nucleosome-free rDNA, tDNA and telomere loci to inhibit genome instability

In the baker's yeast Saccharomyces cerevisiae, NuA4 and SWR1-C, two multisubunit complexes, are involved in histone acetylation and chromatin remodeling, respectively. Eaf1 is the assembly platform subunit of NuA4, Swr1 is the assembly platform and catalytic subunit of SWR1-C, while Swc4, Yaf9, Arp4 and Act1 form a functional module, and is present in both NuA4 and SWR1 complexes. ACT1 and ARP4 are essential for cell survival. Deletion of SWC4, but not YAF9, EAF1 or SWR1 results in a severe growth defect, but the underlying mechanism remains largely unknown. Here, we show that swc4Δ, but not yaf9Δ, eaf1Δ, or swr1Δ cells display defects in DNA ploidy and chromosome segregation, suggesting that the defects observed in swc4Δ cells are independent of NuA4 or SWR1-C integrity. Swc4 is enriched in the nucleosome-free regions (NFRs) of the genome, including characteristic regions of RDN5s, tDNAs and telomeres, independently of Yaf9, Eaf1 or Swr1. In particular, rDNA, tDNA and telomere loci are more unstable and prone to recombination in the swc4Δ cells than in wild-type cells. Taken together, we conclude that the chromatin associated Swc4 protects nucleosome-free chromatin of rDNA, tDNA and telomere loci to ensure genome integrity.

Predicting nucleosome positioning using statistical equilibrium models in budding yeast

We present a protocol using thermodynamic models to predict nucleosome positioning with transcription factors (TFs) and chromatin remodelers. We describe step-by-step approaches to annotate genome-wide nucleosome-depleted regions (NDRs), compute nucleosome and TF occupancy, optimize parameters, and evaluate model performance. These models identify nucleosome-displacing TFs in the budding yeast genome and predict the locations and sizes of NDRs solely based on DNA sequence and TF motifs. The protocol can be applied to all organisms with prior knowledge of TF motifs. For complete details on the use and execution of this protocol, please refer to Kharerin and Bai (2021).1.

Histone H2A monoubiquitination marks are targeted to specific sites by cohesin subunits in Arabidopsis

Histone H2A monoubiquitination (H2Aub1) functions as a conserved posttranslational modification in eukaryotes to maintain gene expression and guarantee cellular identity. Arabidopsis H2Aub1 is catalyzed by the core components AtRING1s and AtBMI1s of polycomb repressive complex 1 (PRC1). Because PRC1 components lack known DNA binding domains, it is unclear how H2Aub1 is established at specific genomic locations. Here, we show that the Arabidopsis cohesin subunits AtSYN4 and AtSCC3 interact with each other, and AtSCC3 binds to AtBMI1s. H2Aub1 levels are reduced in atsyn4 mutant or AtSCC3 artificial microRNA knockdown plants. ChIP-seq assays indicate that most binding events of AtSYN4 and AtSCC3 are associated with H2Aub1 along the genome where transcription is activated independently of H3K27me3. Finally, we show that AtSYN4 binds directly to the G-box motif and directs H2Aub1 to these sites. Our study thus reveals a mechanism for cohesin-mediated recruitment of AtBMI1s to specific genomic loci to mediate H2Aub1.

Beyond Sequence: Internucleosomal Interactions Dominate Array Assembly

The organization of the nucleosome array is a critical component of the chromatin assembly into higher order structure as well as its function. Here, we investigated the contributions of the DNA sequence and internucleosomal interactions on the organization of the nucleosomal arrays in compact structures using atomic force microscopy. We assembled nucleosomes on DNA substrates allowing for the formation of tetranucleosomes. We found that nucleosomes are capable of close positioning with no discernible space between them, even in the case of assembled dinucleosomes. This morphology of the array is in contrast with that observed for arrays assembled with repeats of the nucleosome positioning motifs separated by uniform spacers. Simulated assembly of tetranucleosomes by random placement along the substrates revealed that nucleosome array compaction is promoted by the interaction of the nucleosomes. We developed a theoretical model to account for the role of DNA sequence and internucleosomal interactions in the formation of the nucleosome structures. These findings suggest that, in the chromatin assembly, the affinity of the nucleosomes to the DNA sequence and the strengths of the internucleosomal interactions are the two major factors defining the compactness of the chromatin.

New insights into genome annotation in Podospora anserina through re-exploiting multiple RNA-seq data

BACKGROUND

Publicly available RNA-seq datasets are often underused although being helpful to improve functional annotation of eukaryotic genomes. This is especially true for filamentous fungi genomes which structure differs from most well annotated yeast genomes. Podospora anserina is a filamentous fungal model, which genome has been sequenced and annotated in 2008. Still, the current annotation lacks information about cis-regulatory elements, including promoters, transcription starting sites and terminators, which are instrumental to integrate epigenomic features into global gene regulation strategies.

RESULTS

Here we took advantage of 37 RNA-seq experiments that were obtained in contrasted developmental and physiological conditions, to complete the functional annotation of P. anserina genome. Out of the 10,800 previously annotated genes, 5'UTR and 3'UTR were defined for 7554, among which, 3328 showed differential transcriptional signal starts and/or transcriptional end sites. In addition, alternative splicing events were detected for 2350 genes, mostly due alternative 3'splice sites and 1732 novel transcriptionally active regions (nTARs) in unannotated regions were identified.

CONCLUSIONS

Our study provides a comprehensive genome-wide functional annotation of P. anserina genome, including chromatin features, cis-acting elements such as UTRs, alternative splicing events and transcription of non-coding regions. These new findings will likely improve our understanding of gene regulation strategies in compact genomes, such as those of filamentous fungi. Characterization of alternative transcripts and nTARs paves the way to the discovery of putative new genes, alternative peptides or regulatory non-coding RNAs.

Alternative transcription start sites contribute to acute-stress-induced transcriptome response in human skeletal muscle

Background

More than half of human protein-coding genes have an alternative transcription start site (TSS). We aimed to investigate the contribution of alternative TSSs to the acute-stress-induced transcriptome response in human tissue (skeletal muscle) using the cap analysis of gene expression approach. TSSs were examined at baseline and during recovery after acute stress (a cycling exercise).

Results

We identified 44,680 CAGE TSS clusters (including 3764 first defined) belonging to 12,268 genes and annotated for the first time 290 TSSs belonging to 163 genes. The transcriptome dynamically changes during the first hours after acute stress; the change in the expression of 10% of genes was associated with the activation of alternative TSSs, indicating differential TSSs usage. The majority of the alternative TSSs do not increase proteome complexity suggesting that the function of thousands of alternative TSSs is associated with the fine regulation of mRNA isoform expression from a gene due to the transcription factor-specific activation of various alternative TSSs. We identified individual muscle promoter regions for each TSS using muscle open chromatin data (ATAC-seq and DNase-seq). Then, using the positional weight matrix approach we predicted time course activation of “classic” transcription factors involved in response of skeletal muscle to contractile activity, as well as diversity of less/un-investigated factors.

Conclusions

Transcriptome response induced by acute stress related to activation of the alternative TSSs indicates that differential TSSs usage is an essential mechanism of fine regulation of gene response to stress stimulus. A comprehensive resource of accurate TSSs and individual promoter regions for each TSS in muscle was created. This resource together with the positional weight matrix approach can be used to accurate prediction of TFs in any gene(s) of interest involved in the response to various stimuli, interventions or pathological conditions in human skeletal muscle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00399-8.

ATP binding facilitates target search of SWR1 chromatin remodeler by promoting one-dimensional diffusion on DNA

One-dimensional (1D) target search is a well-characterized phenomenon for many DNA-binding proteins but is poorly understood for chromatin remodelers. Herein, we characterize the 1D scanning properties of SWR1, a conserved yeast chromatin remodeler that performs histone exchange on +1 nucleosomes adjacent to a nucleosome-depleted region (NDR) at gene promoters. We demonstrate that SWR1 has a kinetic binding preference for DNA of NDR length as opposed to gene-body linker length DNA. Using single and dual color single-particle tracking on DNA stretched with optical tweezers, we directly observe SWR1 diffusion on DNA. We found that various factors impact SWR1 scanning, including ATP which promotes diffusion through nucleotide binding rather than ATP hydrolysis. A DNA-binding subunit, Swc2, plays an important role in the overall diffusive behavior of the complex, as the subunit in isolation retains similar, although faster, scanning properties as the whole remodeler. ATP-bound SWR1 slides until it encounters a protein roadblock, of which we tested dCas9 and nucleosomes. The median diffusion coefficient, 0.024 μm2/s, in the regime of helical sliding, would mediate rapid encounter of NDR-flanking nucleosomes at length scales found in cellular chromatin.

Global proteomic analysis of pre-implantational ovine embryos produced in-vitro

The present study was conducted to characterize the major proteome of pre-implantation (D6) ovine embryos produced in vitro. COCs were aspirated from antral follicles (2-6 mm), matured and fertilized in vitro, and cultured until day six. Proteins were extracted separately from three pools of 45 embryos and separately run in SDS-PAGE. Proteins from each pool were individually subjected to in-gel digestion followed by LC-MS/MS. Three "raw. files" and protein lists were produced by Pattern Lab software but only proteins present in all three lists were used for the bioinformatics analyses. https://david.ncifcrf.govThere were 2,262 proteins identified in the 6-day old ovine embryos, including albumin, zona pellucida glycoprotein 2, 3 and 4, peptidyl arginine deiminase 6, actin cytoplasmic 1, gamma-actin 1, pyruvate kinase, heat shock protein 90 and protein disulfide isomerase, among others. Major biological processes linked to the sheep embryo proteome were translation, protein transport and protein stabilization, and molecular functions, defined as ATP binding, oxygen carrier activity and oxygen binding. There were 42 enriched functional clusters according to the 2,147 genes (UniProt database). Ten selected clusters with potential association with embryo development included translation, structural constituent of ribosomes, ribosomes, nucleosomes, structural constituent of the cytoskeleton, microtubule-based process, translation initiation factor activity, regulation of translational initiation, cell body and nucleotide biosynthetic process. The most representative KEEG pathways were ribosome, oxidative phosphorylation, glutathione metabolism, gap junction, mineral absorption, DNA replication and cGMP-PKG signaling pathway. Analyses of functional clusters clearly showed differences associated with the proteome of pre-implantation (D6) sheep embryos generated after in vitro fertilization in comparison with in vivo counterparts (Sanchez et al., 2021; https://doi.org/10.1111/rda.13897), confirming that the quality of in vitro derived blastocysts are unlike those produced in vivo. The present study portrays the first comprehensive overview of the proteome of pre-implantational ovine embryos grown in vitro.

GC content strongly influences the role of poly(dA) in the intrinsic nucleosome positioning in Saccharomyces cerevisiae

The nucleosome is the basic structural element of genomic DNA packaging and plays a role in transcription, replication, and recombination. Poly(dA) tracts are considered major sequence determinants of nucleosome positioning, although their role is not well understood. Here, we show that the homopolymeric character and the low GC content of poly(dA)s play different roles in nucleosome formation. We found that the inherent low GC content of poly(dA) alone can account for the deep and anisotropic nucleosome depletion at structurally and functionally important regions of promoters and origins of replication. We also show that the level of nucleosome occupancy at poly(dA) is strongly related to the local nucleotide background and its high frequency of occurrence in Saccharomyces cerevisiae does not appear merely to be associated with its intrinsic nucleosome-excluding properties. In addition, we show that the GC content alone can predict more than 60% of the in vitro nucleosome map, providing further evidence that the intrinsic nucleosome positioning is more greatly determined by GC content than poly(dA) stretches. Our results are consistent with a model in which poly(dA) stretches act at two distinct levels: first, by its low GC content, which intrinsically contributes to hinder nucleosome formation, and second, by its contiguous runs of dA that selectively drive the recruitment of non-histone proteins with structural and functional roles.

Transposon insertional mutagenesis of diverse yeast strains suggests coordinated gene essentiality polymorphisms

Due to epistasis, the same mutation can have drastically different phenotypic consequences in different individuals. This phenomenon is pertinent to precision medicine as well as antimicrobial drug development, but its general characteristics are largely unknown. We approach this question by genome-wide assessment of gene essentiality polymorphism in 16 Saccharomyces cerevisiae strains using transposon insertional mutagenesis. Essentiality polymorphism is observed for 9.8% of genes, most of which have had repeated essentiality switches in evolution. Genes exhibiting essentiality polymorphism lean toward having intermediate numbers of genetic and protein interactions. Gene essentiality changes tend to occur concordantly among components of the same protein complex or metabolic pathway and among a group of over 100 mitochondrial proteins, revealing molecular machines or functional modules as units of gene essentiality variation. Most essential genes tolerate transposon insertions consistently among strains in one or more coding segments, delineating nonessential regions within essential genes.

Whole-genome methods to define DNA and histone accessibility and long-range interactions in chromatin

Defining the genome-wide chromatin landscape has been a goal of experimentalists for decades. Here we review highlights of these efforts, from seminal experiments showing discontinuities in chromatin structure related to gene activation to extensions of these methods elucidating general features of chromatin related to gene states by exploiting deep sequencing methods. We also review chromatin conformational capture methods to identify patterns in long-range interactions between genomic loci.

Chemically Induced Chromosomal Interaction (CICI) method to study chromosome dynamics and its biological roles

Numerous intra- and inter-chromosomal contacts have been mapped in eukaryotic genomes, but it remains challenging to link these 3D structures to their regulatory functions. To establish the causal relationships between chromosome conformation and genome functions, we  develop a method, Chemically Induced Chromosomal Interaction (CICI), to selectively perturb the chromosome conformation at targeted loci. Using this method, long-distance chromosomal interactions can be induced dynamically between intra- or inter-chromosomal loci pairs, including the ones with very low Hi-C contact frequencies. Measurement of CICI formation time allows us to probe chromosome encounter dynamics between different loci pairs across the cell cycle. We also conduct two functional tests of CICI. We perturb the chromosome conformation near a DNA double-strand break and observe altered donor preference in homologous recombination; we force interactions between early and late-firing DNA replication origins and find no significant changes in replication timing. These results suggest that chromosome conformation plays a deterministic role in homology-directed DNA repair, but not in the establishment of replication timing. Overall, our study demonstrates that CICI is a powerful tool to study chromosome dynamics and 3D genome function.

Methods to selectively manipulate specific long-distance chromosomal interactions are limited. Here the authors develop a method called Chemically Induced Chromosomal Interaction (CICI) to engineer interactions and demonstrate that 3D conformation plays a causal role in establishing donor DNA preference during DNA repair.

Control of Gene Expression via the Yeast CWI Pathway

Living cells exposed to stressful environmental situations can elicit cellular responses that guarantee maximal cell survival. Most of these responses are mediated by mitogen-activated protein kinase (MAPK) cascades, which are highly conserved from yeast to humans. Cell wall damage conditions in the yeast Saccharomyces cerevisiae elicit rescue mechanisms mainly associated with reprogramming specific transcriptional responses via the cell wall integrity (CWI) pathway. Regulation of gene expression by this pathway is coordinated by the MAPK Slt2/Mpk1, mainly via Rlm1 and, to a lesser extent, through SBF (Swi4/Swi6) transcription factors. In this review, we summarize the molecular mechanisms controlling gene expression upon cell wall stress and the role of chromatin structure in these processes. Some of these mechanisms are also discussed in the context of other stresses governed by different yeast MAPK pathways. Slt2 regulates both transcriptional initiation and elongation by interacting with chromatin at the promoter and coding regions of CWI-responsive genes but using different mechanisms for Rlm1- and SBF-dependent genes. Since MAPK pathways are very well conserved in eukaryotic cells and are essential for controlling cellular physiology, improving our knowledge regarding how they regulate gene expression could impact the future identification of novel targets for therapeutic intervention.

Mapping of Nonhomologous End Joining-Mediated Integration Facilitates Genome-Scale Trackable Mutagenesis in Yarrowia lipolytica

Genome-scale mutagenesis, phenotypic screening, and tracking the causal mutations is a powerful approach for genetic analysis. However, classic mutagenesis approaches require extensive effort to identify causal mutations. It is desirable to demonstrate a powerful approach for rapid trackable mutagenesis. Here, we mapped the distribution of nonhomologous end joining (NHEJ)-mediated integration for the first time and demonstrated that it can be used for constructing the genome-scale trackable mutagenesis library in Yarrowia lipolytica. The sequencing of 9.15 × 10⁵ insertions showed that NHEJ-mediated integration inserted DNA randomly across the chromosomes, and the transcriptional regulatory regions exhibited integration preference. The insertions were located in both nucleosome-occupancy regions and nucleosome-free regions. Using NHEJ-mediated integration to construct the genome-scale mutagenesis library, the new targets that improved β-carotene biosynthesis and acetic acid tolerance were identified rapidly. This mutagenesis approach is readily applicable to other organisms with strong NHEJ preference and will contribute to cell factory construction.

Predicting the configuration and energy of DNA in a nucleosome by coarse-grain modelling

We present a novel algorithm which uses a coarse-grained model and an energy minimisation procedure to predict the sequence-dependent DNA configuration in a nucleosome together with its energetic cost.

Reb1, Cbf1, and Pho4 bias histone sliding and deposition away from their binding sites

In transcriptionally active genes, nucleosome positions in promoters are regulated by nucleosome displacing factors (NDFs) and chromatin remodeling enzymes. Depletion of NDFs or the RSC chromatin remodeler shrinks or abolishes the nucleosome depleted regions (NDRs) in promoters, which can suppress gene activation and result in cryptic transcription. Despite their vital cellular functions, how the action of chromatin remodelers may be directly affected by site-specific binding factors like NDFs is poorly understood. Here we demonstrate that two NDFs, Reb1 and Cbf1, can direct both Chd1 and RSC chromatin remodeling enzymes in vitro, stimulating repositioning of the histone core away from their binding sites. Interestingly, although the Pho4 transcription factor had a much weaker effect on nucleosome positioning, both NDFs and Pho4 were able to similarly redirect positioning of hexasomes. In chaperone-mediated nucleosome assembly assays, Reb1 but not Pho4 showed an ability to block deposition of the histone H3/H4 tetramer, but Reb1 did not block addition of the H2A/H2B dimer to hexasomes. Our in vitro results show that NDFs bias the action of remodelers to increase the length of the free DNA in the vicinity of their binding sites. These results suggest that NDFs could directly affect NDR architecture through chromatin remodelers.

Pioneer factors in development and cancer

Transcription factors (TFs) are essential mediators of epigenetic regulation and modifiers of penetrance. Studies from the past decades have revealed a sub-class of TF that is capable of remodeling closed chromatin states through targeting nucleosomal motifs. This pioneer factor (PF) class of chromatin remodeler is ATP independent in its roles in epigenetic initiation, with nucleosome-motif recognition and association with repressive chromatin regions. Increasing evidence suggests that the fundamental properties of PFs can be coopted in human cancers. We explore the role of PFs in the larger context of tissue-specific epigenetic regulation. Moreover, we highlight an emerging class of chimeric PF derived from translocation partners in human disease and PFs associated with rare tumors. In the age of site-directed genome editing and targeted protein degradation, increasing our understanding of PFs will provide access to next-generation therapy for human disease driven from altered transcriptional circuitry.

Proximal-end bias from in-vitro reconstituted nucleosomes and the result on downstream data analysis

The most basic level of eukaryotic gene regulation is the presence or absence of nucleosomes on DNA regulatory elements. In an effort to elucidate in vivo nucleosome patterns, in vitro studies are frequently used. In vitro, short DNA fragments are more favorable for nucleosome formation, increasing the likelihood of nucleosome occupancy. This may in part result from the fact that nucleosomes prefer to form on the terminal ends of linear DNA. This phenomenon has the potential to bias in vitro reconstituted nucleosomes and skew results. If the ends of DNA fragments are known, the reads falling close to the ends are typically discarded. In this study we confirm the phenomenon of end bias of in vitro nucleosomes. We describe a method in which nearly identical libraries, with different known ends, are used to recover nucleosomes which form towards the terminal ends of fragmented DNA. Finally, we illustrate that although nucleosomes prefer to form on DNA ends, it does not appear to skew results or the interpretation thereof.

Linker histone H1.8 inhibits chromatin binding of condensins and DNA topoisomerase II to tune chromosome length and individualization

DNA loop extrusion by condensins and decatenation by DNA topoisomerase II (topo II) are thought to drive mitotic chromosome compaction and individualization. Here, we reveal that the linker histone H1.8 antagonizes condensins and topo II to shape mitotic chromosome organization. In vitro chromatin reconstitution experiments demonstrate that H1.8 inhibits binding of condensins and topo II to nucleosome arrays. Accordingly, H1.8 depletion in Xenopus egg extracts increased condensins and topo II levels on mitotic chromatin. Chromosome morphology and Hi-C analyses suggest that H1.8 depletion makes chromosomes thinner and longer through shortening the average loop size and reducing the DNA amount in each layer of mitotic loops. Furthermore, excess loading of condensins and topo II to chromosomes by H1.8 depletion causes hyper-chromosome individualization and dispersion. We propose that condensins and topo II are essential for chromosome individualization, but their functions are tuned by the linker histone to keep chromosomes together until anaphase.

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.

RoboCOP: jointly computing chromatin occupancy profiles for numerous factors from chromatin accessibility data

Chromatin is a tightly packaged structure of DNA and protein within the nucleus of a cell. The arrangement of different protein complexes along the DNA modulates and is modulated by gene expression. Measuring the binding locations and occupancy levels of different transcription factors (TFs) and nucleosomes is therefore crucial to understanding gene regulation. Antibody-based methods for assaying chromatin occupancy are capable of identifying the binding sites of specific DNA binding factors, but only one factor at a time. In contrast, epigenomic accessibility data like MNase-seq, DNase-seq, and ATAC-seq provide insight into the chromatin landscape of all factors bound along the genome, but with little insight into the identities of those factors. Here, we present RoboCOP, a multivariate state space model that integrates chromatin accessibility data with nucleotide sequence to jointly compute genome-wide probabilistic scores of nucleosome and TF occupancy, for hundreds of different factors. We apply RoboCOP to MNase-seq and ATAC-seq data to elucidate the protein-binding landscape of nucleosomes and 150 TFs across the yeast genome, and show that our model makes better predictions than existing methods. We also compute a chromatin occupancy profile of the yeast genome under cadmium stress, revealing chromatin dynamics associated with transcriptional regulation.

Fast and inexpensive whole-genome sequencing library preparation from intact yeast cells

Through the increase in the capacity of sequencing machines massively parallel sequencing of thousands of samples in a single run is now possible. With the improved throughput and resulting drop in the price of sequencing, the cost and time for preparation of sequencing libraries have become the major bottleneck in large-scale experiments. Methods using a hyperactive variant of the Tn5 transposase efficiently generate libraries starting from cDNA or genomic DNA in a few hours and are highly scalable. For genome sequencing, however, the time and effort spent on genomic DNA isolation limit the practicability of sequencing large numbers of samples. Here, we describe a highly scalable method for preparing high-quality whole-genome sequencing libraries directly from Saccharomyces cerevisiae cultures in less than 3 h at 34 cents per sample. We skip the rate-limiting step of genomic DNA extraction by directly tagmenting lysed yeast spheroplasts and add a nucleosome release step prior to enrichment PCR to improve the evenness of genomic coverage. Resulting libraries do not show any GC bias and are comparable in quality to libraries processed from genomic DNA with a commercially available Tn5-based kit. We use our protocol to investigate CRISPR/Cas9 on- and off-target edits and reliably detect edited variants and shared polymorphisms between strains. Our protocol enables rapid preparation of unbiased and high-quality, sequencing-ready indexed libraries for hundreds of yeast strains in a single day at a low price. By adjusting individual steps of our workflow, we expect that our protocol can be adapted to other organisms.

Transcription at a Distance in the Budding Yeast Saccharomyces cerevisiae

Proper transcriptional regulation depends on the collaboration of multiple layers of control simultaneously. Cells tightly balance cellular resources and integrate various signaling inputs to maintain homeostasis during growth, development and stressors, among other signals. Many eukaryotes, including the budding yeast Saccharomyces cerevisiae, exhibit a non-random distribution of functionally related genes throughout their genomes. This arrangement coordinates the transcription of genes that are found in clusters, and can occur over long distances. In this work, we review the current literature pertaining to gene regulation at a distance in budding yeast.

Comparative analysis and prediction of nucleosome positioning using integrative feature representation and machine learning algorithms

Background

Nucleosome plays an important role in the process of genome expression, DNA replication, DNA repair and transcription. Therefore, the research of nucleosome positioning has invariably received extensive attention. Considering the diversity of DNA sequence representation methods, we tried to integrate multiple features to analyze its effect in the process of nucleosome positioning analysis. This process can also deepen our understanding of the theoretical analysis of nucleosome positioning.

Results

Here, we not only used frequency chaos game representation (FCGR) to construct DNA sequence features, but also integrated it with other features and adopted the principal component analysis (PCA) algorithm. Simultaneously, support vector machine (SVM), extreme learning machine (ELM), extreme gradient boosting (XGBoost), multilayer perceptron (MLP) and convolutional neural networks (CNN) are used as predictors for nucleosome positioning prediction analysis, respectively. The integrated feature vector prediction quality is significantly superior to a single feature. After using principal component analysis (PCA) to reduce the feature dimension, the prediction quality of H. sapiens dataset has been significantly improved.

Conclusions

Comparative analysis and prediction on H. sapiens, C. elegans, D. melanogaster and S. cerevisiae datasets, demonstrate that the application of FCGR to nucleosome positioning is feasible, and we also found that integrative feature representation would be better.

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

The nucleosome is a major modulator of DNA accessibility to other cellular factors. Nucleosome positioning has a critical importance in regulating cell processes such as transcription, replication, recombination or DNA repair. The DNA sequence has an influence on the position of nucleosomes on genomes, although other factors are also implicated, such as ATP-dependent remodelers or competition of the nucleosome with DNA binding proteins. Different sequence motifs can promote or inhibit the nucleosome formation, thus influencing the accessibility to the DNA. Sequence-encoded nucleosome positioning having functional consequences on cell processes can then be selected or counter-selected during evolution. We review the interplay between sequence evolution and nucleosome positioning evolution. We first focus on the different ways to encode nucleosome positions in the DNA sequence, and to which extent these mechanisms are responsible of genome-wide nucleosome positioning in vivo. Then, we discuss the findings about selection of sequences for their nucleosomal properties. Finally, we illustrate how the nucleosome can directly influence sequence evolution through its interactions with DNA damage and repair mechanisms. This review aims to provide an overview of the mutual influence of sequence evolution and nucleosome positioning evolution, possibly leading to complex evolutionary dynamics.

Basis of specificity for a conserved and promiscuous chromatin remodeling protein

Eukaryotic genomes are organized dynamically through the repositioning of nucleosomes. Isw2 is an enzyme that has been previously defined as a genome-wide, nonspecific nucleosome spacing factor. Here, we show that Isw2 instead acts as an obligately targeted nucleosome remodeler in vivo through physical interactions with sequence-specific factors. We demonstrate that Isw2-recruiting factors use small and previously uncharacterized epitopes, which direct Isw2 activity through highly conserved acidic residues in the Isw2 accessory protein Itc1. This interaction orients Isw2 on target nucleosomes, allowing for precise nucleosome positioning at targeted loci. Finally, we show that these critical acidic residues have been lost in the Drosophila lineage, potentially explaining the inconsistently characterized function of Isw2-like proteins. Altogether, these data suggest an 'interacting barrier model,' where Isw2 interacts with a sequence-specific factor to accurately and reproducibly position a single, targeted nucleosome to define the precise border of phased chromatin arrays.

TATA and paused promoters active in differentiated tissues have distinct expression characteristics

Core promoter types differ in the extent to which RNA polymerase II (Pol II) pauses after initiation, but how this affects their tissue-specific gene expression characteristics is not well understood. While promoters with Pol II pausing elements are active throughout development, TATA promoters are highly active in differentiated tissues. We therefore used a genomics approach on late-stage Drosophila embryos to analyze the properties of promoter types. Using tissue-specific Pol II ChIP-seq, we found that paused promoters have high levels of paused Pol II throughout the embryo, even in tissues where the gene is not expressed, while TATA promoters only show Pol II occupancy when the gene is active. The promoter types are associated with different chromatin accessibility in ATAC-seq data and have different expression characteristics in single-cell RNA-seq data. The two promoter types may therefore be optimized for different properties: paused promoters show more consistent expression when active, while TATA promoters have lower background expression when inactive. We propose that tissue-specific genes have evolved to use two different strategies for their differential expression across tissues.

High-Resolution Mapping of Transcription Initiation in the Asexual Stages of Toxoplasma gondii

Toxoplasma gondii is a common parasite of humans and animals, causing life-threatening disease in the immunocompromized, fetal abnormalities when contracted during gestation, and recurrent ocular lesions in some patients. Central to the prevalence and pathogenicity of this protozoan is its ability to adapt to a broad range of environments, and to differentiate between acute and chronic stages. These processes are underpinned by a major rewiring of gene expression, yet the mechanisms that regulate transcription in this parasite are only partially characterized. Deciphering these mechanisms requires a precise and comprehensive map of transcription start sites (TSSs); however, Toxoplasma TSSs have remained incompletely defined. To address this challenge, we used 5'-end RNA sequencing to genomically assess transcription initiation in both acute and chronic stages of Toxoplasma. Here, we report an in-depth analysis of transcription initiation at promoters, and provide empirically-defined TSSs for 7603 (91%) protein-coding genes, of which only 1840 concur with existing gene models. Comparing data from acute and chronic stages, we identified instances of stage-specific alternative TSSs that putatively generate mRNA isoforms with distinct 5' termini. Analysis of the nucleotide content and nucleosome occupancy around TSSs allowed us to examine the determinants of TSS choice, and outline features of Toxoplasma promoter architecture. We also found pervasive divergent transcription at Toxoplasma promoters, clustered within the nucleosomes of highly-symmetrical phased arrays, underscoring chromatin contributions to transcription initiation. Corroborating previous observations, we asserted that Toxoplasma 5' leaders are among the longest of any eukaryote studied thus far, displaying a median length of approximately 800 nucleotides. Further highlighting the utility of a precise TSS map, we pinpointed motifs associated with transcription initiation, including the binding sites of the master regulator of chronic-stage differentiation, BFD1, and a novel motif with a similar positional arrangement present at 44% of Toxoplasma promoters. This work provides a critical resource for functional genomics in Toxoplasma, and lays down a foundation to study the interactions between genomic sequences and the regulatory factors that control transcription in this parasite.

Thermodynamic modeling of genome-wide nucleosome depleted regions in yeast

Nucleosome positioning in the genome is essential for the regulation of many nuclear processes. We currently have limited capability to predict nucleosome positioning in vivo, especially the locations and sizes of nucleosome depleted regions (NDRs). Here, we present a thermodynamic model that incorporates the intrinsic affinity of histones, competitive binding of sequence-specific factors, and nucleosome remodeling to predict nucleosome positioning in budding yeast. The model shows that the intrinsic affinity of histones, at near-saturating histone concentration, is not sufficient in generating NDRs in the genome. However, the binding of a few factors, especially RSC towards GC-rich and poly(A/T) sequences, allows us to predict ~ 66% of genome-wide NDRs. The model also shows that nucleosome remodeling activity is required to predict the correct NDR sizes. The validity of the model was further supported by the agreement between the predicted and the measured nucleosome positioning upon factor deletion or on exogenous sequences introduced into yeast. Overall, our model quantitatively evaluated the impact of different genetic components on NDR formation and illustrated the vital roles of sequence-specific factors and nucleosome remodeling in this process.

Nucleosome is the basic unit of chromatin, containing 147 base-pairs of DNA wrapped around a histone core. The positioning of nucleosomes, i.e., which parts of DNA are inside nucleosome and which parts are nucleosome-free, is highly regulated. In particular, regulatory sequences tend to be exposed in nucleosome-depleted regions (NDRs), and such exposure is crucial for a variety of processes including DNA replication, repair, and gene expression. Here, we used a thermodynamics model to predict nucleosome positioning on the yeast genome. The model shows that the intrinsic sequence preference of histones is not sufficient in generating NDRs. In contrast, binding of a few transcription factors, especially RSC, is largely responsible for NDR formation. Nucleosome remodeling activity is also required in the model to recapitulate the NDR sizes. This model contributes to our understanding of the mechanisms that regulate nucleosome positioning. It can also be used to predict nucleosome positioning in mutant yeast or on novel DNA sequences.

Measuring DNA mechanics on the genome scale

Mechanical deformations of DNA such as bending are ubiquitous and have been implicated in diverse cellular functions1. However, the lack of high-throughput tools to measure the mechanical properties of DNA has limited our understanding of how DNA mechanics influence chromatin transactions across the genome. Here we develop 'loop-seq'-a high-throughput assay to measure the propensity for DNA looping-and determine the intrinsic cyclizabilities of 270,806 50-base-pair DNA fragments that span Saccharomyces cerevisiae chromosome V, other genomic regions, and random sequences. We found sequence-encoded regions of unusually low bendability within nucleosome-depleted regions upstream of transcription start sites (TSSs). Low bendability of linker DNA inhibits nucleosome sliding into the linker by the chromatin remodeller INO80, which explains how INO80 can define nucleosome-depleted regions in the absence of other factors2. Chromosome-wide, nucleosomes were characterized by high DNA bendability near dyads and low bendability near linkers. This contrast increases for deeper gene-body nucleosomes but disappears after random substitution of synonymous codons, which suggests that the evolution of codon choice has been influenced by DNA mechanics around gene-body nucleosomes. Furthermore, we show that local DNA mechanics affect transcription through TSS-proximal nucleosomes. Overall, this genome-scale map of DNA mechanics indicates a 'mechanical code' with broad functional implications.

Nucleosome Positioning around Transcription Start Site Correlates with Gene Expression Only for Active Chromatin State in Drosophila Interphase Chromosomes

We analyzed the whole-genome experimental maps of nucleosomes in Drosophila melanogaster and classified genes by the expression level in S2 cells (RPKM value, reads per kilobase million) as well as the number of tissues in which a gene was expressed (breadth of expression, BoE). Chromatin in 5′-regions of genes we classified on four states according to the hidden Markov model (4HMM). Only the Aquamarine chromatin state we considered as Active, while the rest three states we defined as Non-Active. Surprisingly, about 20/40% of genes with 5′-regions mapped to Active/Non-Active chromatin possessed the minimal/at least modest RPKM and BoE. We found that regardless of RPKM/BoE the genes of Active chromatin possessed the regular nucleosome arrangement in 5′-regions, while genes of Non-Active chromatin did not show respective specificity. Only for genes of Active chromatin the RPKM/BoE positively correlates with the number of nucleosome sites upstream/around TSS and negatively with that downstream TSS. We propose that for genes of Active chromatin, regardless of RPKM value and BoE the nucleosome arrangement in 5′-regions potentiates transcription, while for genes of Non-Active chromatin, the transcription machinery does not require the substantial support from nucleosome arrangement to influence gene expression.

Transcription-dependent enrichment of the yeast FACT complex influences nucleosome dynamics on the RNA polymerase III-transcribed genes

The FACT (FAcilitates Chromatin Transactions) complex influences transcription initiation and enables passage of RNA polymerase (pol) II through gene body nucleosomes during elongation. In the budding yeast, ~280 non-coding RNA genes highly transcribed in vivo by pol III are found in the nucleosome-free regions bordered by positioned nucleosomes. The downstream nucleosome dynamics was found to regulate transcription via controlling the gene terminator accessibility and hence, terminator-dependent pol III recycling. As opposed to the enrichment at the 5'-ends of pol II-transcribed genes, our genome-wide mapping found transcription-dependent enrichment of the FACT subunit Spt16 near the 3'-end of all pol III-transcribed genes. Spt16 physically associates with the pol III transcription complex and shows gene-specific occupancy levels on the individual genes. On the non-tRNA pol III-transcribed genes, Spt16 facilitates transcription by reducing the nucleosome occupany on the gene body. On the tRNA genes, it maintains the position of the nucleosome at the 3' gene-end and affects transcription in gene-specific manner. Under nutritional stress, Spt16 enrichment is abolished in the gene downstream region of all pol III-transcribed genes and reciprocally changed on the induced or repressed pol II-transcribed ESR genes. Under the heat and replicative stress, its occupancy on the pol III-transcribed genes increases significantly. Our results show that Spt16 elicits a differential, gene-specific and stress-responsive dynamics, which provides a novel stress-sensor mechanism of regulating transcription against external stress. By primarily influencing the nucleosomal organization, FACT links the downstream nucleosome dynamics to transcription and environmental stress on the pol III-transcribed genes.

FACT and Ash1 promote long-range and bidirectional nucleosome eviction at the HO promoter

The Saccharomyces cerevisiae HO gene is a model regulatory system with complex transcriptional regulation. Budding yeast divide asymmetrically and HO is expressed only in mother cells where a nucleosome eviction cascade along the promoter during the cell cycle enables activation. HO expression in daughter cells is inhibited by high concentration of Ash1 in daughters. To understand how Ash1 represses transcription, we used a myo4 mutation which boosts Ash1 accumulation in both mothers and daughters and show that Ash1 inhibits promoter recruitment of SWI/SNF and Gcn5. We show Ash1 is also required for the efficient nucleosome repopulation that occurs after eviction, and the strongest effects of Ash1 are seen when Ash1 has been degraded and at promoter locations distant from where Ash1 bound. Additionally, we defined a specific nucleosome/nucleosome-depleted region structure that restricts HO activation to one of two paralogous DNA-binding factors. We also show that nucleosome eviction occurs bidirectionally over a large distance. Significantly, eviction of the more distant nucleosomes is dependent upon the FACT histone chaperone, and FACT is recruited to these regions when eviction is beginning. These last observations, along with ChIP experiments involving the SBF factor, suggest a long-distance loop transiently forms at the HO promoter.

Promoter Architecture and Promoter Engineering in Saccharomyces cerevisiae

Promoters play an essential role in the regulation of gene expression for fine-tuning genetic circuits and metabolic pathways in Saccharomyces cerevisiae (S. cerevisiae). However, native promoters in S. cerevisiae have several limitations which hinder their applications in metabolic engineering. These limitations include an inadequate number of well-characterized promoters, poor dynamic range, and insufficient orthogonality to endogenous regulations. Therefore, it is necessary to perform promoter engineering to create synthetic promoters with better properties. Here, we review recent advances related to promoter architecture, promoter engineering and synthetic promoter applications in S. cerevisiae. We also provide a perspective of future directions in this field with an emphasis on the recent advances of machine learning based promoter designs.