High-resolution mapping of transcription factor binding sites on native chromatin

Multiplexed chromatin immunoprecipitation sequencing for quantitative study of histone modifications and chromatin factors

ChIP-seq is a widely used technique for studying histone post-translational modifications and DNA-binding proteins. DNA fragments associated with a specific protein or histone modification epitope are captured by using antibodies, sequenced and mapped to a reference genome. Albeit versatile and popular, performing many parallel ChIP-seq experiments to compare different conditions, replicates and epitopes is laborious, is prone to experimental variation and does not allow quantitative comparisons unless adequate spike-in chromatin is included. We present a detailed protocol for performing and analyzing a multiplexed quantitative chromatin immunoprecipitation-sequencing experiment (MINUTE-ChIP), in which multiple samples are profiled against multiple epitopes in a single workflow. Multiplexing not only dramatically increases the throughput of ChIP-seq experiments (e.g., profiling 12 samples against multiple histone modifications or DNA-binding proteins in a single experiment), but also enables accurate quantitative comparisons. The protocol consists of four parts: sample preparation (i.e., lysis, chromatin fragmentation and barcoding of native or formaldehyde-fixed material), pooling and splitting of the barcoded chromatin into parallel immunoprecipitation reactions, preparation of next-generation sequencing libraries from input and immunoprecipitated DNA and data analysis using our dedicated analysis pipeline. This pipeline autonomously generates quantitatively scaled ChIP-seq tracks for downstream analysis and visualization, alongside necessary quality control indicators. The entire workflow requires basic knowledge in molecular biology and bioinformatics and can be completed in 1 week. MINUTE-ChIP empowers biologists to perform every ChIP-seq experiment with an appropriate number of replicates and control conditions, delivering more statistically robust, exquisitely quantitative and biologically meaningful results.

DMF-ChIP-seq for Highly Sensitive and Integrated Epigenomic Profiling of Low-Input Cells

Mapping genome-wide DNA-protein interactions (DPIs) provides insights into the epigenetic landscape of complex biological systems and elucidates the mechanisms of epigenetic regulation in biological progress. However, current technologies in DPI profiling still suffer from high cell demands, low detection sensitivity, and large reagent consumption. To address these problems, we developed DMF-ChIP-seq that builds on digital microfluidic (DMF) technology to profile genome-wide DPIs in a highly efficient, cost-effective, and user-friendly way. The entire workflow including cell pretreatment, antibody recognition, pA-Tn5 tagmentation, fragment enrichment, and PCR amplification is programmatically manipulated on a single chip. Leveraging closed submicroliter reaction volumes and a superhydrophobic interface, DMF-ChIP-seq presented higher sensitivity in peak enrichment than other current methods, with high accuracy (Pearson Correlation Coefficient (PCC) > 0.86) and high repeatability (PCC > 0.92). Furthermore, DMF-ChIP-seq was capable of processing the samples with as few as 8 cells while maintaining a high signal-to-noise ratio. By applying DMF-ChIP-seq, H3K27ac histone modification of early embryonic cells during differentiation was profiled for the investigation of epigenomic landscape dynamics. With the benefits of high efficiency and sensitivity in DPI analysis, the system provides great promise in studying epigenetic regulation during various biological processes.

H3K4me3 changes occur in cell wall genes during the development of Fagopyrum tataricum morphogenic and non-morphogenic calli

Epigenetic changes accompany the dynamic changes in the cell wall composition during the development of callus cells. H3K4me3 is responsible for active gene expression and reaction to environmental cues. Chromatin immunoprecipitation (ChIP) is a powerful technique for studying the interplay between epigenetic modifications and the DNA regions of interest. In combination with sequencing, it can provide the genome-wide enrichment of the specific epigenetic mark, providing vital information on its involvement in the plethora of cellular processes. Here, we describe the genome-wide distribution of H3K4me3 in morphogenic and non-morphogenic callus of Fagopyrum tataricum. Levels of H3K4me3 were higher around the transcription start site, in agreement with the role of this mark in transcriptional activation. The global levels of methylation were higher in the non-morphogenic callus, which indicated increased gene activation compared to the morphogenic callus. We also employed ChIP to analyse the changes in the enrichment of this epigenetic mark on the cell wall-related genes in both calli types during the course of the passage. Enrichment of H3K4me3 on cell wall genes was specific for callus type, suggesting that the role of this mark in cell-wall remodelling is complex and involved in many processes related to dedifferentiation and redifferentiation. This intricacy of the cell wall composition was supported by the immunohistochemical analysis of the cell wall epitopes' distribution of pectins and extensins. Together, these data give a novel insight into the involvement of H3K4me3 in the regeneration processes in F. tataricum in vitro callus tissue culture.

A continuum of zinc finger transcription factor retention on native chromatin underlies dynamic genome organization

Transcription factor (TF) residence on chromatin translates into quantitative transcriptional or structural outcomes on genome. Commonly used formaldehyde crosslinking fixes TF-DNA interactions cumulatively and compromises the measured occupancy level. Here we mapped the occupancy level of global or individual zinc finger TFs like CTCF and MAZ, in the form of highly resolved footprints, on native chromatin. By incorporating reinforcing perturbation conditions, we established S-score, a quantitative metric to proxy the continuum of CTCF or MAZ retention across different motifs on native chromatin. The native chromatin-retained CTCF sites harbor sequence features within CTCF motifs better explained by S-score than the metrics obtained from other crosslinking or native assays. CTCF retention on native chromatin correlates with local SUMOylation level, and anti-correlates with transcriptional activity. The S-score successfully delineates the otherwise-masked differential stability of chromatin structures mediated by CTCF, or by MAZ independent of CTCF. Overall, our study established a paradigm continuum of TF retention across binding sites on native chromatin, explaining the dynamic genome organization.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

MotifHub: Detection of trans-acting DNA motif group with probabilistic modeling algorithm

BACKGROUND

Trans-acting factors are of special importance in transcription regulation, which is a group of proteins that can directly or indirectly recognize or bind to the 8-12 bp core sequence of cis-acting elements and regulate the transcription efficiency of target genes. The progressive development in high-throughput chromatin capture technology (e.g., Hi-C) enables the identification of chromatin-interacting sequence groups where trans-acting DNA motif groups can be discovered. The problem difficulty lies in the combinatorial nature of DNA sequence pattern matching and its underlying sequence pattern search space.

METHOD

Here, we propose to develop MotifHub for trans-acting DNA motif group discovery on grouped sequences. Specifically, the main approach is to develop probabilistic modeling for accommodating the stochastic nature of DNA motif patterns.

RESULTS

Based on the modeling, we develop global sampling techniques based on EM and Gibbs sampling to address the global optimization challenge for model fitting with latent variables. The results reflect that our proposed approaches demonstrate promising performance with linear time complexities.

CONCLUSION

MotifHub is a novel algorithm considering the identification of both DNA co-binding motif groups and trans-acting TFs. Our study paves the way for identifying hub TFs of stem cell development (OCT4 and SOX2) and determining potential therapeutic targets of prostate cancer (FOXA1 and MYC). To ensure scientific reproducibility and long-term impact, its matrix-algebra-optimized source code is released at http://bioinfo.cs.cityu.edu.hk/MotifHub.

DisP-seq reveals the genome-wide functional organization of DNA-associated disordered proteins

Intrinsically disordered regions (IDRs) in DNA-associated proteins are known to influence gene regulation, but their distribution and cooperative functions in genome-wide regulatory programs remain poorly understood. Here we describe DisP-seq (disordered protein precipitation followed by DNA sequencing), an antibody-independent chemical precipitation assay that can simultaneously map endogenous DNA-associated disordered proteins genome-wide through a combination of biotinylated isoxazole precipitation and next-generation sequencing. DisP-seq profiles are composed of thousands of peaks that are associated with diverse chromatin states, are enriched for disordered transcription factors (TFs) and are often arranged in large lineage-specific clusters with high local concentrations of disordered proteins and different combinations of histone modifications linked to regulatory potential. We use DisP-seq to analyze cancer cells and reveal how disordered protein-associated islands enable IDR-dependent mechanisms that control the binding and function of disordered TFs, including oncogene-dependent sequestration of TFs through long-range interactions and the reactivation of differentiation pathways upon loss of oncogenic stimuli in Ewing sarcoma.

Native ChIP: Studying the Genome-Wide Distribution of Histone Modifications in Cells and Tissue

For the genome-wide mapping of histone modifications, chromatin immunoprecipitation (ChIP) followed by high-throughput sequencing remains the benchmark method. While crosslinked ChIP can be used for all kinds of targets, native ChIP is predominantly used for strong and direct DNA interactors like histones and their modifications. Here we describe a native ChIP protocol that can be used for cells and tissue material.

Epigenetic regulator RNF20 underlies temporal hierarchy of gene expression to regulate postnatal cardiomyocyte polarization

Differentiated cardiomyocytes (CMs) must undergo diverse morphological and functional changes during postnatal development. However, the mechanisms underlying initiation and coordination of these changes remain unclear. Here, we delineate an integrated, time-ordered transcriptional network that begins with expression of genes for cell-cell connections and leads to a sequence of structural, cell-cycle, functional, and metabolic transitions in mouse postnatal hearts. Depletion of histone H2B ubiquitin ligase RNF20 disrupts this gene network and impairs CM polarization. Subsequently, assay for transposase-accessible chromatin using sequencing (ATAC-seq) analysis confirmed that RNF20 contributes to chromatin accessibility in this context. As such, RNF20 is likely to facilitate binding of transcription factors at the promoters of genes involved in cell-cell connections and actin organization, which are crucial for CM polarization and functional integration. These results suggest that CM polarization is one of the earliest events during postnatal heart development and provide insights into how RNF20 regulates CM polarity and the postnatal gene program.

Functional genomics in Spiralia

Our understanding of the mechanisms that modulate gene expression in animals is strongly biased by studying a handful of model species that mainly belong to three groups: Insecta, Nematoda and Vertebrata. However, over half of the animal phyla belong to Spiralia, a morphologically and ecologically diverse animal clade with many species of economic and biomedical importance. Therefore, investigating genome regulation in this group is central to uncovering ancestral and derived features in genome functioning in animals, which can also be of significant societal impact. Here, we focus on five aspects of gene expression regulation to review our current knowledge of functional genomics in Spiralia. Although some fields, such as single-cell transcriptomics, are becoming more common, the study of chromatin accessibility, DNA methylation, histone post-translational modifications and genome architecture are still in their infancy. Recent efforts to generate chromosome-scale reference genome assemblies for greater species diversity and optimise state-of-the-art approaches for emerging spiralian research systems will address the existing knowledge gaps in functional genomics in this animal group.

High-resolution CTCF footprinting reveals impact of chromatin state on cohesin extrusion dynamics

DNA looping is vital for establishing many enhancer-promoter interactions. While CTCF is known to anchor many cohesin-mediated loops, the looped chromatin fiber appears to predominantly exist in a poorly characterized actively extruding state. To better characterize extruding chromatin loop structures, we used CTCF MNase HiChIP data to determine both CTCF binding at high resolution and 3D contact information. Here we present FactorFinder, a tool that identifies CTCF binding sites at near base-pair resolution. We leverage this substantial advance in resolution to determine that the fully extruded (CTCF-CTCF) state is rare genome-wide with locus-specific variation from ~1-10%. We further investigate the impact of chromatin state on loop extrusion dynamics, and find that active enhancers and RNA Pol II impede cohesin extrusion, facilitating an enrichment of enhancer-promoter contacts in the partially extruded loop state. We propose a model of topological regulation whereby the transient, partially extruded states play active roles in transcription.

Research progress of single-cell sequencing in tuberculosis

Tuberculosis is a major infectious disease caused by Mycobacterium tuberculosis infection. The pathogenesis and immune mechanism of tuberculosis are not clear, and it is urgent to find new drugs, diagnosis, and treatment targets. A useful tool in the quest to reveal the enigmas related to Mycobacterium tuberculosis infection and disease is the single-cell sequencing technique. By clarifying cell heterogeneity, identifying pathogenic cell groups, and finding key gene targets, the map at the single cell level enables people to better understand the cell diversity of complex organisms and the immune state of hosts during infection. Here, we briefly reviewed the development of single-cell sequencing, and emphasized the different applications and limitations of various technologies. Single-cell sequencing has been widely used in the study of the pathogenesis and immune response of tuberculosis. We review these works summarizing the most influential findings. Combined with the multi-molecular level and multi-dimensional analysis, we aim to deeply understand the blank and potential future development of the research on Mycobacterium tuberculosis infection using single-cell sequencing technology.

Simultaneous Mapping of DNA Binding and Nucleosome Positioning with SpLiT-ChEC

The organization of chromatin - including the positions of nucleosomes and the binding of other proteins to DNA - helps define transcriptional profiles in eukaryotic organisms. While techniques like ChIP-Seq and MNase-Seq can map protein-DNA and nucleosome localization separately, assays designed to simultaneously capture nucleosome positions and protein-DNA interactions can produce a detailed picture of the chromatin landscape. Most assays that monitor chromatin organization and protein binding rely on antibodies, which often exhibit nonspecific binding, and/or the addition of bulky adducts to the DNA-binding protein being studied, which can affect their expression and activity. Here, we describe SpyCatcher Linked Targeting of Chromatin Endogenous Cleavage (SpLiT-ChEC), where a 13-amino acid SpyTag peptide, appended to a protein of interest, serves as a highly-specific targeting moiety for in situ enzymatic digestion. The SpyTag/SpyCatcher system forms a covalent bond, linking the target protein and a co-expressed MNase-SpyCatcher fusion construct. SpyTagged proteins are expressed from endogenous loci, whereas MNase-SpyCatcher expression is induced immediately before harvesting cultures. MNase is activated with high concentrations of calcium, which primarily digests DNA near target protein binding sites. By sequencing the DNA fragments released by targeted MNase digestion, we found that this method recovers information on protein binding and proximal nucleosome positioning. SpLiT-ChEC provides precise temporal control that we anticipate can be used to monitor chromatin under various conditions and at distinct points in the cell cycle.

Genome-wide maps of rare and atypical UV photoproducts reveal distinct patterns of damage formation and mutagenesis in yeast chromatin

Ultraviolet (UV) light induces different classes of mutagenic photoproducts in DNA, namely cyclobutane pyrimidine dimers (CPDs), 6-4 photoproducts (6-4PPs), and atypical thymine-adenine photoproducts (TA-PPs). CPD formation is modulated by nucleosomes and transcription factors (TFs), which has important ramifications for Ultraviolet (UV) mutagenesis. How chromatin affects the formation of 6-4PPs and TA-PPs is unclear. Here, we use UV damage endonuclease-sequencing (UVDE-seq) to map these UV photoproducts across the yeast genome. Our results indicate that nucleosomes, the fundamental building block of chromatin, have opposing effects on photoproduct formation. Nucleosomes induce CPDs and 6-4PPs at outward rotational settings in nucleosomal DNA but suppress TA-PPs at these settings. Our data also indicate that DNA binding by different classes of yeast TFs causes lesion-specific hotspots of 6-4PPs or TA-PPs. For example, DNA binding by the TF Rap1 generally suppresses CPD and 6-4PP formation but induces a TA-PP hotspot. Finally, we show that 6-4PP formation is strongly induced at the binding sites of TATA-binding protein (TBP), which is correlated with higher mutation rates in UV-exposed yeast. These results indicate that the formation of 6-4PPs and TA-PPs is modulated by chromatin differently than CPDs and that this may have important implications for UV mutagenesis.

In situ tools for chromatin structural epigenomics

Technological progress over the past 15 years has fueled an explosion in genome-wide chromatin profiling tools that take advantage of low-cost short-read sequencing technologies to map particular chromatin features. Here, we survey the recent development of epigenomic tools that provide precise positions of chromatin proteins genome-wide in intact cells or nuclei. Some profiling tools are based on tethering Micrococcal Nuclease to chromatin proteins of interest in situ, whereas others similarly tether Tn5 transposase to integrate DNA sequencing adapters (tagmentation) and so eliminate the need for library preparation. These in situ cleavage and tagmentation tools have gained in popularity over the past few years, with many protocol enhancements and adaptations for single-cell and spatial chromatin profiling. The application of experimental and computational tools to address problems in gene regulation, eukaryotic development, and human disease are helping to define the emerging field of chromatin structural epigenomics.

The impact of single-cell genomics on the field of mycobacterial infection

Genome sequencing projects of humans and other organisms reinforced that the complexity of biological systems is largely attributed to the tight regulation of gene expression at the epigenome and RNA levels. As a consequence, plenty of technological developments arose to increase the sequencing resolution to the cell dimension creating the single-cell genomics research field. Single-cell RNA sequencing (scRNA-seq) is leading the advances in this topic and comprises a vast array of different methodologies. scRNA-seq and its variants are more and more used in life science and biomedical research since they provide unbiased transcriptomic sequencing of large populations of individual cells. These methods go beyond the previous “bulk” methodologies and sculpt the biological understanding of cellular heterogeneity and dynamic transcriptomic states of cellular populations in immunology, oncology, and developmental biology fields. Despite the large burden caused by mycobacterial infections, advances in this field obtained via single-cell genomics had been comparatively modest. Nonetheless, seminal research publications using single-cell transcriptomics to study host cells infected by mycobacteria have become recently available. Here, we review these works summarizing the most impactful findings and emphasizing the different and recent single-cell methodologies used, potential issues, and problems. In addition, we aim at providing insights into current research gaps and potential future developments related to the use of single-cell genomics to study mycobacterial infection.

Kinetic principles underlying pioneer function of GAGA transcription factor in live cells

How pioneer factors interface with chromatin to promote accessibility for transcription control is poorly understood in vivo. Here, we directly visualize chromatin association by the prototypical GAGA pioneer factor (GAF) in live Drosophila hemocytes. Single-particle tracking reveals that most GAF is chromatin bound, with a stable-binding fraction showing nucleosome-like confinement residing on chromatin for more than 2 min, far longer than the dynamic range of most transcription factors. These kinetic properties require the full complement of GAF's DNA-binding, multimerization and intrinsically disordered domains, and are autonomous from recruited chromatin remodelers NURF and PBAP, whose activities primarily benefit GAF's neighbors such as Heat Shock Factor. Evaluation of GAF kinetics together with its endogenous abundance indicates that, despite on-off dynamics, GAF constitutively and fully occupies major chromatin targets, thereby providing a temporal mechanism that sustains open chromatin for transcriptional responses to homeostatic, environmental and developmental signals.

nMOWChIP-seq: low-input genome-wide mapping of non-histone targets

Genome-wide profiling of interactions between genome and various functional proteins is critical for understanding regulatory processes involved in development and diseases. Conventional assays require a large number of cells and high-quality data on tissue samples are scarce. Here we optimized a low-input chromatin immunoprecipitation followed by sequencing (ChIP-seq) technology for profiling RNA polymerase II (Pol II), transcription factor (TF), and enzyme binding at the genome scale. The new approach produces high-quality binding profiles using 1,000-50,000 cells. We used the approach to examine the binding of Pol II and two TFs (EGR1 and MEF2C) in cerebellum and prefrontal cortex of mouse brain and found that their binding profiles are highly reflective of the functional differences between the two brain regions. Our analysis reveals the potential for linking genome-wide TF or Pol II profiles with neuroanatomical origins of brain cells.

High-resolution mapping demonstrates inhibition of DNA excision repair by transcription factors

DNA base damage arises frequently in living cells and needs to be removed by base excision repair (BER) to prevent mutagenesis and genome instability. Both the formation and repair of base damage occur in chromatin and are conceivably affected by DNA-binding proteins such as transcription factors (TFs). However, to what extent TF binding affects base damage distribution and BER in cells is unclear. Here, we used a genome-wide damage mapping method, N-methylpurine-sequencing (NMP-seq), and characterized alkylation damage distribution and BER at TF binding sites in yeast cells treated with the alkylating agent methyl methanesulfonate (MMS). Our data show that alkylation damage formation was mainly suppressed at the binding sites of yeast TFs ARS binding factor 1 (Abf1) and rDNA enhancer binding protein 1 (Reb1), but individual hotspots with elevated damage levels were also found. Additionally, Abf1 and Reb1 binding strongly inhibits BER in vivo and in vitro, causing slow repair both within the core motif and its adjacent DNA. Repair of ultraviolet (UV) damage by nucleotide excision repair (NER) was also inhibited by TF binding. Interestingly, TF binding inhibits a larger DNA region for NER relative to BER. The observed effects are caused by the TF–DNA interaction, because damage formation and BER can be restored by depletion of Abf1 or Reb1 protein from the nucleus. Thus, our data reveal that TF binding significantly modulates alkylation base damage formation and inhibits repair by the BER pathway. The interplay between base damage formation and BER may play an important role in affecting mutation frequency in gene regulatory regions.

Refactoring transcription factors for metabolic engineering

Due to the ability to regulate target metabolic pathways globally and dynamically, metabolic regulation systems composed of transcription factors have been widely used in metabolic engineering and synthetic biology. This review introduced the categories, action principles, prediction strategies, and related databases of transcription factors. Then, the application of global transcription machinery engineering technology and the transcription factor-based biosensors and quorum sensing systems are overviewed. In addition, strategies for optimizing the transcriptional regulatory tools' performance by refactoring transcription factors are summarized. Finally, the current limitations and prospects of constructing various regulatory tools based on transcription factors are discussed. This review will provide theoretical guidance for the rational design and construction of transcription factor-based metabolic regulation systems.

Three classes of epigenomic regulators converge to hyperactivate the essential maternal gene deadhead within a heterochromatin mini-domain

The formation of a diploid zygote is a highly complex cellular process that is entirely controlled by maternal gene products stored in the egg cytoplasm. This highly specialized transcriptional program is tightly controlled at the chromatin level in the female germline. As an extreme case in point, the massive and specific ovarian expression of the essential thioredoxin Deadhead (DHD) is critically regulated in Drosophila by the histone demethylase Lid and its partner, the histone deacetylase complex Sin3A/Rpd3, via yet unknown mechanisms. Here, we identified Snr1 and Mod(mdg4) as essential for dhd expression and investigated how these epigenomic effectors act with Lid and Sin3A to hyperactivate dhd. Using Cut&Run chromatin profiling with a dedicated data analysis procedure, we found that dhd is intriguingly embedded in an H3K27me3/H3K9me3-enriched mini-domain flanked by DNA regulatory elements, including a dhd promoter-proximal element essential for its expression. Surprisingly, Lid, Sin3a, Snr1 and Mod(mdg4) impact H3K27me3 and this regulatory element in distinct manners. However, we show that these effectors activate dhd independently of H3K27me3/H3K9me3, and that dhd remains silent in the absence of these marks. Together, our study demonstrates an atypical and critical role for chromatin regulators Lid, Sin3A, Snr1 and Mod(mdg4) to trigger tissue-specific hyperactivation within a unique heterochromatin mini-domain.

Multicellular development depends on a tight control of gene expression in each cell type. This relies on the coordinated activities of nuclear proteins that interact with DNA or its histone scaffold to promote or restrict gene transcription. For example, we previously showed that the histone modifying enzymes Lid and Sin3A/Rpd3 are required in Drosophila ovaries for the massive expression of deadhead (dhd), a gene encoding for a thioredoxin that is essential for fertility. In this paper, we have further identified two additional dhd regulators, Snr1 and Mod(mdg4) and dissected the mechanism behind hyperactivation of this gene. Using the epigenomic profiling method Cut&Run with a dedicated data analysis approach, we unexpectedly found that dhd is embedded in an unusual chromatin mini-domain featuring repressive histone modifications H3K27me3 and H3K9me3 and flanked by two regulatory elements. However, we further showed that Lid, Sin3A, Snr1 and Mod(mdg4) behave like obligatory activators of dhd independently of this mini-domain. Our study unveils how multiple broad-acting epigenomic effectors operate in non-canonical manners to ensure a critical and specialized gene activation event. These findings challenge our knowledge on these regulatory mechanisms and their roles in development and pathology.

CUT&RUN Profiling of the Budding Yeast Epigenome

Mapping the epigenome is key to describe the relationship between chromatin landscapes and the control of DNA-based cellular processes such as transcription. Cleavage under targets and release using nuclease (CUT&RUN) is an in situ chromatin profiling strategy in which controlled cleavage by antibody-targeted Micrococcal Nuclease solubilizes specific protein-DNA complexes for paired-end DNA sequencing. When applied to budding yeast, CUT&RUN profiling yields precise genome-wide maps of histone modifications, histone variants, transcription factors, and ATP-dependent chromatin remodelers, while avoiding cross-linking and solubilization issues associated with the most commonly used chromatin profiling technique Chromatin Immunoprecipitation (ChIP). Furthermore, targeted chromatin complexes cleanly released by CUT&RUN can be used as input for a subsequent native immunoprecipitation step (CUT&RUN.ChIP) to simultaneously map two epitopes in single molecules genome-wide. The intrinsically low background and high resolution of CUT&RUN and CUT&RUN.ChIP allows for identification of transient genomic features such as dynamic nucleosome-remodeling intermediates. Starting from cells, one can perform CUT&RUN or CUT&RUN.ChIP and obtain purified DNA for sequencing library preparation in 2 days.

CUT&RUN Identifies Centromeric DNA Regions of Rhodotorula toruloides IFO0880

Rhodotorula toruloides has been increasingly explored as a host for bioproduction of lipids, fatty acid derivatives, and terpenoids. Various genetic tools have been developed, but neither a centromere nor an autonomously replicating sequence (ARS), both necessary elements for stable episomal plasmid maintenance, have yet been reported. In this study, Cleavage Under Targets and Release Using Nuclease (CUT&RUN), a method used for genome-wide mapping DNA-protein interactions, was used to identify R. toruloides IFO0880 genomic regions associated with the centromeric histone H3 protein Cse4, a marker of centromeric DNA. Fifteen putative centromeres ranging from 8 to 19 kb in length were identified and analyzed, and four were tested for, but did not show, ARS activity. These centromeric sequences contained below average GC content, corresponded to transcriptional cold-spots, were primarily nonrepetitive, and shared some vestigial transposon-related sequences but otherwise did not show significant sequence conservation. Future efforts to identify an ARS in this yeast can utilize these centromeric DNA sequences to improve the stability of episomal plasmids derived from putative ARS elements.

Biophysics is reshaping our perception of the epigenome: from DNA-level to high-throughput studies

Epigenetic research holds great promise to advance our understanding of biomarkers and regulatory processes in health and disease. An increasing number of new approaches, ranging from molecular to biophysical analyses, enable identifying epigenetic changes on the level of a single gene or the whole epigenome. The aim of this review is to highlight how the field is shifting from completely molecular-biology-driven solutions to multidisciplinary strategies including more reliance on biophysical analysis tools. Biophysics not only offers technical advancements in imaging or structure analysis but also helps to explore regulatory interactions. New computational methods are also being developed to meet the demand of growing data volumes and their processing. Therefore, it is important to capture these new directions in epigenetics from a biophysical perspective and discuss current challenges as well as multiple applications of biophysical methods and tools. Specifically, we gradually introduce different biophysical research methods by first considering the DNA-level information and eventually higher-order chromatin structures. Moreover, we aim to highlight that the incorporation of bioinformatics, machine learning, and artificial intelligence into biophysical analysis allows gaining new insights into complex epigenetic processes. The gained understanding has already proven useful in translational and clinical research providing better patient stratification options or new therapeutic insights. Together, this offers a better readiness to transform bench-top experiments into industrial high-throughput applications with a possibility to employ developed methods in clinical practice and diagnostics.

Identification of chromatin states during zebrafish gastrulation using CUT&RUN and CUT&Tag

BACKGROUND

Cell fate decisions are governed by interactions between sequence-specific transcription factors and a dynamic chromatin landscape. Zebrafish offer a powerful system for probing the mechanisms that drive these cell fate choices, especially in the context of early embryogenesis. However, technical challenges associated with conventional methods for chromatin profiling have slowed progress toward understanding the exact relationships between chromatin changes, transcription factor binding, and cellular differentiation during zebrafish embryogenesis.

RESULTS

To overcome these challenges, we adapted the chromatin profiling methods Cleavage Under Targets and Release Using Nuclease (CUT&RUN) and CUT&Tag for use in zebrafish and applied these methods to generate high-resolution enrichment maps for H3K4me3, H3K27me3, H3K9me3, RNA polymerase II, and the histone variant H2A.Z using tissue isolated from whole, mid-gastrula stage embryos. Using this data, we identify a subset of genes that may be bivalently regulated during both zebrafish and mouse gastrulation, provide evidence for an evolving H2A.Z landscape during embryo development, and demonstrate the effectiveness of CUT&RUN for detecting H3K9me3 enrichment at repetitive sequences.

CONCLUSIONS

Our results demonstrate the power of combining CUT&RUN and CUT&Tag methods with the strengths of the zebrafish system to define emerging chromatin landscapes in the context of vertebrate embryogenesis.

The glucose-sensing transcription factor MLX balances metabolism and stress to suppress apoptosis and maintain spermatogenesis

Male germ cell (GC) production is a metabolically driven and apoptosis-prone process. Here, we show that the glucose-sensing transcription factor (TF) MAX-Like protein X (MLX) and its binding partner MondoA are both required for male fertility in the mouse, as well as survival of human tumor cells derived from the male germ line. Loss of Mlx results in altered metabolism as well as activation of multiple stress pathways and GC apoptosis in the testes. This is concomitant with dysregulation of the expression of male-specific GC transcripts and proteins. Our genomic and functional analyses identify loci directly bound by MLX involved in these processes, including metabolic targets, obligate components of male-specific GC development, and apoptotic effectors. These in vivo and in vitro studies implicate MLX and other members of the proximal MYC network, such as MNT, in regulation of metabolism and differentiation, as well as in suppression of intrinsic and extrinsic death signaling pathways in both spermatogenesis and male germ cell tumors (MGCTs).

Su(Hw) primes 66D and 7F Drosophila chorion genes loci for amplification through chromatin decondensation

Suppressor of Hairy wing [Su(Hw)] is an insulator protein that participates in regulating chromatin architecture and gene repression in Drosophila. In previous studies we have shown that Su(Hw) is also required for pre-replication complex (pre-RC) recruitment on Su(Hw)-bound sites (SBSs) in Drosophila S2 cells and pupa. Here, we describe the effect of Su(Hw) on developmentally regulated amplification of 66D and 7F Drosophila amplicons in follicle cells (DAFCs), widely used as models in replication studies. We show Su(Hw) binding co-localizes with all known DAFCs in Drosophila ovaries, whereas disruption of Su(Hw) binding to 66D and 7F DAFCs causes a two-fold decrease in the amplification of these loci. The complete loss of Su(Hw) binding to chromatin impairs pre-RC recruitment to all amplification regulatory regions of 66D and 7F loci at early oogenesis (prior to DAFCs amplification). These changes coincide with a considerable Su(Hw)-dependent condensation of chromatin at 66D and 7F loci. Although we observed the Brm, ISWI, Mi-2, and CHD1 chromatin remodelers at SBSs genome wide, their remodeler activity does not appear to be responsible for chromatin decondensation at the 66D and 7F amplification regulatory regions. We have discovered that, in addition to the CBP/Nejire and Chameau histone acetyltransferases, the Gcn5 acetyltransferase binds to 66D and 7F DAFCs at SBSs and this binding is dependent on Su(Hw). We propose that the main function of Su(Hw) in developmental amplification of 66D and 7F DAFCs is to establish a chromatin structure that is permissive to pre-RC recruitment.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

RoboCOP: jointly computing chromatin occupancy profiles for numerous factors from chromatin accessibility data

Chromatin is a tightly packaged structure of DNA and protein within the nucleus of a cell. The arrangement of different protein complexes along the DNA modulates and is modulated by gene expression. Measuring the binding locations and occupancy levels of different transcription factors (TFs) and nucleosomes is therefore crucial to understanding gene regulation. Antibody-based methods for assaying chromatin occupancy are capable of identifying the binding sites of specific DNA binding factors, but only one factor at a time. In contrast, epigenomic accessibility data like MNase-seq, DNase-seq, and ATAC-seq provide insight into the chromatin landscape of all factors bound along the genome, but with little insight into the identities of those factors. Here, we present RoboCOP, a multivariate state space model that integrates chromatin accessibility data with nucleotide sequence to jointly compute genome-wide probabilistic scores of nucleosome and TF occupancy, for hundreds of different factors. We apply RoboCOP to MNase-seq and ATAC-seq data to elucidate the protein-binding landscape of nucleosomes and 150 TFs across the yeast genome, and show that our model makes better predictions than existing methods. We also compute a chromatin occupancy profile of the yeast genome under cadmium stress, revealing chromatin dynamics associated with transcriptional regulation.

The zinc finger protein CLAMP promotes long-range chromatin interactions that mediate dosage compensation of the Drosophila male X-chromosome

Background

Drosophila dosage compensation is an important model system for defining how active chromatin domains are formed. The male-specific lethal dosage compensation complex (MSLc) increases transcript levels of genes along the length of the single male X-chromosome to equalize with that expressed from the two female X-chromosomes. The strongest binding sites for MSLc cluster together in three-dimensional space largely independent of MSLc because clustering occurs in both sexes. CLAMP, a non-sex specific, ubiquitous zinc finger protein, binds synergistically with MSLc to enrich the occupancy of both factors on the male X-chromosome.

Results

Here, we demonstrate that CLAMP promotes the observed three-dimensional clustering of MSLc binding sites. Moreover, the X-enriched CLAMP protein more strongly promotes longer-range three-dimensional interactions on the X-chromosome than autosomes. Genome-wide, CLAMP promotes three-dimensional interactions between active chromatin regions together with other insulator proteins.

Conclusion

Overall, we define how long-range interactions which are modulated by a locally enriched ubiquitous transcription factor promote hyper-activation of the X-chromosome to mediate dosage compensation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00399-3.

G-quadruplexes are transcription factor binding hubs in human chromatin

BACKGROUND

The binding of transcription factors (TF) to genomic targets is critical in the regulation of gene expression. Short, double-stranded DNA sequence motifs are routinely implicated in TF recruitment, but many questions remain on how binding site specificity is governed.

RESULTS

Herein, we reveal a previously unappreciated role for DNA secondary structures as key features for TF recruitment. In a systematic, genome-wide study, we discover that endogenous G-quadruplex secondary structures (G4s) are prevalent TF binding sites in human chromatin. Certain TFs bind G4s with affinities comparable to double-stranded DNA targets. We demonstrate that, in a chromatin context, this binding interaction is competed out with a small molecule. Notably, endogenous G4s are prominent binding sites for a large number of TFs, particularly at promoters of highly expressed genes.

CONCLUSIONS

Our results reveal a novel non-canonical mechanism for TF binding whereby G4s operate as common binding hubs for many different TFs to promote increased transcription.

A plug and play microfluidic platform for standardized sensitive low-input chromatin immunoprecipitation

Epigenetic profiling by chromatin immunoprecipitation followed by sequencing (ChIP-seq) has become a powerful tool for genome-wide identification of regulatory elements, for defining transcriptional regulatory networks, and for screening for biomarkers. However, the ChIP-seq protocol for low-input samples is laborious and time-consuming and suffers from experimental variation, resulting in poor reproducibility and low throughput. Although prototypic microfluidic ChIP-seq platforms have been developed, these are poorly transferable as they require sophisticated custom-made equipment and in-depth microfluidic and ChIP expertise, while lacking parallelization. To enable standardized, automated ChIP-seq profiling of low-input samples, we constructed microfluidic PDMS-based plates capable of performing 24 sensitive ChIP reactions within 30 min of hands-on time and 4.5 h of machine-running time. These disposable plates can be conveniently loaded into a widely available controller for pneumatics and thermocycling. In light of the plug and play (PnP) ChIP plates and workflow, we named our procedure PnP-ChIP-seq. We show high-quality ChIP-seq on hundreds to a few thousand of cells for all six post-translational histone modifications that are included in the International Human Epigenome Consortium set of reference epigenomes. PnP-ChIP-seq robustly detects epigenetic differences on promoters and enhancers between naive and more primed mouse embryonic stem cells (mESCs). Furthermore, we used our platform to generate epigenetic profiles of rare subpopulations of mESCs that resemble the two-cell stage of embryonic development. PnP-ChIP-seq allows nonexpert laboratories worldwide to conveniently run robust, standardized ChIP-seq, whereas its high throughput, consistency, and sensitivity pave the way toward large-scale profiling of precious sample types such as rare subpopulations of cells or biopsies.

Cooperative binding between distant transcription factors is a hallmark of active enhancers

Enhancers harbor binding motifs that recruit transcription factors (TFs) for gene activation. While cooperative binding of TFs at enhancers is known to be critical for transcriptional activation of a handful of developmental enhancers, the extent of TF cooperativity genome-wide is unknown. Here, we couple high-resolution nuclease footprinting with single-molecule methylation profiling to characterize TF cooperativity at active enhancers in the Drosophila genome. Enrichment of short micrococcal nuclease (MNase)-protected DNA segments indicates that the majority of enhancers harbor two or more TF-binding sites, and we uncover protected fragments that correspond to co-bound sites in thousands of enhancers. From the analysis of co-binding, we find that cooperativity dominates TF binding in vivo at the majority of active enhancers. Cooperativity is highest between sites spaced 50 bp apart, indicating that cooperativity occurs without apparent protein-protein interactions. Our findings suggest nucleosomes promoting cooperativity because co-binding may effectively clear nucleosomes and promote enhancer function.

GAGA factor: a multifunctional pioneering chromatin protein

The Drosophila GAGA factor (GAF) is a multifunctional protein implicated in nucleosome organization and remodeling, activation and repression of gene expression, long distance enhancer-promoter communication, higher order chromosome structure, and mitosis. This broad range of activities poses questions about how a single protein can perform so many seemingly different and unrelated functions. Current studies argue that GAF acts as a "pioneer" factor, generating nucleosome-free regions of chromatin for different classes of regulatory elements. The removal of nucleosomes from regulatory elements in turn enables other factors to bind to these elements and carry out their specialized functions. Consistent with this view, GAF associates with a collection of chromatin remodelers and also interacts with proteins implicated in different regulatory functions. In this review, we summarize the known activities of GAF and the functions of its protein partners.

High-Resolution Genome-Wide Occupancy in Candida spp. Using ChEC-seq

To persist in their dynamic human host environments, fungal pathogens must sense and adapt by modulating their gene expression to fulfill their cellular needs. Understanding transcriptional regulation on a global scale would uncover cellular processes linked to persistence and virulence mechanisms that could be targeted for antifungal therapeutics. Infections associated with the yeast Candida albicans, a highly prevalent fungal pathogen, and the multiresistant related species Candida auris are becoming a serious public health threat. To define the set of a gene regulated by a transcriptional regulator in C. albicans, chromatin immunoprecipitation (ChIP)-based techniques, including ChIP with microarray technology (ChIP-chip) or ChIP-DNA sequencing (ChIP-seq), have been widely used. Here, we describe a new set of PCR-based micrococcal nuclease (MNase)-tagging plasmids for C. albicans and other Candida spp. to determine the genome-wide location of any transcriptional regulator of interest using chromatin endogenous cleavage (ChEC) coupled to high-throughput sequencing (ChEC-seq). The ChEC procedure does not require protein-DNA cross-linking or sonication, thus avoiding artifacts related to epitope masking or the hyper-ChIPable euchromatic phenomenon. In a proof-of-concept application of ChEC-seq, we provided a high-resolution binding map of the SWI/SNF chromatin remodeling complex, a master regulator of fungal fitness in C. albicans, in addition to the transcription factor Nsi1 that is an ortholog of the DNA-binding protein Reb1 for which genome-wide occupancy was previously established in Saccharomyces cerevisiae The ChEC-seq procedure described here will allow a high-resolution genomic location definition which will enable a better understanding of transcriptional regulatory circuits that govern fungal fitness and drug resistance in these medically important fungi.IMPORTANCE Systemic fungal infections caused by Candida albicans and the "superbug" Candida auris are becoming a serious public health threat. The ability of these yeasts to cause disease is linked to their faculty to modulate the expression of genes that mediate their escape from the immune surveillance and their persistence in the different unfavorable niches within the host. Comprehensive knowledge on gene expression control of fungal fitness is consequently an interesting framework for the identification of essential infection processes that could be hindered by chemicals as potential therapeutics. Here, we expanded the use of ChEC-seq, a technique that was initially developed in the yeast model Saccharomyces cerevisiae to identify genes that are modulated by a transcriptional regulator, in pathogenic yeasts from the genus Candida This robust technique will allow a better characterization of key gene expression regulators and their contribution to virulence and antifungal resistance in these pathogenic yeasts.

Profiling chromatin regulatory landscape: insights into the development of ChIP-seq and ATAC-seq

Chromatin regulatory landscape plays a critical role in many disease processes and embryo development. Epigenome sequencing technologies such as chromatin immunoprecipitation sequencing (ChIP-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) have enabled us to dissect the pan-genomic regulatory landscape of cells and tissues in both time and space dimensions by detecting specific chromatin state and its corresponding transcription factors. Pioneered by the advancement of chromatin immunoprecipitation-chip (ChIP-chip) technology, abundant epigenome profiling technologies have become available such as ChIP-seq, DNase I hypersensitive site sequencing (DNase-seq), ATAC-seq and so on. The advent of single-cell sequencing has revolutionized the next-generation sequencing, applications in single-cell epigenetics are enriched rapidly. Epigenome sequencing technologies have evolved from low-throughput to high-throughput and from bulk sample to the single-cell scope, which unprecedentedly benefits scientists to interpret life from different angles. In this review, after briefly introducing the background knowledge of epigenome biology, we discuss the development of epigenome sequencing technologies, especially ChIP-seq & ATAC-seq and their current applications in scientific research. Finally, we provide insights into future applications and challenges.

Efficient low-cost chromatin profiling with CUT&Tag

We recently introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling strategy in which antibodies are bound to chromatin proteins in situ in permeabilized nuclei. These antibodies are then used to tether the cut-and-paste transposase Tn5. Activation of the transposase simultaneously cleaves DNA and adds adapters ('tagmentation') for paired-end DNA sequencing. Here, we introduce a streamlined CUT&Tag protocol that suppresses DNA accessibility artefacts to ensure high-fidelity mapping of the antibody-targeted protein and improves the signal-to-noise ratio over current chromatin profiling methods. Streamlined CUT&Tag can be performed in a single PCR tube, from cells to amplified libraries, providing low-cost genome-wide chromatin maps. By simplifying library preparation CUT&Tag requires less than a day at the bench, from live cells to sequencing-ready barcoded libraries. As a result of low background levels, barcoded and pooled CUT&Tag libraries can be sequenced for as little as $25 per sample. This enables routine genome-wide profiling of chromatin proteins and modifications and requires no special skills or equipment.

Genome-wide off-rates reveal how DNA binding dynamics shape transcription factor function

Protein-DNA interactions are dynamic, and these dynamics are an important aspect of chromatin-associated processes such as transcription or replication. Due to a lack of methods to study on- and off-rates across entire genomes, protein-DNA interaction dynamics have not been studied extensively. Here, we determine in vivo off-rates for the Saccharomyces cerevisiae chromatin organizing factor Abf1, at 191 sites simultaneously across the yeast genome. Average Abf1 residence times span a wide range, varying between 4.2 and 33 min. Sites with different off-rates are associated with different functional characteristics. This includes their transcriptional dependency on Abf1, nucleosome positioning and the size of the nucleosome-free region, as well as the ability to roadblock RNA polymerase II for termination. The results show how off-rates contribute to transcription factor function and that DIVORSEQ (Determining In Vivo Off-Rates by SEQuencing) is a meaningful way of investigating protein-DNA binding dynamics genome-wide.

Ecdysone-Induced 3D Chromatin Reorganization Involves Active Enhancers Bound by Pipsqueak and Polycomb

Evidence suggests that Polycomb (Pc) is present at chromatin loop anchors in Drosophila. Pc is recruited to DNA through interactions with the GAGA binding factors GAF and Pipsqueak (Psq). Using HiChIP in Drosophila cells, we find that the psq gene, which has diverse roles in development and tumorigenesis, encodes distinct isoforms with unanticipated roles in genome 3D architecture. The BR-C, ttk, and bab domain (BTB)-containing Psq isoform (Psq^L) colocalizes genome-wide with known architectural proteins. Conversely, Psq lacking the BTB domain (Psq^S) is consistently found at Pc loop anchors and at active enhancers, including those that respond to the hormone ecdysone. After stimulation by this hormone, chromatin 3D organization is altered to connect promoters and ecdysone-responsive enhancers bound by Psq^S. Our findings link Psq variants lacking the BTB domain to Pc-bound active enhancers, thus shedding light into their molecular function in chromatin changes underlying the response to hormone stimulus.

In vivo chromatin organization on native yeast telomeric regions is independent of a cis-telomere loopback conformation

Background

DNA packaging into chromatin regulates all DNA-related processes and at chromosomal ends could affect both essential functions of telomeres: protection against DNA damage response and telomere replication. Despite this primordial role of chromatin, little is known about chromatin organization, and in particular about nucleosome positioning on unmodified subtelomere–telomere junctions in Saccharomyces cerevisiae.

Results

By ChEC experiments and indirect end-labeling, we characterized nucleosome positioning as well as specialized protein–DNA associations on most subtelomere–telomere junctions present in budding yeast. The results show that there is a relatively large nucleosome-free region at chromosome ends. Despite the absence of sequence homologies between the two major classes of subtelomere–telomere junctions (i.e.: Y’-telomeres and X-telomeres), all analyzed subtelomere–telomere junctions show a terminal nucleosome-free region just distally from the known Rap1-covered telomeric repeats. Moreover, previous evidence suggested a telomeric chromatin fold-back structure onto subtelomeric areas that supposedly was implicated in chromosome end protection. The in vivo ChEC method used herein in conjunction with several proteins in a natural context revealed no evidence for such structures in bulk chromatin.

Conclusions

Our study allows a structural definition of the chromatin found at chromosome ends in budding yeast. This definition, derived with direct in vivo approaches, includes a terminal area that is free of nucleosomes, certain positioned nucleosomes and conserved DNA-bound protein complexes. This organization of subtelomeric and telomeric areas however does not include a telomeric cis-loopback conformation. We propose that the observations on such fold-back structures may report rare and/or transient associations and not stable or constitutive structures.

Resolving the 3D Landscape of Transcription-Linked Mammalian Chromatin Folding

Whereas folding of genomes at the large scale of epigenomic compartments and topologically associating domains (TADs) is now relatively well understood, how chromatin is folded at finer scales remains largely unexplored in mammals. Here, we overcome some limitations of conventional 3C-based methods by using high-resolution Micro-C to probe links between 3D genome organization and transcriptional regulation in mouse stem cells. Combinatorial binding of transcription factors, cofactors, and chromatin modifiers spatially segregates TAD regions into various finer-scale structures with distinct regulatory features including stripes, dots, and domains linking promoters-to-promoters (P-P) or enhancers-to-promoters (E-P) and bundle contacts between Polycomb regions. E-P stripes extending from the edge of domains predominantly link co-expressed loci, often in the absence of CTCF and cohesin occupancy. Acute inhibition of transcription disrupts these gene-related folding features without altering higher-order chromatin structures. Our study uncovers previously obscured finer-scale genome organization, establishing functional links between chromatin folding and gene regulation.

RoboCOP: Multivariate State Space Model Integrating Epigenomic Accessibility Data to Elucidate Genome-Wide Chromatin Occupancy

Chromatin is the tightly packaged structure of DNA and protein within the nucleus of a cell. The arrangement of different protein complexes along the DNA modulates and is modulated by gene expression. Measuring the binding locations and level of occupancy of different transcription factors (TFs) and nucleosomes is therefore crucial to understanding gene regulation. Antibody-based methods for assaying chromatin occupancy are capable of identifying the binding sites of specific DNA binding factors, but only one factor at a time. On the other hand, epigenomic accessibility data like ATAC-seq, DNase-seq, and MNase-seq provide insight into the chromatin landscape of all factors bound along the genome, but with minimal insight into the identities of those factors. Here, we present RoboCOP, a multivariate state space model that integrates chromatin information from epigenomic accessibility data with nucleotide sequence to compute genome-wide probabilistic scores of nucleosome and TF occupancy, for hundreds of different factors at once. RoboCOP can be applied to any epigenomic dataset that provides quantitative insight into chromatin accessibility in any organism, but here we apply it to MNase-seq data to elucidate the protein-binding landscape of nucleosomes and 150 TFs across the yeast genome. Using available protein-binding datasets from the literature, we show that our model more accurately predicts the binding of these factors genome-wide.

High Sensitivity Profiling of Chromatin Structure by MNase-SSP

A complete view of eukaryotic gene regulation requires that we accurately delineate how transcription factors (TFs) and nucleosomes are arranged along linear DNA in a sensitive, unbiased manner. Here we introduce MNase-SSP, a single-stranded sequencing library preparation method for nuclease-digested chromatin that enables simultaneous mapping of TF and nucleosome positions. As a proof of concept, we apply MNase-SSP toward the genome-wide, high-resolution mapping of nucleosome and TF occupancy in murine embryonic stem cells (mESCs). Compared with existing MNase-seq protocols, MNase-SSP markedly enriches for short DNA fragments, enabling detection of binding by subnucleosomal particles and TFs, in addition to nucleosomes. From these same data, we identify multiple, sequence-dependent binding modes of the architectural TF Ctcf and extend this analysis to the TF Nrsf/ Rest. Looking forward, we anticipate that single stranded protocol (SSP) adaptations of any protein-DNA interaction mapping technique (e.g., ChIP-exo and CUT&RUN) will enhance the information content of the resulting data.

Ramani et al. describe MNase-SSP, a single-stranded DNA sequencing library preparation method for profiling chromatin structure. MNase-SSP libraries harbor diminished sequence bias and capture shorter DNA fragments compared to traditional MNase-seq libraries. Applying MNase-SSP to murine embryonic stem cells enables simultaneous analysis of nucleosomal, subnucleosomal, and transcription factor-DNA interactions genome-wide.

Genome-wide MNase hypersensitivity assay unveils distinct classes of open chromatin associated with H3K27me3 and DNA methylation in Arabidopsis thaliana

BACKGROUND

Regulation of transcription depends on interactions between cis-regulatory elements (CREs) and regulatory proteins. Active CREs are imbedded in open chromatin that are accessible to nucleases. Several techniques, including DNase-seq, which is based on nuclease DNase I, and ATAC-seq, which is based on transposase Tn5, have been widely used to identify genomic regions associated with open chromatin. These techniques have played a key role in dissecting the regulatory networks in gene expression in both animal and plant species.

RESULTS

We develop a technique, named MNase hypersensitivity sequencing (MH-seq), to identify genomic regions associated with open chromatin in Arabidopsis thaliana. Genomic regions enriched with MH-seq reads are referred as MNase hypersensitive sites (MHSs). MHSs overlap with the majority (~ 90%) of the open chromatin identified previously by DNase-seq and ATAC-seq. Surprisingly, 22% MHSs are not covered by DNase-seq or ATAC-seq reads, which are referred to "specific MHSs" (sMHSs). sMHSs tend to be located away from promoters, and a substantial portion of sMHSs are derived from transposable elements. Most interestingly, genomic regions containing sMHSs are enriched with epigenetic marks, including H3K27me3 and DNA methylation. In addition, sMHSs show a number of distinct characteristics including association with transcriptional repressors. Thus, sMHSs span distinct classes of open chromatin that may not be accessible to DNase I or Tn5. We hypothesize that the small size of the MNase enzyme relative to DNase I or Tn5 allows its access to relatively more condensed chromatin domains.

CONCLUSION

MNase can be used to identify open chromatin regions that are not accessible to DNase I or Tn5. Thus, MH-seq provides an important tool to identify and catalog all classes of open chromatin in plants.

Epigenome Regulation by Dynamic Nucleosome Unwrapping

Gene regulation in eukaryotes requires the controlled access of sequence-specific transcription factors (TFs) to their sites in a chromatin landscape dominated by nucleosomes. Nucleosomes are refractory to TF binding, and often must be removed from regulatory regions. Recent genomic studies together with in vitro measurements suggest that the nucleosome barrier to TF binding is modulated by dynamic nucleosome unwrapping governed by ATP-dependent chromatin remodelers. Genome-wide occupancy and the regulation of subnucleosomal intermediates have gained recent attention with the application of high-resolution approaches for precision mapping of protein-DNA interactions. We summarize here recent findings on nucleosome substructures and TF binding dynamics, and highlight how unwrapped nucleosomal intermediates provide a novel signature of active chromatin.

A native chromatin immunoprecipitation (ChIP) protocol for studying histone modifications in strawberry fruits

BACKGROUND

Covalent modifications of histones and histone variants have great influence on chromatin structure, which is involved in the transcriptional regulation of gene expression. Chromatin immunoprecipitation (ChIP) is a powerful tool for studying in vivo DNA-histone interactions. Strawberry is a model for Rosaceae and non-climacteric fruits, in which histone modifications have been implicated to affect fruit development and ripening. However, a validated ChIP method has not been reported in strawberry, probably due to its high levels of polysaccharides which affect the quality of prepared chromatin and the efficiency of immunoprecipitation.

RESULTS

We describe a native chromatin immunoprecipitation (N-ChIP) protocol suitable for strawberry by optimizing the parameters for nuclei isolation, chromatin extraction, DNA fragmentation and validation analysis using quantitative real-time PCR (qRT-PCR). The qRT-PCR results show that both the active mark H3K36me3 and the silent mark H3K9me2 are efficiently immunoprecipitated for the enriched regions. Compared to X-ChIP (cross-linked chromatin followed by immunoprecipitation), our optimized N-ChIP procedure has a higher signal-to-noise ratio and a lower background for both the active and the silent histone modifications. Furthermore, high-throughput sequencing following N-ChIP demonstrates that nearly 90% of the enriched H3K9/K14ac peaks are overlapped between biological replicates, indicating its remarkable consistency and reproducibility.

CONCLUSIONS

An N-ChIP method suitable for the fleshy fruit tissues of woodland strawberry Fragaria vesca is described in this study. The efficiency and reproducibility of our optimized N-ChIP protocol are validated by both qRT-PCR and high-throughput sequencing. We conclude that N-ChIP is a more suitable method for strawberry fruit tissues relative to X-ChIP, which could be combined with high-throughput sequencing to investigate the impact of histone modifications in strawberry and potentially in other fruits with high content of polysaccharides.

Native internally calibrated chromatin immunoprecipitation for quantitative studies of histone post-translational modifications

Chromatin immunoprecipitation coupled to next-generation sequencing (ChIP-seq) has served as the central method for the study of histone modifications for the past decade. In ChIP-seq analyses, antibodies selectively capture nucleosomes bearing a modification of interest and the associated DNA is then mapped to the genome to determine the distribution of the mark. This approach has several important drawbacks: (i) ChIP interpretation necessitates the assumption of perfect antibody specificity, despite growing evidence that this is often not the case. (ii) Common methods for evaluating antibody specificity in other formats have little or no bearing on specificity within a ChIP experiment. (iii) Uncalibrated ChIP is reported as relative enrichment, which is biologically meaningless outside the experimental reference frame defined by a discrete immunoprecipitation (IP), thus preventing facile comparison across experimental conditions or modifications. (iv) Differential library amplification and loading onto next-generation sequencers, as well as computational normalization, can further compromise quantitative relationships that may exist between samples. Consequently, the researcher is presented with a series of potential pitfalls and is blind to nearly all of them. Here we provide a detailed protocol for internally calibrated ChIP (ICeChIP), a method we recently developed to resolve these problems by spike-in of defined nucleosomal standards within a ChIP procedure. This protocol is optimized for specificity and quantitative power, allowing for measurement of antibody specificity and absolute measurement of histone modification density (HMD) at genomic loci on a biologically meaningful scale enabling unambiguous comparisons. We provide guidance on optimal conditions for next-generation sequencing (NGS) and instructions for data analysis. This protocol takes between 17 and 18 h, excluding time for sequencing or bioinformatic analysis. The ICeChIP procedure enables accurate measurement of histone post-translational modifications (PTMs) genome-wide in mammalian cells as well as Drosophila melanogaster and Caenorhabditis elegans, indicating suitability for use in eukaryotic cells more broadly.

Genome-wide histone modification profiling of inner cell mass and trophectoderm of bovine blastocysts by RAT-ChIP

Chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq) has revolutionized our understanding of chromatin-related biological processes. The method, however, requires thousands of cells and has therefore limited applications in situations where cell numbers are limited. Here we describe a novel method called Restriction Assisted Tagmentation Chromatin Immunoprecipitation (RAT-ChIP) that enables global histone modification profiling from as few as 100 cells. The method is simple, cost-effective and takes a single day to complete. We demonstrate the sensitivity of the method by deriving the first genome-wide maps of histone H3K4me3 and H3K27me3 modifications of inner cell mass and trophectoderm of bovine blastocyst stage embryos.

Engineered Chromatin Remodeling Proteins for Precise Nucleosome Positioning

Regulation of chromatin structure is essential for controlling access of DNA to factors that require association with specific DNA sequences. Here we describe the development and validation of engineered chromatin remodeling proteins (E-ChRPs) for inducing programmable changes in nucleosome positioning by design. We demonstrate that E-ChRPs function both in vitro and in vivo to specifically reposition target nucleosomes and entire nucleosomal arrays. We show that induced, systematic positioning of nucleosomes over yeast Ume6 binding sites leads to Ume6 exclusion, hyperacetylation, and transcriptional induction at target genes. We also show that programmed global loss of nucleosome-free regions at Reb1 targets is generally inhibitory with mildly repressive transcriptional effects. E-ChRPs are compatible with multiple targeting modalities, including the SpyCatcher and dCas9 moieties, resulting in high versatility and enabling diverse future applications. Thus, engineered chromatin remodeling proteins represent a simple and robust means to probe and disrupt DNA-dependent processes in different chromatin contexts.