Chromatin landscape instructs precise transcription factor regulome during embryonic lineage specification

Embryos, originating from fertilized eggs, undergo continuous cell division and differentiation, accompanied by dramatic changes in transcription, translation, and metabolism. Chromatin regulators, including transcription factors (TFs), play indispensable roles in regulating these processes. Recently, the trophoblast regulator TFAP2C was identified as crucial in initiating early cell fate decisions. However, Tfap2c transcripts persist in both the inner cell mass and trophectoderm of blastocysts, prompting inquiry into Tfap2c's function in post-lineage establishment. In this study, we delineate the dynamics of TFAP2C during the mouse peri-implantation stage and elucidate its collaboration with the key lineage regulators CDX2 and NANOG. Importantly, we propose that de novo formation of H3K9me3 in the extraembryonic ectoderm during implantation antagonizes TFAP2C binding to crucial developmental genes, thereby maintaining its lineage identity. Together, these results highlight the plasticity of the chromatin environment in designating the genomic binding of highly adaptable lineage-specific TFs and regulating embryonic cell fates.

The CUT&RUN greenlist: genomic regions of consistent noise are effective normalizing factors for quantitative epigenome mapping

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a recent development for epigenome mapping, but its unique methodology can hamper proper quantitative analyses. As traditional normalization approaches have been shown to be inaccurate, we sought to determine endogenous normalization factors based on the human genome regions of constant nonspecific signal. This constancy was determined by applying Shannon's information entropy, and the set of normalizer regions, which we named the 'Greenlist', was extensively validated using publicly available datasets. We demonstrate here that the greenlist normalization outperforms the current top standards, and remains consistent across different experimental setups, cell lines and antibodies; the approach can even be applied to different species or to CUT&Tag. Requiring no additional experimental steps and no added cost, this approach can be universally applied to CUT&RUN experiments to greatly minimize the interference of technical variation over the biological epigenome changes of interest.

Cleavage Under Targets & Release Using Nuclease (CUT&RUN) of Histone Modifications in Epidermal Stem Cells of Adult Murine Skin

Cleavage Under Targets & Release Using Nuclease (CUT&RUN) has swiftly become the preferred procedure over the past few years for genomic mapping and detecting interactions between chromatin and its bound proteins. CUT&RUN is now being widely used for characterizing the epigenetic landscape in many cell types as it utilizes far less cell numbers when compared to Chromatin Immunoprecipitation-sequencing (ChIP-seq), thereby making it a powerful tool for researchers working with limited material. This protocol has been specifically optimized for detecting histone modifications in fluorescence-activated cell sorting (FACS)-isolated epidermal stem cells from adult mice.

Activity of the pleiotropic drug resistance transcription factors Pdr1p and Pdr3p is modulated by binding site flanking sequences

The transcription factors Pdr1p and Pdr3p regulate pleiotropic drug resistance (PDR) in Saccharomyces cerevisiae via the PDR responsive elements (PDREs) to modulate gene expression. However, the exact mechanisms underlying the differences in their regulons remain unclear. Employing genomic occupancy profiling (CUT&RUN), binding assays, and transcription studies, we characterized the differences in sequence specificity between transcription factors. Findings reveal distinct preferences for core PDRE sequences and the flanking sequences for both proteins. While flanking sequences moderately alter DNA binding affinity, they significantly impact Pdr1/3p transcriptional activity. Notably, both proteins demonstrated the ability to bind half sites, showing potential enhancement of transcription from adjacent PDREs. This insight sheds light on ways Pdr1/3p can differentially regulate PDR.

A histone deacetylase network regulates epigenetic reprogramming and viral silencing in HIV-infected cells

A long-lived latent reservoir of HIV-1-infected CD4 T cells persists with antiretroviral therapy and prevents cure. We report that the emergence of latently infected primary CD4 T cells requires the activity of histone deacetylase enzymes HDAC1/2 and HDAC3. Data from targeted HDAC molecules, an HDAC3-directed PROTAC, and CRISPR-Cas9 knockout experiments converge on a model where either HDAC1/2 or HDAC3 targeting can prevent latency, whereas all three enzymes must be targeted to achieve latency reversal. Furthermore, HDACi treatment targets features of memory T cells that are linked to proviral latency and persistence. Latency prevention is associated with increased H3K9ac at the proviral LTR promoter region and decreased H3K9me3, suggesting that this epigenetic switch is a key proviral silencing mechanism that depends on HDAC activity. These findings support further mechanistic work on latency initiation and eventual clinical studies of HDAC inhibitors to interfere with latency initiation.

Multiomic analysis implicates FOXO4 in genetic regulation of chick lens fiber cell differentiation

A classic model for identification of novel differentiation mechanisms and pathways is the eye lens that consists of a monolayer of quiescent epithelial cells that are the progenitors of a core of mature fully differentiated fiber cells. The differentiation of lens epithelial cells into fiber cells follows a coordinated program involving cell cycle exit, expression of key structural proteins and the hallmark elimination of organelles to achieve transparency. Although multiple mechanisms and pathways have been identified to play key roles in lens differentiation, the entirety of mechanisms governing lens differentiation remain to be discovered. A previous study established that specific chromatin accessibility changes were directly associated with the expression of essential lens fiber cell genes, suggesting that the activity of transcription factors needed for expression of these genes could be regulated through binding access to the identified chromatin regions. Sequence analysis of the identified chromatin accessible regions revealed enhanced representation of the binding sequence for the transcription factor FOXO4 suggesting a direct role for FOXO4 in expression of these genes. FOXO4 is known to regulate a variety of cellular processes including cellular response to metabolic and oxidative stress, cell cycle withdrawal, and homeostasis, suggesting a previously unidentified role for FOXO4 in the regulation of lens cell differentiation. To further evaluate the role of FOXO4 we employed a multiomics approach to analyze the relationship between genome-wide FOXO4 binding, the differentiation-specific expression of key genes, and chromatin accessibility. To better identify active promoters and enhancers we also examined histone modification through analysis of H3K27ac. Specific methods included CUT&RUN (FOXO4 binding and H3K27ac modification), RNA-seq (differentiation state specific gene expression), and ATAC-seq (chromatin accessibility). CUT&RUN identified 20,966 FOXO4 binding sites and 33,921 H3K27ac marked regions across the lens fiber cell genome. RNA-seq identified 956 genes with significantly greater expression levels in fiber cells compared to epithelial cells (log2FC > 0.7, q < 0.05) and 2548 genes with significantly lower expression levels (log2FC < -0.7, q < 0.05). Integrated analysis identified 1727 differentiation-state specific genes that were nearest neighbors to at least one FOXO4 binding site, including genes encoding lens gap junctions (GJA1, GJA3), lens structural proteins (BFSP1, CRYBB1, ASL1), and genes required for lens transparency (HSF4, NRCAM). Multiomics analysis comparing the identified FOXO4 binding sites in published ATAC-seq data revealed that chromatin accessibility was associated with FOXO4-dependent gene expression during lens differentiation. The results provide evidence for an important requirement for FOXO4 in the regulated expression of key genes required for lens differentiation and link epigenetic regulation of chromatin accessibility and H3K27ac histone modification with the function of FOXO4 in controlling lens gene expression during lens fiber cell differentiation.

Cell Type-Specific Whole-Genome Landscape of ΔFOSB Binding in the Nucleus Accumbens After Chronic Cocaine Exposure

BACKGROUND

The ability of neurons to respond to external stimuli involves adaptations of gene expression. Induction of the transcription factor ΔFOSB in the nucleus accumbens, a key brain reward region, is important for the development of drug addiction. However, a comprehensive map of ΔFOSB's gene targets has not yet been generated.

METHODS

We used CUT&RUN (cleavage under targets and release using nuclease) to map the genome-wide changes in ΔFOSB binding in the 2 main types of nucleus accumbens neurons-D1 or D2 medium spiny neurons-after chronic cocaine exposure. To annotate genomic regions of ΔFOSB binding sites, we also examined the distributions of several histone modifications. Resulting datasets were leveraged for multiple bioinformatic analyses.

RESULTS

The majority of ΔFOSB peaks occur outside promoter regions, including intergenic regions, and are surrounded by epigenetic marks indicative of active enhancers. BRG1, the core subunit of the SWI/SNF chromatin remodeling complex, overlaps with ΔFOSB peaks, a finding consistent with earlier studies of ΔFOSB's interacting proteins. Chronic cocaine use induces broad changes in ΔFOSB binding in both D1 and D2 nucleus accumbens medium spiny neurons of male and female mice. In addition, in silico analyses predict that ΔFOSB cooperatively regulates gene expression with homeobox and T-box transcription factors.

CONCLUSIONS

These novel findings uncover key elements of ΔFOSB's molecular mechanisms in transcriptional regulation at baseline and in response to chronic cocaine exposure. Further characterization of ΔFOSB's collaborative transcriptional and chromatin partners specifically in D1 and D2 medium spiny neurons will reveal a broader picture of the function of ΔFOSB and the molecular basis of drug addiction.

The CUT&RUN suspect list of problematic regions of the genome

BACKGROUND

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is an increasingly popular technique to map genome-wide binding profiles of histone modifications, transcription factors, and co-factors. The ENCODE project and others have compiled blacklists for ChIP-seq which have been widely adopted: these lists contain regions of high and unstructured signal, regardless of cell type or protein target, indicating that these are false positives. While CUT&RUN obtains similar results to ChIP-seq, its biochemistry and subsequent data analyses are different. We found that this results in a CUT&RUN-specific set of undesired high-signal regions.

RESULTS

We compile suspect lists based on CUT&RUN data for the human and mouse genomes, identifying regions consistently called as peaks in negative controls. Using published CUT&RUN data from our and other labs, we show that the CUT&RUN suspect regions can persist even when peak calling is performed with SEACR or MACS2 against a negative control and after ENCODE blacklist removal. Moreover, we experimentally validate the CUT&RUN suspect lists by performing reiterative negative control experiments in which no specific protein is targeted, showing that they capture more than 80% of the peaks identified.

CONCLUSIONS

We propose that removing these problematic regions can substantially improve peak calling in CUT&RUN experiments, resulting in more reliable datasets.

The time-resolved genomic impact of Wnt/β-catenin signaling

Wnt signaling orchestrates gene expression via its effector, β-catenin. However, it is unknown whether β-catenin binds its target genomic regions simultaneously and how this impacts chromatin dynamics to modulate cell behavior. Using a combination of time-resolved CUT&RUN against β-catenin, ATAC-seq, and perturbation assays in different cell types, we show that Wnt/β-catenin physical targets are tissue-specific, β-catenin "moves" on different loci over time, and its association to DNA accompanies changing chromatin accessibility landscapes that determine cell behavior. In particular, Wnt/β-catenin progressively shapes the chromatin of human embryonic stem cells (hESCs) as they undergo mesodermal differentiation, a behavior that we define as "plastic." In HEK293T cells, on the other hand, Wnt/β-catenin drives a transient chromatin opening, followed by re-establishment of the pre-stimulation state, a response that we define as "elastic." Future experiments shall assess whether other cell communication mechanisms, in addition to Wnt signaling, are ruled by time, cellular idiosyncrasies, and chromatin constraints. A record of this paper's transparent peer review process is included in the supplemental information.

Single-cell response to Wnt signaling activation reveals uncoupling of Wnt target gene expression

Wnt signaling drives nuclear translocation of β-catenin and its subsequent association with the DNA-bound TCF/LEF transcription factors, which dictate target gene specificity by recognizing Wnt responsive elements across the genome. β-Catenin target genes are therefore thought to be collectively activated upon Wnt pathway stimulation. However, this appears in contrast with the non-overlapping patterns of Wnt target gene expression in several contexts, including early mammalian embryogenesis. Here we followed Wnt target gene expression in human embryonic stem cells after Wnt pathway stimulation at a single-cell resolution. Cells changed gene expression program over time consistent with three key developmental events: i) loss of pluripotency, ii) induction of Wnt target genes, and iii) mesoderm specification. Contrary to our expectation, not all cells displayed equal amplitude of Wnt target gene activation; rather, they distributed in a continuum from strong to weak responders when ranked based on the expression of the target AXIN2. Moreover, high AXIN2 did not always correspond to elevated expression of other Wnt targets, which were activated in different proportions in individual cells. The uncoupling of Wnt target gene expression was also identified in single cell transcriptomics profiling of other Wnt-responding cell types, including HEK293T, murine developing forelimbs, and human colorectal cancer. Our finding underlines the necessity to identify additional mechanisms that explain the heterogeneity of the Wnt/β-catenin-mediated transcriptional outputs in single cells.

A single-cell multiomic analysis of kidney organoid differentiation

Kidney organoids differentiated from pluripotent stem cells are powerful models of kidney development and disease but are characterized by cell immaturity and off-target cell fates. Comparing the cell-specific gene regulatory landscape during organoid differentiation with human adult kidney can serve to benchmark progress in differentiation at the epigenome and transcriptome level for individual organoid cell types. Using single-cell multiome and histone modification analysis, we report more broadly open chromatin in organoid cell types compared to the human adult kidney. We infer enhancer dynamics by cis-coaccessibility analysis and validate an enhancer driving transcription of HNF1B by CRISPR interference both in cultured proximal tubule cells and also during organoid differentiation. Our approach provides an experimental framework to judge the cell-specific maturation state of human kidney organoids and shows that kidney organoids can be used to validate individual gene regulatory networks that regulate differentiation.

High-throughput methods for the analysis of transcription factors and chromatin modifications: Low input, single cell and spatial genomic technologies

Genome-wide analysis of transcription factors and epigenomic features is instrumental to shed light on DNA-templated regulatory processes such as transcription, cellular differentiation or to monitor cellular responses to environmental cues. Two decades of technological developments have led to a rich set of approaches progressively pushing the limits of epigenetic profiling towards single cells. More recently, disruptive technologies using innovative biochemistry came into play. Assays such as CUT&RUN, CUT&Tag and variations thereof show considerable potential to survey multiple TFs or histone modifications in parallel from a single experiment and in native conditions. These are in the path to become the dominant assays for genome-wide analysis of TFs and chromatin modifications in bulk, single-cell, and spatial genomic applications. The principles together with pros and cons are discussed.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

Disruption of the ATXN1-CIC complex reveals the role of additional nuclear ATXN1 interactors in spinocerebellar ataxia type 1

Spinocerebellar ataxia type 1 (SCA1) is a paradigmatic neurodegenerative disease in that it is caused by a mutation in a broadly expressed protein, ATXN1; however, only select populations of cells degenerate. The interaction of polyglutamine-expanded ATXN1 with the transcriptional repressor CIC drives cerebellar Purkinje cell pathogenesis; however, the importance of this interaction in other vulnerable cells remains unknown. Here, we mutated the 154Q knockin allele of Atxn1^154Q/2Q mice to prevent the ATXN1-CIC interaction globally. This normalized genome-wide CIC binding; however, it only partially corrected transcriptional and behavioral phenotypes, suggesting the involvement of additional factors in disease pathogenesis. Using unbiased proteomics, we identified three ATXN1-interacting transcription factors: RFX1, ZBTB5, and ZKSCAN1. We observed altered expression of RFX1 and ZKSCAN1 target genes in SCA1 mice and patient-derived iNeurons, highlighting their potential contributions to disease. Together, these data underscore the complexity of mechanisms driving cellular vulnerability in SCA1.

Efficient Genome-Wide Chromatin Profiling by CUT&RUN with Low Numbers of Muscle Stem Cells

Adult muscle stem cells (MuSCs), also called satellite cells, are situated under the basal lamina of myofibers in skeletal muscles. MuSCs are instrumental for postnatal muscle growth and regeneration of skeletal muscles. Under physiological conditions, the majority of MuSCs is actively maintained in a quiescent state but becomes rapidly activated during muscle regeneration, which is accompanied with massive changes in the epigenome. Moreover, aging, but also pathological conditions, such as in muscle dystrophy, results in profound changes of the epigenome, which can be monitored with different approaches. However, a better understanding of the role of chromatin dynamics in MuSCs and its function for skeletal muscle physiology and disease has been hampered by technical limitations, mostly due to the relatively low number of MuSCs but also due to the strongly condensed chromatin state of quiescent MuSCs. Traditional chromatin immunoprecipitation (ChIP) usually requires large amounts of cells and has several other shortcomings. Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a simple alternative to ChIP for chromatin profiling, providing higher efficiency and better resolution at lower costs. CUT&RUN maps genome-wide chromatin features, including genome-wide localization of transcription factor binding in small numbers of freshly isolated MuSCs, facilitating analysis of different subpopulations of MuSCs. Here we describe an optimized protocol to profile global chromatin in freshly isolated MuSCs using CUT&RUN.

CUT&Tag Using "Stress-Free" Con A-Conjugated Dynabeads®

Epigenome research has employed various methods to identify the genomic location of proteins of interest, such as transcription factors and histone modifications. CUT&Tag is a recently established method used in epigenome research to determine the genomic location of proteins of interest, such as transcription factors and histone modifications. In CUT&Tag method, cells are bound and hold on concanavalin A (con A)-coated magnetic beads, then a Protein-A Tn5 transposase fusion protein cuts the genome and inserts adapter sequences nearby the target protein. Here we describe the updated CUT&Tag procedure using "home-made" con A-conjugated magnetic beads. This method is free of poor suspendability and severe aggregation, hence providing improved sensitivity.

A new CUT&RUN low volume-urea (LoV-U) protocol optimized for transcriptional co-factors uncovers Wnt/β-catenin tissue-specific genomic targets

Upon WNT/β-catenin pathway activation, stabilized β-catenin travels to the nucleus where it associates with the TCF/LEF transcription factors, constitutively bound to genomic Wnt-responsive elements (WREs), to activate target gene transcription. Discovering the binding profile of β-catenin is therefore required to unambiguously assign direct targets of WNT signaling. Cleavage under targets and release using nuclease (CUT&RUN) has emerged as prime technique for mapping the binding profile of DNA-interacting proteins. Here, we present a modified version of CUT&RUN, named LoV-U (low volume and urea), that enables the robust and reproducible generation of β-catenin binding profiles, uncovering direct WNT/β-catenin target genes in human cells, as well as in cells isolated from developing mouse tissues. CUT&RUN-LoV-U outperforms original CUT&RUN when targeting co-factors that do not bind the DNA, can profile all classes of chromatin regulators and is well suited for simultaneous processing of several samples. We believe that the application of our protocol will allow the detection of the complex system of tissue-specific WNT/β-catenin target genes, together with other non-DNA-binding transcriptional regulators that act downstream of ontogenetically fundamental signaling cascades.

GoPeaks: histone modification peak calling for CUT&Tag

Genome-wide mapping of histone modifications is critical to understanding transcriptional regulation. CUT&Tag is a new method for profiling histone modifications, offering improved sensitivity and decreased cost compared with ChIP-seq. Here, we present GoPeaks, a peak calling method specifically designed for histone modification CUT&Tag data. We compare the performance of GoPeaks against commonly used peak calling algorithms to detect histone modifications that display a range of peak profiles and are frequently used in epigenetic studies. We find that GoPeaks robustly detects genome-wide histone modifications and, notably, identifies a substantial number of H3K27ac peaks with improved sensitivity compared to other standard algorithms.

CUT&RUN for Chromatin Profiling in Caenorhabditis elegans

Cleavage under targets and release using nuclease (CUT&RUN) is a recently developed chromatin profiling technique that uses a targeted micrococcal nuclease cleavage strategy to obtain high-resolution binding profiles of protein factors or to map histones with specific post-translational modifications. Due to its high sensitivity, CUT&RUN allows quality binding profiles to be obtained with only a fraction of the starting material and sequencing depth typically required for other chromatin profiling techniques such as chromatin immunoprecipitation. Although CUT&RUN has been widely adopted in multiple model systems, it has rarely been utilized in Caenorhabditis elegans, a model system of great importance to genomic research. Cell dissociation techniques, which are required for this approach, can be challenging in C. elegans due to the toughness of the worm's cuticle and the sensitivity of the cells themselves. Here, we describe a robust CUT&RUN protocol for use in C. elegans to determine the genome-wide localization of protein factors and specific histone marks. With a simple protocol utilizing live, uncrosslinked tissue as the starting material, performing CUT&RUN in worms has the potential to produce physiologically relevant data at a higher resolution than chromatin immunoprecipitation. This protocol involves a simple dissociation step to uniformly permeabilize worms while avoiding sample loss or cell damage, resulting in high-quality CUT&RUN profiles with as few as 100 worms and detectable signal with as few as 10 worms. This represents a significant advancement over chromatin immunoprecipitation, which typically uses thousands or hundreds of thousands of worms for a single experiment. The protocols presented here provide a detailed description of worm growth, sample preparation, CUT&RUN workflow, library preparation for high-throughput sequencing, and a basic overview of data analysis, making CUT&RUN simple and accessible for any worm lab. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Growth and synchronization of C. elegans Basic Protocol 2: Worm dissociation, sample preparation, and optimization Basic Protocol 3: CUT&RUN chromatin profiling Alternate Protocol: Improving CUT&RUN signal using a secondary antibody Basic Protocol 4: CUT&RUN library preparation for Illumina high-throughput sequencing Basic Protocol 5: Basic data analysis using Linux.

Identification of Transcription Factor Binding Sites by Cleavage Under Target and Release Using Nuclease in Zebrafish

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) has emerged as a chromatin profiling strategy that excels traditional methods. Although CUT&RUN has been widely utilized in mammalian cells, its use in the zebrafish is at its early stages. In this study, we have developed a protocol to successfully perform CUT&RUN to map transcription factor (TF) binding sites in embryonic, adult tissues, and FACS-sorted zebrafish cells. We also provide a detailed workflow for the identification of predicted TF binding sites that can be utilized in any animal species. Altogether, our strategy will expand this invaluable tool to the zebrafish community, improving the epigenetic resolution that can be achieved in this model organism.

SEAseq: a portable and cloud-based chromatin occupancy analysis suite

Background

Genome-wide protein-DNA binding is popularly assessed using specific antibody pulldown in Chromatin Immunoprecipitation Sequencing (ChIP-Seq) or Cleavage Under Targets and Release Using Nuclease (CUT&RUN) sequencing experiments. These technologies generate high-throughput sequencing data that necessitate the use of multiple sophisticated, computationally intensive genomic tools to make discoveries, but these genomic tools often have a high barrier to use because of computational resource constraints.

Results

We present a comprehensive, infrastructure-independent, computational pipeline called SEAseq, which leverages field-standard, open-source tools for processing and analyzing ChIP-Seq/CUT&RUN data. SEAseq performs extensive analyses from the raw output of the experiment, including alignment, peak calling, motif analysis, promoters and metagene coverage profiling, peak annotation distribution, clustered/stitched peaks (e.g. super-enhancer) identification, and multiple relevant quality assessment metrics, as well as automatic interfacing with data in GEO/SRA. SEAseq enables rapid and cost-effective resource for analysis of both new and publicly available datasets as demonstrated in our comparative case studies.

Conclusions

The easy-to-use and versatile design of SEAseq makes it a reliable and efficient resource for ensuring high quality analysis. Its cloud implementation enables a broad suite of analyses in environments with constrained computational resources. SEAseq is platform-independent and is aimed to be usable by everyone with or without programming skills. It is available on the cloud at https://platform.stjude.cloud/workflows/seaseq and can be locally installed from the repository at https://github.com/stjude/seaseq.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04588-z.

SIX1 reprograms myogenic transcription factors to maintain the rhabdomyosarcoma undifferentiated state

Rhabdomyosarcoma (RMS) is a pediatric muscle sarcoma characterized by expression of the myogenic lineage transcription factors (TFs) MYOD1 and MYOG. Despite high expression of these TFs, RMS cells fail to terminally differentiate, suggesting the presence of factors that alter their functions. Here, we demonstrate that the developmental TF SIX1 is highly expressed in RMS and critical for maintaining a muscle progenitor-like state. SIX1 loss induces differentiation of RMS cells into myotube-like cells and impedes tumor growth in vivo. We show that SIX1 maintains the RMS undifferentiated state by controlling enhancer activity and MYOD1 occupancy at loci more permissive to tumor growth over muscle differentiation. Finally, we demonstrate that a gene signature derived from SIX1 loss correlates with differentiation status and predicts RMS progression in human disease. Our findings demonstrate a master regulatory role of SIX1 in repression of RMS differentiation via genome-wide alterations in MYOD1 and MYOG-mediated transcription.

Profiling m6A RNA Modifications in Low Amounts of Plant Cells Using Maize Meiocytes

RNA modifications can influence gene expression via multiple aspects such as RNA stability and alternative splicing. The most prominent RNA modification is m6A (N6-methyladenosine). Its profiling from low starting amounts of <100 cells is challenging. We describe here a complete workflow from cell isolation to data analysis that is based on using an RNA CUT&RUN-supported m6A-RIP (RNA immunoprecipitation) procedure and a subsequent adaptor-tagging library synthesis. Male meiocytes isolated from maize anthers were used as a test system to establish the protocol.

Chromatin Profiling of Human Naïve Pluripotent Stem Cells

Chromatin immunoprecipitation combined with high-throughput sequencing (ChIP-sequencing) facilitates the genome-wide mapping of DNA sequences that are enriched for specific chromatin-binding proteins or histone post-translational modifications. More recently developed chromatin profiling methods called Cleavage Under Targets and Release Using Nuclease (CUT&RUN) and Cleavage Under Targets and Tagmentation (CUT&Tag) have adapted the ChIP-sequencing approach to produce similar data from a smaller amount of starting material, and while overcoming many of the conventional drawbacks of ChIP-sequencing. Here, we present detailed protocols for ChIP-seq, CUT&RUN, and CUT&Tag to profile genome-wide protein-DNA interactions in naïve human pluripotent stem cells.

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.

greenCUT&RUN: Efficient Genomic Profiling of GFP-Tagged Transcription Factors and Chromatin Regulators

Genome-wide mapping of transcription factors and chromatin regulators is important to distinguish their direct from indirect effects on gene transcription or chromatin function. Novel approaches for studying their genomic localization under native conditions, such us cleavage under target and release using nuclease (CUT&RUN), offer higher resolution and lower sequencing costs than classical chromatin immunoprecipitation (ChIP) assays, and require fewer cells but they still depend on the availability of high-quality antibodies. Here, we describe detailed and robust protocols for greenCUT&RUN, which is a generic CUT&RUN-based approach for mapping the genome-wide localization of green fluorescent protein (GFP)-tagged factors in intact mammalian cells. The greenCUT&RUN method makes use of a micrococcal nuclease (MNase) coupled to a high affinity nanobody against GFP, which exploits the accessibility of multiple surfaces of the GFP tag, thus eliminating issues of antibody variability and availability. We also provide efficient protocols for the expression and purification of two different GFP nanobodies, which recognize non-overlapping GFP epitopes and can be combined for a further gain in sensitivity and accuracy. Compared to traditional CUT&RUN, genomic localization by greenCUT&RUN reduces handling time and experimental variability. GreenCUT&RUN is a versatile, robust, and universal procedure for surveying the genome-wide localization of GFP-tagged versions of proteins that drive key transcriptional programs and regulate chromatin function. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Standard greenCUT&RUN for GFP-tagged proteins in mammalian cells Alternate Protocol: High-Ca⁺⁺ /low-salt greenCUT&RUN for GFP-tagged histone proteins in mammalian cells Support Protocol: Expression and purification of GFP nanobody-MNase fusion proteins for greenCUT&RUN.

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.

Establishment, maintenance, and recall of inflammatory memory

Known for nearly a century but through mechanisms that remain elusive, cells retain a memory of inflammation that equips them to react quickly and broadly to diverse secondary stimuli. Using murine epidermal stem cells as a model, we elucidate how cells establish, maintain, and recall inflammatory memory. Specifically, we landscape and functionally interrogate temporal, dynamic changes to chromatin accessibility, histone modifications, and transcription factor binding that occur during inflammation, post-resolution, and in memory recall following injury. We unearth an essential, unifying role for the general stress-responsive transcription factor FOS, which partners with JUN and cooperates with stimulus-specific STAT3 to establish memory; JUN then remains with other homeostatic factors on memory domains, facilitating rapid FOS re-recruitment and gene re-activation upon diverse secondary challenges. Extending our findings, we offer a comprehensive, potentially universal mechanism behind inflammatory memory and less discriminate recall phenomena with profound implications for tissue fitness in health and disease.

Integrative Cistromic and Transcriptomic Analyses Identify CREB Target Genes in Cystic Renal Epithelial Cells

BACKGROUND

Genome-wide mapping of transcription factor (TF) binding sites is essential to identify a TF's direct target genes in kidney development and diseases. However, due to the cellular complexity of the kidney and limited numbers of a given cell type, it has been challenging to determine the binding sites of a TF in vivo. cAMP response element-binding protein (CREB) is phosphorylated and hyperactive in autosomal dominant polycystic kidney disease (ADPKD). We focus on CREB as an example to profile genomic loci bound by a TF and to identify its target genes using low numbers of specific kidney cells.

METHODS

Cleavage under targets and release using nuclease (CUT&RUN) assays were performed with Dolichos biflorus agglutinin (DBA)-positive tubular epithelial cells from normal and ADPKD mouse kidneys. Pharmacologic inhibition of CREB with 666-15 and genetic inhibition with A-CREB were undertaken using ADPKD mouse models.

RESULTS

CUT&RUN to profile genome-wide distribution of phosphorylated CREB (p-CREB) indicated correlation of p-CREB binding with active histone modifications (H3K4me3 and H3K27ac) in cystic epithelial cells. Integrative analysis with CUT&RUN and RNA-sequencing revealed CREB direct targets, including genes involved in ribosome biogenesis and protein synthesis. Pharmacologic and genetic inhibition of CREB suppressed cyst growth in ADPKD mouse models.

CONCLUSIONS

CREB promotes cystogenesis by activating ribosome biogenesis genes. CUT&RUN, coupled with transcriptomic analysis, enables interrogation of TF binding and identification of direct TF targets from a low number of specific kidney cells.

ssvQC: an integrated CUT&RUN quality control workflow for histone modifications and transcription factors

OBJECTIVE

Among the different methods to profile the genome-wide patterns of transcription factor binding and histone modifications in cells and tissues, CUT&RUN has emerged as a more efficient approach that allows for a higher signal-to-noise ratio using fewer number of cells compared to ChIP-seq. The results from CUT&RUN and other related sequence enrichment assays requires comprehensive quality control (QC) and comparative analysis of data quality across replicates. While several computational tools currently exist for read mapping and analysis, a systematic reporting of data quality is lacking. Our aims were to (1) compare methods for using frozen versus fresh cells for CUT&RUN and (2) to develop an easy-to-use pipeline for assessing data quality.

RESULTS

We compared a workflow for CUT&RUN with fresh and frozen samples, and present an R package called ssvQC for quality control and comparison of data quality derived from CUT&RUN and other enrichment-based sequence data. Using ssvQC, we evaluate results from different CUT&RUN protocols for transcription factors and histone modifications from fresh and frozen tissue samples. Overall, this process facilitates evaluation of data quality across datasets and permits inspection of peak calling analysis, replicate analysis of different data types. The package ssvQC is readily available at https://github.com/FrietzeLabUVM/ssvQC .

Extensive NEUROG3 occupancy in the human pancreatic endocrine gene regulatory network

Objective

Mice lacking the bHLH transcription factor (TF) Neurog3 do not form pancreatic islet cells, including insulin-secreting beta cells, the absence of which leads to diabetes. In humans, homozygous mutations of NEUROG3 manifest with neonatal or childhood diabetes. Despite this critical role in islet cell development, the precise function of and downstream genetic programs regulated directly by NEUROG3 remain elusive. Therefore, we mapped genome-wide NEUROG3 occupancy in human induced pluripotent stem cell (hiPSC)–derived endocrine progenitors and determined NEUROG3 dependency of associated genes to uncover direct targets.

Methods

We generated a novel hiPSC line (NEUROG3-HA-P2A-Venus) where NEUROG3 is HA-tagged and fused to a self-cleaving fluorescent VENUS reporter. We used the CUT&RUN technique to map NEUROG3 occupancy and epigenetic marks in pancreatic endocrine progenitors (PEP) that were differentiated from this hiPSC line. We integrated NEUROG3 occupancy data with chromatin status and gene expression in PEPs as well as their NEUROG3-dependence. In addition, we investigated whether NEUROG3 binds type 2 diabetes mellitus (T2DM)–associated variants at the PEP stage.

Results

CUT&RUN revealed a total of 863 NEUROG3 binding sites assigned to 1263 unique genes. NEUROG3 occupancy was found at promoters as well as at distant cis-regulatory elements that frequently overlapped within PEP active enhancers. De novo motif analyses defined a NEUROG3 consensus binding motif and suggested potential co-regulation of NEUROG3 target genes by FOXA or RFX transcription factors. We found that 22% of the genes downregulated in NEUROG3^−/− PEPs, and 10% of genes enriched in NEUROG3-Venus positive endocrine cells were bound by NEUROG3 and thus likely to be directly regulated. NEUROG3 binds to 138 transcription factor genes, some with important roles in islet cell development or function, such as NEUROD1, PAX4, NKX2-2, SOX4, MLXIPL, LMX1B, RFX3, and NEUROG3 itself, and many others with unknown islet function. Unexpectedly, we uncovered that NEUROG3 targets genes critical for insulin secretion in beta cells (e.g., GCK, ABCC8/KCNJ11, CACNA1A, CHGA, SCG2, SLC30A8, and PCSK1). Thus, analysis of NEUROG3 occupancy suggests that the transient expression of NEUROG3 not only promotes islet destiny in uncommitted pancreatic progenitors, but could also initiate endocrine programs essential for beta cell function. Lastly, we identified eight T2DM risk SNPs within NEUROG3-bound regions.

Conclusion

Mapping NEUROG3 genome occupancy in PEPs uncovered unexpectedly broad, direct control of the endocrine genes, raising novel hypotheses on how this master regulator controls islet and beta cell differentiation.

•

NEUROG3 CUT&RUN analysis revealed 1263 target genes in human pancreatic endocrine progenitors (PEPs).

•

NEUROG3 binding sites overlap with active chromatin regions in PEPs.

•

1/5 of the genes downregulated in NEUROG3^−/− hESC-derived PEPs are bound by NEUROG3.

•

NEUROG3 targets islet-specific TFs and regulators of insulin secretion.

•

Several T2DM risk alleles lie within NEUROG3-bound regions.

Sox2 controls neural stem cell self-renewal through a Fos-centered gene regulatory network

The Sox2 transcription factor is necessary for the long-term self-renewal of neural stem cells (NSCs). Its mechanism of action is still poorly defined. To identify molecules regulated by Sox2, and acting in mouse NSC maintenance, we transduced, into Sox2-deleted NSC, genes whose expression is strongly downregulated following Sox2 loss (Fos, Jun, Egr2), individually or in combination. Fos alone rescued long-term proliferation, as shown by in vitro cell growth and clonal analysis. Furthermore, pharmacological inhibition by T-5224 of FOS/JUN AP1 complex binding to its targets decreased cell proliferation and expression of the putative target Suppressor of cytokine signaling 3 (Socs3). Additionally, Fos requirement for efficient long-term proliferation was demonstrated by the reduction of NSC clones capable of long-term expansion following CRISPR/Cas9-mediated Fos inactivation. Previous work showed that the Socs3 gene is strongly downregulated following Sox2 deletion, and its re-expression by lentiviral transduction rescues long-term NSC proliferation. Fos appears to be an upstream regulator of Socs3, possibly together with Jun and Egr2; indeed, Sox2 re-expression in Sox2-deleted NSC progressively activates both Fos and Socs3 expression; in turn, Fos transduction activates Socs3 expression. Based on available SOX2 ChIPseq and ChIA-PET data, we propose a model whereby Sox2 is a direct activator of both Socs3 and Fos, as well as possibly Jun and Egr2; furthermore, we provide direct evidence for FOS and JUN binding on Socs3 promoter, suggesting direct transcriptional regulation. These results provide the basis for developing a model of a network of interactions, regulating critical effectors of NSC proliferation and long-term maintenance.

FOS Rescues Neuronal Differentiation of Sox2-Deleted Neural Stem Cells by Genome-Wide Regulation of Common SOX2 and AP1(FOS-JUN) Target Genes

The transcription factor SOX2 is important for brain development and for neural stem cells (NSC) maintenance. Sox2-deleted (Sox2-del) NSC from neonatal mouse brain are lost after few passages in culture. Two highly expressed genes, Fos and Socs3, are strongly downregulated in Sox2-del NSC; we previously showed that Fos or Socs3 overexpression by lentiviral transduction fully rescues NSC's long-term maintenance in culture. Sox2-del NSC are severely defective in neuronal production when induced to differentiate. NSC rescued by Sox2 reintroduction correctly differentiate into neurons. Similarly, Fos transduction rescues normal or even increased numbers of immature neurons expressing beta-tubulinIII, but not more differentiated markers (MAP2). Additionally, many cells with both beta-tubulinIII and GFAP expression appear, indicating that FOS stimulates the initial differentiation of a "mixed" neuronal/glial progenitor. The unexpected rescue by FOS suggested that FOS, a SOX2 transcriptional target, might act on neuronal genes, together with SOX2. CUT&RUN analysis to detect genome-wide binding of SOX2, FOS, and JUN (the AP1 complex) revealed that a high proportion of genes expressed in NSC are bound by both SOX2 and AP1. Downregulated genes in Sox2-del NSC are highly enriched in genes that are also expressed in neurons, and a high proportion of the "neuronal" genes are bound by both SOX2 and AP1.

A functional map of genomic HIF1α-DNA complexes in the eye lens revealed through multiomics analysis

BACKGROUND

During eye lens development the embryonic vasculature regresses leaving the lens without a direct oxygen source. Both embryonically and throughout adult life, the lens contains a decreasing oxygen gradient from the surface to the core that parallels the natural differentiation of immature surface epithelial cells into mature core transparent fiber cells. These properties of the lens suggest a potential role for hypoxia and the master regulator of the hypoxic response, hypoxia-inducible transcription factor 1 (HIF1), in the regulation of genes required for lens fiber cell differentiation, structure and transparency. Here, we employed a multiomics approach combining CUT&RUN, RNA-seq and ATACseq analysis to establish the genomic complement of lens HIF1α binding sites, genes activated or repressed by HIF1α and the chromatin states of HIF1α-regulated genes.

RESULTS

CUT&RUN analysis revealed 8375 HIF1α-DNA binding complexes in the chick lens genome. One thousand one hundred ninety HIF1α-DNA binding complexes were significantly clustered within chromatin accessible regions (χ2 test p < 1 × 10- 55) identified by ATACseq. Formation of the identified HIF1α-DNA complexes paralleled the activation or repression of 526 genes, 116 of which contained HIF1α binding sites within 10kB of the transcription start sites. Some of the identified HIF1α genes have previously established lens functions while others have novel functions never before examined in the lens. GO and pathway analysis of these genes implicate HIF1α in the control of a wide-variety of cellular pathways potentially critical for lens fiber cell formation, structure and function including glycolysis, cell cycle regulation, chromatin remodeling, Notch and Wnt signaling, differentiation, development, and transparency.

CONCLUSIONS

These data establish the first functional map of genomic HIF1α-DNA complexes in the eye lens. They identify HIF1α as an important regulator of a wide-variety of genes previously shown to be critical for lens formation and function and they reveal a requirement for HIF1α in the regulation of a wide-variety of genes not yet examined for lens function. They support a requirement for HIF1α in lens fiber cell formation, structure and function and they provide a basis for understanding the potential roles and requirements for HIF1α in the development, structure and function of more complex tissues.

A unified atlas of CD8 T cell dysfunctional states in cancer and infection

CD8 T cells play an essential role in defense against viral and bacterial infections and in tumor immunity. Deciphering T cell loss of functionality is complicated by the conspicuous heterogeneity of CD8 T cell states described across experimental and clinical settings. By carrying out a unified analysis of over 300 assay for transposase-accessible chromatin sequencing (ATAC-seq) and RNA sequencing (RNA-seq) experiments from 12 studies of CD8 T cells in cancer and infection, we defined a shared differentiation trajectory toward dysfunction and its underlying transcriptional drivers and revealed a universal early bifurcation of functional and dysfunctional T cell states across models. Experimental dissection of acute and chronic viral infection using single-cell ATAC (scATAC)-seq and allele-specific single-cell RNA (scRNA)-seq identified state-specific drivers and captured the emergence of similar TCF1+ progenitor-like populations at an early branch point, at which functional and dysfunctional T cells diverge. Our atlas of CD8 T cell states will facilitate mechanistic studies of T cell immunity and translational efforts.

Binding of the Treslin-MTBP Complex to Specific Regions of the Human Genome Promotes the Initiation of DNA Replication

The processes that control where higher eukaryotic cells initiate DNA replication throughout the genome are not understood clearly. In metazoans, the Treslin-MTBP complex mediates critical final steps in formation of the activated replicative helicase prior to initiation of replication. Here, we map the genome-wide distribution of the MTBP subunit of this complex in human cells. Our results indicate that MTBP binds to at least 30,000 sites in the genome. A majority of these sites reside in regions of open chromatin that contain transcriptional-regulatory elements (e.g., promoters, enhancers, and super-enhancers), which are known to be preferred areas for initiation of replication. Furthermore, many binding sites encompass two genomic features: a nucleosome-free DNA sequence (e.g., G-quadruplex DNA or AP-1 motif) and a nucleosome bearing histone marks characteristic of open chromatin, such as H3K4me2. Taken together, these findings indicate that Treslin-MTBP associates coordinately with multiple genomic signals to promote initiation of replication.

Kumagai and Dunphy show that Treslin-MTBP, activator of the replicative helicase, binds to at least 30,000 sites in the human genome. Many sites contain a nucleosome with active chromatin marks and nucleosome-free DNA (G-quadruplex or AP-1 site). Thus, Treslin-MTBP associates with multiple genomic elements to promote initiation of DNA replication.

Zinc Finger Protein SALL4 Functions through an AT-Rich Motif to Regulate Gene Expression

The zinc finger transcription factor SALL4 is highly expressed in embryonic stem cells, downregulated in most adult tissues, but reactivated in many aggressive cancers. This unique expression pattern makes SALL4 an attractive therapeutic target. However, whether SALL4 binds DNA directly to regulate gene expression is unclear, and many of its targets in cancer cells remain elusive. Here, through an unbiased screen of protein binding microarray (PBM) and cleavage under targets and release using nuclease (CUT&RUN) experiments, we identify and validate the DNA binding domain of SALL4 and its consensus binding sequence. Combined with RNA sequencing (RNA-seq) analyses after SALL4 knockdown, we discover hundreds of new SALL4 target genes that it directly regulates in aggressive liver cancer cells, including genes encoding a family of histone 3 lysine 9-specific demethylases (KDMs). Taken together, these results elucidate the mechanism of SALL4 DNA binding and reveal pathways and molecules to target in SALL4-dependent tumors.

In this paper, Kong et al. elucidate the DNA binding mechanisms of the transcription factor SALL4 and an epigenetic pathway that it regulates. Due to its important role in driving aggressive cancers, better understanding of SALL4 function will lead to strategies to target this protein in cancer.

Conserved Gsx2/Ind homeodomain monomer versus homodimer DNA binding defines regulatory outcomes in flies and mice

How homeodomain proteins gain sufficient specificity to control different cell fates has been a long-standing problem in developmental biology. The conserved Gsx homeodomain proteins regulate specific aspects of neural development in animals from flies to mammals, and yet they belong to a large transcription factor family that bind nearly identical DNA sequences in vitro. Here, we show that the mouse and fly Gsx factors unexpectedly gain DNA binding specificity by forming cooperative homodimers on precisely spaced and oriented DNA sites. High-resolution genomic binding assays revealed that Gsx2 binds both monomer and homodimer sites in the developing mouse ventral telencephalon. Importantly, reporter assays showed that Gsx2 mediates opposing outcomes in a DNA binding site-dependent manner: Monomer Gsx2 binding represses transcription, whereas homodimer binding stimulates gene expression. In Drosophila, the Gsx homolog, Ind, similarly represses or stimulates transcription in a site-dependent manner via an autoregulatory enhancer containing a combination of monomer and homodimer sites. Integrating these findings, we test a model showing how the homodimer to monomer site ratio and the Gsx protein levels defines gene up-regulation versus down-regulation. Altogether, these data serve as a new paradigm for how cooperative homeodomain transcription factor binding can increase target specificity and alter regulatory outcomes.

CUT&RUN detects distinct DNA footprints of RNA polymerase II near the transcription start sites

CUT&RUN is a powerful tool to study protein-DNA interactions in vivo. DNA fragments cleaved by the targeted micrococcal nuclease identify the footprints of DNA-binding proteins on the chromatin. We performed CUT&RUN on human lung carcinoma cell line A549 maintained in a multi-well cell culture plate to profile RNA polymerase II. Long (> 270 bp) DNA fragments released by CUT&RUN corresponded to the bimodal peak around the transcription start sites, as previously seen with chromatin immunoprecipitation. However, we found that short (< 120 bp) fragments identify a well-defined peak localised at the transcription start sites. This distinct DNA footprint of short fragments, which constituted only about 5% of the total reads, suggests the transient positioning of RNA polymerase II before promoter-proximal pausing, which has not been detected in the physiological settings by standard chromatin immunoprecipitation. We showed that the positioning of the large-size-class DNA footprints around the short-fragment peak was associated with the directionality of transcription, demonstrating the biological significance of distinct CUT&RUN footprints of RNA polymerase II.

Chromatin Profiling in Mouse Embryonic Germ Cells by CUT&RUN

Cleavage under targets and release using nuclease (CUT&RUN) allows the chromatin profiling of proteins of interest for which specific antibodies are available. Because it is performed on intact chromatin in situ, CUT&RUN provides exceptional signal over background, making it an ideal choice for chromatin profiling on primary cells available at limited numbers. Here, we describe its application to the profiling of histone post-translational modifications in germ cells isolated from mouse embryos from 12.5 up to 18.5 days postfertilization. This approach can be applied to as low as 100 isolated germ cells, allowing the generation of multiple genome-wide profiles from the cells obtained from a single embryo.

PathSTORM: a road to early cancer detection

Disruption of chromatin structure could enable early carcinogenesis by facilitating malignant transformation. Using stochastic optical reconstruction microscopy optimized for pathological tissue (PathSTORM), we uncovered a gradual decompaction of higher-order chromatin folding through progressive stages of carcinogenesis. We demonstrated potential detection of pre-cancerous genomic architecture not easily discernible by conventional pathology.

Opposing Functions of BRD4 Isoforms in Breast Cancer

Bromodomain-containing protein 4 (BRD4) is a cancer therapeutic target in ongoing clinical trials disrupting primarily BRD4-regulated transcription programs. The role of BRD4 in cancer has been attributed mainly to the abundant long isoform (BRD4-L). Here we show, by isoform-specific knockdown and endogenous protein detection, along with transgene expression, the less abundant BRD4 short isoform (BRD4-S) is oncogenic while BRD4-L is tumor-suppressive in breast cancer cell proliferation and migration, as well as mammary tumor formation and metastasis. Through integrated RNA-seq, genome-wide ChIP-seq, and CUT&RUN association profiling, we identify the Engrailed-1 (EN1) homeobox transcription factor as a key BRD4-S coregulator, particularly in triple-negative breast cancer. BRD4-S and EN1 comodulate the extracellular matrix (ECM)-associated matrisome network, including type II cystatin gene cluster, mucin 5, and cathepsin loci, via enhancer regulation of cancer-associated genes and pathways. Our work highlights the importance of targeted therapies for the oncogenic, but not tumor-suppressive, activity of BRD4.

Genomic methods in profiling DNA accessibility and factor localization

Recent advancements in next-generation sequencing technologies and accompanying reductions in cost have led to an explosion of techniques to examine DNA accessibility and protein localization on chromatin genome-wide. Generally, accessible regions of chromatin are permissive for factor binding and are therefore hotspots for regulation of gene expression; conversely, genomic regions that are highly occupied by histone proteins are not permissive for factor binding and are less likely to be active regulatory regions. Identifying regions of differential accessibility can be useful to uncover putative gene regulatory regions, such as enhancers, promoters, and insulators. In addition, DNA-binding proteins, such as transcription factors that preferentially bind certain DNA sequences and histone proteins that form the core of the nucleosome, play essential roles in all DNA-templated processes. Determining the genomic localization of chromatin-bound proteins is therefore essential in determining functional roles, sequence motifs important for factor binding, and regulatory networks controlling gene expression. In this review, we discuss techniques for determining DNA accessibility and nucleosome positioning (DNase-seq, FAIRE-seq, MNase-seq, and ATAC-seq) and techniques for detecting and functionally characterizing chromatin-bound proteins (ChIP-seq, DamID, and CUT&RUN). These methods have been optimized to varying degrees of resolution, specificity, and ease of use. Here, we outline some advantages and disadvantages of these techniques, their general protocols, and a brief discussion of their development. Together, these complimentary approaches have provided an unparalleled view of chromatin architecture and functional gene regulation.

Peak calling by Sparse Enrichment Analysis for CUT&RUN chromatin profiling

BACKGROUND

CUT&RUN is an efficient epigenome profiling method that identifies sites of DNA binding protein enrichment genome-wide with high signal to noise and low sequencing requirements. Currently, the analysis of CUT&RUN data is complicated by its exceptionally low background, which renders programs designed for analysis of ChIP-seq data vulnerable to oversensitivity in identifying sites of protein binding.

RESULTS

Here we introduce Sparse Enrichment Analysis for CUT&RUN (SEACR), an analysis strategy that uses the global distribution of background signal to calibrate a simple threshold for peak calling. SEACR discriminates between true and false-positive peaks with near-perfect specificity from "gold standard" CUT&RUN datasets and efficiently identifies enriched regions for several different protein targets. We also introduce a web server ( http://seacr.fredhutch.org ) for plug-and-play analysis with SEACR that facilitates maximum accessibility across users of all skill levels.

CONCLUSIONS

SEACR is a highly selective peak caller that definitively validates the accuracy of CUT&RUN for datasets with known true negatives. Its ease of use and performance in comparison with existing peak calling strategies make it an ideal choice for analyzing CUT&RUN data.

Direct and widespread role for the nuclear receptor EcR in mediating the response to ecdysone in Drosophila

The ecdysone pathway was among the first experimental systems employed to study the impact of steroid hormones on the genome. In Drosophila and other insects, ecdysone coordinates developmental transitions, including wholesale transformation of the larva into the adult during metamorphosis. Like other hormones, ecdysone controls gene expression through a nuclear receptor, which functions as a ligand-dependent transcription factor. Although it is clear that ecdysone elicits distinct transcriptional responses within its different target tissues, the role of its receptor, EcR, in regulating target gene expression is incompletely understood. In particular, EcR initiates a cascade of transcription factor expression in response to ecdysone, making it unclear which ecdysone-responsive genes are direct EcR targets. Here, we use the larval-to-prepupal transition of developing wings to examine the role of EcR in gene regulation. Genome-wide DNA binding profiles reveal that EcR exhibits widespread binding across the genome, including at many canonical ecdysone response genes. However, the majority of its binding sites reside at genes with wing-specific functions. We also find that EcR binding is temporally dynamic, with thousands of binding sites changing over time. RNA-seq reveals that EcR acts as both a temporal gate to block precocious entry to the next developmental stage as well as a temporal trigger to promote the subsequent program. Finally, transgenic reporter analysis indicates that EcR regulates not only temporal changes in target enhancer activity but also spatial patterns. Together, these studies define EcR as a multipurpose, direct regulator of gene expression, greatly expanding its role in coordinating developmental transitions.

Profiling of Pluripotency Factors in Single Cells and Early Embryos

Cell fate decisions are governed by sequence-specific transcription factors (TFs) that act in small populations of cells within developing embryos. To understand their functions in vivo, it is important to identify TF binding sites in these cells. However, current methods cannot profile TFs genome-wide at or near the single-cell level. Here we adapt the cleavage under targets and release using nuclease (CUT&RUN) method to profile TFs in low cell numbers, including single cells and individual pre-implantation embryos. Single-cell experiments suggest that only a fraction of TF binding sites are occupied in most cells, in a manner broadly consistent with measurements of peak intensity from multi-cell studies. We further show that chromatin binding by the pluripotency TF NANOG is highly dependent on the SWI/SNF chromatin remodeling complex in individual blastocysts but not in cultured cells. Ultra-low input CUT&RUN (uliCUT&RUN) therefore enables interrogation of TF binding from rare cell populations of particular importance in development or disease.

High-Resolution Chromatin Profiling Using CUT&RUN

Determining the genomic location of DNA-binding proteins is essential to understanding their function. Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a powerful method for mapping protein-DNA interactions at high resolution. In CUT&RUN, a recombinant protein A-microccocal nuclease (pA-MN) fusion is recruited by an antibody targeting the chromatin protein of interest; this can be done with either uncrosslinked or formaldehyde-crosslinked cells. DNA fragments near sites of antibody binding are released from the insoluble bulk chromatin through endonucleolytic cleavage and used to build barcoded DNA-sequencing libraries that can be sequenced in pools of at least 30. Therefore, CUT&RUN provides an alternative to ChIP-seq approaches for mapping chromatin proteins, which typically have relatively high signal-to-noise ratios, while using fewer cells and at a lower cost. Here, we describe the methods for performing CUT&RUN, generating DNA-sequencing libraries, and analyzing the resulting datasets. © 2019 by John Wiley & Sons, Inc.

Automated in situ chromatin profiling efficiently resolves cell types and gene regulatory programs

BACKGROUND

Our understanding of eukaryotic gene regulation is limited by the complexity of protein-DNA interactions that comprise the chromatin landscape and by inefficient methods for characterizing these interactions. We recently introduced CUT&RUN, an antibody-targeted nuclease cleavage method that profiles DNA-binding proteins, histones and chromatin-modifying proteins in situ with exceptional sensitivity and resolution.

RESULTS

Here, we describe an automated CUT&RUN platform and apply it to characterize the chromatin landscapes of human cells. We find that automated CUT&RUN profiles of histone modifications crisply demarcate active and repressed chromatin regions, and we develop a continuous metric to identify cell-type-specific promoter and enhancer activities. We test the ability of automated CUT&RUN to profile frozen tumor samples and find that our method readily distinguishes two pediatric glioma xenografts by their subtype-specific gene expression programs.

CONCLUSIONS

The easy, cost-effective workflow makes automated CUT&RUN an attractive tool for high-throughput characterization of cell types and patient samples.

RSC-Associated Subnucleosomes Define MNase-Sensitive Promoters in Yeast

The classic view of nucleosome organization at active promoters is that two well-positioned nucleosomes flank a nucleosome-depleted region (NDR). However, this view has been recently disputed by contradictory reports as to whether wider (≳150 bp) NDRs instead contain unstable, micrococcal nuclease-sensitive ("fragile") nucleosomal particles. To determine the composition of fragile particles, we introduce CUT&RUN.ChIP, in which targeted nuclease cleavage and release is followed by chromatin immunoprecipitation. We find that fragile particles represent the occupancy of the RSC (remodeling the structure of chromatin) nucleosome remodeling complex and RSC-bound, partially unwrapped nucleosomal intermediates. We also find that general regulatory factors (GRFs) bind to partially unwrapped nucleosomes at these promoters. We propose that RSC binding and its action cause nucleosomes to unravel, facilitate subsequent binding of GRFs, and constitute a dynamic cycle of nucleosome deposition and clearance at the subset of wide Pol II promoter NDRs.

Direct Promoter Repression by BCL11A Controls the Fetal to Adult Hemoglobin Switch

Fetal hemoglobin (HbF, α₂γ₂) level is genetically controlled and modifies severity of adult hemoglobin (HbA, α₂β₂) disorders, sickle cell disease, and β-thalassemia. Common genetic variation affects expression of BCL11A, a regulator of HbF silencing. To uncover how BCL11A supports the developmental switch from γ- to β- globin, we use a functional assay and protein binding microarray to establish a requirement for a zinc-finger cluster in BCL11A in repression and identify a preferred DNA recognition sequence. This motif appears in embryonic and fetal-expressed globin promoters and is duplicated in γ-globin promoters. The more distal of the duplicated motifs is mutated in individuals with hereditary persistence of HbF. Using the CUT&RUN approach to map protein binding sites in erythroid cells, we demonstrate BCL11A occupancy preferentially at the distal motif, which can be disrupted by editing the promoter. Our findings reveal that direct γ-globin gene promoter repression by BCL11A underlies hemoglobin switching.