PMC:1885650 JSONTXT 6 Projects

Real-time PCR mapping of DNaseI-hypersensitive sites using a novel ligation-mediated amplification technique Abstract Mapping sites within the genome that are hypersensitive to digestion with DNaseI is an important method for identifying DNA elements that regulate transcription. The standard approach to locating these DNaseI-hypersensitive sites (DHSs) has been to use Southern blotting techniques, although we, and others, have recently published alternative methods using a range of technologies including high-throughput sequencing and genomic array tiling paths. In this article, we describe a novel protocol to use real-time PCR to map DHS. Advantages of the technique reported here include the small cell numbers required for each analysis, rapid, relatively low-cost experiments with minimal need for specialist equipment. Presented examples include comparative DHS mapping of known TAL1/SCL regulatory elements between human embryonic stem cells and K562 cells. INTRODUCTION Mapping the location of DNaseI-hypersensitive sites (DHSs) remains central to developing our understanding of transcriptional regulation. Elements with a range of transcriptional regulatory functions have been identified initially as DHSs. These include transcriptional enhancers (1,2) and repressors (3,4) as well as chromatin insulators and barrier elements (5,6). A number of techniques have been published recently that permit the mapping of DHSs without the need for Southern blotting (7–13). These include high-throughput sequencing of cloned DNA libraries derived from DNaseI-digested chromatin (8,9), and a number of different approaches that use genomic tiling path arrays to map DHSs (7,11,13). While these approaches have the advantage of covering large genomic regions with a limited number of experiments, they are inherently costly and less applicable to the rapid DHS mapping of specific genomic sites in a range of cell types. An alternative approach that does permit the targeted semi-quantitative DHS mapping of specific loci has used quantitative real-time PCR to map sites from digested DNA (10). This technique depends on the quantification of relative loss of PCR signal observed when PCR primers amplify across regions of digested DNA compared with amplification of undigested DNA. This approach has been reported as showing good sensitivity, but a limitation is the large number of PCR reactions required to quantify the calculated loss of signal with significant certainty. In this article, we present an alternative method to identify DHS using real-time PCR by adapting the protocol we have used to map DHS using genomic tiling path arrays (11). We have previously demonstrated the high sensitivity and specificity of the basic protocol with regard to identifying DHS across large genomic regions (11). Here, we detail the laboratory protocol that permits the rapid comparative mapping of known and candidate DHS between different cell types using real-time PCR. As examples, we include comparative DHS mapping of regulatory elements located across the extended TAL1 (T-cell acute lymphocytic leukaemia-1, also known as SCL (stem cell leukaemia)) locus in human embryonic stem (hES) cells and the leukaemia cell line K562. RESULTS Generating a library of DHSs The basic protocol is outlined in Figure 1. Nuclei are extracted from living cells, then digested on ice for 1 h with a range of DNaseI concentrations, as detailed in the ‘Materials and methods’ sections. Following RNaseA and proteinaseK treatment, the DNA is extracted and run on a 1% agarose gel to check for the size of digested DNA. The gel in Figure 1A shows the DNA from mouse thymocyte nuclei digested with 0, 40 and 120 units of DNaseI. Maximal enrichment at DHS is usually observed with samples that are not over-digested. In the experiment shown, maximal enrichment at a known control DHS was seen with 40 units DNaseI (real-time PCR quantification shown in Figure 2B). The DNA is then blunt-ended using T4 polymerase (Figure 1B) and ligated with an asymmetric double-stranded linker. After extraction, DNA is amplified using a biotinylated linker-specific primer and 35 thermal polymerase cycles, as detailed in ‘Materials and methods’ section. As the linker will ligate to both ends of digested DNA, the amplification will represent a mix of primer extension and PCR, depending on the length of DNA amplified. As has been previously reasoned (13), one double strand of DNA in a region of DNaseI hypersensitivity is more likely to be digested twice within a short distance than non-hypersensitive DNA. This will lead to the preferential amplification of DNA from regions of DHS. The amplified DNA is then extracted using para-magnetic streptavidin beads, which provides a DNA library representative of whole-genome DHS. Agarose gel electrophoresis of the DNA recovered from the beads confirms that the vast majority of these products are between 300 and 500 base pairs in length (Figure 1B). Figure 1. (A) Nuclei are isolated from cells, and aliquots are digested with a range of DNaseI as detailed in ‘Materials and methods’ section. The DNA is extracted and run on 1% agarose gels using gel electrophoresis (right panel). With the example shown, the 40-units sample gave maximal enrichment at a housekeeping promoter using the complete protocol. (B) The digested DNA is blunt-ended with T4 polymerase and ligated to the LP21–25 linker as detailed in ‘Materials and methods’ section. Following amplification with the biotinylated LP25 primer, the extracted DNA template represents a library of whole-genome DHS. When samples of this library are visualised using gel electrophoresis (bottom panel), the majority of products are between 300 and 500 bp in size. Figure 2. (A) The left-hand panels show the Sybr-green real-time PCR profiles (fluorescence versus cycle number) of 5-fold dilutions of quantified mouse genomic DNA standards for two primer pairs within (ppA) and 3′ (ppB) of the Stil promoter. The right-hand panel shows that DNA template derived from the 40 units DNaseI-treated sample amplified nearly six PCR cycles in advance of the 0 DNaseI-treated sample with ppA. No difference in amplification kinetics was seen between the 0 and 40 samples with ppB. (B) Quantification of samples relative to genomic standards using primers for the ‘housekeeping’ Hmbs promoter shows maximal enrichment with 40 units DNaseI treatment for primary mouse thymocytes (left panel) compared with maximal enrichment with 120 units DNaseI for the mouse T-cell line, BW5147 (right panel). Real-time PCR quantification of DNA within the DHS library Figure 2 documents the quantification of specific DNA sites from different mouse-cell-derived DNA libraries. Figure 2A shows the quantification of material from two primer sets, which are within (primer pair A) and 3′ (primer pair B) to the mouse Stil promoter. The left-hand panels of Figure 2A show the Sybr-green real-time fluorescence profiles using serial 5-fold dilutions of quantified mouse genomic DNA standards with primers A and B. A calculated standard curve then permits the quantification of DNA from this sequence in library samples. The right-hand panels of Figure 2A show the quantification of samples from 0 and 40 units DNaseI-treated material using primers A and B. While there is no difference in amplification between the samples using primer pair B, with primer pair A, the 40 units DNaseI sample amplifies ∼5.5 PCR cycles before the 0 units sample. This equates to greater than 40-fold enrichment at the Stil promoter compared with no observed enrichment downstream of the promoter. Figure 2B shows the quantifiable differences in enrichment at the porphobilinogen deaminase (Hmbs) promoter between primary mouse thymocytes and the mouse T-cell line BW5147. This representative primary mouse cell experiment was performed with the DNA shown in Figure 1A, which confirms that with these cells, maximal enrichment is observed from DNA that is not over-digested. The comparison of the thymocytes with the BW5147 cells illustrates another common finding that, in our experience of primary cell experiments, maximal enrichment is often obtained using lower amounts of DNaseI compared with cell lines. Comparison of real-time PCR DHS with Southern blotting As the protocol generates template by primer extending away from a DNaseI-digestion site, as well as PCR amplification between DNaseI-digestion sites, there is the potential for the genomic size of the real-time DHS to be larger than sites revealed by Southern blotting. There is also potential for signal to be lost at the most ‘open’ stretch of DNA, due to over-digestion with DNaseI. These issues are addressed in Figure 3. The upper panel shows a Southern blot of the human STIL (SCL/TAL1 interrupting locus) promoter in K562 cells. The BglII restricted fragment is probed from the 3′ end, as shown in the upper panel. This reveals a central DHS ∼500 bp wide, with a suggestion of weaker hypersensitivity for a few hundred base pairs 5′ and 3′ to the central region. The lower panel shows real-time PCR data from K562 material using 10 primer sets, each generating an amplicon ∼120 bp long. This permits the quantification of enrichments over 1200 bp, centred around the STIL transcription start site. The lower panel represents the mean ± SD enrichment from three independent biological replicates from K562 cells. The black bar denotes the location of the ‘core’ hypersensitive site as defined by Southern blotting. There is good correlation between the location of the DHS between the two techniques. The 5′ extension of enrichment seen in the lower panel appears to reflect the weaker 5′ extension seen in the Southern blot in the upper panel. A dip in enrichment is observed in the lower panel over the previously mapped transcription-factor-binding sites (14), which may represent over-digestion of the most accessible core region of the DHS. Figure 3. The upper panel shows a Southern blot DHS map of the human STIL promoter in K562 cells using a 3′ probe (grey box labelled ‘P’) to analyse the BglII-restricted fragment (cut −4 kb and +4.2 kb relative to the STIL transcription start site). This revealed a core DHS ∼500 bp wide. The lower panel represents the mean ± SD real-time PCR quantification of enrichments from three biological replicates of K562 material relative to known genomic standards, using 10 separate primer sets that together amplify ∼1200 bp centred around the STIL transcription start site. The black bar represents the approximate location of the DHS, as identified by the Southern blot experiment in the upper panel. The real-time PCR amplification profile closely represents the Southern blot profile, although a central dip in real-time PCR signal is observed over the site of transcription factor binding. Quantification of DHS across the TAL1 locus in hES and K562 cells As examples of relative enrichments at different regulatory elements, we present a comparison of enrichments across the extended TAL1 locus from human embryonic stem (hES) cells and K562 cells in Figure 4B. Relative quantification of mRNA expression using real-time PCR shows that both of these cell types express STIL, but only K562 cells express TAL1 (Figure 4A). TAL1 expression is critical to the establishment of haematopoiesis in a developing embryo (15). Post embryogenesis, TAL1 expression is maintained in the non-lymphoid haematopoietic system, although in addition, expression is observed in a range of non-haematopoietic tissues, including endothelium and brain (16–18). Over the past 10 years, we and others have dissected the regulatory elements that direct the expression of TAL1 to different tissues (11,19–31). These include the 3′ haematopoietic stem cell enhancer (+19) (24,30), the 5′ endothelial-haematopoietic enhancer (−4) (31), the 3′ erythroid enhancer (+40 in mouse, +50 in human) (28), and a number of neuronal elements (25,26). The location of the STIL promoter has been previously mapped (14), although no other STIL regulatory elements have yet been identified. Figure 4. (A) Quantitative reverse transcriptase PCR quantification of STIL and TAL1 expression in hES and K562 cells. Expression was normalised to β-actin and plotted as detailed in ‘Materials and methods’ section. Both cell types express STIL, but TAL1 transcripts were not detectable (U) in hES cells. (B) Quantitative real-time PCR DNaseI profiles from hES and K562 cells using nine primer sets across the TAL1 locus. Primer set numbers refer to approximate kilobases (kb) relative to the start of TAL1 exon 1a. Primer sets corresponding to known regulatory elements are highlighted in red. Profiles represent the mean ± SD enrichment of material derived from two independent hES cell biological replicates, and two independent K562 cell cultures. Each bar represents the mean ratio of enrichment obtained from DNaseI-treated samples versus DNaseI-untreated samples for each cell type. With both cell types, no enrichment is seen at the four control sites (−16, +5, +22, +70 kb). K562 cells show maximal enrichments at the TAL1 regulatory elements and the STIL promoter, while hES cells show the highest enrichments at the STIL promoter and +50 enhancer, with minimal enrichment at the other TAL1 regulatory elements. Figure 4B shows the mean ± SD ratio (DNaseI-treated enrichment/DNaseI-untreated enrichment) of quantitative enrichments from two independent hES cell biological replicates and two independent K562 biological replicates, using primer sets from different locations as indicated in the figure. Primers indicated in red correspond to known regulatory elements. The hES cell cultures were maintained as detailed in methods. Both cell types show significant enrichment at the STIL promoter. Four control regions were selected: −16 kb, upstream of TAL1; +5 and +22 kb, within the TAL1 locus; +70 kb, downstream of TAL1. These sequences were chosen as control regions as they show limited homology between a number of species (26) and were considered highly unlikely to represent regulatory elements. All four regions showed no enrichments between DNaseI-treated and -untreated samples in either the hES cells or K562 cells. The most striking differences between hES cells and K562 cells were found at the TAL1 regulatory elements. With K562 cells, there is prominent enrichment at the TAL1 promoter 1a, and the −4, +19 and +50 enhancers. These cells express high levels of TAL1 transcripts, and were originally derived from a patient with blastic transformation of chronic myeloid leukaemia. K562 cells are relatively undifferentiated, although they have a partial erythroid phenotype. This potentially explains the marked enrichment seen at the +19 stem cell enhancer and the +50 erythroid enhancer. This is similar to data obtained from primary mouse erythroid cells (day 14.5 foetal liver), which show the highest enrichments at the mouse +40 enhancer (homologous to the human +50) (data not shown). In contrast, the hES cells show markedly reduced enrichment at the TAL1 regulatory elements when compared with K562 cells. Although there is enrichment at the +50 enhancer, there is minimal enrichment at the −4, promoter 1A and +19 elements. Reduced/absent accessibility of regulatory elements to DNaseI digestion is consistent with the lack of expression of TAL1 in hES cells (Figure 4A). As ES cells differentiate to form embryoid bodies, expression of TAL1 is rapidly switched on, appearing from day 3 of mouse ES cell differentiation (32). The relative accessibility of the +50 enhancer may reflect a poised chromatin state in hES cells, which permits a rapid response to the changing transcription factor environment that accompanies differentiation. DISCUSSION The development of techniques that permit the rapid comparative mapping of DHS between different cell types will greatly facilitate the study of transcriptional regulation in both normal and diseased cells. Recently published high-throughput techniques that map DHS sites using high-throughput sequencing (8,9) and genomic array tiling paths (7,11,13) have clear advantages of scale over more targeted approaches. Major disadvantages of these approaches include cost and lack of focus. This makes them less suitable for many laboratories that want to assess the chromatin accessibility of a number of defined or presumed regulatory elements in a range of cell types. A real-time-PCR-based approach to DHS mapping has, therefore, a number of potential advantages for researchers interested in specific regulatory questions at defined loci. Real-time PCR is relatively inexpensive compared with the large-scale techniques, and permits a rapid, focused DHS analysis of defined regions of DNA from multiple cell types. It also provides flexibility, as any genomic region can be analysed from the DNA library derived from the DNaseI-treated samples by designing further real-time PCR primer sets. We have previously shown that our basic technique for amplifying a library of DHS generates a template representative of known DHS with excellent sensitivity and specificity (11). Although the experiments presented in this paper were each performed using 5 million cells/digestion condition, we have obtained reproducible data using 5-fold fewer cells as starting material. We feel our technique can deliver acceptable specificity and sensitivity for DHS mapping with small numbers of cells, and will therefore be of use to those researchers working with limited numbers of primary cells. An alternative approach using real-time PCR to define DHS has been previously published (10). The two approaches differ in that Dorscher et al. quantify the DHS through the loss of PCR signal obtained from DHS when DNA is digested, whereas our approach uses real-time PCR to quantify a gain of signal observed from DHS. Dorscher et al. report excellent sensitivity using their approach to map DHS. However, the technique depends on large numbers of comparative quantitative real-time PCR reactions across a region in both digested and undigested material, in order to quantify the loss of enrichment. One advantage of our technique published here is that data can be obtained using far fewer quantitative PCR reactions. The technique is highly reproducible, with relatively little variation in quantifiable enrichments observed between different biological replicates. Moreover, we demonstrate tissue specificity, with variable enrichment at known regulatory elements between different cell types. The technique published here permits the rapid comparative analysis of DHS between different cell types from relatively small numbers of cells. It will have potential use for researchers across a broad spectrum of biology for the study of transcriptional regulation in both healthy and diseased tissues. MATERIALS AND METHODS Cell culture Cell lines were maintained in culture as previously described (11). Care was taken to ensure maximum viability of cells when taken for experiments. The primary thymocytes used for Figure 2 were prepared by gentle physical disassociation of a whole thymus gland into PBS supplemented with 2% FCS. Cells were filtered to ensure that a single-cell suspension was taken forward for nuclei isolation. Human embryonic stem cells (hES) were grown in chemically defined media in the presence of Activin and FGF2, as detailed previously (33). In these conditions, hES cells remain homogenously undifferentiated. Nuclei preparation Up to 3 × 107 cells were washed in ice-cold PBS and resuspended in 2 ml of ice-cold cell lysis buffer [300 mM sucrose, 10 mM Tris pH 7.4, 15 mM NaCl, 5 mM MgCl, 0.1 mM EGTA, 60 mM KCl, 0.2% NP-40, 0.5 mM DTT, 0.5 μM spermidine, 1× protease inhibitor (complete, Roche)]. After 5 min, the lysed cells were spun at 500 g for 5 min at 4°C with brakes off. After careful removal of supernatant, the nuclei were gently resuspended in 200 μl of ice-cold reaction buffer (20 μl 10× DNaseI buffer, 4 μl glycerol, 176 μl water) using pipette tips with cut off ends. The nuclei were spun again at 500 g for 5 min at 4°C and, following supernatant removal, were resuspended in 30 μl of reaction buffer per DNaseI condition. For example, if 2 × 107 cells were taken to look at four different conditions (e.g. 0, 20, 60, 120 units DNaseI), the nuclei were resuspended in 120 μl of reaction buffer. Separate 30 μl aliquots were then taken and gently mixed with 70 μl of DNaseI mix (see Table 1) on ice. This made a final digestion volume of 100 μl for each sample, which was left to incubate for 1 h on ice in the cold room. Table 1. DNaseI mixes of 70 μl aliquots were prepared as detailed The resuspended nuclei were added as 30 μl aliquots to make a final digestion volume of 100 μl. After 1 h, 700 μl of nuclei lysis buffer (100 mM tris HCL pH 8, 5 mM EDTA pH 8, 200 mM NaCl, 0.2% SDS) was added to each sample with 50 μg proteinase K. Following gentle mixing with inversion, the lysed samples were incubated at 55°C for 1 h. RNaseA (10 μg (Ambion)) was then added to each sample and further incubated at 37°C for 30 min. DNA was then extracted using standard phenol–chloroform techniques. Care was taken to use cut-off tips and very gentle pipetting to reduce non-specific DNA sheering. Following precipitation, DNA was resuspended in 200 μl of 0.1 TE and quantified using spectophotometry. Samples (1 μg) were analysed using gel-electrophoresis, as shown in Figure 1A. DNA library preparation Following quantification, 7.5 μg of DNA was taken for each sample condition. This was blunt-ended using T4 polymerase (10 μl 10× buffer, 0.5 μl 25 mM dNTP, 3 μl BSA, 1 μl T4 polymerase (3 U/μl NE Biolabs), 85.5 μl of water/DNA in 0.1 TE). The samples were mixed gently on ice, then put at 12°C for 16 min. The reaction was stopped with excess EDTA (4 μl 0.5 M) followed by 75°C for 10 min. The DNA was further extracted using phenol–chloroform and precipitated with ethanol. Glycogen carrier (5 μg) was used at this stage, and the pellet resuspended in 80 μl water. The blunt-ended DNA was then ligated to the LP21–25 linker. A stock of linker was prepared by mixing 80 μl LP21 (100 μM: GAATTCAGATCTCCCGGGTCA) with 80 μl LP25 (100 μM: GCGGTGACCCGGGAGATCTG AATTC) and 240 μl water. The mix was then placed on a 95°C hot block and the power supply turned off. When the block had cooled to room temperature, the LP21–25 linker was aliquoted and frozen. The 80 μl DNA sample was split into two samples for ligation (5 μl 10× buffer, 3 μl LP21–25 linker, 3 μl ligase (NE biolabs 400 U/μl and 39 μl DNA). This was mixed well with a pipette and left at 16°C overnight. Following ligation, the DNA was precipitated with ethanol, and resuspended in 42 μl water. The sample was amplified using Vent exo-polymerase and a biotinylated LP25 primer. The mix (5 μl ThermoPol 10× buffer, 0.5 μl B-LP25 (200 μM), 2 μl Vent exo-(NE Biolabs 2 U/μl), 1 μl dNTP (25 mM) and 41.5 μl DNA) was amplified with the following thermal cycle: (95°C 3 min, 95°C 30 s, 61°C 30 s, 72°C 30 s) × 35 cycles. Following amplification, the biotinylated products were extracted using Dynal streptavidin beads. For each sample, 30 μl of beads (Dyabeads M-270; Dynal biotech) were washed twice in 2× binding buffer (10 mM Tris HCL pH7.5, 1 mM EDTA, 2 M NaCl) using a magnet, and then resuspended in 50 μl of 2× binding buffer/sample. Each post amplification sample (50 μl) was then mixed with 50 μl of resuspended beads and incubated at room temperature on a rotator for 1 h. After binding, the samples were washed twice with 2× binding buffer, then once with 1× TE. After the last wash, each sample was resuspended in 30 μl 0.1 TE. The paired samples that were split before ligation were then pooled, making a 60 μl final aliquot. This was heated at 95°C for 10 min to free the DNA from the beads, and then stored at 4°C. Real-time PCR quantification We have used Stratagene Brilliant SYBR Green QPCR Master Mix as our standard kit for quantification of genomic controls and samples. We have used a range of primer sets. The HMBS primers are: human—Forward; ATGCTGCCTATTTCAAGGTTGT, Reverse; GAATT GGAACATTGCGACAGT, and mouse—Forward; CGGAGTCATGTCCGGTAAC, Reverse; CGACCAA TAGACGACGAGAA. Primers from the TAL1 locus (human and mouse) are available on request. Amplification conditions and PCR mixes were as recommended by Stratagene. Genomic standard curves were calculated for each primer using serial 5-fold dilutions from 50 ng genomic DNA/PCR reaction. For each sample PCR, a 45 μl stock was made using 1.3 μl of sample DNA + 43.7 μl water. This stock was used as 5 μl per real-time PCR reaction, and following amplification, quantified relative to the genomic DNA standard curve, for each primer set. The quantifications shown in Figure 3 represent the absolute quantification of DNA at each primer set for the 5 μl PCR sample. The mRNA expression levels of TAL1 and STIL in Figure 4A were normalised relative to β-actin and plotted relative to mRNA expression from normal human donor CD34 + cells, as described previously (11).