- Short Report
- Open Access
Complexities due to single-stranded RNA during antibody detection of genomic rna:dna hybrids
BMC Research Notesvolume 8, Article number: 127 (2015)
Long genomic R-loops in eukaryotes were first described at the immunoglobulin heavy chain locus switch regions using bisulfite sequencing and functional studies. A mouse monoclonal antibody called S9.6 has been used for immunoprecipitation (IP) to identify R-loops, based on the assumption that it is specific for RNA:DNA over other nucleic acid duplexes. However, recent work has demonstrated that a variable domain of S9.6 binds AU-rich RNA:RNA duplexes with a KD that is only 5.6-fold weaker than for RNA:DNA duplexes. Most IP protocols do not pre-clear the genomic nucleic acid with RNase A to remove free RNA. Fold back of ssRNA can readily generate RNA:RNA duplexes that may bind the S9.6 antibody, and adventitious binding of RNA may also create short RNA:DNA regions. Here we investigate whether RNase A is needed to obtain reliable IP with S9.6.
As our test locus, we chose the most well-documented site for kilobase-long mammalian genomic R-loops, the immunoglobulin heavy chain locus (IgH) class switch regions. The R-loops at this locus can be induced by using cytokines to stimulate transcription from germline transcript promoters. We tested IP using S9.6 with and without various RNase treatments. The RNase treatments included RNase H to destroy the RNA in an RNA:DNA duplex and RNase A to destroy single-stranded (ss) RNA to prevent it from binding S9.6 directly (as duplex RNA) and to prevent the ssRNA from annealing to the genome, resulting in adventitious RNA:DNA hybrids. We find that optimal detection of RNA:DNA duplexes requires removal of ssRNA using RNase A. Without RNase A treatment, known regions of R-loop formation containing RNA:DNA duplexes can not be reliably detected. With RNase A treatment, a signal can be detected over background, but only within a limited 2 or 3-fold range, even with a stable kilobase-long genomic R-loop.
Any use of the S9.6 antibody must be preceded by RNase A treatment to remove free ssRNA that may compete for the S9.6 binding by forming RNA:RNA regions or short, transient RNA:DNA duplexes. Caution should be used when interpreting S9.6 data, and confirmation by independent structural and functional methods is essential.
Physical presence of long genomic R-loops in eukaryotes was first described at the IgH switch regions . This genomic R-loop analysis was based on (a) bisulfite sequencing [1,2] and (b) oligonucleotide hybridization (colony lift hybridization) [3,4], and both were done with and without an RNase H challenge treatment to confirm the RNA:DNA conformation. Moreover, we did not observe any large effects of parallel RNase A treatment in vivo [1,3-5] or in vitro [6-8]. Our original description of kilobase long mammalian genomic R-loops was further built upon and had the advantage of several lines of independent evidence including (a) the large body of IgH switch DNA sequence and recombination junctional sequence information [9,10]; (b) many functional studies of IgH switch region transcription ; (c) concurrent studies of IgH switch region orientation [12,13]; and (d) detailed in vitro biochemical studies of transcription through switch regions [6-8,14-17]. Therefore, these well-documented regions with R-loops are ideal positive control targets of IP using S9.6 antibody.
While S9.6 antibody can recognize RNA:DNA duplexes [18-28], complete characterization of the binding specificity of S9.6 was initially limited to ELISA measurements on its binding to long nucleic acid duplexes . Such ELISA measurements can be complicated by multiple antibodies binding to a single long duplex. This multi-antibody complex would reflect the combination of affinities of multiple antibodies [(KD)n, where n = the number of antibodies bound to a given duplex]. Recent work has shown that a single-chain variable domain of the S9.6 antibody can bind RNA:RNA duplexes with an affinity that is only 5.6-fold weaker than to RNA:DNA duplexes, raising the serious concern that S9.6 can indeed cross react with RNA species . Most of the studies that have used S9.6 to identify R-loops in eukaryotes have not used RNase A to eliminate any artifacts due to free RNA [18-28,30,31], which might be present during cell lysis and/or harvest of the nucleic acid. In addition, free RNA can reanneal with the template DNA during transient breathing of the DNA, and short RNA:DNA hybrids also exist during DNA replication at RNA primer annealing sites. Here we examine the complexities of R-loop analysis when using the S9.6 antibody with and without various RNase treatments. We used the mouse B-cell line, CH12F3.2a , which is able to specifically and efficiently switch to IgA upon cytokine stimulation, thus providing the only widely-accepted extended R-loop as a positive control [33-35].
CH12F3.2a and its derivative cells were cultured in RPMI medium supplemented with 10 % FCS and 50 μM β-mercaptoethanol . For experiments specifying cytokine stimulation, two million healthy CH12F3.2a cells at a density around 1×106 cells/ml were supplemented with anti-CD40 (eBioscience Cat. No. 16-0404-86), IL-4 (R&D Cat. No. 404-ML-010) and TGF- β1 (R&D Cat. No. 240-B-002) for 24 hours.
ATCC HB-8730 hybridoma line (generously provided by Bradley Cairns) was cultured in a CELLine 1000 bioreactor (Satorius Biotech, NY) according to manufacturer’s instructions. Harvested antibody (culture supernatant) was purified on a column packed with Protein G Sepharose 4 Fast Flow (GE Healthcare) equilibrated with 1x phosphate buffered saline (PBS).
S9.6 Immunoprecipitation and qPCR
Genomic DNA from CH12F3.2a cells with and without cytokine stimulation was prepared by overnight proteinase K digestion, phenol-chloroform extraction and ethanol precipitation. Genomic DNA was digested with EcoRI; importantly, RNase A was added at this step to prevent S9.6 antibody binding to RNA species in subsequent steps in the experiments with RNase A treatment (Phillips et al. ). Five microgram of digested genomic DNA were incubated with 5 μg of S9.6 antibody in 400 μl IP buffer (10 mM sodium phosphate pH 7.0, 140 mM NaCl, 0.1% Tween 20) for 2 h at 4°C with rotation. Ten microliter pre-blocked Dynabeads (Invitrogen 10004D) were added to the mixture and gently rotated at 4°C for 2 h. The beads were washed with IP buffer three times before the bead-bound DNA was recovered by proteinase K treatment overnight, phenol-chloroform extraction, and ethanol precipitation. The recovered DNA was then analyzed by quantitative real-time PCR (qPCR) using TaqMan probes. The sequence details of the primers and probes are shown in Additional file 1: Table S1. For the beta-actin gene, the probe for real-time PCR was: ZZ218 FAM-CACCGCAAGTGCTTCTAGGCGGAC-BHQ, and the primers were: ZZ221 - GAGTCCGGCCCCTCCAT and ZZ222 - GGTTTTGTCAAAGAAAGGGTGTAAA. Each IP experiment was done in triplicate, and each qPCR is done in duplicate.
R-loops at the Mammalian IgH Switch Regions
We stimulated CH12F3.2a cells (which are capable of robust cytokine-dependent CSR in culture) for 24 hr before genomic DNA was harvested and digested in preparation for immunoprecipitation (IP) as described previously . The transcription level of the IgH Sα in CH12F3.2a cells increased slightly upon cytokine stimulation . Four different DNA regions were chosen for qPCR analysis after IP. The first two are the IgH Sα region and IgH Sμ region, which we have shown to form R-loops within cells, based on genomic bisulfite sequencing [3,33,38]. IgH Sμ is already highly transcribed without cytokine stimulation; therefore, cytokine stimulation does not increase R-loop formation in this region of the IgH locus. We chose two nearby regions as controls. The first control region is located upstream of the adjacent transcriptional promoters. The second control region is located far downstream (Figure 1; note IgH locus with EcoRI sites and qPCR amplicons). These controls reside within different restriction fragments that are distinct from the Ig class switch regions.
In these IP experiments , we find that the pull down of the DNA fragment containing the Sα region is more than two fold higher with cytokine stimulation (Figure 2, Group 5 versus Group 6, column b) while other regions showed no or very little change (Figure 2, Group 5 versus Group 6, columns a, c, and d). The downstream region, which showed a very small increase of pull down in the stimulated cells, may be due to the R-loop formation in the IgH Sα affecting the efficient digestion of EcoRI site in the DNA fragment. To further confirm that the signals that we detected arose from RNA:DNA hybrids, we treated the DNA samples with RNase H before the IP experiment, an endonuclease that specifically degrades RNA in RNA:DNA duplexes. Upon RNase H treatment, the signals from all regions that we examined decreased 5 to 10-fold (Figure 2, compare Group 5 versus Group 4) indicating that the IP signals are indeed due to RNA:DNA hybrids.
It is surprising that we detected signals clearly above the background within the fragments far upstream and downstream of the alpha locus even though these signals do not change with cytokine stimulation. We have previously documented that the region upstream of the germline transcript promoter for the IgH alpha locus and the downstream region do not harbor any R-loops, regardless of cytokine treatment [1,3,4,33]. We confirmed that this was not merely due to incomplete restriction digestion because the signal is present after multiple rounds of digestion.
We also examined the beta-actin gene in the DNA samples from the same IP experiments. The beta-actin gene is at another location in the genome, and it has not been documented to have any extended R-loop by physical methods. The signal at the beta-actin gene is similar to all regions other than IgH Sα (Figure 2, all Groups, column e), and this signal also does not change with cytokine stimulation. These findings indicate that while the S9.6 antibody binds to well-documented R-loops in the IgH Sα, this antibody may also bind to DNA regions that are clearly known to have no extended R-loops. In addition to binding to bona fide extended R-loops [1,3,4,33], it is possible that the S9.6 antibody can bind to very small RNA:DNA hybrids, such as those formed by RNA primers for DNA replication and the RNA transcripts that may transiently reanneal back to the DNA template during or before genomic DNA is harvested. This may account for the RNase H-sensitivity of these IP signals at locations outside of IgH switch regions (Figure 2, Group 4).
Taken together, the S9.6 antibody appears to have good specificity for binding to RNA:DNA hybrids, but its ability to distinguish extended R-loops versus small RNA:DNA hybrids arising during DNA replication (RNA primer sites) or from other causes is highly uncertain. With the potential high background due to the antibody binding to non-R-loop RNA:DNA hybrids, the dynamic range of the R-loop signal is greatly affected.
Cellular RNA complicates R-loop detection
In all the IP experiments of DNA samples without RNase A treatment regardless of cytokine stimulation, qPCR did not detect any signal above background in any regions examined (Figure 2, Groups 1 to 3). This is also surprising since RNase A does not affect RNA:DNA hybrids. This finding indicates that the vast amount of RNA in the DNA samples can indeed affect S9.6 binding to RNA:DNA hybrids. While the mechanism involved in this interference is undetermined in our experiments, the increase of signals in experiments with RNase A treatment reflects the possibility that RNA species in the DNA may titrate out the S9.6 antibody (Figure 2, Groups 2 and 3 versus Groups 5 and 6). It is possible that the free RNA can fold back on itself to form RNA:RNA duplexes, which is known to compete with RNA:DNA for S9.6 binding . This finding strongly indicates that RNase A pre-treatment is necessary for the S9.6 antibody to clearly detect the presence of RNA:DNA duplexes.
The IgH class switch region R-loops can be up to kilobases in length, and thus provide the largest possible target upon which to test the S9.6 antibody for detection of RNA:DNA duplex in the mammalian genome [1,2]. Despite such a large target, we found that the detection of these kilobase R-loops is less than 3-fold above background with borderline significance even with the essential pre-treatment of genomic nucleic acid with RNase A. We also found that the regions known to have no extended R-loops can also be pulled down by the S9.6 antibody. We suspect that abundant free RNA transcripts can reanneal to the template DNA strands to form RNA:DNA hybrids for short lengths at low levels throughout the genome. Also, RNA primers generated during the course of physiologic cellular DNA replication in the sub-population of proliferating cells in S phase might also have short (~20 bp) RNA:DNA duplexes that become sites of S9.6 binding. It is likely that these two forms of RNA:DNA hybrids compete with the long and stable R-loops, such as at the IgH switch locus, and result in promiscuous S9.6 binding on a genome-wide basis in IP experiments (Additional file 1: Figure S1) These findings suggest that while the S9.6 antibody binds to RNA:DNA hybrids with good specificity, its ability to distinguish extensive R-loops of biological significance and very short RNA:DNA hybrids due to transient RNA reannealing to DNA is uncertain. This not surprising since the S9.6 antibody was generated using short RNA:DNA hybrids. We also found that RNA remaining in the genomic DNA harvested from the cells affects S9.6 binding to RNA:DNA hybrids dramatically. It is quite common that free RNA folds back on itself to form RNA:RNA duplexes, which is known to compete with RNA:DNA for S9.6 binding . This binding competition would potentially contribute to the low level of the IP signal when RNase A is not used and thus RNA:RNA duplexes are abundant. Removal of cellular RNA with RNase A reduces this competition of the antibody, thereby allowing S9.6 to bind more specifically to RNA:DNA hybrids. Based on our experiments, it is clear that identification of biologically important R-loops based solely on S9.6 is insufficient, especially when RNase A is not used prior to the assay. Findings in such studies should always be accompanied by chemical probing [1,2], quantitation of the number of alleles involved in the R-loop formation [3,4], and functional studies showing biological effects [12,34,35]. Without rigorous confirmation using independent methods, conclusions drawn from such studies would be highly questionable.
Availability of supporting data
We will immediately provide anyone with our primary data.
Class switch recombination
Quantitative PCR (also termed real-time PCR)
Yu K, Chedin F, Hsieh C-L, Wilson TE, Lieber MR. R-loops at immunoglobulin class switch regions in the chromosomes of stimulated B cells. Nature Immunol. 2003;4:442–51.
Yu K, Roy D, Huang FT, Lieber MR. Detection and structural analysis of R-loops. Methods Enzymol. 2006;409:316–29.
Huang F-T, Yu K, Balter BB, Selsing E, Oruc Z, Khamlichi AA, et al. Sequence-dependence of chromosomal R-loops at the immunoglobulin heavy chain Smu class switch region. Mol Cell Biol. 2007;27:5921–32.
Huang F-T, Yu K, Hsieh C-L, Lieber MR. The downstream boundary of chromosomal R-loops at murine switch regions: implications for the mechanism of class switch recombination. Proc Natl Acad Sci. 2006;103:5030–5.
Yu K, Lieber MR. Nucleic acid structures and enzymes in the immunoglobulin class switch recombination mechanism. DNA Repair. 2003;2:1163–74.
Roy D, Lieber MR. G clustering is important for the initiation of transcription-induced R-loops in vitro, whereas high G density without clustering is sufficient thereafter. Mol Cell Biol. 2009;29(11):3124–33.
Roy D, Yu K, Lieber MR. Mechanism of R-loop formation at immunoglobulin class switch sequences. Mol Cell Biol. 2008;28(1):50–60.
Roy D, Zhang Z, Lu Z, Hsieh CL, Lieber MR. Competition Between the RNA Transcript and the Nontemplate DNA Strand During R-Loop Formation In Vitro: A Nick Can Serve as a Strong R-loop Initiation Site. Mol Cell Biol. 2010;30:146–59.
Dunnick WA, Hertz GZ, Scappino L, Gritzmacher C. DNA sequence at immunoglobulin switch region recombination sites. Nucl Acid Res. 1993;21(3):365–72.
Davis M, Kim SK, Hood LE. DNA Sequences Mediating Class Switching in α-Immunoglobulins. Science. 1980;209(September):1360–5.
Chaudhuri J, Basu U, Zarrin A, Yan C, Franco S, Perlot T, et al. Evolution of the immunoglobulin heavy chain class switch recombination mechanism. Adv Immunol. 2007;94:157–214.
Shinkura R, Tian M, Khuong C, Chua K, Pinaud E, Alt FW. The influence of transcriptional orientation on endogenous switch region function. Nature Immunol. 2003;4:435–41.
Daniels GA, Lieber MR. Strand-specificity in the transcriptional targeting of recombination at immunoglobulin class switch sequences. Proc Natl Acad Sci U S A. 1995;92:5625–9.
Reaban ME, Griffin JA. Induction of RNA-stabilized DNA conformers by transcription of an immunoglobulin switch region. Nature. 1990;348(22 November):342–4.
Reaban ME, Lebowitz J, Griffin JA. Transcription induces the formation of a stable RNA.DNA hybrid in the immunoglobulin alpha switch region. J Biol Chem. 1994;269:21850–7.
Daniels GA, Lieber MR. RNA:DNA complex formation upon transcription of immunoglobulin switch regions: implications for the mechanism and regulation of class switch recombination. Nucl Acids Res. 1995;23:5006–11.
Yu K, Roy D, Bayramyan M, Haworth IS, Lieber MR. Fine-structure analysis of activation-induced deaminase accessbility to class switch region R-loops. Mol Cell Biol. 2005;25:1730–6.
Stirling PC, Chan YA, Minaker SW, Aristizabal MJ, Barrett I, Sipahimalani P, et al. R-loop-mediated genome instability in mRNA cleavage and polyadenylation mutants. Genes Dev. 2012;26(2):163–75.
Skourti-Stathaki K, Proudfoot NJ, Gromak N. Human senataxin resolves RNA/DNA hybrids formed at transcriptional pause sites to promote Xrn2-dependent termination. Mol Cell. 2011;42(6):794–805.
Wahba L, Amon JD, Koshland D, Vuica-Ross M. RNase H and multiple RNA biogenesis factors cooperate to prevent RNA:DNA hybrids from generating genome instability. Mol Cell. 2011;44(6):978–88.
Sun Q, Csorba T, Skourti-Stathaki K, Proudfoot NJ, Dean C. R-loop stabilization represses antisense transcription at the Arabidopsis FLC locus. Science. 2013;340(6132):619–21.
Mischo HE, Gomez-Gonzalez B, Grzechnik P, Rondon AG, Wei W, Steinmetz L, et al. Yeast Sen1 helicase protects the genome from transcription-associated instability. Mol Cell. 2011;41(1):21–32.
Yang Y, McBride KM, Hensley S, Lu Y, Chedin F, Bedford MT. Arginine Methylation Facilitates the Recruitment of TOP3B to Chromatin to Prevent R Loop Accumulation. Mol Cell. 2014;53(3):484–97.
Pfeiffer V, Crittin J, Grolimund L, Lingner J. The THO complex component Thp2 counteracts telomeric R-loops and telomere shortening. EMBO J. 2013;32(21):2861–71.
Balk B, Maicher A, Dees M, Klermund J, Luke-Glaser S, Bender K, et al. Telomeric RNA-DNA hybrids affect telomere-length dynamics and senescence. Nat Struct Mol Biol. 2013;20(10):1199–205.
Nakama M, Kawakami K, Kajitani T, Urano T, Murakami Y. DNA-RNA hybrid formation mediates RNAi-directed heterochromatin formation. Genes to cells: devoted to molecular & cellular mechanisms. 2012;17(3):218–33.
Ginno PA, Lim YW, Lott PL, Korf I, Chedin F. GC skew at the 5′ and 3′ ends of human genes links R-loop formation to epigenetic regulation and transcription termination. Genome Res. 2013;23(10):1590–600.
Ginno PA, Lott PL, Christensen HC, Korf I, Chedin F. R-loop formation is a distinctive characteristic of unmethylated human CpG island promoters. Mol Cell. 2012;45(6):814–25.
Boguslawski SJ, Smith DE, Michalak MA, Mickelson KE, Yehle CO, Patterson WL, et al. Characterization of monoclonal antibody to DNA.RNA and its application to immunodetection of hybrids. J Immunol Methods. 1986;89(1):123–30.
Phillips DD, Garboczi DN, Singh K, Hu Z, Leppla SH, Leysath CE. The sub-nanomolar binding of DNA-RNA hybrids by the single-chain Fv fragment of antibody S9.6. J Mol Recognit. 2013;26(8):376–81.
El Hage A, French SL, Beyer AL, Tollervey D. Loss of Topoisomerase I leads to R-loop-mediated transcriptional blocks during ribosomal RNA synthesis. Genes Dev. 2010;24(14):1546–58.
Han L, Masani S, Yu K. Overlapping activation-induced cytidine deaminase hotspot motifs in Ig class-switch recombination. Proc Natl Acad Sci U S A. 2011;108(28):11584–9.
Kao YP, Hsieh WC, Hung ST, Huang CW, Lieber MR, Huang FT. Detection and characterization of R-loops at the murine immunoglobulin Salpha region. Mol Immunol. 2013;54(2):208–16.
Zhang ZZ, Pannunzio NR, Han L, Hsieh CL, Yu K, Lieber MR. The Strength of an Ig Switch Region Is Determined by Its Ability to Drive R Loop Formation and Its Number of WGCW Sites. Cell Rep. 2014;8:557–69.
Zhang ZZ, Pannunzio NR, Hsieh CL, Yu K, Lieber MR. The role of G-density in switch region repeats for immunoglobulin class switch recombination. Nucleic Acids Res. 2014;42(21):13186–93.
Han L, Yu K. Altered kinetics of nonhomologous end joining and class switch recombination in ligase IV–deficient B cells. J Exp Med. 2008;205(12):2745–53.
Nakamura M, Kondo S, Sugai M, Nazarea M, Imamura S, Honjo T. High frequency class switching of an IgM+ B lymphoma clone CH12F3 to IgA+ cells. Int Immunol. 1996;8(2):193–201.
Hu Z, Zhang A, Storz G, Gottesman S, Leppla SH. An antibody-based microarray assay for small RNA detection. Nucleic Acids Res. 2006;34(7):e52.
We thank members of the Lieber and Yu labs for comments. This work was supported by NIH grants to MRL and KY.
The authors declare that they have no competing interests.
ZZZ, performed all of the experiments; NRP, purified the S9.6 antibody and performed unpublished gel shifts for S9.6 on defined nucleic acid species; CLH, advised on aspects of the study; KY, grew the S9.6 hybridoma and advised on the use of the CH12 cells; MRL, advised on all aspects of the study. All authors contributed to the writing of the paper. All authors read and approved the final manuscript.
Structure and Origins of Targets of the S9.6 Antibody. In addition to the stable R-loops at mammalian IgH switch regions (A), short R-loop-like structure may also form when RNA transcripts anneal back to the DNA template strand during DNA breathing (B), or when primases synthesize RNA primers during DNA replication (C). RNA with secondary structure may also be recognized by S9.6 (D). Table S1. Real-Time PCR Oligonucleotides.
About this article
- Immunogloblulin heavy chain class switch recombination
- S9.6 monoclonal antibody
- Nucleic acid structure