MNase titration reveals differences between nucleosome occupancy and chromatin accessibility

Chromatin and gene regulation in archaea

The chromatinisation of DNA by nucleoid-associated proteins (NAPs) in archaea 'formats' the genome structure in profound ways, revealing both striking differences and analogies to eukaryotic chromatin. However, the extent to which archaeal NAPs actively regulate gene expression remains poorly understood. The dawn of quantitative chromatin mapping techniques and first NAP-specific occupancy profiles in different archaea promise a more accurate view. A picture emerges where in diverse archaea with very different NAP repertoires chromatin maintains access to regulatory motifs including the gene promoter independently of transcription activity. Our re-analysis of genome-wide occupancy data of the crenarchaeal NAP Cren7 shows that these chromatin-free regions are flanked by increased Cren7 binding across the transcription start site. While bacterial NAPs often form heterochromatin-like regions across islands with xenogeneic genes that are transcriptionally silenced, there is little evidence for similar structures in archaea and data from Haloferax show that the promoters of xenogeneic genes remain accessible. Local changes in chromatinisation causing wide-ranging effects on transcription restricted to one chromosomal interaction domain (CID) in Saccharolobus islandicus hint at a higher-order level of organisation between chromatin and transcription. The emerging challenge is to integrate results obtained at microscale and macroscale, reconciling molecular structure and function with dynamic genome-wide chromatin landscapes.

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

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

Genome wide nucleosome landscape shapes 3D chromatin organization

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

Recent advances in exploring transcriptional regulatory landscape of crops

Crop breeding entails developing and selecting plant varieties with improved agronomic traits. Modern molecular techniques, such as genome editing, enable more efficient manipulation of plant phenotype by altering the expression of particular regulatory or functional genes. Hence, it is essential to thoroughly comprehend the transcriptional regulatory mechanisms that underpin these traits. In the multi-omics era, a large amount of omics data has been generated for diverse crop species, including genomics, epigenomics, transcriptomics, proteomics, and single-cell omics. The abundant data resources and the emergence of advanced computational tools offer unprecedented opportunities for obtaining a holistic view and profound understanding of the regulatory processes linked to desirable traits. This review focuses on integrated network approaches that utilize multi-omics data to investigate gene expression regulation. Various types of regulatory networks and their inference methods are discussed, focusing on recent advancements in crop plants. The integration of multi-omics data has been proven to be crucial for the construction of high-confidence regulatory networks. With the refinement of these methodologies, they will significantly enhance crop breeding efforts and contribute to global food security.

Structural dynamics in chromatin unraveling by pioneer transcription factors

Pioneer transcription factors are proteins with a dual function. First, they regulate transcription by binding to nucleosome-free DNA regulatory elements. Second, they bind to DNA while wrapped around histone proteins in the chromatin and mediate chromatin opening. The molecular mechanisms that connect the two functions are yet to be discovered. In recent years, pioneer factors received increased attention mainly because of their crucial role in promoting cell fate transitions that could be used for regenerative therapies. For example, the three factors required to induce pluripotency in somatic cells, Oct4, Sox2, and Klf4 were classified as pioneer factors and studied extensively. With this increased attention, several structures of complexes between pioneer factors and chromatin structural units (nucleosomes) have been resolved experimentally. Furthermore, experimental and computational approaches have been designed to study two unresolved, key scientific questions: First, do pioneer factors induce directly local opening of nucleosomes and chromatin fibers upon binding? And second, how do the unstructured tails of the histones impact the structural dynamics involved in such conformational transitions? Here we review the current knowledge about transcription factor-induced nucleosome dynamics and the role of the histone tails in this process. We discuss what is needed to bridge the gap between the static views obtained from the experimental structures and the key structural dynamic events in chromatin opening. Finally, we propose that integrating nuclear magnetic resonance spectroscopy with molecular dynamics simulations is a powerful approach to studying pioneer factor-mediated dynamics of nucleosomes and perhaps small chromatin fibers using native DNA sequences.

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.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

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

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

A data-adaptive methods in detecting exogenous methyltransferase accessible chromatin in human genome using nanopore sequencing

MOTIVATION

Identifying chromatin accessibility is one of the key steps in studying the regulation of eukaryotic genomes. The combination of exogenous methyltransferase and nanopore sequencing provides an strategy to identify open chromatin over long genomic ranges at the single-molecule scale. However, endogenous methylation, non-open-chromatin-specific exogenous methylation and base-calling errors limit the accuracy and hinders its application to complex genomes.

RESULTS

We systematically evaluated the impact of these three influence factors, and developed a model-based computational method, methyltransferase accessible genome region finder (MAGNIFIER), to address the issues. By incorporating control data, MAGNIFIER attenuates the three influence factors with data-adaptive comparison strategy. We demonstrate that MAGNIFIER is not only sensitive to identify the open chromatin with much improved accuracy, but also able to detect the chromatin accessibility of repetitive regions that are missed by NGS-based methods. By incorporating long-read RNA-seq data, we revealed the association between the accessible Alu elements and non-classic gene isoforms.

AVAILABILITY AND IMPLEMENTATION

Freely available on web at https://github.com/Goatofmountain/MAGNIFIER.

Skeletal muscle differentiation induces wide-ranging nucleosome repositioning in muscle gene promoters

In a previous report, we demonstrated that Cbx1, PurB and Sp3 inhibited cardiac muscle differentiation by increasing nucleosome density around cardiac muscle gene promoters. Since cardiac and skeletal muscle express many of the same proteins, we asked if Cbx1, PurB and Sp3 similarly regulated skeletal muscle differentiation. In a C2C12 model of skeletal muscle differentiation, Cbx1 and PurB knockdown increased myotube formation. In contrast, Sp3 knockdown inhibited myotube formation, suggesting that Sp3 played opposing roles in cardiac muscle and skeletal muscle differentiation. Consistent with this finding, Sp3 knockdown also inhibited various muscle-specific genes. The Cbx1, PurB and Sp3 proteins are believed to influence gene-expression in part by altering nucleosome position. Importantly, we developed a statistical approach to determine if changes in nucleosome positioning were significant and applied it to understanding the architecture of muscle-specific genes. Through this novel statistical approach, we found that during myogenic differentiation, skeletal muscle-specific genes undergo a set of unique nucleosome changes which differ significantly from those shown in commonly expressed muscle genes. While Sp3 binding was associated with nucleosome loss, there appeared no correlation with the aforementioned nucleosome changes. In summary, we have identified a novel role for Sp3 in skeletal muscle differentiation and through the application of quantifiable MNase-seq have discovered unique fingerprints of nucleosome changes for various classes of muscle genes during myogenic differentiation.

The multiple arrays of a PD1-derived peptide on chromatin specifically bind to PD-L1 and induce doxorubicin resistance in cancer cell lines

Programmed cell death-1 (PD-1) and its ligand 1 (PD-L1) have been extensively studied to block the activation of the PD-1 axis immune checkpoint pathway on T cells. However, recent evidence suggests that PD-1-independent PD-L1 pathways in cancer cells and macrophages are important in cellular proliferation and survival of those cell types but the underlying mechanism is little understood. To recapitulate the binding of multiple PD-1 to the high levels of PD-L1 expression on cancer cell membranes, recombinant histones carrying a PD-1 receptor-derived peptide (PDR) were prepared and used to assemble a PDR-displayed nucleosome array. PDR-displayed chromatin showed physiologically spaced nucleosomes, unaffected by N-terminal histone tails and biotinylated DNA modifications. PDR-displayed chromatin exhibited selective binding to the breast cancer cell line with high PD-L1 expression, MDA-MB-231, where saturated binding occurred within 3 h, comparable to well-known antigen-antibody reactions. Notably, PDR-displayed chromatin enhanced resistance to doxorubicin in PD-L1-overexpressed cancer cells, revealing a PD-L1 level-dependent effect on cell survival upon chemotherapy. Our study sheds light on the potential of PDR-displayed chromatin as a tool to induce chemoresistance in cancer cells without employing anti-PD-L1 antibodies or soluble PD-1 receptors. Further investigation into the underlying signaling pathways and therapeutic applications is warranted, positioning PDR-displayed chromatin as a valuable tool for understanding PD-L1-mediated intracellular signaling pathways in cancer cells with high PD-L1 expression.

Regulation of chromatin architecture by transcription factor binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome to control the transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers at nucleosome resolution to determine TF effect on chromatin architecture and shed light into its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. At the same time, the level of linker histone and tail acetylation, or different linker DNA lengths, regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

Noncanonical Activity of Med4 as a Gatekeeper of Metastasis through Epigenetic Control of Integrin Signaling

Breast cancer metastatic relapse after a latency period, known as metastatic dormancy. Through genetic screening in mice, we identified the mediator complex subunit 4 (Med4) as a novel tumor-cell intrinsic gatekeeper in metastatic reactivation. Med4 downregulation effectively awakened dormant breast cancer cells, prompting macroscopic metastatic outgrowth in the lungs. Med4 depletion results in profound changes in nuclear size and three-dimensional chromatin architecture from compacted to relaxed states in contrast to the canonical function of the Mediator complex. These changes rewire the expression of extracellular matrix proteins, integrins, and signaling components resulting in integrin-mediated mechano-transduction and activation of YAP and MRTF. The assembly of stress fibers pulls on the nuclear membrane and contributes to reinforcing the overall chromatin modifications by Med4 depletion. MED4 gene deletions were observed in patients with metastatic breast cancer, and reduced MED4 expression correlates with worse prognosis, highlighting its significance as a potential biomarker for recurrence.

The transcriptional activator Klf5 recruits p300-mediated H3K27ac for maintaining trophoblast stem cell pluripotency

The effective proliferation and differentiation of trophoblast stem cells (TSCs) is indispensable for the development of the placenta, which is the key to maintaining normal fetal growth during pregnancy. Kruppel-like factor 5 (Klf5) is implicated in the activation of pluripotency gene expression in embryonic stem cells (ESCs), yet its function in TSCs is poorly understood. Here, we showed that Klf5 knockdown resulted in the downregulation of core TSC-specific genes, consequently causing rapid differentiation of TSCs. Consistently, Klf5-depleted embryos lost the ability to establish TSCs in vitro. At the molecular level, Klf5 preferentially occupied the proximal promoter regions and maintained an open chromatin architecture of key TSC-specific genes. Deprivation of Klf5 impaired the enrichment of p300, a major histone acetyl transferase of H3 lysine 27 acetylation (H3K27ac), and further reduced the occupancy of H3K27ac at promoter regions, leading to decreased transcriptional activity of TSC pluripotency genes. Thus, our findings highlight a novel mechanism of Klf5 in regulating the self-renewal and differentiation of TSCs and provide a reference for understanding placental development and improving pregnancy rates.

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.

Analysis of HBV cccDNA Minichromosome Accessibility by MNase-qPCR and High-Throughput Sequencing

The covalently closed circular DNA (cccDNA) of the hepatitis B virus (HBV) is organized as a minichromosome structure in the nucleus of infected hepatocytes and considered the major obstacle to the discovery of a cure for HBV. Until now, no strategies directly targeting cccDNA have been advanced to clinical stages as much is unknown about the accessibility and activity regulation of the cccDNA minichromosome. We have described the method for evaluation of the cccDNA minichromosome accessibility using micrococcal nuclease-quantitative polymerase chain reaction and high-throughput sequencing, which could be useful tools for cccDNA research and HBV cure studies.

Expanding the list of sequence-agnostic enzymes for chromatin conformation capture assays with S1 nuclease

This study presents a novel approach for mapping global chromatin interactions using S1 nuclease, a sequence-agnostic enzyme. We develop and outline a protocol that leverages S1 nuclease's ability to effectively introduce breaks into both open and closed chromatin regions, allowing for comprehensive profiling of chromatin properties. Our S1 Hi-C method enables the preparation of high-quality Hi-C libraries, marking a significant advancement over previously established DNase I Hi-C protocols. Moreover, S1 nuclease's capability to fragment chromatin to mono-nucleosomes suggests the potential for mapping the three-dimensional organization of the genome at high resolution. This methodology holds promise for an improved understanding of chromatin state-dependent activities and may facilitate the development of new genomic methods.

Regulation of Chromatin Architecture by Transcription Factor Binding

Transcription factors (TF) bind to chromatin and regulate the expression of genes. The pair Myc:Max binds to E-box regulatory DNA elements throughout the genome, controlling transcription of a large group of specific genes. We introduce an implicit modeling protocol for Myc:Max binding to mesoscale chromatin fibers to determine TF effect on chromatin architecture and shed light on its mechanism of gene regulation. We first bind Myc:Max to different chromatin locations and show how it can direct fiber folding and formation of microdomains, and how this depends on the linker DNA length. Second, by simulating increasing concentrations of Myc:Max binding to fibers that differ in the DNA linker length, linker histone density, and acetylation levels, we assess the interplay between Myc:Max and other chromatin internal parameters. Third, we study the mechanism of gene silencing by Myc:Max binding to the Eed gene loci. Overall, our results show how chromatin architecture can be regulated by TF binding. The position of TF binding dictates the formation of microdomains that appear visible only at the ensemble level. On the other hand, the presence of linker histone, acetylations, or different linker DNA lengths regulates the concentration-dependent effect of TF binding. Furthermore, we show how TF binding can repress gene expression by increasing fiber folding motifs that help compact and occlude the promoter region. Importantly, this effect can be reversed by increasing linker histone density. Overall, these results shed light on the epigenetic control of the genome dictated by TF binding.

DNA elements tether canonical Polycomb Repressive Complex 1 to human genes

Development of multicellular animals requires epigenetic repression by Polycomb group proteins. The latter assemble in multi-subunit complexes, of which two kinds, Polycomb Repressive Complex 1 (PRC1) and Polycomb Repressive Complex 2 (PRC2), act together to repress key developmental genes. How PRC1 and PRC2 recognize specific genes remains an open question. Here we report the identification of several hundreds of DNA elements that tether canonical PRC1 to human developmental genes. We use the term tether to describe a process leading to a prominent presence of canonical PRC1 at certain genomic sites, although the complex is unlikely to interact with DNA directly. Detailed analysis indicates that sequence features associated with PRC1 tethering differ from those that favour PRC2 binding. Throughout the genome, the two kinds of sequence features mix in different proportions to yield a gamut of DNA elements that range from those tethering predominantly PRC1 or PRC2 to ones capable of tethering both complexes. The emerging picture is similar to the paradigmatic targeting of Polycomb complexes by Polycomb Response Elements (PREs) of Drosophila but providing for greater plasticity.

Different elongation factors distinctly modulate RNA polymerase II transcription in Arabidopsis

Various transcript elongation factors (TEFs) including modulators of RNA polymerase II (RNAPII) activity and histone chaperones tune the efficiency of transcription in the chromatin context. TEFs are involved in establishing gene expression patterns during growth and development in Arabidopsis, while little is known about the genomic distribution of the TEFs and the way they facilitate transcription. We have mapped the genome-wide occupancy of the elongation factors SPT4-SPT5, PAF1C and FACT, relative to that of elongating RNAPII phosphorylated at residues S2/S5 within the carboxyterminal domain. The distribution of SPT4-SPT5 along transcribed regions closely resembles that of RNAPII-S2P, while the occupancy of FACT and PAF1C is rather related to that of RNAPII-S5P. Under transcriptionally challenging heat stress conditions, mutant plants lacking the corresponding TEFs are differentially impaired in transcript synthesis. Strikingly, in plants deficient in PAF1C, defects in transcription across intron/exon borders are observed that are cumulative along transcribed regions. Upstream of transcriptional start sites, the presence of FACT correlates with nucleosomal occupancy. Under stress conditions FACT is particularly required for transcriptional upregulation and to promote RNAPII transcription through +1 nucleosomes. Thus, Arabidopsis TEFs are differently distributed along transcribed regions, and are distinctly required during transcript elongation especially upon transcriptional reprogramming.

Changes in adenoviral chromatin organization precede early gene activation upon infection

Within the virion, adenovirus DNA associates with the virus-encoded, protamine-like structural protein pVII. Whether this association is organized, and how genome packaging changes during infection and subsequent transcriptional activation is currently unclear. Here, we combined RNA-seq, MNase-seq, ChIP-seq, and single genome imaging during early adenovirus infection to unveil the structure- and time-resolved dynamics of viral chromatin changes as well as their correlation with gene transcription. Our MNase mapping data indicates that the adenoviral genome is arranged in precisely positioned nucleoprotein particles with nucleosome-like characteristics, that we term adenosomes. We identified 238 adenosomes that are positioned by a DNA sequence code and protect about 60-70 bp of DNA. The incoming adenoviral genome is more accessible at early gene loci that undergo additional chromatin de-condensation upon infection. Histone H3.3 containing nucleosomes specifically replaces pVII at distinct genomic sites and at the transcription start sites of early genes. Acetylation of H3.3 is predominant at the transcription start sites and precedes transcriptional activation. Based on our results, we propose a central role for the viral pVII nucleoprotein architecture, which is required for the dynamic structural changes during early infection, including the regulation of nucleosome assembly prior to transcription initiation. Our study thus may aid the rational development of recombinant adenoviral vectors exhibiting sustained expression in gene therapy.

Altered cfDNA fragmentation profile in hypomethylated regions as diagnostic markers in breast cancer

BACKGROUND

Breast cancer, the most common malignancy in women worldwide, has been proven to have both altered plasma cell-free DNA (cfDNA) methylation and fragmentation profiles. Nevertheless, simultaneously detecting both of them for breast cancer diagnosis has never been reported. Moreover, although fragmentation pattern of cfDNA is determined by nuclease digestion of chromatin, structure of which may be affected by DNA methylation, whether cfDNA methylation and fragmentation are biologically related or not still remains unclear.

METHODS

Improved cfMeDIP-seq were utilized to characterize both cfDNA methylation and fragmentation profiles in 49 plasma samples from both healthy individuals and patients with breast cancer. The feasibility of using cfDNA fragmentation profile in hypo- and hypermethylated regions as diagnostic markers for breast cancer was evaluated.

RESULTS

Mean size of cfDNA fragments (100-220 bp) mapped to hypomethylated regions decreased more in patients with breast cancer (4.60 bp, 172.33 to 167.73 bp) than in healthy individuals (2.87 bp, 174.54 to 171.67 bp). Furthermore, proportion of short cfDNA fragments (100-150 bp) in hypomethylated regions when compared with it in hypermethylated regions was found to increase more in patients with breast cancer in two independent discovery cohort. The feasibility of using abnormality of short cfDNA fragments ratio in hypomethylated genomic regions for breast cancer diagnosis in validation cohort was evaluated. 7 out of 11 patients were detected as having breast cancer (63.6% sensitivity), whereas no healthy individuals were mis-detected (100% specificity).

CONCLUSION

We identified enriched short cfDNA fragments after 5mC-immunoprecipitation (IP) in patients with breast cancer, and demonstrated the enriched short cfDNA fragments might originated from hypomethylated genomic regions. Furthermore, we proved the feasibility of using differentially methylated regions (DMRs)-dependent cfDNA fragmentation profile for breast cancer diagnosis.

Bioinformatics for wet-lab scientists: practical application in sequencing analysis

BACKGROUND

Genomics data is available to the scientific community after publication of research projects and can be investigated for a multitude of research questions. However, in many cases deposited data is only assessed and used for the initial publication, resulting in valuable resources not being exploited to their full depth. MAIN: A likely reason for this is that many wetlab-based researchers are not formally trained to apply bioinformatic tools and may therefore assume that they lack the necessary experience to do so themselves. In this article, we present a series of freely available, predominantly web-based platforms and bioinformatic tools that can be combined in analysis pipelines to interrogate different types of next-generation sequencing data. Additionally to the presented exemplary route, we also list a number of alternative tools that can be combined in a mix-and-match fashion. We place special emphasis on tools that can be followed and used correctly without extensive prior knowledge in programming. Such analysis pipelines can be applied to existing data downloaded from the public domain or be compared to the results of own experiments.

CONCLUSION

Integrating transcription factor binding to chromatin (ChIP-seq) with transcriptional output (RNA-seq) and chromatin accessibility (ATAC-seq) can not only assist to form a deeper understanding of the molecular interactions underlying transcriptional regulation but will also help establishing new hypotheses and pre-testing them in silico.

Dynamic nucleosome remodeling mediated by YY1 underlies early mouse development

Nucleosome positioning can alter the accessibility of DNA-binding proteins to their cognate DNA elements, and thus its precise control is essential for cell identity and function. Mammalian preimplantation embryos undergo temporal changes in gene expression and cell potency, suggesting the involvement of dynamic epigenetic control during this developmental phase. However, the dynamics of nucleosome organization during early development are poorly understood. In this study, using a low-input MNase-seq method, we show that nucleosome positioning is globally obscure in zygotes but becomes well defined during subsequent development. Down-regulation of the chromatin assembly in embryonic stem cells can partially reverse nucleosome organization into a zygote-like pattern, suggesting a possible link between the chromatin assembly pathway and fuzzy nucleosomes in zygotes. We also reveal that YY1, a zinc finger-containing transcription factor expressed upon zygotic genome activation, regulates the de novo formation of well-positioned nucleosome arrays at the regulatory elements through identifying YY1-binding sites in eight-cell embryos. The YY1-binding regions acquire H3K27ac enrichment around the eight-cell and morula stages, and YY1 depletion impairs the morula-to-blastocyst transition. Thus, our study delineates the remodeling of nucleosome organization and its underlying mechanism during early mouse development.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

OCRFinder: a noise-tolerance machine learning method for accurately estimating open chromatin regions

Open chromatin regions are the genomic regions associated with basic cellular physiological activities, while chromatin accessibility is reported to affect gene expressions and functions. A basic computational problem is to efficiently estimate open chromatin regions, which could facilitate both genomic and epigenetic studies. Currently, ATAC-seq and cfDNA-seq (plasma cell-free DNA sequencing) are two popular strategies to detect OCRs. As cfDNA-seq can obtain more biomarkers in one round of sequencing, it is considered more effective and convenient. However, in processing cfDNA-seq data, due to the dynamically variable chromatin accessibility, it is quite difficult to obtain the training data with pure OCRs or non-OCRs, and leads to a noise problem for either feature-based approaches or learning-based approaches. In this paper, we propose a learning-based OCR estimation approach with a noise-tolerance design. The proposed approach, named OCRFinder, incorporates the ideas of ensemble learning framework and semi-supervised strategy to avoid potential overfitting of noisy labels, which are the false positives on OCRs and non-OCRs. Compared to different noise control strategies and state-of-the-art approaches, OCRFinder achieved higher accuracies and sensitivities in the experiments. In addition, OCRFinder also has an excellent performance in ATAC-seq or DNase-seq comparison experiments.

PGE2 alters chromatin through H2A.Z-variant enhancer nucleosome modification to promote hematopoietic stem cell fate

Prostaglandin E2 (PGE2) and 16,16-dimethyl-PGE2 (dmPGE2) are important regulators of hematopoietic stem and progenitor cell (HSPC) fate and offer potential to enhance stem cell therapies [C. Cutler et al. Blood 122, 3074-3081(2013); W. Goessling et al. Cell Stem Cell 8, 445-458 (2011); W. Goessling et al. Cell 136, 1136-1147 (2009)]. Here, we report that PGE2-induced changes in chromatin at enhancer regions through histone-variant H2A.Z permit acute inflammatory gene induction to promote HSPC fate. We found that dmPGE2-inducible enhancers retain MNase-accessible, H2A.Z-variant nucleosomes permissive of CREB transcription factor (TF) binding. CREB binding to enhancer nucleosomes following dmPGE2 stimulation is concomitant with deposition of histone acetyltransferases p300 and Tip60 on chromatin. Subsequent H2A.Z acetylation improves chromatin accessibility at stimuli-responsive enhancers. Our findings support a model where histone-variant nucleosomes retained within inducible enhancers facilitate TF binding. Histone-variant acetylation by TF-associated nucleosome remodelers creates the accessible nucleosome landscape required for immediate enhancer activation and gene induction. Our work provides a mechanism through which inflammatory mediators, such as dmPGE2, lead to acute transcriptional changes and modify HSPC behavior to improve stem cell transplantation.

Disruption of polyhomeotic polymerization decreases nucleosome occupancy and alters genome accessibility

Chromatin attains its three-dimensional (3D) conformation by establishing contacts between different noncontiguous regions. Sterile Alpha Motif (SAM)-mediated polymerization of the polyhomeotic (PH) protein regulates subnuclear clustering of Polycomb Repressive Complex 1 (PRC1) and chromatin topology. The mutations that perturb the ability of the PH to polymerize, disrupt long-range chromatin contacts, alter Hox gene expression, and lead to developmental defects. To understand the underlying mechanism, we combined the experiments and theory to investigate the effect of this SAM domain mutation on nucleosome occupancy and accessibility on a genome wide scale. Our data show that disruption of PH polymerization because of SAM domain mutation decreases nucleosome occupancy and alters accessibility. Polymer simulations investigating the interplay between distant chromatin contacts and nucleosome occupancy, both of which are regulated by PH polymerization, suggest that nucleosome density increases when contacts between different regions of chromatin are established. Taken together, it appears that SAM domain-mediated PH polymerization biomechanically regulates the organization of chromatin at multiple scales from nucleosomes to chromosomes and we suggest that higher order organization can have a top-down causation effect on nucleosome occupancy.

Rig1 receptor plays a critical role in cardiac reprogramming via YY1 signaling

We discovered that innate immunity plays an important role in the reprogramming of fibroblasts into cardiomyocytes. In this report, we define the role of a novel retinoic acid-inducible gene 1 Yin Yang 1 (Rig1:YY1) pathway. We found that fibroblast to cardiomyocyte reprogramming efficacy was enhanced by specific Rig1 activators. To understand the mechanism of action, we performed various transcriptomic, nucleosome occupancy, and epigenomic approaches. Analysis of the datasets indicated that Rig1 agonists had no effect on reprogramming-induced changes in nucleosome occupancy or loss of inhibitory epigenetic motifs. Instead, Rig1 agonists were found to modulate cardiac reprogramming by promoting the binding of YY1 specifically to cardiac genes. To conclude, these results show that the Rig1:YY1 pathway plays a critical role in fibroblast to cardiomyocyte reprogramming.

Toward cis-regulation in soybean: a 3D genome scope

In eukaryotic cells, 3D genome plays an important role in the regulation of gene spatiotemporal expression, which is essential for the biological and developmental processes in a life cycle. In the past decade, the development of high-throughput technologies greatly enhances our ability to map the 3D genome organization, identifies multiple 3D genome structures, and investigates the functional role of 3D genome organization in gene regulation, which facilitates our understandings of cis-regulatory landscape and biological development. Comparing with the comprehensive analyses of 3D genome in mammals and model plants, the progress in soybean is much less. Future development and application of tools to precisely manipulate 3D genome structure at different levels will significantly strengthen the functional genome study and molecular breeding in soybean. Here, we review the recent progresses in 3D genome study and discuss future directions, which may help to improve soybean 3D functional genome study and molecular breeding.

Factors and Methods for the Detection of Gene Expression Regulation

Gene-expression regulation involves multiple processes and a range of regulatory factors. In this review, we describe the key factors that regulate gene expression, including transcription factors (TFs), chromatin accessibility, histone modifications, DNA methylation, and RNA modifications. In addition, we also describe methods that can be used to detect these regulatory factors.

Characterizing crosstalk in epigenetic signaling to understand disease physiology

Epigenetics, the inheritance of genomic information independent of DNA sequence, controls the interpretation of extracellular and intracellular signals in cell homeostasis, proliferation and differentiation. On the chromatin level, signal transduction leads to changes in epigenetic marks, such as histone post-translational modifications (PTMs), DNA methylation and chromatin accessibility to regulate gene expression. Crosstalk between different epigenetic mechanisms, such as that between histone PTMs and DNA methylation, leads to an intricate network of chromatin-binding proteins where pre-existing epigenetic marks promote or inhibit the writing of new marks. The recent technical advances in mass spectrometry (MS) -based proteomic methods and in genome-wide DNA sequencing approaches have broadened our understanding of epigenetic networks greatly. However, further development and wider application of these methods is vital in developing treatments for disorders and pathologies that are driven by epigenetic dysregulation.

Genome-Wide Identification of Open Chromatin in Plants Using MH-Seq

Functional cis-regulatory elements (CREs) act as precise transcriptional switches for fine-tuning gene transcription. Identification of CREs is critical for understanding regulatory mechanisms of gene expression associated with various biological processes in eukaryotes. It is well known that CREs reside in open chromatin that exhibits hypersensitivity to enzyme cleavage and physical shearing. Currently, high-throughput methodologies, such as DNase-seq, ATAC-seq, and FAIRE-seq, have been widely applied in mapping open chromatin in various eukaryotic genomes. More recently, differential MNase (micrococcal nuclease) treatment has been successfully employed to map open chromatin in addition to profiling nucleosome landscape in both mammalian and plant species. We have developed a MNase hypersensitivity sequencing (MH-seq) technique in plants. The MH-seq procedure includes plant nuclei fixation and purification, differential treatments of purified nuclei with MNase, specific recovery of MNase-trimmed small DNA fragments within 20~100 bp in length, and MH-seq library construction followed by Illumina sequencing and data analysis. MH-seq has been successfully applied for global identification of open chromatin in both Arabidopsis thaliana and maize. It has been proven to be an attractive alternative for profiling open chromatin. Thus, MH-seq is expected to be valuable in probing chromatin accessibility on a genome-wide scale for other plants with sequenced genomes. Moreover, MHS data allow to implement footprinting assays to unveil binding sites of transcription factors.

Profiling Accessible Chromatin and Nucleosomes in the Mammalian Genome

Genomic DNA wraps around core histones to form nucleosomes, which provides steric constraints on how transcription factors (TFs) can interact with gene regulatory sequences. It is increasingly apparent that well-positioned, accessible nucleosomes are an inherent feature of active enhancers and can facilitate cooperative TF binding, referred to as nucleosome-mediated cooperativity. Thus, profiling chromatin and nucleosome properties (accessibility, positioning, and occupancy) on the genome is crucial to understand cell-type-specific gene regulation. Here we describe a simplified protocol to profile accessible nucleosomes in the mammalian genome using low-level and high-level micrococcal nuclease (MNase) digestion followed by genome-wide sequencing.

Comprehensive computational analysis of epigenetic descriptors affecting CRISPR-Cas9 off-target activity

BACKGROUND

A common issue in CRISPR-Cas9 genome editing is off-target activity, which prevents the widespread use of CRISPR-Cas9 in medical applications. Among other factors, primary chromatin structure and epigenetics may influence off-target activity.

METHODS

In this work, we utilize crisprSQL, an off-target database, to analyze the effect of 19 epigenetic descriptors on CRISPR-Cas9 off-target activity. Termed as 19 epigenetic features/scores, they consist of 6 experimental epigenetic and 13 computed nucleosome organization-related features. In terms of novel features, 15 of the epigenetic scores are newly considered. The 15 newly considered scores consist of 13 freshly computed nucleosome occupancy/positioning scores and 2 experimental features (MNase and DRIP). The other 4 existing scores are experimental features (CTCF, DNase I, H3K4me3, RRBS) commonly used in deep learning models for off-target activity prediction. For data curation, MNase was aggregated from existing experimental nucleosome occupancy data. Based on the sequence context information available in crisprSQL, we also computed nucleosome occupancy/positioning scores for off-target sites.

RESULTS

To investigate the relationship between the 19 epigenetic features and off-target activity, we first conducted Spearman and Pearson correlation analysis. Such analysis shows that some computed scores derived from training-based models and training-free algorithms outperform all experimental epigenetic features. Next, we evaluated the contribution of all epigenetic features in two successful machine/deep learning models which predict off-target activity. We found that some computed scores, unlike all 6 experimental features, significantly contribute to the predictions of both models. As a practical research contribution, we make the off-target dataset containing all 19 epigenetic features available to the research community.

CONCLUSIONS

Our comprehensive computational analysis helps the CRISPR-Cas9 community better understand the relationship between epigenetic features and CRISPR-Cas9 off-target activity.

Chromatin accessibility: methods, mechanisms, and biological insights

Access to DNA is a prerequisite to the execution of essential cellular processes that include transcription, replication, chromosomal segregation, and DNA repair. How the proteins that regulate these processes function in the context of chromatin and its dynamic architectures is an intensive field of study. Over the past decade, genome-wide assays and new imaging approaches have enabled a greater understanding of how access to the genome is regulated by nucleosomes and associated proteins. Additional mechanisms that may control DNA accessibility in vivo include chromatin compaction and phase separation - processes that are beginning to be understood. Here, we review the ongoing development of accessibility measurements, we summarize the different molecular and structural mechanisms that shape the accessibility landscape, and we detail the many important biological functions that are linked to chromatin accessibility.

DNA Sequence-Dependent Properties of Nucleosome Positioning in Regions of Distinct Chromatin States in Mouse Embryonic Stem Cells

Chromatin architecture is orchestrated, and plays crucial roles during the developmental process by regulating gene expression. In embryonic stem cells (ESCs), three types of chromatin states, including active, repressive and poised states, were previously identified and characterized with specific chromatin modification marks and different transcription activity, but it is largely unknown how nucleosomes are organized in these chromatin states. In this study, by using a DNA deformation energy model, we investigated the sequence-dependent nucleosome organization within the chromatin states in mouse ESCs. The results revealed that: (1) compared with poised genes, active genes are characterized with a higher level of nucleosome occupancy around their transcription start sites (TSS) and transcription termination sites (TTS), and both types of genes do not have a nucleosome-depleted region at their TTS, contrasting with the MNase-seq based result; (2) based on our previous DNA bending energy model, we developed an improved model capable of predicting both rotational positioning and nucleosome occupancy determined by a chemical mapping approach; (3) DNA bending-energy-based analyses demonstrated that the fragile nucleosomes positioned at both gene ends could be explained largely by enhanced rotational positioning signals encoded in DNA, but nucleosome phasing around the TSS of active genes was not determined by sequence preference; (4) the nucleosome occupancy landscape around the binding sites of some developmentally important transcription factors known to bind with different chromatin contexts, was also successfully predicted; (5) the difference of nucleosome occupancy around the TSS between CpG-rich and CpG-poor promoters was partly captured by our sequence-dependent model. Taken together, by developing an improved deformation-energy-based model, we revealed some sequence-dependent properties of the nucleosome arrangements in regions of distinct chromatin states in mouse ESCs.

Transcription factor-nucleosome dynamics from plasma cfDNA identifies ER-driven states in breast cancer

Genome-wide binding profiles of estrogen receptor (ER) and FOXA1 reflect cancer state in ER+ breast cancer. However, routine profiling of tumor transcription factor (TF) binding is impractical in the clinic. Here, we show that plasma cell-free DNA (cfDNA) contains high-resolution ER and FOXA1 tumor binding profiles for breast cancer. Enrichment of TF footprints in plasma reflects the binding strength of the TF in originating tissue. We defined pure in vivo tumor TF signatures in plasma using ER+ breast cancer xenografts, which can distinguish xenografts with distinct ER states. Furthermore, state-specific ER-binding signatures can partition human breast tumors into groups with significantly different ER expression and mortality. Last, TF footprints in human plasma samples can identify the presence of ER+ breast cancer. Thus, plasma TF footprints enable minimally invasive mapping of the regulatory landscape of breast cancer in humans and open vast possibilities for clinical applications across multiple tumor types.

Chromatin Dynamics in Digestive System Cancer: Commander and Regulator

Digestive system tumors have a poor prognosis due to complex anatomy, insidious onset, challenges in early diagnosis, and chemoresistance. Epidemiological statistics has verified that digestive system tumors rank first in tumor-related death. Although a great number of studies are devoted to the molecular biological mechanism, early diagnostic markers, and application of new targeted drugs in digestive system tumors, the therapeutic effect is still not satisfactory. Epigenomic alterations including histone modification and chromatin remodeling are present in human cancers and are now known to cooperate with genetic changes to drive the cancer phenotype. Chromatin is the carrier of genetic information and consists of DNA, histones, non-histone proteins, and a small amount of RNA. Chromatin and nucleosomes control the stability of the eukaryotic genome and regulate DNA processes such as transcription, replication, and repair. The dynamic structure of chromatin plays a key role in this regulatory function. Structural fluctuations expose internal DNA and thus provide access to the nuclear machinery. The dynamic changes are affected by various complexes and epigenetic modifications. Variation of chromatin dynamics produces early and superior regulation of the expression of related genes and downstream pathways, thereby controlling tumor development. Intervention at the chromatin level can change the process of cancer earlier and is a feasible option for future tumor diagnosis and treatment. In this review, we introduced chromatin dynamics including chromatin remodeling, histone modifications, and chromatin accessibility, and current research on chromatin regulation in digestive system tumors was also summarized.

A study of strong nucleosomes in the human genome

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, nucleosomes are highly dynamic and susceptible to experimental conditions, resulting in extreme variability across nucleosome maps, which complicates the generation of accurate nucleosome organization data. We mapped nucleosomes from different individuals using improved MNase-seq. The improvements included setting different digestion levels (low, medium, high) and naked DNA correction to remove the noise caused by experimental manipulation and comparing maps to obtain the accurate position and occupancy of strong nucleosomes (SNs) in the whole genome. In addition, the characteristics of SNs were further excavated. SNs were enriched in Alu elements and near the centromere of Chr12. SNs contain some specific sequences, and the GC content of SNs is different from that of dynamic nucleosomes. The findings suggest that nucleosome location in the genome and the DNA sequence may affect nucleosome stability.

•

Naked DNA correction improved the accuracy of nucleosome map in partial digestion

•

Level of MNase digestion has effects on nucleosome organization

•

A type of strong nucleosomes (SNs) exist across different nucleosome maps

•

Nucleosome stability may be related to its location and the DNA sequence

Nucleosome-Omics: A Perspective on the Epigenetic Code and 3D Genome Landscape

Genetic information is loaded on chromatin, which involves DNA sequence arrangement and the epigenetic landscape. The epigenetic information including DNA methylation, nucleosome positioning, histone modification, 3D chromatin conformation, and so on, has a crucial impact on gene transcriptional regulation. Out of them, nucleosomes, as basal chromatin structural units, play an important central role in epigenetic code. With the discovery of nucleosomes, various nucleosome-level technologies have been developed and applied, pushing epigenetics to a new climax. As the underlying methodology, next-generation sequencing technology has emerged and allowed scientists to understand the epigenetic landscape at a genome-wide level. Combining with NGS, nucleosome-omics (or nucleosomics) provides a fresh perspective on the epigenetic code and 3D genome landscape. Here, we summarized and discussed research progress in technology development and application of nucleosome-omics. We foresee the future directions of epigenetic development at the nucleosome level.

Mammalian PERIOD2 regulates H2A.Z incorporation in chromatin to orchestrate circadian negative feedback

Mammalian circadian oscillators are built on a feedback loop in which the activity of the transcription factor CLOCK-BMAL1 is repressed by the PER-CRY complex. Here, we show that murine Per^-/- fibroblasts display aberrant nucleosome occupancy around transcription start sites (TSSs) and at promoter-proximal and distal CTCF sites due to impaired histone H2A.Z deposition. Knocking out H2A.Z mimicked the Per null chromatin state and disrupted cellular rhythms. We found that endogenous mPER2 complexes retained CTCF as well as the specific H2A.Z-deposition chaperone YL1-a component of the ATP-dependent remodeler SRCAP and p400-TIP60 complex. While depleting YL1 or mutating chaperone-binding sites on H2A.Z lengthened the circadian period, H2A.Z deletion abrogated BMAL1 chromatin recruitment and promoted its proteasomal degradation. We propose that a PER2-mediated H2A.Z deposition pathway (1) compacts CLOCK-BMAL1 binding sites to establish negative feedback, (2) organizes circadian chromatin landscapes using CTCF and (3) bookmarks genomic loci for BMAL1 binding to impinge on the positive arm of the subsequent cycle.

Chromatin accessibility profiling by ATAC-seq

The assay for transposase-accessible chromatin using sequencing (ATAC-seq) provides a simple and scalable way to detect the unique chromatin landscape associated with a cell type and how it may be altered by perturbation or disease. ATAC-seq requires a relatively small number of input cells and does not require a priori knowledge of the epigenetic marks or transcription factors governing the dynamics of the system. Here we describe an updated and optimized protocol for ATAC-seq, called Omni-ATAC, that is applicable across a broad range of cell and tissue types. The ATAC-seq workflow has five main steps: sample preparation, transposition, library preparation, sequencing and data analysis. This protocol details the steps to generate and sequence ATAC-seq libraries, with recommendations for sample preparation and downstream bioinformatic analysis. ATAC-seq libraries for roughly 12 samples can be generated in 10 h by someone familiar with basic molecular biology, and downstream sequencing analysis can be implemented using benchmarked pipelines by someone with basic bioinformatics skills and with access to a high-performance computing environment.

Phosphorylation of the FACT histone chaperone subunit SPT16 affects chromatin at RNA polymerase II transcriptional start sites in Arabidopsis

The heterodimeric histone chaperone FACT, consisting of SSRP1 and SPT16, contributes to dynamic nucleosome rearrangements during various DNA-dependent processes including transcription. In search of post-translational modifications that may regulate the activity of FACT, SSRP1 and SPT16 were isolated from Arabidopsis cells and analysed by mass spectrometry. Four acetylated lysine residues could be mapped within the basic C-terminal region of SSRP1, while three phosphorylated serine/threonine residues were identified in the acidic C-terminal region of SPT16. Mutational analysis of the SSRP1 acetylation sites revealed only mild effects. However, phosphorylation of SPT16 that is catalysed by protein kinase CK2, modulates histone interactions. A non-phosphorylatable version of SPT16 displayed reduced histone binding and proved inactive in complementing the growth and developmental phenotypes of spt16 mutant plants. In plants expressing the non-phosphorylatable SPT16 version we detected at a subset of genes enrichment of histone H3 directly upstream of RNA polymerase II transcriptional start sites (TSSs) in a region that usually is nucleosome-depleted. This suggests that some genes require phosphorylation of the SPT16 acidic region for establishing the correct nucleosome occupancy at the TSS of active genes.

Superstructure Detection in Nucleosome Distribution Shows Common Pattern within a Chromosome and within the Genome

Nucleosome positioning plays an important role in crucial biological processes such as replication, transcription, and gene regulation. It has been widely used to predict the genome’s function and chromatin organisation. So far, the studies of patterns in nucleosome positioning have been limited to transcription start sites, CTCFs binding sites, and some promoter and loci regions. The genome-wide organisational pattern remains unknown. We have developed a theoretical model to coarse-grain nucleosome positioning data in order to obtain patterns in their distribution. Using hierarchical clustering on the auto-correlation function of this coarse-grained nucleosome positioning data, a genome-wide clustering is obtained for Candida albicans. The clustering shows the existence beyond hetero- and eu-chromatin inside the chromosomes. These non-trivial clusterings correspond to different nucleosome distributions and gene densities governing differential gene expression patterns. Moreover, these distribution patterns inside the chromosome appeared to be conserved throughout the genome and within species. The pipeline of the coarse grain nucleosome positioning sequence to identify underlying genomic organisation used in our study is novel, and the classifications obtained are unique and consistent.

Unveiling RCOR1 as a rheostat at transcriptionally permissive chromatin

RCOR1 is a known transcription repressor that recruits and positions LSD1 and HDAC1/2 on chromatin to erase histone methylation and acetylation. However, there is currently an incomplete understanding of RCOR1’s range of localization and function. Here, we probe RCOR1’s distribution on a genome-wide scale and unexpectedly find that RCOR1 is predominantly associated with transcriptionally active genes. Biochemical analysis reveals that RCOR1 associates with RNA Polymerase II (POL-II) during transcription and deacetylates its carboxy-terminal domain (CTD) at lysine 7. We provide evidence that this non-canonical RCOR1 activity is linked to dampening of POL-II productive elongation at actively transcribing genes. Thus, RCOR1 represses transcription in two ways—first, via a canonical mechanism by erasing transcriptionally permissive histone modifications through associating with HDACs and, second, via a non-canonical mechanism that deacetylates RNA POL-II’s CTD to inhibit productive elongation. We conclude that RCOR1 is a transcription rheostat.

The classical neuronal-gene corepressor RCOR1/CoREST is paradoxically enriched in transcriptionally active chromatin. Here the authors show RCOR1 is recruited during promoter-proximal pausing and negatively regulates the nascent-transcript synthesis. They also show that an RCOR1-LSD1- HDAC1 complex removes lysine acetylation from RNA polymerase II to repress transcription.

Hemogen /BRG1 cooperativity modulates promoter and enhancer activation during erythropoiesis

Hemogen, also known as EDAG, is a hematopoietic tissue-specific gene that regulates the proliferation and differentiation of hematopoietic cells. However, the mechanism underlying hemogen function in erythropoiesis is unknown. We found that depletion of hemogen in human CD34+ erythroid progenitor cells and HUDEP2 cells significantly reduced the expression of genes associated with heme and hemoglobin synthesis, supporting a positive role of hemogen in erythroid maturation. In human K562 cells, hemogen antagonized the occupancy of co-repressors NuRD complex and facilitated LDB1 complex-mediated chromatin looping. Hemogen recruited SWI/SNF complex ATPase BRG1 as a co-activator to regulate nucleosome accessibility and H3K27ac enrichment for promoter and enhancer activity. To ask if hemogen/BRG1 cooperativity is conserved in mammalian systems, we generated hemogen KO/KI mice and investigated hemogen/BRG1 function in murine erythropoiesis. Loss of hemogen in E12.5-E16.5 fetal liver cells impeded erythroid differentiation through reducing the production of mature erythroblasts. ChIP-seq in WT and hemogen KO animal revealed BRG1 is largely dependent on hemogen to regulate chromatin accessibility at erythroid gene promoters and enhancers. In summary, hemogen/BRG1 interaction in mammals is essential for fetal erythroid maturation and hemoglobin production through its active role in promoter and enhancer activity and chromatin organization.

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.

Integrative epigenomic profiling reveal AP-1 is a key regulator in intrahepatich cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is a malignant tumor with poor prognosis while its mechanisms of pathogenesis remain elusive. In this study, we performed systemic epigenomic and transcriptomic profiling via MNase-seq, ChIP-seq and RNA-seq in normal cholangiocyte and ICC cell lines. We showed that active histone modifications (H3K4me3, H3K4me1 and H3K27ac) were less enriched on cancer-related genes in ICC cell lines compared to control. The region of different histone modification patterns is enrichment in sites of AP-1 motif. Subsequent analysis showed that ICC had different nucleosome occupancy in differentially expressed genes compared to a normal cell line. Furthermore, we found that AP-1 plays a key role in ICC and regulates ICC-related genes through its AP-1 binding site. This study is the first report showing the global features of histone modification, transcript, and nucleosome profiles in ICC; we also show that the transcription factor AP-1 might be a key target gene in ICC.

Genome-wide cancer-specific chromatin accessibility patterns derived from archival processed xenograft tumors

Chromatin accessibility states that influence gene expression and other nuclear processes can be altered in disease. The constellation of transcription factors and chromatin regulatory complexes in cells results in characteristic patterns of chromatin accessibility. The study of these patterns in tissues has been limited because existing chromatin accessibility assays are ineffective for archival formalin-fixed, paraffin-embedded (FFPE) tissues. We have developed a method to efficiently extract intact chromatin from archival tissue via enhanced cavitation with a nanodroplet reagent consisting of a lipid shell with a liquid perfluorocarbon core. Inclusion of nanodroplets during the extraction of chromatin from FFPE tissues enhances the recovery of intact accessible and nucleosome-bound chromatin. We show that the addition of nanodroplets to the chromatin accessibility assay formaldehyde-assisted isolation of regulatory elements (FAIRE), does not affect the accessible chromatin signal. Applying the technique to FFPE human tumor xenografts, we identified tumor-relevant regions of accessible chromatin shared with those identified in primary tumors. Further, we deconvoluted non-tumor signal to identify cellular components of the tumor microenvironment. Incorporation of this method of enhanced cavitation into FAIRE offers the potential for extending chromatin accessibility to clinical diagnosis and personalized medicine, while also enabling the exploration of gene regulatory mechanisms in archival samples.