Preparation of RNA samples
Timing: 1h
- The quality of the starting RNA samples is critical for the success of the PEAT strategy. We recommend TRIzol reagent to isolate total RNA from higher eukaryotes (both animal and plant). For yeast cells, hot-phenol approach1 tends to give the most reproducible results. To avoid potential biases that may incur during library construction, it is preferable to start with ~150 μg total RNA if the sample is easy to be obtained. For low-abundance samples, it is possible to reduce the amount to 10 μg. In this protocol, we use Drosophila embryos to demonstrate the PEAT strategy. One can adapt the protocol to construct PEAT libraries for other species as well
RNA cleanup
Timing: 1 h
RNeasy Mini kit (QIAGEN) is used to remove genomic DNA contamination in RNA preparation. On-column DNase I digestion step is included according to the manufacturer’s protocol Starting with ~150 μg total RNA, the final product is eluted with 80 μl nuclease-free water.
(Optional Step) The integrity of the purified RNA is checked with1.2% formaldehyde agarose gel (FA gel). Alternatively, Agilent RNA Nano Labchip can be used to determine the RNA quality. Proceed to the next step if the RNA samples pass quality control.
PolyA+ RNA selection with Dynabeads (dT)25
Timing: 3 h
Perform two rounds of polyA+ selection with Dynabeads (dT)25 according to manufacturer’s protocol. Use all of the purified total RNAs from the previous step and elute the polyA+ RNA with 30 μl elution buffer (10mM Tris-HCl, pH7.5).
Lithium salt resulting from polyA+ selection step might affect downstream enzymatic activity. RNeasy MinRlute kit is therefore used to remove residual lithium by following the manufacturer’s protocol. We recommend to elute the column twice, each with 40 μl nuclease-free water. The resulting products are then quantified with a nanodrop spectrophotometer. From 150 μg total RNA, the typical yield of purified polyA+ RNA is 2-3 μg.
Bacterial Alkaline Phosphatase (BAP) treatment
Timing: 2.5 h
- Preparing the following reaction in an RNase-free 1.7ml tube:
RNA from step 5 80 μl
BAP enzyme (Takara) 6 μl
10x BAP buffer (Takara) 10 μl
RNasin ribonuclease inhibitor (Promega) 2.5 μl
Nuclease-free water 1.5 μl
Total reaction volume 100 μl
- Incubate the reaction at 37 °C for 45 min. Add 100 μl nuclease-free water to the reaction. Extract once with 200 μl phenol/chloroform and then with 200 μl chloroform. Perform ethanol precipitation by adding 20 μl 3 M NaOAc (pH 5.2), 2 μl GlycoBlue and 500 μl ethanol. Wash the pellet once with 70% ethanol, followed by adding 85.5 μl nuclease-free water to dissolve RNA.
Critical step: Completely removal of the phosphate group on the uncapped RNA molecules before TAP treatment is critical for reducing background noise.
Tobacco Acid Pyrophosphatase (TAP) treatment
Timing: 2.5 h
- Preparing the following reaction in a RNase-free 1.7 ml tube:
BAP treated RNA from step 7 85.5 μl
10x TAP buffer (Epicentre) 10 μl
TAP enzyme (Epicentre) 2 μl
RNasin ribonuclease inhibitor (Promega) 2.5 μl
Total reaction volume 100 μl
- Incubate the reaction at 37 °C for 1 h. Add 100 μl nuclease-free water to the reaction. Extract once with 200 μl phenol/chloroform and then with 200 μl chloroform. Perform ethanol precipitation by adding 20 μl 3 M NaOAc (pH 5.2), 2 μl GlycoBlue and 500 μl ethanol. Wash the pellet once with 70% ethanol. Dissolve the BAP and TAP treated RNA in 13 μl nuclease-free water.
Ligation with 5’ chimeric linker
Timing: 3 h (or overnight)
- Combing the following in a 1.7ml RNase-free tube:
BAT/TAP treated polyA+ RNA from step 9 13 μl
30 μM PAGE purified Mme_CAP_RDhyb (IDT) 2 μl
10x T4 RNA ligase buffer (NEB) 10 μl
40% PEG8000 (0.2 μM filtered and autoclaved) 62.5 μl
RNasin ribonuclease inhibitor (Promega) 2.5 μl
T4 RNA ligase1 (NEB) 10 μl
Total reaction volume 100 μl
Critical step: The ratio of linker:RNA should be larger than 10:1 to ensure efficient ligation. The sequence of Mme_CAP_RDhyb: 5’-CTC AAG CTT CTA ACG ATG TAC GCT CG rArGrUrCrCrArArC-3’, ‘r’ denotes ribonucleotide.
- Mix the reaction well by pipetting and incubate at room temperature for at least 3 h (or overnight). Add 100 μl nuclease-free water to the reaction and extract once with 200 μl phenol/chloroform and a second time with 200 μl chloroform.
Removal of excess linkers
Timing: 0.5 h
- The supernatant from Step 11 is further purified with RNeasy MinElute kit (QIAGEN) to remove excess linkers by following the manufacturer’s protocol. Elute linker-ligated RNA with 20 μl nuclease-free water.
Reverse transcription with random primer
Timing: 2 h
- Prepare the reaction in a PCR tube:
RNA sample from previous step 20 μl
10mM dNTP, each (Bioline) 2 μl
10uM Mme_RTN6 (IDT) 2 μl
Incubate the reaction for 5 min at 65 °C, then chill on ice for ≥ 1 min. Add the following:
5x first strand buffer (Invitrogen) 8 μl
0.1M DTT (Invitrogen) 2 μl
RNasin ribonuclease inhibitor (Promega) 2 μl
20x Actinomycin D (USB) 2 μl
SuperScript III reverse transcriptase (Invitrogen) 2 μl
Total reaction volume 40 μl
The sequence of Mme_RTN6: 5’-GCG GCT GAA GAC GGC CTA TCC GAC NNN NNN-3’
Mix and incubate the reaction at 25 °C for 5 min, 50 °C for 1 h, 70 °C 15 min. Hold at 10 °C.
Add 1 μl RNase H (2 Units, Invitrogen) to the reaction and incubate at 37 °C for 20 min. This treatment helps to remove RNA strand from RNA:DNA hybrid. Use ZYMO DNA clean & concentrator-5 kit to purify the samples following manufacturer’s protocol. Elute with 20 μl nuclease-free water.
TIP: Add at least 7-fold of the binding buffer to efficiently recover single-stranded cDNA.
Quality control step 1
RT-PCR is used to evaluate the efficiency of linker ligation.
Timing: 2 h
- Prepare the following in a PCR tube:
Purified cDNA from step 15 0.25 μl
5x HF buffer (NEB) 4 μl
10 μM Mme_F (IDT) 1 μl
10 μM Mme_R or fly_GAPDH1_R (IDT) 1ul
10 mM dNTP, each (Bioline) 0.4 μl
Phusion high fidelity DNA polymerase (NEB) 0.2 μl
Nuclease-free water 13.15 μl
Total reaction volume 20 μl
Perform the thermal cycling below: 98 °C for 30 sec; 30 cycles of 98 °C for 10 sec, 67 °C for 10 sec and 72 °C for 90 sec; 72 °C for 10 min.
Sequences of PCR primers: Mme_F: 5’-CTC AAG CTT CTA ACG ATG TAC GCT CGA-3’; Mme_R: 5’-GCG GCT GAA GAC GGC CTA TCC-3’; Fly_GAPDH1_R: 5’-GGG CCG AGA TGA TGA CCT TCT T-3’
Take 3 μl RT-PCR product to run 1.5% agarose gel. QC is passed if positive control (Mme_F + fly_GAPDH1_R) shows specific band, experimental control (Mme_F + Mme_R) shows strong smear and single primer amplification (Mme_F only or Mme_R only) shows no or weak smear.
Amplification of cDNA by low-cycle PCR
Timing: 1 h
- Prepare the following in a UV treated PCR tube:
cDNA from step 19 16 μl
5x HF buffer (NEB) 10 μl
10 μM Mme_F (IDT) 2.5 μl
10 μM Mme_R (IDT) 2.5 μl
10 mM dNTP, each (Bioline) 1 μl
Phusion HF DNA polymerase (NEB) 0.5 μl
Nuclease-free water 17.5 μl
Total reaction volume 50 μl
- Perform the thermal cycling: 98 °C for 30 sec; 5 cycles of 98 °C for 10 sec, 67 °C for 10 se and 72 °C for 90 sec; 72 °C for 10 min. If starting with less amount of RNA, more cycles can be added (up to 15 cycles).
Exo I digestion to remove excess primers
Timing: 1.5 h
- Add 5 μl Exonuclease I (Exo I) to the 50 μl PCR reaction. Incubate at 37 °C for 45 min, 80 °C for 20 min, and hold the reaction at 10 °C. Purify the samples with ZYMO Clean & Concentrator-5 kit. Elute with 10 μl nuclease-free water.
TIP: To avoid airborne DNA contamination, use fresh binding buffer and wash buffer of the ZYMO kit. Or filter the buffer with 0.2 μm filter unit and treated with UV.
Circularization of PCR product
Timing: 3 h
- Combine the following in a UV treated PCR tube:
Purified PCR product from step 21 10 μl
10x Ampligase buffer (Epicentre) 3 μl
OptiKinase (USB) 1.5 μl
100 mM ATP (VWR) 1.5 μl
100 mM DTT (Invitrogen) 0.3 μl
Nuclease-free water 11.3 μl
Reaction volume at this point 27.6 μl
Incubate the reaction at 37 °C for 30 min, then 95 °C for 2 min to heat inactivate OptiKinase. Hold the reaction at 4 °C.
Add the following two components in PCR hood to avoid potential contamination:
10 μM TSS_Col3_short (IDT) 0.9 μl
Ampligase (Epicentre) 1.5 μl
Total reaction volume 30 μl
Mix the reaction and perform thermal cycling as following: 5 cycles (30 sec at 95 °C, 2 min at 68 °C, 1 min at 55 °C, 5 min at 60 °C); 5 cycles (30 sec at 95 °C, 2 min at 65 °C, 1 min at 55 °C, 5 min at 60 °C). Hold the reaction at 10 °C.
Sequence of TSS_Col3_short: 5’-GCC GTC TTC AGC CGC CTCA AGC TTC TAA CGA TGT ACG-3’
Critical step: Any circular DNA contamination in this step will be amplified in the following RCA step and should be avoided2. Use Ampligase buffer, 100mM ATP and 100mM DTT which were aliquoted and kept at -20°C. UV treats tubes and nuclease-free water to remove potential DNA contaminations.
Exonuclease digestion
Timing: 1.5 h
- Add 3 μl Exo I (NEB) and 0.6 μl Exo III (NEB) to the 30 μl circularization reaction. Incubate 45 min at 37 °C, 20 min at 80 °C, and hold at 10 °C. This step eliminates DNA fragments in linear form, and only single-stranded circular DNA remain after the digestion. Directly use the reaction for rolling circle amplification (RCA).
Rolling circle amplification (RCA)
Timing: 17 h
- Combine the following in four UV treated PCR tubes:
Circular DNA from previous step 2 μl
25 mM dNTP, each (Epicentre) 0.8 μl
10 mg/ml BSA (NEB) 0.4 μl
10 x Phi29 buffer (NEB) 2 μl
100 μM N6T2 (IDT) 2 μl
Phi29 high fidelity DNA polymerase (NEB) 1 μl
DMSO, UV-treated (Stratagene) 2 μl
H2O (Ambion, UV-treated) 9.8 μl
Total reaction volume 20 μl
TIP: Before adding Phi29 DNA polymerase, heat the reaction at 80 °C for 2 min and snap cool down on ice for 2 min to reduce bias of random primer binding. N6T2: 5’-NNNNNN-3’; * = phosphothiol.
Tips: Perform multiple reactions to reduce amplification bias. We typically prepare four reactions for each sample.
- Mix the reactions by pipetting and incubate them as following: 10 °C for 10 min, 28 °C for 16 h, 65 °C for 10 min, and hold at 10 °C. Combine the four tubes of same sample after RCA.
Quality control step 2
This QC step is to determine whether the circular DNAs are specifically amplified using a built-in XhoI site in the 5’ adaptor. This step also can help to estimate the distribution of the circular molecules.
Timing: 1.5 h
- QC by XhoI digestion. Prepare the following reaction:
RCA product 2 μl
10 x buffer 2 (NEB) 2 μl
100 x BSA (NEB) 0.2 μl
XhoI (NEB) 1 μl
H2O 14.8 μl
Total reaction volume 20 μl
- Incubate the reaction at 37 °C for 1 h. Take 10 μl reaction to run in a 1.5% agarose gel. Results passed QC: Majority of the RCA products were digested by XhoI and generate smear ranging from 200 bp to 1000 bp and peaked at ~400 bp.
Precipitation of RCA products
Timing: 1.5 h
- For each sample, combine 4 tubes of RCA products together. Add 120 μl nuclease-free water to bring a total volume of 200 μl. Extract once with phenol/chloroform then with chloroform followed by ethanol precipitation. Qualify the RCA products by Nanodrop spectrophotometer. The typical yield is around 60 μg.
Generation of ditag by MmeI digestion
Timing: 2.5 h
- Prepare the following reaction in a 1.7 ml tube:
RCA product from previous step 16 μg
10x buffer 4 (NEB) 18 μl
SAM, dilute to 500 μM using 1x NEB buffer 4 20 μl
MmeI, 2U/μl (NEB) 10 μl
Nuclease-free water add to 200 μl
Total reaction volume 200 μl
TIP: The original concentration of SAM is 32 mM. Use 1x buffer 4 (NEB) to dilute SAM.
- Incubate the reaction at 37 °C for 30 min. MmeI digestion enables the release of small fragment (94bp) which containing ditag3. Concentrate the reaction to 12 μl with two DNA clean & concentrator-5 columns (ZYMO). Load 6 μl of purified DNA to one lane (one sample has two lanes) in 6% acrylamide gel. We use 18 x 16 cm acylamide gel to obtain better separation of MmeI digested product (94 bp) from background. Run the vertical gel at 200 V for 1.5 h. Compare the intensity of bands with the DNA standard (100 bp band of Invitrogen’s 10 bp DNA ladder), we estimate that the desired DNA fragment is ~7.2 pmol.
Purification of MmeI digested small fragment
Timing: overnight
- Use a new Razor Blade (VWR) to cut the desired band (~94 bp) which contains both TSS sequence and 3’ downstream tag. Elute each gel slice with 400 μl gel elution buffer (0.1% SDS, 0.2 NaCl, 10 mM MgOAc2). Rotate at least 6 h at room temperature or overnight at 4 °C. Phenol/chloroform extraction followed by ethanol precipitation. Resuspend the sample in 20 μl nuclease-free water.
Ligation with paired-end linkers
Timing: overnight
- Combine the following in a PCR tube:
Purified DNA from previous step 20 μl
12.5 μM PE2_AN2 1.5 μl
12.5 μM PE2_BN2 1.5 μl
T4 DNA ligase buffer (NEB) 5 μl
40% PEG8000 (filtered and autoclaved) 20 μl
T4 DNA ligase, 2000U/ul (NEB) 2 μl
Total reaction volume 50 μl
Incubate the reaction overnight at 16 °C.
TIP: Mix the reaction for at least 20 times to ensure the even distribution of PEG8000.
Linker sequences:
PE2_AN2:
5’-ACCGAGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCTNN-3’
3’-TGGCTCTAGATGTGAGAAAGGGATGTGCTGCGAGAAGGCTAGA-p-5’
PE2_BN2:
5’-p-AGATCGGAAGAGCGGTTCAGCAGGAATGCCGAGACCGATCT-3’
3’-NNTCTAGCCTTCTCGCCAAGTCGTCCTTACGGCTCTGGCTAGA-5’
‘p’ denotes phosphate group; ‘N’ represents any nucleotide.
- Phenol/chloroform extraction followed by concentrate with ZYMO Clean & Concentrator-5 kit. Elute the ligation product with 10 μl nuclease-free water.
Final library amplification by PCR
Timing: 1.5 h
- Prepare the following in a PCR tube:
Linker ligated DNA from previous step 10 μl
5x HF buffer (NEB) 10 μl
1 μM PE2_AFN2_short (IDT) 2.5 μl
10 μM PE2_A_short (IDT) 2.5 μl
1 μM PE2_BRN2_short (IDT) 2.5 μl
10 μM PE2_B_short (IDT) 2.5 μl
10 mM dNTP, each (Bioline) 1 μl
Phusion HF DNA polymerase (NEB) 0.5 μl
Nuclease-free water 18.5 μl
Total reaction volume 50 μl
Sequence of primers:
PE2_AFN2_short: 5’-AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACA-3’
PE2_BRN2_short: 5’-CAAGCAGAAGACGGCATACGAGATCGGTCTCGGCATTCCT-3’
PE2_A_short: 5’-AATGATACGGCGACCACCGAGA-3’
PE2_B_short: 5’-CAAGCAGAAGACGGCATACGAGA-3’
- Run the following PCR program: incubate for 30 sec at 98 °C, followed by a 2 cycles of (10 sec and 98 °C, 10 sec at 66 °C, 30 sec at 72°C) and 10 cycles of (10 sec and 98 °C, 10 sec at 69 °C, 30 sec at 72°C), additional extension for 10 min at 72 °C and then hold at 10 °C.
Quality control step 3
Timing: 0.5 h
- Take 3 μl of PCR product to run a 1.5% agarose gel. One should see a specific band of 213 bp. If very strong non-specific amplification occurs, perform size selection (~180 bp) after linker ligation.
TIP: A library with good quality should have a relatively strong band (~213 bp) compared to the primer dimers.
Size selection in an acrylamide gel
Timing: overnight
- Use ZYMO DNA clean & concentrator-5 kit to purify the remaining 47 μl PCR product and elute with 6 μl nuclease-free water. Load the purified sample to an 8% acrylamide gel and run at 1 w for 1 hr. Stain the acrylamide gel with EtBr for 2 min and cut the ~213 bp band under UV transilluminator. Elute the gel slice with 400 μl gel elution buffer (0.1% SDS, 0.3M NaOAc, pH5.2) and rotate at room temperature for at least 6 h or at 4 °C overnight. Purify the final library with ZYMO DNA clean & concentrator-5 kit. Load multiple times for mixture more than 700 μl. Elute with 30 μl nuclease-free water. Quantify the final library by Qubit fluoromerter according to manufacturer’s protocol. The concentration of our sample (fly 0-24h embryo) is 1.3 ng/μl.
TIP: This size selection step is to remove non-specific amplification. It helps to increase the quality of paired-end library.
Quality control step 4
We use Sanger sequencing to evaluate the quality of the final TSS library.
Timing: two days
- Use 5 μl out of 30 μl final library for the QC by Sanger sequencing. We obtained 20 qualified clones; 18 of the clones could be mapped to annotated Drosophila TSSs (± 250nt).
Paired-end sequencing by Illumina GAIIx
Timing: 3~4 days
- The final library passed QC was subjected to paired-end sequencing by Illumina Genome Analyzer IIx.