We thank all the authors of the original ERRBS report. We thank Mame Fall for technical assistance. We acknowledge the Weill Cornell Medical College Epigenomics Core for technical services and assistance. The work was supported by a Sass Foundation Judah Folkman Fellowship, an NCI K08CA169055 and ASH-AMFDP12005 to FGB, NIH R01HG006798 and R01NS076465, funding from the Irma T. Hirschl and Monique Weill-Caulier Charitable Trusts and STARR Consortium (I7-A765) to CEM, and an LLS SCORE grant (7006-13) to AMM.

NOTE: Institutional review board approval was obtained at Weill Cornell Medical College (protocol number 0805009783) and this study was performed in accordance with the Helsinki protocol.
NOTE: Please consult material safety data sheets for relevant materials before use (indicated throughout the protocol with &#8220;CAUTION&#8221;). Several of the reagents used are toxic and appropriate safety measures are advised (personal protective equipment and fume hood).
NOTE: All steps are performed at room temperature unless otherwise indicated throughout the protocol
1. Preparation and digestion of genomic DNA
<ol>
	<li>Prepare 50 ng of high quality genomic DNA (&#62; 40 kilobases in size for human DNA) as starting material in 50 microliter (&#181;l) of DNase-free water. Quantify the DNA using a fluorescence-based quantitation assay per manufacturer&#8217;s recommendations.</li>
	<li>Mix 50 ng of DNA with 2 &#181;l of MspI (100,000 units/milliliter), 10 &#181;l of appropriate 10X reaction buffer, and DNase-free water to bring the reaction to a total volume of 100 &#181;l.</li>
	<li>Incubate the reaction at 37 &#176;C in a thermal cycler or water bath for at least 18 hours.</li>
	<li>DNA extraction and precipitation
	<ol>
		<li>Purify the digested DNA using a phenol-chloroform extraction.
		<ol>
			<li>Add 200 &#956;l of 10 millimolar (mM) tris(hydroxymethyl)aminomethane (Tris-Cl) pH 8.0 buffer to each reaction from step 1.3 to bring the volume to 300 &#956;l.</li>
			<li>Add 150 &#181;l of Tris- Ethylenediaminetetraacetic acid (TE) -saturated phenol and 150 &#181;l of chloroform mix (1:1; CAUTION) to the DNA in a chemical hood and vortex briefly.</li>
			<li>Centrifuge at 20,800 x g in a microcentrifuge for 10 min at room temperature.</li>
			<li>Transfer top aqueous phase (approximately 300 &#181;l) from last step into a new tube.</li>
		</ol>
		</li>
		<li>Precipitate the digested DNA using ethanol precipitation.
		<ol>
			<li>Add 1 &#181;l of glycogen (20 milligrams/milliliter), 1/10 volume (30 &#181;l) of 3 M pH 5.2 sodium acetate and 2.5 volumes (750 &#181;l) of room temperature 100% ethanol. Mix by vortexing at high speed.</li>
			<li>Centrifuge at 4 &#176;C for 45-60 min at 20,800 x g in a microcentrifuge.</li>
			<li>Remove ethanol quickly by inverting and dragging each tube over a paper towel with a quick flick of the wrist. Ensure that the DNA pellet remains on the side wall of the tube.</li>
			<li>Add 600 &#181;l of room temperature 70% ethanol to each tube. Visualize the pellet as it dislodges from the tube side. Mix gently by inverting the tube five times.</li>
			<li>Centrifuge at 4 &#176;C for 45-60 min at 20,800 x g in a microcentrifuge.</li>
			<li>Remove the ethanol by quickly inverting and dragging each tube over a paper towel with a quick flick of the wrist.</li>
			<li>Remove as much ethanol as possible by carefully removing any residual volume. Use a vacuum apparatus with a DNase-free 10 &#181;l pipette tip at the end, and allow the pellet to air dry if necessary to fully eliminate any further residual ethanol. 
			NOTE: Avoid over drying the pellet since that will reduce the solubility upon resuspension in the next step.</li>
			<li>Resuspend the DNA pellet into 30 &#181;l of 10 mM Tris-Cl, pH 8.5 (Buffer EB). 
			NOTE: To avoid re-annealing of the &#8220;sticky&#8221; CG overhangs created by the MspI digestion, it is necessary to complete the protocol through the ligation step on the same day (protocol step 4).</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
2. End-repair
<ol>
	<li>Transfer the 30 &#181;l of MspI digested DNA (from step protocol 1.4.2.8) into a PCR tube on ice (4 &#176;C) and add the reagents from Table 1.</li>
	<li>Incubate the end-repair reaction for 30 min at 20 &#176;C in a thermal cycler with the heat lid on.</li>
	<li>Purify the DNA products to remove unincorporated dNTPs and other reaction reagents using a commercial column-based kit that binds of DNA in high-salt buffer and elutes DNA in low-salt buffer conditions per manufacturer&#8217;s recommendations. Elute the DNA products in 32 &#181;l of Buffer EB.</li>
</ol>
3. A-tailing
<ol>
	<li>Transfer the 32 &#181;l of DNA solution from step 2.3 into a PCR tube on ice (4 &#176;C) and add the reagents listed in Table 2.</li>
	<li>Incubate the A-tailing reaction for 30 min at 37 &#176;C in a thermal cycler with the heat lid on.</li>
	<li>Purify the DNA products to remove unincorporated dATPs and other reaction reagents using a commercial column-based kit that binds of DNA in high-salt buffer and elutes DNA in low-salt buffer conditions per manufacturer&#8217;s recommendations. Elute in 10 &#181;l of Buffer EB.</li>
</ol>
4. Adapter ligation
<ol>
	<li>Transfer the 10 &#181;l of A-tailed DNA into a PCR tube on ice (4 &#176;C)</li>
	<li>Add the ligation reaction reagents to the DNA as detailed in Table 3 and adapters as detailed in Table 4.</li>
	<li>Incubate the ligation reaction overnight at 16 &#176;C in a thermal cycler with the heat lid on.</li>
	<li>Purify the ligation products using a solid-phase reversible immobilization (SPRI) bead isolation protocol per manufacturer&#8217;s recommendations. For example, with Agencourt AMPure XP beads, use a 1.8X ratio of bead volume to sample volume (90 &#181;l of beads for the 50 &#181;l ligation reaction). Elute into 30 &#181;l of DNase-free water.
	<ol>
		<li>Ligated DNA can be stored at -20 &#176;C before proceeding with size selection.</li>
	</ol>
	</li>
</ol>
5. Size Selection
NOTE: Follow either section 5.1 for an automated size selection protocol or section 5.2 for a manual gel extraction size selection protocol. For samples with 25 ng or more of input DNA, an automated size selection protocol (using an instrument such as the Pippin Prep) can be used. Manual gel extraction is necessary for low DNA input amounts of 5-10 ng.
<ol>
	<li>Size selection using Pippin Prep for samples with 25 ng or higher of input DNA (Figure 2A) .
	<ol>
		<li>Create a new protocol on the Pippin Prep for size selection.</li>
		<li>Select &#8220;2% DF Marker L&#8221; as the Cassette. Click the &#8220;Use internal standards&#8221; button. Verify that the &#8220;Ref lane&#8221; numbers match the lane numbers.</li>
		<li>Select &#8220;Range&#8221; as the collection mode for each lane. Enter 135 under &#34;BP Start&#34;, 410 under &#34;BP End&#34;, and 240 under &#34;BP Pause&#34; for each lane.</li>
		<li>Save the protocol.</li>
		<li>Add 10 &#181;l of marker L to each 30 &#181;l sample from the ligation reaction (step 4.3), bringing the total volume to 40 &#181;l. Follow standard Pippin Prep protocol to prepare the gel cassette and instrument.</li>
		<li>Remove 40 &#181;l of electrophoresis buffer from the sample wells, and load a 40 &#181;l sample from step 5.1.5 into each of the 5 wells of the dye-free gel cassette. 
		NOTE: Dye-free, agarose gel cassettes are critical to size selecting adapter-ligated DNA fragments as the presence of ethidium bromide can alter the migration properties of the forked-adapter bound fragments.</li>
		<li>Select the protocol created in step 5.1.1, and start the run.</li>
		<li>For each lane used collect 40 &#181;l from the elution module when the run pauses at the 240 base pair (bp) point (lower library fraction: 135-240 bp). The instrument will pause individually for each lane. Repeat steps 5.1.5-5.1.7 for each lane as it pauses at 240 bp.</li>
		<li>Wash the elution module with 40 &#181;l of fresh electrophoresis buffer by pipetting up and down three times. Discard the 40 &#181;l of wash buffer. Repeat the wash step two additional times and remove all residual liquid from the elution module. Washing the elution module will decrease the amount of lower fraction DNA carried over into the higher fraction collection.</li>
		<li>Add 40 &#181;l of fresh electrophoresis buffer and re-seal the elution module and resume the run.</li>
		<li>Collect 40 &#181;l from each elution module into a new tube when the instrument indicates that the elution is complete. This elution contains the higher library fraction (240-410 bp). If running multiple gel cassettes, store size-selected samples at 4 &#176;C before proceeding with bisulfite conversion.</li>
	</ol>
	</li>
	<li>Size selection usingmanual gel extraction (Figure 2B)
	<ol>
		<li>Prepare a 1.5% agarose gel with 0.2 &#181;g/ml ethidium bromide (CAUTION).</li>
		<li>Load 50 bp and 100 bp ladders prepared with loading buffer in adjacent wells on both sides of the gel.</li>
		<li>Add 2 &#956;l of 6X Orange G loading dye to the cleaned-up ligation products. Load the entire volume of each sample into individual wells, skipping at least one well between samples to avoid cross-contamination.</li>
		<li>Run the gel at 3.5 volts per centimeter until ladder is fully separated and the Orange G dye runs to the bottom of the gel (minimum of one hour)</li>
		<li>Slice off the ladders using a clean razor blade and visualize them on a UV transilluminator (CAUTION). Do not expose the samples to UV light. Mark the 150 bp, 250 bp, and 400 bp bands of the ladders using a razor blade or pipette tips. Marking the ladder bands will allow the excision of samples without exposing them to UV light.</li>
		<li>Align the ladder slices with the rest of the gel. Using the marked bands as reference, excise one slice containing the 150-250 bp library fraction (lower) and another slice containing the 250-400 bp library fraction (higher). Place each of the two slices into different tubes.</li>
		<li>Repeat the excision for all samples prepared using a clean razor blade for each gel lane used.</li>
		<li>Purify the lower and higher library fractions using a gel extraction kit following manufacturer&#8217;s protocols. Elute each sample extracted into 40 &#181;l of EB. To ensure efficient PCR amplification of both the lower and higher library fractions in step 7, maintain the two fractions as independent samples for the bisulfite conversion step.</li>
	</ol>
	</li>
</ol>
6. Bisulfite Conversion
<ol>
	<li>Set up a bisulfite conversion control per manufacturer&#8217;s recommendations (e.g., Universal Methylated Human DNA Standard kit). Treat the control in the same manner as the samples throughout the bisulfite conversion protocol.</li>
	<li>Perform bisulfite conversion using a commercial kit protocol per manufacturer&#8217;s recommendations. Use a commercial kit that has limited template degradation and DNA loss during treatment and cleanup, while resulting in high bisulfite conversion rates of the DNA. If using the EZ DNA Methylation Kit, perform the protocol per manufacturer&#8217;s recommendations with the following exception: incubate the samples in a thermal cycler using the following protocol: 55 cycles: 95 &#176;C for 30 sec, 50 &#176;C for 15 min. Hold at 4 &#176;C.</li>
	<li>Elute the samples in 40 &#181;l of DNase free water. Proceed with the enrichment PCR (step 7) on the same day that the bisulfite conversion is completed. Bisulfite-treated DNA is AU-rich and single-stranded, which reduces its stability.</li>
	<li>Confirm efficient bisulfite conversion of the control used via Sanger sequencing with non-methylated cytosines greater than 99% converted.</li>
</ol>
7. Enrichment PCR
<ol>
	<li>Prepare a PCR master mix using the reagents in Table 5 per library fraction (optimized for the use of FastStart Taq DNA Polymerase).</li>
	<li>Add the PCR master mix to each 40 &#181;l bisulfite-converted library fraction. Mix by pipetting. Divide the 200 &#181;l reaction into four PCR tubes with 50 &#181;l each.</li>
	<li>Amplify the reactions in a thermal cycler with the following protocol: Set the heat lid to 100 oC. Initialize with a step of 94 &#176;C for 5 min. Run 18 cycles of denaturing, annealing and extension/elongation steps: 94 &#176;C for 20 seconds followed by 65 &#176;C for 30 seconds followed by 72 &#176;C for 1 min. Run a final extension/elongation step of 72 &#176;C for 3 min and hold at 4 &#176;C. 
	NOTE: 18 cycles of PCR is recommended for low input material quantities (less than 10 ng) and for first time users of the protocol. The protocol can be adjusted for a lower number of PCR cycles (as low as 14 cycles for 50 ng of input DNA; see Table 6).</li>
</ol>
8. Purify PCR reactions
<ol>
	<li>Combine each 4-reaction set (four 50 &#181;l PCR reactions per library fraction).</li>
	<li>Purify the PCR products using a SPRI bead approach (or other approach which can remove unincorporated primers and other PCR reaction reagents) per manufacturer&#8217;s recommendations. The following steps have been optimized for Agencourt AMPure XP.
	<ol>
		<li>Add 1.7X volume of SPRI beads (340 &#956;l for 200 &#956;l of PCR product) to the lower fraction amplification products.</li>
		<li>Add 1.1X volume of SPRI beads (220 &#956;l for 200 &#956;l of PCR product) to the higher fraction amplification products.</li>
		<li>Purify PCR products per manufacturer&#8217;s recommendations using 800 &#181;l of 70% ethanol for wash steps.</li>
		<li>Elute into 50 &#181;l of DNase-free water and add 1M Tris buffer to bring each library to 10 mM Tris-Cl, pH 8.5. Store libraries at -20 &#176;C.</li>
	</ol>
	</li>
</ol>
9. Library Quality Control
<ol>
	<li>Quantify libraries using a fluorescence-based quantitation assay selective for double-stranded DNA per manufacturer&#8217;s recommendations. Spectrophotometry-based measurements are not reliable. Expected concentrations measure 10 - 50 ng/&#181;l for the lower library fraction and 3 - 15 ng/&#181;l for the higher library fraction.</li>
	<li>Assess library sizes and quality using a bioanalyzer instrument and the High Sensitivity DNA Kit. Visualize library products and determine the average size of each library fraction. The lower library fraction typically has an average size between 180 and 210 bp. The higher fraction typically has an average size between 280 and 310 bp.</li>
</ol>
10. Prepare libraries for sequencing
<ol>
	<li>Calculate the molarity of each library fraction as follows:
	<ol>
		<li>Size each library fraction using the trace obtained in step 9.2. For example, if using a bioanalyzer, use the &#8220;region&#8221; feature, and cover the beginning and end of the library fraction assessed.</li>
		<li>Record the average size of each library fraction in bp.</li>
		<li>Calculate molarity of each library fraction (Nanomolar concentration of DNA) using the following formula: Nanomolar concentration of DNA = [(ng/1000)/(bp*660)]*10^9 where ng is the concentration expressed in ng/&#956;l (as measured in step 9.1) , bp is the average size of the library fraction, and 660 is the weight of a standard double stranded DNA base pair. For example: with a library that is 326 bp in size, at a concentration of 14.2 ng/&#956;l, the Nanomolar concentration of DNA is 66 nM.</li>
	</ol>
	</li>
	<li>Prepare 10 &#181;l of a 2 nM solution of each library fraction. Dilute the library with DNase-free water.</li>
	<li>Pool the library fractions by combining 10 &#181;l of the 2 nM lower fraction with 10 &#181;l of the 2 nM higher fraction for each sample library prepared. This pool is the final ERRBS library for sequencing.</li>
	<li>Sequence the ERRBS library
	<ol>
		<li>Load the libraries with a goal of optimal densities of 600,000 &#8211; 650,000 clusters per millimeter squared (suggested loading concentration of 7-8 picomolar).</li>
		<li>Sequence the libraries using a minimum of single-read 51 cycles sequencing on a HiSeq 2500 sequencer in high output mode.</li>
		<li>Use a designated control lane (see discussion for rationale).</li>
	</ol>
	</li>
</ol>
11. Data Analysis
NOTE: Please refer to the Supplemental code files 1 and 2 for full details of commands and scripts recommended for use.
<ol>
	<li>Convert the base call files (.bcl files) to individual FASTQ files for each sample using the software provided by the sequencer&#8217;s manufacturer (example: CASAVA v1.8.2 for Illumina). If a sequencing core facility is used, this step may be provided and/or performed as an automated step on-board the sequencing computer.</li>
	<li>Filter the sequencing reads for reads that pass quality filtering (see Supplemental code file 1 for details). The CASAVA FASTQ file contains both reads that pass quality filtering and reads that do not pass quality filtering. Use a custom script that utilizes the &#60;is filtered&#62; element of the sequence identifier and keeps reads that have passed quality filtering.</li>
	<li>Trim adapter sequences from 3&#39; end of the sequence reads in the filtered FASTQ files using a software tool that can remove adapter sequences, such as Flexbar 37. The minimum overlap length (-ao) set to 6, the minimum read length to keep after adapter removal (--m) set to 21, the cutoff for the number of allowed mismatches (-at) set to 2, and defaults for all other parameters. As an alternate to Flexbar, any software which can remove adapters can be used in step 11.3 (examples: cutadapt38, trimmomatic39).</li>
	<li>Use Bismark41 to align the filtered, adapter trimmed sequence reads to the bisulfite converted reference human genome hg19 (whole genome alignment approach) and determine the methylation context for each cytosine. Bismark is a customized short read mapping tool that aligns bisulfite treated reads to a bisulfite converted genome (where all &#8216;C&#8217;s are converted to &#8216;T&#8217;s) and returns methylation calls for cytosines in CpG, CHG, and CHH context. Typically set the seed length (-l) to the read length for alignment accuracy and use defaults for all other parameters.</li>
	<li>Sort aligned reads first by chromosome, then start position, and finally strand, after alignment is complete.</li>
	<li>Use custom scripts (see Supplemental code file 2 for script commands) to iterate over the sorted methylation calls output by Bismark to compute the percent methylation scores of bisulfite converted Cytosines (T&#8217;s; representing unmethylated C&#8217;s) and non-converted C&#8217;s (representing methylated C&#8217;s) for each cytosine methylation context (CpG, CHG, CHH), retaining only the cytosines that have at least phred quality score of 20 and have at least 10x coverage. 
	NOTE: The outputs are a methylation score file for each cytosine context with columns corresponding to: the position, strand, coverage, percent cytosine, and percent thymine, and for CHG and CHH context, a column for the next base (the H).</li>
	<li>Use a custom script (see Supplemental code file 2 for script commands) to compute the conversion rates, from the methylation scores, and output the total of other C&#8217;s considered (CHG and CHH context), average conversion rate, and median conversion rate. Compute conversion rates for both strands independently as well as summarized over the entire library. Mean conversion rate is the fraction of C&#8217;s (converted and unconverted) in non-CpG context of the total number of C&#8217;s. Cytosines in CHG or CHH context are typically unmethylated and therefore present as thymine in the sequencing data.</li>
	<li>Generate a BAM file of the aligned reads using Bismark&#8217;s bismark2SAM_v5_xm.pl41 and SAMTOOLS42, and generate a wiggle format file using a custom script to convert the CpG methylation calls. The post-processed output can be converted to other formats, such as a bedgraph, and viewed in a genome browser such as the UCSC Genome Browser43 or IGV44.</li>
</ol>

The authors have no conflicts of interest to disclose.

The protocol presented yields base-pair resolution data of cytosine methylation at biologically-relevant genomic regions. The protocol as written is optimized for 50 ng of starting material, however, it can be adapted to handle a range of input material (5 ng or more)26. This will require adjustments of some of the protocol steps as seen in Table 6. The ERRBS libraries are amenable to paired end sequencing and further genomic coverage can also be accomplished by sequencing reads longer than 51 cycles. Multiplexed sequencing will offer a lower cost protocol per sample, however, this will result in reduced coverage per CpG site represented in the data (Figure 5 and Table 8), and will not yield sufficient depth of coverage to perform analyses which require high coverage per CpG site (e.g. as described by Landan et al.33). Finally, this protocol (or any bisulfite-based protocol) cannot distinguish between methyl-cytosine and hydroxymethyl-cytosine46,47. However, the data generated can be integrated with other protocol results48,49 to delineate the different modifications, and other cytosine modifications recently reported50, should they be of interest.
High quality libraries will appear as shown in Figure 3A-C, and once pooled for sequencing yields a trace as shown in Figure 3G (red trace) representing equal molar contributions from both library fractions. Library preparation failure can result from any step during the procedure. If degraded DNA is processed it will result in libraries that are not enriched in MspI fragments and hence in low CpG coverage using the sequencing parameters described in this protocol. If an enzyme is non-functional or inadvertently excluded from one of the reactions, the protocol will not yield the expected library. If the ligation reaction is inefficient, adaptors are at a higher concentration than expected, and/or the primers concentration used is a limiting reagent for the final amplification steps, library failure can occur. Excess adapters (seen as a peaks at ~150 bp in bioanalyzer results; Figure 3D-F) in the library will also interfere with sequencing due to the indiscriminate clustering of both the library and excess adapters. While such a library may sequence apparently normally, a significant portion of the reads will be merely adapter sequences. If excess adapters are observed in a library, it is best to repeat the library preparation if material is available using optimal input material to adapter quantity ratios. Finally, to ensure efficient PCR amplification of the libraries, the lower and higher library fractions are maintained as separate samples throughout the bisulfite conversion and PCR enrichment steps. Failure to do so yields differential efficiency of amplification during the PCR reaction of higher and lower fractions (as seen in Figure 3G blue trace) and the potential for unequal representation of the respective genomic loci covered in each library fraction during sequencing. The user may opt to include a quantitative PCR step immediately after the bisulfite conversion for further titration of optimal PCR cycles needed to amplify the libraries being generated.
ERRBS library preparation protocol has several key steps in which specific reagents are recommended. At the end-repair step, the use of a four-nucleotide dNTP mix allows for end-repair of any products not containing the CG overhang, such as those resulting from MspI enzymatic star activity and sheared DNA fragments present in the original DNA sample. This results in improved CpG representation in the results. At the ligation step it is critical to use a high concentration ligase (2,000,000 units/ml) and methylated adapters to ensure that the ligation reaction is efficient and that the bisulfite conversion does not influence the adapter sequences essential for accurate data alignment. At the PCR step, using a polymerase capable of amplifying bisulfite-treated GC-rich DNA fragments is necessary for high specificity. Finally, to ensure elimination of excess adapters and primers, SPRI bead purification (for example: Agencourt AMPure XP) is recommended rather than column based assays for ligation and PCR product isolations.
In order to generate high quality data, it is important to ensure efficient bisulfite conversion. The control presented offers the user the ability to determine conversion efficiency prior to sequencing. As an alternative, a non-human DNA such as lambda DNA can be used as an internal control (spike-in). Due to the differences in species, this type of a control can be directly included in downstream sequencing (e.g. as used by Yu, et al.34). However, if the spike-in is utilized, it cannot be used to determine conversion efficiency prior to library sequencing unless uniquely amplified and independently sequenced prior to library sequencing. The conversion rates determined are based on the methylation status at non-CpG sites. This may not be appropriate for use in the context of high cytosine methylation in non-CpG context (for example embryonic stem cells) and parallel samples or other means of assessing for conversion efficiency can be utilized for this purpose.
There are a few caveats to address that are unique to the sequencing of ERRBS libraries. The first three bases of the library fractions sequenced are nearly uniformly non-random due to the MspI recognition cut site (C^CGG; see Figure 4B, C). This results in the potential for significant data loss due to low quality reads resulting from poor cluster localization in spite of apparent high cluster density during sequencing. To overcome this barrier, include a high complexity library in an independent lane (PhiX control or other library type) as a dedicated control lane. High complexity libraries have ends containing a balanced representation of A, C, T and G in the first four bases sequenced. Suitable control lanes include libraries such as RNA-seq, ChIP-seq, whole genome sequencing, or a control offered by the sequencing machine manufacturer (e.g. PhiX Control v3). When designated as a control lane for the respective sequencing run, it can serve as the basis for the matrix generation which is utilized during the first four bases of sequencing to detect cluster positions. The higher quality reads captured will raise the mean coverage per CpG site by 5.2 (n = 4). Alternatively, this technical difficulty can also be overcome using a dark sequencing approach as previously described23. Other sequencing criteria follow standard operating procedures per manufacturer&#8217;s protocols. Finally, the coverage per CpG chosen for data analysis will be guided by the user and in part by the biological questions of interest. 10x coverage threshold affords a high coverage analysis approach, however this threshold can be lowered should that be of interest.
A full discussion of ERRBS data analysis is beyond the scope of this article, however, differentially methylated cytosines and regions can be determined using open source tools 31,51-53. Additional analysis considerations and approaches have been well-described54,55, and the reader is encouraged to search the literature for tools most appropriate to the analysis planned.
Compared to other published methods, ERRBS offers a four-day protocol which when performed as described yields high rates of reproducibility. It has been validated compared to the gold standard MassARRAY EpiTYPER26, is cost-effective for high coverage data, and is adaptable for various input material amounts (favorable for clinical sample processing and other cell types of low frequency) and sequencing approaches. It offers base-pair resolution at biologically relevant loci and can be used in integrative analyses with other techniques profiling genome-wide transcription factor binding, chromatin remodeling, epigenetic marks and other cytosine modifications of interest. ERRBS data use in such studies can contribute to a comprehensive molecular approach and allow for high dimensional analyses in the study of biological models and human disease.

DNA methylation at cytosine (5-methylcytosine) is an epigenetic mark critical in mammalian cells for a variety of biological processes, including but not limited to imprinting, X chromosome inactivation, development, and regulation of gene expression1-8. The study of DNA methylation patterns in malignant and other disorders has determined disease specific patterns and contributed to the understanding of disease pathogenesis and potential biomarker discoveries9-17.&#160;There are many protocols that interrogate the epigenome for DNA methylation status. These can be divided into affinity-based, restriction enzyme-based, and bisulfite conversion-based assays that utilize microarray or sequencing platforms downstream. Furthermore, there are a few protocols that bridge these general categories including, but not limited to, Combined Bisulfite Restriction Analysis18 and Reduced Representation Bisulfite Sequencing (RRBS19).
RRBS was originally described by Meissner et al.19,20. The protocol introduced a step to enrich GC-rich genomic regions followed by bisulfite sequencing, which resulted in quantitative base-pair resolution data that is cost effective21,22. The GC-rich regions are targeted by the MspI (C^CGG) restriction enzyme, and cytosine methylation is resolved by bisulfite conversion of cytosines (deamination of unmodified cytosines to uracil), followed by polymerase chain reaction (PCR) amplification. RRBS covered the majority of gene promoters and CpG islands in a fraction of the sequencing required for a whole genome; however RRBS had limited coverage of CpG shores and other intergenic regions of biological relevance. Several groups have published updated RRBS protocols since the original report that improve upon the methodology and resultant coverage of these genomic regions23-25. Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) includes library preparation modifications and an alternate data alignment approach26 when compared to RRBS. ERRBS resulted in a higher number of CpGs represented in the data generated and increased coverage of all genomic regions interrogated26. This method has been used to resolve DNA methylation patterns in human patient and other animal specimens26-30.
The ERRBS protocol described offers details on all steps needed for completion and data was generated using representative human DNA (samples were obtained from previously reported, de-identified patient samples31, and a CD34+ bone marrow sample from a normal human donor). The protocol includes an automated size selection process, which reduces the processing time per sample and allows for increased accuracy in library size selection. The protocol combines a series of established molecular biology techniques. High molecular weight DNA is digested with a methylation-insensitive restriction enzyme (MspI) followed by end-repair, A-tailing, and ligation of methylated adapters. Size selection of the GC-rich fragments is followed by bisulfite conversion and PCR amplification prior to sequencing. Bisulfite conversion has been previously described32 and detailed review of data analysis and applications is beyond the scope of this paper, however recommendations and references are included for the readers&#8217; use. The protocol can be performed over four days and is amenable to small input (50 ng or less) material amounts. The protocol as described yields data with high coverage per CpG site sufficient not only for differential methylation site and region determinations but also for epigenetic polymorphism detection as described by Landan, et al.33.

<table border="0" cellpadding="0" cellspacing="0" width="1216">
	<tbody>
		<tr height="42">
			<td height="42" style="height:42px;width:269px;">Name of Reagent/ Equipment</td>
			<td style="width:153px;">Company</td>
			<td style="width:104px;">Catalog Number</td>
			<td style="width:193px;">Comments/Description</td>
			<td style="width:497px;">URL</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">MspI</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">R0106M</td>
			<td style="width:193px;">100,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/r0106-mspi</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">NEBuffer 2</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">B7002S</td>
			<td style="width:193px;">Reaction buffer for MspI enzyme; protocol step 1.2</td>
			<td style="width:497px;">https://www.neb.com/products/b7002-nebuffer-2</td>
		</tr>
		<tr height="140">
			<td height="140" style="height:140px;width:269px;">Phenol solution</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">P4557</td>
			<td style="width:193px;">Equilibrated with 10&#160;mM Tris HCl, pH&#160;8.0; see safety and handling instructions at http://www.sigmaaldrich.com/catalog/product/sigma/p4557</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/p4557</td>
		</tr>
		<tr height="100">
			<td height="100" style="height:100px;width:269px;">Chloroform</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">C2432</td>
			<td style="width:193px;">See safety and handling instructions at http://www.sigmaaldrich.com/catalog/product/sial/c2432</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sial/c2432</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Glycogen</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">G1767</td>
			<td style="width:193px;">19-22 mg/ml</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/g1767</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">NaOAc</td>
			<td>Sigma-Aldrich</td>
			<td>S7899</td>
			<td style="width:193px;">3M pH 5.2</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/s7899</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Ethanol</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">E7023</td>
			<td style="width:193px;">200 proof, for molecular biology</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sial/e7023</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Buffer EB</td>
			<td style="width:153px;">Qiagen</td>
			<td style="width:104px;">19086</td>
			<td style="width:193px;">10 mM Tris-Cl, pH 8.5</td>
			<td style="width:497px;">http://www.qiagen.com/products/catalog/lab-essentials-and-accessories/buffer-eb</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">tris(hydroxymethyl)aminomethane (Tris)</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">T1503</td>
			<td style="width:193px;">prepare a 1M pH 8.5 solution</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/t1503</td>
		</tr>
		<tr height="60">
			<td height="60" style="height:60px;">Tris- Ethylenediaminetetraacetic acid (TE)</td>
			<td style="width:153px;">Sigma-Aldrich</td>
			<td style="width:104px;">T9285</td>
			<td style="width:193px;">Dilute to 1X buffer solution per manufacturer&#39;s recommendations</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/t9285</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">T4 DNA Ligase Reaction Buffer</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">B0202S</td>
			<td style="width:193px;">10X concentration</td>
			<td style="width:497px;">https://www.neb.com/products/b0202-t4-dna-ligase-reaction-buffer</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Deoxynucleotide triphosphate (dNTP) Solution Mix</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">N0447L</td>
			<td style="width:193px;">10 mM each nucleotide</td>
			<td style="width:497px;">https://www.neb.com/products/n0447-deoxynucleotide-dntp-solution-mix</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">T4 DNA Polymerase</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">M0203L</td>
			<td style="width:193px;">3,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/m0203-t4-dna-polymerase</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">DNA Polymerase I, Large (Klenow) Fragment</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">M0210L</td>
			<td style="width:193px;">5,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/m0210-dna-polymerase-i-large-klenow-fragment</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">T4 Polynucleotide Kinase</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">M0201L</td>
			<td style="width:193px;">10,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/m0201-t4-polynucleotide-kinase</td>
		</tr>
		<tr height="60">
			<td height="60" style="height:60px;width:269px;">QIAquick PCR Purification Kit</td>
			<td style="width:153px;">Qiagen</td>
			<td style="width:104px;">28104</td>
			<td style="width:193px;">Used for DNA product purification in protocol step 2.3</td>
			<td style="width:497px;">http://www.qiagen.com/products/catalog/sample-technologies/dna-sample-technologies/dna-cleanup/qiaquick-pcr-purification-kit</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">2&#39;-deoxyadenosine 5&#39;-triphosphate (dATP)</td>
			<td style="width:153px;">Promega</td>
			<td style="width:104px;">U1201</td>
			<td style="width:193px;">100 mM</td>
			<td style="width:497px;">http://www.promega.com/products/pcr/routine-pcr/deoxynucleotide-triphosphates-_dntps_/</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Klenow Fragment (3&#39;&#8594;5&#39; exo-)</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">M0212L</td>
			<td style="width:193px;">5,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/m0212-klenow-fragment-3-5-exo</td>
		</tr>
		<tr height="60">
			<td height="60" style="height:60px;width:269px;">MinElute PCR Purification Kit</td>
			<td style="width:153px;">Qiagen</td>
			<td style="width:104px;">28004</td>
			<td style="width:193px;">Used for DNA product purification in protocol step 3.3</td>
			<td style="width:497px;">http://www.qiagen.com/products/catalog/sample-technologies/dna-sample-technologies/dna-cleanup/minelute-pcr-purification-kit</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">T4 DNA Ligase</td>
			<td style="width:153px;">New England Biolabs</td>
			<td style="width:104px;">M0202M</td>
			<td style="width:193px;">2,000,000 units/ml</td>
			<td style="width:497px;">https://www.neb.com/products/m0202-t4-dna-ligase</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Methylation Adapter Oligo Kit</td>
			<td style="width:153px;">Illumina</td>
			<td style="width:104px;">ME-100-0010</td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="120">
			<td height="120" style="height:120px;width:269px;">Agencourt AMPure XP</td>
			<td style="width:153px;">Beckman Coulter</td>
			<td style="width:104px;">A63881</td>
			<td style="width:193px;">Used in protocol sections that implement magnetic bead purification steps (steps 4.3 and 8.2). Equilibrate to room temperature before use</td>
			<td style="width:497px;">https://www.beckmancoulter.com/wsrportal/wsrportal/wsr/research-and-discovery/products-and-services/nucleic-acid-sample-preparation/agencourt-ampure-xp-pcr-purification/index.htm?i=A63882#2/10//0/25/1/0/asc/2/A63882///0/1//0/</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">Pippin Prep Gel Cassettes, 2% Agarose, dye-free</td>
			<td style="width:153px;">Sage Science</td>
			<td style="width:104px;">CDF2010</td>
			<td style="width:193px;">with internal standards</td>
			<td style="width:497px;">http://store.sagescience.com/s.nl/it.A/id.1036/.f</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">Certified Low Range Ultra Agarose</td>
			<td>Bio-Rad</td>
			<td>161-3106</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.bio-rad.com/en-us/sku/161-3106-certified-low-range-ultra-agarose</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Tris-Borate-EDTA (TBE) buffer</td>
			<td>Sigma-Aldrich</td>
			<td>T4415</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/t4415</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Ethidium bromide solution</td>
			<td>Sigma-Aldrich</td>
			<td>E1510</td>
			<td style="width:193px;">10 mg/ml</td>
			<td style="width:497px;">http://www.sigmaaldrich.com/catalog/product/sigma/e1510</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">50 bp DNA Ladder</td>
			<td>NEB</td>
			<td>N3236S</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">https://www.neb.com/products/n3236-50-bp-dna-ladder</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">100 bp DNA Ladder</td>
			<td>NEB</td>
			<td>N3231S</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">https://www.neb.com/products/n3231-100-bp-dna-ladder</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Gel Loading Dye, Orange (6X)</td>
			<td>NEB</td>
			<td>B7022S</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">https://www.neb.com/products/b7022-gel-loading-dye-orange-6x</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">Scalpel&#160;Blade&#160;No.&#160;11</td>
			<td>Fisher Scientific</td>
			<td>3120030</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.fishersci.com/ecomm/servlet/fsproductdetail?position=content&#38;tab=Items&#38;productId=11876776</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">QIAquick Gel Extraction Kit</td>
			<td>Qiagen</td>
			<td>28704</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.qiagen.com/products/catalog/sample-technologies/dna-sample-technologies/dna-cleanup/qiaquick-gel-extraction-kit</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">EZ DNA Methylation Kit</td>
			<td style="width:153px;">Zymo Research</td>
			<td style="width:104px;">D5001</td>
			<td style="width:193px;">Used in protocol step 6.2</td>
			<td style="width:497px;">http://www.zymoresearch.com/epigenetics/dna-methylation/bisulfite-conversion/ez-dna-methylation-kits/ez-dna-methylation-kit</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">EZ DNA Methylation-Lightning Kit</td>
			<td style="width:153px;">Zymo Research</td>
			<td style="width:104px;">D5030</td>
			<td style="width:193px;">Alternative for step 6.2</td>
			<td style="width:497px;">http://www.zymoresearch.com/epigenetics/dna-methylation/bisulfite-conversion/ez-dna-methylation-lightning-kit</td>
		</tr>
		<tr height="60">
			<td height="60" style="height:60px;width:269px;">Universal Methylated Human DNA Standard</td>
			<td style="width:153px;">Zymo Research</td>
			<td>D5011</td>
			<td style="width:193px;">Used as bisulfite conversion control</td>
			<td style="width:497px;">http://www.zymoresearch.com/epigenetics/dna-methylation/methylated-dna-standards/universal-methylated-human-dna-standard</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">FastStart High Fidelity PCR System</td>
			<td style="width:153px;">Roche</td>
			<td style="width:104px;">03553426001</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://lifescience.roche.com/shop/products/faststart-high-fidelity-pcr-system#tab-0</td>
		</tr>
		<tr height="80">
			<td height="80" style="height:80px;width:269px;">Qubit dsDNA High Sensitivity Assay Kit</td>
			<td style="width:153px;">Life Technologies</td>
			<td style="width:104px;">Q32854</td>
			<td style="width:193px;">A fluorescence-based DNA quantitation assay; used in protocol steps 1.1, 9.1 and 10.1</td>
			<td style="width:497px;">http://www.lifetechnologies.com/order/catalog/product/Q32854</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">DynaMag-2 Magnet</td>
			<td style="width:153px;">Life Technologies</td>
			<td>12321D</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">https://www.lifetechnologies.com/order/catalog/product/12321D</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">High Sensitivity DNA Kit</td>
			<td style="width:153px;">Agilent Technologies</td>
			<td style="width:104px;">5067-4626</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.genomics.agilent.com/en/Bioanalyzer-DNA-RNA-Kits/High-Sensitivity-DNA-Analysis-Kits/?cid=AG-PT-105&#38;tabId=AG-PR-1069</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;width:269px;">2100 Bioanalyzer</td>
			<td style="width:153px;">Agilent Technologies</td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.genomics.agilent.com/en/Bioanalyzer-System/2100-Bioanalyzer-Instruments/?cid=AG-PT-106&#38;tabId=AG-PR-1001</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">PhiX Control v3</td>
			<td style="width:153px;">Illumina</td>
			<td>FC-110-3001</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.illumina.com/products/phix_control_v3.ilmn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">HiSeq 2500</td>
			<td style="width:153px;">Illumina</td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.illumina.com/systems/hiseq_2500_1500.ilmn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Pippin Prep</td>
			<td style="width:153px;">Sage Science</td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.sagescience.com/products/pippin-prep/</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Qubit 2.0 Fluorometer</td>
			<td>Life Technologies</td>
			<td>Q32872</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.lifetechnologies.com/order/catalog/product/Q32872</td>
		</tr>
		<tr height="40">
			<td height="40" style="height:40px;">TruSeq SR Cluster Kit v3-cBot-HS</td>
			<td>Illumina</td>
			<td>GD-401-3001</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.illumina.com/products/truseq_sr_cluster_kit_v3-cbot-hs.ilmn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">TruSeq SBS Kit v3-HS</td>
			<td>Illumina</td>
			<td>FC-401-3002</td>
			<td style="width:193px;"></td>
			<td style="width:497px;">http://www.illumina.com/products/truseq_sbs_kit_v3-hs.ilmn</td>
		</tr>
		<tr height="80">
			<td height="80" style="height:80px;">TruSeq RNA Sample prep</td>
			<td>Illumina</td>
			<td>RS-122-2001</td>
			<td style="width:193px;">Barcoded adapters used for multiplexing libraries; See Supplemental file for multiplexing protocol</td>
			<td style="width:497px;">http:/http://www.illumina.com/products/truseq-rna-access-kit.ilmn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Microcentrifuge</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Vortex Mixer</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Dry Block Heater</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Thermal Cycler</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Water Bath</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Gel electrophoresis system</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Electrophoresis power supply</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">Gel doc</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;width:269px;">UV or blue light transilluminator</td>
			<td></td>
			<td></td>
			<td style="width:193px;"></td>
			<td style="width:497px;"></td>
		</tr>
	</tbody>
</table>

enhanced reduced representation bisulfite sequencing for assessment of dna methylation at base pair resolution

DNA methylation pattern mapping is heavily studied in normal and diseased tissues. A variety of methods have been established to interrogate the cytosine methylation patterns in cells. Reduced representation of whole genome bisulfite sequencing was developed to detect quantitative base pair resolution cytosine methylation patterns at GC-rich genomic loci. This is accomplished by combining the use of a restriction enzyme followed by bisulfite conversion. Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) increases the biologically relevant genomic loci covered and has been used to profile cytosine methylation in DNA from human, mouse and other organisms. ERRBS initiates with restriction enzyme digestion of DNA to generate low molecular weight fragments for use in library preparation. These fragments are subjected to standard library construction for next generation sequencing. Bisulfite conversion of unmethylated cytosines prior to the final amplification step allows for quantitative base resolution of cytosine methylation levels in covered genomic loci. The protocol can be completed within four days. Despite low complexity in the first three bases sequenced, ERRBS libraries yield high quality data when using a designated sequencing control lane. Mapping and bioinformatics analysis is then performed and yields data that can be easily integrated with a variety of genome-wide platforms. ERRBS can utilize small input material quantities making it feasible to process human clinical samples and applicable in a range of research applications. The video produced demonstrates critical steps of the ERRBS protocol.

Watch this Scientific Journal Video about Enhanced Reduced Representation Bisulfite Sequencing for Assessment of DNA Methylation at Base Pair Resolution at JoVE.com

Enhanced Reduced Representation Bisulfite Sequencing for Assessment of DNA Methylation at Base Pair Resolution

Enhanced Reduced Representation Bisulfite Sequencing is a method for the preparation of sequencing libraries for DNA methylation analysis based on restriction enzyme digestion combined with cytosine bisulfite conversion. This protocol requires 50 ng of starting material and yields base pair resolution data at GC-rich genomic regions.

DNA methylation pattern mapping is heavily studied in normal and diseased tissues. A variety of methods have been established to interrogate the cytosine ...

enhanced-reduced-representation-bisulfite-sequencing-for-assessment

Research

JoVE Journal

Biology

25.2K Views.  Weill Cornell Medical College. Enhanced Reduced Representation Bisulfite Sequencing is a method for the preparation of sequencing libraries for DNA methylation analysis based on restriction enzyme digestion combined with cytosine bisulfite conversion. This protocol requires 50 ng of starting material and yields base pair resolution data at GC-rich genomic regions.