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.

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

In vivo Imaging of Intact Drosophila Larvae at Sub-cellular Resolution

Recent improvements in optical imaging, genetically encoded fluorophores and genetic tools allowing efficient establishment of desired transgenic animal lines have enabled biological processes to be studied in the context of a living, and in some instances even behaving, organism. In this protocol we will describe how to anesthetize intact Drosophila larvae, using the volatile anesthetic desflurane, to follow the development and plasticity of synaptic populations at sub-cellular resolution1-3. While other useful methods to anesthetize Drosophila melanogaster larvae have been previously described4,5,6,7,8, the protocol presented herein demonstrates significant improvements due to the following combined key features: (1) A very high degree of anesthetization; even the heart beat is arrested allowing for lateral resolution of up to 150 nm1, (2) a high survival rate of &gt; 90% per anesthetization cycle, permitting the recording of more than five time-points over a period of hours to days2 and (3) a high sensitivity enabling us in 2 instances to study the dynamics of proteins expressed at physiological levels. In detail, we were able to visualize the postsynaptic glutamate receptor subunit GluR-IIA expressed via the endogenous promoter1 in stable transgenic lines and the exon trap line FasII-GFP1. (4) In contrast to other methods4,7 the larvae can be imaged not only alive, but also intact (i.e. non-dissected) allowing observation to occur over a number of days1. The accompanying video details the function of individual parts of the in vivo imaging chamber2,3, the correct mounting of the larvae, the anesthetization procedure, how to re-identify specific positions within a larva and the safe removal of the larvae from the imaging chamber.

In vivo Imaging of Intact Drosophila Larvae at Sub-cellular Resolution

This protocol describes a reliable method for anesthetization and imaging of intact Drosophila melanogaster larvae. We have utilized the volatile anesthetic desflurane to allow for repetitive imaging at sub-cellular resolution and re-identification of structures for up to a few days1.

Recent improvements in optical imaging, genetically encoded fluorophores and genetic tools allowing efficient establishment of desired transgenic animal ...

Determination of DNA Methylation of Imprinted Genes in Arabidopsis Endosperm

Arabidopsis thaliana is an excellent model organism for studying epigenetic mechanisms. One of the reasons is the loss-of-function null mutant of DNA methyltransferases is viable, thus providing a system to study how loss of DNA methylation in a genome affects growth and development. Imprinting refers to differential expression of maternal and paternal alleles and plays an important role in reproduction development in both mammal and plants. DNA methylation is critical for determining whether the maternal or paternal alleles of an imprinted gene is expressed or silenced. In flowering plants, there is a double fertilization event in reproduction: one sperm cell fertilizes the egg cell to form embryo and a second sperm fuses with the central cell to give rise to endosperm. Endosperm is the tissue where imprinting occurs in plants. MEDEA, a SET domain Polycomb group gene, and FWA, a transcription factor regulating flowering, are the first two genes shown to be imprinted in endosperm and their expression is controlled by DNA methylation and demethylation in plants. In order to determine imprinting status of a gene and methylation pattern in endosperm, we need to be able to isolate endosperm first. Since seed is tiny in Arabidopsis, it remains challenging to isolate Arabidopsis endosperm and examine its methylation. In this video protocol, we report how to conduct a genetic cross, to isolate endosperm tissue from seeds, and to determine the methylation status by bisulfite sequencing.

Determination of DNA Methylation of Imprinted Genes in Arabidopsis Endosperm

Imprinting is a phenomenon in plant and mammal reproduction. DNA methylation plays an important role in mechanisms of imprinting. Isolating endosperm and determining methylation status of imprinted genes in Arabidopsis can be difficult. In this protocol, we describe how to isolate endosperm and determine methylation by bisulfite sequencing.

Arabidopsis thaliana is an excellent model organism for studying epigenetic mechanisms. One of the reasons is the loss-of-function null mutant of DNA ...

DNA Methylation: Bisulphite Modification and Analysis

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene regulation is complex, but essentially involves modifications to the DNA itself or the proteins with which DNA associates. The predominant epigenetic modification of DNA in mammalian genomes is methylation of cytosine nucleotides (5-MeC). DNA methylation provides instruction to gene expression machinery as to where and when the gene should be expressed. The primary target sequence for DNA methylation in mammals is 5'-CpG-3' dinucleotides (Figure 1). CpG dinucleotides are not uniformly distributed throughout the genome, but are concentrated in regions of repetitive genomic sequences and CpG "islands" commonly associated with gene promoters (Figure 1). DNA methylation patterns are established early in development, modulated during tissue specific differentiation and disrupted in many disease states including cancer. To understand the biological role of DNA methylation and its role in human disease, precise, efficient and reproducible methods are required to detect and quantify individual 5-MeCs.
This protocol for bisulphite conversion is the "gold standard" for DNA methylation analysis and facilitates identification and quantification of DNA methylation at single nucleotide resolution. The chemistry of cytosine deamination by sodium bisulphite involves three steps (Figure 2). (1) Sulphonation: The addition of bisulphite to the 5-6 double bond of cytosine (2) Hydrolic Deamination: hydrolytic deamination of the resulting cytosine-bisulphite derivative to give a uracil-bisulphite derivative (3) Alkali Desulphonation: Removal of the sulphonate group by an alkali treatment, to give uracil. Bisulphite preferentially deaminates cytosine to uracil in single stranded DNA, whereas 5-MeC, is refractory to bisulphite-mediated deamination. Upon PCR amplification, uracil is amplified as thymine while 5-MeC residues remain as cytosines, allowing methylated CpGs to be distinguished from unmethylated CpGs by presence of a cytosine "C" versus thymine "T" residue during sequencing.
DNA modification by bisulphite conversion is a well-established protocol that can be exploited for many methods of DNA methylation analysis. Since the detection of 5-MeC by bisulphite conversion was first demonstrated by Frommer et al.1 and Clark et al.2, methods based around bisulphite conversion of genomic DNA account for the majority of new data on DNA methylation. Different methods of post PCR analysis may be utilized, depending on the degree of specificity and resolution of methylation required. Cloning and sequencing is still the most readily available method that can give single nucleotide resolution for methylation across the DNA molecule.

The gold standard for DNA methylation analysis is genomic sequencing of bisulphite converted DNA. This method takes advantage of the increased sensitivity of cytosine compared with 5-methylcytosine (5-MeC) to bisulphite deamination under acidic conditions. Unmethylated cytosines can be distinguished from methylated cytosines after PCR amplification of the target genomic DNA.

Epigenetics describes the heritable changes in gene function that occur independently to the DNA sequence. The molecular basis of epigenetic gene ...

Reduced-gravity Environment Hardware Demonstrations of a Prototype Miniaturized Flow Cytometer and Companion Microfluidic Mixing Technology

Until recently, astronaut blood samples were collected in-flight, transported to earth on the Space Shuttle, and analyzed in terrestrial laboratories. If humans are to travel beyond low Earth orbit, a transition towards space-ready, point-of-care (POC) testing is required. Such testing needs to be comprehensive, easy to perform in a reduced-gravity environment, and unaffected by the stresses of launch and spaceflight. Countless POC devices have been developed to mimic laboratory scale counterparts, but most have narrow applications and few have demonstrable use in an in-flight, reduced-gravity environment. In fact, demonstrations of biomedical diagnostics in reduced gravity are limited altogether, making component choice and certain logistical challenges difficult to approach when seeking to test new technology. To help fill the void, we are presenting a modular method for the construction and operation of a prototype blood diagnostic device and its associated parabolic flight test rig that meet the standards for flight-testing onboard a parabolic flight, reduced-gravity aircraft. The method first focuses on rig assembly for in-flight, reduced-gravity testing of a flow cytometer and a companion microfluidic mixing chip. Components are adaptable to other designs and some custom components, such as a microvolume sample loader and the micromixer may be of particular interest. The method then shifts focus to flight preparation, by offering guidelines and suggestions to prepare for a successful flight test with regard to user training, development of a standard operating procedure (SOP), and other issues. Finally, in-flight experimental procedures specific to our demonstrations are described.

Spaceflight blood diagnostics need innovation. Few demonstrations have been published illustrating in-flight, reduced-gravity health diagnostic technology. Here we present a method for construction and operation of a parabolic flight test rig for a prototype point-of-care flow-cytometry design, with components and preparation strategies adaptable to other setups.

Until recently, astronaut blood samples were collected in-flight, transported to earth on the Space Shuttle, and analyzed in terrestrial laboratories. If ...

In vitro Methylation Assay to Study Protein Arginine Methylation

Protein arginine methylation is one of the most abundant post-translational modifications in the nucleus. Protein arginine methylation can be identified and/or determined via proteomic approaches, and/or immunoblotting with methyl-arginine specific antibodies. However, these techniques sometimes can be misleading and often provide false positive results. Most importantly, these techniques cannot provide direct evidence in support of the PRMT substrate specificity. In vitro methylation assays, on the other hand, are useful biochemical assays, which are sensitive, and consistently reveal if the identified proteins are indeed PRMT substrates. A typical in vitro methylation assay includes purified, active PRMTs, purified substrate and a radioisotope labeled methyl donor (S-adenosyl-L-[methyl-3H] methionine). Here we describe a step-by-step protocol to isolate catalytically active PRMT1, a ubiquitously expressed PRMT family member. The methyl transferase activities of the purified PRMT1 were later tested on Ras-GTPase activating protein binding protein 1 (G3BP1), a known PRMT substrate, in the presence of S-adenosyl-L-[methyl-3H] methionine as the methyl donor. This protocol can be employed not only for establishing the methylation status of novel physiological PRMT1 substrates, but also for understanding the basic mechanism of protein arginine methylation.

In vitro Methylation Assay to Study Protein Arginine Methylation

Protein arginine methylation, catalyzed by a class of enzymes viz., protein arginine methyl transferases (PRMTs), is the process of enzymatic addition of methyl group(s) to arginines within proteins. The in vitro methylation assay is the most dependable tool for assessing the methylation status of known or novel PRMT substrates.

Protein arginine methylation is one of the most abundant post-translational modifications in the nucleus. Protein arginine methylation can be identified ...

Correlating Gene-specific DNA Methylation Changes with Expression and Transcriptional Activity of Astrocytic KCNJ10 (Kir4.1)

DNA methylation serves to regulate gene expression through the covalent attachment of a methyl group onto the C5 position of a cytosine in a cytosine-guanine dinucleotide. While DNA methylation provides long-lasting and stable changes in gene expression, patterns and levels of DNA methylation are also subject to change based on a variety of signals and stimuli. As such, DNA methylation functions as a powerful and dynamic regulator of gene expression. The study of neuroepigenetics has revealed a variety of physiological and pathological states that are associated with both global and gene-specific changes in DNA methylation. Specifically, striking correlations between changes in gene expression and DNA methylation exist in neuropsychiatric and neurodegenerative disorders, during synaptic plasticity, and following CNS injury. However, as the field of neuroepigenetics continues to expand its understanding of the role of DNA methylation in CNS physiology, delineating causal relationships in regards to changes in gene expression and DNA methylation are essential. Moreover, in regards to the larger field of neuroscience, the presence of vast region and cell-specific differences requires techniques that address these variances when studying the transcriptome, proteome, and epigenome. Here we describe FACS sorting of cortical astrocytes that allows for subsequent examination of a both RNA transcription and DNA methylation. Furthermore, we detail a technique to examine DNA methylation, methylation sensitive high resolution melt analysis (MS-HRMA) as well as a luciferase promoter assay. Through the use of these combined techniques one is able to not only explore correlative changes between DNA methylation and gene expression, but also directly assess if changes in the DNA methylation status of a given gene region are sufficient to affect transcriptional activity.

Correlating Gene-specific DNA Methylation Changes with Expression and Transcriptional Activity of Astrocytic KCNJ10 Kir4.1

DNA methylation is capable of maintaining stable levels of gene expression as well as allowing for dynamic changes in gene expression in response to a variety of stimuli. We detail techniques that allow the study of gene-specific changes in DNA methylation and the effect of these changes on gene expression.

DNA methylation serves to regulate gene expression through the covalent attachment of a methyl group onto the C5 position of a cytosine in a ...

Targeted DNA Methylation Analysis by Next-generation Sequencing

The role of epigenetic processes in the control of gene expression has been known for a number of years. DNA methylation at cytosine residues is of particular interest for epigenetic studies as it has been demonstrated to be both a long lasting and a dynamic regulator of gene expression. Efforts to examine epigenetic changes in health and disease have been hindered by the lack of high-throughput, quantitatively accurate methods. With the advent and popularization of next-generation sequencing (NGS) technologies, these tools are now being applied to epigenomics in addition to existing genomic and transcriptomic methodologies. For epigenetic investigations of cytosine methylation where regions of interest, such as specific gene promoters or CpG islands, have been identified and there is a need to examine significant numbers of samples with high quantitative accuracy, we have developed a method called Bisulfite Amplicon Sequencing (BSAS). This method combines bisulfite conversion with targeted amplification of regions of interest, transposome-mediated library construction and benchtop NGS. BSAS offers a rapid and efficient method for analysis of up to 10 kb of targeted regions in up to 96 samples at a time that can be performed by most research groups with basic molecular biology skills. The results provide absolute quantitation of cytosine methylation with base specificity. BSAS can be applied to any genomic region from any DNA source. This method is useful for hypothesis testing studies of target regions of interest as well as confirmation of regions identified in genome-wide methylation analyses such as whole genome bisulfite sequencing, reduced representation bisulfite sequencing, and methylated DNA immunoprecipitation sequencing.

Bisulfite amplicon sequencing (BSAS) is a method for quantifying cytosine methylation in targeted genomic regions of interest. This method uses bisulfite conversion paired with PCR amplification of target regions prior to next-generation sequencing to produce absolute quantitation of DNA methylation at a base-specific level.

The role of epigenetic processes in the control of gene expression has been known for a number of years. DNA methylation at cytosine residues is of ...

Amplification, Next-generation Sequencing, and Genomic DNA Mapping of Retroviral Integration Sites

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of diverse libraries of retroviral integration sites using ligation-mediated PCR (LM-PCR) amplification and next-generation sequencing (NGS), (2) mapping the genomic location of each virus-host junction using BEDTools, and (3) analyzing the data for statistical relevance. Genomic DNA extracted from infected cells is fragmented by digestion with restriction enzymes or by sonication. After suitable DNA end-repair, double-stranded linkers are ligated onto the DNA ends, and semi-nested PCR is conducted using primers complementary to both the long terminal repeat (LTR) end of the virus and the ligated linker DNA. The PCR primers carry sequences required for DNA clustering during NGS, negating the requirement for separate adapter ligation. Quality control (QC) is conducted to assess DNA fragment size distribution and adapter DNA incorporation prior to NGS. Sequence output files are filtered for LTR-containing reads, and the sequences defining the LTR and the linker are cropped away. Trimmed host cell sequences are mapped to a reference genome using BLAT and are filtered for minimally 97% identity to a unique point in the reference genome. Unique integration sites are scrutinized for adjacent nucleotide (nt) sequence and distribution relative to various genomic features. Using this protocol, integration site libraries of high complexity can be constructed from genomic DNA in three days. The entire protocol that encompasses exogenous viral infection of susceptible tissue culture cells to integration site analysis can therefore be conducted in approximately one to two weeks. Recent applications of this technology pertain to longitudinal analysis of integration sites from HIV-infected patients.

We describe a protocol for amplifying retroviral integration sites from the genomic DNA of infected cells, sequencing the amplified virus-host junctions, and then mapping these sequences to a reference genome. We also describe techniques to quantify the distribution of integration sites relative to various genomic annotations using BEDTools.

Retroviruses exhibit signature integration preferences on both the local and global scales. Here, we present a detailed protocol for (1) generation of ...

The ChIP-exo Method: Identifying Protein-DNA Interactions with Near Base Pair Precision

Chromatin immunoprecipitation (ChIP) is an indispensable tool in the fields of epigenetics and gene regulation that isolates specific protein-DNA interactions. ChIP coupled to high throughput sequencing (ChIP-seq) is commonly used to determine the genomic location of proteins that interact with chromatin. However, ChIP-seq is hampered by relatively low mapping resolution of several hundred base pairs and high background signal. The ChIP-exo method is a refined version of ChIP-seq that substantially improves upon both resolution and noise. The key distinction of the ChIP-exo methodology is the incorporation of lambda exonuclease digestion in the library preparation workflow to effectively footprint the left and right 5&#39; DNA borders of the protein-DNA crosslink site. The ChIP-exo libraries are then subjected to high throughput sequencing. The resulting data can be leveraged to provide unique and ultra-high resolution insights into the functional organization of the genome. Here, we describe the ChIP-exo method that we have optimized and streamlined for mammalian systems and next-generation sequencing-by-synthesis platform.

Here, we present a protocol to achieve near base pair resolution of protein-DNA interactions. This is obtained by exonuclease treatment of DNA fragments selectively enriched by chromatin immunoprecipitation (ChIP-exo) followed by high throughput sequencing.

Chromatin immunoprecipitation (ChIP) is an indispensable tool in the fields of epigenetics and gene regulation that isolates specific protein-DNA ...

Generating Transposon Insertion Libraries in Gram-Negative Bacteria for High-Throughput Sequencing

Transposon sequencing (Tn-seq) is a powerful method that combines transposon mutagenesis and massive parallel sequencing to identify genes and pathways that contribute to bacterial fitness under a wide range of environmental conditions. Tn-seq applications are extensive and have not only enabled examination of genotype-phenotype relationships at an organism level but also at the population, community and systems levels. Gram-negative bacteria are highly associated with antimicrobial resistance phenotypes, which has increased incidents of antibiotic treatment failure. Antimicrobial resistance is defined as bacterial growth in the presence of otherwise lethal antibiotics. The &#8220;last-line&#8221; antimicrobial colistin is used to treat Gram-negative bacterial infections. However, several Gram-negative pathogens, including Acinetobacter baumannii can develop colistin resistance through a range of molecular mechanisms, some of which were characterized using Tn-seq. Furthermore, signal transduction pathways that regulate colistin resistance vary within Gram-negative bacteria. Here we propose an efficient method of transposon mutagenesis in A. baumannii that streamlines generation of a saturating transposon insertion library and amplicon library construction by eliminating the need for restriction enzymes, adapter ligation, and gel purification. The methods described herein will enable in-depth analysis of molecular determinants that contribute to A. baumannii fitness when challenged with colistin. The protocol is also applicable to other Gram-negative ESKAPE pathogens, which are primarily associated with drug resistant hospital-acquired infections.

We describe a method to generate saturating transposon mutant libraries in Gram-negative bacteria and subsequent preparation of DNA amplicon libraries for high-throughput sequencing. As an example, we focus on the ESKAPE pathogen, Acinetobacter baumannii, but this protocol is amenable to a wide range of Gram-negative organisms.

Transposon sequencing (Tn-seq) is a powerful method that combines transposon mutagenesis and massive parallel sequencing to identify genes and pathways ...