Best practices for the visualization, mapping, and manipulation of R-loops.

R-loops represent an abundant class of large non-B DNA structures in genomes. Even though they form transiently and at modest frequencies, interfering with R-loop formation or dissolution has significant impacts on genome stability. Addressing the mechanism(s) of R-loop-mediated genome destabilization requires a precise characterization of their distribution in genomes. A number of independent methods have been developed to visualize and map R-loops, but their results are at times discordant, leading to confusion. Here, we review the main existing methodologies for R-loop mapping and assess their limitations as well as the robustness of existing datasets. We offer a set of best practices to improve the reproducibility of maps, hoping that such guidelines could be useful for authors and referees alike. Finally, we propose a possible resolution for the apparent contradictions in R-loop mapping outcomes between antibody-based and RNase H1-based mapping approaches.

Exploration and analysis of R-loop mapping data with RLBase.

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. In 2012, Ginno et al. introduced the first R-loop mapping method. Since that time, dozens of R-loop mapping studies have been conducted, yielding hundreds of publicly available datasets. Current R-loop databases provide only limited access to these data. Moreover, no web tools for analyzing user-supplied R-loop datasets have yet been described. In our recent work, we reprocessed 810 R-loop mapping samples, building the largest R-loop data resource to date. We also defined R-loop consensus regions and developed a framework for R-loop data analysis. Now, we introduce RLBase, a user-friendly database that provides the capability to (i) explore hundreds of public R-loop mapping datasets, (ii) explore R-loop consensus regions, (iii) analyze user-supplied data and (iv) download standardized and reprocessed datasets. RLBase is directly accessible via the following URL: https://gccri.bishop-lab.uthscsa.edu/shiny/rlbase/.

Prevalent, Dynamic, and Conserved R-Loop Structures Associate with Specific Epigenomic Signatures in Mammals.

R-loops are three-stranded nucleic acid structures formed upon annealing of an RNA strand to one strand of duplex DNA. We profiled R-loops using a high-resolution, strand-specific methodology in human and mouse cell types. R-loops are prevalent, collectively occupying up to 5% of mammalian genomes. R-loop formation occurs over conserved genic hotspots such as promoter and terminator regions of poly(A)-dependent genes. In most cases, R-loops occur co-transcriptionally and undergo dynamic turnover. Detailed epigenomic profiling revealed that R-loops associate with specific chromatin signatures. At promoters, R-loops associate with a hyper-accessible state characteristic of unmethylated CpG island promoters. By contrast, terminal R-loops associate with an enhancer- and insulator-like state and define a broad class of transcription terminators. Together, this suggests that the retention of nascent RNA transcripts at their site of expression represents an abundant, dynamic, and programmed component of the mammalian chromatin that affects chromatin patterning and the control of gene expression.

Genomic patterns of transcription-replication interactions in mouse primary B cells.

Conflicts between transcription and replication machinery are a potent source of replication stress and genome instability; however, no technique currently exists to identify endogenous genomic locations prone to transcription-replication interactions. Here, we report a novel method to identify genomic loci prone to transcription-replication interactions termed transcription-replication immunoprecipitation on nascent DNA sequencing, TRIPn-Seq. TRIPn-Seq employs the sequential immunoprecipitation of RNA polymerase 2 phosphorylated at serine 5 (RNAP2s5) followed by enrichment of nascent DNA previously labeled with bromodeoxyuridine. Using TRIPn-Seq, we mapped 1009 unique transcription-replication interactions (TRIs) in mouse primary B cells characterized by a bimodal pattern of RNAP2s5, bidirectional transcription, an enrichment of RNA:DNA hybrids, and a high probability of forming G-quadruplexes. TRIs are highly enriched at transcription start sites and map to early replicating regions. TRIs exhibit enhanced Replication Protein A association and TRI-associated genes exhibit higher replication fork termination than control transcription start sites, two marks of replication stress. TRIs colocalize with double-strand DNA breaks, are enriched for deletions, and accumulate mutations in tumors. We propose that replication stress at TRIs induces mutations potentially contributing to age-related disease, as well as tumor formation and development.

Quality-controlled R-loop meta-analysis reveals the characteristics of R-loop consensus regions.

R-loops are three-stranded nucleic acid structures formed from the hybridization of RNA and DNA. While the pathological consequences of R-loops have been well-studied to date, the locations, classes, and dynamics of physiological R-loops remain poorly understood. R-loop mapping studies provide insight into R-loop dynamics, but their findings are challenging to generalize. This is due to the narrow biological scope of individual studies, the limitations of each mapping modality, and, in some cases, poor data quality. In this study, we reprocessed 810 R-loop mapping datasets from a wide array of biological conditions and mapping modalities. From this data resource, we developed an accurate R-loop data quality control method, and we reveal the extent of poor-quality data within previously published studies. We then identified a set of high-confidence R-loop mapping samples and used them to define consensus R-loop sites called 'R-loop regions' (RL regions). In the process, we identified a stark divergence between RL regions detected by S9.6 and dRNH-based mapping methods, particularly with respect to R-loop size, location, and colocalization with RNA binding factors. Taken together, this work provides a much-needed method to assess R-loop data quality and offers novel context regarding the differences between dRNH- and S9.6-based R-loop mapping approaches.

Arginine methylation facilitates the recruitment of TOP3B to chromatin to prevent R loop accumulation.

Tudor domain-containing protein 3 (TDRD3) is a major methylarginine effector molecule that reads methyl-histone marks and facilitates gene transcription. However, the underlying mechanism by which TDRD3 functions as a transcriptional coactivator is unknown. We identified topoisomerase IIIB (TOP3B) as a component of the TDRD3 complex. TDRD3 serves as a molecular bridge between TOP3B and arginine-methylated histones. The TDRD3-TOP3B complex is recruited to the c-MYC gene promoter primarily by the H4R3me2a mark, and the complex promotes c-MYC gene expression. TOP3B relaxes negative supercoiled DNA and reduces transcription-generated R loops in vitro. TDRD3 knockdown in cells increases R loop formation at the c-MYC locus, and Tdrd3 null mice exhibit elevated R loop formation at this locus in B cells. Tdrd3 null mice show significantly increased c-Myc/Igh translocation, a process driven by R loop structures. By reducing negative supercoiling and resolving R loops, TOP3B promotes transcription, protects against DNA damage, and reduces the frequency of chromosomal translocations.

Interplay between DNA sequence and negative superhelicity drives R-loop structures.

R-loops are abundant three-stranded nucleic-acid structures that form in cis during transcription. Experimental evidence suggests that R-loop formation is affected by DNA sequence and topology. However, the exact manner by which these factors interact to determine R-loop susceptibility is unclear. To investigate this, we developed a statistical mechanical equilibrium model of R-loop formation in superhelical DNA. In this model, the energy involved in forming an R-loop includes four terms-junctional and base-pairing energies and energies associated with superhelicity and with the torsional winding of the displaced DNA single strand around the RNA:DNA hybrid. This model shows that the significant energy barrier imposed by the formation of junctions can be overcome in two ways. First, base-pairing energy can favor RNA:DNA over DNA:DNA duplexes in favorable sequences. Second, R-loops, by absorbing negative superhelicity, partially or fully relax the rest of the DNA domain, thereby returning it to a lower energy state. In vitro transcription assays confirmed that R-loops cause plasmid relaxation and that negative superhelicity is required for R-loops to form, even in a favorable region. Single-molecule R-loop footprinting following in vitro transcription showed a strong agreement between theoretical predictions and experimental mapping of stable R-loop positions and further revealed the impact of DNA topology on the R-loop distribution landscape. Our results clarify the interplay between base sequence and DNA superhelicity in controlling R-loop stability. They also reveal R-loops as powerful and reversible topology sinks that cells may use to nonenzymatically relieve superhelical stress during transcription.

Emerging roles for R-loop structures in the management of topological stress.

R-loop structures are a prevalent class of alternative non-B DNA structures that form during transcription upon invasion of the DNA template by the nascent RNA. R-loops form universally in the genomes of organisms ranging from bacteriophages, bacteria, and yeasts to plants and animals, including mammals. A growing body of work has linked these structures to both physiological and pathological processes, in particular to genome instability. The rising interest in R-loops is placing new emphasis on understanding the fundamental physicochemical forces driving their formation and stability. Pioneering work in Escherichia coli revealed that DNA topology, in particular negative DNA superhelicity, plays a key role in driving R-loops. A clear role for DNA sequence was later uncovered. Here, we review and synthesize available evidence on the roles of DNA sequence and DNA topology in controlling R-loop formation and stability. Factoring in recent developments in R-loop modeling and single-molecule profiling, we propose a coherent model accounting for the interplay between DNA sequence and DNA topology in driving R-loop structure formation. This model reveals R-loops in a new light as powerful and reversible topological stress relievers, an insight that significantly expands the repertoire of R-loops' potential biological roles under both normal and aberrant conditions.

The extruded non-template strand determines the architecture of R-loops.

Three-stranded R-loop structures have been associated with genomic instability phenotypes. What underlies their wide-ranging effects on genome stability remains poorly understood. Here we combined biochemical and atomic force microscopy approaches with single molecule R-loop footprinting to demonstrate that R-loops formed at the model Airn locus in vitro adopt a defined set of three-dimensional conformations characterized by distinct shapes and volumes, which we call R-loop objects. Interestingly, we show that these R-loop objects impose specific physical constraints on the DNA, as revealed by the presence of stereotypical angles in the surrounding DNA. Biochemical probing and mutagenesis experiments revealed that the formation of R-loop objects at Airn is dictated by the extruded non-template strand, suggesting that R-loops possess intrinsic sequence-driven properties. Consistent with this, we show that R-loops formed at the fission yeast gene sum3 do not form detectable R-loop objects. Our results reveal that R-loops differ by their architectures and that the organization of the non-template strand is a fundamental characteristic of R-loops, which could explain that only a subset of R-loops is associated with replication-dependent DNA breaks.

R-loop formation is a distinctive characteristic of unmethylated human CpG island promoters.

CpG islands (CGIs) function as promoters for approximately 60% of human genes. Most of these elements remain protected from CpG methylation, a prevalent epigenetic modification associated with transcriptional silencing. Here, we report that methylation-resistant CGI promoters are characterized by significant strand asymmetry in the distribution of guanines and cytosines (GC skew) immediately downstream from their transcription start sites. Using innovative genomics methodologies, we show that transcription through regions of GC skew leads to the formation of long R loop structures. Furthermore, we show that GC skew and R loop formation potential is correlated with and predictive of the unmethylated state of CGIs. Finally, we provide evidence that R loop formation protects from DNMT3B1, the primary de novo DNA methyltransferase in early development. Altogether, these results suggest that protection from DNA methylation is a built-in characteristic of the DNA sequence of CGI promoters that is revealed by the cotranscriptional formation of R loop structures.

Nascent Connections: R-Loops and Chromatin Patterning.

RNA molecules, such as long noncoding RNAs (lncRNAs), have critical roles in regulating gene expression, chromosome architecture, and the modification states of chromatin. Recent developments suggest that RNA also influences gene expression and chromatin patterns through the interaction of nascent transcripts with their DNA template via the formation of co-transcriptional R-loop structures. R-loop formation over specific, conserved, hotspots occurs at thousands of genes in mammalian genomes and represents an important and dynamic feature of mammalian chromatin. Here, focusing primarily on mammalian systems, I describe the accumulating connections and possible mechanisms linking R-loop formation and chromatin patterning. The possible contribution of aberrant R-loops to pathological conditions is also discussed.

Comprehensive Proteomic Analysis of Mesenchymal Stem Cell Exosomes Reveals Modulation of Angiogenesis via Nuclear Factor-KappaB Signaling.

Mesenchymal stem cells (MSC) are known to facilitate healing of ischemic tissue related diseases through proangiogenic secretory proteins. Recent studies further show that MSC derived exosomes function as paracrine effectors of angiogenesis, however, the identity of which components of the exosome proteome responsible for this effect remains elusive. To address this we used high-resolution isoelectric focusing coupled liquid chromatography tandem mass spectrometry, an unbiased high throughput proteomics approach to comprehensively characterize the proteinaceous contents of MSCs and MSC derived exosomes. We probed the proteome of MSCs and MSC derived exosomes from cells cultured under expansion conditions and under ischemic tissue simulated conditions to elucidate key angiogenic paracrine effectors present and potentially differentially expressed in these conditions. In total, 6,342 proteins were identified in MSCs and 1,927 proteins in MSC derived exosomes, representing to our knowledge the first time these proteomes have been probed comprehensively. Multilayered analyses identified several putative paracrine effectors of angiogenesis present in MSC exosomes and increased in expression in MSCs exposed to ischemic tissue-simulated conditions; these include platelet derived growth factor, epidermal growth factor, fibroblast growth factor, and most notably nuclear factor-kappaB (NFkB) signaling pathway proteins. NFkB signaling was identified as a key mediator of MSC exosome induced angiogenesis in endothelial cells by functional in vitro validation using a specific inhibitor. Collectively, the results of our proteomic analysis show that MSC derived exosomes contain a robust profile of angiogenic paracrine effectors, which have potential for the treatment of ischemic tissue-related diseases.

GC skew is a conserved property of unmethylated CpG island promoters across vertebrates.

GC skew is a measure of the strand asymmetry in the distribution of guanines and cytosines. GC skew favors R-loops, a type of three stranded nucleic acid structures that form upon annealing of an RNA strand to one strand of DNA, creating a persistent RNA:DNA hybrid. Previous studies show that GC skew is prevalent at thousands of human CpG island (CGI) promoters and transcription termination regions, which correspond to hotspots of R-loop formation. Here, we investigated the conservation of GC skew patterns in 60 sequenced chordates genomes. We report that GC skew is a conserved sequence characteristic of the CGI promoter class in vertebrates. Furthermore, we reveal that promoter GC skew peaks at the exon 1/ intron1 junction and that it is highly correlated with gene age and CGI promoter strength. Our data also show that GC skew is predictive of unmethylated CGI promoters in a range of vertebrate species and that it imparts significant DNA hypomethylation for promoters with intermediate CpG densities. Finally, we observed that terminal GC skew is conserved for a subset of vertebrate genes that tend to be located significantly closer to their downstream neighbors, consistent with a role for R-loop formation in transcription termination.

Transcription-associated R-loop formation across the human FMR1 CGG-repeat region.

Expansion of a trinucleotide (CGG) repeat element within the 5' untranslated region (5'UTR) of the human FMR1 gene is responsible for a number of heritable disorders operating through distinct pathogenic mechanisms: gene silencing for fragile X syndrome (>200 CGG) and RNA toxic gain-of-function for FXTAS (∼ 55-200 CGG). Existing models have focused almost exclusively on post-transcriptional mechanisms, but co-transcriptional processes could also contribute to the molecular dysfunction of FMR1. We have observed that transcription through the GC-rich FMR1 5'UTR region favors R-loop formation, with the nascent (G-rich) RNA forming a stable RNA:DNA hybrid with the template DNA strand, thereby displacing the non-template DNA strand. Using DNA:RNA (hybrid) immunoprecipitation (DRIP) of genomic DNA from cultured human dermal fibroblasts with both normal (∼ 30 CGG repeats) and premutation (55

GC skew at the 5' and 3' ends of human genes links R-loop formation to epigenetic regulation and transcription termination.

Strand asymmetry in the distribution of guanines and cytosines, measured by GC skew, predisposes DNA sequences toward R-loop formation upon transcription. Previous work revealed that GC skew and R-loop formation associate with a core set of unmethylated CpG island (CGI) promoters in the human genome. Here, we show that GC skew can distinguish four classes of promoters, including three types of CGI promoters, each associated with unique epigenetic and gene ontology signatures. In particular, we identify a strong and a weak class of CGI promoters and show that these loci are enriched in distinct chromosomal territories reflecting the intrinsic strength of their protection against DNA methylation. Interestingly, we show that strong CGI promoters are depleted from the X chromosome while weak CGIs are enriched, a property consistent with the acquisition of DNA methylation during dosage compensation. Furthermore, we identify a third class of CGI promoters based on its unique GC skew profile and show that this gene set is enriched for Polycomb group targets. Lastly, we show that nearly 2000 genes harbor GC skew at their 3' ends and that these genes are preferentially located in gene-dense regions and tend to be closely arranged. Genomic profiling of R-loops accordingly showed that a large proportion of genes with terminal GC skew form R-loops at their 3' ends, consistent with a role for these structures in permitting efficient transcription termination. Altogether, we show that GC skew and R-loop formation offer significant insights into the epigenetic regulation, genomic organization, and function of human genes.

R-loop formation at Snord116 mediates topotecan inhibition of Ube3a-antisense and allele-specific chromatin decondensation.

Prader-Willi syndrome (PWS) and Angelman syndrome (AS) are oppositely imprinted autism-spectrum disorders with known genetic bases, but complex epigenetic mechanisms underlie their pathogenesis. The PWS/AS locus on 15q11-q13 is regulated by an imprinting control region that is maternally methylated and silenced. The PWS imprinting control region is the promoter for a one megabase paternal transcript encoding the ubiquitous protein-coding Snrpn gene and multiple neuron-specific noncoding RNAs, including the PWS-related Snord116 repetitive locus of small nucleolar RNAs and host genes, and the antisense transcript to AS-causing ubiquitin ligase encoding Ube3a (Ube3a-ATS). Neuron-specific transcriptional progression through Ube3a-ATS correlates with paternal Ube3a silencing and chromatin decondensation. Interestingly, topoisomerase inhibitors, including topotecan, were recently identified in an unbiased drug screen for compounds that could reverse the silent paternal allele of Ube3a in neurons, but the mechanism of topotecan action on the PWS/AS locus is unknown. Here, we demonstrate that topotecan treatment stabilizes the formation of RNA:DNA hybrids (R loops) at G-skewed repeat elements within paternal Snord116, corresponding to increased chromatin decondensation and inhibition of Ube3a-ATS expression. Neural precursor cells from paternal Snord116 deletion mice exhibit increased Ube3a-ATS levels in differentiated neurons and show a reduced effect of topotecan compared with wild-type neurons. These results demonstrate that the AS candidate drug topotecan acts predominantly through stabilizing R loops and chromatin decondensation at the paternally expressed PWS Snord116 locus. Our study holds promise for targeted therapies to the Snord116 locus for both AS and PWS.

Accelerated DNA replication fork speed due to loss of R-loops in myelodysplastic syndromes with SF3B1 mutation.

Myelodysplastic syndromes (MDS) with mutated SF3B1 gene present features including a favourable outcome distinct from MDS with mutations in other splicing factor genes SRSF2 or U2AF1. Molecular bases of these divergences are poorly understood. Here we find that SF3B1-mutated MDS show reduced R-loop formation predominating in gene bodies associated with intron retention reduction, not found in U2AF1- or SRSF2-mutated MDS. Compared to erythroblasts from SRSF2- or U2AF1-mutated patients, SF3B1-mutated erythroblasts exhibit augmented DNA synthesis, accelerated replication forks, and single-stranded DNA exposure upon differentiation. Importantly, histone deacetylase inhibition using vorinostat restores R-loop formation, slows down DNA replication forks and improves SF3B1-mutated erythroblast differentiation. In conclusion, loss of R-loops with associated DNA replication stress represents a hallmark of SF3B1-mutated MDS ineffective erythropoiesis, which could be used as a therapeutic target.

DNMT3L modulates significant and distinct flanking sequence preference for DNA methylation by DNMT3A and DNMT3B in vivo.

The DNTM3A and DNMT3B de novo DNA methyltransferases (DNMTs) are responsible for setting genomic DNA methylation patterns, a key layer of epigenetic information. Here, using an in vivo episomal methylation assay and extensive bisulfite methylation sequencing, we show that human DNMT3A and DNMT3B possess significant and distinct flanking sequence preferences for target CpG sites. Selection for high or low efficiency sites is mediated by the base composition at the -2 and +2 positions flanking the CpG site for DNMT3A, and at the -1 and +1 positions for DNMT3B. This intrinsic preference reproducibly leads to the formation of specific de novo methylation patterns characterized by up to 34-fold variations in the efficiency of DNA methylation at individual sites. Furthermore, analysis of the distribution of signature methylation hotspot and coldspot motifs suggests that DNMT flanking sequence preference has contributed to shaping the composition of CpG islands in the human genome. Our results also show that the DNMT3L stimulatory factor modulates the formation of de novo methylation patterns in two ways. First, DNMT3L selectively focuses the DNA methylation machinery on properly chromatinized DNA templates. Second, DNMT3L attenuates the impact of the intrinsic DNMT flanking sequence preference by providing a much greater boost to the methylation of poorly methylated sites, thus promoting the formation of broader and more uniform methylation patterns. This study offers insights into the manner by which DNA methylation patterns are deposited and reveals a new level of interplay between members of the de novo DNMT family.

A SPOT on the chromatin landscape? Histone peptide arrays as a tool for epigenetic research.

Post-translational modifications of histones serve as docking sites and signals for effector proteins and chromatin-remodeling enzymes, thereby influencing many fundamental cellular processes. Nevertheless, there are huge gaps in the knowledge of which proteins read and write the 'histone code'. Several techniques have been used to decipher complex histone-modification patterns. However, none is entirely satisfactory owing to the inherent limitations of in vitro studies of histones, such as deficits in the knowledge of the proteins involved, and the associated difficulties in the consistent and quantitative generation of histone marks. An alternative technique that could prove to be a useful tool in the study of the histone code is the use of synthetic peptide arrays (SPOT blot analysis) as a screening approach to characterize macromolecules that interact with specific covalent modifications of histone tails.

Recognition of multivalent histone states associated with heterochromatin by UHRF1 protein.

Histone modifications and DNA methylation represent two layers of heritable epigenetic information that regulate eukaryotic chromatin structure and gene activity. UHRF1 is a unique factor that bridges these two layers; it is required for maintenance DNA methylation at hemimethylated CpG sites, which are specifically recognized through its SRA domain and also interacts with histone H3 trimethylated on lysine 9 (H3K9me3) in an unspecified manner. Here we show that UHRF1 contains a tandem Tudor domain (TTD) that recognizes H3 tail peptides with the heterochromatin-associated modification state of trimethylated lysine 9 and unmodified lysine 4 (H3K4me0/K9me3). Solution NMR and crystallographic data reveal the TTD simultaneously recognizes H3K9me3 through a conserved aromatic cage in the first Tudor subdomain and unmodified H3K4 within a groove between the tandem subdomains. The subdomains undergo a conformational adjustment upon peptide binding, distinct from previously reported mechanisms for dual histone mark recognition. Mutant UHRF1 protein deficient for H3K4me0/K9me3 binding shows altered localization to heterochromatic chromocenters and fails to reduce expression of a target gene, p16(INK4A), when overexpressed. Our results demonstrate a novel recognition mechanism for the combinatorial readout of histone modification states associated with gene silencing and add to the growing evidence for coordination of, and cross-talk between, the modification states of H3K4 and H3K9 in regulation of gene expression.