- Methodology article
- Open Access
Rapid high-throughput analysis of DNaseI hypersensitive sites using a modified Multiplex Ligation-dependent Probe Amplification approach
BMC Genomics volume 10, Article number: 412 (2009)
Mapping DNaseI hypersensitive sites is commonly used to identify regulatory regions in the genome. However, currently available methods are either time consuming and laborious, expensive or require large numbers of cells. We aimed to develop a quick and straightforward method for the analysis of DNaseI hypersensitive sites that overcomes these problems.
We have developed a modified Multiplex Ligation-dependent Probe Amplification (MLPA) approach for the identification and analysis of genomic regulatory regions. The utility of this approach was demonstrated by simultaneously analysing 20 loci from the ENCODE project for DNaseI hypersensitivity in a range of different cell lines. We were able to obtain reproducible results with as little as 5 × 104 cells per DNaseI treatment. Our results broadly matched those previously reported by the ENCODE project, and both technical and biological replicates showed high correlations, indicating the sensitivity and reproducibility of this method.
This new method will considerably facilitate the identification and analysis of DNaseI hypersensitive sites. Due to the multiplexing potential of MLPA (up to 50 loci can be examined) it is possible to analyse dozens of DNaseI hypersensitive sites in a single reaction. Furthermore, the high sensitivity of MLPA means that fewer than 105 cells per DNaseI treatment can be used, allowing the discovery and analysis of tissue specific regulatory regions without the need for pooling. This method is quick and easy and results can be obtained within 48 hours after harvesting of cells or tissues. As no special equipment is required, this method can be applied by any laboratory interested in the analysis of DNaseI hypersensitive regions.
Open chromatin is a characteristic of genomic loci with regulatory functions. These regions are preferentially digested by DNaseI , and the identification of DNaseI hypersensitive sites is frequently used to identify and analyse regulatory regions such as promoters, enhancers and silencers [2, 3]. However, currently available methods have significant limitations. A commonly used approach involves Southern blotting, but this is time consuming, usually requires radioactivity and is limited to short stretches of DNA. Several PCR-based methods have been described [4, 5], but these do not readily allow multiplexing. Recent reports of large scale analysis of DNaseI hypersensitive sites have used either microarrays [6–9] or deep sequencing [9, 10]. Whilst valuable for genome wide analysis, the costs involved are a limiting factor for many applications (such as comparing different developmental stages or tissues). Another disadvantage of those methods is that they usually require many millions of cells. For ex vivo studies, this might require extensive pooling of tissues, meaning that these methods are not suitable for all applications.
Multiplex Ligation-dependent Probe Amplification (MLPA) was originally developed to detect deletions and duplications in genomic DNA , and has become popular in diagnostic settings for a range of disorders [12, 13]. It has since been modified for several other applications as well, including methylation analysis , mRNA expression analysis , identifying copy number variation in normal populations [16, 17], genotyping of mouse models [18, 19] and measuring the efficiency of Cre-mediated recombination in mouse models . The principle advantages of this method are the sensitivity and multiplexing potential. It can be used to analyse up to 50 genomic loci with as little as 20 ng genomic DNA in a single reaction. Furthermore, the only equipment that is required is a thermocycler and DNA sequencer, readily available to most researchers.
We describe here a quick and straightforward protocol for analysing DNaseI hypersensitive sites. This is demonstrated by the analysis of 20 different loci in a single reaction, based on data published by the ENCODE consortium .
Figure 1 outlines the protocol used in this study. To examine the general feasibility of our approach, we designed 11 probes to cover randomly chosen DNaseI hypersensitive sites in HeLa cells as published by the ENCODE consortium . In addition, nine probes were designed in regions that showed no evidence of DNaseI hypersensitivity. To be able to cover larger genomic regions and to give greater flexibility in probe design we employed the recently developed extension MLPA . All 20 probes were combined in a single mix, and could be differentiated from each other on the basis of length (the final product length range was 94-207 bp).
Nuclei isolation and DNaseI treatment
We tested different methods for the isolation of nuclei, with the Ne-Per Kit (Thermo Scientific) found to be most suitable. Although intended for the isolation and separation of nuclear and cytoplasmic proteins, it proved to be useful for the isolation of nuclei. Another obstacle to overcome was the isolation of the genomic DNA. This approach requires efficient recovery of very long as well as very short fragments of genomic DNA. Unfortunately, DNA obtained using the standard method for this task, phenol-chloroform purification, is not suitable for use with MLPA, as traces of phenol are known to interfere with the enzymatic reactions. Instead we used the HighPure PCR Purification Kit (Roche), which, according to the manufacturer, is suitable for the isolation of fragments of up to 50 kb. As expected, in our experiments the recovery of DNA from the untreated and therefore most intact sample, was often least efficient, however, the quality and quantity of the DNA was still sufficient for MLPA analysis. Figure 2 shows a typical example of digested DNA obtained after treating isolated nuclei with increasing amounts of DNaseI and purification.
Analysis of DNaseI hypersensitive sites in HeLa cells
We next used the same probe mix to analyse the DNaseI hypersensitivity of DNA within intact nuclei of HeLa cells. Figure 3A shows typical peak patterns obtained using nuclei aliquots digested with increasing amounts of DNaseI. Several peaks show significantly decreased peak height with increased DNaseI concentration, while others remain virtually unchanged. Nine of 11 probes targeting previously described DNaseI hypersensitive sites in HeLa cells show a clear decrease in normalized peak heights (defined as < 75% of the equivalent peak in undigested DNA), whereas all nine probes targeting non-sensitive loci show no significant decrease (figure 3B). These results were highly reproducible, with the technical and biological replicates giving an r2 > 0.9.
We also tested the same probe mix on DNA directly isolated from HeLa cells. As can be seen in Figure 3C, there are no significant changes in normalized peak heights with increasing DNaseI digestion. This was expected as these DNA samples should have no chromatin structure, and all cuts by DNaseI should therefore occur in a random fashion.
Comparison of results using other cell lines
We then analysed the same 20 loci for DNaseI hypersensitivity in other cell lines. For this comparison we used HEK293, C28 and S97 cells. While the results for C28 cells were very similar to those of HeLa cells (data not shown), we could identify differences in sensitivity for several probes in HEK293 and S97 cells. As shown in Figure 4, some probes identifying hypersensitive sites in HeLa cells stay unchanged in the other cell lines (probes S9 in HEK293 cells, S6 in S97 cells) and two of the probes not showing the expected sensitivity in HeLa cells (S3 and S4) do show sensitivity in HEK293 cells.
Determining the minimum cell number
We also tried to determine the minimum cell number that could be used with this approach. As we were able to obtain robust results and sufficient DNA for several replicates from aliquots containing nuclei from about 2.5 × 105 cells, we estimated the lower limit of required cells to be around 5 × 104 (theoretically ~300 ng DNA). To confirm this number, we used aliquots of 5 × 104 HeLa cells and treated them as described above and compared the results with those obtained from experiments using 2.5 × 105 cells per aliquot. Figure 5 shows the comparison of results from 5 × 104 and 2.5 × 105 cells (r2 = 0.97).
We describe here a simple technique that allows the rapid analysis of many DNaseI hypersensitive sites using little starting material. For a proof of principle we chose loci that had been analysed by both array analysis and deep sequencing as part of the ENCODE project .
We used 20 probes in one reaction (11 sensitive, 9 non-sensitive in HeLa cells according to ENCODE data). Nine of 11 of the probes targeted to hypersensitive sites showed the expected decrease in peak height for HeLa cells. The fact that two sensitive sites did not show the expected drop in peak heights is probably due to different growth conditions, differences in passage numbers or to our batch of HeLa cells being different from the one used in the ENCODE study rather than to limitations of our method. This is supported by the high correlation observed for both technical and biological replicates (r2 > 0.9) and the fact that those two probes showed sensitivity in HEK293 cells.
We also examined three other cell lines, in order to examine the cell line specificity of our experimental approach. We could show clear differences in hypersensitivity for some loci compared to HeLa cells. Again, the conclusion that these are cell type-specific differences is supported by the concordance of replicate experiments (r2 > 0.9).
To determine suitable DNaseI digestion conditions for our method, a range of different DNaseI concentrations and incubation conditions were tested. Incubating the digestion reactions with the indicated DNaseI amounts at room temperature for 20 minutes gave more reproducible results than, for example, using less DNaseI and incubating at 37°C or incubations on ice (data not shown). By analysing the digested DNA on an agarose gel it was possible to identify the most appropriate conditions. Most of our experiments were carried out using nuclei from 2.5 × 105 cells per DNaseI digestion aliquot, providing sufficient DNA for several technical replicates. However, in the case of HeLa cells this number could be reduced to only 5 × 104 cells, still providing highly reproducible results as shown by an r2 of 0.97 when comparing the results of the experiments with both cell numbers. Furthermore, we found that two different DNaseI concentrations and an undigested control are sufficient to analyse hypersensitive sites, which reduces the amount of total starting material required to about 1.5 × 105 cells.
Although we have used 20 probes in this study, there is the potential to significantly increase the degree of multiplexing. Reducing the interprobe spacing to 3 bp would allow up to 40 loci to be examined. Indeed, by designing probes that can be labelled with different fluorophores [22, 23] it would be possible to analyse > 100 loci in a single reaction.
We describe here a rapid and accurate method for assaying DNaseI hypersensitive sites. In contrast to genome-wide approaches such as deep sequencing or microarray analysis, we consider the primary strength of this approach to be when < 100 genomic loci are being analysed. As all loci are analysed in a single reaction, relatively little starting material is required. We have shown that < 105 cells per DNaseI treatment can be used, which will allow the study of e.g. embryonic organ development without extensive pooling. The protocol is straightforward and results can be easily obtained from many samples within 48 hours. Finally, using this method makes it possible to also detect heterozygous and homozygous deletions and duplications of the examined regions, which is an artifact that is known to occur in cultured cells.
HEK293 (human embryonic kidney cells), C28 (chondrogenic cells), HeLa (cervical cancer cells) and S97 (dermal fibroblasts) were grown in DMEM containing 10% FBS and supplemented with L-glutamine in T-25 flasks and incubated at 37°C containing 5% CO2 in a humidified atmosphere.
Cell harvest and isolation of nuclei
After reaching 95-100% of confluence, cells were washed with PBS and harvested using 1 ml of 0.025% Trypsin-EDTA. After 5-10 min incubation at 37°C, cells were washed with 2 ml of 10% FBS containing DMEM and carefully resuspended in a 1.5 ml Eppendorf tube in cold PBS and placed on ice. To limit changes in chromatin structure during treatment, nuclei isolation was performed using 2.5 × 105 to 106 cells as soon as possible.
For isolation of nuclei, the NE-PER Kit (Thermo Scientific) was used with the following modifications. To prevent nuclei from excessive clumping and releasing DNA, twice the recommended volumes of the solutions were used and resuspension of cells/nuclei was carried out by carefully pipetting up and down rather than vortexing. Centrifugation steps were carried out at 250 × g at 4°C.
DNaseI treatment of nuclei
Isolated nuclei were washed in 500 μl cold DNaseI buffer containing 2% glycerol and carefully resuspended in 75 μl DNaseI buffer containing 2% glycerol. 25 μl aliquots of the nuclei suspension were added to 2 ml Eppendorf tubes containing one control without DNaseI and increasing amounts (0.5 - 2 units) of DNaseI (Promega) in 50 μl DNaseI buffer containing 2% glycerol. The solutions were mixed by carefully flicking the tubes and incubated for 20 min at room temperature (23°C).
Isolation and purification of DNaseI treated genomic DNA
DNaseI treated nuclei were lysed by adding 250 μl of nuclei lysis buffer (200 mM NaCl, 150 mM Tris HCl pH 8, 10 mM EDTA pH 8, 0.2% SDS) containing 50 μg of Proteinase K and incubated for 45 min at 55°C. 20 μg of RNaseA (Sigma) was added and incubated for 30 min at 37°C.
The DNaseI treated DNA was isolated and purified using the HighPure PCR Purification Kit (Roche) according to the manufacturer's instructions. DNA was eluted with 50 μl elution buffer provided with the kit. Purified DNA was checked for quality and degree of DNaseI digestion by agarose gel electrophoresis. DNA concentrations were measured using the NanoDrop.
MLPA and fragment analysis
The MLPA probe mix was created by combining each MLPA oligonucleotide in table 1 to a final concentration of 4 nM. The extension MLPA was performed essentially as described . 100-200 ng of DNaseI-treated DNA of each sample was denatured at 98°C for 5 min, and allowed to cool to room temperature. 1.5 μl MLPA buffer (MRC-Holland, the Netherlands) and 1.5 μl MLPA probe mix were added to each sample, denatured at 95°C for 1 min then hybridized overnight at 60°C. The following morning a ligation and elongation step was carried out at 54°C by adding 3 μl ligase buffer A, 3 μl ligase buffer B, 1 μl Ligase-65 (all MRC Holland, the Netherlands), 2 mM dNTPs, 1 U Stoffel Taq polymerase (Applied Biosystems) and H2O to a final volume of 40 μl. After 15-20 min the ligase was inactivated by heating to 95°C for 5 min. The PCR reaction consisted of 5 μl of ligated mix added to 2 μl SALSA PCR buffer, 1 μl SALSA enzyme buffer, 1 μl SALSA PCR primer, 0.25 μl SALSA polymerase (all MRC-Holland, the Netherlands) and 15.75 μl H2O. The PCR reaction started with a 1 min denaturation step at 95°C, followed by 35-38 cycles of amplification (each cycle consisting of 95°C for 30 sec, 58°C for 30 sec, 72°C for 30 sec) and concluded with 20 min at 72°C.
Fragment analysis was performed on an ABI3130 capillary sequencer. Peak data was extracted using GeneMarker software (Soft Genetics) and exported to Excel (Microsoft) for further analysis.
MLPA data analysis
Basic data analysis was performed as described . Peak data for the naked DNA experiment was normalised to the average of all probes in the examined samples, with all other reactions being normalized to all non-sensitive probes. To assist in visualisation of the results, normalized ratios of all untreated samples were set to 1.
Gross DS, Garrard WT: Nuclease hypersensitive sites in chromatin. Annual review of biochemistry. 1988, 57: 159-197. 10.1146/annurev.bi.57.070188.001111.
Wu C, Wong YC, Elgin SC: The chromatin structure of specific genes: II. Disruption of chromatin structure during gene activity. Cell. 1979, 16 (4): 807-814. 10.1016/0092-8674(79)90096-5.
Burgess-Beusse B, Farrell C, Gaszner M, Litt M, Mutskov V, Recillas-Targa F, Simpson M, West A, Felsenfeld G: The insulation of genes from external enhancers and silencing chromatin. Proceedings of the National Academy of Sciences of the United States of America. 2002, 99 (Suppl 4): 16433-16437. 10.1073/pnas.162342499.
Follows GA, Janes ME, Vallier L, Green AR, Gottgens B: Real-time PCR mapping of DNaseI-hypersensitive sites using a novel ligation-mediated amplification technique. Nucleic acids research. 2007, 35 (8): e56-10.1093/nar/gkm108.
Dorschner MO, Hawrylycz M, Humbert R, Wallace JC, Shafer A, Kawamoto J, Mack J, Hall R, Goldy J, Sabo PJ, et al: High-throughput localization of functional elements by quantitative chromatin profiling. Nature methods. 2004, 1 (3): 219-225. 10.1038/nmeth721.
Crawford GE, Davis S, Scacheri PC, Renaud G, Halawi MJ, Erdos MR, Green R, Meltzer PS, Wolfsberg TG, Collins FS: DNase-chip: a high-resolution method to identify DNase I hypersensitive sites using tiled microarrays. Nature methods. 2006, 3 (7): 503-509. 10.1038/nmeth888.
Sabo PJ, Kuehn MS, Thurman R, Johnson BE, Johnson EM, Cao H, Yu M, Rosenzweig E, Goldy J, Haydock A, et al: Genome-scale mapping of DNase I sensitivity in vivo using tiling DNA microarrays. Nature methods. 2006, 3 (7): 511-518. 10.1038/nmeth890.
Follows GA, Dhami P, Gottgens B, Bruce AW, Campbell PJ, Dillon SC, Smith AM, Koch C, Donaldson IJ, Scott MA, et al: Identifying gene regulatory elements by genomic microarray mapping of DNaseI hypersensitive sites. Genome research. 2006, 16 (10): 1310-1319. 10.1101/gr.5373606.
Boyle AP, Davis S, Shulha HP, Meltzer P, Margulies EH, Weng Z, Furey TS, Crawford GE: High-resolution mapping and characterization of open chromatin across the genome. Cell. 2008, 132 (2): 311-322. 10.1016/j.cell.2007.12.014.
Crawford GE, Holt IE, Whittle J, Webb BD, Tai D, Davis S, Margulies EH, Chen Y, Bernat JA, Ginsburg D, et al: Genome-wide mapping of DNase hypersensitive sites using massively parallel signature sequencing (MPSS). Genome research. 2006, 16 (1): 123-131. 10.1101/gr.4074106.
Schouten JP, McElgunn CJ, Waaijer R, Zwijnenburg D, Diepvens F, Pals G: Relative quantification of 40 nucleic acid sequences by multiplex ligation-dependent probe amplification. Nucleic acids research. 2002, 30 (12): e57-10.1093/nar/gnf056.
Lalic T, Vossen RH, Coffa J, Schouten JP, Guc-Scekic M, Radivojevic D, Djurisic M, Breuning MH, White SJ, den Dunnen JT: Deletion and duplication screening in the DMD gene using MLPA. Eur J Hum Genet. 2005, 13 (11): 1231-1234. 10.1038/sj.ejhg.5201465.
Kozlowski P, Roberts P, Dabora S, Franz D, Bissler J, Northrup H, Au KS, Lazarus R, Domanska-Pakiela D, Kotulska K, et al: Identification of 54 large deletions/duplications in TSC1 and TSC2 using MLPA, and genotype-phenotype correlations. Human genetics. 2007, 121 (3-4): 389-400. 10.1007/s00439-006-0308-9.
Nygren AO, Ameziane N, Duarte HM, Vijzelaar RN, Waisfisz Q, Hess CJ, Schouten JP, Errami A: Methylation-specific MLPA (MS-MLPA): simultaneous detection of CpG methylation and copy number changes of up to 40 sequences. Nucleic acids research. 2005, 33 (14): e128-10.1093/nar/gni127.
Eldering E, Spek CA, Aberson HL, Grummels A, Derks IA, de Vos AF, McElgunn CJ, Schouten JP: Expression profiling via novel multiplex assay allows rapid assessment of gene regulation in defined signalling pathways. Nucleic acids research. 2003, 31 (23): e153-10.1093/nar/gng153.
Young JM, Endicott RM, Parghi SS, Walker M, Kidd JM, Trask BJ: Extensive copy-number variation of the human olfactory receptor gene family. American journal of human genetics. 2008, 83 (2): 228-242. 10.1016/j.ajhg.2008.07.005.
White SJ, Vissers LE, Geurts van Kessel A, de Menezes RX, Kalay E, Lehesjoki AE, Giordano PC, Vosse van de E, Breuning MH, Brunner HG, et al: Variation of CNV distribution in five different ethnic populations. Cytogenetic and genome research. 2007, 118 (1): 19-30. 10.1159/000106437.
Notini AJ, Li R, Western PS, Sinclair AH, White SJ: Rapid and reliable determination of transgene zygosity in mice by multiplex ligation-dependent probe amplification. Transgenic Res. 2009,
Hewitt CA, Carmichael CL, Wilkins EJ, Cannon PZ, Pritchard MA, Scott HS: Multiplex ligation-dependent probe amplification (MLPA) genotyping assay for mouse models of down syndrome. Front Biosci. 2007, 12: 3010-3016. 10.2741/2291.
Leonhard WN, Roelfsema JH, Lantinga-van Leeuwen IS, Breuning MH, Peters DJ: Quantification of Cre-mediated recombination by a novel strategy reveals a stable extra-chromosomal deletion-circle in mice. BMC biotechnology. 2008, 8: 18-10.1186/1472-6750-8-18.
Birney E, Stamatoyannopoulos JA, Dutta A, Guigo R, Gingeras TR, Margulies EH, Weng Z, Snyder M, Dermitzakis ET, Thurman RE, et al: Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature. 2007, 447 (7146): 799-816. 10.1038/nature05874.
Harteveld CL, Voskamp A, Phylipsen M, Akkermans N, den Dunnen JT, White SJ, Giordano PC: Nine unknown rearrangements in 16p13.3 and 11p15.4 causing alpha- and beta-thalassaemia characterised by high resolution multiplex ligation-dependent probe amplification. Journal of medical genetics. 2005, 42 (12): 922-931. 10.1136/jmg.2005.033597.
White SJ, Vink GR, Kriek M, Wuyts W, Schouten J, Bakker B, Breuning MH, den Dunnen JT: Two-color multiplex ligation-dependent probe amplification: detecting genomic rearrangements in hereditary multiple exostoses. Human mutation. 2004, 24 (1): 86-92. 10.1002/humu.20054.
White SJ, Breuning MH, den Dunnen JT: Detecting copy number changes in genomic DNA: MAPH and MLPA. Methods in cell biology. 2004, 75: 751-768. full_text.
We thank Dr. Craig Smith for providing HeLa and HEK293 cells, Dr. Peter Kannu for providing C28 cells and Dr. Laura Zamurs for providing S97 cells. This work was supported by the National Health and Medical Research Council [491293; 546478 S.W., 334314 A.S.].
TO and SJW conceived the experimental design. TO, SE, WNL, AHS, and SJW participated in acquisition, analysis, and interpretation of the data. TO and SJW drafted the manuscript. TO, SE, WNL, AHS, and SJW participated in critical revision and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Ohnesorg, T., Eggers, S., Leonhard, W.N. et al. Rapid high-throughput analysis of DNaseI hypersensitive sites using a modified Multiplex Ligation-dependent Probe Amplification approach. BMC Genomics 10, 412 (2009). https://doi.org/10.1186/1471-2164-10-412
- HeLa Cell
- Hypersensitive Site
- Encode Project
- Normalize Peak Height
- Encode Consortium