Dec 08, 2021

Public workspaceATAC-STARR-seq V.2

  • 1Vanderbilt University
Icon indicating open access to content
QR code linking to this content
Protocol CitationTyler Hansen, Emily Hodges 2021. ATAC-STARR-seq. protocols.io https://dx.doi.org/10.17504/protocols.io.b2nuqdewVersion created by Tyler Hansen
Manuscript citation:
Hansen, T. & Hodges, E. Simultaneous profiling of regulatory activity, chromatin accessibility, and transcription factor occupancy with ATAC-STARR-seq. In Revision. 2020.
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License,  which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it’s working
Created: December 06, 2021
Last Modified: December 08, 2021
Protocol Integer ID: 55732
Disclaimer
DISCLAIMER – FOR INFORMATIONAL PURPOSES ONLY; USE AT YOUR OWN RISK

The protocol content here is for informational purposes only and does not constitute legal, medical, clinical, or safety advice, or otherwise; content added to protocols.io is not peer reviewed and may not have undergone a formal approval of any kind. Information presented in this protocol should not substitute for independent professional judgment, advice, diagnosis, or treatment. Any action you take or refrain from taking using or relying upon the information presented here is strictly at your own risk. You agree that neither the Company nor any of the authors, contributors, administrators, or anyone else associated with protocols.io, can be held responsible for your use of the information contained in or linked to this protocol or any of our Sites/Apps and Services.
Abstract
Massively parallel reporter assays test the capacity of putative cis-regulatory elements (CREs) to drive transcription on a genome-wide scale. In nearly all cases, chromatin accessibility is necessary to drive activity, so most CREs are inactive due to chromatin context rather than intrinsic DNA sequence properties. Here, we combined assay for transposase-accessible chromatin (ATAC-seq) with self-transcribing active regulatory region sequencing (STARR-seq) to selectively assay the regulatory potential of nucleosome-free DNA genome-wide. Our approach enabled high-resolution testing of ~50 million unique DNA fragments tiling ~101,000 accessible chromatin regions in human lymphoblastoid cells. To illustrate the application of our approach, we show that 30% of all accessible regions contain an activator, a silencer or both. Benchmarking against standard ATAC-seq, our approach faithfully captures chromatin accessibility and transcription factor (TF) footprints with high signal-to-noise. Integrating three layers of genomic information (accessibility, TF occupancy, and activity) provided by ATAC-STARR-seq, we stratified active and silent CREs by the presence of several TF footprints and show that CREs with specific TF combinations are associated with distinct gene regulatory pathways. Altogether, these data highlight the power of ATAC-STARR-seq to comprehensively investigate the regulatory landscape of the human genome from a single DNA source.
Materials
*Oligo sequences available in Hansen & Hodges NAR 2021.

* Plasmids
  • hSTARR-Seq_ORI (Addgene #99296)
  • pcDNA3.1-eGFP (Addgene #13031)

*Special Equipment:
  • Gene Pulser Xcell Microbial System (Bio-Rad, #1652662) – the Cortez Lab has one.
  • Neon™ Transfection System (Invitrogen, #MPK5000)
  • 4150 TapeStation System (Agilent, #G2992AA)
  • DynaMag™-PCR Magnet (Invitrogen, #492025)
  • DynaMag™-2 Magnet (Invitrogen, #12321D)

*Buffer Recipes:
  • STE Buffer - (10 mM Tris-HCl pH 8.0, 50 mM NaCl, 1mM EDTA, pH 8.0)
  • Tn5 Storage Buffer - (50 mM Hepes pH 7.2, 100 mM NaCl, 0.1 mM EDTA, 1mM DTT,
0.1% Triton-X 100, 50% Glycerol)
  • ATAC-Resuspension Buffer - (10mM Tris-HCl pH 7.4, 10mM NaCl, 3mM MgCl2)
  • 5x Tris-DMF - (50 mM Tris-HCl pH 7.5, 25 mM MgCl2, 50% DMF). Pipette DMF with a glass pipette.
  • Flow cytometry buffer – (1xPBS + 1%BSA)
  • Recovery Media – RPMI 1640, 20% FBS
  • Maintenance Media – RPMI 1640, 15% FBS, 2mM GlutaMax, 1X Penn/Strep
  • Binding Buffer – 20 mM Tris-HCl, pH 7.5, 1.0 M LiCl, 2 mM EDTA
  • Wash Buffer A – 10 mM Tris-HCl, pH 7.5 ,0.15 M LiCl, 1 mM EDTA 10 mM Tris-HCl, pH 7.5

*Reagents and consumables:
  • 10% IGEPAL CA-630 – (Roche cat# 11332473001) Also called NP-40 substitute, IGEPAL CA-630 is supplied at 10%. Store at 4°C.
  • Digitonin (Cayman Chemical, #14952) – Prepared Digitonin at 2% in DMSO (20mg/mL). Dilute 1:1 with water to make a 1% (100x) stock solution. Avoid more than 5 freeze thaw cycles of 2% solution. 1% can be kept at -20°C for up to 6 monthsand does not freeze at -20°C.
  • Tween-20 – (Roche cat# 11332465001) Tween-20 is supplied at 10%. Store at 4°C.
  • DNA Clean & Concentrate-5 (Zymo Research, #D4004)
  • NEBNext® High-Fidelity 2X PCR Master Mix (M0541L)
  • High Sensitivity D5000 ScreenTape (Agilent, #5067-5592)
  • High Sensitivity D5000 Reagents (Agilent, #5067-5593)
  • SPRIselect (Beckman-Coulter, #B23317)
  • DynaMag™-PCR Magnet (Invitrogen, #492025)
  • NEBNext® Ultra™ II Q5® Master Mix (NEB, #M0544L)
  • NEBuilder® HiFi DNA Assembly Master Mix (NEB, #E2621L)
  • GlycoBlue™ Coprecipitant (15 mg/mL) (Invitrogen, # AM9515)
  • MegaX DH10B™ T1R Electrocomp™ Cells (Invitrogen, # C640003)
  • Gene Pulser/MicroPulser Electroporation Cuvettes, 0.1 cm gap (Bio-Rad, #1652089), these can be reused—clean with 70% EtOH and dry overnight.
  • ZymoPURE™ II Plasmid Gigaprep Kit (Zymo Research, #D4204)
  • Neon™ Transfection System 100 µL Kit (Invitrogen, #MPK10096)
  • Propidium iodide solution (1mg/mL in water) (Sigma-Aldrich, #P4864)
  • TRIzol™ Reagent and Phasemaker™ Tubes Complete System (Invitrogen™, #A33251)
  • ZymoPURE™ II Plasmid Midiprep Kit (Zymo Research, #D4200)
  • Oligo d(T)25 Magnetic Beads (NEB, #S1419S)
  • DNaseI (NEB, #M0303S)
  • Zymo Research RNA Clean & Concentrator-25 kit (Zymo Research, #R1018)
  • PrimeScript™ Reverse Transcriptase (Takara, #2680)
  • KAPA dNTP Mix (10 mM each) (Roche, #KK1017)
  • RNaseOUT™ Recombinant Ribonuclease Inhibitor (Invitrogen, #10777019). We have also used RiboLock RNase Inhibitor (Thermo Scientific, #EO0382).
  • RNase A, DNase and protease-free (10 mg/mL) (Thermo Scientific, #EN0531)
Anneling Oligos
Anneling Oligos
Dilute oligos to 100 µM in STE buffer.
Make Adaptor Mixes in 200μL PCR tubes:

Mix (Adaptor A) [50µL]:
5µL Tn5MERV oligo (100µM)
5µL Tn5_1 Oligo (100µM)
40µL H2O

Mix (Adaptor B) [50µL]:
5µL Tn5MERV oligo (100µM)
5µL Tn5_2_ME_Comp Oligo (100µM)
40µL H2O
Anneal Adaptor solutions separately with thermocycler conditions as follows:
95°C, 3 min
65°C, 3 min
Ramp to 24°C, -1ºC per cycle, 30 s
Hold at 24°C forever
After annealing mix:
50µL Adaptor solution A
50µL Adaptor solution B
100µL Glycerol (≥99.5%)
200 µL Adaptor mixture, which is now 5 µM.
Transposome Assembly
Transposome Assembly
Mix the following: 10 µL 5µM Adaptors Mix 10 µL Tn5 (3mg/mL in Tn5 storage buffer) 20 µL Total Reaction Volume
Let stand at room temp for 30-60min, then immediately place on ice. Assembled transposome may be stored at -20°C for several months.
ATAC Reaction
ATAC Reaction
Turn on all equipment!

Set the microcentrifuge to 4ºC & the thermomixer to 37ºC. Each of these takes 15-30min to get to their desired temp.

Count cells in duplicate using a hemocytometer. Record cell count.

Note: It is important that cell viability is greater than 85%. If less than 85%, I strongly recommend trying again when viability is much better.
In two tubes, pellet 200,000 viable cells in each tube at 500 RCF at 4°C for 10 min in a fixed angle centrifuge.

Use low-bind Eppendorf tubes to prevent pellets from sticking along the back of the tube (this is much more important for reproducibility than you might think).
While waiting, make lysis buffers.
ATAC-Resuspension Buffer (RSB) containing 0.1% NP40, 0.1% Tween-20, and 0.01% digitonin (lysis buffer):

For 110μL, mix 1.1μL 10% Tween-20,1.1μL 10% IGEPAL CA-630, 1.1μL 1% digitonin, 106.7μL 1x RSB.
ATAC-RSB containing only 0.1% Tween-20 (wash out buffer):

For 2.2mL, mix 220μL 10% Tween-20 and 2mL 1x RSB.

Keep on ice.
Aspirate supernatant and resuspend each cell pellet in 50µL cold ATAC-Resuspension Buffer (RSB) containing 0.1% IGEPAL CA-630, 0.1% Tween-20, and 0.01% Digitonin and pipette up and down 5 times.

Note: The volume used here is a modification to the Omin-ATAC protocol. In my hands, 50μL for this number of cells works just as well as 200μL and is much more manageable.
Incubate on ice for 3 minutes.
Wash out lysis with 1 ml of cold ATAC-RSB containing ONLY 0.1% Tween-20 and invert tube 3 times to mix.

Note: The volume used here is a modification to the Omin-ATAC protocol. In my hands, 1mL for this number of cells works just as well as 4mL and is much more manageable.
Immediately pellet nuclei at 500 RCF for 10min at 4°C in a fixed angle centrifuge. While waiting, make transposition mix:

Aspirate all supernatant, carefully avoiding visible cell pellet, using two pipetting steps (aspirate down to 100 ul with a p1000 pipette and remove final 100 ul with a p200 pipette).
Resuspend cell pellets in 180µL of transposition mix by pipetting up and down 10 times. Combine like samples for a total volume of 360μL.
Add 40µL assembled Tn5 transposase and quickly vortex to mix well. This is the equivalent to 5μL/reaction.
Divide resuspension into 8-1.5mL tubes (50μL/tube).
Incubate reaction at 37°C for 30min in a thermomixer with 1000 RPM mixing.
After incubation, pool the eight like-reactions together in a 15mL conical tube containing 2mL DNA binding Buffer.
Purify the small-isolated DNA fragments using the Zymo C+C kit. Run all 2.4mL through one column, adding in 800μL increments. Elute in 22µl pre-warmed elution buffer. Incubate on column for 5min.

Can stop here and store at -20°C.
PCR Amplification
PCR Amplification
Make PCR mastermix:

Note: do not use a hot-start polymerase here. The reaction requires polymerase activity in the initial extension step.


Run PCR with the following thermocycler parameters:

72 ºC, 5 min
98ºC, 30 sec
4 cycles of
98ºC, 10 sec
62ºC, 30 sec
72ºC, 1 min
Final extension 72ºC, 2 min
Purify PCR Product using Zymo Clean and Concentrator Kit. Mix 250μL Binding Buffer with 50μL PCR products. Elute in 32µl pre-warmed elution buffer. Incubate on column for 5min.
Analyze by running 2uL on High-Sensitivity D5000 Tapestation. Determine size distribution and molarity.

If it looks good, move on to cloning. Need between 0.04pmols and 0.27pmols total—no more, no less. If there are too many large fragments, perform the SPRI bead size selection below.
SPRI bead size-selection
SPRI bead size-selection
*Note: Perform a 0.6X ratio right-side selection to remove fragments larger than 500bp. The volumes listed below are for that ratio specifically. Reference the SPRISelect User Manual if this needs to be altered.

Thoroughly shake the SPRIselect bottle to resuspend the SPRI beads.
Add 18μL of SPRIselect to the sample. Pipette up and down and vortex for at least 5 seconds to mix thoroughly.

Note: Insufficient mixing of sample and SPRIselect will lead to inconsistent size selection results.
Place the reaction vessel on DynaMag-PCR magnet and allow the SPRI beads to settle to the magnet for ~1min.
Transfer the clear supernatant, which contains the Right-Side Size Selected sample, to a new reaction vessel. The reaction vessel with the remaining beads can be discarded.
Add 36μL of SPRIselect to the Right-Side Size Selected sample. Pipette up and down and vortex for at least 5 seconds to mix thoroughly.
Place the reaction vessel on an appropriate magnetic stand or plate and allow the SPRI beads to settle to the magnet for ~1min.
Remove and discard the clear supernatant.
With the reaction vessel still on the magnet, add 180 μL of freshly-made 85% ethanol and incubate at RT for 30 seconds. Remove and discard ALL of the ethanol supernatant.
Let the beads dry but not over dry.
Once cracks are starting to form in the pellet, elute by removing the reaction vessel from the magnet and adding 12μL of 10mM Tris-HCl pH 8.0, elution buffer, or nuclease-free water.
Pipette up and down and vortex for at least 5 seconds to mix thoroughly.
Place the reaction vessel on an appropriate magnetic stand or plate and allow the SPRI beads to settle to the magnet for ~1min.
Transfer the eluate (size selected sample) to an appropriate storage vessel.

Note: It is very important to avoid pipetting beads over to the new tube. They can cause the tape station to run inaccurately.
Analyze size-selected sample by running on 2uL on High-Sensitivity D5000 Tapestation and analyze for concentration and size range.
Linearization of hSTARR-Seq_ORI
Linearization of hSTARR-Seq_ORI
*Note: at this step, you will generate the vector required for the gibson cloning reaction. This linearized vector contains the homology arms for gibson. In the NAR paper we used an oligo to add the N504 barcode to the plasmid, this is not required for most purposes. For this reason, I’ve denoted the reverse primer as N5XX/universal because you can add whatever index you want OR not add an index at all, just like the forward primer.


PCR amplify vector backbone from the hSTARR-Seq_ORI plasmid (Addgene #99296) with the mixture and thermocycler conditions below.

98ºC, 30 sec
2 cycles of
98ºC, 10 sec
59ºC, 30 sec
72ºC, 90 sec
28 cycles of
98ºC, 10 sec
72ºC, 2min
Final extension 72ºC, 2 min
Purify PCR Product using Zymo Clean and Concentrator – 5 Kit. Elute in 20μL elution buffer.
Determine concentration via Nanodrop.
Analyze by running 0.1ng on 1% agarose gel with EtBr stain to determine purity.
Clone fragments into linearized hSTARR-Seq_ORI via Gibson Cloning
Clone fragments into linearized hSTARR-Seq_ORI via Gibson Cloning
*This part of the protocol is intentionally flexible in order to add the entire volume of tagments to these gibson reactions. In this step, perform four Gibson reactions per sample at a 1:2 rato of vector:insert.

Note: this used to be a 1:3 ratio but was changed to 1:2 after learning that 1:2 would require fewer PCR cycles and is the recommended ratio for this enzyme cocktail.

The insert (tagments) should be at a total amount of 0.04pmols and 0.27pmols. Calculate the reaction molarity as follows:


In 4 PCR tubes, combine the volumes of insert and vector as written above.
Bring up to a total volume of 5μL with water.

Note: The total volume cannot exceed 5μL.
In a separate PCR tube, preparea negative control where you substitute the inserts for water (only 1 reaction).
Add 5μL of NEBuilder Hifi DNA Assembly 2x MasterMix and mix well.
Incubate at 50°C for 1hr.
Ethanol Precipitation of Gibson Reactions
Ethanol Precipitation of Gibson Reactions
*This is needed to prevent arching in the Electroporation step. DO NOT SKIP.

Pool all four like-reactions together into a low-bind Eppendorf tube.
Adjust volume to 222.5μL by adding 182.5μL elution buffer (212.5μL to negative control).
Add 25μL 3M NaAc pH5.2 and Vortex
Add 2.5μL GlycoBlue and Vortex.
Add 750μL cold (-20°C) 100% EtOH and vortex.
Incubate at -20°C for 30-120min.
Spin at full speed for 30min @ 4°C
Aspirate supernatant
Wash with 750μL cold 70% EtOH. Spin at full speed for 15min @ 4°C.
Aspirate supernatant and repeat wash.
Dry pellet at room temp until dry.
Resuspend pellet in 10μL 10mM Tris-HCl pH 8.0. Pipette up/down until reconstituted. Can stop here, store @ -20°C.
Electroporation of STARR-seq Plasmid Library into Bacterial Cells
Electroporation of STARR-seq Plasmid Library into Bacterial Cells
*In addition to the negative control, do a positive control using the supplied pUC19. For the controls, only do one reaction and use refrozen cells if possible.

Pre-cool pipette tips, GenePulser cuvettes, and 1.5mL tubes at -20°C. Prewarm 2L LB+amp flasks and SOC @ 37°C.
Add all 10μL of each gibson sample to a pre-cooled tube (2.5μL for controls). Keep on ice.
Thaw MegaX DH10B cells on ice. Each tube contains ~120μL.
Add 80μL of thawed cells to each 10μL gibson sample and 20μL to each 2.5μL control using pre-cooled tips. Keep on ice.
Pipette SLOWLY up/down to mix with pre-cooled tips. Then add 20μL (max volume) to a GenePulser 0.1cm cuvette—may need to use 3-4 cuvettes for each sample, but only one for each control. Keep cuvettes on ice.

Note: make sure the volume is touching both metal plates, otherwise it will arc (the best technique here is to pipette into the corner of the cuvette at an angle so the cells reach the bottom). To help, gently tap to the bottom.
Electroporate with “exponential decay” pulse type at 2kV, 200Ω, and 25μF - (Gap distance is 0.1cm for brown-capped cuvettes).
Immediately add 1mL pre-warmed SOC media to the 20μL DNA/bacteria mixture in the cuvette.
Combine like-reactions into a round bottom 14mL culture tube.
Recover at 37°C for 1 hour in shaker incubator.
After recovery make a dilution series:

Plate 100μL of each dilution onto a prewarmed LB+ampicillin(100μg/mL) plate. Incubate at 37°C for ~16hrs.
With the remaining undiluted culture, add directly to 1L pre-warmed liquid LB+ ampicillin(100μg/mL) media. Shake at 37°C O/N 225 RPM. Do not do this for the controls!

DO NOT FORGET TO ADD THE AMPICILLIN TO THE 1L CULTURE!!!
Quantification of colony counts to estimate number of transformants
Quantification of colony counts to estimate number of transformants
Count colonies on plates after overnight growth either manually or via automated means (imaging and particle counting with imageJ)

The equations below assume a plating of 100µL from a 1mL dilution preparation:
(.01 colony count) * 1000 = # colonies in 1mL of recovery media
(.001 colony count) * 10000 = # colonies in 1mL of recovery media
(.0001 colony count) * 100000 = # colonies in 1mL of recovery media
Purification of Large Scale Plasmid Library Growth
Purification of Large Scale Plasmid Library Growth
Record OD600.
Pellet 1L culture at 3400xg for 10min.
Proceed with Zymo Plasmid GigaPrep Kit protocol according to manufacturer directions.
Elute with 2mL ZymoPURE Elution Buffer. Nanodrop for purity and concentration. Store DNA at -20°C.
Low-Read Sequencing of Plasmid Library to Determine Insert Complexity
Low-Read Sequencing of Plasmid Library to Determine Insert Complexity
*Before proceeding to electroporation, it is important to determine whether the library has sufficient complexity. To do this, perform limited-cycle PCR to add illumina adapters and amplify the insert off the plasmid. Then submit the sample sequencing and target 15 million reads. Must be PE, but the cycle number is irrelevant as long as reads can be mapped.

* To reduce rxn bias effect, perform two PCRs for each sample.

Prepare a mastermix:

Aliquot 50μL mix into PCR tubes.
Vortex and centrifuge, briefly. Run the following PCR reaction:

98ºC, 30 sec
8 cycles of
98ºC, 10 sec
68ºC, 30 sec
72ºC, 30 sec
Final extension 72ºC, 2 min
Pool like-reactions and then purify using Zymo Clean and Concentrate - 5 kit – elute in 30μL elution buffer (warm buffer to 55°C or 37°C prior to elution and incubate on column for 5min).
Determine concentration and size distribution via HS D5000 ScreenTape on the Agilent Tapestation.

Concentration should be greater than 5nM but not much higher. Anything greater than 50nM is significantly overamplified.
Submit samples for paired end sequencing. In the NAR paper, we requested 25M reads at PE150 on NovaSeq.
Electroporate ATAC-STARR Libraries into Mammalian Cells
Electroporate ATAC-STARR Libraries into Mammalian Cells
*Each transfection reaction has a max capacity of 5 x106 cells and 10μL DNA. Do 20 reactions to reach 100x106 cells per sample.

*To assess transfection efficiency, conduct an eGFP electroporation and no DNA electroporation (as unstained control). For this, only do one reaction of 5 million cells. The eGFP plasmid is ~2x the size of the ATAC-STARR library plasmids, so to match the required molarity, increase the concentration 2-fold. This control should be done for each individual cell type used in the experiment.

Prepare DNA samples:
Add 110μg ATAC-STARR library DNA into a 1.5 mL Eppendorf tube in a max volume of 220μL [This equates to 5μg DNA per 5 million cells].
For GFP control, add 11μg of pcDNA3-eGFP plasmid (endotoxin-free) to a 1.5mL Eppendorf tube in a max volume of 11μL. [This equates to 10μg DNA per million cells].
For no DNA control, add 11μL water to a 1.5mL Eppendorf tube.
Prepare Suspension Cells:
Day 1: Perform cell count (with hemocytometer)
Per sample, add 50 mL of warm recovery medium (without antibiotics) into four T-75 culture flasks. Incubate at 37°C until use.
For each GFP or noDNA control reaction, add 10mL of media into a T-12.5 culture flasks.
Transfer 121 million cells to as many 50 mL falcon tube as necessary and centrifuge at 300xg for 5 minutes. Aspirate supernatant. Wash and resuspend with 45mL of DPBS or PBS, combine pellets. Centrifuge for 300xg for 5 min.
Set up the Neon Pipette Station:

a. Fill a fresh neon tube with 3 mL of pre-warmed Buffer E2. The electrode on the side of the tube should be completely submerged in buffer.
b. Insert tube into the pipette station until you hear a click.
c. Replace neon tube between samples—only need to change between experimental and controls (can use the same tube for water and GFP if doing water first).
Resuspend cells in 2178μL Buffer R/T. (R for GM12878, GM11831, and SUDHL-6 and T for LCL8664). Gently pipette for a single cell suspension. (Do not leave cells resuspended in Buffer R/T for longer than 15-30 min).
Electroporation:
Add 1980μL cells to each ATAC-STARR library DNA and gently mix.
Add 99μL cells to GFP plasmid. Gently mix.
Add 99μL cells to water. Gently mix.
Electroporate 20 ATAC-STARR reactions and both control reactions (100μL at a time) using the following conditions:

100μL Neon Tip, 1100V, 30ms, 2 pulses.

If the sample arcs, eject into waste and grab a new tip—19 good reactions are better than 19 good + 1 very dead.

Note: Each tip can be used 5 times. Reuse the tip only if doing an aliquot of the same sample. After use, save the tip and it can be regenerated to use for another application. Regenerated tips are somewhat less efficient and more prone to arcing than fresh tips, so do not use regenerated tips in this protocol. That said, these tips are quite expensive at ~$20 each, so keep that in mind when using.
After each electroporation, remove the pipette from the station and expel contents into one of the prepared T-75 (or T-12.5 for control reactions) flasks from step 6 above. Repeat electroporation with remaining aliquots and expel 5 aliquots into each T-75 flask containing 50mL media—cell concentration will be ~500,000 cells/mL. Gently swirl flask after addition of each transfected aliquot. Place in incubator at 37°C for 24 hours.
Check transformation efficiency via flow cytometry
Check transformation efficiency via flow cytometry
***We will use PI as a viability dye; PI will stain dead cells***

After ~24 hours, for both GFP and No DNA transfections, prepare cells for flow:
Spin down cells in 15mL conical at 300xg, 5min.
Remove supernatant and resuspend in 5mL 1x PBS or DPBS.
Spin at 300xg, 5min. Remove Supernatant.
Resuspend in 1000μL Flow Cytometry Buffer (1%BSA in 1xPBS) to achieve a cell density of ~1.25x106 cells/mL.
Split each sample into two aliquots of 500μL each. Pass through the flow cytometry tube cap filter to obtain single-cell suspensions.
Stain one from each with propidium iodine:

Make a 1:1 dilution of 1μg/μL PI solution with 1xPBS+1%BSA.

Add 1.25μL of diluted PI to the desired tubes. Vortex.
Analyze samples with 3-laser or 5-laser Flow Cytometer. Use no DNA as unstained compensation. No DNA+PI and GFP noPI as compensation controls. Analyze viable GFP+PI cells for GFP/no GFP.
Harvest RNA and plasmid DNA from ATAC-STARR transfected cells
Harvest RNA and plasmid DNA from ATAC-STARR transfected cells
Exactly 24 hours after electroporation (counting from last transfection added to media), mix all four T-75 flasks—200mL should all fit in one.
Take 10 µl sample from each flask and perform cell count/assess viability.
Divide into two 100mL volumes and process each separately.

Note: one is for RNA harvest, the other is for isolated plasmid DNA harvest. Process the RNA aliquot first as it is the most time-sensitive.
Before processing each, spin down in 50mL conical tubes, 300xg 5min. Remove sup and wash with 25mL 1x PBS or 1x DPBS. Spin down cells at 300xg, 5min and remove supernatant.
RNA Extraction:

Note: Trizol and Chloroform steps should be done in a chemical fume hood.
Add 5mL Trizol per sample (this equates to 1ml Trizol per 10e6 cells). Pipette up/down to homogenize. Note: can stop here and store at -20°C overnight.
Centrifuge 5 Invitrogen™ Phasemaker™ Tubes at 12,000xg for 30s.
Aliquot five volumes of 1mL into Invitrogen™ Phasemaker™ Tubes.
Incubate for 5 min at RT. Note: can stop here and store at -20°c overnight.
Add 0.2mL chloroform to each 1mL aliquot and mix vigorously by inversion.
Incubate for 2-3min at RT.
Spin at 12,000xg for 15min @4°C.
Transfer the aqueous layer to a 1.5mL Eppendorf tube—try to get ~600μL from each.
Add 2μL Glycoblue co-precipitant (15mg/ml) to each. Do not combine like-reactions yet.
Add 0.5mL isopropanol to each aqueous phase.
Incubate at RT for 10min.
Spin for 10min at 12,000xg @4°C.
Check for a white gel-like precipitate. If there, remove sup and continue with protocol. If no pellet spin again.
Resuspend Centrifuge 5 Invitrogen™ Phasemaker™ Tubesat 12,000xg for 30s. pellet in 1mL 75% EtOH. RNA can be stored in this solution at -20°C for up to a year.
Vortex briefly, and centrifuge 5min at 7500xg @4°C. Completely discard supernatant. It helps to use a P10 after briefly spinning.
Air-dry pellet for 5-10min.
Combine by suspending all 5 pellets in 100μL nuclease-free water. Incubate @ 55°C for 10min.
Determine RNA yield with Nanodrop (expected yield: 2.5-1.5μg/μL, 260/280 > 2, 260/230 > 1.8). RNA can be stored at -80°C overnight.
Plasmid DNA extraction:

*Use ZymoPURE II midiprep to isolate plasmid DNA
Follow manufactures protocol for kit. Add 8mL P1 buffer directly to cell pellet.
Elute in 50μL pre-warmed 10mM Tris pH 8.0. Incubate 5min on column.
Determine concentration and purity via nanodrop (expect ~300-400ng/μL). Store at -20°C. These samples will be amplified in a PCR reaction in the NGS PCR step.
Poly(A)+ RNA enrichment with oligo(dT)25 magnetic beads
Poly(A)+ RNA enrichment with oligo(dT)25 magnetic beads
*The magnetic beads are supplied at 5mg/mL and each mg can bind up to 5μg Poly(A)+ RNA. We will assume that Poly(A)+ RNA makes up about 5% (this is at the higher end because we are accounting for both mRNAs and the Reporter RNAs). Therefore, for 100μg Total RNA, we expect 5μg to be Poly(A)+ RNA. To capture all 5μg, we need 1mg of magnetic beads. Therefore, the ratio is 1μg Total RNA:10μg beads.

Determine the amount of beads to wash:


Chill Wash Buffer B at 4°C.
Wash magnetic beads:
Resuspend magnetic beads in the vial (vortex for 30s).







Transfer the desired amount to a 1.5mL Eppendorf tube (use more than one tube if volume > 750μL).
Add an equal volume of Binding Buffer and resuspend.
Place tube on DynaMag-2 magnet and incubate for 1min. Remove supernatant.
Remove tube from magnet and resuspend beads in the same volume of binding buffer used in step 98.3
Isolate Poly(A)+ RNA from Total RNA:
Add 100μL binding buffer to each 100μL Total RNA sample.
Heat to 65°C for 2min to disrupt secondary structures.
Immediately place on ice. Incubate 1min.
Add the indicated amount of washed beads to add to the 200μL Total RNA solution.

Note: Gently handle the beads for the remaining steps of the protocol.
Mix well and incubate on a rotator for 5min at 4°C.
Place tube on the magnet and incubate for ~30s. Remove supernatant.
Remove tube from magnet and add the same volume used in step iv of Washing Buffer B. Mix by pipetting up/down gently.
Rotate for 2-3min at 4°C.
Incubate on magnet for ~30s and remove supernatant.
Repeat wash for a total of 4 washes.
Add 82μL of 10mM Tris-HCl pH 7.5 to the beads.
Gently mix. Incubate at 75°C for 2min.
Place on magnet and quickly transfer eluate to a 1.5mL tube. Determine concentration via nanodrop (should be ~500-250ng/μL). Note: Can stop here and store at -70°C.
DNaseI digestion
DNaseI digestion
*20 units are sufficient to treat ~100μg RNA

Per sample, add 10μL DNaseI Reaction buffer and 10μL DNase I (RNase-free)(20 units).
Mix well and incubate @ 37°C for 10min.
Purify Poly(A)+ RNA with Zymo RNA Clean & Concentrator – 25 kit.

Volumes: 200μL RNA binding buffer, 300μL 100% EtOH

Elute in 60μL Nuclease-Free water.
Nanodrop to determine concentration.
Reverse Transcription
Reverse Transcription
For each 20μL reaction, the input maximum is 1μg Poly(A)+ RNA. We expect to have ~25μg, so scale up 25-fold. Do ten 50μL reactions.

Anneal primer to template RNA:
Make MasterMix following the recipe below:

Aliquot 25μL into ten 0.2mL PCR tubes.
Heat the RNA-primer mix at 65°C for 5min and then immediately incubate on ice for at least 1 min.
Prepare RT reaction mix:
Vortex and briefly centrifuge the 5x PrimeScript Buffer
Make MasterMix following the recipe below:

Add 25μL of the MasterMix to each of the 10 RNA-Primer reactions. Mix well.
Incubate at 42°C for 60min.
To inactivate the reaction, heat at 70°C for 15min and then cool on ice.
RNaseA Treatment
RNaseA Treatment
Pool like reactions for a total volume of ~500μL
Add 1μL RNAseA (10mg/mL) and incubate 37°C for 1hr.

Note: final concentration is ~20μg/mL in low salt conditions [KCl] = 75mM and [MgCl2] = 8mM. At low salt concentrations, RNase A cleaves single-stranded and double-stranded RNA as well the RNA strand in RNA-DNA hybrids.
Purify cDNA using Zymo DNA Clean and Concentrate-25 kit:

a.Add 3.5mL Binding Buffer (7:1). Mix in a 15mL conical. Spin through column in 800μL increments.
b.Wash with 200μL wash buffer twice.
c.Elute in 28μL elution buffer.
Nanodrop to determine yield and purity. Expected yield is ~10-50 ng/μL
Next-Generation sequencing PCR
Next-Generation sequencing PCR
* To reduce rxn bias effect, perform two PCRs for each sample.

Prepare a mastermix:

Aliquot 39μL mix into PCR tubes and add 10μL of sample to each.
Add 1μL of the correct 50μM P7 Nextera primer (N7XX).
Vortex and centrifuge, briefly. Run the following PCR reaction:

98ºC, 30 sec
8 or 12 cycles of àplasmid DNA = 8 cycles, reporter RNA = 13 cycles
98ºC, 10 sec
68ºC, 30 sec
72ºC, 30 sec
Final extension 72ºC, 2 min
Pool like-reactions and then purify cDNA using Zymo Clean and Concentrate -5 kit – elute in 30μL elution buffer (warm buffer to 55°C or 37°C prior to elution and incubate on column for 5min).
Determine concentration and size distribution via HS D5000 ScreenTape on the Agilent Tapestation. Concentration should be greater than 5nM but not much higher. Anything greater than 50nM is significantly overamplified.

Note: you will see the template and initial amplification products in the Tapestation, these products are not capable of generating reads on the sequencer.
Submit samples for paired end sequencing. In the NAR paper, we requested 50M reads for the DNA samples and 75M for the RNA samples at PE150 on NovaSeq.