We gratefully acknowledge the UND Genomics Core facility for outstanding technical assistance.
This work was funded by the National Institutes of Health [P20GM104360 to M.T., P20 GM104360 to A.D.] and start-up funds provided by the University of North Dakota School of Medicine and Health Sciences, Department of Biomedical Sciences [to M.T.].

1. Preparations before beginning the experiment
<ol>
	<li>Prepare lysis buffer stock. 
	NOTE: The optimal nuclear isolation buffer can be different for each cell type. We recommend testing both the hypotonic buffer used in the original paper8 and a CSK buffer12,13 for each cell type using trypan blue staining.
	<ol>
		<li>To prepare hypotonic buffer (Greenleaf buffer), mix 10 mM Tris-HCl pH 7.4, 10 mM NaCl, 3 mM MgCl2, and 0.1% NP-40.</li>
		<li>To prepare CSK buffer, mix 10 mM PIPES pH 6.8, 100 mM NaCl, 300 mM sucrose, 3 mM MgCl2, and 0.1% Triton X-100.</li>
	</ol>
	</li>
	<li>Ensure that the cell culture conditions are optimal; examine the cells under a light microscope to make sure that there are no elevated levels of apoptosis.</li>
	<li>Turn on the centrifuge to let it reach 4 &#176;C.</li>
</ol>
2. Cell harvest
<ol>
	<li>Aspirate media from a culture of cells at &#60;80% confluency. Cells are typically grown in a 35 mm or 60 mm dish, based on the number of cells needed, treatments, etc.</li>
	<li>Wash the cells with 1 mL or 3 mL of PBS.</li>
	<li>Add 0.5 mL or 1 mL of Trypsin and incubate for 2-3 min (depending on the cell types). Carefully use a 1 mL pipette to mix well.</li>
	<li>Add 1 mL or 2 mL of media to the plate and mix well. Transfer to a 15 mL conical tube. Centrifuge the cells in a 15 mL tube for 3 min at 300 x g. 
	NOTE: A swinging bucket rotor is recommended to minimize cell loss.</li>
	<li>Aspirate the media and resuspend the cells&#160;in 1 mL or 2 mL of medium. 
	NOTE: It is critical to prepare a homogenous single-cell solution. For tissues, a Dounce homogenizer can be used to homogenize tissues and minimize large cell clumps or aggregates. Some tissues require filtration (e.g., cell strainers) and/or FACS cell sorting to isolate viable cells and obtain a single cell solution.</li>
	<li>Count the cells using a cell counter. 
	NOTE: In this study, an automated cell counter (Table of Materials) was used.</li>
	<li>Centrifuge the required cell volume (e.g., a volume containing 1 x 106 cells based on cell count) at 300 x g for 3 min at room temperature (RT).</li>
	<li>Resuspend the cell pellet containing 1 x 106 cells in 1 mL of cold PBS. 
	NOTE: If cell numbers are less than 1 x 106&#160;cells, decrease the cold PBS volume to prepare the cell solution at the same concentration (1 x 106&#160;cells/mL).</li>
</ol>
3. Cell lysis
<ol>
	<li>Transfer 25 &#181;L (= 2.5 x 104 cells) to a new 1.7 mL microcentrifuge tube. Centrifuge&#160;for 5 min at 500 x g at 4 &#176;C and carefully discard the supernatant 
	NOTE: A swinging bucket rotor is recommended to minimize cell loss. Leave about 3 &#181;L of supernatant to ensure that the cells are not being discarded.</li>
	<li>Add 25 &#181;L of lysis buffer and resuspend the cells by gentle pipetting up and down. Incubate on ice for 5 min</li>
	<li>Centrifuge for 5 min at 500 x g at 4 &#176;C and carefully&#160;discard the supernatant. 
	NOTE: A swinging bucket rotor is recommended to minimize cell loss. Leave about 3 &#181;L of supernatant to ensure that the cells are not being discarded.</li>
</ol>
4. Tn5 tagmentation
<ol>
	<li>Immediately resuspend the pellet in 25 &#181;L of Tn5 reaction mixture. The Tn5 reaction mixture is as follows: 12.5 &#181;L of 2x Tagmentation DNA buffer (TD buffer), 1.25-5 &#181;L of Tn5 transposase, and 11.25-7.5 &#181;L of nuclease-free water. 
	NOTE: The standard concentration of Tn5 is 2.5 &#181;L for 2.5 x 104 cells.</li>
	<li>Incubate for 30 min at 37 &#176;C. 
	&#8203;NOTE: Mix the solution every 10 min.</li>
</ol>
5. DNA purification
NOTE: Purification is required before amplification. DNA purification is done using the MinElute PCR purification kit (Table of Materials).
<ol>
	<li>Add 5 &#181;L of 3 M sodium acetate (pH 5.2).</li>
	<li>Immediately add 125 &#181;L of Buffer PB and mix well.</li>
	<li>Follow the PCR purification kit protocol starting at step 2.
	<ol>
		<li>Place the spin column in a 2 mL collection tube.</li>
		<li>Apply the sample to the spin column and centrifuge for 1 min at 17,900 x g at RT.</li>
		<li>Discard the flow-through and place the spin column back in the same collection tube.</li>
		<li>Add 750 &#181;L of Buffer PE to the spin column and centrifuge for 1 min at 17,900 x g at RT.</li>
		<li>Discard the flow-through and place the spin column back in the same collection tube.</li>
		<li>Centrifuge the spin column for 5 min to dry the membrane completely.</li>
		<li>Discard the flow-through and place the spin column in a new 1.7 mL microcentrifuge tube.</li>
		<li>Elute the DNA fragments with 10 &#181;L of Buffer EB.</li>
		<li>Let the column stand for 1 min at RT.</li>
		<li>Centrifuge for 1 min at 17,900 x g at RT to elute the DNA.</li>
	</ol>
	</li>
</ol>
6. PCR amplification
NOTE: PCR amplification of transposed (tagmented) DNAs is necessary for sequencing. Nextera kit adapters (Table of Materials) were used in this example. The primers used in this study are listed in Table 1.
<ol>
	<li>Set up the following PCR reaction in 200 &#181;L PCR tubes (50 &#181;L for each sample). The PCR reaction mixture is as follows: 10 &#181;L of eluted DNA (step 5.), 2.5 &#181;L of 25 mM Adapter 1, 2.5 &#181;L of 25 mM Adapter 2, 25 &#181;L of 2x PCR master mix, and 10 &#181;L of nuclease-free water</li>
	<li>Perform PCR amplification program part 1 (Table 2). 
	NOTE: For initial amplification, the same conditions are used for all samples using a standard thermal cycler.</li>
	<li>Real-time qPCR (total of 14.5 &#181;L per sample)
	<ol>
		<li>Using 5 &#181;L of the product from the PCR amplification part 1 (step 6.2.), set up the following reaction, this time including SYBR gold and detection in a real-time PCR machine.</li>
		<li>Prepare the PCR reaction mixture as follows: 5 &#181;L of PCR product from step 6.2., 0.75 &#181;L of SYBR gold (1000x diluted), 5 &#181;L of 2x PCR master mix, and 3.75 &#181;L of nuclease-free water. 
		&#8203;NOTE: The precise cycle numbers for library amplification before reaching saturation for each sample are determined by qPCR (Table 3) to reduce the GC and size bias in PCR10. SYBR Gold was diluted with nuclease-free water. To make the diluted solution, 1 &#181;L of SYBR Gold (stock concentration 10,000x) was added to&#160;999 &#181;L of nuclease-free water. The run time of RT-qPCR mentioned in Table 3 is about 60 min.</li>
	</ol>
	</li>
	<li>Determine the required number of additional PCR cycles using the run parameters from step 6.3.
	<ol>
		<li>Number of cycles = 1/4 maximum (saturated) fluorescence intensity (typically 3 or 4 PCR cycles)</li>
	</ol>
	</li>
	<li>PCR amplification program part 2 (volume = 45 &#181;L; Table 4): Once you determine the cycle number, set up PCRs with additional cycles calculated at step 6.4.1. For instance, if the cycle number calculated is 3, set up 3 cycles in the program below. The total cycles would be 8 (5 in step 6.2., plus 3 additional cycles in step 6.5.)</li>
</ol>
7. Beads purification
NOTE: Here, AMPure XP (Table of Materials) beads were used.
<ol>
	<li>To the PCR products amplified in the previous step, add 150 &#181;L of beads at RT. 
	NOTE: Homemade AMPure beads14 can be used in this process.</li>
	<li>Incubate for 15 min at RT.</li>
	<li>Place the tubes on a magnetic stand for 5 min.</li>
	<li>Carefully remove the supernatant.</li>
	<li>Wash the beads with 200 &#181;L of 80% ethanol. 
	NOTE: 80% ethanol should be prepared fresh.</li>
	<li>Repeat step 7.4.</li>
	<li>Remove the ethanol completely.</li>
	<li>Air dry the samples for 10 min at RT.</li>
	<li>Resuspend with 50 &#181;L of elution buffer (10 mM Tris-HCL pH 8.0).</li>
	<li>Repeat the beads purification (steps 7.1.-7.9.) to further minimize primer contamination.</li>
</ol>
8. DNA concentration and quality check
<ol>
	<li>Measure the DNA concentration by using a nucleic acid quantification kit (Table of Materials).</li>
	<li>Check the amplified DNA fragments by gel electrophoresis using SYBR Gold (0.5x TBE, 1.5% agarose gel) or an automated electrophoresis system (e.g., TapeStation, bioanalyzer).</li>
</ol>
9. Sequencing
<ol>
	<li>Perform high-throughput sequencing using the amplified DNA samples12,13. 
	&#8203;NOTE: Amplified DNA samples are ready for high-throughput sequencing. To retain nucleosomal DNA fragment information, paired-end sequencing is recommended over single-end sequencing. Both ends of the DNA fragments indicate where the transposase bind. For mapping open chromatin regions, 30 million reads/sample is typically sufficient.</li>
</ol>
10. Data analysis
<ol>
	<li>Data filtration based on base call quality: Set the minimum average base quality score to 20 (99% accuracy).</li>
	<li>Adapter trimming: Remove the adapter sequences using an appropriate tool. 
	NOTE: The Trim Galore wrapper tool (http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/) was used to remove adapter sequences. For ATAC-seq libraries prepared by Illumina Tn5, the following sequence is used as an adapter sequence: CTGTCTCTTATACACATCT. The minimum sequence length to recognize adapter contamination is set to 5.</li>
	<li>Genome mapping: Map the reads to the appropriate genome with Bowtie using the following parameters &#34;-I 0 -X 2000 -m 1&#34;15. An example of mapping quality from MDA-MB-231 basal breast cancer cells is below: 
	Reads with at least one reported alignment: 52.79% 
	Reads that failed to align: 13.10% 
	Reads with alignments suppressed due to -m: 34.11%</li>
	<li>Deduplication: Mark the PCR duplicates by Picard tools16.</li>
	<li>Generating a genome coverage file for data visualization: Convert the deduplicated paired-end reads to single-end read fragments. Genome coverage tracks are obtained using genomeCoverageBed from bedtools17.</li>
</ol>

<ol>
	<li>Klemm, S. L., Shipony, Z., Greenleaf, W. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chromatin+accessibility+and+the+regulatory+epigenome.">Chromatin accessibility and the regulatory epigenome.</a> Nature Reviews Genetics. 20 (4), 207-220 (2019).</li><li>Liu, Y. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Clinical+implications+of+chromatin+accessibility+in+human+cancers.">Clinical implications of chromatin accessibility in human cancers.</a> Oncotarget. 11 (18), 1666-1678 (2020).</li><li>Sanghi, 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=Chromatin+accessibility+associates+with+protein-RNA+correlation+in+human+cancer.">Chromatin accessibility associates with protein-RNA correlation in human cancer.</a> Nature Communications. 12 (1), 5732 (2021).</li><li>Corces, M. 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=The+chromatin+accessibility+landscape+of+primary+human+cancers.">The chromatin accessibility landscape of primary human cancers.</a> Science. 362 (6413), (2018).</li><li>Boyle, A. P., 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=High-resolution+mapping+and+characterization+of+open+chromatin+across+the+genome.">High-resolution mapping and characterization of open chromatin across the genome.</a> Cell. 132 (2), 311-322 (2008).</li><li>Giresi, P. G., Kim, J., McDaniell, R. M., Iyer, V. R., Lieb, J. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=FAIRE+(Formaldehyde-Assisted+Isolation+of+Regulatory+Elements)+isolates+active+regulatory+elements+from+human+chromatin.">FAIRE (Formaldehyde-Assisted Isolation of Regulatory Elements) isolates active regulatory elements from human chromatin.</a> Genome Research. 17 (6), 877-885 (2007).</li><li>Pardo, C. E., Nabilsi, N. H., Darst, R. P., Kladde, M. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Integrated+DNA+methylation+and+chromatin+structural+analysis+at+single-molecule+resolution.">Integrated DNA methylation and chromatin structural analysis at single-molecule resolution.</a> Methods in Molecular Biology. 1288, 123-141 (2015).</li><li>Buenrostro, J. D., Giresi, P. G., Zaba, L. C., Chang, H. Y., Greenleaf, W. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transposition+of+native+chromatin+for+fast+and+sensitive+epigenomic+profiling+of+open+chromatin,+DNA-binding+proteins+and+nucleosome+position.">Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position.</a> Nature Methods. 10 (12), 1213-1218 (2013).</li><li>Song, L., 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=Open+chromatin+defined+by+DNaseI+and+FAIRE+identifies+regulatory+elements+that+shape+cell-type+identity.">Open chromatin defined by DNaseI and FAIRE identifies regulatory elements that shape cell-type identity.</a> Genome Research. 21 (10), 1757-1767 (2011).</li><li>Buenrostro, J. D., Wu, B., Chang, H. Y., Greenleaf, W. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=ATAC-seq:+A+method+for+assaying+chromatin+accessibility+genome-wide.">ATAC-seq: A method for assaying chromatin accessibility genome-wide.</a> Current Protocols in Molecular Biology. , 1 (2015).</li><li>Mehrmohamadi, M., Sepehri, M. H., Nazer, N., Norouzi, M. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+comparative+overview+of+epigenomic+profiling+methods.">A comparative overview of epigenomic profiling methods.</a> Frontiers in Cell and Developmental Biology. 9, 714687 (2021).</li><li>Fujiwara, S., Baek, S., Varticovski, L., Kim, S., Hager, G. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=High+quality+ATAC-Seq+data+recovered+from+cryopreserved+breast+cell+lines+and+tissue.">High quality ATAC-Seq data recovered from cryopreserved breast cell lines and tissue.</a> Scientific Reports. 9 (1), 516 (2019).</li><li>Takaku, 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=GATA3-dependent+cellular+reprogramming+requires+activation-domain+dependent+recruitment+of+a+chromatin+remodeler.">GATA3-dependent cellular reprogramming requires activation-domain dependent recruitment of a chromatin remodeler.</a> Genome Biology. 17, 36 (2016).</li><li>Rohland, N., Reich, D. C. o. s. t. -. e. f. f. e. c. t. i. v. e. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=high-throughput+DNA+sequencing+libraries+for+multiplexed+target+capture.">high-throughput DNA sequencing libraries for multiplexed target capture.</a> Genome Research. 22 (5), 939-946 (2012).</li><li>Langmead, B., Trapnell, C., Pop, M., Salzberg, S. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ultrafast+and+memory-efficient+alignment+of+short+DNA+sequences+to+the+human+genome.">Ultrafast and memory-efficient alignment of short DNA sequences to the human genome.</a> Genome Biology. 10 (3), 25 (2009).</li><li>. Picard Toolkit Broad Institute, GitHub Repository Available from: <a target="_blank" href="https://broadinstitute.github.io/picard/">https://broadinstitute.github.io/picard/</a> (2022)</li><li>Quinlan, A. R., Hall, I. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=BEDTools:+A+flexible+suite+of+utilities+for+comparing+genomic+features.">BEDTools: A flexible suite of utilities for comparing genomic features.</a> Bioinformatics. 26 (6), 841-842 (2010).</li><li>Strober, W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Trypan+blue+exclusion+test+of+cell+viability.">Trypan blue exclusion test of cell viability.</a> Current Protocols in Immunology. 111, 1-3 (2015).</li><li>Zhou, B., 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=INO80+governs+superenhancer-mediated+oncogenic+transcription+and+tumor+growth+in+melanoma.">INO80 governs superenhancer-mediated oncogenic transcription and tumor growth in melanoma.</a> Genes &amp; Development. 30 (12), 1440-1453 (2016).</li><li>Takaku, 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=GATA3+zinc+finger+2+mutations+reprogram+the+breast+cancer+transcriptional+network.">GATA3 zinc finger 2 mutations reprogram the breast cancer transcriptional network.</a> Nature Communications. 9 (1), 1059 (2018).</li><li>Oldfield, A. 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=NF-Y+controls+fidelity+of+transcription+initiation+at+gene+promoters+through+maintenance+of+the+nucleosome-depleted+region.">NF-Y controls fidelity of transcription initiation at gene promoters through maintenance of the nucleosome-depleted region.</a> Nature Communications. 10 (1), 3072 (2019).</li><li>Langer, L. F., Ward, J. M., Archer, T. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Tumor+suppressor+SMARCB1+suppresses+super-enhancers+to+govern+hESC+lineage+determination.">Tumor suppressor SMARCB1 suppresses super-enhancers to govern hESC lineage determination.</a> Elife. 8, 45672 (2019).</li><li>Elrod, N. 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=The+integrator+complex+attenuates+promoter-proximal+transcription+at+protein-coding+genes.">The integrator complex attenuates promoter-proximal transcription at protein-coding genes.</a> Molecular Cell. 76 (5), 738-752 (2019).</li><li>Corces, M. 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=An+improved+ATAC-seq+protocol+reduces+background+and+enables+interrogation+of+frozen+tissues.">An improved ATAC-seq protocol reduces background and enables interrogation of frozen tissues.</a> Nature Methods. 14 (10), 959-962 (2017).</li><li>Porter, J. 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=Global+inhibition+with+specific+activation:+How+p53+and+MYC+redistribute+the+transcriptome+in+the+DNA+double-strand+break+response.">Global inhibition with specific activation: How p53 and MYC redistribute the transcriptome in the DNA double-strand break response.</a> Molecular Cell. 67 (6), 1013-1025 (2017).</li><li>Dhasarathy, A., Kajita, M., Wade, P. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+transcription+factor+snail+mediates+epithelial+to+mesenchymal+transitions+by+repression+of+estrogen+receptor-alpha.">The transcription factor snail mediates epithelial to mesenchymal transitions by repression of estrogen receptor-alpha.</a> Molecular Endocrinology. 21 (12), 2907-2918 (2007).</li><li>McCarthy, M. T., O'Callaghan, C. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=PeaKDEck:+A+kernel+density+estimator-based+peak+calling+program+for+DNaseI-seq+data.">PeaKDEck: A kernel density estimator-based peak calling program for DNaseI-seq data.</a> Bioinformatics. 30 (9), 1302-1304 (2014).</li><li>Grandi, F. C., Modi, H., Kampman, L., Corces, M. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Chromatin+accessibility+profiling+by+ATAC-seq.">Chromatin accessibility profiling by ATAC-seq.</a> Nature Protocols. 7 (6), 1518-1552 (2022).</li><li>Heinz, S., 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=Simple+combinations+of+lineage-determining+transcription+factors+prime+cis-regulatory+elements+required+for+macrophage+and+B+cell+identities.">Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities.</a> Molecular Cell. 38 (4), 576-589 (2010).</li><li>Hanahan, D., Weinberg, 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=Hallmarks+of+cancer:+The+next+generation.">Hallmarks of cancer: The next generation.</a> Cell. 144 (5), 646-674 (2011).</li><li>Garcia-Martinez, L., Zhang, Y., Nakata, Y., Chan, H. L., Morey, L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Epigenetic+mechanisms+in+breast+cancer+therapy+and+resistance.">Epigenetic mechanisms in breast cancer therapy and resistance.</a> Nature Communications. 12 (1), 1786 (2021).</li><li>Bhattacharya, 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=The+calcium+channel+proteins+ORAI3+and+STIM1+mediate+TGF-&#946;+induced+Snai1+expression.">The calcium channel proteins ORAI3 and STIM1 mediate TGF-&#946; induced Snai1 expression.</a> Oncotarget. 9 (50), 29468-29483 (2018).</li><li>Schep, A. 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=Structured+nucleosome+fingerprints+enable+high-resolution+mapping+of+chromatin+architecture+within+regulatory+regions.">Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions.</a> Genome Research. 25 (11), 1757-1770 (2015).</li></ol>

The authors declare that there are no relevant or material financial interests that relate to the research described in this paper.

ATAC-seq has been widely used for mapping open and active chromatin regions. Cancer cell progression is frequently driven by genetic alterations and epigenetic reprogramming, resulting in altered chromatin accessibility and gene expression. An example of this reprogramming is seen during the epithelial-to-mesenchymal transition (EMT) and its reverse process, mesenchymal-to-epithelial transition (MET), which are known to be key cellular reprogramming processes during tumor metastasis30. Another example is the acquisition of drug resistance against hormone therapies is often observed in luminal breast tumors31. ATAC-seq is useful to monitor such cell phenotypic alterations at the chromatin level and can be used to predict a cell of origin and cancer subtypes4.
To precisely measure chromatin remodeling by ATAC-seq, the establishment of a reproducible and consistent protocol is necessary. In this manuscript, a standard strategy for ATAC-seq library optimization is described. The ATAC-seq data from MDA-MB-231 and NMuMG cell lines are used as examples of optimization processes. NMuMG cells have been used to study EMT, and MDA-MB-231 cells have been used to study its reverse process, MET. These cellular models have been previously used to study epigenetic alterations during EMT and MET13,32. The use of an appropriate buffer for cell lysis and Tn5 concentration is the key component of successful ATAC-seq library preparation. In addition to these components, a high percentage (&#62;90%) of cell viability, appropriate concentrations of detergent in the lysis buffer, and the preparation of the single-cell suspension are also very important. When Tn5 digestion efficiency is significantly different across samples, it is challenging to correct data variations. This is due to the lack of internal and external controls to quantitatively assess digestion efficiency. Cellular characteristics or properties might also be different between the cells from the control group versus the cells from the test group. Therefore, careful consideration of the experimental conditions, including the choice of lysis buffer, is necessary to obtain high-quality ATAC-seq data from both experimental groups (Figure 5 and Figure 6).
Besides open chromatin profiling, ATAC-seq has been used to identify transcription factor footprints and nucleosome positioning8,33. It is important to note that such information is mostly derived from open chromatin regions (Figure 1F). Therefore, results would be biased toward open chromatin regions. Nevertheless, ATAC-seq is a powerful tool for studying chromatin architecture. This revolutionary tool has more recently been adopted for single-cell genomics and continues to contribute to improving our understanding of gene regulation and chromatin biology.
DATA AVAILABILITY: 
The data presented in this protocol are available at Gene Expression Omnibus under Accession Number GSE202791. The ATAC-seq data from GSE72141 and GSE99479 were also used in the analysis.

Chromatin accessibility is a key requirement for the regulation of gene expression on a genome-wide scale1. Changes in chromatin accessibility are frequently associated with several disease states, including cancer2,3,4. Over the years, numerous techniques have been developed to enable researchers to probe the chromatin landscape by mapping regions of chromatin accessibility. Some of them include DNase-seq (DNase I hypersensitive sites sequencing)5, FAIRE-seq (formaldehyde-assisted isolation of regulatory elements)6, MAPit (methyltransferase accessibility protocol for individual templates)7, and the focus of this paper, ATAC-seq (assay for transposase-accessible chromatin)8. DNase-seq maps accessible regions by employing a key feature of DNase, namely the preferential digestion of naked DNA free from histones and other proteins such as transcription factors5. FAIRE-seq, similar to ChIP-seq, utilizes formaldehyde crosslinking and sonication, except no immunoprecipitation is involved, and the nucleosome-free regions are isolated by phenol-chloroform extraction6. The MAPit method uses a GC methyltransferase to probe chromatin structure at single-molecule resolution7. ATAC-seq relies on the hyperactive transposase, Tn58. The Tn5 transposase preferentially binds to open chromatin regions and inserts sequencing adapters into accessible regions. Tn5 operates through a DNA-mediated &#34;cut and paste&#34; mechanism, whereby the transposase preloaded with adapters binds to open chromatin sites, cuts DNA, and ligates the adapters8. Tn5 bound regions are recovered by PCR amplification using primers that anneal to these adapters. FAIRE-seq and DNase-seq require a large amount of starting material (~100,000 cells to 225,000 cells) and a separate library preparation step before sequencing9. On the other hand, the ATAC-seq protocol is relatively simple and requires a small number of cells (&#60;50,000 cells)10. Unlike the FAIRE-seq and DNase-seq techniques, the sequencing library preparation of ATAC-seq is relatively easy, as the isolated DNA sample is already being tagged with the sequencing adapters by Tn5. Therefore, only the PCR amplification step with appropriate primers is needed to complete the library preparation, and the prior processing steps such as end-repair and adapter ligation need not be performed, thus saving time11. Secondly, ATAC-seq avoids the need for bisulfite conversion, cloning, and amplification with region-specific primers required for MAPit7. Due to these advantages, ATAC-seq has become a hugely popular method for defining open chromatin regions. Although the ATAC-seq method is simple, multiple steps require optimization to obtain high-quality and reproducible data. This manuscript discusses optimization procedures for standard ATAC-seq library preparation, especially highlighting three parameters: (1) lysis buffer composition, (2) Tn5 transposase concentration, and (3) cell number. In addition, this paper provides example data from the optimization conditions using both cancerous and non-cancerous adherent epithelial cells.

<table><tbody><tr><td>1.5 mL microcentrifuge tubes</td><td>USA Scientific</td><td>1615-5500</td><td>Natural</td></tr><tr><td>10 &amp;micro;L XL TipOne tips</td><td>USA Scientific</td><td>1120-3810</td><td>Filtered and low-retention</td></tr><tr><td>100 &amp;micro;L XL TipOne RPT tips</td><td>USA Scientific</td><td>1182-1830</td><td>Filtered and low-retention</td></tr><tr><td>100 &amp;micro;L XL TipOne tips</td><td>USA Scientific</td><td>1120-1840</td><td>Filtered and low-retention. Beveled Grade</td></tr><tr><td>15 mL Conical Centrifuge Tubes</td><td>Corning</td><td>352096</td><td></td></tr><tr><td>20 &amp;micro;L TipOne RPT tips</td><td>USA Scientific</td><td>1183-1810</td><td>Filtered and low-retention</td></tr><tr><td>200 &amp;micro;L TipOne RPT tips</td><td>USA Scientific</td><td>1180-8810</td><td>Filtered and low-retention</td></tr><tr><td>50 mL Centrifuge Tubes</td><td>Fisherbrand</td><td>06-443-19</td><td></td></tr><tr><td>Agarose</td><td>ThermoFisher Scientific</td><td>YBP136010</td><td>Genetic Analysis Grade</td></tr><tr><td>All the cell lines used in this study are obtained from ATCC</td><td>ATCC</td><td></td><td></td></tr><tr><td>Allegra X-30R Centrifuge</td><td>Beckman Coulter</td><td>364658</td><td>SX2415</td></tr><tr><td>AMPure XP beads</td><td>Beckman Coulter</td><td>A63881</td><td>Bead purification kit</td></tr><tr><td>CellDrop Cell Counter</td><td>DeNovix</td><td>CellDrop FL</td><td>Cell counter</td></tr><tr><td>EDTA</td><td>MilliporeSigma</td><td>EDS</td><td>BioUltra, anhydrous, &amp;ge;99% (titration)</td></tr><tr><td>EGTA</td><td>MilliporeSigma</td><td>E3889</td><td></td></tr><tr><td>Ethanol 100%</td><td>ThermoFisher Scientific</td><td>AC615100020</td><td>Anhydrous; Fisher Scientific - Decon Labs Sterilization Products</td></tr><tr><td>Fetal Bovine Serum - TET Tested</td><td>R&amp;amp;D Systems</td><td>S10350</td><td>Triple 0.1 &amp;micro;m filtered</td></tr><tr><td>Gibco DMEM 1x</td><td>ThermoFisher Scientific</td><td>11965092</td><td>[+] 4.5 g/L D-glucose; [+] L-Glutamine; [-] Sodium pyruvate</td></tr><tr><td>Gibco PBS 1x</td><td>ThermoFisher Scientific</td><td>10010023</td><td>pH 7.4</td></tr><tr><td>Gibco Trypsin-EDTA 1x</td><td>ThermoFisher Scientific</td><td>25200056</td><td>(0.25%), phenol red</td></tr><tr><td>Glycerol</td><td>IBI Scientific</td><td>56-81-5</td><td></td></tr><tr><td>Glycine</td><td>MilliporeSigma</td><td>G8898</td><td></td></tr><tr><td>HCl</td><td>MilliporeSigma</td><td>H1758</td><td></td></tr><tr><td>HEPES</td><td>MilliporeSigma</td><td>H3375</td><td></td></tr><tr><td>Invitrogen Qubit Fluorometer</td><td>ThermoFisher Scientific</td><td>Q32857</td><td></td></tr><tr><td>MgCl&lt;sub&gt;2&lt;/sub&gt;</td><td>MilliporeSigma</td><td>M3634</td><td></td></tr><tr><td>MinElute PCR Purification kit</td><td>Qiagen</td><td>28004</td><td>DNA purification kit</td></tr><tr><td>NaCl</td><td>IBI Scientific</td><td>7647-14-5</td><td></td></tr><tr><td>NaOH</td><td>MilliporeSigma</td><td>S8045</td><td>BioXtra, &amp;ge;98% (acidimetric), pellets (anhydrous)</td></tr><tr><td>NEBNext High-Fidelity 2x PCR Master Mix</td><td>New England Biolabs</td><td>M0541</td><td></td></tr><tr><td>Nextera DNA Sample Preparation Kit</td><td>Illumina</td><td>FC-121-1030</td><td>2x TD and Tn5 Transposase</td></tr><tr><td>NP - 40 (IGEPAL CA-630)</td><td>MilliporeSigma</td><td>I8896</td><td>for molecular biology</td></tr><tr><td>PCR Detection System</td><td>BioRad</td><td>1855484</td><td>CFX384 Real-Time System. C1000 Touch Thermal Cycler</td></tr><tr><td>PIPES</td><td>MilliporeSigma</td><td>P1851</td><td>BioPerformance Certified, suitable for cell culture</td></tr><tr><td>Qubit dsDNA HS Assay kit</td><td>ThermoFisher Scientific</td><td>Q32854</td><td>Invitrogen; Nucleic acid quantitation kit</td></tr><tr><td>Quibit Assay Tubes</td><td>ThermoFisher Scientific</td><td>Q32856</td><td>Invitrogen</td></tr><tr><td>SDS</td><td>MilliporeSigma</td><td>L3771</td><td></td></tr><tr><td>Sodium Acetate</td><td>Homemade</td><td>-</td><td>pH 5.2</td></tr><tr><td>Sucrose</td><td>IBI Scientific</td><td>57-50-1</td><td></td></tr><tr><td>SYBR Gold</td><td>ThermoFisher Scientific</td><td>S11494</td><td></td></tr><tr><td>SYBR Green Supermix, 1.25 mL</td><td>BioRad</td><td>1708882</td><td></td></tr><tr><td>T100 Thermal Cycler</td><td>BioRad</td><td>1861096</td><td></td></tr><tr><td>TempAssure 0.2 mL PCR 8-Tube Strips</td><td>USA Scientific</td><td>1402-4700</td><td>Flex-free, natural, polypropylene</td></tr><tr><td>TempPlate 384-WELL PCR PLATE</td><td>USA Scientific</td><td>1438-4700</td><td>Single notch. Natural polypropylene</td></tr><tr><td>Tris Base</td><td>MilliporeSigma</td><td>648311</td><td>ULTROL Grade</td></tr><tr><td>Triton x-100</td><td>IBI Scientific</td><td>9002-93-1</td><td></td></tr><tr><td>TrueSeq Dual Index Sequencing Primer Kit</td><td>Illumina</td><td>PE-121-1003</td><td>paired-end</td></tr><tr><td>Trypan Blue Stain</td><td>ThermoFisher Scientific</td><td>Q32851</td><td></td></tr><tr><td>Tween-20</td><td>MilliporeSigma</td><td>P7949</td><td>BioXtra, viscous liquid</td></tr><tr><td>Water</td><td>MilliporeSigma</td><td>W3500</td><td>sterile-filtered, BioReagent, suitable for cell culture</td></tr></tbody></table>

atac-seq optimization for cancer epigenetics research

The assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) probes deoxyribonucleic acid (DNA) accessibility using the hyperactive Tn5 transposase. Tn5 cuts and ligates adapters for high-throughput sequencing within accessible chromatin regions. In eukaryotic cells, genomic DNA is packaged into chromatin, a complex of DNA, histones, and other proteins, which acts as a physical barrier to the transcriptional machinery. In response to extrinsic signals, transcription factors recruit chromatin remodeling complexes to enable access to the transcriptional machinery for gene activation. Therefore, identifying open chromatin regions is useful when monitoring enhancer and gene promoter activities during biological events such as cancer progression. Since this protocol is easy to use and has a low cell input requirement, ATAC-seq has been widely adopted to define open chromatin regions in various cell types, including cancer cells. For successful data acquisition, several parameters need to be considered when preparing ATAC-seq libraries. Among them, the choice of cell lysis buffer, the titration of the Tn5 enzyme, and the starting volume of cells are crucial for ATAC-seq library preparation in cancer cells. Optimization is essential for generating high-quality data. Here, we provide a detailed description of the ATAC-seq optimization methods for epithelial cell types.

To obtain successful and high-quality ATAC-seq data, it is important to optimize the experimental conditions. ATAC-seq library preparation can be separated into the five major steps (Figure 1), namely cell lysis, tagmentation (fragmentation and adapter insertion by Tn5), genomic DNA purification, PCR amplification, and data analysis. As an initial process, the cell lysis (nuclear isolation) buffer must be first optimized for each cell type. Either the hypotonic buffer described in the original ATAC-seq paper8 or the CSK buffer12,13 was used to lyse breast cancer cells. Trypan blue staining can be used to confirm nuclei isolation18. For ATAC-seq library preparation in human cancer cell lines such as MDA-MB-231, T47D, MCF7, and A375 cells, CSK buffer has been used by multiple groups12,13,19,20. CSK buffer was also used for ATAC-seq library preparation in other cell types such as embryonic stem cells21,22 and Drosophila S2 cells23. For the lysis buffer optimization, Trypan blue staining is useful to assess the cell lysis efficiency. When NMuMG, non-cancerous mouse mammary gland epithelial cells, were treated with the hypotonic buffer8 (see step 1) and CSK buffer, higher cell lysis efficiency was observed in CSK-treated cells (Figure 2). For frozen tissues, Omni-ATAC has been shown to improve ATAC-seq signals24. In the Omni-ATAC, a buffer containing 0.1% NP40, 0.1% Tween-20, and 0.01% Digitonin is used for cell lysis. Although the Digitonin-based buffer is known to decrease the fraction of mitochondrial DNA, the signal-to-noise ratios and TSS read counts from ATAC-seq with CSK buffer were shown to be better than those from Omni-ATAC in breast cancer cells12. Therefore, the rest of this manuscript is mainly focused on the results from ATAC-seq data with CSK buffer.
Following the determination of the best cell lysis buffer (nuclear isolation) conditions, it is important to optimize the ratios between input cell number and Tn5 transposase concentration12,25. Typically, the Tn5 concentration can be titrated by varying the volume of Tn5 from 1.25 µL, 2.5 µL, to 5 µL in a total volume of 25 µL reaction mixture. Alternatively, input cell numbers can be changed (e.g., from 10,000 to 100,000), while maintaining the Tn5 enzyme unchanged at 2.5 µL to optimize the Tn5 tagmentation reaction. To evaluate the quality of ATAC-seq libraries and the efficiency of Tn5-induced chromatin digestion, the library PCR products can be analyzed by agarose gel electrophoresis. Figure 3 shows examples of successful ATAC-seq libraries. Typically, PCR amplification of 8 or 9 cycles (total PCR cycles including initial PCR) results in 3-6 ng/µL ATAC-seq library concentration. As Tn5 preferentially "attacks" nucleosome-free regions at open chromatin (Figure 1), the nucleosome ladders should be obvious on the gel image (Figure 3). Ethidium bromide or SYBR Gold can be used to detect amplified DNAs on the agarose gel.
Next, qPCR analysis using ATAC-seq libraries can be used to briefly check the quality of the libraries and the fragment enrichment at open chromatin regions such as promoters of active genes (Figure 4). Primers can be designed to generate &lt;80 bp amplicons using promoters of known active genes as positive controls and known "closed" (inactive) chromatin regions as negative controls (Figure 4A). Among the tested conditions, higher enrichment was observed in the T47D ATAC-seq libraries with CSK buffer (Figure 4B). As expected, we found more enrichment using primers K and L at the Estrogen Receptor alpha (ESR1) gene promoter region relative to a closed region upstream of ESR1 with primer C (Figure 4B)26. Figure 5 shows examples of ATAC-seq data with different Tn5 concentrations (25,000 cells were used). In this case, higher background noise was detected in the lowest (1.25 µL of Tn5) Tn5 condition (Figure 5, vertical bars). ATAC-seq data from 2.5 µL or 5 µL of Tn5 showed similar levels of ATAC-seq signals at open chromatin regions with lower background signals. Considering the potential over-digestion by Tn5 at the higher concentration, it can be concluded that 2.5 µL of Tn5 is suitable for this cell type.
We also performed the ATAC-seq using different cell numbers in NMuMG cells, frequently used to study the epithelial-mesenchymal transition (Figure 6). To optimize the starting volume, four different cell numbers (25,000, 50,000, 75,000 and 100,000) were used for the ATAC-seq library optimization. The cells were lysed with CSK buffer, and nuclei were incubated with 2.5 µL of Tn5 transposase in 25 µL of the reaction mixture. To explore the efficacy of Tn5 digestion and nucleosomal ladders, the size distribution of inserted DNAs was bioinformatically analyzed (Figure 6A). When lower cell numbers were used, nucleosome-free DNA fragments (&lt;100 bp) were more enriched in the sequencing data. On the other hand, the fraction of mono-nucleosomal DNA fragments (175-225 bp) was less in the low cell input samples (25,000 cells and 50,000 cells) compared to higher cell input conditions (75,000 cells and 100,000 cells). Although differential DNA fragment patterns were observed between the samples, input cell number appears to have minimal impact on the enrichment of ATAC-seq signals at promoters and other open chromatin regions (Figure 6B). To further investigate the impact of input cell numbers on downstream analysis, ATAC-seq peaks were defined. The ATAC-seq peaks were determined by the PeaKDEck peak calling program27 using the following parameters: -sig 0.0001 -bin 300 -back 3,000 -npBack 2,500,000. From ~4.9 million uniquely mapped reads, 38,878, 43,832, 41,509, and 45,530 peaks were observed in ATAC-seq data from 25,000, 50,000, 75,000, and 100,000 cell inputs, respectively. The 100,000-cell input condition showed the highest number of peaks. Since the Transcription Start Site (TSS) enrichment score has been used to assess the quality of ATAC-seq data, the TSS score of each ATAC-seq condition was calculated using the GitHub program: Jiananlin/TSS_enrichment_score_calculation28. The TSS scores for the 25,000, 50,000, 75,000, and 100,000 cell input conditions were 9.03, 9,02, 8.82, and 8.28, respectively (the TSS enrichment score is considered ideal if it is greater than 728). This indicated a gradual decrease in the TSS enrichment score with increasing cell numbers. While the largest peak number was observed in the 100,000 cell input condition, the 25,000 cell input condition gave a better TSS score. Peak annotation analysis by Homer29 indicated that the majority of the increased peaks are categorized as promoter distal peaks (Figure 6C). Based on these results, it can be concluded that the higher cell number inputs can produce a higher number of peaks and more frequently detect non-promoter peaks compared to the lower cell number inputs in this condition.
<img alt="ATAC-seq workflow diagram: Chromatin accessibility analysis using Tn5 transposase and sequencing results." class="xfigimg" src="/files/ftp_upload/64242/64242fig01.jpg" data-altupdatedat="1747911167000" data-alt="Figure 1"> 
Figure 1: Outline of ATAC-seq method. (A) ATAC-seq measures chromatin accessibility using a hyperactive Tn5 transposase preloaded with forward and reverse adapters. The various steps in the process include (B) cell lysis (including optimization of lysis buffer composition, titration of Tn5, and the number of cells used); (C) transposition (binding of Tn5 to open/accessible chromatin); (D) tagmentation of chromatin; (E) DNA fragment purification and amplification of libraries, and (F) sequencing of libraries and data analysis. Fragment length distribution analysis of ATAC-seq libraries typically shows an initial peak around ~50 bp (nucleosome-free regions) and another peak at ~200 bp (mono-nucleosome). Therefore, without computational or experimental size filtration, ATAC-seq signals contain both nucleosome-free and nucleosomal regions in open chromatin. Created with BioRender.com. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig01large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Cell microscopy in hypotonic vs. CSK buffer; comparative analysis, 50 µm scale, microscopy image." class="xfigimg" src="/files/ftp_upload/64242/64242fig02.jpg" data-altupdatedat="1747911167000" data-alt="Figure 2"> 
Figure 2: Comparison of cell lysis efficiency using trypan blue. (A) NMuMG cells were treated with the hypotonic buffer8 and stained with trypan blue. (B) NMuMG cells were treated with CSK buffer and stained with trypan blue. The scale bars indicate 50 µm. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig02large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Gel electrophoresis results; DNA purification by magnetic beads, agarose gels comparison, automated electrophoresis." class="xfigimg" src="/files/ftp_upload/64242/64242fig03.jpg" data-altupdatedat="1747911167000" data-alt="Figure 3"> 
Figure 3: Amplified DNA fragment analysis of ATAC-seq libraries. ATAC-seq libraries were prepared from MDA-MB-231 breast cancer cells. (A) After PCR amplification, PCR products were analyzed on a 0.5x TBE, 1.5% agarose gel before and after beads purification. Primers are mostly removed by the initial beads purification. DNA band patterns from (B) 1x TAE, 1% agarose gel electrophoresis and (C) an automated electrophoresis are also indicated. SYBR Gold was used to visualize the DNA fragments shown in A and B. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig03large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="ATAC-seq data analysis at ESR1 locus; diagram of sequencing peaks and enrichment bar chart results." class="xfigimg" src="/files/ftp_upload/64242/64242fig04.jpg" data-altupdatedat="1747911167000" data-alt="Figure 4"> 
Figure 4: Enrichment analysis of ATAC-seq libraries. (A) Genome browser track showing the ESR1 locus. The genome coverage of the ATAC-seq data in T47D cells20 is shown. Primers from a previous publication were used26, and their target regions are highlighted in yellow. (B) Bar graph depicting primer amplicon enrichment in ATAC-seq libraries. ATAC-seq libraries were prepared by the hypotonic buffer or CSK buffer. Different concentrations of Tn5 were also used to generate ATAC-seq libraries. 1 ng of each library was used to perform qPCR. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig04large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="ATAC-seq data comparison with varying Tn5 concentrations; genome accessibility, histogram, ENCODE data." class="xfigimg" src="/files/ftp_upload/64242/64242fig05.jpg" data-altupdatedat="1747911167000" data-alt="Figure 5"> 
Figure 5: Genome browser visualization for ATAC-seq optimization. Different amounts of Tn5 transposase were added during library preparation. The 1.25 µL of Tn5 condition shows higher background signals (highlighted in yellow), while the other two conditions look similar. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig05large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="DNA fragment size distribution, peak analysis; graphs show read counts, genomic regions, peak numbers." class="xfigimg" src="/files/ftp_upload/64242/64242fig06.jpg" data-altupdatedat="1747911167000" data-alt="Figure 6"> 
Figure 6: DNA fragment size distribution of ATAC-seq libraries. (A) ATAC-seq libraries were prepared by using the indicated cell numbers. The ATAC-seq data from the 25,000 cell input condition showed the highest enrichment of nucleosome-free fragments, while the 100,000 cell input condition presented the highest mono-nucleosomal DNAs. (B) Genome browser tracks showing ATAC-seq data from different cell numbers. (C) Homer peak annotation analysis. Each peak was classified as a promoter-proximal or distal peak by Homer29. <a href="https://www-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/64242/64242fig06large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Name</td>
			<td>Sequence</td>
		</tr>
		<tr>
			<td>ESR1 C Forward</td>
			<td>TGGTGACTCATATTTGAACAAGCC</td>
		</tr>
		<tr>
			<td>ESR1 C Reverse</td>
			<td>CTCCTCCGTTGAATGTGTCTCC</td>
		</tr>
		<tr>
			<td>ESR1 K Forward</td>
			<td>TGTGGCTGGCTGCGTATG</td>
		</tr>
		<tr>
			<td>ESR1 K Reverse</td>
			<td>TGTCTCTCTTTCTGTTTGATTCCC</td>
		</tr>
		<tr>
			<td>ESR1 L Forward</td>
			<td>TGTGCCTGGAGTGATGTTTAAG</td>
		</tr>
		<tr>
			<td>ESR1 L Reverse</td>
			<td>CATTACAAAGGTGCTGGAGGAC</td>
		</tr>
	</tbody>
</table>
Table 1: List of primers used in this study.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Steps</td>
			<td>Temperature</td>
			<td>Time</td>
			<td>Cycle</td>
		</tr>
		<tr>
			<td>Gap filling</td>
			<td>72 °C</td>
			<td>5 min</td>
			<td>1</td>
		</tr>
		<tr>
			<td>Initial Denaturation</td>
			<td>98 °C</td>
			<td>30 s</td>
			<td>1</td>
		</tr>
		<tr>
			<td>Denaturation</td>
			<td>98 °C</td>
			<td>10 s</td>
			<td rowspan="3">5</td>
		</tr>
		<tr>
			<td>Annealing</td>
			<td>63 °C</td>
			<td>30 s</td>
		</tr>
		<tr>
			<td>Extension</td>
			<td>72 °C</td>
			<td>1 min</td>
		</tr>
		<tr>
			<td>Hold</td>
			<td>4 °C</td>
			<td>forever</td>
			<td></td>
		</tr>
	</tbody>
</table>
Table 2: PCR amplification program part 1.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Steps</td>
			<td>Temperature</td>
			<td>Cycle</td>
		</tr>
		<tr>
			<td>Initial Denaturation</td>
			<td>98 °C 30 s</td>
			<td>1</td>
		</tr>
		<tr>
			<td>Denaturation</td>
			<td>98 °C 10 s</td>
			<td rowspan="4">20</td>
		</tr>
		<tr>
			<td>Annealing</td>
			<td>63 °C 30 s</td>
		</tr>
		<tr>
			<td>Extension</td>
			<td>72 °C 1 min</td>
		</tr>
		<tr>
			<td>Plate read</td>
			<td></td>
		</tr>
	</tbody>
</table>
Table 3: qPCR settings.
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td>Steps</td>
			<td>Temperature</td>
			<td>Time</td>
			<td>Cycle</td>
		</tr>
		<tr>
			<td>Initial Denaturation</td>
			<td>98 °C</td>
			<td>30 s</td>
			<td>1</td>
		</tr>
		<tr>
			<td>Denaturation</td>
			<td>98 °C</td>
			<td>10 s</td>
			<td rowspan="3">Total cycles calculated at step 6.4</td>
		</tr>
		<tr>
			<td>Annealing</td>
			<td>63 °C</td>
			<td>30 s</td>
		</tr>
		<tr>
			<td>Extension</td>
			<td>72 °C</td>
			<td>1 min</td>
		</tr>
		<tr>
			<td>Hold</td>
			<td>4 °C</td>
			<td>forever</td>
			<td></td>
		</tr>
	</tbody>
</table>
Table 4: PCR amplification program part 2.

Watch this Scientific Journal Video about ATAC-Seq Optimization for Cancer Epigenetics Research at JoVE.com

ATAC-Seq Optimization for Cancer Epigenetics Research

ATAC-seq is a DNA sequencing method that uses the hyperactive mutant transposase, Tn5, to map changes in chromatin accessibility mediated by transcription factors. ATAC-seq enables the discovery of the molecular mechanisms underlying phenotypic alterations in cancer cells. This protocol outlines optimization procedures for ATAC-seq in epithelial cell types, including cancer cells.

The assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) probes deoxyribonucleic acid (DNA) accessibility using the ...

atac-seq-optimization-for-cancer-epigenetics-research

Nanyang Technological University

Research

JoVE Journal

Cancer Research

4.1K Views.  University of North Dakota School of Medicine and Health Sciences. ATAC-seq stelt ons in staat om de status van de open chromatinegebieden in het genoom te identificeren, die vaak wordt gewijzigd in de ziektetoestanden in vergelijking met de normale toestand. Minder celinvoer, een vereenvoudigd protocol van Tn5-vertering, gevolgd door isolatie van gefragmenteerd DNA en bibliotheekvoorbereiding door PCR-amplificatie en kortere experimentele tijd maken het een voordeligere techniek. Vergelijking van het veranderde epigenetische landschap door ATAC-seq tussen n...

Video: ATAC-Seq Optimization for Cancer Epigenetics Research

Mapping Genome-wide Accessible Chromatin in Primary Human T Lymphocytes by ATAC-Seq

Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method used for the identification of open (accessible) regions of chromatin. These regions represent regulatory DNA elements (e.g., promoters, enhancers, locus control regions, insulators) to which transcription factors bind. Mapping the accessible chromatin landscape is a powerful approach for uncovering active regulatory elements across the genome. This information serves as an unbiased approach for discovering the network of relevant transcription factors and mechanisms of chromatin structure that govern gene expression programs. ATAC-seq is a robust and sensitive alternative to DNase I hypersensitivity analysis coupled with next-generation sequencing (DNase-seq) and formaldehyde-assisted isolation of regulatory elements (FAIRE-seq) for genome-wide analysis of chromatin accessibility and to the sequencing of micrococcal nuclease-sensitive sites (MNase-seq) to determine nucleosome positioning. We present a detailed ATAC-seq protocol optimized for human primary immune cells i.e. CD4+ lymphocytes (T helper 1 (Th1) and Th2 cells). This comprehensive protocol begins with cell harvest, then describes the molecular procedure of chromatin tagmentation, sample preparation for next-generation sequencing, and also includes methods and considerations for the computational analyses used to interpret the results. Moreover, to save time and money, we introduced quality control measures to assess the ATAC-seq library prior to sequencing. Importantly, the principles presented in this protocol allow its adaptation to other human immune and non-immune primary cells and cell lines. These guidelines will also be useful for laboratories which are not proficient with next-generation sequencing methods.

Assay for Transposase-Accessible Chromatin coupled with high-throughput sequencing (ATAC-seq) is a genome-wide method to uncover accessible chromatin. This is a step-by-step ATAC-seq protocol, from molecular to the final computational analysis, optimized for human lymphocytes (Th1/Th2). This protocol can be adopted by researchers without prior experience in next-generation sequencing methods.

Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) is a method used for the identification of open (accessible) regions ...

Genetics

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture 4C-seq

The identification of physical interactions between genes and regulatory elements is challenging but has been facilitated by chromosome conformation capture methods. This modification to the 4C-seq protocol mitigates PCR bias by minimizing over-amplification of PCR templates and maximizes the mappability of reads by incorporating an addition restriction enzyme digest step.

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and ...

ATAC-seq Assay with Low Mitochondrial DNA Contamination from Primary Human CD4+ T Lymphocytes

ATAC-seq has become a widely used methodology in the study of epigenetics due to its rapid and simple approach to mapping genome-wide accessible chromatin. In this paper we present an improved ATAC-seq protocol that reduces contaminating mitochondrial DNA reads. While previous ATAC-seq protocols have struggled with an average of 50% contaminating mitochondrial DNA reads, the optimized&#160;lysis buffer introduced in this protocol reduces mitochondrial DNA contamination to an average of 3%. This improved ATAC-seq protocol allows for a near 50% reduction in the sequencing cost. We demonstrate how these high-quality ATAC-seq libraries can be prepared from activated CD4+ lymphocytes, providing step-by-step instructions from CD4+ lymphocyte isolation from whole blood through data analysis. This improved ATAC-seq protocol has been validated in a wide range of cell types and will be of immediate use to researchers studying chromatin accessibility.

Here, we present a protocol to perform an assay for transposase-accessible chromatin sequencing (ATAC-seq) on activated CD4+ human lymphocytes. The protocol has been modified to minimize contaminating mitochondrial DNA reads from 50% to 3% through the introduction of a new lysis buffer.

ATAC-seq has become a widely used methodology in the study of epigenetics due to its rapid and simple approach to mapping genome-wide accessible ...

Adipocyte-Specific ATAC-Seq with Adipose Tissues Using Fluorescence-Activated Nucleus Sorting

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a robust technique that enables genome-wide chromatin accessibility profiling. This technique has been useful for understanding the regulatory mechanisms of gene expression in a range of biological processes. Although ATAC-seq has been modified for different types of samples, there have not been effective modifications of ATAC-seq methods for adipose tissues. Challenges with adipose tissues include the complex cellular heterogeneity, large lipid content, and high mitochondrial contamination. To overcome these problems, we have developed a protocol that allows adipocyte-specific ATAC-seq by employing fluorescence-activated nucleus sorting with adipose tissues from the transgenic reporter Nuclear tagging and Translating Ribosome Affinity Purification (NuTRAP) mouse. This protocol produces high-quality data with minimal wasted sequencing reads while reducing the amount of nucleus input and reagents. This paper provides detailed step-by-step instructions for the ATAC-seq method validated for the use of adipocyte nuclei isolated from mouse adipose tissues. This protocol will aid in the investigation of chromatin dynamics in adipocytes upon diverse biological stimulations, which will allow for novel biological insights.

We present a protocol for assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) specifically on adipocytes using nucleus sorting with adipose tissues isolated from transgenic reporter mice with nuclear fluorescence labeling.

Assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) is a robust technique that enables genome-wide chromatin ...

CATCH-UP: A High-Throughput Upstream-Pipeline for Bulk ATAC-Seq and ChIP-Seq Data

Assay for transposase-accessible chromatin (ATAC) and chromatin immunoprecipitation (ChIP), coupled with next-generation sequencing (NGS), have revolutionized the study of gene regulation. A lack of standardization in the analysis of the highly dimensional datasets generated by these techniques has made reproducibility difficult to achieve, leading to discrepancies in the published, processed data. Part of this problem is due to the diverse range of bioinformatic tools available for the analysis of these types of data. Secondly, a number of different bioinformatic tools are required sequentially to convert raw data into a fully processed and interpretable output, and these tools require varying levels of computational skills. Furthermore, there are many options for quality control that are not uniformly employed during data processing. We address these issues with a complete assay for transposase-accessible chromatin sequencing (ATAC-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) upstream pipeline (CATCH-UP), an easy-to-use, Python-based pipeline for the analysis of bulk ChIP-seq and ATAC-seq datasets from raw fastq files to visualizable bigwig tracks and peaks calls. This pipeline is simple to install and run, requiring minimal computational knowledge. The pipeline is modular, scalable, and parallelizable on various computing infrastructures, allowing for easy reporting of methodology to enable reproducible analysis of novel or published datasets.

ATAC-seq and ChIP-seq allow detailed investigation of gene regulation; however, processing these data types is challenging and often inconsistent between research groups. We present CATCH-UP: an easy-to-use computational pipeline that allows standardized and reproducible data processing and analysis of new and published ATAC/ChIP-seq datasets.

Assay for transposase-accessible chromatin (ATAC) and chromatin immunoprecipitation (ChIP), coupled with next-generation sequencing (NGS), have ...

In Vitro Selection of Engineered Transcriptional Repressors for Targeted Epigenetic Silencing

Gene inactivation is instrumental to study gene function and represents a promising strategy for the treatment of a broad range of diseases. Among traditional technologies, RNA interference suffers from partial target abrogation and the requirement for life-long treatments. In contrast, artificial nucleases can impose stable gene inactivation through induction of a DNA double strand break (DSB), but recent studies are questioning the safety of this approach. Targeted epigenetic editing via engineered transcriptional repressors (ETRs) may represent a solution, as a single administration of specific ETR combinations can lead to durable silencing without inducing DNA breaks.
ETRs are proteins containing a programmable DNA-binding domain (DBD) and effectors from naturally occurring transcriptional repressors. Specifically, a combination of three ETRs equipped with the KRAB domain of human ZNF10, the catalytic domain of human DNMT3A and human DNMT3L, was shown to induce heritable repressive epigenetic states on the ETR-target gene. The hit-and-run nature of this platform, the lack of impact on the DNA sequence of the target, and the possibility to revert to the repressive state by DNA demethylation on demand, make epigenetic silencing a game-changing tool. A critical step is the identification of the proper ETRs&#39; position on the target gene to maximize on-target and minimize off-target silencing. Performing this step in the final ex vivo or in vivo preclinical setting can be cumbersome.
Taking the CRISPR/catalytically dead Cas9 system as a paradigmatic DBD for ETRs, this paper describes a protocol consisting of the in vitro screen of guide RNAs (gRNAs) coupled to the triple-ETR combination for efficient on-target silencing, followed by evaluation of the genome-wide specificity profile of top hits. This allows for reduction of the initial repertoire of candidate gRNAs to a short list of promising ones, whose complexity is suitable for their final&#160;evaluation in the therapeutically relevant setting of interest.

In Vitro Selection of Engineered Transcriptional Repressors for Targeted Epigenetic Silencing

Here, we present a protocol for the in vitro selection of engineered transcriptional repressors (ETRs) with high, long-term, stable, on-target silencing efficiency and low genome-wide, off-target activity. This workflow allows for reducing an initial, complex repertoire of candidate ETRs to a short list, suitable for further evaluation in therapeutically relevant settings.

Gene inactivation is instrumental to study gene function and represents a promising strategy for the treatment of a broad range of diseases. Among ...

Bioengineering

Detecting Somatic Genetic Alterations in Tumor Specimens by Exon Capture and Massively Parallel Sequencing

Efforts to detect and investigate key oncogenic mutations have proven valuable to facilitate the appropriate treatment for cancer patients. The establishment of high-throughput, massively parallel "next-generation" sequencing has aided the discovery of many such mutations. To enhance the clinical and translational utility of this technology, platforms must be high-throughput, cost-effective, and compatible with formalin-fixed paraffin embedded (FFPE) tissue samples that may yield small amounts of degraded or damaged DNA. Here, we describe the preparation of barcoded and multiplexed DNA libraries followed by hybridization-based capture of targeted exons for the detection of cancer-associated mutations in fresh frozen and FFPE tumors by massively parallel sequencing. This method enables the identification of sequence mutations, copy number alterations, and select structural rearrangements involving all targeted genes. Targeted exon sequencing offers the benefits of high throughput, low cost, and deep sequence coverage, thus conferring high sensitivity for detecting low frequency mutations.

We describe the preparation of barcoded DNA libraries and subsequent hybridization-based exon capture for detection of key cancer-associated mutations in clinical tumor specimens by massively parallel "next generation" sequencing. Targeted exon sequencing offers the benefits of high throughput, low cost, and deep sequence coverage, thus yielding high sensitivity for detecting low frequency mutations.

Efforts  to detect and investigate key oncogenic mutations have proven valuable to  facilitate the appropriate treatment for cancer patients. The ...

Biology

Multiplexed Analysis of Retinal Gene Expression and Chromatin Accessibility Using scRNA-Seq and scATAC-Seq

Powerful next generation sequencing techniques offer robust and comprehensive analysis to investigate how retinal gene regulatory networks function during development and in disease states. Single-cell RNA sequencing allows us to comprehensively profile gene expression changes observed in retinal development and disease at a cellular level, while single-cell ATAC-Seq allows analysis of chromatin accessibility and transcription factor binding to be profiled at similar resolution. Here the use of these techniques in the developing retina is described, and MULTI-Seq is demonstrated, where individual samples are labeled with a modified oligonucleotide-lipid complex, enabling researchers to both increase the scope of individual experiments and substantially reduce costs.

Here, the authors showcase the utility of MULTI-seq for phenotyping and subsequent paired scRNA-seq and scATAC-seq in characterizing the transcriptomic and chromatin accessibility profiles in retina.

Powerful next generation sequencing techniques offer robust and comprehensive analysis to investigate how retinal gene regulatory networks function during ...

Neuroscience

Scalable ATAC-seq Protocol for Chromatin Accessibility in Rodent Neutrophils

ATAC-Seq Library Preparation of Murine Bone Marrow-Derived Neutrophils

Assay for Transposase-Accessible Chromatin with sequencing (ATAC-seq) is a powerful, high-throughput technique for assessing chromatin accessibility and understanding epigenomic regulation. Neutrophils, as a crucial leukocyte type in immune responses, undergo substantial chromatin architectural changes during differentiation and activation, which significantly impact the gene expression necessary for their functions. ATAC-seq has been instrumental in uncovering key transcription factors in neutrophil maturation, revealing pathogen-specific epigenomic signatures, and identifying therapeutic targets for autoimmune diseases. However, neutrophils' sensitivity to the external milieu complicates high-quality ATAC-seq data production. Here, we propose a scalable protocol for preparing ATAC-seq libraries from rodent bone marrow-derived neutrophils, featuring improved immunomagnetic separation to ensure optimal cell viability and high-quality libraries. The vital elements impacting the library quality and optimization principles for methodological extension are discussed in detail. This protocol will support the researchers who are willing to study the chromatin architecture and epigenomic reprogramming of neutrophils, advancing studies in basic and clinical immunology.

This manuscript outlines a protocol for preparing an ATAC-seq library of neutrophils from murine bone marrow, aiming to guarantee optimal neutrophil viability and high library quality. It offers step-by-step instructions on BMC preparation, immunomagnetic sorting, and library construction, serving as a valuable guide, especially for newcomers to studying neutrophils.

Assay for Transposase-Accessible Chromatin with sequencing (ATAC-seq) is a powerful, high-throughput technique for assessing chromatin accessibility and ...

Immunology and Infection

3' End Sequencing Library Preparation with A-seq2

Studies in the last decade have revealed a complex and dynamic variety of pre-mRNA cleavage and polyadenylation reactions. mRNAs with long 3&#39; untranslated regions (UTRs) are generated in differentiated cells whereas proliferating cells preferentially express transcripts with short 3&#39;UTRs. We describe the A-seq protocol, now at its second version, which was developed to map polyadenylation sites genome-wide and study the regulation of pre-mRNA 3&#39; end processing. Also this current protocol takes advantage of the polyadenylate (poly(A)) tails that are added during the biogenesis of most mammalian mRNAs to enrich for fully processed mRNAs. A DNA adaptor with deoxyuracil at its fourth position allows the precise processing of mRNA 3&#39; end fragments for sequencing. Not including the cell culture and the overnight ligations, the protocol requires about 8 h hands-on time. Along with it, an easy-to-use software package for the analysis of the derived sequencing data is provided. A-seq2 and the associated analysis software provide an efficient and reliable solution to the mapping of pre-mRNA 3&#39; ends in a wide range of conditions, from 106 or fewer cells.

This protocol describes a method for mapping pre-mRNA 3&#39; end processing sites.

ATAC-Seq Optimization for Cancer Epigenetics Research

Samenvatting

Meer video's verkennen