The authors thank Colleen Van Kirk for mouse retina and cerebellum samples, and Peter Gregory for figure generation. The authors also thank Dr. Allison Gillaspy and the Laboratory for Molecular Biology and Cytometry Research core for access to the MiSeq. This work was supported by NIH grants DA029405, EY021716, and AG026607 and the Donald W. Reynolds Foundation.

1. Nucleic Acid Isolation from Tissue
NOTE: Co-isolation of DNA and RNA from the same tissue sample offers the opportunity to collect DNA for epigenetic analysis and paired RNA for gene expression analysis. The example given here is a form of column purification from experimental tissue but alternative approaches for different sample types or other isolation methods can be employed. The primary requirement is for highly purified nucleic acid without contaminating protein or organic solvents.
<ol>
	<li>Select a piece of fresh-frozen or fresh tissue for mechanical homogenization and place on ice. Store or transfer tissue to a 2.0 ml round-bottom microcentrifuge tube.</li>
	<li>Add one 5 mm stainless steel bead to each tube containing tissue. Add an appropriate amount of lysis buffer, based on manufacturer&#8217;s suggestions of cell number or tissue weight.</li>
	<li>Homogenize the tissue using a bead mill mechanical homogenization at 30 Hz for 30 sec. 
	NOTE: The frequency and duration of mechanical homogenization is dependent on the tissue type. Softer tissues will completely homogenize using the recommended settings, while tougher tissues require setting optimization.</li>
	<li>Carry out the DNA and RNA isolation using silica spin columns at room temperature following manufacturer&#8217;s protocol.
	<ol>
		<li>Elute RNA in 50 &#181;l of RNase-free water and place on ice. Use RNA for mRNA expression analysis by qPCR.</li>
		<li>Elute DNA in 100 &#181;l of TE buffer. Repeat elution off DNA column using the eluent to increase DNA concentration and yield. Use DNA for targeted methylation quantitation.</li>
	</ol>
	</li>
	<li>Assess the quality of DNA and RNA using a standard spectrophotometer for absorbance at 230 nm, 260 nm, 280 nm, and 320 nm. 
	NOTE: A260/A280 ratio should be ~1.8-2 for high quality DNA. A260/A280 ratio should be ~2 for high quality RNA. The presence of excess absorbance at 230 nm indicates contaminants from the isolation procedure.</li>
	<li>Quality check the RNA further by using a capillary electrophoresis chip according to manufacturer&#8217;s instructions. 
	NOTE: High quality RNA will have a RNA integrity number (RIN) &#62;8. The presence of extra peaks in the electropherogram besides rRNA peaks indicate sample contamination.</li>
	<li>Quantify DNA and RNA by fluorescent assay according to manufacturer&#8217;s protocol. 
	NOTE: Using a fluorometric assay is more sensitive and specific for quantifying nucleic acids than spectrophotometry alone.</li>
	<li>Store isolated DNA and RNA at -80 &#176;C.</li>
</ol>
2. Target Identification and Primer Design
<ol>
	<li>Determine the genomic region(s) of interest to be analyzed by BSAS based on the appropriate reference genome sequence. Gene promoters of differentially expressed genes are common targets of methylation quantitation.</li>
	<li>Prepare an in silico bisulfite-converted reference genome for later alignment of sequencing reads by altering the .fasta sequence file in a text editor. In the 5&#8217;-3&#8217; orientation, replace non-CpG cytosine&#8217;s with thymine&#8217;s (Figure 2).</li>
	<li>Design primer sets to amplify regions of interest from bisulfite converted DNA.
	<ol>
		<li>Select and copy the unconverted region of interest, in the 5&#8217;-3&#8217; orientation, into a bisulfite-specific PCR primer design program (Figure 3). 
		NOTE: An optimal bisulfite-PCR amplicon length is 250-400 bp per amplicon as bisulfite treatment fragments DNA and it is difficult to amplify large &#62;400 bp regions. If, for example, a region of 1 kb is of interest, multiple primer pairs can be designed to cover this region. Design amplicons to be of equal bp size for easy pooling. Avoid amplicon lengths &#60;250 bp because these may be of insufficient length for library generation. An optimal primer length for BSAS is &#62; 20 bp per primer.</li>
		<li>Select a Tm range of 55-65 &#176;C and a max Tm difference between forward and reverse primers of 1-2 &#176;C.</li>
		<li>Select primer pairs that best cover the region of interest. Do not select primers that contain CpG sites or are directly adjacent to CpG sites, as this will cause bias in the PCR reactions. Use standard PCR primers that can be ordered from any number of academic or commercial providers.</li>
		<li>Reconstitute lyophilized primers in RNase/DNase free water at a working stock of 100 &#181;M.</li>
	</ol>
	</li>
	<li>Store PCR primers at -20 &#176;C. Dilute 100 &#181;M primer stock to 10 &#181;M working stock for PCR reactions.</li>
</ol>
3. Bisulfite Specific PCR Optimization
NOTE: For bisulfite conversion of genomic DNA, a number of different commercial kits for bisulfite conversion are available. Select the kit or protocol that best suits the planned experiment.
<ol>
	<li>Use between 200 ng to 2 &#181;g of genomic DNA. For optimization experiments, run multiple conversion reactions in order to have sufficient bisulfite converted DNA for multiple BS PCR reactions.</li>
	<li>Use small (e.g., 10 &#181;l) elution volumes of bisulfite converted DNA to maintain high DNA concentrations. 
	NOTE: 1-2 &#181;l of bisulfite converted DNA is sufficient for optimization PCR reactions if 1 &#181;g genomic DNA was used for conversion input.</li>
	<li>Amplify bisulfite converted DNA by bisulfite specific PCR. BS-PCR reactions require using a Taq polymerase capable of amplifying bisulfite converted DNA.
	<ol>
		<li>Assemble the following reaction for target amplification optimization of a single amplicon. For multiple samples, assemble reactions in a 96-well PCR plate. 
		25 &#181;l 2x reaction buffer 
		0.5 &#181;l dNTP Mix 
		5 &#181;l 10 &#181;M Forward Primer (Final 1 &#181;M) 
		5 &#181;l 10 &#181;M Reverse Primer (Final 1 &#181;M) 
		2 &#181;l template bisulfite converted DNA 
		0.4 &#181;l (5 U/&#181;l) DNA polymerase 
		12.1 &#181;l DNase free water (Final reaction volume 50 &#181;l)</li>
		<li>Seal PCR plate with the appropriate adhesive or heat-sealed film.</li>
		<li>Place reaction in an appropriate thermal cycler and use the following cycling conditions with a heated lid:
		<ol>
			<li>Perform an initial denaturation at 95 &#176;C for 10 min.</li>
			<li>Denature at 95 &#176;C for 30 sec.</li>
			<li>Anneal for 30 sec at the specific Tm of primers being used. Start at an annealing temperature a few &#176;C below the Tm of the primer for optimization reactions. 
			NOTE: Best annealing temperatures can also be determined using a gradient thermal cycler to test a range of temperatures.</li>
			<li>Perform an extension at 72 &#176;C for 30 sec. Lengthen the extension time for longer amplicons.</li>
			<li>Repeat steps 3.3.3.2-3.3.3.4 for 35 cycles total for initial optimization. Higher cycle numbers might be necessary but should generally be avoided to prevent clonal amplifications that will create reaction artifacts.</li>
			<li>Perform a final extension at 72 &#176;C for 7 min.</li>
			<li>Hold reactions at 4 &#176;C.</li>
		</ol>
		</li>
		<li>Visualize amplicons by PAGE electrophoresis. Alternatively, use a capillary electrophoresis DNA chip according to manufacturer&#8217;s protocol.
		<ol>
			<li>Mix 24 &#181;l of the PCR reaction with 6 &#181;l of 5x loading dye on ice and vortex.</li>
			<li>Make 1 L of 1x TBE running buffer by mixing 200 ml of 5x TBE running buffer and 800 ml of ddH2O. Dilute DNA ladder for a final concentration of 0.5 &#181;g per well in ddH2O and loading dye.</li>
			<li>Load wells with 30 &#181;l of PCR reactions or ladder.</li>
			<li>Run the gel at 200 V for ~45 min. Stop the gel run when the first dye front reaches the end of the gel.</li>
			<li>Stain the gel with 5 &#181;g/ml ethidium bromide by mixing 5 &#181;l of 10 mg/ml in 100 ml of ddH2O. Cover the entire gel and let incubate for ~5 min.</li>
			<li>Image the gel via a UV light box or multimodal imager at an excitation wavelength of 482 nm.</li>
			<li>In reference to the ladder, size the PCR amplicons and determine specificity of the reaction by looking for multiple bands other than the expected amplicon size (Figure 4A and C).</li>
		</ol>
		</li>
		<li>Once optimal conditions for BS-PCR have been determined, perform BS-PCR on experimental samples.
		<ol>
			<li>Purify amplicons from remaining primers and other reaction components. Perform silica column purification, SPRI bead, or gel-based clean-up and quantification.</li>
			<li>Elute amplicons in 30 &#181;l TE buffer.</li>
			<li>Quantify amplicons using fluorometric assay using manufacturer&#8217;s protocols.</li>
			<li>If there are multiple amplicons per sample, pool them at equal weight per volume amounts and store at -20 &#176;C.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
4. NGS Library Preparation and Library QC
NOTE: BSAS NGS library preparation uses a simplified transposome-mediated protocol with dual-indexing (See Materials List for details on library preparation kit selection). This method provides an extremely rapid and high-throughput route for library construction that has been validated and shows little to no bias in bisulfite sequencing applications16,19. Dual-indexing allows for a higher multiplexing of samples than single indexing. Other library preparation styles may be possible but have not been tested.
<ol>
	<li>Prepare NGS libraries from BS PCR amplicons in a 96-well PCR plate. Dilute amplicons to 0.2 ng/&#181;l for a final input mass of 1 ng, or 5 &#181;l of the 0.2 ng/&#181;l dilution. Perform a clean-up of generated libraries using SPRI beads. Use 50-90 &#181;l SPRI beads for isolating libraries with this protocol. The volume of SPRI beads depends on target size of library and cycle number of desired sequencing reactions. 
	NOTE: Use a microplate shaker or plate vortexer to resuspend beads.
	<ol>
		<li>Elute the libraries off the beads and transfer to a new 96-well plate and seal with heat-sealing film for storage at -20 &#176;C.</li>
	</ol>
	</li>
	<li>Determine the size and concentration of libraries (Figure 4B and D).
	<ol>
		<li>Run libraries on a capillary electrophoresis DNA sizing chip according to manufacturer&#8217;s protocols.
		<ol>
			<li>Use a high sensitivity capillary electrophoresis DNA assay to determine the average size of detected library peaks and an estimation of molarity according to manufacturer&#8217;s protocols.</li>
		</ol>
		</li>
		<li>Quantify libraries by qPCR, using primers designed against the adapter sequences in the generated libraries and use standards with a known concentration.
		<ol>
			<li>Run qPCR reactions on a real-time instrument using the absolute quantitation settings according to manufacturer&#8217;s protocols.</li>
			<li>Use the size of libraries and the molar quantitation from qPCR to calculate the molar concentrations of the libraries.</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>Dilute libraries to 4 nM using 10 mM Tris-Cl with 0.1% Tween-20 at pH 8.5. Store libraries in a 96-well PCR plate sealed with heat-sealing film at -20 &#176;C.</li>
</ol>
5. Sequencing Libraries Using a Benchtop Sequencer
<ol>
	<li>Thaw out and prepare reagents according to manufacturer&#8217;s instructions.</li>
	<li>Thaw 4 nM libraries on ice.</li>
	<li>Make 1 ml of 0.2 N NaOH in a 1.5 ml microcentrifuge tube.</li>
	<li>Pool 4 nM libraries in equal volume amounts, if multiple libraries are to be used.</li>
	<li>Dilute and denature pooled libraries according to manufacturer&#8217;s instructions to the desired final molar concentration.
	<ol>
		<li>Keep diluted and denature pooled library on ice until ready to load on reagent cartridge.</li>
	</ol>
	</li>
	<li>Generate a sample sheet containing sample names and corresponding index information and specify to generate .fastq files only.
	<ol>
		<li>Load sample sheet in the corresponding sample sheet folder on the sequencer.</li>
	</ol>
	</li>
	<li>Add a total of 600 &#181;l diluted and denatured library to the reagent cartridge in the corresponding well.</li>
	<li>Run the sequencing reactions to completion.</li>
</ol>
6. Methylation Quantitation Data Analysis
NOTE: There are multiple programs available for NGS data analysis including both commercial and open source. Details on running programs can be found with each specific package. General instructions are given at each step along with the specific commands for this software package. (See Materials List for details on software used in this protocol).
<ol>
	<li>Import .fastq.gz files into appropriate sequencing analysis pipeline retaining quality scores (Qscores) for use in read trimming. (For this example program select &#8216;Import&#8217; and select the sequencing method used to generate reads. Select .fastq.gz read files to import, and retain Qscores from .fastq files for read trimming applications by deselecting &#8216;Discard quality scores&#8217; under &#8216;General Options&#8217;. Select a location for the imported reads and select &#8216;Finish&#8217;.)
	<ol>
		<li>Trim and retain reads using the following settings in a NGS data analysis pipeline. (For this program select &#8216;Trim Sequences&#8217; under &#8216;NGS Core Tools&#8217;, and select the reads to be trimmed.)
		<ol>
			<li>Retain only reads containing &#8805; Q30 scores. (For this program select &#8216;Trim using quality scores&#8217; and set the limit to 0.001.)</li>
			<li>Retain only reads containing &#8804; 1 ambiguous nucleotide. (For this program select &#8216;Trim ambiguous nucleotides&#8217; and set the &#8216;maximum number of ambiguities&#8217; to 1. Select a location to save trimmed reads and select &#8216;Finish&#8217;.)</li>
		</ol>
		</li>
		<li>Map trimmed reads to the in silico bisulfite converted gene reference from 2.2. (For the example software select &#8216;Map Reads to Reference&#8217; under &#8216;NGS Core Tools&#8217;. Select trimmed reads to be mapped and select &#8216;next&#8217;. Select the converted reference sequence and select &#8216;no masking&#8217; under &#8216;reference masking&#8217;. Select &#8216;next&#8217;.)
		<ol>
			<li>Assign a maximum penalty for mismatches, insertions, and deletions. Set parameters such that 100% of the read length is required to map with 90% sequence identity. For example software, assign a score of 3 for mismatches, insertions, and deletions. Assign a score of 1 for minimum fraction of read length mapped to reference, and assign a score of 0.9 for minimum fraction of identity between mapped read and reference.</li>
			<li>Discard or ignore non-specific reads mapped to reference. These reads can map equally well to multiple regions. For example software select &#8216;Ignore&#8217; under &#8216;Non-specific match handling&#8217; and select &#8216;next&#8217;.</li>
			<li>For alignment output generate each sample as an independent file. For example software, select &#8216;Create stand-alone read mappings&#8217; under &#8216;Output options&#8217; and &#8216;Save&#8217; under &#8216;Result handling&#8217;. Select a location to save the read mapping and select &#8216;Finish&#8217;.</li>
		</ol>
		</li>
		<li>Run a variant calling workflow on the mapped reads. For the example software, select &#8216;Probabilistic Variant Detection&#8217; under &#8216;Resequencing Analysis&#8217;, select the read mapping to be analyzed and select &#8216;next&#8217;. Select &#8216;Ignore non-specific matches&#8217; under &#8216;Read filters&#8217; and select &#8216;next&#8217;.
		<ol>
			<li>Ensure that the sequence is at a depth of &#8805; 1000 X by setting the &#8216;minimum coverage&#8217; under &#8216;Significance&#8217; to 1,000.</li>
			<li>Ensure that the variant call is present in both forward and reverse reads. For example software, select &#8216;Require presence in both forward and reverse stands&#8217; under &#8216;Variant filters&#8217; and select &#8216;next&#8217;.</li>
			<li>Export aligned, variant called data as an annotated table. For example software, select &#8216;Create annotated table&#8217; under &#8216;Output options&#8217; and select a location to save the variant table file. Select &#8216;Finish&#8217;.</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>Calculate 5-mC percentage using the variant frequencies at annotated CpG sites in the region of interest from the variant table file. The frequency of cytosine at a reference C in CG is the 5-mC percentage. 
	NOTE: When using variant calling algorithms for highly methylated CpG&#8217;s, replace CpG cytosines with thymines. This will identify the mapped cytosine&#8217;s as variants, resulting in cytosine frequency at CpG&#8217;s.</li>
	<li>Determine bisulfite conversion efficiency (bsc) by first calculating the percentage of C&#8217;s not in CpG sites in the reference (C). Multiply the percentage of non CpG C&#8217;s by number of total T&#8217;s mapped (Tmp), and divide by total C&#8217;s mapped at the reference T&#8217;s (CmpT). 100 minus the frequency calculated is the bisulfite conversion efficiency (bsc) of that library (Equation 1).</li>
</ol>
<img alt="Equation 1" fo:content-width="1.5in" src="/files/ftp_upload/52488/52488eq1.jpg" />
Equation 1. Bisulfite conversion efficiency. 
NOTE: While mammalian genomes contain very low amounts of non-CpG cytosine methylation (with the exception potentially of stem cells20), plant genomes contain a significant amount. In order to determine bisulfite conversion efficiency in these reference genomes, a suitable route is to sequence a portion of the mitochondrial genome or chloroplast genome or known unmethylated genes present in plant genomes21 and calculating conversion efficiency as described above. Bisulfite conversion efficiency should be &#8805;98% in most cases, with the unconverted C&#8217;s potentially arising from non CpG methylation.

<ol>
	<li>Hotchkiss, R. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+quantitative+separation+of+purines,+pyrimidines,+and+nucleosides+by+paper+chromatography.">The quantitative separation of purines, pyrimidines, and nucleosides by paper chromatography.</a> The Journal of Biological Chemistry. 175 (1), 315-332 (1948).</li><li>Bird, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DNA+methylation+patterns+and+epigenetic+memory.">DNA methylation patterns and epigenetic memory.</a> Genes &amp; Development. 16 (1), 6-21 (2002).</li><li>Smith, Z. D., Meissner, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DNA+methylation:+roles+in+mammalian+development.">DNA methylation: roles in mammalian development.</a> Nature reviews. Genetics. 14 (3), 204-220 (2013).</li><li>Heard, E., Martienssen, R. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transgenerational+epigenetic+inheritance:+myths+and+mechanisms.">Transgenerational epigenetic inheritance: myths and mechanisms.</a> Cell. 157 (1), 95-109 (2014).</li><li>Baylin, S. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DNA+methylation+and+gene+silencing+in+cancer.">DNA methylation and gene silencing in cancer.</a> Nature Clinical Practice. Oncology. 2, S4-S11 (2005).</li><li>Baylin, S. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+cancer+epigenome:+its+origins,+contributions+to+tumorigenesis,+and+translational+implications.">The cancer epigenome: its origins, contributions to tumorigenesis, and translational implications.</a> Proceedings of the American Thoracic Society. 9 (2), 64-65 (2012).</li><li>Beck, S., Rakyan, V. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+methylome:+approaches+for+global+DNA+methylation+profiling.">The methylome: approaches for global DNA methylation profiling.</a> Trends in Genetics : TIG. 24 (5), 231-237 (2008).</li><li>Laird, P. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Principles+and+challenges+of+genomewide+DNA+methylation+analysis.">Principles and challenges of genomewide DNA methylation analysis.</a> Nature reviews. Genetics. 11 (3), 191-203 (2010).</li><li>Mikeska, T., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Optimization+of+quantitative+MGMT+promoter+methylation+analysis+using+pyrosequencing+and+combined+bisulfite+restriction+analysis.">Optimization of quantitative MGMT promoter methylation analysis using pyrosequencing and combined bisulfite restriction analysis.</a> The Journal of Molecular Diagnostics : JMD. 9 (3), 368-381 (2007).</li><li>Kreutz, M., Hochstein, N., Kaiser, J., Narz, F., Peist, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Pyrosequencing:+powerful+and+quantitative+sequencing+technology.+Current+Protocols+In.">Pyrosequencing: powerful and quantitative sequencing technology. Current Protocols In.</a> Molecular Biology / edited by Frederick M. Ausubel ... [et al.]. 104 (Unit 7 15), (2013).</li><li>Dikow, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantification+of+the+methylation+status+of+the+PWS/AS+imprinted+region:+comparison+of+two+approaches+based+on+bisulfite+sequencing+and+methylation-sensitive+MLPA.">Quantification of the methylation status of the PWS/AS imprinted region: comparison of two approaches based on bisulfite sequencing and methylation-sensitive MLPA.</a> Molecular And Cellular Probes. 21 (3), 208-215 (2007).</li><li>Parrish, R. R., Day, J. J., Lubin, F. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Direct+bisulfite+sequencing+for+examination+of+DNA+methylation+with+gene+and+nucleotide+resolution+from+brain+tissues.">Direct bisulfite sequencing for examination of DNA methylation with gene and nucleotide resolution from brain tissues.</a> Current Protocols In Neuroscience / editorial board, Jacqueline N. Crawley ... [et al.]. 7 (Unit 7 24), (2012).</li><li>Shapiro, R., Servis, R. E., Welcher, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Reactions+of+uracil+and+cytosine+derivatives+with+sodium+bisulfite.+A+specific+deamination+method.">Reactions of uracil and cytosine derivatives with sodium bisulfite. A specific deamination method.</a> Journal of the American Chemical Society. 92, 422-424 (1970).</li><li>Hayatsu, H., Wataya, Y., Kazushige, K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+addition+of+sodium+bisulfite+to+uracil+and+to+cytosine.">The addition of sodium bisulfite to uracil and to cytosine.</a> Journal of the American Chemical Society. 92 (3), 724-726 (1970).</li><li>Wang, R. Y., Gehrke, C. W., Ehrlich, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comparison+of+bisulfite+modification+of+5-methyldeoxycytidine+and+deoxycytidine+residues.">Comparison of bisulfite modification of 5-methyldeoxycytidine and deoxycytidine residues.</a> Nucleic Acids Research. 8 (20), 4777-4790 (1980).</li><li>Masser, D. R., Berg, A. S., Focused Freeman, W. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=high+accuracy+5-methylcytosine+quantitation+with+base+resolution+by+benchtop+next-generation+sequencing.">high accuracy 5-methylcytosine quantitation with base resolution by benchtop next-generation sequencing.</a> Epigenetics &amp; Chromatin. 6 (1), 33 (2013).</li><li>Caruccio, N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Preparation+of+next-generation+sequencing+libraries+using+Nextera+technology:+simultaneous+DNA+fragmentation+and+adaptor+tagging+by+in+vitro+transposition.">Preparation of next-generation sequencing libraries using Nextera technology: simultaneous DNA fragmentation and adaptor tagging by in vitro transposition.</a> Methods in Molecular Biology. 733, 241-255 (2011).</li><li>Quail, M. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+tale+of+three+next+generation+sequencing+platforms:+comparison+of+Ion+Torrent,+Pacific+Biosciences+and+Illumina+MiSeq+sequencers.">A tale of three next generation sequencing platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq sequencers.</a> BMC Genomics. 13, 341 (2012).</li><li>Adey, A., Shendure, J. U. l. t. r. a. -. l. o. w. -. i. n. p. u. t. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=tagmentation-based+whole-genome+bisulfite+sequencing.">tagmentation-based whole-genome bisulfite sequencing.</a> Genome Research. 22 (6), 1139-1143 (2012).</li><li>Smallwood, S. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Single-cell+genome-wide+bisulfite+sequencing+for+assessing+epigenetic+heterogeneity.">Single-cell genome-wide bisulfite sequencing for assessing epigenetic heterogeneity.</a> Nature Methods. 11 (8), 817-820 (2014).</li><li>Wang, J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Universal+endogenous+gene+controls+for+bisulphite+conversion+in+analysis+of+plant+DNA+methylation.">Universal endogenous gene controls for bisulphite conversion in analysis of plant DNA methylation.</a> Plant Methods. 7, 39 (2011).</li><li>Lister, R., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Human+DNA+methylomes+at+base+resolution+show+widespread+epigenomic+differences.">Human DNA methylomes at base resolution show widespread epigenomic differences.</a> Nature. 462 (7271), 315-322 (2009).</li><li>Ivanov, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=In-solution+hybrid+capture+of+bisulfite-converted+DNA+for+targeted+bisulfite+sequencing+of+174+ADME+genes.">In-solution hybrid capture of bisulfite-converted DNA for targeted bisulfite sequencing of 174 ADME genes.</a> Nucleic Acids Research. 41 (6), e72 (2013).</li><li>Booth, M. J., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Oxidative+bisulfite+sequencing+of+5-methylcytosine+and+5-hydroxymethylcytosine.">Oxidative bisulfite sequencing of 5-methylcytosine and 5-hydroxymethylcytosine.</a> Nature Protocols. 8 (10), 1841-1851 (2013).</li></ol>

The authors have nothing to disclose.

BSAS allows for focused, accurate DNA methylation quantitation with high throughput capabilities in both number of samples and number of targets. Additionally, this library generation using transposome-mediated tagmentation requires less input DNA, is quicker and contains fewer steps than traditional shearing and subsequent adapter ligation techniques. These improvements over existing methods allow for precise base-level DNA methylation analysis of any target of interest. Moreover, BSAS provides a new method for confirmation of regions of interest (differentially methylation regions (DMRs), CpG islands, regulatory regions) that have been identified using whole genome bisulfite22 or capture23 approaches. Using orthogonal confirmation approaches and more quantitatively precise methods improves the analytical rigor of epigenomic studies. The high read depth achieved with BSAS allows for accurate quantitation in a narrow confidence interval, resulting in precise quantitation16. Additionally, the ability to run many samples in a single experiment can increase the sample size for confirmations of whole genome methods, which will increase the statistical power; a weakness of many current epigenetic studies. Importantly, BSAS provides a base-specific CpG methylation quantitation, which allows for the generation of testable hypotheses for epigenetic research. Hypotheses such as increased methylation at a specific response element will result in decreased binding of a specific transcription factor. BSAS, combined with paired mRNA expression data, provides a method for obtaining both epigenetic regulation and mRNA expression within a single piece of tissue or cell population. Paired measurements of regulation and expression also increase the scientific rigor and impact of the findings.
Critical steps within the BSAS protocol include primer design, amplicon optimization, and library generation/QC. PCR bias is introduced if primers are not designed properly and amplify at different efficiencies8. The use of methylation standards, such as enzymatically generated 100% and 0% cytosine methylated standards, can be used to determine the degree, if any, of methylation quantitation bias, and is a proper methodological control when quantifying methylation16. Similarly, care should be taken when optimizing PCR reactions for generating a single high-quality amplicon. If no amplicon is present using the suggested guidelines detailed in the protocol, optimization steps include increasing PCR cycle numbers or bsDNA input amount. If there are multiple PCR products, including primer-dimers, optimization steps include decreasing bsDNA input amount, decreasing the PCR primer molar input concentrations, increasing annealing temperatures, or decreasing PCR cycle number. One or a combination of these optimization procedures yields sufficient results to generate a single high-quality amplicon. Alternatively, redesigning primers is a good approach for primer sets that do not yield a single amplicon. For those regions where redesign is not suitable or possible, degenerate primer sets including both cytosine&#8217;s and thymine&#8217;s at CpG sites in forward primers and adenine&#8217;s and guanine&#8217;s at CpG sites in reverse primers may work. However, these primer sets need further optimization to ensure no PCR bias is occurring, a process outside the scope of this protocol. Quality libraries are extremely important for success. Ensure that QC measures are carried out to determine size, quality and quantity of libraries before sequencing. Sequencing reactions are prone to failure with input of low-quality libraries. Bias in methylation quantitation can be mitigated by ensuring that methylation frequencies are present in both forward and reverse sequencing reads. Additionally, quality scores of bases aligning to CpG sites should be &#8805;Q30 for high confidence in quantitation. While others have shown previously little to no bias in tagmentation reactions of bisulfite converted DNA19, ensuring the metrics described above allows for confirmation of low bias methylation quantitation.
BSAS currently measures methylation at CpG sites, however, future expansions of this method will allow paired measurements of CpG 5-hmC with the addition of experimental steps distinguishing between 5-mC and 5-hmC such as introducing potassium perruthenate before bisulfite conversion24. This will increase the impact of BSAS utility by allowing for paired methylation, both 5-mC and 5-hmC, and expression data in order to determine gene expression regulatory roles of DNA methylation. Transposome-mediated library preparation of PCR amplicons does result in decreased sequencing depth at the ends of amplified regions. Therefore, regions of high interest should be designed to reside in the middle of amplicons. However, because BSAS generates sequencing read depths of &#62;1,000 X, the quantitative accuracy of 5-mC measurements near the ends of amplified regions remains high.

It has been over a half-century since the first report of naturally occurring DNA modifications in the form of cytosine methylation1. Cytosine methylation is a modified cytosine nucleotide base where the 5th carbon on the base ring has a methyl group (5-mC), most frequently occurring in the CpG dinucleotide motif in mammalian genomes. The functional presence of 5-mC in gene promoters is generally associated with transcriptional repression, while the absence is associated with transcriptional activity2.
Tremendous advances have been made in our understanding of the role of DNA methylation in development3, transgeneration propagation of epigenetic profiles4, cancer pathogenesis5,6, and a number of other research areas. Many of these advances have come in the last few years as new methodologies have been developed to profile DNA methylation7. With the advent and popularization of next-generation sequencing (NGS) techniques there are now a number of methods to profile DNA methylation across a whole genome or large regions of a genome8. These approaches open the possibility of discovery studies that quantify DNA methylation patterns and differences in DNA methylation. In addition to these hypothesis-generating approaches, there is a significant need for targeted/hypothesis testing DNA methylation quantitation methods.
Pyrosequencing, methylation specific PCR, and direct Sanger sequencing of bisulfite converted DNA have been the most used methods for analysis of targeted regions (i.e., a promoter region of a single gene or a CpG island)9-12. All of these methods rely on bisulfite conversion, a chemical deamination reaction where unmethylated cytosine&#8217;s are converted into uracil&#8217;s, while methylated cytosine&#8217;s remain intact13,14. PCR amplification of bisulfite-converted DNA results in replacement of uracil&#8217;s with thymine15 allowing the differential methylation to read out as a base difference. While highly useful, the limitations to these methods include low quantitative accuracy, short read length, and low sample throughput. To address these limitations we developed a method we have termed Bisulfite Amplicon Sequencing (BSAS)16. The goal of this method is to be able to analyze target genomic regions of interest in a large number of samples with high quantitative accuracy.
BSAS uses elements of existing methods (bisulfite conversion and region-specific PCR amplification) and combines them with simple next-generation library construction (transposome-mediated)17 and benchtop NGS18 (Figure 1). This approach provides for a more quantitative and higher throughput method to examine the cytosine methylation in any region of interest. Here we present the method in detail and describe approaches for troubleshooting and quality checking the BSAS process.

<table border="0" cellpadding="0" cellspacing="0" width="840">
	<tbody>
		<tr height="20">
			<td height="20" style="height:20px;width:304px;">Name</td>
			<td style="width:255px;">Company</td>
			<td style="width:97px;">Cat #</td>
			<td style="width:185px;">Comments</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Agilent 2100 Bioanlyzer</td>
			<td>Agilent</td>
			<td>G29393AA</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Agilent RNA 6000 Nano kit</td>
			<td>Agilent</td>
			<td>5067-1511</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Agilent High Sensitivity DNA kit</td>
			<td>Agilent</td>
			<td>5067-4626</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Typhoon 9200</td>
			<td>Amersham/GE Healthcare Life Sciences</td>
			<td></td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Agencourt AMPure XP - PCR Purification</td>
			<td>Beckman Coulter</td>
			<td>A63880</td>
			<td>5 mL</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">PowerPac Basic Power Supply</td>
			<td>Bio Rad</td>
			<td>164-5050</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">twin.tech PCR plate 96&#160;</td>
			<td>Eppendorf</td>
			<td>951020362</td>
			<td>semi-skirted</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Heatsealing Film</td>
			<td>Eppendorf</td>
			<td>0030127838</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Heatsealing Foil</td>
			<td>Eppendorf</td>
			<td>0030127854</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Heat Sealer</td>
			<td>Eppendorf</td>
			<td>5390000.024</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">0.1-10ul Tips</td>
			<td>Eppendorf</td>
			<td>0030073.002</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">2-200ul Tips</td>
			<td>Eppendorf</td>
			<td>0030073.045</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">50-1000ul Tips</td>
			<td>Eppendorf</td>
			<td>0030073.100</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">MixMate&#160;</td>
			<td>Eppendorf</td>
			<td>5353000.014</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">1.5 mL Microcentrifuge Tubes</td>
			<td>Eppendorf</td>
			<td>0030121.023</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Centrifuge 5424</td>
			<td>Eppendorf</td>
			<td>5424000.010</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">2.0 mL Microcentrifuge Tube</td>
			<td>Eppendorf</td>
			<td>022363352</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Nextera XT DNA Sample Preparation Kit</td>
			<td>Illumina</td>
			<td>FC-131-1024</td>
			<td>24 rxn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Nextera XT Index Kit</td>
			<td>Illumina</td>
			<td>FC-131-1001</td>
			<td>24 indexes</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">MiSeq Desktop Sequencer</td>
			<td>Illumina</td>
			<td></td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">MiSeq Reagent Kit v2</td>
			<td>Illumina</td>
			<td>MS-102-2002</td>
			<td>300cycle</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">KAPA SYBR FAST Universal qPCR kit</td>
			<td>KAPA Biosystems</td>
			<td>KK4824</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Quant-iT PicoGreen dsDNA Assay kit</td>
			<td>Life Technologies</td>
			<td>P7589</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Quant-iT RiboGreen RNA Assay kit</td>
			<td>Life Technologies</td>
			<td>R11490</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">ProFlex 96-well PCR system</td>
			<td>Life Technologies</td>
			<td>4484075</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Novex 6% TBE DNA gel</td>
			<td>Life Technologies</td>
			<td>EC6261BOX</td>
			<td>10-well</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Novex TBE Running Buffer</td>
			<td>Life Technologies</td>
			<td>LC6675</td>
			<td>5X (1 L)</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Novex Hi-Density TBE Sample Buffer</td>
			<td>Life Technologies</td>
			<td>LC6678</td>
			<td>5X (10 mL)</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">1 Kb Plus DNA Ladder</td>
			<td>Life Technologies</td>
			<td>10787-018</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Xcell SureLock Mini-Cell</td>
			<td>Life Technologies</td>
			<td>EI0001</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">UltraPure 10mg/mL Ethidium Bromide</td>
			<td>Life Technologies</td>
			<td>15585-011</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">7900HT Fast Real-Time PCR system</td>
			<td>Life Technologies</td>
			<td>4329001</td>
			<td>384well block</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">DynaMag-96 Side Skirted</td>
			<td>Life Technologies</td>
			<td>12027</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">SpectroMax M2 Multi-Mode Microplate Reader</td>
			<td>Molecular Devices/VWR</td>
			<td>89429-532</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">AllPrep DNA/RNA Mini Kit</td>
			<td>Qiagen</td>
			<td>80204</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Stainless Steel Beads</td>
			<td>Qiagen</td>
			<td>69989</td>
			<td>5 mm</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">CLC Genomics Workbench</td>
			<td>Qiagen</td>
			<td></td>
			<td>Software for methylation pipeline</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">QIAQuick PCR Purification kit</td>
			<td>Qiagen</td>
			<td>28104</td>
			<td>50 rxn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">TissueLyser II</td>
			<td>Retsch/Qiagen</td>
			<td>85300</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Nuclease-Free Water</td>
			<td>Sigma</td>
			<td>W4502</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">TRIS hydrochloride</td>
			<td>Sigma</td>
			<td>PHG0002-100G</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Tween-20</td>
			<td>Sigma</td>
			<td>P9416-50ML</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">NanoDrop 2000</td>
			<td>Thermo Scientific</td>
			<td></td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">384-well Full-Skirted Plate, Standard</td>
			<td>Thermo Scientific</td>
			<td>AB-1384</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">Absolute qPCR Seal</td>
			<td>Thermo Scientific</td>
			<td>AB-1170</td>
			<td></td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">EZ DNA Methylation-Lightning kit</td>
			<td>Zymo Research</td>
			<td>D5030</td>
			<td>50 rxn</td>
		</tr>
		<tr height="20">
			<td height="20" style="height:20px;">ZymoTaq DNA Polymerase</td>
			<td>Zymo Research</td>
			<td>E2001</td>
			<td>50 rxn</td>
		</tr>
	</tbody>
</table>

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.

BSAS reads properly aligned to converted reference sequence will resemble Figure 5. CpG dinucleotides can clearly be identified and methylation states can be estimated by observing base calls at CpG sites in the mapped reads. For example, with methylation controls, 0% methylation will result in all the reads mapped to CpG sites containing T&#8217;s (Figure 5A). 100% methylation controls will result in all the reads mapped reads to CpG sites containing C&#8217;s (Figure 5B).
Quantitation of cytosine frequency (C/C+T) at CpG sites yields the methylation frequency in the original sample. A representative standard curve generated from whole genome methylation controls (n = 3/methylation ratio) shows linearity of methylation quantitation as well as precision in quantitation at each methylation ratio (Figure 6). As an example of this method in total, RNA and DNA were co-isolated from mouse cerebellum and retina (n=4/group). Rhodopsin expression, selectively expressed in retinal tissue, was measured using qPCR. Expression was only detected in retina (Figure 7A). Using BSAS, CpG methylation levels were quantified in the rhodopsin promoter region. Cumulative methylation levels across the promoter region were &#62;80% in the cerebellum compared to &#60;15% in the retina (p &#60; 0.001, parametric t-test) (Figure 7B). BSAS methylation quantitation benefits from site-specific methylation quantitation, and methylation levels can be compared on a CpG-specific basis across any given genomic region. CpG methylation levels across the rhodopsin promoter were significantly higher in cerebellar samples when compared to retina samples (p &#60;0.001, parametric t-test at each CpG site) (Figure 7C).
<img alt="Figure 1" src="/files/ftp_upload/52488/52488fig1highres.jpg" /> 
Figure 1: Schematic of BSAS method. In this method, genomic DNA is bisulfite converted to modify unmethylation cytosines to uracils. Subsequently, during PCR these uracils are changed to thymines. PCR amplification is directed against regions of interest and highly enriches for just these sequences. The resulting PCR amplicons are made into dual-indexed libraries through a simple tagmentation process. Subsequently the libraries are sequenced on a benchtop next generation sequencer and the sequencing reads are mapped to the in silico converted reference sequence and percent methylation of cytosine&#8217;s is determined in a base-specific manner.
<img alt="Figure 2" src="/files/ftp_upload/52488/52488fig2highres.jpg" /> 
Figure 2: In silico bisulfite conversion of region of interest. Any region of interest from any genomic reference can be selected for in silico bisulfite conversion. Non-CpG cytosines are replaced with thymine&#8217;s in the reference sequence. CpG cytosines remain cytosines in the bisulfite converted reference.
<img alt="Figure 3" src="/files/ftp_upload/52488/52488fig3highres.jpg" /> 
Figure 3: Primer placement. Bisulfite PCR primers are designed against an in silico bisulfite converted reference sequence. Primers are to be designed so they do not overlap CG dinucleotides or designed adjacent to CG dinucleotides.
<img alt="Figure 4" src="/files/ftp_upload/52488/52488fig4highres.jpg" /> 
Figure 4: Examples of high and poor quality PCR amplicons and sequencing libraries. Representative electropherogram traces and gel show (A) ideal amplicon for transposome-mediated library generation with a single high concentration product. (B) Library generated from this amplicon shows high concentration and even size distribution. (C) Poor quality PCR amplicons can be low in size, concentration, and contain multiple products and/or primer dimers. These poor quality amplicons will lead to (D) failed library generation with low concentration and low size distribution.
<img alt="Figure 5" src="/files/ftp_upload/52488/52488fig5highres.jpg" /> 
Figure 5: Examples of sequencing outputs from methylation controls. (A) In silico converted reference sequencing with CpG sites highlighted in red. 0% methylation control mapped reads to the selected region showing thymines (greed highlight) mapping at CpG sites. (B) 100% methylation control mapped reads show cytosine&#8217;s (blue highlight) mapping at CpG sites.
<img alt="Figure 6" src="/files/ftp_upload/52488/52488fig6highres.jpg" /> 
Figure 6: Methylation quantitation across a range of methylation controls. Whole genome methylation controls mixed at methylation ratios (n = 3/ratio) from 0% to 100% were quantified using BSAS. Plotting quantitation of expected versus quantified demonstrates linearity of methylation control quantitation. Across all controls, there was high precision in methylation quantitation.
<img alt="Figure 7" src="/files/ftp_upload/52488/52488fig7highres.jpg" /> 
Figure 7: Example of paired DNA methylation and gene expression analysis from tissue samples. (A) Rhodopsin relative mRNA expression from mouse cerebellar and retinal tissue. (B) Rhodopsin average promoter methylation quantified by BSAS between cerebellar and retinal DNA (***p &#60;0.001, parametric t-test). C) CpG site-specific methylation levels across rhodopsin promoter (-244 to +6, relative to transcription start site [TSS]) in mouse cerebellum and retina DNA.

Watch this Scientific Journal Video about Targeted DNA Methylation Analysis by Next-generation Sequencing at JoVE.com

Targeted DNA Methylation Analysis by Next-generation 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 ...

targeted-dna-methylation-analysis-by-next-generation-sequencing

Nanyang Technological University

Research

JoVE Journal

Biology

36.9K Views.  University of Oklahoma College of Medicine. 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.

Video: Targeted DNA Methylation Analysis by Next-generation Sequencing

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 ...

Fluorescence-microscopy Screening and Next-generation Sequencing: Useful Tools for the Identification of Genes Involved in Organelle Integrity

This protocol describes a fluorescence microscope-based screening of Arabidopsis seedlings and describes how to map recessive mutations that alter the subcellular distribution of a specific tagged fluorescent marker in the secretory pathway. Arabidopsis is a powerful biological model for genetic studies because of its genome size, generation time, and conservation of molecular mechanisms among kingdoms. The array genotyping as an approach to map the mutation in alternative to the traditional method based on molecular markers is advantageous because it is relatively faster and may allow the mapping of several mutants in a really short time frame. This method allows the identification of proteins that can influence the integrity of any organelle in plants. Here, as an example, we propose a screen to map genes important for the integrity of the endoplasmic reticulum (ER). Our approach, however, can be easily extended to other plant cell organelles (for example see1,2), and thus represents an important step toward understanding the molecular basis governing other subcellular structures.

A fundamental quest in cell biology is to define the mechanisms that underlie the identity of the organelles that make eukaryotic cells. Here we propose a method to identify the genes responsible for the morphological and functional integrity of plant organelles using fluorescence microscopy and next-generation sequencing tools.

This protocol describes a fluorescence microscope-based screening of Arabidopsis seedlings and describes how to map recessive mutations that alter the ...

Next-generation Sequencing of 16S Ribosomal RNA Gene Amplicons

One of the major questions in microbial ecology is &#8220;who is there?&#8221; This question can be answered using various tools, but one of the long-lasting gold standards is to sequence 16S ribosomal RNA (rRNA) gene amplicons generated by domain-level PCR reactions amplifying from genomic DNA. Traditionally, this was performed by cloning and Sanger (capillary electrophoresis) sequencing of PCR amplicons. The advent of next-generation sequencing has tremendously simplified and increased the sequencing depth for 16S rRNA gene sequencing. The introduction of benchtop sequencers now allows small labs to perform their 16S rRNA sequencing in-house in a matter of days. Here, an approach for 16S rRNA gene amplicon sequencing using a benchtop next-generation sequencer is detailed. The environmental DNA is first amplified by PCR using primers that contain sequencing adapters and barcodes. They are then coupled to spherical particles via emulsion PCR. The particles are loaded on a disposable chip and the chip is inserted in the sequencing machine after which the sequencing is performed. The sequences are retrieved in fastq format, filtered and the barcodes are used to establish the sample membership of the reads. The filtered and binned reads are then further analyzed using publically available tools. An example analysis where the reads were classified with a taxonomy-finding algorithm within the software package Mothur is given. The method outlined here is simple, inexpensive and straightforward and should help smaller labs to take advantage from the ongoing genomic revolution.

Characterizing microbial community has been a longstanding goal in environmental microbiology. Next-generation sequencing methods now allow for the characterization of microbial communities at an unprecedented depth with minimal cost and labor. We detail here our approach to sequence bacterial 16S ribosomal RNA genes using a benchtop sequencer.

One of the major questions in microbial ecology is &#8220;who is there?&#8221; This question can be answered using various tools, but one of the ...

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.

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 ...

Fecal Glucocorticoid Analysis: Non-invasive Adrenal Monitoring in Equids

Adrenal activity can be assessed in the equine species by analysis of feces for corticosterone metabolites. During a potentially aversive situation, corticotrophin releasing hormone (CRH) is released from the hypothalamus in the brain. This stimulates the release of adrenocorticotrophic hormone (ACTH) from the pituitary gland, which in turn stimulates release of glucocorticoids from the adrenal gland. In horses the glucocorticoid corticosterone is responsible for several adaptations needed to support equine flight behaviour and subsequent removal from the aversive situation. Corticosterone metabolites can be detected in the feces of horses and assessment offers a non-invasive option to evaluate long term patterns of adrenal activity. Fecal assessment offers advantages over other techniques that monitor adrenal activity including blood plasma and saliva analysis. The non-invasive nature of the method avoids sampling stress which can confound results. It also allows the opportunity for repeated sampling over time and is ideal for studies in free ranging horses. This protocol describes the enzyme linked immunoassay (EIA) used to assess feces for corticosterone, in addition to the associated biochemical validation.

Adrenal activity can be assessed in the equine species by analysis of feces for corticosterone metabolites. The method offers a non-invasive option to assess long term patterns in both domestic and free ranging horses. This protocol describes the enzyme linked immunoassay involved and the associated biochemical validation.

Adrenal activity can be assessed in the equine species by analysis of feces for corticosterone metabolites. During a potentially aversive situation, ...

DamID-seq: Genome-wide Mapping of Protein-DNA Interactions by High Throughput Sequencing of Adenine-methylated DNA Fragments

The DNA adenine methyltransferase identification (DamID) assay is a powerful method to detect protein-DNA interactions both locally and genome-wide. It is an alternative approach to chromatin immunoprecipitation (ChIP). An expressed fusion protein consisting of the protein of interest and the E. coli DNA adenine methyltransferase can methylate the adenine base in GATC motifs near the sites of protein-DNA interactions. Adenine-methylated DNA fragments can then be specifically amplified and detected. The original DamID assay detects the genomic locations of methylated DNA fragments by hybridization to DNA microarrays, which is limited by the availability of microarrays and the density of predetermined probes. In this paper, we report the detailed protocol of integrating high throughput DNA sequencing into DamID (DamID-seq). The large number of short reads generated from DamID-seq enables detecting and localizing protein-DNA interactions genome-wide with high precision and sensitivity. We have used the DamID-seq assay to study genome-nuclear lamina (NL) interactions in mammalian cells, and have noticed that DamID-seq provides a high resolution and a wide dynamic range in detecting genome-NL interactions. The DamID-seq approach enables probing NL associations within gene structures and allows comparing genome-NL interaction maps with other functional genomic data, such as ChIP-seq and RNA-seq.

We describe herein an assay by coupling DNA adenine methyltransferase identification (DamID) to high throughput sequencing (DamID-seq). This improved method provides a higher resolution and a wider dynamic range, and allows analyzing DamID-seq data in conjunction with other high throughput sequencing data such as ChIP-seq, RNA-seq, etc.

The DNA adenine methyltransferase identification (DamID) assay is a powerful method to detect protein-DNA interactions both locally and genome-wide. It is ...

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 ...

Targeted RNA Sequencing Assay to Characterize Gene Expression and Genomic Alterations

RNA sequencing (RNAseq) is a versatile method that can be utilized to detect and characterize gene expression, mutations, gene fusions, and noncoding RNAs. Standard RNAseq requires 30 - 100 million sequencing reads and can include multiple RNA products such as mRNA and noncoding RNAs. We demonstrate how targeted RNAseq (capture) permits a focused study on selected RNA products using a desktop sequencer. RNAseq capture can characterize unannotated, low, or transiently expressed transcripts that may otherwise be missed using traditional RNAseq methods. Here we describe the extraction of RNA from cell lines, ribosomal RNA depletion, cDNA synthesis, preparation of barcoded libraries, hybridization and capture of targeted transcripts and multiplex sequencing on a desktop sequencer. We also outline the computational analysis pipeline, which includes quality control assessment, alignment, fusion detection, gene expression quantification and identification of single nucleotide variants. This assay allows for targeted transcript sequencing to characterize gene expression, gene fusions, and mutations.

We describe a targeted RNA sequencing-based method that includes preparation of indexed cDNA libraries, hybridization and capture with custom probes and data analysis to interrogate selected transcripts for gene expression, mutations, and gene fusions. Targeted RNAseq permits cost-effective, rapid evaluation of selected transcripts on a desktop sequencer.

RNA sequencing (RNAseq) is a versatile method that can be utilized to detect and characterize gene expression, mutations, gene fusions, and noncoding ...

Interactome-Seq: A Protocol for Domainome Library Construction, Validation and Selection by Phage Display and Next Generation Sequencing

Folding reporters are proteins with easily identifiable phenotypes, such as antibiotic resistance, whose folding and function is compromised when fused to poorly folding proteins or random open reading frames. We have developed a strategy where, by using TEM-1 &#946;-lactamase (the enzyme conferring ampicillin resistance) on a genomic scale, we can select collections of correctly folded protein domains from the coding portion of the DNA of any intronless genome. The protein fragments obtained by this approach, the so called &#34;domainome&#34;, will be well expressed and soluble, making them suitable for structural/functional studies.
By cloning and displaying the &#34;domainome&#34; directly in a phage display system, we have showed that it is possible to select specific protein domains with the desired binding properties (e.g., to other proteins or to antibodies), thus providing essential experimental information for gene annotation or antigen identification.
The identification of the most enriched clones in a selected polyclonal population can be achieved by using novel next-generation sequencing technologies (NGS). For these reasons, we introduce deep sequencing analysis of the library itself and the selection outputs to provide complete information on diversity, abundance and precise mapping of each of the selected fragment. The protocols presented here show the key steps for library construction, characterization, and validation.

The protocols described allow the construction, characterization and selection (against the target of choice) of a &#34;domainome&#34; library made from any DNA source. This is achieved by a research pipeline that combines different technologies: phage display, a folding reporter and next generation sequencing with a web tool for data analysis.

Folding reporters are proteins with easily identifiable phenotypes, such as antibiotic resistance, whose folding and function is compromised when fused to ...