Bisulfite sequencing of small DNA/cell samples (PROT35)
Sascha Tierling, Sabine Reither, Jörn Walter
Saarland University, FR 8.3 Biosciences
Dept. of Genetics/Epigenetics
66123 Saarbrücken, Germany
Email feedback to:
IntroductionSequencing of sodium-bisulfite modified genomic DNA originally introduced by M. Frommer (Frommer et al., 1992) is a widely used “gold standard” method for DNA-methylation analysis. Since this method relies on a harsh chemical treatment of DNA it causes a lot of DNA damage and hence a dramatic loss of high quality DNA for PCR amplification and further analysis. In the meantime several commercial kits are available for this procedure which work reasonably well when starting with large amounts of DNA.
Here we describe a protocol for small numbers of cells and little DNA which requires some specific handling. The protocol is based on a strategy originally introduced by our lab (Hajkova et al., 2002) using agarose embedded DNA. This physical trapping helps to avoid DNA loss during the various incubation steps while maintaining a good bisulphite conversion rate. We will introduce two alternative procedures to perform bisulphite treatment of agarose embedded small DNA aliquots or cells and guide through some generally critical points in the bisulphite reaction and primer design.
We also include tips for the process of data processing after sequencing which is facilitated by a new and very useful software tool (BiQ Analyzer). This tool is freely available and allows rapid and reproducible processing and evaluation of bisulphite sequencing data. It generates standardized table output formats allowing direct database integration.
Two alternative protocols (A and B) for agarose embedded DNA of cell preparation
In the following we will describe two protocols for bisulphite treatment of small aliquots of mammalian cells/DNA. Both procedures only differ in the first step, i.e. the preparation of the DNA for bisulphite treatment. We routinely use protocol A for collected mammalian cells in the range >100 and B for very scarce material (such as zygotes/early mouse embryos), i.e. few to 100 cells.
Protocol A: Preparation of genomic DNA from samples with more than 100 cells per sample
- Mix the collected mammalian cells (>100) with 50 µl of solution A (25 mM EDTA, 75 mM NaCl), 50 µl of solution B (10 mM EDTA, 10 mM Tri-HCl, pH 8.0, 1% SDS), add 2 µl of Proteinase K (20 mg/ml) and 5 µg of yeast t-RNA (use 1.5 ml Eppendorf tubes)
- Incubate 1h at 37°C
- Add 50 µl of phenol and 50 µl of chloroform/isoamylalcohol (24:1)
- Mix the samples gently by inverting the tubes
- Centrifuge the samples at 13,000 rpm for 10-15 min
- Transfer the supernatant into fresh tubes and add 100 µl of chloroform/isoamylalcohol (24:1)
- Mix the samples gently by inverting the tubes
- Centrifuge 5 min at 13,000 rpm
- Transfer the supernatant into fresh tubes and add another 5 µg of t-RNA and 1/10 volume of 3M sodium acetate (pH5.5), then add 350 µl of 100% ethanol
- Precipitate the DNA at -20°C over night or longer
- Centrifuge the samples for 30 min at 13,000 rpm at 4°C
- Remove the ethanol and let the pellets dry at room temperature
- Dissolve the pellet in 4 µl ddH2O over night (do not dissolve the pellets in larger volumes since the DNA is later embedded into one single agarose bead – highest yield; dissolving the pellet over night is necessary to bring the DNA entirely into solution)
- Add 0.4 µl of freshly prepared 2 M NaOH and 10 µl of molten 2% LMP agarose (comment 1)
- Overlay with 200-500 µl of heavy mineral oil (white, Aldrich)
- Boil for 10 min in a water bath
- Cool down on ice until the bead is solid
- Proceed with bisulphite treatment
Protocol B: Bisulfite treatment in agarose embedded cells
This Protocol should be used as an alternative if you want to avoid precipitation steps, e.g. when working with very low amounts of cells/DNA; few to 100 cells or equivalent DNA content (up to 300 pg).
- Add 15 µl of 2% LMP agarose to the sample and overlay with 100-300 µl of heavy mineral oil (note: total volume of the Cell suspension should not exceed 5µl)
- Boil the samples for 10 min in a water bath to destroy the cells
- Cool the samples for at least 10 min on ice until the bead is solid
- Remove the oil (be careful in case if the bead is still not solid)
- Add 400 µl of solution A (25mM EDTA, 75 mM NaCl) and 50 µl of Proteinase K (20mg/ml). Centrifuge briefly to make sure that there is no residual oil between the bead and the solution
- Incubate at 50°C overnight
- Carefully remove the solution as complete as possible
- Wash the samples 2x with 400 µl 1xTE for 15 min.
- Rinse the bead with 400 µl of 0.3M NaOH
- Wash the sample 2x with 400 µl of 0.3M NaOH
- Wash the samples 2x with 1xTE
- Remove the TE entirely and store the bead at 4°C (no longer than a couple of days)
- Proceed with bisulphite treatment
Agarose embedded bisulfite modification
Agarose beads formed according to protocols A or B are used for the following procedure:
- Prepare bisulfite solutions freshly as follows: mix 1.9 g sodium metabisulfite (Merck no. 1.06528.0500) with 2.5 ml ddH2O and add 750 µl of freshly prepared 2 M NaOH, dissolve the powder by strong vortexing (it takes time, if the powder does not dissolve it can be done at 50°C). Dissolve 55 mg of hydroquinon (Sigma) in 500 µl ddH2O at 50°C. Mix both solutions when chemicals are dissolved completely. (comment 2)
- Add 400 µl of freshly prepared bisulfite solution to each bead (as prepared according to the protocol A or B. Overlay with mineral oil (following protocol B). (When adding the solution to beads from protocol A make sure the beads are in contact with the solution and there is no residual oil between the bead and the solution; brief spin down helps)
- Transfer the samples to 50°C and incubate for 3.5 h
- Cool the samples briefly on ice (makes the bead better visible and “harder” for solution exchange)
- Remove the bisulfite solution and rinse the bead with 400 µl of 1x TE at RT (put 400 µl of 1x TE onto the bead and remove the TE immediately to avoid starting diffusion)
- Wash the samples 2x with 400 µl 1xTE for 15 min at RT
- Rinse the bead with 400 µl of 0.3M NaOH at RT
- Wash the sample 2x with 400 µl of 0.3M NaOH at RT
- Wash the samples 2x with 1xTE at RT
- Remove the TE entirely and store the bead at 4°C (can be stored for 2 weeks max.)
- Proceed with Bisulphite PCR amplification (Bis-PCR)
Bisulphite PCR amplification, cloning and sequencing
For Bisulphite PCR-amplification we regularly use the HotFirePol from Solis BioDyne, Tartu, Estonia. It combines the characteristics of a Hot Start Taq with a good performance und a high cost-efficiency. We experienced that some amplicons may work better with HotStarTaq from Qiagen (has to be tested).
Our standard Bis-PCR protocol runs as follows:
Melt the agarose bead at 70-80°C and transfer 3 µl of the molten agarose into the PCR tube containing the following ingredients (kept at RT or on ice):
- 5 µl 10xbuffer (Solis BioDyne)
- 5 µl MgCL2 (25 mM)
- 4 µl dNTPs (2.5 mM each)
- 1 µl forward primer (10 µM)
- 1 µl reverse primer (10 µM)
- 0.6 µl HotFirePol (5 U/µl)
- 30.4 µl ddH2O
Transfer into a preheated (95°C) PCR block and start the following program:
- 97°C 15 min (required for HotFirePol)
- 40 cycles
- 95°C 1 min
- 52°C 1 min
- 72°C 1min30sec
Load 5 µl on a 1.2% agarose gel and purify the PCR product by gel extraction (using JetsorbTM from Genomed) or precipitation with SureClean from Bioline, store rest at -20°C. (comment 3)
Cloning and “colony” PCR
We usually use 3 µl of the purified PCR product for pGEM-T T/A cloning, following the recommendations of the manufacturer, and plate on LB agar containing 100 mg/ml Amp, 0.04% X-GAL, 0.2 mM IPTG. Colonies are picked after overnight incubation with a sterile tooth pick. Briefly dip the tooth pick into a PCR-tube containing 30 ml of the following Colony-PCR Mix :
- 3 µl 10x Standard PCR buffer (containing 20mM MgCl)
- 2.4 µl dNTPs (2.5 mM each)
- 1 µl colony PCR forward primer (10 µM)
- 1 µl colony PCR reverse primer (10 µM)
- 0.5 µl Taq polymerase ( U/µl)
- 22.1 µl ddH2O
For amplification we use the Colony-PCR program:
- 95°C 5 min
- 30 cycles of:
- 95°C 30 sec
- 65°C 30 sec
- 72°C 1min
- 72°C 5 min
Comments to cloning and colony PCR: We use pGEM-T specific “colony PCR” forward (5´ gct att acg cca gct ggc gaa agg ggg atg tg 3´) and reverse (5´ ccc cag gct tta cac ttt atg ctt ccg gct cg 3´) primers which are located in the vector approx. 150 bp upstream and downstream of the insert hence enlarging the insert by 300bp. The “colony PCR fragments” include the universal M13 forward and reverse primer sequences which make them amendable for Sanger sequencing.
For cloning we routinely use the pGEM-T cloning kit from Promega, but PCR products may behave different with respect to cloning efficiency, so that it might be occasionally useful test other vector systems. We even experienced that some fragments are not clonable at all into a particular vector. We do not find much of a difference with respect to the bacterial strain used (such as Sure, Top10). Different ligation times and temperatures may also influence the cloning efficiency. Best results we obtained incubating the ligation mix at 16°C over night.
The particular sequence composition of bisulphite converted DNA often results in repeated structures and extended homopolymer stretches. This might influence the orientation of the cloned fragment and occasionally the maintenance of the cloned fragments in bacteria (instabilities). We recommend to pick clones from fresh overnight plates. (comment 4)
Sequencing of bisulphite clones
We perform the sequencing following standard Sanger-procedures (Sanger et al., 1977) on a Beckmann-Coulter sequencer using the DTCS QickStart Sequencing kit. In principle all standard sequencing platforms can be used such as ABI 3100, 3730, Amersham Megabase etc. For sequencing we usually prepare (“purify”) the colony-PCR templates by PEG precipitation (protocol available on request) or use gel purification. For a standard sequencing reaction about 100-200ng of template is needed. We mostly use standard sequencing primers (M13 forward: 5´ gtt ttc cca gtc acg acg 3´, M13 reverse: agg aaa cag cta tga cca t 3´) for the reaction but occasionally perform reactions with template specific primers.
One of the most critical issues for successful bisulphite experiments is to carefully select the regions suitable for bisulphite analysis and the design of accurately positioned and matching primers. Our experience is that there is no good specific prediction program for the design of bisulphite PCR primers available. Hand made selections (with some “experience”) work best. (comment 5) Selected primers should however be examined with software such as Oligo 6.0; Vector NTI etc. manually for similarities in melting temperatures, secondary structures and primer dimer formation. We usually vary the oligo length to adjust annealing performance and avoid secondary structure effects. Primers work best between 25-28 bp in length. Both primers should have an optimal Tm at around 65°C which means an effective annealing temperature of 52°C in the PCR reaction. (comment 6) As a threshold we recommend the primers do not differ more than 1-2°C in their Tm’s. The melting temperature of the PCR product should also not differ too much from the melting temperature of the primers.
Primers should be designed within originally cytosine rich but CpG depleted regions. One should try to include at least 4-5 cytosines preferably at the 3´end of the primer which will be converted to T’s. This enhances the “specificity” to truly amplify bisulphite converted DNA sequences and to “eliminate” not fully converted material. CpG positions within primers should be avoided - our experience with the use of “wobble” primers (to cover CpG positions with either C or T) were not very positive, i.e. we observed preferential amplification of methylated templates.
In older protocols (Hajkova et al., 2002) we recommended nested PCR approaches to enhance sensitivity and specificity. The current protocol yields better DNA quality and we mostly find the nested approach to be unnecessary. Avoiding nested amplifications has advantages because it reduces the risk of cross contamination and of preferential primary template amplification (amplification bias). According to our experience in almost all cases one round of PCR (max. 40 cycles) is sufficient, even if small amount of material is used (analysis of single copy genes from 20 cells, i.e. 60 pg genomic DNA). (comment 7) (figure 1)
BiQ-Analysis of Bisulphite-sequence data (clones)
The conventional procedure to analyze bisulfite-modified DNA is to amplify the region of interest, to clone the PCR product and analyze single clones by Sanger-sequencing. By this high resolution single chromosome profiles are obtained. A big issue in this process is the evaluation (quality) and quantification of sequencing profiles. To facilitate the analysis we have developed the BiQ Analyzer program (Bock et al., 2005, which runs on PC and Mac’s and is freely available at http://biq-analyzer.bioinf.mpi-sb.mpg.de/). The program uses ClustalW to automatically align the single amplicon sequences to the original genomic sequence. It calculates sequence identity, conversion rate and methylation degree for each clone and provides a colour coded alignment (see figure 2a). Sequences below 80% sequence identity to the original sequence and below 90% conversion rate are marked and can be excluded from analysis. The data are exported in excel type tables and automatically provide a graphical representation (“Lollipop” graphs) (figure 2b) and a bar diagram calculating the methylation degree for every CpG position in all clones are given as output files (figure 2c). Whole methylation of the amplicon is calculated by the following formula: no. of methylated sites/ no. of all analyzed sites. BiQ Analyzer saves time and provides standardized data and output formats. BiQ still requires some manually quality control of the processing - the program sometimes does not exclude “bad” sequence traces that show very low sequence quality.
The statistical significance of the bisulphite results is very much dependent on the number of sequences analyzed. This number can vary from experiment to experiment depending on the methylation variation observed. We recommend analyzing at least some 15-20 clones for each amplicon, ideally derived from two independent bisulphite treatments to obtain a first reasonable estimate. Ideally the data obtained by bisulphite clone sequencing should be confirmed by a different method since cloning bias (i.e. preferred cloning of unmethylated or methylated sequences) and preferential amplification may play a role (see also Clark et al., 2006).
Brief overview of other methods for sequence analysis
Besides labor and cost intensive single clone analysis by Sanger sequencing it is also possible to analyze the average methylation of CpG positions within the PCR amplicon by “direct” sequencing of the bisulphite-PCR product. Here pyrosequencing has been shown to be very efficient and quantitative but reaching only few CpG positions (Tost et al., 2007). Chromatograms from direct sequencing by the Sanger-method are rather difficult to quantitate without sophisticated software, i.e. such an analysis requires software such as ESME (Lewin et al., 2004; see also http://www.epigenome.org/), which are not publicly available. If only single or few CpG position shall be interrogated SNuPE based methods such as SIRPH (El-Maarri et al., 2002, El-Maarri, 2004) or real time PCR approaches such as MethyLight (Trinh et al., 2001) can be recommended since they are less laborious and once established very reliable and quantitative.
SIRPH as an alternative to sequencing
In case the exact CpG positions to be quantified for methylation are known in a given fragment, SIRPH (and other SNuPE-based methods) offer(s) a cheaper and less laborious possibility to obtain average methylation values. We are using this method routinely for large scale sample comparison. Briefly, a primer is placed next to the CpG position to be interrogated. Extension by ddCTP (or ddGTP on the other strand) indicated methylation, extension by dTTP an unmethylated state. The ratio of incorporated nucleotides (extended primers) can be examined simultaneously on a DHPLC WAVETM system from Transgenomic (using DNASep column). T and C extended primers are simply separated according to their hydrophobicity and can be distinguished from unextended primers by length. By measuring the height or area of the obtained peaks the ratio C/C+T is calculated. Further information on the method is available on request.
Although SIRPH is a simple, quick and cost-efficient method it can only analyze one or two CpG positions simultaneously. Therefore it is not suitable for the analysis of complex methylation patterns. Moreover to really quantify the methylation for a given position requires the establishment of a reference curve – since the C and T incorporations at a specific position may vary form position to position. Additionally, conversion rate (i.e. the bisulphite quality check) will have to be tested using an additional primer or partial sequencing of clones.
Reviewed by: Jörg Tost Centre National de Génotypage, CEA / Institut de Génomique, 91057 Evry Cedex, France.
- Supplier should be specified. We are using SeaPlaqueTM low melting agarose from FMC.
- We have experienced no difference in performance if the solution stays slightly clouded and contains some undissolved crystals.
- A simple and quick alternative for purification consists of a centrifugation step of the PCR product through a Multiscreen plate (Millipore) filled with P100 Bio-Gel (Bio-Rad) acrylamide suspension.
- This is an important point and we have frequently experienced problems with instability of inserts. Shorter incubation times and even lower concentrations of antibiotics might help to overcome this problem. In some cases a slightly lower incubation temperature 36°C might have beneficial effects.
- We have good success rates using freely available design software such as MethPrimer, which is based on the widely used Primer3 algorithm or MethylPrimer express from Applied Biosystems.
- We use an annealing temperature of 60°C, which in our hand reflects quite well the real temperature when using the Qiagen HotStar Taq.
- We do also not recommend any nested PCR approaches.
Fig. 1 An example of proper (green) and improper (red) primer positioning. Highlighted in green is the best position to place the primer. This region contains originally 7 cytosines where at least four of them are located at the 3´ part of the primer. The region highlighted in red shows a position where in total 5 originally cytosines are included, but they are all located in the 5´part of the primer. In this case it might happen that not fully converted sequences are amplified as well.
Fig. 2a Example of a BiQ generated ClustlW alignment representing the methylation of 7 CpG positions in single clones; unconverted Cs are highlighted in green.
Fig. 2b “Lollipop” scheme representing unmethylated (white circles) and methylated (black circles) CpG positions; this diagram is generated based on the ClustlW alignment including the relative CpG distances.
Fig. 2c Bar diagram showing the numbers of methylated and unmethylated CpG positions; the unmethylated portion is shown in blue, the methylated in yellow; not present designates CpGs that were not analysable due to sequencing errors; the distance between the single CpG positions is given in bp.
- Frommer M, et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. (1992) Proc Natl Acad Sci U S A. :1827-31.
- Hajkova, P., El-Maarri, O., Engemann, S., Oswald, J., Olek, A. and Walter, J. (2000) DNA-methylation analysis by the bisulfite-assisted genomic sequencing method Methods Mol. Biol. 200: 143-154.
- Sanger, F., Nicklen, S. and Coulson, A.R. (1977) DNA sequencing with chain-terminating inhibitors Proc Natl Acad Sci USA 74: 5463-5467.
- Bock, C., Reither, S., Mikeska, T., Paulsen, M., Walter, J. and Lengauer, T. (2005) BiQ Analyzer: visualization and quality control for DNA methylation data from bisulfite sequencing Bioinformatics 21: 4067-4068.
- Lewin, J., Schmitt, A.O., Adorjan, P., Hildmann, T. and Piepenbrock, C. (2004) Quantitative DNA methylation analysis based on four-dye trace data from direct sequencing of PCR amplificates Bioinformatics 20: 3005-3012.
- El-Maarri, O., Herbiniaux, U., Walter, J. and Oldenburg, J. (2002) A rapid, quantitative, non-radioactive bisulfite-SNuPE- IP RP HPLC assay for methylation analysis at specific CpG sites Nucleic Acids Res. 30: e25.
- El-Maarri, O. (2004) SIRPH analysis: SNuPE with IP-RP-HPLC for quantitative measurements of DNA methylation at specific CpG sites Methods Mol. Biol. 287: 195-205.
- Trinh, B.N., Long, T.I. and Laird, P.W. (2001) DNA methylation analysis by MethyLight technology Methods 25: 456-462.
- Clark SJ, Statham A, Stirzaker C, Molloy PL, Frommer M. (2006) DNA methylation: bisulphite modification and analysis Nat Protoc: 2353-64.
- Tost J, Gut IG.DNA methylation analysis by pyrosequencing. Nat Protoc. 2007;2(9):2265-75.