Genome-scale identification of nucleosome positions in S. cerevisiae

H2A.Z deposition by the SWR complex is stimulated by polyadenine DNA sequences in nucleosomes

The variant histone H2A.Z is deposited into nucleosomes immediately downstream of promoters, where it plays a critical role in transcription. The site-specific deposition of H2A.Z is catalyzed by the SWR complex, a conserved chromatin remodeler with affinity for promoter-proximal nucleosome-depleted regions (NDRs) and histone acetylation. By comparing the genomic distribution of H2A.Z in wild-type and SWR-deficient cells, we found that SWR is also responsible for depositing H2A.Z at thousands of non-canonical sites not directly linked to NDRs or histone acetylation. To understand the targeting mechanism of H2A.Z, we presented SWR to a library of canonical nucleosomes isolated from yeast and analyzed the preferred substrates. Our results revealed that SWR preferentially deposited H2A.Z into a subset of endogenous H2A.Z sites, which are overrepresented by polyadenine tracts on the top strands of the DNA duplex at the nucleosomal entry-exit sites. Insertion of polyadenine sequences into recombinant nucleosomes near the outgoing H2A-H2B dimer enhanced SWR's affinity for the nucleosomal substrate and increased its H2A.Z insertion activity. These findings suggest that the genome encodes sequence-based information that facilitates remodeler-mediated targeting of H2A.Z.

Identifying associations between short tandem repeat sequences and gene expression in yeast reveals specific repeated motifs encoding transcriptional regulatory proteins

Tandem repeat sequences (TRs), a class of repetitive genomic elements, are broadly distributed in both coding and non-coding regions. Investigating the relationship between sequences and function is essential for understanding the genome. Saccharomyces cerevisiae serves as a vital model organism and is widely used as an engineered strain. Although the transcriptional regulatory functions of TRs in the promoters of S.cerevisiae have been elucidated, our understanding of their roles within coding sequences (CDS) remains limited. In this study, we integrate RNA-seq, ChIP-seq, ATAC-seq, Hi-C, and Micro-C data from S.cerevisiae to analyze the types and distribution of TRs, and their impact on gene expression. Our results indicate that genes containing short tandem repeats (STRs) in their CDS exhibit lower expression levels. Epigenetic analysis reveals that these regions are characterized by high levels of repressive histone modifications and low levels of activating marks, with reduced chromatin accessibility and fewer chromatin interactions. Furthermore, trinucleotide and hexanucleotide repeated motifs of STR are found primarily enriched in genes encoding transcriptional regulatory proteins. This study provides new insights into the functions and characteristics of STRs in the CDS of S.cerevisiae. The identification of key STR motifs offers potential targets for the design of transcriptional regulatory elements.

Thermodynamics of nucleosome breathing and positioning

Nucleosomes are fundamental units of chromatin in which a length of genomic DNA is wrapped around a histone octamer spool in a left-handed superhelix. Large-scale nucleosome maps show a wide distribution of DNA wrapping lengths, which in some cases are tens of base pairs (bp) shorter than the 147 bp canonical wrapping length observed in nucleosome crystal structures. Here, we develop a thermodynamic model that assumes a constant free energy cost of unwrapping a nucleosomal bp. Our model also incorporates linker DNA-short DNA segments between neighboring nucleosomes imposed by the folding of nucleosome arrays into chromatin fibers and other higher-order chromatin structures. We use this model to study nucleosome positioning and occupancy in the presence of nucleosome "breathing"-partial unwrapping and rewrapping of nucleosomal DNA due to interactions with the neighboring particles. We find that, as the unwrapping cost per bp and the chemical potential are varied, the nucleosome arrays are characterized by three distinct states, with low, intermediate, and high densities. The transition between the latter two states proceeds through an equiprobable state in which all nucleosome wrapping lengths are equally likely. We study the equiprobable state theoretically using a mean-field approach, obtaining an excellent agreement with numerical simulations. Finally, we use our model to reproduce S. cerevisiae nucleosome occupancy profiles observed in the vicinity of transcription start sites, as well as genome-wide distributions of nucleosome wrapping lengths. Overall, our results highlight the key role of partial nucleosome unwrapping in shaping the genome-wide patterns of nucleosome positioning and occupancy.

Footprint-C reveals transcription factor modes in local clusters and long-range chromatin interactions

The proximity ligation-based Hi-C and derivative methods are the mainstream tools to study genome-wide chromatin interactions. These methods often fragment the genome using enzymes functionally irrelevant to the interactions per se, restraining the efficiency in identifying structural features and the underlying regulatory elements. Here we present Footprint-C, which yields high-resolution chromatin contact maps built upon intact and genuine footprints protected by transcription factor (TF) binding. When analyzed at one-dimensional level, the billions of chromatin contacts from Footprint-C enable genome-wide analysis at single footprint resolution, and reveal preferential modes of local TF co-occupancy. At pairwise contact level, Footprint-C exhibits higher efficiency in identifying chromatin structural features when compared with other Hi-C methods, segregates chromatin interactions emanating from adjacent TF footprints, and uncovers multiway interactions involving different TFs. Altogether, Footprint-C results suggest that rich regulatory modes of TF may underlie both local residence and distal chromatin interactions, in terms of TF identity, valency, and conformational configuration.

RNA Structure: Past, Future, and Gene Therapy Applications

First believed to be a simple intermediary between the information encoded in deoxyribonucleic acid and that functionally displayed in proteins, ribonucleic acid (RNA) is now known to have many functions through its abundance and intricate, ubiquitous, diverse, and dynamic structure. About 70-90% of the human genome is transcribed into protein-coding and noncoding RNAs as main determinants along with regulatory sequences of cellular to populational biological diversity. From the nucleotide sequence or primary structure, through Watson-Crick pairing self-folding or secondary structure, to compaction via longer distance Watson-Crick and non-Watson-Crick interactions or tertiary structure, and interactions with RNA or other biopolymers or quaternary structure, or with metabolites and biomolecules or quinary structure, RNA structure plays a critical role in RNA's lifecycle from transcription to decay and many cellular processes. In contrast to the success of 3-dimensional protein structure prediction using AlphaFold, RNA tertiary and beyond structures prediction remains challenging. However, approaches involving machine learning and artificial intelligence, sequencing of RNA and its modifications, and structural analyses at the single-cell and intact tissue levels, among others, provide an optimistic outlook for the continued development and refinement of RNA-based applications. Here, we highlight those in gene therapy.

In vivo assembly of complete eukaryotic nucleosomes and (H3-H4)-only non-canonical nucleosomal particles in the model bacterium Escherichia coli

As a fundamental unit for packaging genomic DNA into chromatin, the eukaryotic nucleosome core comprises a canonical octamer with two copies for each histone, H2A, H2B, H3, and H4, wrapped around with 147 base pairs of DNA. While H3 and H4 share structure-fold with archaeal histone-like proteins, the eukaryotic nucleosome core and the complete nucleosome (the core plus H1 histone) are unique to eukaryotes. To explore whether the eukaryotic nucleosome can assemble in prokaryotes and to reconstruct the possible route for its emergence in eukaryogenesis, we developed an in vivo system for assembly of nucleosomes in the model bacterium, Escherichia coli, and successfully reconstituted the core nucleosome, the complete nucleosome, and unexpectedly the non-canonical (H3-H4)4 octasome. The core and complete nucleosomes assembled in E. coli exhibited footprints typical of eukaryotic hosts after in situ micrococcal nuclease digestion. Additionally, they caused condensation of E. coli nucleoid. We also demonstrated the stable formation of non-canonical (H3-H4)2 tetrasome and (H3-H4)4 octasomes in vivo, which are suggested to be 'fossil complex' that marks the intermediate in the progressive development of eukaryotic nucleosome. The study presents a unique platform in a bacterium for in vivo assembly and studying the properties of non-canonical variants of nucleosome.

Seasonal and circadian rhythms of clerodane diterpenes and glycosylated flavonoids in two varieties of Casearia sylvestris Sw. (Salicaceae)

Casearia sylvestris Sw. (Salicaceae) is noted for its morphological and chemical plasticity and pharmacological properties. The present study investigates two of its varieties: C. sylvestris var. sylvestris, predominant in dense and humid forests and ecotones and characterized by clerodane diterpenes; and C. sylvestris var. lingua, mainly found in xeric and open savannah areas and containing phenolic compounds. Despite their comprehensive chemical profiles, the dynamics of clerodane diterpenes and glycosylated flavonoids remain unknown. This study thus aimed to describe seasonal and circadian variations in their content in the leaves of the two varieties. The relative contents of five diterpenes and three glycosylated flavonoids were monitored monthly, every 3 h for 48 h, over 1 year via high-performance liquid chromatography coupled to diode array detection (HPLC-UV-DAD). The differential expression of photosynthetic proteins (Rubisco and photosystem II) was analyzed by Western blotting. The contents of both chemical classes decreased during the reproductive stage, though the prevalence of diterpenes in var. sylvestris and flavonoids in var. lingua remained unchanged; furthermore, even when the plants are grown under the same geographic and environmental conditions, Rubisco expression in var. lingua is twice that of var. sylvestris. In var. lingua, photosystem II proteins are 10 % less expressed. The study reveals the circadian and seasonal fluctuations and, thus, prevalence of the two main compound classes in the examined varieties. The expression of the investigated photosynthetic proteins provides insights into the two varieties, supporting the prevalence of var. lingua in Cerrado areas and var. sylvestris in Atlantic Forest areas.

An integrated machine-learning model to predict nucleosome architecture

We demonstrate that nucleosomes placed in the gene body can be accurately located from signal decay theory assuming two emitters located at the beginning and at the end of genes. These generated wave signals can be in phase (leading to well defined nucleosome arrays) or in antiphase (leading to fuzzy nucleosome architectures). We found that the first (+1) and the last (-last) nucleosomes are contiguous to regions signaled by transcription factor binding sites and unusual DNA physical properties that hinder nucleosome wrapping. Based on these analyses, we developed a method that combines Machine Learning and signal transmission theory able to predict the basal locations of the nucleosomes with an accuracy similar to that of experimental MNase-seq based methods.

Screening of broad-host expression promoters for shuttle expression vectors in non-conventional yeasts and bacteria

BACKGROUND

Non-conventional yeasts and bacteria gain significance in synthetic biology for their unique metabolic capabilities in converting low-cost renewable feedstocks into valuable products. Improving metabolic pathways and increasing bioproduct yields remain dependent on the strategically use of various promoters in these microbes. The development of broad-spectrum promoter libraries with varying strengths for different hosts is attractive for biosynthetic engineers.

RESULTS

In this study, five Yarrowia lipolytica constitutive promoters (yl.hp4d, yl.FBA1in, yl.TEF1, yl.TDH1, yl.EXP1) and five Kluyveromyces marxianus constitutive promoters (km.PDC1, km.FBA1, km.TEF1, km.TDH3, km.ENO1) were selected to construct promoter-reporter vectors, utilizing α-amylase and red fluorescent protein (RFP) as reporter genes. The promoters' strengths were systematically characterized across Y. lipolytica, K. marxianus, Pichia pastoris, Escherichia coli, and Corynebacterium glutamicum. We discovered that five K. marxianus promoters can all express genes in Y. lipolytica and that five Y. lipolytica promoters can all express genes in K. marxianus with variable expression strengths. Significantly, the yl.TEF1 and km.TEF1 yeast promoters exhibited their adaptability in P. pastoris, E. coli, and C. glutamicum. In yeast P. pastoris, the yl.TEF1 promoter exhibited substantial expression of both amylase and RFP. In bacteria E. coli and C. glutamicum, the eukaryotic km.TEF1 promoter demonstrated robust expression of RFP. Significantly, in E. coli, The RFP expression strength of the km.TEF1 promoter reached ∼20% of the T7 promoter.

CONCLUSION

Non-conventional yeast promoters with diverse and cross-domain applicability have great potential for developing innovative and dynamic regulated systems that can effectively manage carbon flux and enhance target bioproduct synthesis across diverse microbial hosts.

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.

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here, we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

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.

Energy-driven genome regulation by ATP-dependent chromatin remodellers

The packaging of DNA into chromatin in eukaryotes regulates gene transcription, DNA replication and DNA repair. ATP-dependent chromatin remodelling enzymes (re)arrange nucleosomes at the first level of chromatin organization. Their Snf2-type motor ATPases alter histone-DNA interactions through a common DNA translocation mechanism. Whether remodeller activities mainly catalyse nucleosome dynamics or accurately co-determine nucleosome organization remained unclear. In this Review, we discuss the emerging mechanisms of chromatin remodelling: dynamic remodeller architectures and their interactions, the inner workings of the ATPase cycle, allosteric regulation and pathological dysregulation. Recent mechanistic insights argue for a decisive role of remodellers in the energy-driven self-organization of chromatin, which enables both stability and plasticity of genome regulation - for example, during development and stress. Different remodellers, such as members of the SWI/SNF, ISWI, CHD and INO80 families, process (epi)genetic information through specific mechanisms into distinct functional outputs. Combinatorial assembly of remodellers and their interplay with histone modifications, histone variants, DNA sequence or DNA-bound transcription factors regulate nucleosome mobilization or eviction or histone exchange. Such input-output relationships determine specific nucleosome positions and compositions with distinct DNA accessibilities and mediate differential genome regulation. Finally, remodeller genes are often mutated in diseases characterized by genome dysregulation, notably in cancer, and we discuss their physiological relevance.

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.

Incorporation of the histone variant H2A.Z counteracts gene silencing mediated by H3K27 trimethylation in Fusarium fujikuroi

BACKGROUND

Fusarium fujikuroi is a pathogen of rice causing diverse disease symptoms such as 'bakanae' or stunting, most likely due to the production of various natural products (NPs) during infection. Fusaria have the genetic potential to synthesize a plethora of these compounds with often diverse bioactivity. The capability to synthesize NPs exceeds the number of those being produced by far, implying a gene regulatory network decisive to induce production. One such regulatory layer is the chromatin structure and chromatin-based modifications associated with it. One prominent example is the exchange of histones against histone variants such as the H2A variant H2A.Z. Though H2A.Z already is well studied in several model organisms, its regulatory functions are not well understood. Here, we used F. fujikuroi as a model to explore the role of the prominent histone variant FfH2A.Z in gene expression within euchromatin and facultative heterochromatin.

RESULTS

Through the combination of diverse '-omics' methods, we show the global distribution of FfH2A.Z and analyze putative crosstalks between the histone variant and two prominent histone marks, i.e., H3K4me3 and H3K27me3, important for active gene transcription and silencing, respectively. We demonstrate that, if FfH2A.Z is positioned at the + 1-nucleosome, it poises chromatin for gene transcription, also within facultative heterochromatin. Lastly, functional characterization of FfH2A.Z overexpression and depletion mutants revealed that FfH2A.Z is important for wild type-like fungal development and secondary metabolism.

CONCLUSION

In this study, we show that the histone variant FfH2A.Z is a mark of positive gene transcription and acts independently of the chromatin state most likely through the stabilization of the + 1-nucleosome. Furthermore, we demonstrate that FfH2A.Z depletion does not influence the establishment of both H3K27me3 and H3K4me3, thus indicating no crosstalk between FfH2A.Z and both histone marks. These results highlight the manifold functions of the histone variant FfH2A.Z in the phytopathogen F. fujikuroi, which are distinct regarding gene transcription and crosstalk with the two prominent histone marks H3K27me3 and H3K4me3, as proposed for other model organisms.

Dynamic 1D search and processive nucleosome translocations by RSC and ISW2 chromatin remodelers

Eukaryotic gene expression is linked to chromatin structure and nucleosome positioning by ATP-dependent chromatin remodelers that establish and maintain nucleosome-depleted regions (NDRs) near transcription start sites. Conserved yeast RSC and ISW2 remodelers exert antagonistic effects on nucleosomes flanking NDRs, but the temporal dynamics of remodeler search, engagement, and directional nucleosome mobilization for promoter accessibility are unknown. Using optical tweezers and two-color single-particle imaging, we investigated the Brownian diffusion of RSC and ISW2 on free DNA and sparse nucleosome arrays. RSC and ISW2 rapidly scan DNA by one-dimensional hopping and sliding, respectively, with dynamic collisions between remodelers followed by recoil or apparent co-diffusion. Static nucleosomes block remodeler diffusion resulting in remodeler recoil or sequestration. Remarkably, both RSC and ISW2 use ATP hydrolysis to translocate mono-nucleosomes processively at ~30 bp/s on extended linear DNA under tension. Processivity and opposing push-pull directionalities of nucleosome translocation shown by RSC and ISW2 shape the distinctive landscape of promoter chromatin.

Structural insights into histone exchange by human SRCAP complex

Histone variant H2A.Z is found at promoters and regulates transcription. The ATP-dependent chromatin remodeler SRCAP complex (SRCAP-C) promotes the replacement of canonical histone H2A-H2B dimer with H2A.Z-H2B dimer. Here, we determined structures of human SRCAP-C bound to H2A-containing nucleosome at near-atomic resolution. The SRCAP subunit integrates a 6-subunit actin-related protein (ARP) module and an ATPase-containing motor module. The ATPase-associated ARP module encircles half of the nucleosome along the DNA and may restrain net DNA translocation, a unique feature of SRCAP-C. The motor module adopts distinct nucleosome binding modes in the apo (nucleotide-free), ADP-bound, and ADP-BeFx-bound states, suggesting that ATPase-driven movement destabilizes H2A-H2B by unwrapping the entry DNA and pulls H2A-H2B out of nucleosome through the ZNHIT1 subunit. Structure-guided chromatin immunoprecipitation sequencing analysis confirmed the requirement of H2A-contacting ZNHIT1 in maintaining H2A.Z occupancy on the genome. Our study provides structural insights into the mechanism of H2A-H2A.Z exchange mediated by SRCAP-C.

Dynamic 1D Search and Processive Nucleosome Translocations by RSC and ISW2 Chromatin Remodelers

Eukaryotic gene expression is linked to chromatin structure and nucleosome positioning by ATP-dependent chromatin remodelers that establish and maintain nucleosome-depleted regions (NDRs) near transcription start-sites. Conserved yeast RSC and ISW2 remodelers exert antagonistic effects on nucleosomes flanking NDRs, but the temporal dynamics of remodeler search, engagement and directional nucleosome mobilization for promoter accessibility are unknown. Using optical tweezers and 2-color single-particle imaging, we investigated the Brownian diffusion of RSC and ISW2 on free DNA and sparse nucleosome arrays. RSC and ISW2 rapidly scan DNA by one-dimensional hopping and sliding respectively, with dynamic collisions between remodelers followed by recoil or apparent co-diffusion. Static nucleosomes block remodeler diffusion resulting in remodeler recoil or sequestration. Remarkably, both RSC and ISW2 use ATP hydrolysis to translocate mono-nucleosomes processively at ~30 bp/sec on extended linear DNA under tension. Processivity and opposing push-pull directionalities of nucleosome translocation shown by RSC and ISW2 shape the distinctive landscape of promoter chromatin.

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.

Genome-wide mapping and cryo-EM structural analyses of the overlapping tri-nucleosome composed of hexasome-hexasome-octasome moieties

The nucleosome is a fundamental unit of chromatin in which about 150 base pairs of DNA are wrapped around a histone octamer. The overlapping di-nucleosome has been proposed as a product of chromatin remodeling around the transcription start site, and previously found as a chromatin unit, in which about 250 base pairs of DNA continuously bind to the histone core composed of a hexamer and an octamer. In the present study, our genome-wide analysis of human cells suggests another higher nucleosome stacking structure, the overlapping tri-nucleosome, which wraps about 300-350 base-pairs of DNA in the region downstream of certain transcription start sites of actively transcribed genes. We determine the cryo-electron microscopy (cryo-EM) structure of the overlapping tri-nucleosome, in which three subnucleosome moieties, hexasome, hexasome, and octasome, are associated by short connecting DNA segments. Small angle X-ray scattering and coarse-grained molecular dynamics simulation analyses reveal that the cryo-EM structure of the overlapping tri-nucleosome may reflect its structure in solution. Our findings suggest that nucleosome stacking structures composed of hexasome and octasome moieties may be formed by nucleosome remodeling factors around transcription start sites for gene regulation.

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.

MINI-AC: inference of plant gene regulatory networks using bulk or single-cell accessible chromatin profiles

Gene regulatory networks (GRNs) represent the interactions between transcription factors (TF) and their target genes. Plant GRNs control transcriptional programs involved in growth, development, and stress responses, ultimately affecting diverse agricultural traits. While recent developments in accessible chromatin (AC) profiling technologies make it possible to identify context-specific regulatory DNA, learning the underlying GRNs remains a major challenge. We developed MINI-AC (Motif-Informed Network Inference based on Accessible Chromatin), a method that combines AC data from bulk or single-cell experiments with TF binding site (TFBS) information to learn GRNs in plants. We benchmarked MINI-AC using bulk AC datasets from different Arabidopsis thaliana tissues and showed that it outperforms other methods to identify correct TFBS. In maize, a crop with a complex genome and abundant distal AC regions, MINI-AC successfully inferred leaf GRNs with experimentally confirmed, both proximal and distal, TF-target gene interactions. Furthermore, we showed that both AC regions and footprints are valid alternatives to infer AC-based GRNs with MINI-AC. Finally, we combined MINI-AC predictions from bulk and single-cell AC datasets to identify general and cell-type specific maize leaf regulators. Focusing on C4 metabolism, we identified diverse regulatory interactions in specialized cell types for this photosynthetic pathway. MINI-AC represents a powerful tool for inferring accurate AC-derived GRNs in plants and identifying known and novel candidate regulators, improving our understanding of gene regulation in plants.

Nucleic acid sequence contributes to remodeler-mediated targeting of histone H2A.Z

The variant histone H2A.Z is inserted into nucleosomes immediately downstream of promoters and is important for transcription. The site-specific deposition of H2A.Z is catalyzed by SWR, a conserved chromatin remodeler with affinity for promoter-proximal nucleosome depleted regions (NDRs) and histone acetylation. By comparing the genomic distribution of H2A.Z in wild-type and SWR-deficient cells, we found that SWR is also responsible for depositing H2A.Z at thousands of non-canonical sites not directly linked to NDRs or histone acetylation. To understand the targeting mechanism of H2A.Z, we presented SWR with a library of nucleosomes isolated from yeast and characterized those preferred by SWR. We found that SWR prefers nucleosomes associated with intergenic over coding regions, especially when polyadenine tracks are present. Insertion of polyadenine sequences into recombinant nucleosomes near the H2A-H2B binding site stimulated the H2A.Z insertion activity of SWR. Therefore, the genome is encoded with information contributing to remodeler-mediated targeting of H2A.Z.

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

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

Examining chromatin heterogeneity through PacBio long-read sequencing of M.EcoGII methylated genomes: an m6A detection efficiency and calling bias correcting pipeline

Recent studies have combined DNA methyltransferase footprinting of genomic DNA in nuclei with long-read sequencing, resulting in detailed chromatin maps for multi-kilobase stretches of genomic DNA from one cell. Theoretically, nucleosome footprints and nucleosome-depleted regions can be identified using M.EcoGII, which methylates adenines in any sequence context, providing a high-resolution map of accessible regions in each DNA molecule. Here we report PacBio long-read sequence data for budding yeast nuclei treated with M.EcoGII and a bioinformatic pipeline which corrects for three key challenges undermining this promising method. First, detection of m6A in individual DNA molecules by the PacBio software is inefficient, resulting in false footprints predicted by random gaps of seemingly unmethylated adenines. Second, there is a strong bias against m6A base calling as AT content increases. Third, occasional methylation occurs within nucleosomes, breaking up their footprints. After correcting for these issues, our pipeline calculates a correlation coefficient-based score indicating the extent of chromatin heterogeneity within the cell population for every gene. Although the population average is consistent with that derived using other techniques, we observe a wide range of heterogeneity in nucleosome positions at the single-molecule level, probably reflecting cellular chromatin dynamics.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Improving genome-wide mapping of nucleosomes in Trypanosome cruzi

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

Chromatin analysis of adult pluripotent stem cells reveals a unique stemness maintenance strategy

Many highly regenerative organisms maintain adult pluripotent stem cells throughout their life, but how the long-term maintenance of pluripotency is accomplished is unclear. To decipher the regulatory logic of adult pluripotent stem cells, we analyzed the chromatin organization of stem cell genes in the planarian Schmidtea mediterranea. We identify a special chromatin state of stem cell genes, which is distinct from that of tissue-specific genes and resembles constitutive genes. Where tissue-specific promoters have detectable transcription factor binding sites, the promoters of stem cell-specific genes instead have sequence features that broadly decrease nucleosome binding affinity. This genic organization makes pluripotency-related gene expression the default state in these cells, which is maintained by the activity of chromatin remodelers ISWI and SNF2 in the stem cells.

Measuring occupancies of the nucleosome and nucleosome-interacting factors in vivo in Saccharomyces cerevisiae genome-wide

Nucleosomes are the repeating units of chromatin. The presence of nucleosomes poses a major impediment to all DNA-dependent processes. As a result, access to DNA in chromatin is dynamically regulated by many factors, including ATP-dependent chromatin remodeling complexes. Digestion of chromatin by micrococcal nuclease (MNase) followed by chromatin immunoprecipitation (ChIP) and sequencing can be leveraged to determine nucleosome occupancy, positioning, and the ability of chromatin interacting factors to alter chromatin accessibility. Here we describe the procedure for performing MNase and MNase ChIP-seq in detail.

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.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Prediction of G4 formation in live cells with epigenetic data: a deep learning approach

G-quadruplexes (G4s) are secondary structures abundant in DNA that may play regulatory roles in cells. Despite the ubiquity of the putative G-quadruplex-forming sequences (PQS) in the human genome, only a small fraction forms G4 structures in cells. Folded G4, histone methylation and chromatin accessibility are all parts of the complex cis regulatory landscape. We propose an approach for prediction of G4 formation in cells that incorporates epigenetic and chromatin accessibility data. The novel approach termed epiG4NN efficiently predicts cell-specific G4 formation in live cells based on a local epigenomic snapshot. Our results confirm the close relationship between H3K4me3 histone methylation, chromatin accessibility and G4 structure formation. Trained on A549 cell data, epiG4NN was then able to predict G4 formation in HEK293T and K562 cell lines. We observe the dependency of model performance with different epigenetic features on the underlying experimental condition of G4 detection. We expect that this approach will contribute to the systematic understanding of correlations between structural and epigenomic feature landscape.

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.

Nucleosome Remodeling at the Yeast PHO8 and PHO84 Promoters without the Putatively Essential SWI/SNF Remodeler

Chromatin remodeling by ATP-dependent remodeling enzymes is crucial for all genomic processes, like transcription or replication. Eukaryotes harbor many remodeler types, and it is unclear why a given chromatin transition requires more or less stringently one or several remodelers. As a classical example, removal of budding yeast PHO8 and PHO84 promoter nucleosomes upon physiological gene induction by phosphate starvation essentially requires the SWI/SNF remodeling complex. This dependency on SWI/SNF may indicate specificity in remodeler recruitment, in recognition of nucleosomes as remodeling substrate or in remodeling outcome. By in vivo chromatin analyses of wild type and mutant yeast under various PHO regulon induction conditions, we found that overexpression of the remodeler-recruiting transactivator Pho4 allowed removal of PHO8 promoter nucleosomes without SWI/SNF. For PHO84 promoter nucleosome removal in the absence of SWI/SNF, an intranucleosomal Pho4 site, which likely altered the remodeling outcome via factor binding competition, was required in addition to such overexpression. Therefore, an essential remodeler requirement under physiological conditions need not reflect substrate specificity, but may reflect specific recruitment and/or remodeling outcomes.

Recycling of modified H2A-H2B provides short-term memory of chromatin states

Chromatin landscapes are disrupted during DNA replication and must be restored faithfully to maintain genome regulation and cell identity. The histone H3-H4 modification landscape is restored by parental histone recycling and modification of new histones. How DNA replication impacts on histone H2A-H2B is currently unknown. Here, we measure H2A-H2B modifications and H2A.Z during DNA replication and across the cell cycle using quantitative genomics. We show that H2AK119ub1, H2BK120ub1, and H2A.Z are recycled accurately during DNA replication. Modified H2A-H2B are segregated symmetrically to daughter strands via POLA1 on the lagging strand, but independent of H3-H4 recycling. Post-replication, H2A-H2B modification and variant landscapes are quickly restored, and H2AK119ub1 guides accurate restoration of H3K27me3. This work reveals epigenetic transmission of parental H2A-H2B during DNA replication and identifies cross talk between H3-H4 and H2A-H2B modifications in epigenome propagation. We propose that rapid short-term memory of recycled H2A-H2B modifications facilitates restoration of stable H3-H4 chromatin states.

Structures of transcription preinitiation complex engaged with the +1 nucleosome

The preinitiation complex (PIC) assembles on promoters of protein-coding genes to position RNA polymerase II (Pol II) for transcription initiation. Previous structural studies revealed the PIC on different promoters, but did not address how the PIC assembles within chromatin. In the yeast Saccharomyces cerevisiae, PIC assembly occurs adjacent to the +1 nucleosome that is located downstream of the core promoter. Here we present cryo-EM structures of the yeast PIC bound to promoter DNA and the +1 nucleosome located at three different positions. The general transcription factor TFIIH engages with the incoming downstream nucleosome and its translocase subunit Ssl2 (XPB in human TFIIH) drives the rotation of the +1 nucleosome leading to partial detachment of nucleosomal DNA and intimate interactions between TFIIH and the nucleosome. The structures provide insights into how transcription initiation can be influenced by the +1 nucleosome and may explain why the transcription start site is often located roughly 60 base pairs upstream of the dyad of the +1 nucleosome in yeast.

Effects of Defective Unloading and Recycling of PCNA Revealed by the Analysis of ELG1 Mutants

Timely and complete replication of the genome is essential for life. The PCNA ring plays an essential role in DNA replication and repair by contributing to the processivity of DNA polymerases and by recruiting proteins that act in DNA replication-associated processes. The ELG1 gene encodes a protein that works, together with the Rfc2-5 subunits (shared by the replication factor C complex), to unload PCNA from chromatin. While ELG1 is not essential for life, deletion of the gene has strong consequences for the stability of the genome, and elg1 mutants exhibit sensitivity to DNA damaging agents, defects in genomic silencing, high mutation rates, and other striking phenotypes. Here, we sought to understand whether all the roles attributed to Elg1 in genome stability maintenance are due to its effects on PCNA unloading, or whether they are due to additional functions of the protein. By using a battery of mutants that affect PCNA accumulation at various degrees, we show that all the phenotypes measured correlate with the amount of PCNA left at the chromatin. Our results thus demonstrate the importance of Elg1 and of PCNA unloading in promoting proper chromatin structure and in maintaining a stable genome.

Are the epigenetic changes predictive of therapeutic efficacy for psychiatric disorders? A translational approach towards novel drug targets

The etiopathogenesis of mental disorders is not fully understood and accumulating evidence support that clinical symptomatology cannot be assigned to a single gene mutation, but it involves several genetic factors. More specifically, a tight association between genes and environmental risk factors, which could be mediated by epigenetic mechanisms, may play a role in the development of mental disorders. Several data suggest that epigenetic modifications such as DNA methylation, post-translational histone modification and interference of microRNA (miRNA) or long non-coding RNA (lncRNA) may modify the severity of the disease and the outcome of the therapy. Indeed, the study of these mechanisms may help to identify patients particularly vulnerable to mental disorders and may have potential utility as biomarkers to facilitate diagnosis and treatment of psychiatric disorders. This article summarizes the most relevant preclinical and human data showing how epigenetic modifications can be central to the therapeutic efficacy of antidepressant and/or antipsychotic agents, as possible predictor of drugs response.

Structural mechanism of extranucleosomal DNA readout by the INO80 complex

The nucleosomal landscape of chromatin depends on the concerted action of chromatin remodelers. The INO80 remodeler specifically places nucleosomes at the boundary of gene regulatory elements, which is proposed to be the result of an ATP-dependent nucleosome sliding activity that is regulated by extranucleosomal DNA features. Here, we use cryo-electron microscopy and functional assays to reveal how INO80 binds and is regulated by extranucleosomal DNA. Structures of the regulatory A-module bound to DNA clarify the mechanism of linker DNA binding. The A-module is connected to the motor unit via an HSA/post-HSA lever element to chemomechanically couple the motor and linker DNA sensing. Two notable sites of curved DNA recognition by coordinated action of the four actin/actin-related proteins and the motor suggest how sliding by INO80 can be regulated by extranucleosomal DNA features. Last, the structures clarify the recruitment of YY1/Ies4 subunits and reveal deep architectural similarities between the regulatory modules of INO80 and SWI/SNF complexes.

Identification and analysis of putative tRNA genes in baculovirus genomes

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

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

Background

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

Results

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

Conclusions

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

Supplementary Information

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

Manipulating chromatin architecture in C. elegans

BACKGROUND

Nucleosome-mediated chromatin compaction has a direct effect on the accessibility of trans-acting activators and repressors to DNA targets and serves as a primary regulatory agent of genetic expression. Understanding the nature and dynamics of chromatin is fundamental to elucidating the mechanisms and factors that epigenetically regulate gene expression. Previous work has shown that there are three types of canonical sequences that strongly regulate nucleosome positioning and thus chromatin accessibility: putative nucleosome-positioning elements, putative nucleosome-repelling sequences, and homopolymeric runs of A/T. It is postulated that these elements can be used to remodel chromatin in C. elegans. Here we show the utility of such elements in vivo, and the extreme efficacy of a newly discovered repelling sequence, PRS-322.

RESULTS

In this work, we show that it is possible to manipulate nucleosome positioning in C. elegans solely using canonical and putative positioning sequences. We have not only tested previously described sequences such as the Widom 601, but also have tested additional nucleosome-positioning sequences: the Trifonov sequence, putative repelling sequence-322 (PRS-322), and various homopolymeric runs of A and T nucleotides.

CONCLUSIONS

Using each of these types of putative nucleosome-positioning sequences, we demonstrate their ability to alter the nucleosome profile in C. elegans as evidenced by altered nucleosome occupancy and positioning in vivo. Additionally, we show the effect that PRS-322 has on nucleosome-repelling and chromatin remodeling.

FRET-FISH probes chromatin compaction at individual genomic loci in single cells

Chromatin compaction is a key biophysical property that influences multiple DNA transactions. Lack of chromatin accessibility is frequently used as proxy for chromatin compaction. However, we currently lack tools for directly probing chromatin compaction at individual genomic loci. To fill this gap, here we present FRET-FISH, a method combining fluorescence resonance energy transfer (FRET) with DNA fluorescence in situ hybridization (FISH) to probe chromatin compaction at select loci in single cells. We first validate FRET-FISH by comparing it with ATAC-seq, demonstrating that local compaction and accessibility are strongly correlated. FRET-FISH also detects expected differences in compaction upon treatment with drugs perturbing global chromatin condensation. We then leverage FRET-FISH to study local chromatin compaction on the active and inactive X chromosome, along the nuclear radius, in different cell cycle phases, and during increasing passage number. FRET-FISH is a robust tool for probing local chromatin compaction in single cells.

Distinct chromosomal “niches” in the genome of Saccharomyces cerevisiae provide the background for genomic innovation and shape the fate of gene duplicates

Nearly one third of Saccharomyces cerevisiae protein coding sequences correspond to duplicate genes, equally split between small-scale duplicates (SSD) and whole-genome duplicates (WGD). While duplicate genes have distinct properties compared to singletons, to date, there has been no systematic analysis of their positional preferences. In this work, we show that SSD and WGD genes are organized in distinct gene clusters that occupy different genomic regions, with SSD being more peripheral and WGD more centrally positioned close to centromeric chromatin. Duplicate gene clusters differ from the rest of the genome in terms of gene size and spacing, gene expression variability and regulatory complexity, properties that are also shared by singleton genes residing within them. Singletons within duplicate gene clusters have longer promoters, more complex structure and a higher number of protein–protein interactions. Particular chromatin architectures appear to be important for gene evolution, as we find SSD gene-pair co-expression to be strongly associated with the similarity of nucleosome positioning patterns. We propose that specific regions of the yeast genome provide a favourable environment for the generation and maintenance of small-scale gene duplicates, segregating them from WGD-enriched genomic domains. Our findings provide a valuable framework linking genomic innovation with positional genomic preferences.

The TRIPLE PHD FINGERS proteins are required for SWI/SNF complex-mediated +1 nucleosome positioning and transcription start site determination in Arabidopsis

Eukaryotes have evolved multiple ATP-dependent chromatin remodelers to shape the nucleosome landscape. We recently uncovered an evolutionarily conserved SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin remodeler complex in plants reminiscent of the mammalian BAF subclass, which specifically incorporates the MINUSCULE (MINU) catalytic subunits and the TRIPLE PHD FINGERS (TPF) signature subunits. Here we report experimental evidence that establishes the functional relevance of TPF proteins for the complex activity. Our results show that depletion of TPF triggers similar pleiotropic phenotypes and molecular defects to those found in minu mutants. Moreover, we report the genomic location of MINU2 and TPF proteins as representative members of this SWI/SNF complex and their impact on nucleosome positioning and transcription. These analyses unravel the binding of the complex to thousands of genes where it modulates the position of the +1 nucleosome. These targets tend to produce 5′-shifted transcripts in the tpf and minu mutants pointing to the participation of the complex in alternative transcription start site usage. Interestingly, there is a remarkable correlation between +1 nucleosome shift and 5′ transcript length change suggesting their functional connection. In summary, this study unravels the function of a plant SWI/SNF complex involved in +1 nucleosome positioning and transcription start site determination.

Linking transcriptional silencing with chromatin remodeling, folding, and positioning in the nucleus

Chromatin organization is important for many DNA-templated processes in eukaryotic cells such as replication and transcription. Recent studies have uncovered the capacity of epigenetic modifications, phase separation, and nuclear architecture and spatial positioning to regulate chromatin organization in both plants and animals. Here, we provide an overview of the recent progress made in understanding how chromatin is organized within the nucleus at both the local and global levels with respect to the regulation of transcriptional silencing in plants. To be concise while covering important mechanisms across a range of scales, we focus on how epigenetic modifications and chromatin remodelers alter local chromatin structure, how liquid-liquid phase separation physically separates broader chromatin domains into distinct droplets, and how nuclear positioning affects global chromatin organization.

Transcriptional regulator Taf14 binds DNA and is required for the function of transcription factor TFIID in the absence of histone H2A.Z

The transcriptional regulator Taf14 is a component of multiple protein complexes involved in transcription initiation and chromatin remodeling in yeast cells. Although Taf14 is not required for cell viability, it becomes essential in conditions where the formation of the transcription preinitiation complex is hampered. The specific role of Taf14 in mediating transcription initiation and preinitiation complex formation is unclear. Here, we explored its role in the general transcription factor IID by mapping Taf14 genetic and proteomic interactions and found that it was needed for the function of the complex if Htz1, the yeast homolog of histone H2A.Z, was absent from chromatin. Dissecting the functional domains of Taf14 revealed that the linker region between the YEATS and ET domains was required for cell viability in the absence of Htz1 protein. We further show that the linker region of Taf14 interacts with DNA. We propose that providing additional DNA binding capacity might be a general role of Taf14 in the recruitment of protein complexes to DNA and chromatin.

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

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

Mononucleotide A-repeats may Play a Regulatory Role in Endothermic Housekeeping Genes

Background:

Coding and non-coding short tandem repeats (STRs) facilitate a great diversity of phenotypic traits. The imbalance of mononucleotide A-repeats around transcription start sites (TSSs) was found in 3 mammals: H. sapiens, M. musculus, and R. norvegicus.

Principal Findings:

We found that the imbalance pattern originated in some vertebrates. A similar pattern was observed in mammals and birds, but not in amphibians and reptiles. We proposed that the enriched A-repeats upstream of TSSs is a novel hallmark of endotherms or warm-blooded animals. Gene ontology analysis indicates that the primary function of upstream A-repeats involves metabolism, cellular transportation, and sensory perception (smell and chemical stimulus) through housekeeping genes.

Conclusions:

Upstream A-repeats may play a regulatory role in the metabolic process of endothermic animals.

NucleoMap: A computational tool for identifying nucleosomes in ultra-high resolution contact maps

Although poorly positioned nucleosomes are ubiquitous in the eukaryotic genome, they are difficult to identify with existing nucleosome identification methods. Recently available enhanced high-throughput chromatin conformation capture techniques such as Micro-C, DNase Hi-C, and Hi-CO characterize nucleosome-level chromatin proximity, probing the positions of mono-nucleosomes and the spacing between nucleosome pairs at the same time, enabling nucleosome profiling in poorly positioned regions. Here we develop a novel computational approach, NucleoMap, to identify nucleosome positioning from ultra-high resolution chromatin contact maps. By integrating nucleosome read density, contact distances, and binding preferences, NucleoMap precisely locates nucleosomes in both prokaryotic and eukaryotic genomes and outperforms existing nucleosome identification methods in both precision and recall. We rigorously characterize genome-wide association in eukaryotes between the spatial organization of mono-nucleosomes and their corresponding histone modifications, protein binding activities, and higher-order chromatin functions. We also find evidence of two tetra-nucleosome folding structures in human embryonic stem cells and analyze their association with multiple structural and functional regions. Based on the identified nucleosomes, nucleosome contact maps are constructed, reflecting the inter-nucleosome distances and preserving the contact distance profiles in original contact maps.