Nucleosomes play a dual role in regulating transcription dynamics.

Transcription has a mechanical component, as the translocation of the transcription machinery or RNA polymerase (RNAP) on DNA or chromatin is dynamically coupled to the chromatin torsion. This posits chromatin mechanics as a possible regulator of eukaryotic transcription, however, the modes and mechanisms of this regulation are elusive. Here, we first take a statistical mechanics approach to model the torsional response of topology-constrained chromatin. Our model recapitulates the experimentally observed weaker torsional stiffness of chromatin compared to bare DNA and proposes structural transitions of nucleosomes into chirally distinct states as the driver of the contrasting torsional mechanics. Coupling chromatin mechanics with RNAP translocation in stochastic simulations, we reveal a complex interplay of DNA supercoiling and nucleosome dynamics in governing RNAP velocity. Nucleosomes play a dual role in controlling the transcription dynamics. The steric barrier aspect of nucleosomes in the gene body counteracts transcription via hindering RNAP motion, whereas the chiral transitions facilitate RNAP motion via driving a low restoring torque upon twisting the DNA. While nucleosomes with low dissociation rates are typically transcriptionally repressive, highly dynamic nucleosomes offer less of a steric barrier and enhance the transcription elongation dynamics of weakly transcribed genes via buffering DNA twist. We use the model to predict transcription-dependent levels of DNA supercoiling in segments of the budding yeast genome that are in accord with available experimental data. The model unveils a paradigm of DNA supercoiling-mediated interaction between genes and makes testable predictions that will guide experimental design.

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.

Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins.

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

In silico design of DNA sequences for in vivo nucleosome positioning.

The computational design of synthetic DNA sequences with designer in vivo properties is gaining traction in the field of synthetic genomics. We propose here a computational method which combines a kinetic Monte Carlo framework with a deep mutational screening based on deep learning predictions. We apply our method to build regular nucleosome arrays with tailored nucleosomal repeat lengths (NRL) in yeast. Our design was validated in vivo by successfully engineering and integrating thousands of kilobases long tandem arrays of computationally optimized sequences which could accommodate NRLs much larger than the yeast natural NRL (namely 197 and 237 bp, compared to the natural NRL of ∼165 bp). RNA-seq results show that transcription of the arrays can occur but is not driven by the NRL. The computational method proposed here delineates the key sequence rules for nucleosome positioning in yeast and should be easily applicable to other sequence properties and other genomes.

Genome wide nucleosome landscape shapes 3D chromatin organization.

The hierarchical chromatin organization begins with formation of nucleosomes, which fold into chromatin domains punctuated by boundaries and ultimately chromosomes. In a hierarchal organization, lower levels shape higher levels. However, the dependence of higher-order 3D chromatin organization on the nucleosome-level organization has not been studied in cells. We investigated the relationship between nucleosome-level organization and higher-order chromatin organization by perturbing nucleosomes across the genome by deleting Imitation SWItch (ISWI) and Chromodomain Helicase DNA-binding (CHD1) chromatin remodeling factors in budding yeast. We find that changes in nucleosome-level properties are accompanied by changes in 3D chromatin organization. Short-range chromatin contacts up to a few kilo-base pairs decrease, chromatin domains weaken, and boundary strength decreases. Boundary strength scales with accessibility and moderately with width of nucleosome-depleted region. Change in nucleosome positioning seems to alter the stiffness of chromatin, which can affect formation of chromatin contacts. Our results suggest a biomechanical "bottom-up" mechanism by which nucleosome distribution across genome shapes 3D chromatin organization.

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.

MNase-Seq Analysis for Identifying Stress-Altered Nucleosome Occupancy in Plants.

Nucleosome occupancy plays an important role in chromatin compaction, affecting biological processes by hampering the binding of cis-acting elements such as transcription factors, RNA polymerase machinery, and coregulatory. Accessible regions allow for cis-acting elements to bind DNA and regulate transcription. Here, we detail our protocol to profile nucleosome occupancy and chromatin structure dynamics under drought stress at the genome-wide scale using micrococcal nuclease (MNase) digestion. Combining variable MNase concentration treatments and high-throughput sequencing, we investigate the changes in the overall chromatin state using bread wheat samples from an exemplary drought experiment.

The PCNA-Pol δ complex couples lagging strand DNA synthesis to parental histone transfer for epigenetic inheritance.

Inheritance of epigenetic information is critical for maintaining cell identity. The transfer of parental histone H3-H4 tetramers, the primary carrier of epigenetic modifications on histone proteins, represents a crucial yet poorly understood step in the inheritance of epigenetic information. Here, we show the lagging strand DNA polymerase, Pol δ, interacts directly with H3-H4 and that the interaction between Pol δ and the sliding clamp PCNA regulates parental histone transfer to lagging strands, most likely independent of their roles in DNA synthesis. When combined, mutations at Pol δ and Mcm2 that compromise parental histone transfer result in a greater reduction in nucleosome occupancy at nascent chromatin than mutations in either alone. Last, PCNA contributes to nucleosome positioning on nascent chromatin. On the basis of these results, we suggest that the PCNA-Pol δ complex couples lagging strand DNA synthesis to parental H3-H4 transfer, facilitating epigenetic inheritance.

C2c: Predicting Micro-C from Hi-C.

MOTIVATION: High-resolution Hi-C data, capable of detecting chromatin features below the level of Topologically Associating Domains (TADs), significantly enhance our understanding of gene regulation. Micro-C, a variant of Hi-C incorporating a micrococcal nuclease (MNase) digestion step to examine interactions between nucleosome pairs, has been developed to overcome the resolution limitations of Hi-C. However, Micro-C experiments pose greater technical challenges compared to Hi-C, owing to the need for precise MNase digestion control and higher-resolution sequencing. Therefore, developing computational methods to derive Micro-C data from existing Hi-C datasets could lead to better usage of a large amount of existing Hi-C data in the scientific community and cost savings. RESULTS: We developed C2c ("high" or upper case C to "micro" or lower case c), a computational tool based on a residual neural network to learn the mapping between Hi-C and Micro-C contact matrices and then predict Micro-C contact matrices based on Hi-C contact matrices. Our evaluation results show that the predicted Micro-C contact matrices reveal more chromatin loops than the input Hi-C contact matrices, and more of the loops detected from predicted Micro-C match the promoter-enhancer interactions. Furthermore, we found that the mutual loops from real and predicted Micro-C better match the ChIA-PET data compared to Hi-C and real Micro-C loops, and the predicted Micro-C leads to more TAD-boundaries detected compared to the Hi-C data. The website URL of C2c can be found in the Data Availability Statement.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF.

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Self-assembling viral histones are evolutionary intermediates between archaeal and eukaryotic nucleosomes.

Nucleosomes are DNA-protein complexes composed of histone proteins that form the basis of eukaryotic chromatin. The nucleosome was a key innovation during eukaryotic evolution, but its origin from histone homologues in Archaea remains unclear. Viral histone repeats, consisting of multiple histone paralogues within a single protein, may reflect an intermediate state. Here we examine the diversity of histones encoded by Nucleocytoviricota viruses. We identified 258 histones from 168 viral metagenomes with variable domain configurations including histone singlets, doublets, triplets and quadruplets, the latter comprising the four core histones arranged in series. Viral histone repeats branch phylogenetically between Archaea and eukaryotes and display intermediate functions in Escherichia coli, self-assembling into eukaryotic-like nucleosomes that stack into archaeal-like oligomers capable of impacting genomic activity and condensing DNA. Histone linkage also facilitates nucleosome formation, promoting eukaryotic histone assembly in E. coli. These data support the hypothesis that viral histone repeats originated in stem-eukaryotes and that nucleosome evolution proceeded through histone repeat intermediates.

Triple lysine and nucleosome-binding motifs of the viral IE19 protein are required for human cytomegalovirus S-phase infections.

Herpesvirus genomes are maintained as extrachromosomal plasmids within the nuclei of infected cells. Some herpesviruses persist within dividing cells, putting the viral genome at risk of being lost to the cytoplasm during mitosis because karyokinesis (nuclear division) requires nuclear envelope breakdown. Oncogenic herpesviruses (and papillomaviruses) avoid genome loss during mitosis by tethering their genomes to cellular chromosomes, thereby ensuring viral genome uptake into newly formed nuclei. These viruses use viral proteins with DNA- and chromatin-binding capabilities to physically link viral and cellular genomes together in a process called tethering. The known viral tethering proteins of human papillomavirus (E2), Epstein-Barr virus (EBNA1), and Kaposi's sarcoma-associated herpesvirus (LANA) each contain two independent domains required for genome tethering, one that binds sequence specifically to the viral genome and another that binds to cellular chromatin. This latter domain is called a chromatin tethering domain (CTD). The human cytomegalovirus UL123 gene encodes a CTD that is required for the virus to productively infect dividing fibroblast cells within the S phase of the cell cycle, presumably by tethering the viral genome to cellular chromosomes during mitosis. The CTD-containing UL123 gene product that supports S-phase infections is the IE19 protein. Here, we define two motifs in IE19 required for S-phase infections: an N-terminal triple lysine motif and a C-terminal nucleosome-binding motif within the CTD.IMPORTANCEThe IE19 protein encoded by human cytomegalovirus (HCMV) is required for S-phase infection of dividing cells, likely because it tethers the viral genome to cellular chromosomes, thereby allowing them to survive mitosis. The mechanism through which IE19 tethers viral genomes to cellular chromosomes is not understood. For human papillomavirus, Epstein-Barr virus, and Kaposi's sarcoma-associated herpesvirus, viral genome tethering is required for persistence (latency) and pathogenesis (oncogenesis). Like these viruses, HCMV also achieves latency, and it modulates the properties of glioblastoma multiforme tumors. Therefore, defining the mechanism through which IE19 tethers viral genomes to cellular chromosomes may help us understand, and ultimately combat or control, HCMV latency and oncomodulation.

FACT maintains chromatin architecture and thereby stimulates RNA polymerase II pausing during transcription in vivo.

Facilitates chromatin transcription (FACT) is a histone chaperone that supports transcription through chromatin in vitro, but its functional roles in vivo remain unclear. Here, we analyze the in vivo functions of FACT with the use of multi-omics analysis after rapid FACT depletion from human cells. We show that FACT depletion destabilizes chromatin and leads to transcriptional defects, including defective promoter-proximal pausing and elongation, and increased premature termination of RNA polymerase II. Unexpectedly, our analysis revealed that promoter-proximal pausing depends not only on the negative elongation factor (NELF) but also on the +1 nucleosome, which is maintained by FACT.

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.

General regulatory factors exert differential effects on nucleosome sliding activity of the ISW1a complex.

BACKGROUND: Chromatin dynamics is deeply involved in processes that require access to DNA, such as transcriptional regulation. Among the factors involved in chromatin dynamics at gene regulatory regions are general regulatory factors (GRFs). These factors contribute to establishment and maintenance of nucleosome-depleted regions (NDRs). These regions are populated by nucleosomes through histone deposition and nucleosome sliding, the latter catalyzed by a number of ATP-dependent chromatin remodeling complexes, including ISW1a. It has been observed that GRFs can act as barriers against nucleosome sliding towards NDRs. However, the relative ability of the different GRFs to hinder sliding activity is currently unknown. RESULTS: Considering this, we performed a comparative analysis for the main GRFs, with focus in their ability to modulate nucleosome sliding mediated by ISW1a. Among the GRFs tested in nucleosome remodeling assays, Rap1 was the only factor displaying the ability to hinder the activity of ISW1a. This effect requires location of the Rap1 cognate sequence on linker that becomes entry DNA in the nucleosome remodeling process. In addition, Rap1 was able to hinder nucleosome assembly in octamer transfer assays. Concurrently, Rap1 displayed the highest affinity for and longest dwell time from its target sequence, compared to the other GRFs tested. Consistently, through bioinformatics analyses of publicly available genome-wide data, we found that nucleosome occupancy and histone deposition in vivo are inversely correlated with the affinity of Rap1 for its target sequences in the genome. CONCLUSIONS: Our findings point to DNA binding affinity, residence time and location at particular translational positions relative to the nucleosome core as the key features of GRFs underlying their roles played in nucleosome sliding and assembly.

Nucleosomal DNA has topological memory.

One elusive aspect of the chromosome architecture is how it constrains the DNA topology. Nucleosomes stabilise negative DNA supercoils by restraining a DNA linking number difference (∆Lk) of about -1.26. However, whether this capacity is uniform across the genome is unknown. Here, we calculate the ∆Lk restrained by over 4000 nucleosomes in yeast cells. To achieve this, we insert each nucleosome in a circular minichromosome and perform Topo-seq, a high-throughput procedure to inspect the topology of circular DNA libraries in one gel electrophoresis. We show that nucleosomes inherently restrain distinct ∆Lk values depending on their genomic origin. Nucleosome DNA topologies differ at gene bodies (∆Lk = -1.29), intergenic regions (∆Lk = -1.23), rDNA genes (∆Lk = -1.24) and telomeric regions (∆Lk = -1.07). Nucleosomes near the transcription start and termination sites also exhibit singular DNA topologies. Our findings demonstrate that nucleosome DNA topology is imprinted by its native chromatin context and persists when the nucleosome is relocated.

DNA methylation-based high-resolution mapping of long-distance chromosomal interactions in nucleosome-depleted regions.

3C-based methods have significantly advanced our understanding of 3D genome organization. However, it remains a formidable task to precisely capture long-range chromosomal interactions between individual loci, such as those between promoters and distal enhancers. Here, we present Methyltransferase Targeting-based chromosome Architecture Capture (MTAC), a method that maps the contacts between a target site (viewpoint) and the rest of the genome in budding yeast with high resolution and sensitivity. MTAC detects hundreds of intra- and inter-chromosomal interactions within nucleosome-depleted regions (NDRs) that cannot be captured by 4C, Hi-C, or Micro-C. By applying MTAC to various viewpoints, we find that (1) most long-distance chromosomal interactions detected by MTAC reflect tethering by the nuclear pore complexes (NPCs), (2) genes co-regulated by methionine assemble into inter-chromosomal clusters near NPCs upon activation, (3) mediated by condensin, the mating locus forms a highly specific interaction with the recombination enhancer (RE) in a mating-type specific manner, and (4) correlation of MTAC signals among NDRs reveal spatial mixing and segregation of the genome. Overall, these results demonstrate MTAC as a powerful tool to resolve fine-scale long-distance chromosomal interactions and provide insights into the 3D genome organization.

Chromatin regulates alternative polyadenylation via the RNA polymerase II elongation rate.

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

Evolutionary Dynamics of Chromatin Structure and Duplicate Gene Expression in Diploid and Allopolyploid Cotton.

Polyploidy is a prominent mechanism of plant speciation and adaptation, yet the mechanistic understandings of duplicated gene regulation remain elusive. Chromatin structure dynamics are suggested to govern gene regulatory control. Here, we characterized genome-wide nucleosome organization and chromatin accessibility in allotetraploid cotton, Gossypium hirsutum (AADD, 2n = 4X = 52), relative to its two diploid parents (AA or DD genome) and their synthetic diploid hybrid (AD), using DNS-seq. The larger A-genome exhibited wider average nucleosome spacing in diploids, and this intergenomic difference diminished in the allopolyploid but not hybrid. Allopolyploidization also exhibited increased accessibility at promoters genome-wide and synchronized cis-regulatory motifs between subgenomes. A prominent cis-acting control was inferred for chromatin dynamics and demonstrated by transposable element removal from promoters. Linking accessibility to gene expression patterns, we found distinct regulatory effects for hybridization and later allopolyploid stages, including nuanced establishment of homoeolog expression bias and expression level dominance. Histone gene expression and nucleosome organization are coordinated through chromatin accessibility. Our study demonstrates the capability to track high-resolution chromatin structure dynamics and reveals their role in the evolution of cis-regulatory landscapes and duplicate gene expression in polyploids, illuminating regulatory ties to subgenomic asymmetry and dominance.

asteRIa enables robust interaction modeling between chromatin modifications and epigenetic readers.

Chromatin, the nucleoprotein complex consisting of DNA and histone proteins, plays a crucial role in regulating gene expression by controlling access to DNA. Chromatin modifications are key players in this regulation, as they help to orchestrate DNA transcription, replication, and repair. These modifications recruit epigenetic 'reader' proteins, which mediate downstream events. Most modifications occur in distinctive combinations within a nucleosome, suggesting that epigenetic information can be encoded in combinatorial chromatin modifications. A detailed understanding of how multiple modifications cooperate in recruiting such proteins has, however, remained largely elusive. Here, we integrate nucleosome affinity purification data with high-throughput quantitative proteomics and hierarchical interaction modeling to estimate combinatorial effects of chromatin modifications on protein recruitment. This is facilitated by the computational workflow asteRIa which combines hierarchical interaction modeling, stability-based model selection, and replicate-consistency checks for a stable estimation of Robust Interactions among chromatin modifications. asteRIa identifies several epigenetic reader candidates responding to specific interactions between chromatin modifications. For the polycomb protein CBX8, we independently validate our results using genome-wide ChIP-Seq and bisulphite sequencing datasets. We provide the first quantitative framework for identifying cooperative effects of chromatin modifications on protein binding.