Senataxin and RNase H2 act redundantly to suppress genome instability during class switch recombination.

Class switch recombination generates distinct antibody isotypes critical to a robust adaptive immune system, and defects are associated with autoimmune disorders and lymphomagenesis. Transcription is required during class switch recombination to recruit the cytidine deaminase AID-an essential step for the formation of DNA double-strand breaks-and strongly induces the formation of R loops within the immunoglobulin heavy-chain locus. However, the impact of R loops on double-strand break formation and repair during class switch recombination remains unclear. Here, we report that cells lacking two enzymes involved in R loop removal-senataxin and RNase H2-exhibit increased R loop formation and genome instability at the immunoglobulin heavy-chain locus without impacting its transcriptional activity, AID recruitment, or class switch recombination efficiency. Senataxin and RNase H2-deficient cells also exhibit increased insertion mutations at switch junctions, a hallmark of alternative end joining. Importantly, these phenotypes were not observed in cells lacking senataxin or RNase H2B alone. We propose that senataxin acts redundantly with RNase H2 to mediate timely R loop removal, promoting efficient repair while suppressing AID-dependent genome instability and insertional mutagenesis.

The immune system is a complex network of cells and molecules, which helps to protect the body from invaders. The adaptive immune system can recognise millions of assailants, kill them, and 'learn' from this experience to mount an even quicker defence the next time the body is infected. To achieve this level of protection, specific immune cells, called B cells, divide when they come into contact with a molecule from a foreign particle, the antigen. The cloned B cells then produce millions of protective proteins, the antibodies, which patrol the blood stream and tag harmful particles for destruction. An antibody resembles a Y-shaped structure that contains a 'variable' region, which gives it the specificity to interact with an antigen, and a 'constant' region, which interacts with components of the immune system and determines the mechanisms used to destroy a pathogen. Based on the constant region, antibodies can be divided into five main classes. B cells are able to switch their production from one antibody class to another in an event known as class switch recombination, by making changes to the constant region. They do this by cutting out a portion of the genes for the constant region from their DNA and fusing the remaining DNA. The resulting antibodies still recognise the same target, but interact with different components of the immune system, ensuring that all the body's forces are mobilised. R-loops are temporary structures that form when a cell 'reads' the instructions in its DNA to make proteins. R-loops provide physical support by anchoring the transcription template to the DNA. They help control the activity of genes, but if they stay on the DNA for too long they could interfere with any form of. DNA repair ­ including the cutting and fusing mechanisms during class switch recombination. To find out more about this process, Zhao et al. used B-cells from mice lacking two specific proteins that usually help to remove R-loops. Without these proteins, the B cells generated more R-loops than normal. Nevertheless, the B-cells were able to undergo class switch recombination, even though their chromosomes showed large areas of DNA damage, and DNA sections that had been repaired contained several mistakes. Errors that occur during class switch recombination have been linked to immune disorders and B cell cancers. The study of Zhao et al. shows that even if R-loops do not affect some processes in B cells, they could still impact the overall health of their DNA. A next step would be to test if an inability to remove R-loops could indeed play a role in immune disorders and B-cell cancers.

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.

PRMT1-dependent regulation of RNA metabolism and DNA damage response sustains pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive cancer that has remained clinically challenging to manage. Here we employ an RNAi-based in vivo functional genomics platform to determine epigenetic vulnerabilities across a panel of patient-derived PDAC models. Through this, we identify protein arginine methyltransferase 1 (PRMT1) as a critical dependency required for PDAC maintenance. Genetic and pharmacological studies validate the role of PRMT1 in maintaining PDAC growth. Mechanistically, using proteomic and transcriptomic analyses, we demonstrate that global inhibition of asymmetric arginine methylation impairs RNA metabolism, which includes RNA splicing, alternative polyadenylation, and transcription termination. This triggers a robust downregulation of multiple pathways involved in the DNA damage response, thereby promoting genomic instability and inhibiting tumor growth. Taken together, our data support PRMT1 as a compelling target in PDAC and informs a mechanism-based translational strategy for future therapeutic development.Statement of significancePDAC is a highly lethal cancer with limited therapeutic options. This study identified and characterized PRMT1-dependent regulation of RNA metabolism and coordination of key cellular processes required for PDAC tumor growth, defining a mechanism-based translational hypothesis for PRMT1 inhibitors.

Recognition of RNA by the S9.6 antibody creates pervasive artifacts when imaging RNA:DNA hybrids.

The S9.6 antibody is broadly used to detect RNA:DNA hybrids but has significant affinity for double-stranded RNA. The impact of this off-target RNA binding activity has not been thoroughly investigated, especially in the context of immunofluorescence microscopy. We report that S9.6 immunofluorescence signal observed in fixed human cells arises predominantly from ribosomal RNA, not RNA:DNA hybrids. S9.6 staining was unchanged by pretreatment with the RNA:DNA hybrid-specific nuclease RNase H1, despite verification in situ that S9.6 recognized RNA:DNA hybrids and that RNase H1 was active. S9.6 staining was, however, significantly sensitive to RNase T1, which specifically degrades RNA. Additional imaging and biochemical data indicate that the prominent cytoplasmic and nucleolar S9.6 signal primarily derives from ribosomal RNA. Importantly, genome-wide maps obtained by DNA sequencing after S9.6-mediated DNA:RNA immunoprecipitation (DRIP) are RNase H1 sensitive and RNase T1 insensitive. Altogether, these data demonstrate that imaging using S9.6 is subject to pervasive artifacts without pretreatments and controls that mitigate its promiscuous recognition of cellular RNAs.

CARM1 Inhibition Enables Immunotherapy of Resistant Tumors by Dual Action on Tumor Cells and T Cells.

A number of cancer drugs activate innate immune pathways in tumor cells but unfortunately also compromise antitumor immune function. We discovered that inhibition of CARM1, an epigenetic enzyme and cotranscriptional activator, elicited beneficial antitumor activity in both cytotoxic T cells and tumor cells. In T cells, Carm1 inactivation substantially enhanced their antitumor function and preserved memory-like populations required for sustained antitumor immunity. In tumor cells, Carm1 inactivation induced a potent type 1 interferon response that sensitized resistant tumors to cytotoxic T cells. Substantially increased numbers of dendritic cells, CD8 T cells, and natural killer cells were present in Carm1-deficient tumors, and infiltrating CD8 T cells expressed low levels of exhaustion markers. Targeting of CARM1 with a small molecule elicited potent antitumor immunity and sensitized resistant tumors to checkpoint blockade. Targeting of this cotranscriptional regulator thus offers an opportunity to enhance immune function while simultaneously sensitizing resistant tumor cells to immune attack. SIGNIFICANCE: Resistance to cancer immunotherapy remains a major challenge. Targeting of CARM1 enables immunotherapy of resistant tumors by enhancing T-cell functionality and preserving memory-like T-cell populations within tumors. CARM1 inhibition also sensitizes resistant tumor cells to immune attack by inducing a tumor cell-intrinsic type 1 interferon response.This article is highlighted in the In This Issue feature, p. 1861.

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.

Bisulfite treatment and single-molecule real-time sequencing reveal D-loop length, position, and distribution.

Displacement loops (D-loops) are signature intermediates formed during homologous recombination. Numerous factors regulate D-loop formation and disruption, thereby influencing crucial aspects of DNA repair, including donor choice and the possibility of crossover outcome. While D-loop detection methods exist, it is currently unfeasible to assess the relationship between D-loop editors and D-loop characteristics such as length and position. Here, we developed a novel in vitro assay to characterize the length and position of individual D-loops with near base-pair resolution and deep coverage, while also revealing their distribution in a population. Non-denaturing bisulfite treatment modifies the cytosines on the displaced strand of the D-loop to uracil, leaving a permanent signature for the displaced strand. Subsequent single-molecule real-time sequencing uncovers the cytosine conversion patch as a D-loop footprint. The D-loop Mapping Assay is widely applicable with different substrates and donor types and can be used to study factors that influence D-loop properties.

Ultra-deep Coverage Single-molecule R-loop Footprinting Reveals Principles of R-loop Formation.

R-loops are a prevalent class of non-B DNA structures that have been associated with both positive and negative cellular outcomes. DNA:RNA immunoprecipitation (DRIP) approaches based on the anti-DNA:RNA hybrid S9.6 antibody revealed that R-loops form dynamically over conserved genic hotspots. We have developed an orthogonal approach that queries R-loops via the presence of long stretches of single-stranded DNA on their looped-out strand. Nondenaturing 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 'footprints' at near nucleotide resolution in a strand-specific manner on long single DNA molecules and at ultra-deep coverage. Single-molecule R-loop footprinting coupled with PacBio sequencing (SMRF-seq) revealed a strong agreement between S9.6-based and bisulfite-based R-loop mapping and confirmed that R-loops form over genic hotspots, including gene bodies and terminal gene regions. Based on the largest single-molecule R-loop dataset to date, we show that individual R-loops form nonrandomly, defining discrete sets of overlapping molecular clusters that pileup through larger R-loop zones. R-loops most often map to intronic regions and their individual start and stop positions do not match with intron-exon boundaries, reinforcing the model that they form cotranscriptionally from unspliced transcripts. SMRF-seq further established that R-loop distribution patterns are not simply driven by intrinsic DNA sequence features but most likely also reflect DNA topological constraints. Overall, DRIP-based and SMRF-based approaches independently provide a complementary and congruent view of R-loop distribution, consolidating our understanding of the principles underlying R-loop formation.

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.

DNA Topoisomerase I differentially modulates R-loops across the human genome.

BACKGROUND: Co-transcriptional R-loops are abundant non-B DNA structures in mammalian genomes. DNA Topoisomerase I (Top1) is often thought to regulate R-loop formation owing to its ability to resolve both positive and negative supercoils. How Top1 regulates R-loop structures at a global level is unknown. RESULTS: Here, we perform high-resolution strand-specific R-loop mapping in human cells depleted for Top1 and find that Top1 depletion results in both R-loop gains and losses at thousands of transcribed loci, delineating two distinct gene classes. R-loop gains are characteristic for long, highly transcribed, genes located in gene-poor regions anchored to Lamin B1 domains and in proximity to H3K9me3-marked heterochromatic patches. R-loop losses, by contrast, occur in gene-rich regions overlapping H3K27me3-marked active replication initiation regions. Interestingly, Top1 depletion coincides with a block of the cell cycle in G0/G1 phase and a trend towards replication delay. CONCLUSIONS: Our findings reveal new properties of Top1 in regulating R-loop homeostasis in a context-dependent manner and suggest a potential role for Top1 in modulating the replication process via R-loop formation.

The Affinity of the S9.6 Antibody for Double-Stranded RNAs Impacts the Accurate Mapping of R-Loops in Fission Yeast.

R-loops, which result from the formation of stable DNA:RNA hybrids, can both threaten genome integrity and act as physiological regulators of gene expression and chromatin patterning. To characterize R-loops in fission yeast, we used the S9.6 antibody-based DRIPc-seq method to sequence the RNA strand of R-loops and obtain strand-specific R-loop maps at near nucleotide resolution. Surprisingly, preliminary DRIPc-seq experiments identified mostly RNase H-resistant but exosome-sensitive RNAs that mapped to both DNA strands and resembled RNA:RNA hybrids (dsRNAs), suggesting that dsRNAs form widely in fission yeast. We confirmed in vitro that S9.6 can immuno-precipitate dsRNAs and provide evidence that dsRNAs can interfere with its binding to R-loops. dsRNA elimination by RNase III treatment prior to DRIPc-seq allowed the genome-wide and strand-specific identification of genuine R-loops that responded in vivo to RNase H levels and displayed classical features associated with R-loop formation. We also found that most transcripts whose levels were altered by in vivo manipulation of RNase H levels did not form detectable R-loops, suggesting that prolonged manipulation of R-loop levels could indirectly alter the transcriptome. We discuss the implications of our work in the design of experimental strategies to probe R-loop functions.

Telomeres in ICF syndrome cells are vulnerable to DNA damage due to elevated DNA:RNA hybrids.

DNA:RNA hybrids, nucleic acid structures with diverse physiological functions, can disrupt genome integrity when dysregulated. Human telomeres were shown to form hybrids with the lncRNA TERRA, yet the formation and distribution of these hybrids among telomeres, their regulation and their cellular effects remain elusive. Here we predict and confirm in several human cell types that DNA:RNA hybrids form at many subtelomeric and telomeric regions. We demonstrate that ICF syndrome cells, which exhibit short telomeres and elevated TERRA levels, are enriched for hybrids at telomeric regions throughout the cell cycle. Telomeric hybrids are associated with high levels of DNA damage at chromosome ends in ICF cells, which are significantly reduced with overexpression of RNase H1. Our findings suggest that abnormally high TERRA levels in ICF syndrome lead to accumulation of telomeric hybrids that, in turn, can result in telomeric dysfunction.

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.

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.

Genome-wide DNA hypomethylation and RNA:DNA hybrid accumulation in Aicardi-Goutières syndrome.

Aicardi-Goutières syndrome (AGS) is a severe childhood inflammatory disorder that shows clinical and genetic overlap with systemic lupus erythematosus (SLE). AGS is thought to arise from the accumulation of incompletely metabolized endogenous nucleic acid species owing to mutations in nucleic acid-degrading enzymes TREX1 (AGS1), RNase H2 (AGS2, 3 and 4), and SAMHD1 (AGS5). However, the identity and source of such immunogenic nucleic acid species remain undefined. Using genome-wide approaches, we show that fibroblasts from AGS patients with AGS1-5 mutations are burdened by excessive loads of RNA:DNA hybrids. Using MethylC-seq, we show that AGS fibroblasts display pronounced and global loss of DNA methylation and demonstrate that AGS-specific RNA:DNA hybrids often occur within DNA hypomethylated regions. Altogether, our data suggest that RNA:DNA hybrids may represent a common immunogenic form of nucleic acids in AGS and provide the first evidence of epigenetic perturbations in AGS, furthering the links between AGS and SLE.

Sequencing the transcriptome of milk production: milk trumps mammary tissue.

BACKGROUND: Studies of normal human mammary gland development and function have mostly relied on cell culture, limited surgical specimens, and rodent models. Although RNA extracted from human milk has been used to assay the mammary transcriptome non-invasively, this assay has not been adequately validated in primates. Thus, the objectives of the current study were to assess the suitability of lactating rhesus macaques as a model for lactating humans and to determine whether RNA extracted from milk fractions is representative of RNA extracted from mammary tissue for the purpose of studying the transcriptome of milk-producing cells. RESULTS: We confirmed that macaque milk contains cytoplasmic crescents and that ample high-quality RNA can be obtained for sequencing. Using RNA sequencing, RNA extracted from macaque milk fat and milk cell fractions more accurately represented RNA from mammary epithelial cells (cells that produce milk) than did RNA from whole mammary tissue. Mammary epithelium-specific transcripts were more abundant in macaque milk fat, whereas adipose or stroma-specific transcripts were more abundant in mammary tissue. Functional analyses confirmed the validity of milk as a source of RNA from milk-producing mammary epithelial cells. CONCLUSIONS: RNA extracted from the milk fat during lactation accurately portrayed the RNA profile of milk-producing mammary epithelial cells in a non-human primate. However, this sample type clearly requires protocols that minimize RNA degradation. Overall, we validated the use of RNA extracted from human and macaque milk and provided evidence to support the use of lactating macaques as a model for human lactation.

Inactive DNMT3B splice variants modulate de novo DNA methylation.

Inactive DNA methyltransferase (DNMT) 3B splice isoforms are associated with changes in DNA methylation, yet the mechanisms by which they act remain largely unknown. Using biochemical and cell culture assays, we show here that the inactive DNMT3B3 and DNMT3B4 isoforms bind to and regulate the activity of catalytically competent DNMT3A or DNMT3B molecules. DNMT3B3 modestly stimulated the de novo methylation activity of DNMT3A and also counteracted the stimulatory effects of DNMT3L, therefore leading to subtle and contrasting effects on activity. DNMT3B4, by contrast, significantly inhibited de novo DNA methylation by active DNMT3 molecules, most likely due to its ability to reduce the DNA binding affinity of co-complexes, thereby sequestering them away from their substrate. Immunocytochemistry experiments revealed that in addition to their effects on the intrinsic catalytic function of active DNMT3 enzymes, DNMT3B3 and DNMT34 drive distinct types of chromatin compaction and patterns of histone 3 lysine 9 tri-methylation (H3K9me3) deposition. Our findings suggest that regulation of active DNMT3 members through the formation of co-complexes with inactive DNMT3 variants is a general mechanism by which DNMT3 variants function. This may account for some of the changes in DNA methylation patterns observed during development and disease.