High-Throughput Single-Molecule R-loop Footprinting Reveals Principles of R-loop Formation

R-loops are a prevalent class of non-B DNA structures that form during transcription upon reannealing of the nascent RNA to the template DNA strand. R-loops have been profiled using the S9.6 antibody to immunoprecipitate DNA:RNA hybrids. S9.6-based DNA:RNA immunoprecipitation (DRIP) techniques revealed that R-loops form dynamically over conserved genic hotspots. We developed an orthogonal profiling methodology that queries R-loops via the presence of long stretches of single-stranded DNA on the looped-out strand. Non-denaturing sodium bisulfite treatment catalyzes the conversion of unpaired cytosines to uracils, creating permanent genetic tags for the position of an R-loop. Long read, single-molecule PacBio sequencing allows the identification of R-loop 9footprints9 at near nucleotide resolution in a strand-specific manner on single DNA molecules and at ultra-deep coverage. Single-molecule R-loop footprinting (SMRF-seq) revealed a strong agreement between S9.6- and bisulfite-based R-loop mapping and confirmed that R-loops form from unspliced transcripts over genic hotspots. Using the largest single-molecule R-loop dataset to date, we show that individual R-loops generate overlapping sets of molecular clusters that pile-up through larger R-loop-prone zones. SMRF-seq further established that R-loop distribution patterns are driven by both intrinsic DNA sequence features and DNA topological constraints, revealing the principles of R-loop formation.