Nucleosome organization in the Drosophila genome

Emerging Approaches to Profile Accessible Chromatin from Formalin-Fixed Paraffin-Embedded Sections

Nucleosomes are non-uniformly distributed across eukaryotic genomes, with stretches of 'open' chromatin strongly associated with transcriptionally active promoters and enhancers. Understanding chromatin accessibility patterns in normal tissue and how they are altered in pathologies can provide critical insights to development and disease. With the advent of high-throughput sequencing, a variety of strategies have been devised to identify open regions across the genome, including DNase-seq, MNase-seq, FAIRE-seq, ATAC-seq, and NicE-seq. However, the broad application of such methods to FFPE (formalin-fixed paraffin-embedded) tissues has been curtailed by the major technical challenges imposed by highly fixed and often damaged genomic material. Here, we review the most common approaches for mapping open chromatin regions, recent optimizations to overcome the challenges of working with FFPE tissue, and a brief overview of a typical data pipeline with analysis considerations.

Three-step mechanism of promoter escape by RNA polymerase II

The transition from transcription initiation to elongation is highly regulated in human cells but remains incompletely understood at the structural level. In particular, it is unclear how interactions between RNA polymerase II (RNA Pol II) and initiation factors are broken to enable promoter escape. Here, we reconstitute RNA Pol II promoter escape in vitro and determine high-resolution structures of initially transcribing complexes containing 8-, 10-, and 12-nt ordered RNAs and two elongation complexes containing 14-nt RNAs. We suggest that promoter escape occurs in three major steps. First, the growing RNA displaces the B-reader element of the initiation factor TFIIB without evicting TFIIB. Second, the rewinding of the transcription bubble coincides with the eviction of TFIIA, TFIIB, and TBP. Third, the binding of DSIF and NELF facilitates TFIIE and TFIIH dissociation, establishing the paused elongation complex. This three-step model for promoter escape fills a gap in our understanding of the initiation-elongation transition of RNA Pol II transcription.

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.

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.

Interplay between the transcription preinitiation complex and the +1 nucleosome

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

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.

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.

Identifying promoter sequence architectures via a chunking-based algorithm using non-negative matrix factorisation

Core promoters are stretches of DNA at the beginning of genes that contain information that facilitates the binding of transcription initiation complexes. Different functional subsets of genes have core promoters with distinct architectures and characteristic motifs. Some of these motifs inform the selection of transcription start sites (TSS). By discovering motifs with fixed distances from known TSS positions, we could in principle classify promoters into different functional groups. Due to the variability and overlap of architectures, promoter classification is a difficult task that requires new approaches. In this study, we present a new method based on non-negative matrix factorisation (NMF) and the associated software called seqArchR that clusters promoter sequences based on their motifs at near-fixed distances from a reference point, such as TSS. When combined with experimental data from CAGE, seqArchR can efficiently identify TSS-directing motifs, including known ones like TATA, DPE, and nucleosome positioning signal, as well as novel lineage-specific motifs and the function of genes associated with them. By using seqArchR on developmental time courses, we reveal how relative use of promoter architectures changes over time with stage-specific expression. seqArchR is a powerful tool for initial genome-wide classification and functional characterisation of promoters. Its use cases are more general: it can also be used to discover any motifs at near-fixed distances from a reference point, even if they are present in only a small subset of sequences.

HapX-mediated H2B deub1 and SreA-mediated H2A.Z deposition coordinate in fungal iron resistance

Plant pathogens are challenged by host-derived iron starvation or excess during infection, but the mechanism through which pathogens counteract iron stress is unclear. Here, we found that Fusarium graminearum encounters iron excess during the colonization of wheat heads. Deletion of heme activator protein X (FgHapX), siderophore transcription factor A (FgSreA) or both attenuated virulence. Further, we found that FgHapX activates iron storage under iron excess by promoting histone H2B deubiquitination (H2B deub1) at the promoter of the responsible gene. Meanwhile, FgSreA is shown to inhibit genes mediating iron acquisition during iron excess by facilitating the deposition of histone variant H2A.Z and histone 3 lysine 27 trimethylation (H3K27 me3) at the first nucleosome after the transcription start site. In addition, the monothiol glutaredoxin FgGrx4 is responsible for iron sensing and control of the transcriptional activity of FgHapX and FgSreA via modulation of their enrichment at target genes and recruitment of epigenetic regulators, respectively. Taken together, our findings elucidated the molecular mechanisms for adaptation to iron excess mediated by FgHapX and FgSreA during infection in F. graminearum and provide novel insights into regulation of iron homeostasis at the chromatin level in eukaryotes.

SEM: sized-based expectation maximization for characterizing nucleosome positions and subtypes

Genome-wide nucleosome profiles are predominantly characterized using MNase-seq, which involves extensive MNase digestion and size selection to enrich for mono-nucleosome-sized fragments. Most available MNase-seq analysis packages assume that nucleosomes uniformly protect 147bp DNA fragments. However, some nucleosomes with atypical histone or chemical compositions protect shorter lengths of DNA. The rigid assumptions imposed by current nucleosome analysis packages ignore variation in nucleosome lengths, potentially blinding investigators to regulatory roles played by atypical nucleosomes. To enable the characterization of different nucleosome types from MNase-seq data, we introduce the Size-based Expectation Maximization (SEM) nucleosome calling package. SEM employs a hierarchical Gaussian mixture model to estimate the positions and subtype identity of nucleosomes from MNase-seq fragments. Nucleosome subtypes are automatically identified based on the distribution of protected DNA fragment lengths at nucleosome positions. Benchmark analysis indicates that SEM is on par with existing packages in terms of standard nucleosome-calling accuracy metrics, while uniquely providing the ability to characterize nucleosome subtype identities. Using SEM on a low-dose MNase H2B MNase-ChIP-seq dataset from mouse embryonic stem cells, we identified three nucleosome types: short-fragment nucleosomes, canonical nucleosomes, and di-nucleosomes. The short-fragment nucleosomes can be divided further into two subtypes based on their chromatin accessibility. Interestingly, the subset of short-fragment nucleosomes in accessible regions exhibit high MNase sensitivity and display distribution patterns around transcription start sites (TSSs) and CTCF peaks, similar to the previously reported "fragile nucleosomes". These SEM-defined accessible short-fragment nucleosomes are found not just in promoters, but also in enhancers and other regulatory regions. Additional investigations reveal their co-localization with the chromatin remodelers Chd6, Chd8, and Ep400. In summary, SEM provides an effective platform for distinguishing various nucleosome subtypes, paving the way for future exploration of non-standard nucleosomes.

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.

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

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

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

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

Transcriptome profiles and chromatin states in mouse androgenetic haploid embryonic stem cells

Haploid embryonic stem cells (haESCs) are derived from the inner cell mass of the haploid blastocyst, containing only one set of chromosomes. Extensive and accurate chromatin remodelling occurs during haESC derivation, but the intrinsic transcriptome profiles and chromatin structure of haESCs have not been fully explored. We profiled the transcriptomes, nucleosome positioning, and key histone modifications of four mouse haESC lines, and compared these profiles with those of other closely-related stem cell lines, MII oocytes, round spermatids, sperm, and mouse embryonic fibroblasts. haESCs had transcriptome profiles closer to those of naïve pluripotent stem cells. Consistent with the one X chromosome in haESCs, Xist was repressed, indicating no X chromosome inactivation. haESCs and ESCs shared a similar global chromatin structure. However, a nucleosome depletion region was identified in 2056 promoters in ESCs, which was absent in haESCs. Furthermore, three characteristic spatial relationships were formed between transcription factor motifs and nucleosomes in both haESCs and ESCs, specifically in the linker region, on the nucleosome central surface, and nucleosome borders. Furthermore, the chromatin state of 4259 enhancers was off in haESCs but active in ESCs. Functional annotation of these enhancers revealed enrichment in regulation of the cell cycle, a predominantly reported mechanism of haESC self-diploidization. Notably, the transcriptome profiles and chromatin structure of haESCs were highly preserved during passaging but different from those of differentiated cell types.

Aclarubicin stimulates RNA polymerase II elongation at closely spaced divergent promoters

Anthracyclines are a class of widely prescribed anticancer drugs that disrupt chromatin by intercalating into DNA and enhancing nucleosome turnover. To understand the molecular consequences of anthracycline-mediated chromatin disruption, we used Cleavage Under Targets and Tagmentation (CUT&Tag) to profile RNA polymerase II during anthracycline treatment in Drosophila cells. We observed that treatment with the anthracycline aclarubicin leads to elevated levels of RNA polymerase II and changes in chromatin accessibility. We found that promoter proximity and orientation affect chromatin changes during aclarubicin treatment, as closely spaced divergent promoter pairs show greater chromatin changes when compared to codirectionally oriented tandem promoters. We also found that aclarubicin treatment changes the distribution of noncanonical DNA G-quadruplex structures both at promoters and at G-rich pericentromeric repeats. Our work suggests that the cancer-killing activity of aclarubicin is driven by the disruption of nucleosomes and RNA polymerase II.

The NELF pausing checkpoint mediates the functional divergence of Cdk9

Promoter-proximal pausing by RNA Pol II is a rate-determining step in gene transcription that is hypothesized to be a prominent point at which regulatory factors act. The pausing factor NELF is known to induce and stabilize pausing, but not all kinds of pausing are NELF-mediated. Here, we find that NELF-depleted Drosophila melanogaster cells functionally recapitulate the NELF-independent pausing we previously observed in fission yeast (which lack NELF). Critically, only NELF-mediated pausing establishes a strict requirement for Cdk9 kinase activity for the release of paused Pol II into productive elongation. Upon inhibition of Cdk9, cells with NELF efficiently shutdown gene transcription, while in NELF-depleted cells, defective, non-productive transcription continues unabated. By introducing a strict checkpoint for Cdk9, the evolution of NELF was likely critical to enable increased regulation of Cdk9 in higher eukaryotes, as Cdk9 availability can be restricted to limit gene transcription without inducing wasteful, non-productive transcription.

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

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

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

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

Actin-related protein 6 facilitates proneural protein-induced gene activation for rapid neural differentiation

Neurogenesis is initiated by basic helix-loop-helix proneural proteins. Here, we show that Actin-related protein 6 (Arp6), a core component of the H2A.Z exchange complex SWR1, interacts with proneural proteins and is crucial for efficient onset of proneural protein target gene expression. Arp6 mutants exhibit reduced transcription in sensory organ precursors (SOPs) downstream of the proneural protein patterning event. This leads to retarded differentiation and division of SOPs and smaller sensory organs. These phenotypes are also observed in proneural gene hypomorphic mutants. Proneural protein expression is not reduced in Arp6 mutants. Enhanced proneural gene expression fails to rescue retarded differentiation in Arp6 mutants, suggesting that Arp6 acts downstream of or in parallel with proneural proteins. H2A.Z mutants display Arp6-like retardation in SOPs. Transcriptomic analyses demonstrate that loss of Arp6 and H2A.Z preferentially decreases expression of proneural protein-activated genes. H2A.Z enrichment in nucleosomes around the transcription start site before neurogenesis correlates highly with greater activation of proneural protein target genes by H2A.Z. We propose that upon proneural protein binding to E-box sites, H2A.Z incorporation around the transcription start site allows rapid and efficient activation of target genes, promoting rapid neural differentiation.

The Histone Variant H2A.Z C-Terminal Domain Has Locus-Specific Differential Effects on H2A.Z Occupancy and Nucleosome Localization

The incorporation of histone variant H2A.Z into nucleosomes creates specialized chromatin domains that regulate DNA-templated processes, such as gene transcription. In Saccharomyces cerevisiae, the diverging H2A.Z C terminus is thought to provide the H2A.Z exclusive functions. To elucidate the roles of this H2A.Z C terminus genome-wide, we used derivatives in which the C terminus was replaced with the corresponding region of H2A (ZA protein), or the H2A region plus a transcriptional activating peptide (ZA-rII'), with the intent of regenerating the H2A.Z-dependent regulation globally. The distribution of these H2A.Z derivatives indicates that the H2A.Z C-terminal region is crucial for both maintaining the occupation level of H2A.Z and the proper positioning of targeted nucleosomes. Interestingly, the specific contribution on incorporation efficiency versus nucleosome positioning varies enormously depending on the locus analyzed. Specifically, the role of H2A.Z in global transcription regulation relies on its C-terminal region. Remarkably, however, this mostly involves genes without a H2A.Z nucleosome in the promoter. Lastly, we demonstrate that the main chaperone complex which deposits H2A.Z to gene regulatory region (SWR1-C) is necessary to localize all H2A.Z derivatives at their specific loci, indicating that the differential association of these derivatives is not due to impaired interaction with SWR1-C. IMPORTANCE We provide evidence that the Saccharomyces cerevisiae C-terminal region of histone variant H2A.Z can mediate its special function in performing gene regulation by interacting with effector proteins and chaperones. These functional interactions allow H2A.Z not only to incorporate to very specific gene regulatory regions, but also to facilitate the gene expression process. To achieve this, we used a chimeric protein which lacks the native H2A.Z C-terminal region but contains an acidic activating region, a module that is known to interact with components of chromatin-remodeling entities and/or transcription modulators. We reasoned that because this activating region can fulfill the role of the H2A.Z C-terminal region, at least in part, the role of the latter would be to interact with these activating region targets.

Contribution of the histone variant H2A.Z to expression of responsive genes in plants

The histone variant H2A.Z plays a critical role in chromatin-based processes such as transcription, replication, and repair in eukaryotes. Although many H2A.Z-associated processes and features are conserved in plants and animals, a distinguishing feature of plant chromatin is the enrichment of H2A.Z in the bodies of genes that exhibit dynamic expression, particularly in response to differentiation and the environment. Recent work sheds new light on the plant machinery that enables dynamic changes in H2A.Z enrichment and identifies additional chromatin-based pathways that contribute to transcriptional properties of H2A.Z-enriched chromatin. In particular, analysis of a variety of responsive loci reveals a repressive role for H2A.Z in expression of responsive genes and identifies roles for SWR1 and INO80 chromatin remodelers in enabling dynamic regulation of H2A.Z levels and transcription. These studies lay the groundwork for understanding how this ancient histone variant is harnessed by plants to enable responsive and dynamic gene expression (Graphical Abstract).

Aclarubicin stimulates RNA polymerase II elongation at closely spaced divergent promoters

Anthracyclines are a class of widely prescribed anti-cancer drugs that disrupt chromatin by intercalating into DNA and enhancing nucleosome turnover. To understand the molecular consequences of anthracycline-mediated chromatin disruption, we utilized CUT&Tag to profile RNA polymerase II during anthracycline treatment in Drosophila cells. We observed that treatment with the anthracycline aclarubicin leads to elevated levels of elongating RNA polymerase II and changes in chromatin accessibility. We found that promoter proximity and orientation impacts chromatin changes during aclarubicin treatment, as closely spaced divergent promoter pairs show greater chromatin changes when compared to codirectionally-oriented tandem promoters. We also found that aclarubicin treatment changes the distribution of non-canonical DNA G-quadruplex structures both at promoters and at G-rich pericentromeric repeats. Our work suggests that the anti-cancer activity of aclarubicin is driven by the effects of nucleosome disruption on RNA polymerase II, chromatin accessibility and DNA structures.

Databases and prospects of dynamic gene regulation in eukaryotes: A mini review

In eukaryotes, dynamic regulation enables DNA polymerases to catalyze a variety of RNA products in spatial and temporal patterns. Dynamic gene expression is regulated by transcription factors (TFs) and epigenetics (DNA methylation and histone modification). The applications of biochemical technology and high-throughput sequencing enhance the understanding of mechanisms of these regulations and affected genomic regions. To provide a searchable platform for retrieving such metadata, numerous databases have been developed based on the integration of genome-wide maps (e.g., ChIP-seq, whole-genome bisulfite sequencing, RNA-seq, ATAC-seq, DNase-seq, and MNase-seq data) and functionally genomic annotation. In this mini review, we summarize the main functions of TF-related databases and outline the prevalent approaches used in inferring epigenetic regulations, their associated genes, and functions. We review the literature on crosstalk between TF and epigenetic regulation and the properties of non-coding RNA regulation, which are challenging topics that promise to pave the way for advances in database development.

CFA: An explainable deep learning model for annotating the transcriptional roles of cis-regulatory modules based on epigenetic codes

Metazoa gene expression is controlled by modular DNA segments called cis-regulatory modules (CRMs). CRMs can convey promoter/enhancer/insulator roles, generating additional regulation layers in transcription. Experiments for understanding CRM roles are low-throughput and costly. Large-scale CRM function investigation still depends on computational methods. However, existing in silico tools only recognize enhancers or promoters exclusively, thus accumulating errors when considering CRM promoter/enhancer/insulator roles altogether. Currently, no algorithm can concurrently consider these CRM roles. In this research, we developed the CRM Function Annotator (CFA) model. CFA provides complete CRM transcriptional role labeling based on epigenetic profiling interpretation. We demonstrated that CFA achieves high performance (test macro auROC/auPRC = 94.1%/90.3%) and outperforms existing tools in promoter/enhancer/insulator identification. CFA is also inspected to recognize explainable epigenetic codes consistent with previous findings when labeling CRM roles. By considering the higher-order combinations of the epigenetic codes, CFA significantly reduces false-positive rates in CRM transcriptional role annotation. CFA is available at https://github.com/cobisLab/CFA/.

Histone variants and chromatin structure, update of advances

Histone proteins are highly conserved among all eukaryotes. They have two important functions in the cell: to package the genomic DNA and to regulate gene accessibility. Fundamental to these functions is the ability of histone proteins to interact with DNA and to form the nucleoprotein complex called chromatin. One of the mechanisms the cells use to regulate chromatin and gene expression is through replacing canonical histones with their variants at specific loci to achieve functional consequence. Recent cryo-electron microscope (cryo-EM) studies of chromatin containing histone variants reveal new details that shed light on how variant-specific features influence the structures and functions of chromatin. In this article, we review the current state of knowledge on histone variants biochemistry and discuss the implication of these new structural information on histone variant biology and their functions in transcription.

Structures of +1 nucleosome-bound PIC-Mediator complex

RNA polymerase II-mediated eukaryotic transcription starts with the assembly of the preinitiation complex (PIC) on core promoters. The +1 nucleosome is well positioned about 40 base pairs downstream of the transcription start site (TSS) and is commonly known as a barrier of transcription. The +1 nucleosome-bound PIC-Mediator structures show that PIC-Mediator prefers binding to T40N nucleosome located 40 base pairs downstream of TSS and contacts T50N but not the T70N nucleosome. The nucleosome facilitates the organization of PIC-Mediator on the promoter by binding TFIIH subunit p52 and Mediator subunits MED19 and MED26 and may contribute to transcription initiation. PIC-Mediator exhibits multiple nucleosome-binding patterns, supporting a structural role of the +1 nucleosome in the coordination of PIC-Mediator assembly. Our study reveals the molecular mechanism of PIC-Mediator organization on chromatin and underscores the significance of the +1 nucleosome in regulating transcription initiation.

Developmental and housekeeping transcriptional programs in Drosophila require distinct chromatin remodelers

Gene transcription is a highly regulated process in all animals. In Drosophila, two major transcriptional programs, housekeeping and developmental, have promoters with distinct regulatory compatibilities and nucleosome organization. However, it remains unclear how the differences in chromatin structure relate to the distinct regulatory properties and which chromatin remodelers are required for these programs. Using rapid degradation of core remodeler subunits in Drosophila melanogaster S2 cells, we demonstrate that developmental gene transcription requires SWI/SNF-type complexes, primarily to maintain distal enhancer accessibility. In contrast, wild-type-level housekeeping gene transcription requires the Iswi and Ino80 remodelers to maintain nucleosome positioning and phasing at promoters. These differential remodeler dependencies relate to different DNA-sequence-intrinsic nucleosome affinities, which favor a default ON state for housekeeping but a default OFF state for developmental gene transcription. Overall, our results demonstrate how different transcription-regulatory strategies are implemented by DNA sequence, chromatin structure, and remodeler activity.

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.

Chromatin structure meets cryo-EM: Dynamic building blocks of the functional architecture

Chromatin is a dynamic molecular complex composed of DNA and proteins that package the DNA in the nucleus of eukaryotic cells. The basic structural unit of chromatin is the nucleosome core particle, composed of ~150 base pairs of genomic DNA wrapped around a histone octamer containing two copies each of four histones, H2A, H2B, H3, and H4. Individual nucleosome core particles are connected by short linker DNAs, forming a nucleosome array known as a beads-on-a-string fiber. Higher-order structures of chromatin are closely linked to nuclear events such as replication, transcription, recombination, and repair. Recently, a variety of chromatin structures have been determined by single-particle cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET), and their structural details have provided clues about the chromatin architecture functions in the cell. In this review, we highlight recent cryo-EM structural studies of a fundamental chromatin unit to clarify the functions of chromatin.

Chromatin remodeler Dmp18 regulates apoptosis by controlling H2Av incorporation in Drosophila imaginal disc development

Programmed Cell Death (PCD) or apoptosis is a highly conserved biological process and plays essential roles both in the development and stress context. In Drosophila, expression of pro-apoptotic genes, including reaper (rpr), head involution defective (hid), grim, and sickle (skl), is sufficient to induce cell death. Here, we demonstrate that the chromatin remodeler Dmp18, the homolog of mammalian Znhit1, plays a crucial role in regulating apoptosis in eye and wing development. We showed that loss of Dmp18 disrupted eye and wing development, up-regulated transcription of pro-apoptotic genes, and induced apoptosis. Inhibition of apoptosis suppressed the eye defects caused by Dmp18 deletion. Furthermore, loss of Dmp18 disrupted H2Av incorporation into chromatin, promoted H3K4me3, but reduced H3K27me3 modifications on the TSS regions of pro-apoptotic genes. These results indicate that Dmp18 negatively regulates apoptosis by mediating H2Av incorporation and histone H3 modifications at pro-apoptotic gene loci for transcriptional regulation. Our study uncovers the role of Dmp18 in regulating apoptosis in Drosophila eye and wing development and provides insights into chromatin remodeling regulating apoptosis at the epigenetic levels.

Stochastic Variation in DNA Methylation Modulates Nucleosome Occupancy and Alternative Splicing in Arabidopsis thaliana

Plants use complex gene regulatory mechanisms to overcome diverse environmental challenges. For instance, cold stress induces rapid and massive transcriptome changes via alternative splicing (AS) to confer cold tolerance in plants. In mammals, mounting evidence suggests chromatin structure can regulate co-transcriptional AS. Recent evidence also supports co-transcriptional regulation of AS in plants, but how dynamic changes in DNA methylation and the chromatin structure influence the AS process upon cold stress remains poorly understood. In this study, we used the DNA methylation inhibitor 5-Aza-2′-Deoxycytidine (5-aza-dC) to investigate the role of stochastic variations in DNA methylation and nucleosome occupancy in modulating cold-induced AS, in Arabidopsis thaliana (Arabidopsis). Our results demonstrate that 5-aza-dC derived stochastic hypomethylation modulates nucleosome occupancy and AS profiles of genes implicated in RNA metabolism, plant hormone signal transduction, and of cold-related genes in response to cold stress. We also demonstrate that cold-induced remodelling of DNA methylation regulates genes involved in amino acid metabolism. Collectively, we demonstrate that sudden changes in DNA methylation via drug treatment can influence nucleosome occupancy levels and modulate AS in a temperature-dependent manner to regulate plant metabolism and physiological stress adaptation.

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

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

The Novel Protease Activities of JMJD5–JMJD6–JMJD7 and Arginine Methylation Activities of Arginine Methyltransferases Are Likely Coupled

The surreptitious discoveries of the protease activities on arginine-methylated targets of a subfamily of Jumonji domain-containing family including JMJD5, JMJD6, and JMJD7 pose several questions regarding their authenticity, function, purpose, and relations with others. At the same time, despite several decades of efforts and massive accumulating data regarding the roles of the arginine methyltransferase family (PRMTs), the exact function of this protein family still remains a mystery, though it seems to play critical roles in transcription regulation, including activation and inactivation of a large group of genes, as well as other biological activities. In this review, we aim to elucidate that the function of JMJD5/6/7 and PRMTs are likely coupled. Besides roles in the regulation of the biogenesis of membrane-less organelles in cells, they are major players in regulating stimulating transcription factors to control the activities of RNA Polymerase II in higher eukaryotes, especially in the animal kingdom. Furthermore, we propose that arginine methylation by PRMTs could be a ubiquitous action marked for destruction after missions by a subfamily of the Jumonji protein family.

Large‐scale analysis of Drosophila core promoter function using synthetic promoters

The core promoter plays a central role in setting metazoan gene expression levels, but how exactly it “computes” expression remains poorly understood. To dissect its function, we carried out a comprehensive structure–function analysis in Drosophila. First, we performed a genome‐wide bioinformatic analysis, providing an improved picture of the sequence motifs architecture. We then measured synthetic promoters’ activities of ~3,000 mutational variants with and without an external stimulus (hormonal activation), at large scale and with high accuracy using robotics and a dual luciferase reporter assay. We observed a strong impact on activity of the different types of mutations, including knockout of individual sequence motifs and motif combinations, variations of motif strength, nucleosome positioning, and flanking sequences. A linear combination of the individual motif features largely accounts for the combinatorial effects on core promoter activity. These findings shed new light on the quantitative assessment of gene expression in metazoans.

Polycomb repressive complex 1 initiates and maintains tailless repression in Drosophila embryo

Maternally-deposited morphogens specify the fates of embryonic cells via hierarchically regulating the expression of zygotic genes that encode various classes of developmental regulators. Once the cell fates are determined, Polycomb-group proteins frequently maintain the repressed state of the genes. This study investigates how Polycomb-group proteins repress the expression of tailless, which encodes a developmental regulator in Drosophila embryo. Previous studies have shown that maternal Tramtrack69 facilitates maternal GAGA-binding factor and Heat shock factor binding to the torso response element (tor-RE) to initiate tailless repression in the stage-4 embryo. Chromatin-immunoprecipitation and genetic-interaction studies exhibit that maternally-deposited Polycomb repressive complex 1 (PRC1) recruited by the tor-RE-associated Tramtrack69 represses tailless expression in the stage-4 embryo. A noncanonical Polycomb-group response element (PRE) is mapped to the tailless proximal region. High levels of Bric-a-brac, Tramtrack, and Broad (BTB)-domain proteins are fundamental for maintaining tailless repression in the stage-8 to -10 embryos. Trmtrack69 sporadically distributes in the linear BTB-domain oligomer, which recruits and retains a high level of PRC1 near the GCCAT cluster for repressing tll expression in the stage-14 embryos. Disrupting the retention of PRC1 decreases the levels of PRC1 and Pleiohomeotic protein substantially on the PRE and causes tailless derepression in the stage-14 embryo. Furthermore, the retained PRC1 potentially serves as a second foundation for assembling the well-characterized polymer of the Sterile alpha motif domain in Polyhomeotic protein, which compacts chromatin to maintain the repressed state of tailless in the embryos after stage 14.

regCNN: identifying Drosophila genome-wide cis-regulatory modules via integrating the local patterns in epigenetic marks and transcription factor binding motifs

Transcription regulation in metazoa is controlled by the binding events of transcription factors (TFs) or regulatory proteins on specific modular DNA regulatory sequences called cis-regulatory modules (CRMs). Understanding the distributions of CRMs on a genomic scale is essential for constructing the metazoan transcriptional regulatory networks that help diagnose genetic disorders. While traditional reporter-assay CRM identification approaches can provide an in-depth understanding of functions of some CRM, these methods are usually cost-inefficient and low-throughput. It is generally believed that by integrating diverse genomic data, reliable CRM predictions can be made. Hence, researchers often first resort to computational algorithms for genome-wide CRM screening before specific experiments. However, current existing in silico methods for searching potential CRMs were restricted by low sensitivity, poor prediction accuracy, or high computation time from TFBS composition combinatorial complexity. To overcome these obstacles, we designed a novel CRM identification pipeline called regCNN by considering the base-by-base local patterns in TF binding motifs and epigenetic profiles. On the test set, regCNN shows an accuracy/auROC of 84.5%/92.5% in CRM identification. And by further considering local patterns in epigenetic profiles and TF binding motifs, it can accomplish 4.7% (92.5%–87.8%) improvement in the auROC value over the average value-based pure multi-layer perceptron model. We also demonstrated that regCNN outperforms all currently available tools by at least 11.3% in auROC values. Finally, regCNN is verified to be robust against its resizing window hyperparameter in dealing with the variable lengths of CRMs. The model of regCNN can be downloaded athttp://cobisHSS0.im.nuk.edu.tw/regCNN/.

DeNOPA: decoding nucleosome positions sensitively with sparse ATAC-seq data

As the basal bricks, the dynamics and arrangement of nucleosomes orchestrate the higher architecture of chromatin in a fundamental way, thereby affecting almost all nuclear biology processes. Thanks to its rather simple protocol, assay for transposase-accessible chromatin using sequencing (ATAC)-seq has been rapidly adopted as a major tool for chromatin-accessible profiling at both bulk and single-cell levels; however, to picture the arrangement of nucleosomes per se remains a challenge with ATAC-seq. In the present work, we introduce a novel ATAC-seq analysis toolkit, named decoding nucleosome organization profile based on ATAC-seq data (deNOPA), to predict nucleosome positions. Assessments showed that deNOPA outperformed state-of-the-art tools with ultra-sparse ATAC-seq data, e.g. no more than 0.5 fragment per base pair. The remarkable performance of deNOPA was fueled by the short fragment reads, which compose nearly half of sequenced reads in the ATAC-seq libraries and are commonly discarded by state-of-the-art nucleosome positioning tools. However, we found that the short fragment reads enrich information on nucleosome positions and that the linker regions were predicted by reads from both short and long fragments using Gaussian smoothing. Last, using deNOPA, we showed that the dynamics of nucleosome organization may not directly couple with chromatin accessibility in the cis-regulatory regions when human cells respond to heat shock stimulation. Our deNOPA provides a powerful tool with which to analyze the dynamics of chromatin at nucleosome position level with ultra-sparse ATAC-seq data.

Structural basis of chromatin regulation by histone variant H2A.Z

The importance of histone variant H2A.Z in transcription regulation has been well established, yet its mechanism-of-action remains enigmatic. Conflicting evidence exists in support of both an activating and a repressive role of H2A.Z in transcription. Here we report cryo-electron microscopy (cryo-EM) structures of nucleosomes and chromatin fibers containing H2A.Z and those containing canonical H2A. The structures show that H2A.Z incorporation results in substantial structural changes in both nucleosome and chromatin fiber. While H2A.Z increases the mobility of DNA terminus in nucleosomes, it simultaneously enables nucleosome arrays to form a more regular and condensed chromatin fiber. We also demonstrated that H2A.Z’s ability to enhance nucleosomal DNA mobility is largely attributed to its characteristic shorter C-terminus. Our study provides the structural basis for H2A.Z-mediated chromatin regulation, showing that the increase flexibility of the DNA termini in H2A.Z nucleosomes is central to its dual-functions in chromatin regulation and in transcription.

High contiguity de novo genome assembly and DNA modification analyses for the fungus fly, Sciara coprophila, using single-molecule sequencing

BACKGROUND

The lower Dipteran fungus fly, Sciara coprophila, has many unique biological features that challenge the rule of genome DNA constancy. For example, Sciara undergoes paternal chromosome elimination and maternal X chromosome nondisjunction during spermatogenesis, paternal X elimination during embryogenesis, intrachromosomal DNA amplification of DNA puff loci during larval development, and germline-limited chromosome elimination from all somatic cells. Paternal chromosome elimination in Sciara was the first observation of imprinting, though the mechanism remains a mystery. Here, we present the first draft genome sequence for Sciara coprophila to take a large step forward in addressing these features.

RESULTS

We assembled the Sciara genome using PacBio, Nanopore, and Illumina sequencing. To find an optimal assembly using these datasets, we generated 44 short-read and 50 long-read assemblies. We ranked assemblies using 27 metrics assessing contiguity, gene content, and dataset concordance. The highest-ranking assemblies were scaffolded using BioNano optical maps. RNA-seq datasets from multiple life stages and both sexes facilitated genome annotation. A set of 66 metrics was used to select the first draft assembly for Sciara. Nearly half of the Sciara genome sequence was anchored into chromosomes, and all scaffolds were classified as X-linked or autosomal by coverage.

CONCLUSIONS

We determined that X-linked genes in Sciara males undergo dosage compensation. An entire bacterial genome from the Rickettsia genus, a group known to be endosymbionts in insects, was co-assembled with the Sciara genome, opening the possibility that Rickettsia may function in sex determination in Sciara. Finally, the signal level of the PacBio and Nanopore data support the presence of cytosine and adenine modifications in the Sciara genome, consistent with a possible role in imprinting.

Epigenetic differences in an identical genetic background modulate alternative splicing in A. thaliana

How stable and temperature-dependent variations in DNA methylation and nucleosome occupancy influence alternative splicing (AS) remains poorly understood in plants. To answer this, we generated transcriptome, whole-genome bisulfite, and MNase sequencing data for an epigenetic Recombinant Inbred Line (epiRIL) of A. thaliana at normal and cold temperature. For comparative analysis, the same data sets for the parental ecotype Columbia (Col-0) were also generated, whereas for DNA methylation, previously published high confidence methylation profiles of Col-0 were used. Significant epigenetic differences in an identical genetic background were observed between Col-0 and epiRIL lines under normal and cold temperatures. Our transcriptome data revealed that differential DNA methylation and nucleosome occupancy modulate expression levels of many genes and AS in response to cold. Collectively, DNA methylation and nucleosome levels exhibit characteristic patterns around intron-exon boundaries at normal and cold conditions, and any perturbation in them, in an identical genetic background is sufficient to modulate AS in Arabidopsis.

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

Background

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

Results

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

Conclusions

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

Coupling between Sequence-Mediated Nucleosome Organization and Genome Evolution

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

Genome information processing by the INO80 chromatin remodeler positions nucleosomes

The fundamental molecular determinants by which ATP-dependent chromatin remodelers organize nucleosomes across eukaryotic genomes remain largely elusive. Here, chromatin reconstitutions on physiological, whole-genome templates reveal how remodelers read and translate genomic information into nucleosome positions. Using the yeast genome and the multi-subunit INO80 remodeler as a paradigm, we identify DNA shape/mechanics encoded signature motifs as sufficient for nucleosome positioning and distinct from known DNA sequence preferences of histones. INO80 processes such information through an allosteric interplay between its core- and Arp8-modules that probes mechanical properties of nucleosomal and linker DNA. At promoters, INO80 integrates this readout of DNA shape/mechanics with a readout of co-evolved sequence motifs via interaction with general regulatory factors bound to these motifs. Our findings establish a molecular mechanism for robust and yet adjustable +1 nucleosome positioning and, more generally, remodelers as information processing hubs that enable active organization and allosteric regulation of the first level of chromatin.

DNA sequence preferences or statistical positioning of histones has not explained genomic patterns of nucleosome organisation in vivo. Here, the authors establish DNA shape/mechanics as key elements that have evolved together with binding sites of DNA sequence-specific barriers so that such information directs nucleosome positioning by chromatin remodelers.

Analysis of Local Chromatin States Reveals Gene Transcription Potential during Mouse Neural Progenitor Cell Differentiation

Chromatin dynamics play a critical role in cell fate determination and maintenance by regulating the expression of genes essential for development and differentiation. In mouse embryonic stem cells (mESCs), maintenance of pluripotency coincides with a poised chromatin state containing active and repressive histone modifications. However, the structural features of poised chromatin are largely uncharacterized. By adopting mild time-course MNase-seq with computational analysis, the low-compact chromatin in mESCs is featured in two groups: one in more open regions, corresponding to an active state, and the other enriched with bivalent histone modifications, considered the poised state. A parameter called the chromatin opening potential index (COPI) is also devised to quantify the transcription potential based on the dynamic changes of MNase-seq signals at promoter regions. Use of COPI provides effective prediction of gene activation potential and, more importantly, reveals a few developmental factors essential for mouse neural progenitor cell (NPC) differentiation.

Multiple roles of H2A.Z in regulating promoter chromatin architecture in human cells

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

The nuclear cap-binding complex as choreographer of gene transcription and pre-mRNA processing

In this review, Rambout and Maquat discuss known roles of the nuclear cap-binding complex (CBC) during the transcription of genes that encode proteins, stitching together past studies from diverse groups to describe the continuum of CBC-mediated checks and balances in eukaryotic cells.

The largely nuclear cap-binding complex (CBC) binds to the 5′ caps of RNA polymerase II (RNAPII)-synthesized transcripts and serves as a dynamic interaction platform for a myriad of RNA processing factors that regulate gene expression. While influence of the CBC can extend into the cytoplasm, here we review the roles of the CBC in the nucleus, with a focus on protein-coding genes. We discuss differences between CBC function in yeast and mammals, covering the steps of transcription initiation, release of RNAPII from pausing, transcription elongation, cotranscriptional pre-mRNA splicing, transcription termination, and consequences of spurious transcription. We describe parameters known to control the binding of generic or gene-specific cofactors that regulate CBC activities depending on the process(es) targeted, illustrating how the CBC is an ever-changing choreographer of gene expression.

Computational methods and next-generation sequencing approaches to analyze epigenetics data: Profiling of methods and applications

Epigenetics is mainly comprised of features that regulate genomic interactions thereby playing a crucial role in a vast array of biological processes. Epigenetic mechanisms such as DNA methylation and histone modifications influence gene expression by modulating the packaging of DNA in the nucleus. A plethora of studies have emphasized the importance of analyzing epigenetics data through genome-wide studies and high-throughput approaches, thereby providing key insights towards epigenetics-based diseases such as cancer. Recent advancements have been made towards translating epigenetics research into a high throughput approach such as genome-scale profiling. Amongst all, bioinformatics plays a pivotal role in achieving epigenetics-related computational studies. Despite significant advancements towards epigenomic profiling, it is challenging to understand how various epigenetic modifications such as chromatin modifications and DNA methylation regulate gene expression. Next-generation sequencing (NGS) provides accurate and parallel sequencing thereby allowing researchers to comprehend epigenomic profiling. In this review, we summarize different computational methods such as machine learning and other bioinformatics tools, publicly available databases and resources to identify key modifications associated with epigenetic machinery. Additionally, the review also focuses on understanding recent methodologies related to epigenome profiling using NGS methods ranging from library preparation, different sequencing platforms and analytical techniques to evaluate various epigenetic modifications such as DNA methylation and histone modifications. We also provide detailed information on bioinformatics tools and computational strategies responsible for analyzing large scale data in epigenetics.

Regulation of Promoter Proximal Pausing of RNA Polymerase II in Metazoans

Regulation of transcription is a tightly choreographed process. The establishment of RNA polymerase II promoter proximal pausing soon after transcription initiation and the release of Pol II into productive elongation are key regulatory processes that occur in early elongation. We describe the techniques and tools that have become available for the study of promoter proximal pausing and their utility for future experiments. We then provide an overview of the factors and interactions that govern a multipartite pausing process and address emerging questions surrounding the mechanism of RNA polymerase II's subsequent advancement into the gene body. Finally, we address remaining controversies and future areas of study.

Basis of specificity for a conserved and promiscuous chromatin remodeling protein

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