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.

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.

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.

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.

High-resolution, strand-specific R-loop mapping via S9.6-based DNA-RNA immunoprecipitation and high-throughput sequencing.

R-loops are prevalent three-stranded non-B DNA structures composed of an RNA-DNA hybrid and a single strand of DNA. R-loops are implicated in various basic nuclear processes, such as class-switch recombination, transcription termination and chromatin patterning. Perturbations in R-loop metabolism have been linked to genomic instability and have been implicated in human disorders, including cancer. As a consequence, the accurate mapping of these structures has been of increasing interest in recent years. Here, we describe two related immunoprecipitation-based methods for mapping R-loop structures: basic DRIP-seq (DNA-RNA immunoprecipitation followed by high-throughput DNA sequencing), an easy, robust, but resolution-limited technique; and DRIPc-seq (DNA-RNA immunoprecipitation followed by cDNA conversion coupled to high-throughput sequencing), a high-resolution and strand-specific iteration of the method that permits accurate R-loop mapping genome wide. Briefly, after gentle DNA extraction and restriction digestion with a cocktail of enzymes, R-loop structures are immunoprecipitated with the anti-RNA-DNA hybrid S9.6 antibody. Compared with DRIP-seq, in which the immunoprecipitated DNA is directly sequenced, DRIPc-seq permits the recovery of the RNA moiety of R-loops, and these RNA strands are subjected to strand-specific RNA sequencing (RNA-seq) analysis. DRIPc-seq can be performed in 5 d and can be applied to any cell type, provided sufficient starting material can be collected. Accurately mapping R-loop distribution in various cell lines and under varied conditions is essential to understanding the formation, roles and dynamic resolution of these important structures.

Co-transcriptional R-loops are the main cause of estrogen-induced DNA damage.

The hormone estrogen (E2) binds the estrogen receptor to promote transcription of E2-responsive genes in the breast and other tissues. E2 also has links to genomic instability, and elevated E2 levels are tied to breast cancer. Here, we show that E2 stimulation causes a rapid, global increase in the formation of R-loops, co-transcriptional RNA-DNA products, which in some instances have been linked to DNA damage. We show that E2-dependent R-loop formation and breast cancer rearrangements are highly enriched at E2-responsive genomic loci and that E2 induces DNA replication-dependent double-strand breaks (DSBs). Strikingly, many DSBs that accumulate in response to E2 are R-loop dependent. Thus, R-loops resulting from the E2 transcriptional response are a significant source of DNA damage. This work reveals a novel mechanism by which E2 stimulation leads to genomic instability and highlights how transcriptional programs play an important role in shaping the genomic landscape of DNA damage susceptibility.

Histone methylation is mechanistically linked to DNA methylation at imprinting control regions in mammals.

Mono-allelic expression of imprinted genes from either the paternal or the maternal allele is mediated by imprinting control regions (ICRs), which are epigenetically marked in an allele-specific fashion. Although, in somatic cells, these epigenetic marks comprise both DNA methylation and histone methylation, the relationship between these two modifications in imprint acquisition and maintenance remains unclear. To address this important question, we analyzed histone modifications at ICRs in mid-gestation embryos that were obtained from Dnmt3L(-/-) females, in which DNA methylation imprints at ICRs are not established during oogenesis. The absence of maternal DNA methylation imprints in these conceptuses led to a marked decrease and loss of allele-specificity of the repressive H3K9me3, H4K20me3 and H2A/H4R3me2 histone modifications, providing the first evidence of a mechanistic link between DNA and histone methylation at ICRs. The existence of this relationship was strengthened by the observation that when DNA methylation was still present at the Snrpn and Peg3 ICRs in some of the progeny of Dnmt3L(-/-) females, these ICRs were associated with the usual patterns of histone methylation. Combined, our data establish that DNA methylation is involved in the acquisition and/or maintenance of histone methylation at ICRs.