3DRAM-seq enables joint epigenome profiling of spatial genome organization, chromatin accessibility and DNA methylation at high resolution
Florian Noack, Boyan Bonev
Abstract
3DRAM-seq ( 3D genome, R NA, A ccessibility and M ethylation sequencing) is a novel multiomic method that simultaneously interrogates spatial genome organization, chromatin accessibility, DNA methylation and gene expression at high resolution. It utilized enzymatic treatment of bulk fixed nuclei with the GpC methyltransferase M.CviPI prior to restriction digestion, biotin fill-in and proximity ligation as in Hi-C. The resulting 3C DNA is then sheared to ~550bp, followed by bisulfite conversion and library preparation based on the Accel-NGS® Methyl-Seq DNA but using streptavidin beads after the adaptor ligation step to enrich for fragments containing informative contacts. 3DRAM-seq can be performed with as few as 200 000 cells and can be combined with immunoFACS to profile the epigenome landscape in a cell-type-specific manner.
Steps
3DRAM-seq: Methylation controls
NOTE: Control DNA has to prepared only once and can be reused.
To prepare GpC methylated control DNA, mix 10µl of CpG methylated pUC19 DNA (Zymo Research, Cat. N.: D5017) with 10µl of unmethylated lambda DNA (Promega, Cat. N: D1521).
Perform GpC methylation by mixing 20µl of the pUC19/lambda DNA mix with 30µl of the GpC methylation mix I and incubate for 1h @ 37 ˚C followed by 20 min. at 65 degree.
Add 40ul (0.8x) of room temperature Ampure XP beads and mix well by pipetting up and down 10 times. Incubate at RT for 5min.
Separate on a magnet and remove the clear solution.
Keeping the beads on the magnet, wash twice with 200μl of 80% ethanol without mixing.
Leave the beads on the magnet for 5 minutes to allow the remaining ethanol to evaporate.
Resuspend the beads with 16µl nuclease-free water, incubate for 10 minutes at RT, pellet the beads on a magnetic-rack and transfer the DNA containing supernatant into a Covaris microTUBE-15 AFA Beads Screw-Cap tubes (Covaris, Cat. N: 520145).
Shear the DNA to a targeted size of 550bp using following parameter on a Covaris S220 sonicator:
- Peak Incident Power (W) 18
- Duty Factor 20%
- Cycles per Burst 50
- Treatment Time (s) 22
Perform biotin fill-in reaction by mixing:
15µl sheared DNA
75µl steril water
10µl 10x DpnII buffer (New England Biolabs, Cat. N: R0543)
20µl Biotin fill-in master mix (see Materials)
Incubate at 23°C for 4 hours (900 RPM mixing; 10 sec. every 5 min.).
Perform a 1x AmpureXP (Agencourt, Cat. N: A63881) purification by adding 120 µl of room temperature AmpureXP beads to the samples followed by a incubation of 10 minutes at RT under slow rotation.
Pellet the beads on a magnatic-rack, remove the supernatant and wash the beads twice with freshly prepared 80% Ethanol.
After the last wash remove as much as possible of the 80% Ethanol and air-dry the beads for approximately 5 minutes at RT.
Resuspend the beads with 50µl nuclease-free water, incubate for 10 minutes at RT, pellet the beads on a magnetic-rack and transfer the DNA containing supernatant into a new tube.
Quantify the yield using the Qubit dsDNA High Sensitivity Assay (ThermoFisher, Cat. N. Q32851). Control DNA can be stored at -20˚C.
Cell fixation
Resuspend dissociated cells in PBS to reach a maximal cell concentration of 2x106 cells/ml and add freshly prepared 2% Formaldehyde solution to reach a final concentration of 1%. Incubate for 10 min. at RT with slow rotation.
Add 2.0M glycine solution (Invitrogen, Cat. N: 15527013) to a final concentration of 0.2M to quench the reaction. Incubate at room temperature for 5 minutes at rotating wheel.
Centrifuge for 5 minutes at 500xG at 4°C. Discard supernatant into an appropriate collection container.
Resuspend cells in 1ml of cold 1x PBS with 0.5% BSA (ThermoFisher, Cat. N: AM2618) and count the cell number.
Fixed cells are pelleted (2500xg for 5min. at 4˚C) and can be either used directly for 3DRAM-seq, ImmunoFAC-sorted or snap frozen in liquid nitrogen for storage at -80˚C.
Please note that snap freezing leads to degradation of RNA, therefore cells for RNA isolation should be directly processed after fixation or ImmunoFACS.
RNA-seq library construction
NOTE: this protocol was optimized for ~2x104 cells. 4 cells.
Pelleted cells are resuspended with 95µl nuclease free water and RNA is isolated using the Quick-RNA FFPE Miniprep (Zymo Research, Cat. N: R1008) according the manual instructions starting from the tissue digestion step.
To minimize the elution volumne Zymo-Spin IC Columns (Zymo Research, Cat. N: C1004-250) can be used.
Purified RNA is quantified using Qubit RNA HS Assay Kit (ThermoFisher, Cat. N.: Q32852) and quality of the RNA is measured using the Bioanalyzer High Sensitivity RNA 6000 Pico Kit (Agilent, Cat. N.: 5067-1513).
100ng of high quality RNA was used for RNA library generation using the NEBNext® Single Cell/Low Input RNA Library Prep Kit (New England Biolabs, Cat. N.: E6420) according to the manual instruction.
3DRAM-seq: Cell lysis and GpC methylation
NOTE: this protocol was optimized for ~2x105 cells. 5 cells.
Pellet fixed cells for 5min. at 2500xg (4˚C) and carefully resuspend them with 190µl of ice-cold Lysis Buffer.
If cells were frozen, thaw the cell pellet first on ice.
Incubate on ice for exactly 10 minutes and centrifuge at 2500xg for 5 minutes at 4°C. and remove the supernatant.
Wash the nuclei pellet by adding 250µl GpC wash buffer followed by a centrifugation at 2500xg for 5 minutes (4°C) and the removal of the supernatant.
Carefully resuspend nuclei with 100µl GpC reaction buffer.
Incubate for 3 hours at 37 ˚C with slight shaking of the tubes (500 RPM) and subsidize the reaction every hour with 1µl 32mM SAM (New England Biolabs, Cat. N: M0227S) and 2µl of M.CviPI (New England Biolabs, Cat. N: M0227S).
3DRAM-seq: DpnII Digestion and biotin fill-in
Wash nuclei once with 190µl MethylHiC wash buffer. Centrifuge at 2500xg for 5 minutes and remove as much of the supernatant as possible without disturbing the nuclei pellet.
Gently resuspend nuclei pellet in 40μl of 0.5% SDS (Invitrogen, Cat. N: AM9823) and incubate at 62°C for 10 minutes. Afterwards place on room temperature.
Add first 112μl of water and then 20μl of freshly made 10% Triton X-100 solution (Sigma Aldrich, Cat. N: X100) to quench the SDS. Mix well by carefully pipetting up-down (avoiding excessive foaming). Incubate at 37°C for 15 minutes on a thermomixer with a rotationspeed of 600 rpm.
Add 20μl of 10X DpnII buffer (New England Biolabs, Cat. N: R0543) and 8µl (400U) of DpnII restriction enzyme (New England Biolabs, Cat. N: R0543) and digest chromatin overnight at 37°C with a rotationspeed of 600 rpm.
Next day incubate samples at 62°C for 20 minutes, then cool the samples to room temperature.
To fill in the restriction overhangs and mark the DNA ends with biotin, add 40μl of Biotin fill-in master mix . Mix by pipetting carefully and incubate for 4 hours at 23°C (900 rpm mixing; 10 sec every 5 min).
MethylHiC: Proximity ligation and decrosslinking
Add 720μl of the Ligation master mix. Mix by inverting and incubate at least for 6h at 16°C (900 rpm mixing, 10 sec every 30 min). Alternativly the ligation reaction can be also performed overnight.
Degrade proteins by adding 10μl of 20mg/ml Proteinase K (New England Biolabs, Cat. N: P8107) and 96μl of 10% SDS (Invitrogen, Cat. N: AM9823) and incubate at 55°C for 30 minutes.
Add 106μl of 5M NaCl (ThermoFisher, Cat. N: AM9760G) and incubate at 68°C overnight (900 rpm 10sec, every 5 min.).
Cool tubes at room temperature, add 1:1 volume Phenol/Chloroform/Isoamylalcohol (Invitrogen, Cat. N. 15593031), vortex for at least 30sec. and spin with maximal speed for 10minutes at 4˚C.
Distribute the upper aqueous phase into two 2ml tubes and add to each tube 2x volumes of ice cold 100% ethanol, 1/10th of the volume 3M NaAc (Ambion, Cat. N: AM9740) and 1µl Glycogen (ThermoFisher, Cat. N.: 10814010). Store the mix at -20˚C at least for 1 hour or -80˚C for 10-15 minutes.
Centrifuge at maximal speed for 20 minutes at 4°C. Carefully remove the supernatant by pipetting.
Wash with 800μl freshly prepared 70% cold ethanol followed by centrifuge at maximal speed for 5 minutes at 4°C.
Wash each tube with 200µl freshly prepared 70% cold ethanol and pool the samples together in one 0.5ml tube (total volume ~ 400µl, make sure the glycogen pellets are transfered). Centrifuge at maximal speed for 5 minutes at 4°C.
Remove as much as possible of the ethanol and airdry the pellet. Dissolve pellet in 16ul of 10mM Tris buffer pH 7.5 (ThermoFisher, Cat. N: 15567027) and incubate at 37°C for 15 minutes to fully dissolve DNA.
3DRAM-seq: DNA shearing (Covaris)
Transfer 16µl of the sample into a Covaris microTUBE-15 AFA Beads Screw-Cap tube (Covaris, Cat. N: 520145) and sheare DNA to a targeted size of 550bp using following parameter on a Covaris S220 sonicator:
- Peak Incident Power (W) 18
- Duty Factor 20%
- Cycles per Burst 50
- Treatment Time (s) 22
Quantify the yield using the Qubit dsDNA High Sensitivity Assay.
3DRAM-seq: Biotin removal
Mix 15µl of the sheared sample with 5µl of the Endrepair master mix.
Incubate for 4h at 20˚C followed by a heat inactivation for 20 minutes at 75˚C. Sample can be stored at -20˚C.
3DRAM-seq: Bisulfite Conversion
NOTE: for the highest yield the libraries should be prepared directly after the bisulfite conversion.
Prior the bisulfite conversion add roughly 0.01% GpC methylated-biotinylated control DNA to the sample. Proceed with the bisulfite conversion using the EZ DNA Methylation-Gold Kit (Zymo Research) according to the manual instructions until the elution step.
Elute the bisulfite converted DNA in 31µl low EDTA TE water provided with the Accel-NGS® Methyl-Seq DNA Library Kit (Swift Bioscience, Cat. N: 30024; now xGen Methyl-Seq DNA Library Prep IDT, Cat. N.: 10009860).
3DRAM-seq: Library construction and Biotin pulldown
Per sample prepare 25µl Dynabeads MyOne Streptavidin T1 beads (ThermoFisher, Cat. N: 65602) by washing them with 400µl 1x Tween washing buffer and resuspend the beads in 60µl of 2x Binding buffer.
Use 30µl of the bisulfite converted DNA for the library preparation using the Accel-NGS® Methyl-Seq DNA Library Kit with the double reaction volume according to the manual instruction until the ligation step (step 16: 15 minutes at 25˚C).
Please note that one reaction volume can be used if less than 100ng of DNA was used for the bisulfite conversion.
After the ligation resuspend samples with 30 µl steril water and subsequently add 60µl of the washed MyOne Streptavidin T1 beads in 2x Binding buffer. Incubate for 15 minutes at room temperature under slow rotation.
Wash beads 4 times with 400µl 1x Tween washing buffer for 2 minutes at 55˚C.
After the last wash remove as much as possible from the Tween washing buffer and wash beads twice with 400 µl sterile water. Thereafter resuspend beads in 50µl nuclease free water.
Final amplification is performed with the with the EpiMark Hot Start Taq (New England Biolabs, Cat. N: M0490S) in 5 separate reactions in order to archive a higher complexity. To do so pipette 10µl of streptavidin bead solution in 5 separate PCR tubes and add 15µl of the Library amplification mix containing sample specific Methyl-Seq Set A Indexing Primers (Swift Bioscience, Cat. N: 36024; now IDT, Cat. N. 10009965 or 10005975).
Amplify library using the following PCR program: 95˚C 30s; {95˚C 15s, 61˚C 30s, 68 ˚C 80s} x10-11; 68˚C 5min; Hold at 10˚C.
After the PCR is finished, pellet the beads and pool the library containing supernatant of different reactions of a sample together (~125µl total volume).
Perform a 0.65x AmpureXP (Agencourt, Cat. N: A63881) purification by adding 81.25µl of room temperature AmpureXP beads to the samples followed by a incubation of 10 minutes at RT under slow rotation.
Please note that a <0.6x AmpureXP beads can lead to the loss of the whole library !
Pellet the beads on a magnatic-rack, remove the supernatant and wash the beads twice with freshly prepared 80% Ethanol.
After the last wash remove as much as possible of the 80% Ethanol and air-dry the beads for approximately 5 minutes at RT.
Resuspend the beads with 15µl nuclease-free water, incubate for 10 minutes at RT, pellet the beads on a magnetic-rack and transfer the DNA containing supernatant into a new tube. Quantify the yield using the Qubit dsDNA High Sensitivity Assay (ThermoFisher, Cat. N. Q32851)