Determinants of nucleosome positioning

The importance of DNA sequence for nucleosome positioning in transcriptional regulation

Nucleosome positioning is a key factor for transcriptional regulation. Nucleosomes regulate the dynamic accessibility of chromatin and interact with the transcription machinery at every stage. Influences to steer nucleosome positioning are diverse, and the according importance of the DNA sequence in contrast to active chromatin remodeling has been the subject of long discussion. In this study, we evaluate the functional role of DNA sequence for all major elements along the process of transcription. We developed a random forest classifier based on local DNA structure that assesses the sequence-intrinsic support for nucleosome positioning. On this basis, we created a simple data resource that we applied genome-wide to the human genome. In our comprehensive analysis, we found a special role of DNA in mediating the competition of nucleosomes with cis-regulatory elements, in enabling steady transcription, for positioning of stable nucleosomes in exons, and for repelling nucleosomes during transcription termination. In contrast, we relate these findings to concurrent processes that generate strongly positioned nucleosomes in vivo that are not mediated by sequence, such as energy-dependent remodeling of chromatin.

Motorized chain models of the ideal chromosome

An array of motor proteins consumes chemical energy in setting up the architectures of chromosomes. Here, we explore how the structure of ideal polymer chains is influenced by two classes of motors. The first class which we call "swimming motors" acts to propel the chromatin fiber through three-dimensional space. They represent a caricature of motors such as RNA polymerases. Previously, they have often been described by adding a persistent flow onto Brownian diffusion of the chain. The second class of motors, which we call "grappling motors" caricatures the loop extrusion processes in which segments of chromatin fibers some distance apart are brought together. We analyze these models using a self-consistent variational phonon approximation to a many-body Master equation incorporating motor activities. We show that whether the swimming motors lead to contraction or expansion depends on the susceptibility of the motors, that is, how their activity depends on the forces they must exert. Grappling motors in contrast to swimming motors lead to long-ranged correlations that resemble those first suggested for fractal globules and that are consistent with the effective interactions inferred by energy landscape analyses of Hi-C data on the interphase chromosome.

nucMACC: An MNase-seq pipeline to identify structurally altered nucleosomes in the genome

Micrococcal nuclease sequencing is the state-of-the-art method for determining chromatin structure and nucleosome positioning. Data analysis is complex due to the AT-dependent sequence bias of the endonuclease and the requirement for high sequencing depth. Here, we present the nucleosome-based MNase accessibility (nucMACC) pipeline unveiling the regulatory chromatin landscape by measuring nucleosome accessibility and stability. The nucMACC pipeline represents a systematic and genome-wide approach for detecting unstable ("fragile") nucleosomes. We have characterized the regulatory nucleosome landscape in Drosophila melanogaster, Saccharomyces cerevisiae, and mammals. Two functionally distinct sets of promoters were characterized, one associated with an unstable nucleosome and the other being nucleosome depleted. We show that unstable nucleosomes present intermediate states of nucleosome remodeling, preparing inducible genes for transcriptional activation in response to stimuli or stress. The presence of unstable nucleosomes correlates with RNA polymerase II proximal pausing. The nucMACC pipeline offers unparalleled precision and depth in nucleosome research and is a valuable tool for future nucleosome studies.

Single-Macromolecule Studies of Eukaryotic Genomic Maintenance

Genomes are self-organized and self-maintained as long, complex macromolecules of chromatin. The inherent heterogeneity, stochasticity, phase separation, and chromatin dynamics of genome operation make it challenging to study genomes using ensemble methods. Various single-molecule force-, fluorescent-, and sequencing-based techniques rooted in different disciplines have been developed to fill critical gaps in the capabilities of bulk measurements, each providing unique, otherwise inaccessible, insights into the structure and maintenance of the genome. Capable of capturing molecular-level details about the organization, conformational changes, and packaging of genetic material, as well as processive and stochastic movements of maintenance factors, a single-molecule toolbox provides an excellent opportunity for collaborative research to understand how genetic material functions in health and malfunctions in disease. In this review, we discuss novel insights brought to genomic sciences by single-molecule techniques and their potential to continue to revolutionize the field-one molecule at a time.

RNA-directed DNA methylation mutants reduce histone methylation at the paramutated maize booster1 enhancer

Paramutation is the transfer of mitotically and meiotically heritable silencing information between two alleles. With paramutation at the maize (Zea mays) booster1 (b1) locus, the low-expressed B' epiallele heritably changes the high-expressed B-I epiallele into B' with 100% frequency. This requires specific tandem repeats and multiple components of the RNA-directed DNA methylation pathway, including the RNA-dependent RNA polymerase (encoded by mediator of paramutation1, mop1), the second-largest subunit of RNA polymerase IV and V (NRP(D/E)2a, encoded by mop2), and the largest subunit of RNA Polymerase IV (NRPD1, encoded by mop3). Mutations in mop genes prevent paramutation and release silencing at the B' epiallele. In this study, we investigated the effect of mutations in mop1, mop2, and mop3 on chromatin structure and DNA methylation at the B' epiallele, and especially the regulatory hepta-repeat 100 kb upstream of the b1 gene. Mutations in mop1 and mop3 resulted in decreased repressive histone modifications H3K9me2 and H3K27me2 at the hepta-repeat. Associated with this decrease were partial activation of the hepta-repeat enhancer function, formation of a multi-loop structure, and elevated b1 expression. In mop2 mutants, which do not show elevated b1 expression, H3K9me2, H3K27me2 and a single-loop structure like in wild-type B' were retained. Surprisingly, high CG and CHG methylation levels at the B' hepta-repeat remained in all three mutants, and CHH methylation was low in both wild type and mutants. Our results raise the possibility of MOP factors mediating RNA-directed histone methylation rather than RNA-directed DNA methylation at the b1 locus.

Acute depletion of BRG1 reveals its primary function as an activator of transcription

The mammalian SWI/SNF-like BAF complexes play critical roles during animal development and pathological conditions. Previous gene deletion studies and characterization of human gene mutations implicate that the complexes both repress and activate a large number of genes. However, the direct function of the complexes in cells remains largely unclear due to the relatively long-term nature of gene deletion or natural mutation. Here we generate a mouse line by knocking in the auxin-inducible degron tag (AID) to the Smarca4 gene, which encodes BRG1, the essential ATPase subunit of the BAF complexes. We show that the tagged BRG1 can be efficiently depleted by osTIR1 expression and auxin treatment for 6 to 10 h in CD4 + T cells, hepatocytes, and fibroblasts isolated from the knock-in mice. The acute depletion of BRG1 leads to decreases in nascent RNAs and RNA polymerase II binding at a large number of genes, which are positively correlated with the loss of BRG1. Further, these changes are correlated with diminished accessibility at DNase I Hypersensitive Sites (DHSs) and p300 binding. The acute BRG1 depletion results in three major patterns of nucleosome shifts leading to narrower nucleosome spacing surrounding transcription factor motifs and at enhancers and transcription start sites (TSSs), which are correlated with loss of BRG1, decreased chromatin accessibility and decreased nascent RNAs. Acute depletion of BRG1 severely compromises the Trichostatin A (TSA) -induced histone acetylation, suggesting a substantial interplay between the chromatin remodeling activity of BRG1 and histone acetylation. Our data suggest BRG1 mainly plays a direct positive role in chromatin accessibility, RNAPII binding, and nascent RNA production by regulating nucleosome positioning and facilitating transcription factor binding to their target sites.

Genome-wide nucleosome and transcription factor responses to genetic perturbations reveal chromatin-mediated mechanisms of transcriptional regulation

Epigenetic mechanisms contribute to gene regulation by altering chromatin accessibility through changes in transcription factor (TF) and nucleosome occupancy throughout the genome. Despite numerous studies focusing on changes in gene expression, the intricate chromatin-mediated regulatory code remains largely unexplored on a comprehensive scale. We address this by employing a factor-agnostic, reverse-genetics approach that uses MNase-seq to capture genome-wide TF and nucleosome occupancies in response to the individual deletion of 201 transcriptional regulators in Saccharomyces cerevisiae, thereby assaying nearly one million mutant-gene interactions. We develop a principled approach to identify and quantify chromatin changes genome-wide, observing differences in TF and nucleosome occupancy that recapitulate well-established pathways identified by gene expression data. We also discover distinct chromatin signatures associated with the up- and downregulation of genes, and use these signatures to reveal regulatory mechanisms previously unexplored in expression-based studies. Finally, we demonstrate that chromatin features are predictive of transcriptional activity and leverage these features to reconstruct chromatin-based transcriptional regulatory networks. Overall, these results illustrate the power of an approach combining genetic perturbation with high-resolution epigenomic profiling; the latter enables a close examination of the interplay between TFs and nucleosomes genome-wide, providing a deeper, more mechanistic understanding of the complex relationship between chromatin organization and transcription.

Circadian regulation of stereotypic chromatin conformations at enhancers

Cooperation between the circadian transcription factor (TF) CLOCK:BMAL1 and other TFs at cis-regulatory elements (CREs) is critical to daily rhythms of transcription. Yet, the modalities of this cooperation are unclear. Here, we analyzed the co-binding of multiple TFs on single DNA molecules in mouse liver using single molecule footprinting (SMF). We found that SMF reads clustered in stereotypic chromatin states that reflect distinguishable organization of TFs and nucleosomes, and that were remarkably conserved between all samples. DNA protection at CLOCK:BMAL1 binding motif (E-box) varied between CREs, from E-boxes being solely bound by CLOCK:BMAL1 to situations where other TFs competed with CLOCK:BMAL1 for E-box binding. SMF also uncovered CLOCK:BMAL1 cooperative binding at E-boxes separated by 250 bp, which structurally altered the CLOCK:BMAL1-DNA interface. Importantly, we discovered multiple nucleosomes with E-boxes at entry/exit sites that were removed upon CLOCK:BMAL1 DNA binding, thereby promoting the formation of open chromatin states that facilitate DNA binding of other TFs and that were associated with rhythmic transcription. These results demonstrate the utility of SMF for studying how CLOCK:BMAL1 and other TFs regulate stereotypical chromatin states at CREs to promote transcription.

Determinant of m6A regional preference by transcriptional dynamics

N6-Methyladenosine (m6A) is the most abundant chemical modification occurring on eukaryotic mRNAs, and has been reported to be involved in almost all stages of mRNA metabolism. The distribution of m6A sites is notably asymmetric along mRNAs, with a strong preference toward the 3' terminus of the transcript. How m6A regional preference is shaped remains incompletely understood. In this study, by performing m6A-seq on chromatin-associated RNAs, we found that m6A regional preference arises during transcription. Nucleosome occupancy is remarkedly increased in the region downstream of m6A sites, suggesting an intricate interplay between m6A methylation and nucleosome-mediated transcriptional dynamics. Notably, we found a remarkable slowdown of Pol-II movement around m6A sites. In addition, inhibiting Pol-II movement increases nearby m6A methylation levels. By analyzing massively parallel assays for m6A, we found that RNA secondary structures inhibit m6A methylation. Remarkably, the m6A sites associated with Pol-II pausing tend to be embedded within RNA secondary structures. These results suggest that Pol-II pausing could affect the accessibility of m6A motifs to the methyltransferase complex and subsequent m6A methylation by mediating RNA secondary structure. Overall, our study reveals a crucial role of transcriptional dynamics in the formation of m6A regional preference.

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.

SWI/SNF complex-mediated chromatin remodeling in Candida glabrata promotes immune evasion

Immune evasion is critical for fungal virulence. However, how the human opportunistic pathogen Candida glabrata (Cg) accomplishes this is unknown. Here, we present the first genome-wide nucleosome map of the macrophage-internalized Cg consisting of ∼12,000 dynamic and 70,000 total nucleosomes. We demonstrate that CgSnf2 (SWI/SNF chromatin remodeling complex-ATPase subunit)-mediated chromatin reorganization in macrophage-internalized Cg upregulates and downregulates the immunosuppressive seven-gene mannosyltransferase-cluster (CgMT-C) and immunostimulatory cell surface adhesin-encoding EPA1 gene, respectively. Consistently, EPA1 overexpression and CgMT-C deletion elevated IL-1β (pro-inflammatory cytokine) production and diminished Cg proliferation in macrophages. Further, Cgsnf2Δ had higher Epa1 surface expression, and evoked increased IL-1β secretion, and was killed in macrophages. Akt-, p38-, NF-κB- or NLRP3 inflammasome-inhibition partially reversed increased IL-1β secretion in Cgsnf2Δ-infected macrophages. Importantly, macrophages responded to multiple Candida pathogens via NF-κB-dependent IL-1β production, underscoring NF-κB signaling's role in fungal diseases. Altogether, our findings directly link the nucleosome positioning-based chromatin remodeling to fungal immunomodulatory molecule expression.

Sequence Dependence in Nucleosome Dynamics

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 μs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Cold-induced deposition of bivalent H3K4me3-H3K27me3 modification and nucleosome depletion in Arabidopsis

Epigenetic regulation of gene expression plays a crucial role in plant development and environmental adaptation. The H3K4me3 and H3K27me3 have not only been discovered in the regulation of gene expression in multiple biological processes but also in responses to abiotic stresses in plants. However, evidence for the presence of both H3K4me3 and H3K27me3 on the same nucleosome is sporadic. Cold-induced deposition of bivalent H3K4me3-H3K27me3 modifications and nucleosome depletion over a considerable number of active genes is documented in potato tubers and provides clues on an additional role of the bivalent modifications. Limited by the available information of genes encoding PcG/TrxG proteins as well as their corresponding mutants in potatoes, the molecular mechanism underlying the cold-induced deposition of the bivalent mark remains elusive. In this study, we found a similar deposition of the bivalent H3K4me3-H3K27me3 mark over 2129 active genes in cold-treated Arabidopsis Col-0 seedlings. The expression levels of the bivalent mark-associated genes tend to be independent of bivalent modification levels. However, these genes were associated with greater chromatin accessibility, presumably to provide a distinct chromatin environment for gene expression. In mutants clf28 and lhp1, failure to deposit H3K27me3 in active genes upon cold treatment implies that the CLF is potentially involved in cold-induced deposition of H3K27me3, with assistance from LHP1. Failure to deposit H3K4me3 during cold treatment in atx1-2 suggests a regulatory role of ATX1 in the deposition of H3K4me3. In addition, we observed a cold-induced global reduction in nucleosome occupancy, which is potentially mediated by LHP1 in an H3K27me3-dependent manner.

Chromatin region binning of gene expression for improving embryo cell subtype identification

Mammalian embryonic development is a complex process, characterized by intricate spatiotemporal dynamics and distinct chromatin preferences. However, the quick diversification in early embryogenesis leads to significant cellular diversity and the sparsity of scRNA-seq data, posing challenges in accurately determining cell fate decisions. In this study, we introduce a chromatin region binning method using scChrBin, designed to identify chromatin regions that elucidate the dynamics of embryonic development and lineage differentiation. This method transforms scRNA-seq data into a chromatin-based matrix, leveraging genomic annotations. Our results showed that the scChrBin method achieves high accuracy, with 98.0% and 89.2% on two single-cell embryonic datasets, demonstrating its effectiveness in analyzing complex developmental processes. We also systematically and comprehensively analysis of these key chromatin binning regions and their associated genes, focusing on their roles in lineage and stage development. The perspective of chromatin region binning method enables a comprehensive analysis of transcriptome data at the chromatin level, allowing us to unveil the dynamic expression of chromatin regions across temporal and spatial development. The tool is available as an application at https://github.com/liameihao/scChrBin.

Transcriptional repression across mitosis: mechanisms and functions

Transcription represents a central aspect of gene expression with RNA polymerase machineries (RNA Pol) driving the synthesis of RNA from DNA template molecules. In eukaryotes, a total of three RNA Pol enzymes generate the plethora of RNA species and RNA Pol II is the one transcribing all protein-coding genes. A high number of cis- and trans-acting factors orchestrates RNA Pol II-mediated transcription by influencing the chromatin recruitment, activation, elongation, and/or termination steps. The levels of DNA accessibility, defining open-euchromatin versus close-heterochromatin, delimits RNA Pol II activity as well as the encounter with other factors acting on chromatin such as the DNA replication or DNA repair machineries. The stage of the cell cycle highly influences RNA Pol II activity with mitosis representing the major challenge. In fact, there is a massive inhibition of transcription during the mitotic entry coupled with chromatin dissociation of most of the components of the transcriptional machinery. Mitosis, as a consequence, highly compromises the transcriptional memory and the perpetuation of cellular identity. Once mitosis ends, transcription levels immediately recover to define the cell fate and to safeguard the proper progression of daughter cells through the cell cycle. In this review, we evaluate our current understanding of the transcriptional repression associated with mitosis with a special focus on the molecular mechanisms involved, on the potential function behind the general repression, and on the transmission of the transcriptional machinery into the daughter cells. We finally discuss the contribution that errors in the inheritance of the transcriptional machinery across mitosis might play in stem cell aging.

Interplay between the transcription preinitiation complex and the +1 nucleosome

Eukaryotic transcription starts with the assembly of a preinitiation complex (PIC) on core promoters. Flanking this region is the +1 nucleosome, the first nucleosome downstream of the core promoter. While this nucleosome is rich in epigenetic marks and plays a key role in transcription regulation, how the +1 nucleosome interacts with the transcription machinery has been a long-standing question. Here, we summarize recent structural and functional studies of the +1 nucleosome in complex with the PIC. We specifically focus on how differently organized promoter-nucleosome templates affect the assembly of the PIC and PIC-Mediator on chromatin and result in distinct transcription initiation.

Collisions of RNA polymerases behind the replication fork promote alternative RNA splicing in newly replicated chromatin

DNA replication produces a global disorganization of chromatin structure that takes hours to be restored. However, how these chromatin rearrangements affect the regulation of gene expression and the maintenance of cell identity is not clear. Here, we use ChOR-seq and ChrRNA-seq experiments to analyze RNA polymerase II (RNAPII) activity and nascent RNA synthesis during the first hours after chromatin replication in human cells. We observe that transcription elongation is rapidly reactivated in nascent chromatin but that RNAPII abundance and distribution are altered, producing heterogeneous changes in RNA synthesis. Moreover, this first wave of transcription results in RNAPII blockages behind the replication fork, leading to changes in alternative splicing. Altogether, our results deepen our understanding of how transcriptional programs are regulated during cell division and uncover molecular mechanisms that explain why chromatin replication is an important source of gene expression variability.

A single fiber view of the nucleosome organization in eukaryotic chromatin

Eukaryotic cells are thought to arrange nucleosomes into extended arrays with evenly spaced nucleosomes phased at genomic landmarks. Here we tested to what extent this stereotypic organization describes the nucleosome organization in Saccharomyces cerevisiae using Fiber-Seq, a long-read sequencing technique that maps entire nucleosome arrays on individual chromatin fibers in a high throughput manner. With each fiber coming from a different cell, Fiber-Seq uncovers cell-to-cell heterogeneity. The long reads reveal the nucleosome architecture even over repetitive DNA such as the ribosomal DNA repeats. The absolute nucleosome occupancy, a parameter that is difficult to obtain with conventional sequencing approaches, is a direct readout of Fiber-Seq. We document substantial deviations from the stereotypical nucleosome organization with unexpectedly long linker DNAs between nucleosomes, gene bodies missing entire nucleosomes, cell-to-cell heterogeneity in nucleosome occupancy, heterogeneous phasing of arrays and irregular nucleosome spacing. Nucleosome array structures are indistinguishable throughout the gene body and with respect to the direction of transcription arguing against transcription promoting array formation. Acute nucleosome depletion destroyed most of the array organization indicating that nucleosome remodelers cannot efficiently pack nucleosomes under those conditions. Given that nucleosomes are cis-regulatory elements, the cell-to-cell heterogeneity uncovered by Fiber-Seq provides much needed information to understand chromatin structure and function.

Epigenetic remodelling under hypoxia

Hypoxia is intrinsic to tumours and contributes to malignancy and metastasis while hindering the efficiency of existing treatments. Epigenetic mechanisms play a crucial role in the regulation of hypoxic cancer cell programs, both in the initial phases of sensing the decrease in oxygen levels and during adaptation to chronic lack of oxygen. During the latter, the epigenetic regulation of tumour biology intersects with hypoxia-sensitive transcription factors in a complex network of gene regulation that also involves metabolic reprogramming. Here, we review the current literature on the epigenetic control of gene programs in hypoxic cancer cells. We highlight common themes and features of such epigenetic remodelling and discuss their relevance for the development of therapeutic strategies.

A genome-wide comprehensive analysis of nucleosome positioning in yeast

In eukaryotic cells, the one-dimensional DNA molecules need to be tightly packaged into the spatially constraining nucleus. Folding is achieved on its lowest level by wrapping the DNA around nucleosomes. Their arrangement regulates other nuclear processes, such as transcription and DNA repair. Despite strong efforts to study nucleosome positioning using Next Generation Sequencing (NGS) data, the mechanism of their collective arrangement along the gene body remains poorly understood. Here, we classify nucleosome distributions of protein-coding genes in Saccharomyces cerevisiae according to their profile similarity and analyse their differences using functional Principal Component Analysis. By decomposing the NGS signals into their main descriptive functions, we compared wild type and chromatin remodeler-deficient strains, keeping position-specific details preserved whilst considering the nucleosome arrangement as a whole. A correlation analysis with other genomic properties, such as gene size and length of the upstream Nucleosome Depleted Region (NDR), identified key factors that influence the nucleosome distribution. We reveal that the RSC chromatin remodeler-which is responsible for NDR maintenance-is indispensable for decoupling nucleosome arrangement within the gene from positioning outside, which interfere in rsc8-depleted conditions. Moreover, nucleosome profiles in chd1Δ strains displayed a clear correlation with RNA polymerase II presence, whereas wild type cells did not indicate a noticeable interdependence. We propose that RSC is pivotal for global nucleosome organisation, whilst Chd1 plays a key role for maintaining local arrangement.

Single molecule MATAC-seq reveals key determinants of DNA replication origin efficiency

Stochastic origin activation gives rise to significant cell-to-cell variability in the pattern of genome replication. The molecular basis for heterogeneity in efficiency and timing of individual origins is a long-standing question. Here, we developed Methylation Accessibility of TArgeted Chromatin domain Sequencing (MATAC-Seq) to determine single-molecule chromatin accessibility of four specific genomic loci. MATAC-Seq relies on preferential modification of accessible DNA by methyltransferases combined with Nanopore-Sequencing for direct readout of methylated DNA-bases. Applying MATAC-Seq to selected early-efficient and late-inefficient yeast replication origins revealed large heterogeneity of chromatin states. Disruption of INO80 or ISW2 chromatin remodeling complexes leads to changes at individual nucleosomal positions that correlate with changes in their replication efficiency. We found a chromatin state with an accessible nucleosome-free region in combination with well-positioned +1 and +2 nucleosomes as a strong predictor for efficient origin activation. Thus, MATAC-Seq identifies the large spectrum of alternative chromatin states that co-exist on a given locus previously masked in population-based experiments and provides a mechanistic basis for origin activation heterogeneity during eukaryotic DNA replication. Consequently, our single-molecule chromatin accessibility assay will be ideal to define single-molecule heterogeneity across many fundamental biological processes such as transcription, replication, or DNA repair in vitro and ex vivo.

The impact of the sequence-dependent physical properties of DNA on chromatin dynamics

The local mechanical properties of DNA depend on local sequence. Here we review recent genomic, structural, and computational efforts at deciphering the "mechanical code", i.e., the mapping between sequence and mechanics. We then discuss works that suggest how evolution has exploited the mechanical code to control the energetics of DNA-deforming biological processes such as nucleosome organization, transcription factor binding, DNA supercoiling, gene regulation, and 3D chromatin organization. As a whole, these recent works suggest that DNA sequence in diverse organisms can encode regulatory information governing diverse processes via the mechanical code.

Histone N-tails modulate sequence-specific positioning of nucleosomes

The precise mechanisms governing sequence-dependent positioning of nucleosomes on DNA remain unknown in detail. Existing algorithms, taking into account the sequence-dependent deformability of DNA and its interactions with the histone globular domains, predict rotational setting of only 65% of human nucleosomes mapped in vivo. To uncover novel factors responsible for the nucleosome positioning, we analyzed potential involvement of the histone N-tails in this process. To this aim, we reconstituted the H2A/H4 N-tailless nucleosomes on human BRCA1 DNA (~100 kb) and compared their positions and sequences with those of the wild-type nucleosomes. In the case of H2A tailless nucleosomes, the AT content of DNA sequences is changed locally at superhelical location (SHL) ±4, while maintaining the same rotational setting as their wild-type counterparts. Conversely, the H4 tailless nucleosomes display widespread changes of the AT content near SHL ±1 and SHL ±2, where the H4 N-tails interact with DNA. Furthermore, a substantial number of H4 tailless nucleosomes exhibit rotational setting opposite to that of the wild-type nucleosomes. Thus, our findings strongly suggest that the histone N-tails are operative in selection of nucleosome positions, which may have wide-ranging implications for epigenetic modulation of chromatin states.

Mathematical Models of HIV-1 Dynamics, Transcription, and Latency

HIV-1 latency is a major barrier to curing infections with antiretroviral therapy and, consequently, to eliminating the disease globally. The establishment, maintenance, and potential clearance of latent infection are complex dynamic processes and can be best described with the help of mathematical models followed by experimental validation. Here, we review the use of viral dynamics models for HIV-1, with a focus on applications to the latent reservoir. Such models have been used to explain the multi-phasic decay of viral load during antiretroviral therapy, the early seeding of the latent reservoir during acute infection and the limited inflow during treatment, the dynamics of viral blips, and the phenomenon of post-treatment control. Finally, we discuss that mathematical models have been used to predict the efficacy of potential HIV-1 cure strategies, such as latency-reversing agents, early treatment initiation, or gene therapies, and to provide guidance for designing trials of these novel interventions.

Poly(dA:dT) Tracts Differentially Modulate Nucleosome Remodeling Activity of RSC and ISW1a Complexes, Exerting Tract Orientation-Dependent and -Independent Effects

The establishment and maintenance of nucleosome-free regions (NFRs) are prominent processes within chromatin dynamics. Transcription factors, ATP-dependent chromatin remodeling complexes (CRCs) and DNA sequences are the main factors involved. In Saccharomyces cerevisiae, CRCs such as RSC contribute to chromatin opening at NFRs, while other complexes, including ISW1a, contribute to NFR shrinking. Regarding DNA sequences, growing evidence points to poly(dA:dT) tracts as playing a direct role in active processes involved in nucleosome positioning dynamics. Intriguingly, poly(dA:dT)-tract-containing NFRs span asymmetrically relative to the location of the tract by a currently unknown mechanism. In order to obtain insight into the role of poly(dA:dT) tracts in nucleosome remodeling, we performed a systematic analysis of their influence on the activity of ISW1a and RSC complexes. Our results show that poly(dA:dT) tracts differentially affect the activity of these CRCs. Moreover, we found differences between the effects exerted by the two alternative tract orientations. Remarkably, tract-containing linker DNA is taken as exit DNA for nucleosome sliding catalyzed by RSC. Our findings show that defined DNA sequences, when present in linker DNA, can dictate in which direction a remodeling complex has to slide nucleosomes and shed light into the mechanisms underlying asymmetrical chromatin opening around poly(dA:dT) tracts.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

Genome-wide single-molecule analysis of long-read DNA methylation reveals heterogeneous patterns at heterochromatin that reflect nucleosome organisation

High-throughput sequencing technology is central to our current understanding of the human methylome. The vast majority of studies use chemical conversion to analyse bulk-level patterns of DNA methylation across the genome from a population of cells. While this technology has been used to probe single-molecule methylation patterns, such analyses are limited to short reads of a few hundred basepairs. DNA methylation can also be directly detected using Nanopore sequencing which can generate reads measuring megabases in length. However, thus far these analyses have largely focused on bulk-level assessment of DNA methylation. Here, we analyse DNA methylation in single Nanopore reads from human lymphoblastoid cells, to show that bulk-level metrics underestimate large-scale heterogeneity in the methylome. We use the correlation in methylation state between neighbouring sites to quantify single-molecule heterogeneity and find that heterogeneity varies significantly across the human genome, with some regions having heterogeneous methylation patterns at the single-molecule level and others possessing more homogeneous methylation patterns. By comparing the genomic distribution of the correlation to epigenomic annotations, we find that the greatest heterogeneity in single-molecule patterns is observed within heterochromatic partially methylated domains (PMDs). In contrast, reads originating from euchromatic regions and gene bodies have more ordered DNA methylation patterns. By analysing the patterns of single molecules in more detail, we show the existence of a nucleosome-scale periodicity in DNA methylation that accounts for some of the heterogeneity we uncover in long single-molecule DNA methylation patterns. We find that this periodic structure is partially masked in bulk data and correlates with DNA accessibility as measured by nanoNOMe-seq, suggesting that it could be generated by nucleosomes. Our findings demonstrate the power of single-molecule analysis of long-read data to understand the structure of the human methylome.

Mechanical codes of chemical-scale specificity in DNA motifs

In gene transcription, certain sequences of double-stranded (ds)DNA play a vital role in nucleosome positioning and expression initiation. That dsDNA is deformed to various extents in these processes leads us to ask: Could the genomic DNA also have sequence specificity in its chemical-scale mechanical properties? We approach this question using statistical machine learning to determine the rigidity between DNA chemical moieties. What emerges for the polyA, polyG, TpA, and CpG sequences studied here is a unique trigram that contains the quantitative mechanical strengths between bases and along the backbone. In a way, such a sequence-dependent trigram could be viewed as a DNA mechanical code. Interestingly, we discover a compensatory competition between the axial base-stacking interaction and the transverse base-pairing interaction, and such a reciprocal relationship constitutes the most discriminating feature of the mechanical code. Our results also provide chemical-scale understanding for experimental observables. For example, the long polyA persistence length is shown to have strong base stacking while its complement (polyAc) exhibits high backbone rigidity. The mechanical code concept enables a direct reading of the physical interactions encoded in the sequence which, with further development, is expected to shed new light on DNA allostery and DNA-binding drugs.

Direct observation of coordinated assembly of individual native centromeric nucleosomes

Eukaryotic chromosome segregation requires the kinetochore, a megadalton-sized machine that forms on specialized centromeric chromatin containing CENP-A, a histone H3 variant. CENP-A deposition requires a chaperone protein HJURP that targets it to the centromere, but it has remained unclear whether HJURP has additional functions beyond CENP-A targeting and why high AT DNA content, which disfavors nucleosome assembly, is widely conserved at centromeres. To overcome the difficulties of studying nucleosome formation in vivo, we developed a microscopy assay that enables direct observation of de novo centromeric nucleosome recruitment and maintenance with single molecule resolution. Using this assay, we discover that CENP-A can arrive at centromeres without its dedicated centromere-specific chaperone HJURP, but stable incorporation depends on HJURP and additional DNA-binding proteins of the inner kinetochore. We also show that homopolymer AT runs in the yeast centromeres are essential for efficient CENP-A deposition. Together, our findings reveal requirements for stable nucleosome formation and provide a foundation for further studies of the assembly and dynamics of native kinetochore complexes.

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.

Recent development of selective inhibitors targeting the HDAC6 as anti-cancer drugs: Structure, function and design

HDAC6, a member of the histone deacetylase family, mainly is a cytosolic protein and regulates cell growth by acting on non-histone substrates, such as α -tubulin, cortactin, heat shock protein HSP90, programmed death 1 (PD-1) and programmed death ligand 1 (PD-L1), that are closely related to the proliferation, invasion, immune escape and angiogenesis of cancer tissues. The approved drugs targeting the HDACs are all pan-inhibitors and have many side effects due to their lack of selectivity. Therefore, development of selective inhibitors of HDAC6 has attracted much attention in the field of cancer therapy. In this review, we will summarize the relationship between HDAC6 and cancer, and discuss the design strategies of HDAC6 inhibitors for cancer treatment in recent years.

A half century of exploring DNA excision repair in chromatin

DNA in eukaryotic cells is packaged into the compact and dynamic structure of chromatin. This packaging is a double-edged sword for DNA repair and genomic stability. Chromatin restricts the access of repair proteins to DNA lesions embedded in nucleosomes and higher order chromatin structures. However, chromatin also serves as a signaling platform in which post-translational modifications of histones and other chromatin-bound proteins promote lesion recognition and repair. Similarly, chromatin modulates the formation of DNA damage, promoting or suppressing lesion formation depending on the chromatin context. Therefore, the modulation of DNA damage and its repair in chromatin is crucial to our understanding of the fate of potentially mutagenic and carcinogenic lesions in DNA. Here, we survey many of the landmark findings on DNA damage and repair in chromatin over the last 50 years (i.e., since the beginning of this field), focusing on excision repair, the first repair mechanism studied in the chromatin landscape. For example, we highlight how the impact of chromatin on these processes explains the distinct patterns of somatic mutations observed in cancer genomes.

Role of Nucleosome Sliding in the Protein Target Search for Covered DNA Sites

Associations of transcription factors (TFs) with specific sites on DNA initiate major cellular processes. But DNA in eukaryotic cells is covered by nucleosomes which prevent TFs from binding. However, nucleosome structures on DNA are not static and exhibit breathing and sliding. We develop a theoretical framework to investigate the effect of nucleosome sliding on a protein target search. By analysis of a discrete-state stochastic model of nucleosome sliding, search dynamics are explicitly evaluated. It is found that for long sliding lengths the target search dynamics are faster for normal TFs that cannot enter the nucleosomal DNA. But for more realistic short sliding lengths, the so-called pioneer TFs, which can invade nucleosomal DNA, locate specific sites faster. It is also suggested that nucleosome breathing, which is a faster process, has a stronger effect on protein search dynamics than that of nucleosome sliding. Theoretical arguments to explain these observations are presented.

Transcription and FACT facilitate the restoration of replication-coupled chromatin assembly defects

Genome duplication occurs through the coordinated action of DNA replication and nucleosome assembly at replication forks. Defective nucleosome assembly causes DNA lesions by fork breakage that need to be repaired. In addition, it causes a loss of chromatin integrity. These chromatin alterations can be restored, even though the mechanisms are unknown. Here, we show that the process of chromatin restoration can deal with highly severe chromatin defects induced by the absence of the chaperones CAF1 and Rtt106 or a strong reduction in the pool of available histones, and that this process can be followed by analyzing the topoisomer distribution of the 2µ plasmid. Using this assay, we demonstrate that chromatin restoration is slow and independent of checkpoint activation, whereas it requires the action of transcription and the FACT complex. Therefore, cells are able to "repair" not only DNA lesions but also chromatin alterations associated with defective nucleosome assembly.

Investigating pioneer factor activity and its coordination with chromatin remodelers using integrated synthetic oligo assay

Chromatin accessibility is regulated by pioneer factors (PFs) and chromatin remodelers (CRs). Here, we present a protocol, based on integrated synthetic oligonucleotide libraries in yeast, to systematically interrogate the nucleosome-displacing activities of PFs and their coordination with CRs. We describe steps for designing oligonucleotide sequences, constructing yeast libraries, measuring nucleosome configurations, and data analyses. This approach potentially can be adapted for use in higher eukaryotes to investigate the activities of many types of chromatin-associated factors. For complete details on the use and execution of this protocol, please refer to Yan et al.,1 and Chen et al.2.

A complex network of transcription factors and epigenetic regulators involved in bovine leukemia virus transcriptional regulation

Bovine Leukemia Virus (BLV) is the etiological agent of enzootic bovine leukosis, a disease characterized by the neoplastic proliferation of B cells in cattle. While most European countries have introduced efficient eradication programs, BLV is still present worldwide and no treatment is available. A major feature of BLV infection is the viral latency, which enables the escape from the host immune system, the maintenance of a persistent infection and ultimately the tumoral development. BLV latency is a multifactorial phenomenon resulting in the silencing of viral genes due to genetic and epigenetic repressions of the viral promoter located in the 5' Long Terminal Repeat (5'LTR). However, viral miRNAs and antisense transcripts are expressed from two different proviral regions, respectively the miRNA cluster and the 3'LTR. These latter transcripts are expressed despite the viral latency affecting the 5'LTR and are increasingly considered to take part in tumoral development. In the present review, we provide a summary of the experimental evidence that has enabled to characterize the molecular mechanisms regulating each of the three BLV transcriptional units, either through cis-regulatory elements or through epigenetic modifications. Additionally, we describe the recently identified BLV miRNAs and antisense transcripts and their implications in BLV-induced tumorigenesis. Finally, we discuss the relevance of BLV as an experimental model for the closely related human T-lymphotropic virus HTLV-1.

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.

HSV-1 Infection Induces a Downstream Shift of Promoter-Proximal Pausing for Host Genes

Herpes simplex virus 1 (HSV-1) infection exerts a profound shutoff of host gene expression at multiple levels. Recently, HSV-1 infection was reported to also impact promoter-proximal RNA polymerase II (Pol II) pausing, a key step in the eukaryotic transcription cycle, with decreased and increased Pol II pausing observed for activated and repressed genes, respectively. Here, we demonstrate that HSV-1 infection induces more complex alterations in promoter-proximal pausing than previously suspected for the vast majority of cellular genes. While pausing is generally retained, it is shifted to more downstream and less well-positioned sites for most host genes. The downstream shift of Pol II pausing was established between 1.5 and 3 h of infection, remained stable until at least 6 hours postinfection, and was observed in the absence of ICP22. The shift in Pol II pausing does not result from alternative de novo transcription initiation at downstream sites or read-in transcription originating from disruption of transcription termination of upstream genes. The use of downstream secondary pause sites associated with +1 nucleosomes was previously observed upon negative elongation factor (NELF) depletion. However, downstream shifts of Pol II pausing in HSV-1 infection were much more pronounced than observed upon NELF depletion. Thus, our study reveals a novel aspect in which HSV-1 infection fundamentally reshapes host transcriptional processes, providing new insights into the regulation of promoter-proximal Pol II pausing in eukaryotic cells. IMPORTANCE This study provides a genome-wide analysis of changes in promoter-proximal polymerase II (Pol II) pausing on host genes induced by HSV-1 infection. It shows that standard measures of pausing, i.e., pausing indices, do not properly capture the complex and unsuspected alterations in Pol II pausing occurring in HSV-1 infection. Instead of a reduction of pausing with increased elongation, as suggested by pausing index analysis, HSV-1 infection leads to a shift of pausing to downstream and less well-positioned sites than in uninfected cells for the majority of host genes. Thus, HSV-1 infection fundamentally reshapes a key regulatory step at the beginning of the host transcriptional cycle on a genome-wide scale.

SeATAC: a tool for exploring the chromatin landscape and the role of pioneer factors

Assay for Transposase-Accessible Chromatin with sequencing (ATAC-seq) reveals chromatin accessibility across the genome. Currently, no method specifically detects differential chromatin accessibility. Here, SeATAC uses a conditional variational autoencoder model to learn the latent representation of ATAC-seq V-plots and outperforms MACS2 and NucleoATAC on six separate tasks. Applying SeATAC to several pioneer factor-induced differentiation or reprogramming ATAC-seq datasets suggests that induction of these factors not only relaxes the closed chromatin but also decreases chromatin accessibility of 20% to 30% of their target sites. SeATAC is a novel tool to accurately reveal genomic regions with differential chromatin accessibility from ATAC-seq data.

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.

Direct observation of coordinated assembly of individual native centromeric nucleosomes

Eukaryotic chromosome segregation requires the kinetochore, a megadalton-sized machine that forms on specialized centromeric chromatin containing CENP-A, a histone H3 variant. CENP-A deposition requires a chaperone protein HJURP that targets it to the centromere, but it has remained unclear whether HJURP has additional functions beyond CENP-A targeting and why high AT DNA content, which disfavors nucleosome assembly, is widely conserved at centromeres. To overcome the difficulties of studying nucleosome formation in vivo, we developed a microscopy assay that enables direct observation of de novo centromeric nucleosome recruitment and maintenance with single molecule resolution. Using this assay, we discover that CENP-A can arrive at centromeres without its dedicated centromere-specific chaperone HJURP, but stable incorporation depends on HJURP and additional DNA-binding proteins of the inner kinetochore. We also show that homopolymer AT runs in the yeast centromeres are essential for efficient CENP-A deposition. Together, our findings reveal requirements for stable nucleosome formation and provide a foundation for further studies of the assembly and dynamics of native kinetochore complexes.

Structural basis of transcription reduction by a promoter-proximal +1 nucleosome

At active human genes, the +1 nucleosome is located downstream of the RNA polymerase II (RNA Pol II) pre-initiation complex (PIC). However, at inactive genes, the +1 nucleosome is found further upstream, at a promoter-proximal location. Here, we establish a model system to show that a promoter-proximal +1 nucleosome can reduce RNA synthesis in vivo and in vitro, and we analyze its structural basis. We find that the PIC assembles normally when the edge of the +1 nucleosome is located 18 base pairs (bp) downstream of the transcription start site (TSS). However, when the nucleosome edge is located further upstream, only 10 bp downstream of the TSS, the PIC adopts an inhibited state. The transcription factor IIH (TFIIH) shows a closed conformation and its subunit XPB contacts DNA with only one of its two ATPase lobes, inconsistent with DNA opening. These results provide a mechanism for nucleosome-dependent regulation of transcription initiation.

Systematic evaluation of chromatin immunoprecipitation sequencing to study histone occupancy in dormancy transitions of grapevine buds

The regulation of DNA accessibility by histone modification has emerged as a paradigm of developmental and environmental programming. Chromatin immunoprecipitation followed by sequencing (ChIP-seq) is a versatile tool to investigate in vivo protein-DNA interaction and has enabled advances in mechanistic understanding of physiologies. The technique has been successfully demonstrated in several plant species and tissues; however, it has remained challenging in woody tissues, in particular complex structures such as perennating buds. Here we developed a ChIP method specifically for mature dormant buds of grapevine (Vitis vinifera cv. Cabernet Sauvignon). Each step of the protocol was systematically optimized, including crosslinking, chromatin extraction, sonication and antibody validation. Analysis of histone H3-enriched DNA was performed to evaluate the success of the protocol and identify occupancy of histone H3 along grapevine bud chromatin. To our best knowledge, this is the first ChIP experiment protocol optimized for the grapevine bud system.

Ageing-associated changes in transcriptional elongation influence longevity

Physiological homeostasis becomes compromised during ageing, as a result of impairment of cellular processes, including transcription and RNA splicing1-4. However, the molecular mechanisms leading to the loss of transcriptional fidelity are so far elusive, as are ways of preventing it. Here we profiled and analysed genome-wide, ageing-related changes in transcriptional processes across different organisms: nematodes, fruitflies, mice, rats and humans. The average transcriptional elongation speed (RNA polymerase II speed) increased with age in all five species. Along with these changes in elongation speed, we observed changes in splicing, including a reduction of unspliced transcripts and the formation of more circular RNAs. Two lifespan-extending interventions, dietary restriction and lowered insulin-IGF signalling, both reversed most of these ageing-related changes. Genetic variants in RNA polymerase II that reduced its speed in worms5 and flies6 increased their lifespan. Similarly, reducing the speed of RNA polymerase II by overexpressing histone components, to counter age-associated changes in nucleosome positioning, also extended lifespan in flies and the division potential of human cells. Our findings uncover fundamental molecular mechanisms underlying animal ageing and lifespan-extending interventions, and point to possible preventive measures.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

Global transcriptome and gene co-expression network analyses reveal regulatory and non-additive effects of drought and heat stress in grapevine

Despite frequent co-occurrence of drought and heat stress, the molecular mechanisms governing plant responses to these stresses in combination have not often been studied. This is particularly evident in non-model, perennial plants. We conducted large scale physiological and transcriptome analyses to identify genes and pathways associated with grapevine response to drought and/or heat stress during stress progression and recovery. We identified gene clusters with expression correlated to leaf temperature and water stress and five hub genes for the combined stress co-expression network. Several differentially expressed genes were common to the individual and combined stresses, but the majority were unique to the individual or combined stress treatments. These included heat-shock proteins, mitogen-activated kinases, sugar metabolizing enzymes, and transcription factors, while phenylpropanoid biosynthesis and histone modifying genes were unique to the combined stress treatment. Following physiological recovery, differentially expressed genes were found only in plants under heat stress, both alone and combined with drought. Taken collectively, our results suggest that the effect of the combined stress on physiology and gene expression is more severe than that of individual stresses, but not simply additive, and that epigenetic chromatin modifications may play an important role in grapevine responses to combined drought and heat stress.

DeepBend: An interpretable model of DNA bendability

The bendability of genomic DNA impacts chromatin packaging and protein-DNA binding. However, we do not have a comprehensive understanding of the motifs influencing DNA bendability. Recent high-throughput technologies such as Loop-Seq offer an opportunity to address this gap but the lack of accurate and interpretable machine learning models still remains. Here we introduce DeepBend, a convolutional neural network model with convolutions designed to directly capture the motifs underlying DNA bendability and their periodic occurrences or relative arrangements that modulate bendability. DeepBend consistently performs on par with alternative models while giving an extra edge through mechanistic interpretations. Besides confirming the known motifs of DNA bendability, DeepBend also revealed several novel motifs and showed how the spatial patterns of motif occurrences influence bendability. DeepBend's genome-wide prediction of bendability further showed how bendability is linked to chromatin conformation and revealed the motifs controlling the bendability of topologically associated domains and their boundaries.

Loss of the Ash2l subunit of histone H3K4 methyltransferase complexes reduces chromatin accessibility at promoters

Changes in gene expression programs are intimately linked to cell fate decisions. Post-translational modifications of core histones contribute to control gene expression. Methylation of lysine 4 of histone H3 (H3K4) correlates with active promoters and gene transcription. This modification is catalyzed by KMT2 methyltransferases, which require interaction with 4 core subunits, WDR5, RBBP5, ASH2L and DPY30, for catalytic activity. Ash2l is necessary for organismal development and for tissue homeostasis. In mouse embryo fibroblasts (MEFs), Ash2l loss results in gene repression, provoking a senescence phenotype. We now find that upon knockout of Ash2l both H3K4 mono- and tri-methylation (H3K4me1 and me3, respectively) were deregulated. In particular, loss of H3K4me3 at promoters correlated with gene repression, especially at CpG island promoters. Ash2l loss resulted in increased loading of histone H3 and reduced chromatin accessibility at promoters, accompanied by an increase of repressing and a decrease of activating histone marks. Moreover, we observed altered binding of CTCF upon Ash2l loss. Lost and gained binding was noticed at promoter-associated and intergenic sites, respectively. Thus, Ash2l loss and reduction of H3K4me3 correlate with altered chromatin accessibility and transcription factor binding. These findings contribute to a more detailed understanding of mechanistic consequences of H3K4me3 loss and associated repression of gene transcription and thus of the observed cellular consequences.