RNA Polymerase II hypertranscription in cancer FFPE samples.

Hypertranscription is common in human cancers and predicts poor prognosis. However detection of hypertranscription is indirect, relying on accurately quantifying mRNA levels and estimating cell numbers. Previously, we introduced FFPE-CUTAC, a genome-wide method for mapping RNA Polymerase II (RNAPII) in formalin-fixed paraffin-embedded (FFPE) sections. Here we use FFPE-CUTAC to demonstrate genome-wide hypertranscription both in transgene-driven mouse gliomas and in assorted human tumors at active regulatory elements and replication-coupled histone genes with reduced mitochondrial DNA abundance. FFPE-CUTAC identified RNAPII-bound regulatory elements shared among diverse cancers and readily categorized human tumors despite using very small samples and low sequencing depths. Remarkably, RNAPII FFPE-CUTAC identified de novo and precisely mapped HER2 amplifications punctuated by likely selective sweeps including genes encoding direct positive regulators of RNAPII itself. Our results demonstrate that FFPE-CUTAC measurements of hypertranscription and classifications of tumors using small sections provides an affordable and sensitive genome-wide strategy for personalized medicine.

Epigenomic analysis of formalin-fixed paraffin-embedded samples by CUT&Tag.

For more than a century, formalin-fixed paraffin-embedded (FFPE) sample preparation has been the preferred method for long-term preservation of biological material. However, the use of FFPE samples for epigenomic studies has been difficult because of chromatin damage from long exposure to high concentrations of formaldehyde. Previously, we introduced Cleavage Under Targeted Accessible Chromatin (CUTAC), an antibody-targeted chromatin accessibility mapping protocol based on CUT&Tag. Here we show that simple modifications of our CUTAC protocol either in single tubes or directly on slides produce high-resolution maps of paused RNA Polymerase II at enhancers and promoters using FFPE samples. We find that transcriptional regulatory element differences produced by FFPE-CUTAC distinguish between mouse brain tumors and identify and map regulatory element markers with high confidence and precision, including microRNAs not detectable by RNA-seq. Our simple workflows make possible affordable epigenomic profiling of archived biological samples for biomarker identification, clinical applications and retrospective studies.

A giant virus genome is densely packaged by stable nucleosomes within virions.

The two doublet histones of Marseillevirus are distantly related to the four eukaryotic core histones and wrap 121 base pairs of DNA to form remarkably similar nucleosomes. By permeabilizing Marseillevirus virions and performing genome-wide nuclease digestion, chemical cleavage, and mass spectrometry assays, we find that the higher-order organization of Marseillevirus chromatin fundamentally differs from that of eukaryotes. Marseillevirus nucleosomes fully protect DNA within virions as closely abutted 121-bp DNA-wrapped cores without linker DNA or phasing along genes. Likewise, we observed that nucleosomes reconstituted onto multi-copy tandem repeats of a nucleosome-positioning sequence are tightly packed. Dense promiscuous packing of fully wrapped nucleosomes rather than "beads on a string" with genic punctuation represents a distinct mode of DNA packaging by histones. We suggest that doublet histones have evolved for viral genome protection and may resemble an early stage of histone differentiation leading to the eukaryotic octameric nucleosome.

Simplified Epigenome Profiling Using Antibody-tethered Tagmentation.

We previously introduced Cleavage Under Targets & Tagmentation (CUT&Tag), an epigenomic profiling method in which antibody tethering of the Tn5 transposase to a chromatin epitope of interest maps specific chromatin features in small samples and single cells. With CUT&Tag, intact cells or nuclei are permeabilized, followed by successive addition of a primary antibody, a secondary antibody, and a chimeric Protein A-Transposase fusion protein that binds to the antibody. Addition of Mg++ activates the transposase and inserts sequencing adapters into adjacent DNA in situ. We have since adapted CUT&Tag to also map chromatin accessibility by simply modifying the transposase activation conditions when using histone H3K4me2, H3K4me3, or Serine-5-phosphorylated RNA Polymerase II antibodies. Using these antibodies, we redirect the tagmentation of accessible DNA sites to produce chromatin accessibility maps with exceptionally high signal-to-noise and resolution. All steps from nuclei to amplified sequencing-ready libraries are performed in single PCR tubes using non-toxic reagents and inexpensive equipment, making our simplified strategy for simultaneous chromatin profiling and accessibility mapping suitable for the lab, home workbench, or classroom.

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation.

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA ('tagmentation'). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here, we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

Cells keep their DNA tidy by wrapping it into structures called nucleosomes. Each of these structures contains a short section of DNA wound around a cluster of proteins called histones. Not only do nucleosomes keep the genetic code organized, they also control whether the proteins that can switch genes on or off have access to the DNA. When genes turn on, the nucleosomes unwrap, exposing sections of genetic code called 'gene regulatory elements'. These elements attract the proteins that help read and copy nearby genes so the cell can make new proteins. Determining which regulatory elements are exposed at any given time can provide useful information about what is happening inside a cell, but the procedure can be expensive. The most popular way to map which regulatory elements are exposed is using a technique called Assay for Transposase-Accessible Chromatin using sequencing, or ATAC-seq for short. The 'transposase' in the acronym is an enzyme that cuts areas of DNA that are not wound around histones and prepares them for detection by DNA sequencing. Unfortunately, the data from ATAC-seq are often noisy (there are random factors that produce a signal that is detected but is not a 'real' result), so more sequencing is required to differentiate between real signal and noise, increasing the expense of ATAC-seq experiments. Furthermore, although ATAC-seq can identify unspooled sections of DNA, it cannot provide a direct connection between active genes and unwrapped DNA. To find the link between unspooled DNA and active genes, Henikoff et al. adapted a technique called CUT&Tag. Like ATAC-seq, it also uses transposases to cut the genome, but it allows more control over where the cuts occur. When genes are switched on, the proteins reading them leave chemical marks on the histones they pass. CUT&Tag attaches a transposase to a molecule that recognizes and binds to those marks. This allowed Henikoff et al. to guide the transposases to unspooled regions of DNA bordering active genes. The maps of gene regulatory elements produced using this method were the same as the best ATAC-seq maps. And, because the transposases could only access gaps near active genes, the data provided evidence that genes switching on leads to regulatory elements in the genome unwrapping. This new technique is simple enough that Henikoff et al. were able to perform it from home on the countertop of a laundry room. By tethering the transposases to histone marks it was possible to detect unspooled DNA that was active more efficiently than with ATAC-seq. This lowers laboratory costs by reducing the cost of DNA sequencing, and may also improve the detection of gaps between nucleosomes in single cells.

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.

Histone deposition pathways determine the chromatin landscapes of H3.1 and H3.3 K27M oncohistones.

Lysine 27-to-methionine (K27M) mutations in the H3.1 or H3.3 histone genes are characteristic of pediatric diffuse midline gliomas (DMGs). These oncohistone mutations dominantly inhibit histone H3K27 trimethylation and silencing, but it is unknown how oncohistone type affects gliomagenesis. We show that the genomic distributions of H3.1 and H3.3 oncohistones in human patient-derived DMG cells are consistent with the DNAreplication-coupled deposition of histone H3.1 and the predominant replication-independent deposition of histone H3.3. Although H3K27 trimethylation is reduced for both oncohistone types, H3.3K27M-bearing cells retain some domains, and only H3.1K27M-bearing cells lack H3K27 trimethylation. Neither oncohistone interferes with PRC2 binding. Using Drosophila as a model, we demonstrate that inhibition of H3K27 trimethylation occurs only when H3K27M oncohistones are deposited into chromatin and only when expressed in cycling cells. We propose that oncohistones inhibit the H3K27 methyltransferase as chromatin patterns are being duplicated in proliferating cells, predisposing them to tumorigenesis.

Improved CUT&RUN chromatin profiling tools.

Previously, we described a novel alternative to chromatin immunoprecipitation, CUT&RUN, in which unfixed permeabilized cells are incubated with antibody, followed by binding of a protein A-Micrococcal Nuclease (pA/MNase) fusion protein (Skene and Henikoff, 2017). Here we introduce three enhancements to CUT&RUN: A hybrid protein A-Protein G-MNase construct that expands antibody compatibility and simplifies purification, a modified digestion protocol that inhibits premature release of the nuclease-bound complex, and a calibration strategy based on carry-over of E. coli DNA introduced with the fusion protein. These new features, coupled with the previously described low-cost, high efficiency, high reproducibility and high-throughput capability of CUT&RUN make it the method of choice for routine epigenomic profiling.

CUT&Tag for efficient epigenomic profiling of small samples and single cells.

Many chromatin features play critical roles in regulating gene expression. A complete understanding of gene regulation will require the mapping of specific chromatin features in small samples of cells at high resolution. Here we describe Cleavage Under Targets and Tagmentation (CUT&Tag), an enzyme-tethering strategy that provides efficient high-resolution sequencing libraries for profiling diverse chromatin components. In CUT&Tag, a chromatin protein is bound in situ by a specific antibody, which then tethers a protein A-Tn5 transposase fusion protein. Activation of the transposase efficiently generates fragment libraries with high resolution and exceptionally low background. All steps from live cells to sequencing-ready libraries can be performed in a single tube on the benchtop or a microwell in a high-throughput pipeline, and the entire procedure can be performed in one day. We demonstrate the utility of CUT&Tag by profiling histone modifications, RNA Polymerase II and transcription factors on low cell numbers and single cells.

Targeted in situ genome-wide profiling with high efficiency for low cell numbers.

Cleavage under targets and release using nuclease (CUT&RUN) is an epigenomic profiling strategy in which antibody-targeted controlled cleavage by micrococcal nuclease releases specific protein-DNA complexes into the supernatant for paired-end DNA sequencing. As only the targeted fragments enter into solution, and the vast majority of DNA is left behind, CUT&RUN has exceptionally low background levels. CUT&RUN outperforms the most widely used chromatin immunoprecipitation (ChIP) protocols in resolution, signal-to-noise ratio and depth of sequencing required. In contrast to ChIP, CUT&RUN is free of solubility and DNA accessibility artifacts and has been used to profile insoluble chromatin and to detect long-range 3D contacts without cross-linking. Here, we present an improved CUT&RUN protocol that does not require isolation of nuclei and provides high-quality data when starting with only 100 cells for a histone modification and 1,000 cells for a transcription factor. From cells to purified DNA, CUT&RUN requires less than a day at the laboratory bench and requires no specialized skills.

A unique chromatin complex occupies young α-satellite arrays of human centromeres.

The intractability of homogeneous α-satellite arrays has impeded understanding of human centromeres. Artificial centromeres are produced from higher-order repeats (HORs) present at centromere edges, although the exact sequences and chromatin conformations of centromere cores remain unknown. We use high-resolution chromatin immunoprecipitation (ChIP) of centromere components followed by clustering of sequence data as an unbiased approach to identify functional centromere sequences. We find that specific dimeric α-satellite units shared by multiple individuals dominate functional human centromeres. We identify two recently homogenized α-satellite dimers that are occupied by precisely positioned CENP-A (cenH3) nucleosomes with two ~100-base pair (bp) DNA wraps in tandem separated by a CENP-B/CENP-C-containing linker, whereas pericentromeric HORs show diffuse positioning. Precise positioning is largely maintained, whereas abundance decreases exponentially with divergence, which suggests that young α-satellite dimers with paired ~100-bp particles mediate evolution of functional human centromeres. Our unbiased strategy for identifying functional centromeric sequences should be generally applicable to tandem repeat arrays that dominate the centromeres of most eukaryotes.

The budding yeast Centromere DNA Element II wraps a stable Cse4 hemisome in either orientation in vivo.

In budding yeast, a single cenH3 (Cse4) nucleosome occupies the ∼120-bp functional centromere, however conflicting structural models for the particle have been proposed. To resolve this controversy, we have applied H4S47C-anchored cleavage mapping, which reveals the precise position of histone H4 in every nucleosome in the genome. We find that cleavage patterns at centromeres are unique within the genome and are incompatible with symmetrical structures, including octameric nucleosomes and (Cse4/H4)2 tetrasomes. Centromere cleavage patterns are compatible with a precisely positioned core structure, one in which each of the 16 yeast centromeres is occupied by oppositely oriented Cse4/H4/H2A/H2B hemisomes in two rotational phases within the population. Centromere-specific hemisomes are also inferred from distances observed between closely-spaced H4 cleavages, as predicted from structural modeling. Our results indicate that the orientation and rotational position of the stable hemisome at each yeast centromere is not specified by the functional centromere sequence. DOI: http://dx.doi.org/10.7554/eLife.01861.001.

The CentO satellite confers translational and rotational phasing on cenH3 nucleosomes in rice centromeres.

Plant and animal centromeres comprise megabases of highly repeated satellite sequences, yet centromere function can be specified epigenetically on single-copy DNA by the presence of nucleosomes containing a centromere-specific variant of histone H3 (cenH3). We determined the positions of cenH3 nucleosomes in rice (Oryza sativa), which has centromeres composed of both the 155-bp CentO satellite repeat and single-copy non-CentO sequences. We find that cenH3 nucleosomes protect 90-100 bp of DNA from micrococcal nuclease digestion, sufficient for only a single wrap of DNA around the cenH3 nucleosome core. cenH3 nucleosomes are translationally phased with 155-bp periodicity on CentO repeats, but not on non-CentO sequences. CentO repeats have an ∼10-bp periodicity in WW dinucleotides and in micrococcal nuclease cleavage, providing evidence for rotational phasing of cenH3 nucleosomes on CentO and suggesting that satellites evolve for translational and rotational stabilization of centromeric nucleosomes.

"Point" centromeres of Saccharomyces harbor single centromere-specific nucleosomes.

The "point" centromere of budding yeast is genetically defined by an ≈ 125-bp sequence. Recent fluorescence measurements of kinetochore clusters have suggested that this sequence specifies multiple centromere histone 3 (CenH3) nucleosomes. However, high-resolution mapping demonstrates that there is only one CenH3 nucleosome per centromere, providing biochemical confirmation of the point centromere model.

Tripartite organization of centromeric chromatin in budding yeast.

The centromere is the genetic locus that organizes the proteinaceous kinetochore and is responsible for attachment of the chromosome to the spindle at mitosis and meiosis. In most eukaryotes, the centromere consists of highly repetitive DNA sequences that are occupied by nucleosomes containing the CenH3 histone variant, whereas in budding yeast, a ∼120-bp centromere DNA element (CDE) that is sufficient for centromere function is occupied by a single right-handed histone variant CenH3 (Cse4) nucleosome. However, these in vivo observations are inconsistent with in vitro evidence for left-handed octameric CenH3 nucleosomes. To help resolve these inconsistencies, we characterized yeast centromeric chromatin at single base-pair resolution. Intact particles containing both Cse4 and H2A are precisely protected from micrococcal nuclease over the entire CDE of all 16 yeast centromeres in both solubilized chromatin and the insoluble kinetochore. Small DNA-binding proteins protect CDEI and CDEIII and delimit the centromeric nucleosome to the ∼80-bp CDEII, only enough for a single DNA wrap. As expected for a tripartite organization of centromeric chromatin, loss of Cbf1 protein, which binds to CDEI, both reduces the size of the centromere-protected region and shifts its location toward CDEIII. Surprisingly, Cse4 overproduction caused genome-wide misincorporation of nonfunctional CenH3-containing nucleosomes that protect ∼135 base pairs and are preferentially enriched at sites of high nucleosome turnover. Our detection of two forms of CenH3 nucleosomes in the yeast genome, a singly wrapped particle at the functional centromere and octamer-sized particles on chromosome arms, reconcile seemingly conflicting in vivo and in vitro observations.

Epigenome characterization at single base-pair resolution.

We have combined standard micrococcal nuclease (MNase) digestion of nuclei with a modified protocol for constructing paired-end DNA sequencing libraries to map both nucleosomes and subnucleosome-sized particles at single base-pair resolution throughout the budding yeast genome. We found that partially unwrapped nucleosomes and subnucleosome-sized particles can occupy the same position within a cell population, suggesting dynamic behavior. By varying the time of MNase digestion, we have been able to observe changes that reflect differential sensitivity of particles, including the eviction of nucleosomes. To characterize DNA-binding features of transcription factors, we plotted the length of each fragment versus its position in the genome, which defined the minimal protected region of each factor. This process led to the precise mapping of protected and exposed regions at and around binding sites, and also determination of the degree to which they are flanked by phased nucleosomes and subnucleosome-sized particles. Our protocol and mapping method provide a general strategy for epigenome characterization, including nucleosome phasing and dynamics, ATP-dependent nucleosome remodelers, and transcription factors, from a single-sequenced sample.

Identification of functional elements and regulatory circuits by Drosophila modENCODE.

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project.

We systematically generated large-scale data sets to improve genome annotation for the nematode Caenorhabditis elegans, a key model organism. These data sets include transcriptome profiling across a developmental time course, genome-wide identification of transcription factor-binding sites, and maps of chromatin organization. From this, we created more complete and accurate gene models, including alternative splice forms and candidate noncoding RNAs. We constructed hierarchical networks of transcription factor-binding and microRNA interactions and discovered chromosomal locations bound by an unusually large number of transcription factors. Different patterns of chromatin composition and histone modification were revealed between chromosome arms and centers, with similarly prominent differences between autosomes and the X chromosome. Integrating data types, we built statistical models relating chromatin, transcription factor binding, and gene expression. Overall, our analyses ascribed putative functions to most of the conserved genome.

H2A.Z nucleosomes enriched over active genes are homotypic.

Nucleosomes that contain the histone variant H2A.Z are enriched around transcriptional start sites, but the mechanistic basis for this enrichment is unknown. A single octameric nucleosome can contain two H2A.Z histones (homotypic) or one H2A.Z and one canonical H2A (heterotypic). To elucidate the function of H2A.Z, we generated high-resolution maps of homotypic and heterotypic Drosophila H2A.Z (H2Av) nucleosomes. Although homotypic and heterotypic H2A.Z nucleosomes mapped throughout most of the genome, homotypic nucleosomes were enriched and heterotypic nucleosomes were depleted downstream of active promoters and intron-exon junctions. The distribution of homotypic H2A.Z nucleosomes resembled that of classical active chromatin and showed evidence of disruption during transcriptional elongation. Both homotypic H2A.Z nucleosomes and classical active chromatin were depleted downstream of paused polymerases. Our results suggest that H2A.Z enrichment patterns result from intrinsic structural differences between heterotypic and homotypic H2A.Z nucleosomes that follow disruption during transcriptional elongation.

Genome-wide kinetics of nucleosome turnover determined by metabolic labeling of histones.

Nucleosome disruption and replacement are crucial activities that maintain epigenomes, but these highly dynamic processes have been difficult to study. Here, we describe a direct method for measuring nucleosome turnover dynamics genome-wide. We found that nucleosome turnover is most rapid over active gene bodies, epigenetic regulatory elements, and replication origins in Drosophila cells. Nucleosomes turn over faster at sites for trithorax-group than polycomb-group protein binding, suggesting that nucleosome turnover differences underlie their opposing activities and challenging models for epigenetic inheritance that rely on stability of histone marks. Our results establish a general strategy for studying nucleosome dynamics and uncover nucleosome turnover differences across the genome that are likely to have functional importance for epigenome maintenance, gene regulation, and control of DNA replication.