Rapid genome-scale mapping of chromatin accessibility in tissue
Epigenetics & Chromatin volume 5, Article number: 10 (2012)
The challenge in extracting genome-wide chromatin features from limiting clinical samples poses a significant hurdle in identification of regulatory marks that impact the physiological or pathological state. Current methods that identify nuclease accessible chromatin are reliant on large amounts of purified nuclei as starting material. This complicates analysis of trace clinical tissue samples that are often stored frozen. We have developed an alternative nuclease based procedure to bypass nuclear preparation to interrogate nuclease accessible regions in frozen tissue samples.
Here we introduce a novel technique that specifically identifies Tissue Accessible Chromatin (TACh). The TACh method uses pulverized frozen tissue as starting material and employs one of the two robust endonucleases, Benzonase or Cyansase, which are fully active under a range of stringent conditions such as high levels of detergent and DTT. As a proof of principle we applied TACh to frozen mouse liver tissue. Combined with massive parallel sequencing TACh identifies accessible regions that are associated with euchromatic features and accessibility at transcriptional start sites correlates positively with levels of gene transcription. Accessible chromatin identified by TACh overlaps to a large extend with accessible chromatin identified by DNase I using nuclei purified from freshly isolated liver tissue as starting material. The similarities are most pronounced at highly accessible regions, whereas identification of less accessible regions tends to be more divergence between nucleases. Interestingly, we show that some of the differences between DNase I and Benzonase relate to their intrinsic sequence biases and accordingly accessibility of CpG islands is probed more efficiently using TACh.
The TACh methodology identifies accessible chromatin derived from frozen tissue samples. We propose that this simple, robust approach can be applied across a broad range of clinically relevant samples to allow demarcation of regulatory elements of considerable prognostic significance.
The coupling of next-generation sequencing methodologies with classical enzymatic chromatin digestion approaches (eg. DNase I and MNase) has provided global, high-resolution information on chromatin features such as nucleosome positioning and chromatin accessibility [1, 2]. DNase-Seq has emerged as a powerful genome-wide tool to identify and characterize chromatin transitions (DNase I hypersensitive sites or DHS) in regulatory regions across a range of biological processes and cell lines [1, 3–5]. Current methodologies to generate genome-wide DHS data start with large numbers of fresh cells from which sizeable quantities of high quality nuclei need to be generated . This has restricted application of the DNase I approach to cells in culture and freshly isolated tissue and has prevented the use of this technology on clinical tissue samples that are typically scarce, stored frozen and from which nuclei are difficult to obtain.
The use of DNase I as a probe of chromatin accessibility was initially described three decades ago [7, 8], where partial digestion of chromatin with DNase I indentified DHS in promoters of heat shock genes. The nucleosomal steric hindrance of DNase I access to DNA is similar to that of transcription factors (TF) and the degree of DNase I accessibility correlates with the level of TF occupancy . Thus DNase I hypersensitivity is considered an independent, unbiased probe for TF accessibility. Today, the primary commercial sources of DNase I are either recombinant proteins or purified enzymes from bovine pancreas. However, DNase I is known to have drawbacks that restrict its broad-spectrum use. Specifically, it is known to have a narrow effective concentration window, requiring fastidious titration of enzyme and chromatin to obtain optimal partial digestion [6, 7]. Moreover, DNase I is inhibited by high levels of actin  and because the actin content varies amongst cell and tissue types, each biological sample must be individually titrated. Preparation of nuclei prior to digestion reduces actin contamination, while in situ NP40 lysis leaving adherent nuclei also allows DNase I treatment of cells grown as a monolayer .
In order to define accessible chromatin compartments in samples derived from frozen tissue, where cell numbers are unknown and nuclear preparation is problematic, we hypothesized that digestion of chromatin with more robust nucleases that function in a coarse environment as well as, over broad concentration ranges, might enable the determination of chromatin accessibility in tissue samples. Benzonase, a recombinant endonuclease derived from Serratia marcescens, is composed of two identical subunits of 30 kDa, requires divalent cations for full activity and is generally used to clear cellular protein extracts of DNA and RNA prior to downstream analysis . It digests DNA and RNA efficiently under a range of conditions to nucleotides of 2–5 base pairs. Although the enzyme is able to cleave DNA at all positions, it has been reported to have a relative preference for GC rich regions over dA/dT tracts . Cyanase is a less described non-Serratia recombinant endonuclease that, like Benzonase, efficiently digests DNA and RNA under harsh conditions. Here we introduce a novel technique that combines rapid processing of frozen tissue using Benzonase and Cyanase to specifically identify Tissue Accessible Chromatin (TACh) from frozen specimens. Accessible regions identified with TACh correlate with features of euchromatin and levels of transcription, suggesting that these accessible regions are indeed regulatory. We propose that TACh will be a valuable tool to identify the physiological or pathological regulatory features of chromatin from clinical materials.
Benzonase and Cyanase as probes for chromatin accessibility
Accessible chromatin has traditionally been identified by DNase I digestion of chromatin using nuclei as starting material (Figure 1A, left panel). Although nuclei can be very efficiently purified from cell lines and fresh tissue within one to two hours, such purification requires disassociation of cells, and washing by centrifugation, conditions that could modify signaling to the nucleus or allow leaching of chromatin-bound components, potentially altering nuclear structures. Extracting nuclei from frozen tissue samples is even more cumbersome and complicated. Thus, in order to minimize the time between the snap freezing of tissue and enzymatic digestion, we have developed a method that avoids nuclear preparation and uses a different endonucleaese, Benzonase or Cyanase, to digest accessible chromatin-embedded DNA (Figure 1A, right panel).
To set a standard for the fidelity of Benzonase and Cyanase as a probe for chromatin accessibility, we initially performed a conventional nuclease hypersensitivity assay using cultured cells. Human promyelocytic leukemia cells (HL-60) grown in suspension were isolated, resuspended hypotonic buffer and incubated with increasing concentrations of Benzonase and Cyanase. Accessible regions at the c-myc promoter were compared using indirect-end labeling and Southern blotting as previously described . We show that Benzonase and Cyanase yielded the same pattern of hypersensitive regions expected for DNase I , demonstrating that Benzonase and Cyanase are useful probes for chromatin accessibility (Figure 1B).
Identification of accessible chromatin in frozen tissue
To test whether Benzonase and Cyanase can interrogate nuclease accessible regions in chromatin from frozen tissue, whole livers from C57BL/6 mice were isolated and frozen immediately in liquid nitrogen. We initially compared different procedures to prepare frozen tissues amenable for nuclease treatment without disrupting chromatin integrity. We found that rapid pulverization of frozen tissue into a fine powder prior to digestion results in the best signal to noise ratio (data not shown). Pulverization was performed on dry ice with equipment pre-cooled in liquid nitrogen and pulverized tissue was stored as aliquots (Figure 1A, right panel). For digestion, pulverized tissue was directly resuspended in a hypotonic digestion buffer and subsequently incubated with Benzonase or Cyanase at different concentrations. DNA fragments from chromatin digested with 0.25U/ml, 1U/ml and 4U/ml of Benzonase or Cyanase were isolated as previously described , sequenced to a depth of 20–30 million uniquely aligning tags (Additional file 1: Figure S 1) and accessible regions were identified as previously described . At all three enzyme concentrations, Benzonase and Cyanase revealed a robust set of nuclease hypersensitive regions in the genome as exemplified by the tyrosine aminotransferase (Tat) gene, a highly expressed liver-specific gene (Figure 2A). Reflecting the usage of frozen tissue, a larger portion of tags generated from TACh aligns with the mitochondrial genome compared to tags generated by DNase I digestion of chromatin from nuclei (Additional file 1: Figure S 1). Genome-wide, the tag density of hotspots identified with 0.25U of Benzonase resembled the tag density of hotspots identified with 1U of Benzonase, correlation coefficient of 0.951 (Figure 2B), suggesting that a fourfold increase in enzyme concentration identifies the same spectrum of hotspots. When the enzyme concentration was increased an additional fourfold to 4U/ml, although the most intense hotspots were reduced in intensity the overall correlation was still 82% with 1U/ml enzyme (Figure 2C). Similar patterns were seen using Cyanase (Additional file 1: Figure S 2) and remarkably at the different enzyme concentrations both enzymes performed very similarly (Additional file 1: Figure S 3). When data was combined from all three concentrations of Benzonase and Cyanase, each identified ~50,000 hotspots with remarkably similar tag densities (Figure 2D, correlation coefficient of 0.973) and an 87% overlap (Figure 2E). Thus in contrast to the narrow concentration windows of DNase I needed for optimal digestion , the hotspot patterns obtained with Benzonase-Cyanase were robustly conserved over a 16-fold range in enzyme concentration. This is a notable advantage for the use of Benzonase or Cyanase with frozen tissue or cell samples when exact cell counts are unavailable.
Correlation of Benzonase hotspots with euchromatin and TSS of active genes
To verify that the TACh procedure identifies accessible regions associated with regulatory elements of gene expression, we mapped the distribution of hotspots in distal upstream regions, promoters, introns, exons and downstream regions, and correlated hotspot intensity with previously mapped histone modifications  and nucleosome occupancy  in mouse liver tissue (Figure 3). The hotspots with the highest tag densities were found primarily at promoters, whereas the weaker hotspots located mainly in distal upstream and intronic regions similar to enhancers and other regulatory elements (Figure 3A). In agreement, the most intense hotspots were flanked by the promoter-specific H3K4me3 histone modification compared to less intense hotspots (Figure 3B). Moreover the most intense hotspots were also the most sensitive to MNase digestion, suggesting that these regions are either nucleosome –free or occupied by highly mobile nucleosomes flanked by H3K4me3-modified nucleosomes (Figure 3C). H3K4me1, present at promoters as well as enhancers, was enriched at both strong and weak Benzonase hotspots (Figure 3D), while H3K27me3, associated with heterochromatic regions, was deficient at Benzonase hotspots (Figure 3E). Thus Benzonase accessibility is associated with euchromatic features, demonstrating that the TACh method identifies accessible regulatory regions of the genome from frozen tissue.
Transcriptional start sites (TSS) of active genes are occupied by highly mobile nucleosomes  and are thus highly accessible to DNase I. In agreement, more than 90% of genes producing more than 16 transcripts were marked by Benzonase and Cyanase hotspots at the TSS (Figure 4A); conversely, only 30% of TSSs of inactive genes contained Benzonase-Cyanase hotspots (Figure 4A). Moreover, active genes had an overall increase in Benzonase and Cyanase accessibility at TSSs, compared to less active or silent genes (Figure 4B). Furthermore, when TSSs were binned into deciles according to the abundance of their gene transcripts, measured by previously published RNA-seq data , a positive correlation of gene transcription with the degree of Benzonase and Cyanase accessibility was observed (Figure 4C).
Benzonase and Cyanase accessible regions overlap with DNase I hotspots
To validate that accessible regions identified by the TACh are indeed bona fide nuclease hypersensitive sites, we mapped DNase I accessible regions using nuclei purified from fresh liver tissue (Figure 1A, left panel). Benzonase, Cyanase and DNase I accessible regions were largely similar at the Tat gene locus. However, we also observed features unique to each nuclease (Figure 5A). Using identical parameters to identify hotspots we detected ~63,000 DNase I hotspots (Figure 5B) which combined with the Benzonase and Cyanase data, identified a total of ~76,000 hotspots. Of these 28% was unique to DNase I, 52% was shared among the three enzymes and 20% was unique to Benzonase-Cyanase (Figure 5B). Parsing nuclease hotspots into quartiles according to tag density, we observed that 62% of the weakest DNase I hotspots were unique whereas 97% of the strongest hotspots overlapped with Benzonase-Cyanase hotpots (Figure 5C). Likewise ~50% of the least intense Benzonase and Cyanase hotspots were unique while close to all of the most intense hotspots overlapped with DNase I hotspots (Figure 5D-E). This suggests that most of highly accessible regions are identified by all enzymes whereas less accessible regions may be unique to particular nucleases. Alternatively many of these less accessible unique regions may have their origin in background digestion by the nucleases and may not be significant. Furthermore Dnase I unique hotspots were preferentially found at introns and distal regions in contrast to Benzonase-Cyanase hotspots which were enriched at promoters (Figure 5F).
Sequence bias for endonucleases
The variation observed between identified hotspots by the nucleases might be explained by the intrinsic methodological differences between TACh and the DNase I-based assays. Specifically, TACh is performed in intact cells with minimum manipulation prior to digestion, while the DNase I assay is performed on nuclei that take at least an hour to process. Alternatively, differences between DNase I, Benzonase and Cyanase can be a consequence of sequence specificity for DNA recognition and cleavage by each of the endonucleases. Benzonase preferentially digests dsDNA enriched for Gs and Cs  while DNase I prefers Ts . In agreement with the base-specificity explanation, Benzonase and Cyanase unique hotspots at the Tat loci overlapped with a GC-rich CpG island proximal to the Marveld3 gene, whereas DNase I unique hotspots overlapped with low-GC regions (Figure 5A). To explore sequence selectivity for cleavage genome-wide, we analyzed the sequence immediately upstream and downstream of all tags sequenced after digestion with DNase I or Benzonase. As shown in Figure 6A, the sequence tags yielded by Benzonase digestion were enriched for Gs at their 5’ ends, whereas the tags produced by DNase I digestion were enriched for 5’ Ts, suggesting that Benzonase-Cyanase preferentially cleaved at accessible DNA regions with high GC content and DNase I at accessible regions with high AT content. In agreement, the hotspots unique to Benzonase-Cyanase had higher overall GC content compared to surrounding regions (Figure 6D) or DNase I unique hotspots (Figure 6B-C). In contrast, DNase I unique hotspots had higher AT content than either neighboring regions or Benzonase-Cyanase hotspots (Additional file 1: Figure S 4). Common hotspots identified by all three enzymes had intermediate GC contents (Figure 6C). Consistent with the preference of Benzonase-Cyanase for high GC content regions, about 23% of hotspots uniquely identified by Benzonase and Cyanase were within CpG islands, whereas less than 1% of hotspots unique to DNase I mapped to CpG islands (Figure 6E). Benzonase-Cyanase specific hotspots were found to a greater extent within promoters (Figure 5F) and the average GC content of these unique Benzonase-Cyanase accessible promoters was 61% with 57% of them overlapping with CpG islands.
Nuclease accessible sites are maintained by the targeted recruitment of remodeling complexes by transcription factors . The transcription factor FOXA2 has been shown to bind AT rich DNA motifs , whereas SREBP has been suggested to bind regions together with SP1  with a resulting preference for GC-rich motifs embedded in CpG islands. Accordingly, the tags sequenced following DNase I digestion were recovered at higher rates from FOXA2 binding sites compared to the tags sequenced using Benzonase-Cyanase digestion (Figure 6F). In contrast, tags from Benzonase-Cyanase digestion were recovered more efficiently at SREBP binding sites (Figure 6G), emphasizing that the differences in base-selectivity at cleavage sites for these enzymes translates into differentially efficient capacities to interrogate the compartments associated with different transcription factors. Though largely similar, Benzonase-Cyanase and DNase I possess unique features that widen the capacity for genome-wide interrogation of chromatin accessibility.
Identification of nuclease accessible sites is a powerful approach to annotate regulatory regions of the genome. In cell lines and fresh tissue where nuclei can be isolated with high efficiency, DNase I has been used for decades to probe for accessibility. However DNase I is inhibited by actin and DNase I can therefore not be used on frozen tissue samples where nuclei cannot be purified in sufficient amounts while maintaining nuclear integrity. Here, we show the use of two other nucleases, Benzonase and Cyanase, in probing chromatin accessibility in frozen whole tissue samples. The digestion patterns of chromatin with Benzonase and Cyanase are remarkably similar and the identified accessible regions correlate with low nucleosomal occupancy, epigenetic marks of euchromatin and levels of transcription of proximal genes, validating the use of TACh for identification of regulatory elements in the genome. Comparisons with DNase I accessibility in liver tissue display a significant overlap between DNase I accessible sites and accessible regions identified by TACh. The similarities are most pronounced at highly accessible regions, whereas identification of less accessible regions tends to be more divergent between nucleases. Interestingly, we show that some of the differences between DNase I and Benzonase relate to their intrinsic sequence biases, with Benzonase and Cyanase preferring GC-rich sequences and disfavoring AT-rich sequences.
Currently, formaldehyde-assisted isolation of regulatory elements (FAIRE) is the only reported method to identify regulatory elements in frozen tissue samples ; however, FAIRE is an ambiguous approach that relies on the density of proteins at regions of crosslinked chromatin, where less dense regions are identified by FAIRE . As a consequence, FAIRE signals depend greatly on crosslink efficiencies of the tissue material, sonication efficiency and organic extraction steps. FAIRE shows some degree of overlap with DNase I accessible regions, however, certain regions are FAIRE specific and many regions are not identified by FAIRE . Since the TACh procedure like the DNase I assay relies on nuclease accessibility of DNA within chromatin, TACh is a more reliable procedure that is less dependent on extrinsic variables.
For clinical samples, TACh provides several advantages compared to current nuclease based methods. 1) TACh uses frozen tissues, which eliminates the need for immediate processing of freshly acquired tissue and is, therefore, compatible with the acquisition of samples under routine clinical conditions. 2) The ability to store whole tissue, fragments or pulverized powders provides flexibility and the powder ensures easy generation of matched aliquots that can be used for additional genomic experiments. 3) The ability to interrupt the procedure at multiple steps makes the whole process more amenable to clinical situations and laboratories. 4) The quick processing and the use of whole cells are likely to minimize nuclear damage and the loss of dissociated chromatin components. 5) Because the chromatin structure is preserved, it is likely that the same pulverized tissue can serve as the starting material for ChIP-Seq interrogations of the same tissue.
Thus TACh will make it possible to mine the accessible genomes from banks of frozen specimens, including biopsies and resected material from diseased and healthy human tissue, and is likely to assist in understanding the pathophysiology of a number of disease states. In particular, since TACh efficiently identifies accessible CpG islands, the method may be used to characterize chromatin structure rearrangements during progression of cancer, which are often associated with abnormal DNA methylation at CpG island rich promoters, leading to deregulation of numerous genes . Identification of genomic regions that specifically change accessibility during tumorigenesis may have significant prognostic value.
The described TACh methodology is a robust method for highly sensitive and comprehensive detection of accessible chromatin in samples derived from frozen tissue. The robustness and quick processing time of the assay provides feasible analysis of multiple tissue biopsies and we propose that application of TACh on clinically derived tissue material will provide knowledge on changes in chromatin accessability during progression of diseases.
HL-60 (ATCC) suspension cells were cultured in RPMI + 10% FBS in CO2 at 37°C.
Male C57BL/6 mice, 8-week old purchased from Harlam Sprague Dawley (Frederick, MD), were maintained according to the NIH guidelines with a 14-h light, 10-h dark cycle and free access to food and water. Mice were killed by decapitation or cervical dislocation and livers were immediately sectioned into 5-10 mm cubes and frozen in liquid nitrogen or processed for nuclei purification. All animal procedures were approved by the Animal Users and Care Committee, NICHD and NCI, NIH.
TACh on whole cells
Cells were collected by centrifugation, washed twice in ice cold cellular wash buffer (CWB) (20 mM TrisHCl pH 7.5, 137 mM NaCl, 1 mM EDTA , 10 mM sodium butyrate, 10 mM sodium orthovanadate, 2 mM sodium fluoride, protease inhibitor cocktial (Roche)) and resuspended in (40 million cells/ml) hypertonic lysis buffer (HLB) (20 mM TrisHCl pH 7.5, 2 mM EDTA, 1 mM EGTA, 0.5% glycerol, 20 mM sodium butyrate, 2 mM sodium orthovanadate, 4 mM sodium fluoride, protease inhibitor cocktial (Roche)). Cells were distributed in 500ul aliquots in 1.5 ml tubes and following addition of 500ul of nuclease digestion (ND) buffer (40 mM TrisHCl pH 8.0, 6 mM MgCl2, 0.3% NP-40, 1% Glycerol) containing a 3-fold dilutions (from 0.125 units/ml to 6U units/ml) of Benzonase (Calbiochem/EMD) or Cyanase (RiboSolutions, Inc). This was mixed gently and incubated for 3 minutes at 37°C. Reactions were terminated by the addition of EDTA (10 mM final) and SDS (0.75% final). Proteinase K was added to a final concentration of 0.5 mg/ml and incubated overnight at 45°C. DNA was purified by Phenol/Chloroform/isoamyl (PCI) extraction, ethanol precipitated and processed for Southern blotting.
Purified nuclease digested DNA was digested with HindIII overnight and then PCI extracted and ethanol precipitated. Precipitated DNA was SacI digested, RNase treated, extracted by PCI and ethanol precipitated. Digested DNA was separated on an agarose gel and transferred to a nylon membrane and incubated with a biotinylated c-myc probe  overnight at 50°C. Membranes were washed, incubated with HRP-conjugated streptavidin solution (kit from KPL, Inc.) and developed using chemiluminescent reagents (ECL Plus from Amersham/GE HealthCare).
TACh on frozen liver tissue
Livers were rapidly removed from mice and frozen immediately in liquid nitrogen. Frozen liver fragments were pulverized using a stainless steel pulverizer (Bessman), pre-chilled in liquid nitrogen. Seven hundred mg of frozen (the procedure may be done on a smaller amount of starting material), pulverized tissue was transferred to a pre-cooled 15 ml tube and suspended in 4 ml HLB (30 mM TrisHCl pH8.0, 2 mM EDTA, 2 mM EGTA, 20 mM sodium butyrate, 2 mM sodium orthovanadate, 4 mM sodium fluoride, protease inhibitor cocktail (Roche)). The suspension was mixed with a 3 cc syringe fitted with a 19 gauge needle and HLB was added to a final volume of 7 ml. The resulting suspension was passed through a 19 gauge needle 5 times, followed by a 22 gauge needle 5 times and finally a 23 gauge 10 times, using a 3 cc syringe. The suspension was distributed in 500ul aliquots in 1.5 ml tubes using a 23 gauge needle and incubated on ice for 5 minutes. 500ul of ND buffer (30 mM TrisHCl pH 8.0, 14 mM MgCl2 0.5% NP-40, 0.2% fatty acid free BSA) containing a 2-fold dilutions of Benzonase (Calbiochem/EMD) or Cyanase (RiboSolutions, Inc) was mixed gently with the 500ul tissue solution and incubated for 3 minutes at 37°C. Reactions were terminated by the addition of a final concentration of 50 mM EDTA and 0.1% SDS. 50ul RNaseA/RNaseT1 (Ambion) was added and the reaction was incubated over night at 40°C. SDS was brought to a final concentration of 0.75% and incubated for 2 hours at 45°C, after which Proteinase K was added (0.8ug/ml final) and incubated overnight at 45°C. DNA fragments of 100-500 bp from a chromatin digestion were purified over sucrose gradients  and precipitated in 0.1 volume NaAc and 0.7 volume isopropanol.
DNase I digestion of chromatin from liver tissue
Livers was removed from mice and homogenized in 8 ml/g tissue of low sucrose buffer (250 mM sucrose, 15 mM Tris–HCl pH 8.0, 15 mM NaCl, 60 mM KCl, 1 mM EDTA, 0.5 mM EGTA, 1 mM spermidine, protease inhibitors (Roche)) using a Type B dounce. Crude nuclear pellets were washed once in low sucrose buffer and resuspended in 9 volumes of high sucrose buffer (2 M sucrose, 15 mM Tris–HCl pH 8.0, 15 mM NaCl, 60 mM KCl, 1 mM EDTA, 0.5 mM EGTA, 1 mM spermidine, protease inhibitor cocktail (Roche)). The nuclear suspension was aliquoted in 2.0 ml tubes and centrifuged at 16,000xg for 30 minutes at 4°C. Nuclear pellets were combined in buffer (15 mM Tris–HCl pH 8.0, 15 mM NaCl, 60 mM KCl, 1 mM EDTA, 0.5 mM EGTA, 1 mM Spermidine, protease inhibitor cocktail (Roche)) and DNase I digestions were performed with 20 million nuclei as previously described (. DNA fragments of 100-500 bp from a chromatin digestion with 60U/ml DNase I (Sigma) were purified using sucrose gradients  and precipitated in 0.1 volume NaAc and 0.7 volume isopropanol.
Sequencing and data analysis
DNA was sequenced using Illumina GII sequencer at the Advanced Technology Center, NCI (Rockville, Md). Sequenced DNA was aligned to the mouse genome (mm9) using Eland or Bowtie . Hotspots were identified as previously described [1, 4, 28] using a Fdr of 0% and a tag density threshold at the mode. Hotspots were identified either from tag libraries generated from individual concentrations of Benzonase, Cyanase or DNase I or tag libraries pooled from different concentrations of enzymes or replicates of the same concentration. Bioinformatic analysis was performed using HOMER , Galaxy  and tools described in [4, 31]. Sequence data is accessible at GEO: GSE39982.
Previous published data used for analysis
RNA-seq ; SREBP-2 ChIP-seq peaks ; FOXA2 ChIP-seq peaks ; MNase-seq tag library (GEO: GSM717558) ; H3K4me1 and H3K4me3 ChIP-seq tag library (SRA008281) ; H3K27me2 ChIP-seq tag library (GEO: GSM751034).
ChIP combined with massive parallel sequencing
DNase I digestion of chromatin combined with massive parallel sequencing
Formaldehyde assisted identification of regulatory elements
Forkhead box A2
Histone3 lysin 4 mono-methylation
Histone3 lysine 4 tri-methylation
Histone3 lysine 27 tri-methylation
Massive parallel sequencing of RNA
Sterol Regulatory Element-Binding Protein
Tissue accessible chromatin
Transcriptional start site.
John S, Sabo PJ, Thurman RE, Sung MH, Biddie SC, Johnson TA: Chromatin accessibility pre-determines glucocorticoid receptor binding patterns. Nat Genet. 2011, 43: 264-268. 10.1038/ng.759.
Valouev A, Johnson SM, Boyd SD, Smith CL, Fire AZ, Sidow A: Determinants of nucleosome organization in primary human cells. Nature. 2011, 474: 516-520. 10.1038/nature10002.
Boyle AP, Davis S, Shulha HP, Meltzer P, Margulies EH, Weng Z: High-resolution mapping and characterization of open chromatin across the genome. Cell. 2008, 132: 311-322. 10.1016/j.cell.2007.12.014.
Siersbaek R, Nielsen R, John S, Sung MH, Baek S, Loft A: Adipogenic development is associated with extensive early remodeling of the chromatin landscape and establishment of transcription factor 'hotspots'. EMBO J. 2011, 30: 1459-1472. 10.1038/emboj.2011.65.
Ernst J, Kheradpour P, Mikkelsen TS, Shoresh N, Ward LD, Epstein CB: Mapping and analysis of chromatin state dynamics in nine human cell types. Nature. 2011, 473: 43-49. 10.1038/nature09906.
Song L, Crawford GE: DNase-seq: a high-resolution technique for mapping active gene regulatory elements across the genome from mammalian cells. Cold Spring Harb Protoc 2010. 2010
Wu C, Bingham PM, Livak KJ, Holmgren R, Elgin SC: The chromatin structure of specific genes: I. Evidence for higher order domains of defined DNA sequence. Cell. 1979, 16: 797-806. 10.1016/0092-8674(79)90095-3.
Wu C: The 5' ends of Drosophila heat shock genes in chromatin are hypersensitive to DNase I. Nature. 1980, 286 (5776): 854-860. 10.1038/286854a0.
Kaplan T, Li XY, Sabo PJ, Thomas S, Stamatoyannopoulos JA, Biggin MD: Quantitative Models of the Mechanisms That Control Genome-Wide Patterns of Transcription Factor Binding during Early Drosophila Development. PLoS Genet. 2011, 7: e1001290-10.1371/journal.pgen.1001290.
Lazarides E, Lindberg U: Actin is the naturally occurring inhibitor of deoxyribonuclease I. Proc Natl Acad Sci U S A. 1974, 71: 4742-4746. 10.1073/pnas.71.12.4742.
Stewart AF, Reik A, Schutz G: A simpler and better method to cleave chromatin with DNase 1 for hypersensitive site analyses. Nucleic Acids Res. 1991, 19: 3157-10.1093/nar/19.11.3157.
Benedik MJ, Strych U: Serratia marcescens and its extracellular nuclease. FEMS Microbiol Lett. 1998, 165: 1-13. 10.1111/j.1574-6968.1998.tb13120.x.
Meiss G, Friedhoff P, Hahn M, Gimadutdinow O, Pingoud A: Sequence preferences in cleavage of dsDNA and ssDNA by the extracellular Serratia marcescens endonuclease. Biochemistry. 1995, 34: 11979-11988. 10.1021/bi00037a040.
Michelotti GA, Michelotti EF, Pullner A, Duncan RC, Eick D, Levens D: Multiple single-stranded cis elements are associated with activated chromatin of the human c-myc gene in vivo. Mol Cell Biol. 1996, 16: 2656-2669.
Hoffman BG, Robertson G, Zavaglia B, Beach M, Cullum R, Lee S: Locus co-occupancy, nucleosome positioning, and H3K4me1 regulate the functionality of FOXA2-, HNF4A-, and PDX1-bound loci in islets and liver. Genome Res. 2010, 20: 1037-1051. 10.1101/gr.104356.109.
Li Z, Schug J, Tuteja G, White P, Kaestner KH: The nucleosome map of the mammalian liver. Nat Struct Mol Biol. 2011, 18: 742-746. 10.1038/nsmb.2060.
Deal RB, Henikoff JG, Henikoff S: Genome-wide kinetics of nucleosome turnover determined by metabolic labeling of histones. Science. 2010, 328: 1161-1164. 10.1126/science.1186777.
Mortazavi A, Williams B, Mccue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.
Dingwall C, Lomonossoff GP, Laskey RA: High sequence specificity of micrococcal nuclease. Nucleic Acids Res. 1981, 9: 2659-2673. 10.1093/nar/9.12.2659.
Biddie SC, John S, Sabo PJ, Thurman RE, Johnson TA, Schiltz RL: Transcription factor AP1 potentiates chromatin accessibility and glucocorticoid receptor binding. Mol Cell. 2011, 43: 145-155. 10.1016/j.molcel.2011.06.016.
Lupien M, Eeckhoute J, Meyer CA, Wang Q, Zhang Y, Li W: FoxA1 translates epigenetic signatures into enhancer-driven lineage-specific transcription. Cell. 2008, 132: 958-970. 10.1016/j.cell.2008.01.018.
Seo YK, Jeon TI, Chong HK, Biesinger J, Xie X, Osborne TF: Genome-wide localization of SREBP-2 in hepatic chromatin predicts a role in autophagy. Cell Metab. 2011, 13: 367-375. 10.1016/j.cmet.2011.03.005.
Gaulton K, Nammo T, Pasquali L, Simon J, Giresi PG, Fogarty M: A map of open chromatin in human pancreatic islets. Nat Genet. 2010, 42: 255-259. 10.1038/ng.530.
Giresi PG, Kim J, McDaniell RM, Iyer VR, Lieb JD: FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements) isolates active regulatory elements from human chromatin. Genome Res. 2007, 17: 877-885. 10.1101/gr.5533506.
Song L, Zhang Z, Grasfeder LL, Boyle AP, Giresi PG, Lee BK: Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity. Genome Res. 2011, 21: 1757-67. 10.1101/gr.121541.111.
Baylin SB, Jones PA: A decade of exploring the cancer epigenome - biological and translational implications. Nat Rev Cancer. 2011, 11: 726-734. 10.1038/nrc3130.
Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10: R25-10.1186/gb-2009-10-3-r25.
Sabo PJ, Kuehn MS, Thurman R, Johnson BE, Johnson EM, Cao H: Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays. Nat Methods. 2006, 3: 511-518. 10.1038/nmeth890.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P: Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010, 38: 576-589. 10.1016/j.molcel.2010.05.004.
Goecks J, Nekrutenko A, Taylor J: Galaxy: a comprehensive approach for supporting accessible, reproducible, and transparent computational research in the life sciences. Genome Biol. 2010, 11: R86-10.1186/gb-2010-11-8-r86.
Wiench M, John S, Baek S, Johnson TA, Sung MH, Escobar T: DNA methylation status predicts cell-type-specific enhancer activity. EMBO J. 2011, 30: 3028-3039. 10.1038/emboj.2011.210.
Wederell ED, Bilenky M, Cullum R, Thiessen N, Dagpinar M, Delaney A: Global analysis of in vivo Foxa2-binding sites in mouse adult liver using massively parallel sequencing. Nucleic Acids Res. 2008, 36: 4549-4564. 10.1093/nar/gkn382.
LG was supported by a research grant from the Lundbeck Foundation. This work was supported by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.
The authors declare that no competing interests exist.
LG carried out the DNase-seq in liver, participated in TACh-seq in liver, analyzed sequencing data and drafted the manuscript. RB developed the TACh method and performed TACh in cell lines and liver tissue, SJ participated in DNase-seq and TACh-seq in liver and helped to draft the manuscript. SB analyzed sequencing data. HJC helped develop the TACh method. YL and GA performed all the mouse work. CO and GLH participated in the study design and GLH helped to the draft of the manuscript. DL designed and coordinated the study and helped draft the manuscript. All authors read and approved the manuscript.
Lars Grøntved, Russell Bandle, Sam John contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1. Table of aligned sequence tags. Figure S2: Correlation of sequence tag densities from chromatin digested with different concentration of Cyananse. Figure S3: Correlation of sequence tag densities from chromatin digested with Benzonase and Cyanase. Figure S4: Distribution of AT percentage surrounding DNase I, Benzonase and Cyanase hotspots. (PDF 1 MB)
Authors’ original submitted files for images
About this article
Cite this article
Grøntved, L., Bandle, R., John, S. et al. Rapid genome-scale mapping of chromatin accessibility in tissue. Epigenetics & Chromatin 5, 10 (2012). https://doi.org/10.1186/1756-8935-5-10