Chemical manipulation of m1A mediates its detection in human tRNA.

N 1-methyl adenosine (m1A) is a widespread RNA modification present in tRNA, rRNA, and mRNA. m1A modification sites in tRNAs are evolutionarily conserved and its formation on tRNA is catalyzed by methyltransferase TRMT61A and TRMT6 complex. m1A promotes translation initiation and elongation. Due to its positive charge under physiological conditions, m1A can notably modulate RNA structure. It also blocks Watson-Crick-Franklin base-pairing and causes mutation and truncation during reverse transcription. Several misincorporation-based high-throughput sequencing methods have been developed to sequence m1A. In this study, we introduce a reduction-based m1A sequencing (red-m1A-seq). We report that NaBH4 reduction of m1A can improve the mutation and readthrough rates using commercially available RT enzymes to give a better positive signature, while alkaline-catalyzed Dimroth rearrangement can efficiently convert m1A to m6A to provide good controls, allowing the detection of m1A with higher sensitivity and accuracy. We applied red-m1A-seq to sequence human small RNA, and we not only detected all the previously reported tRNA m1A sites, but also new m1A sites in mt-tRNAAsn-GTT and 5.8S rRNA.

5-Hydroxymethylcytosine Profiling of Cell-Free DNA Identifies Bivalent Genes That Are Prognostic of Survival in High-Risk Neuroblastoma.

PURPOSE: Neuroblastoma is the most common extracranial solid tumor in childhood. We previously showed that circulating cell-free DNA (cfDNA) and tumor biopsy derived 5-hydroxymethylcytosime (5-hmC) profiles identified patients with neuroblastoma who experienced subsequent relapse. Here, we hypothesized that 5-hmC modifications selectively enriched in cfDNA compared with tumor biopsy samples would identify epigenetic changes associated with aggressive tumor behavior and identify novel biomarkers of outcome in patients with high-risk neuroblastoma. METHODS: 5-hmC profiles from cfDNA (n = 64) and tumor biopsies (n = 48) were compared. Two neuroblastoma cell lines underwent chromatin immunoprecipitation followed by sequencing (ChIP-Seq) for H3K27me3, H3K4me3, and H3K27ac; kethoxal-associated single-stranded DNA sequencing; hmC-Seal for 5-hmC; and RNA-sequencing (RNA-Seq). Genes enriched for both H3K27me3 and H3K4me3 in the included cell lines were defined as bivalent. Using bivalent genes defined in vitro, a bivalent signature was established in three publicly available cohorts of patients with neuroblastoma through gene set variation analysis. Differences between tumors with high or low bivalent signatures were assessed by the Kaplan-Meier method and Cox proportional hazards models. RESULTS: In cfDNA compared with tumor biopsy derived 5-hmC profiles, we found increased 5-hmC deposition on Polycomb Repressive Complex 2 target genes, a finding previously described in the context of bivalent genes. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogeneous neuroblastoma cell lines. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in patients with high-risk neuroblastoma. CONCLUSION: Low expression of bivalent genes is a biomarker of worse outcome in patients with high-risk neuroblastoma.

PTPN2 copper-sensing rapidly relays copper level fluctuations into EGFR/CREB activation and associated CTR1 transcriptional repression.

Fluxes in human intra- and extracellular copper levels recently garnered attention for roles in cellular signaling, including affecting levels of the signaling molecule cyclic adenosine monophosphate (cAMP). We herein applied an unbiased temporal evaluation of the whole-genome transcriptional activities modulated by fluctuations in copper levels to identify the copper sensor proteins responsible for driving these activities. We found that fluctuations in physiologically-relevant copper levels rapidly modulate EGFR/MAPK/ERK signal transduction and activation of the transcription factor cAMP response element-binding protein (CREB). Both intracellular and extracellular assays support Cu 1+ inhibition of the EGFR-phosphatase PTPN2 (and potentially the homologous PTPN1)-via direct ligation to the PTPN2 active site cysteine side chain-as the underlying mechanism of copper-stimulated EGFR signal transduction activation. Depletion of copper represses this signaling pathway. We additionally show i ) copper supplementation drives transcriptional repression of the copper importer CTR1 and ii ) CREB activity is inversely correlated with CTR1 expression. In summary, our study reveals PTPN2 as a physiological copper sensor and defines a regulatory mechanism linking feedback control of copper-stimulated MAPK/ERK/CREB-signaling and CTR1 expression, thereby uncovering a previously unrecognized link between copper levels and cellular signal transduction.

KAS-Analyzer: a novel computational framework for exploring KAS-seq data.

Motivation: Kethoxal-assisted ssDNA sequencing (KAS-seq) is rapidly gaining popularity as a robust and effective approach to study the nascent dynamics of transcriptionally engaged RNA polymerases through profiling of genome-wide single-stranded DNA (ssDNA). Its latest variant, spKAS-seq, a strand-specific version of KAS-seq, has been developed to map genome-wide R-loop structures by detecting imbalances of ssDNA on two strands. However, user-friendly, open-source computational tools tailored for KAS-seq data are still lacking. Results: Here, we introduce KAS-Analyzer, the first comprehensive computational framework aimed at streamlining and enhancing the analysis and interpretation of KAS-seq and spKAS-seq data. In addition to standard analyses, KAS-Analyzer offers many novel tools specifically designed for KAS-seq data, including, but not limited to: calculation of transcription-related metrics, identification of single-stranded transcribing (SST) enhancers, high-resolution mapping of R-loops, and differential RNA polymerase activity analysis. We provided a detailed overview of KAS-seq data and its diverse applications through the implementation of KAS-Analyzer. Using the example time-course KAS-seq datasets, we further showcase the robust capabilities of KAS-Analyzer for investigating dynamic transcriptional regulatory programs in response to UVB radiation. Availability and implementation: KAS-Analyzer is available at https://github.com/Ruitulyu/KAS-Analyzer.

An engineered hypercompact CRISPR-Cas12f system with boosted gene-editing activity.

Compact CRISPR-Cas systems offer versatile treatment options for genetic disorders, but their application is often limited by modest gene-editing activity. Here we present enAsCas12f, an engineered RNA-guided DNA endonuclease up to 11.3-fold more potent than its parent protein, AsCas12f, and one-third of the size of SpCas9. enAsCas12f shows higher DNA cleavage activity than wild-type AsCas12f in vitro and functions broadly in human cells, delivering up to 69.8% insertions and deletions at user-specified genomic loci. Minimal off-target editing is observed with enAsCas12f, suggesting that boosted on-target activity does not impair genome-wide specificity. We determine the cryo-electron microscopy (cryo-EM) structure of the AsCas12f-sgRNA-DNA complex at a resolution of 2.9 Å, which reveals dimerization-mediated substrate recognition and cleavage. Structure-guided single guide RNA (sgRNA) engineering leads to sgRNA-v2, which is 33% shorter than the full-length sgRNA, but with on par activity. Together, the engineered hypercompact AsCas12f system enables robust and faithful gene editing in mammalian cells.

5-hydroxymethylcytosine profiling of cell-free DNA identifies bivalent genes that are prognostic of survival in high-risk neuroblastoma.

Neuroblastoma is the most common extra-cranial solid tumor in childhood and epigenetic dysregulation is a key driver of this embryonal disease. In cell-free DNA from neuroblastoma patients with high-risk disease, we found increased 5-hydroxymethylcytosine (5-hmC) deposition on Polycomb Repressive Complex 2 (PRC2) target genes, a finding previously described in the context of bivalent genes. As bivalent genes, defined as genes bearing both activating (H3K4me3) and repressive (H3K27me3) chromatin modifications, have been shown to play an important role in development and cancer, we investigated the potential role of bivalent genes in maintaining a de-differentiated state in neuroblastoma and their potential use as a biomarker. We identified 313 genes that bore bivalent chromatin marks, were enriched for mediators of neuronal differentiation, and were transcriptionally repressed across a panel of heterogenous neuroblastoma cell lines. Through gene set variance analysis, we developed a clinically implementable bivalent signature. In three distinct clinical cohorts, low bivalent signature was significantly and independently associated with worse clinical outcome in high-risk neuroblastoma patients. Thus, low expression of bivalent genes is a biomarker of ultra-high-risk disease and may represent a therapeutic opportunity in neuroblastoma.

Exon architecture controls mRNA m6A suppression and gene expression.

N6-methyladenosine (m6A) is the most abundant messenger RNA (mRNA) modification and plays crucial roles in diverse physiological processes. Using a massively parallel assay for m6A (MPm6A), we discover that m6A specificity is globally regulated by suppressors that prevent m6A deposition in unmethylated transcriptome regions. We identify exon junction complexes (EJCs) as m6A suppressors that protect exon junction-proximal RNA within coding sequences from methylation and regulate mRNA stability through m6A suppression. EJC suppression of m6A underlies multiple global characteristics of mRNA m6A specificity, with the local range of EJC protection sufficient to suppress m6A deposition in average-length internal exons but not in long internal and terminal exons. EJC-suppressed methylation sites colocalize with EJC-suppressed splice sites, which suggests that exon architecture broadly determines local mRNA accessibility to regulatory complexes.

spKAS-seq reveals R-loop dynamics using low-input materials by detecting single-stranded DNA with strand specificity.

R-loops affect transcription and genome stability. Dysregulation of R-loops is related to human diseases. Genome-wide R-loop mapping typically uses the S9.6 antibody or inactive ribonuclease H, both requiring a large number of cells with varying results observed depending on the approach applied. Here, we present strand-specific kethoxal-assisted single-stranded DNA (ssDNA) sequencing (spKAS-seq) to map R-loops by taking advantage of the presence of a ssDNA in the triplex structure. We show that spKAS-seq detects R-loops and their dynamics at coding sequences, enhancers, and other intergenic regions with as few as 50,000 cells. A joint analysis of R-loops and chromatin-bound RNA binding proteins (RBPs) suggested that R-loops can be RBP binding hotspots on the chromatin.

KAS-seq: genome-wide sequencing of single-stranded DNA by N3-kethoxal-assisted labeling.

Transcription and its dynamics are crucial for gene expression regulation. However, very few methods can directly read out transcriptional activity with low-input material and high temporal resolution. This protocol describes KAS-seq, a robust and sensitive approach for capturing genome-wide single-stranded DNA (ssDNA) profiles using N3-kethoxal-assisted labeling. We developed N3-kethoxal, an azido derivative of kethoxal that reacts with deoxyguanosine bases of ssDNA in live cells within 5-10 min at 37 °C, allowing the capture of dynamic changes. Downstream biotinylation of labeled DNA occurs via copper-free click chemistry. Altogether, the KAS-seq procedure involves N3-kethoxal labeling, DNA isolation, biotinylation, fragmentation, affinity pull-down, library preparation, sequencing and bioinformatics analysis. The pre-library construction labeling and enrichment can be completed in as little as 3-4 h and is applicable to both animal tissue and as few as 1,000 cultured cells. Our recent study shows that ssDNA signals measured by KAS-seq simultaneously reveal the dynamics of transcriptionally engaged RNA polymerase (Pol) II, transcribing enhancers, RNA Pol I and Pol III activities and potentially non-canonical DNA structures with high analytical sensitivity. In addition to the experimental protocol, we also introduce here KAS-pipe, a user-friendly integrative data analysis pipeline for KAS-seq.

Tumour suppressor TET2 safeguards enhancers from aberrant DNA methylation and epigenetic reprogramming in ERα-positive breast cancer cells.

Aberrant DNA methylation is an epigenetic hallmark of malignant tumours. The DNA methylation level is regulated by not only DNA methyltransferases (DNMTs) but also Ten-Eleven Translocation (TET) family proteins. However, the exact role of TET genes in breast cancer remains controversial. Here, we uncover that the ERα-positive breast cancer patients with high TET2 mRNA expression had better overall survival rates. Consistently, knockout of TET2 promotes the tumorigenesis of ERα-positive MCF7 breast cancer cells. Mechanistically, TET2 loss leads to aberrant DNA methylation (gain of 5mC) at a large proportion of enhancers, accompanied by significant reduction in H3K4me1 and H3K27ac enrichment. By analysing the epigenetically reprogrammed enhancers, we identify oestrogen responsive element (ERE) as one of the enriched motifs of transcriptional factors. Importantly, TET2 loss impairs 17beta-oestradiol (E2)-induced transcription of the epigenetically reprogrammed EREs-associated genes through attenuating the binding of ERα. Taken together, these findings shed light on our understanding of the epigenetic mechanisms underlying the enhancer reprogramming during breast cancer pathogenesis.

Regulation of TET2 gene expression and 5mC oxidation in breast cancer cells by estrogen signaling.

Estrogen signaling plays important roles in diverse physiological and pathophysiological processes. However, the relationship between estrogen signaling and epigenetic regulation is not fully understood. Here, we explored the effect of estrogen signaling on the expression of Ten-Eleven Translocation (TET) family genes and DNA hydroxylmethylation in estrogen receptor alpha positive (ERα+) breast cancer cells. By analyzing the RNA-seq data, we identified TET2 as an estradiol (E2)-responsive gene in ERα+ MCF7 cells. RT-qPCR and Western blot analyses confirmed that both the mRNA and protein levels of TET2 gene were upregulated in MCF7 cells by E2 treatment. ChIP-seq and qPCR analyses showed that the enrichment of ERα and H3K27ac on the upstream regulatory regions of TET2 gene was increased in MCF7 cells upon E2 treatment. Moreover, E2 treatment also led to a significant increase in the global 5-hydroxymethylcytosine (5hmC) level, while knockout of TET2 abolished such E2-induced 5hmC increase. Conversely, treatment with ICI 182780, a potent and selective estrogen receptor degrader (SERD), inhibited TET2 gene expression and down-regulated the 5hmC level in MCF7 cells. Taken together, our study identified an ERα/TET2/5hmC epigenetic pathway, which may participate in the estrogen-associated physiological and pathophysiological processes.

NSD2 dimethylation at H3K36 promotes lung adenocarcinoma pathogenesis.

The etiological role of NSD2 enzymatic activity in solid tumors is unclear. Here we show that NSD2, via H3K36me2 catalysis, cooperates with oncogenic KRAS signaling to drive lung adenocarcinoma (LUAD) pathogenesis. In vivo expression of NSD2E1099K, a hyperactive variant detected in individuals with LUAD, rapidly accelerates malignant tumor progression while decreasing survival in KRAS-driven LUAD mouse models. Pathologic H3K36me2 generation by NSD2 amplifies transcriptional output of KRAS and several complementary oncogenic gene expression programs. We establish a versatile in vivo CRISPRi-based system to test gene functions in LUAD and find that NSD2 loss strongly attenuates tumor progression. NSD2 knockdown also blocks neoplastic growth of PDXs (patient-dervived xenografts) from primary LUAD. Finally, a treatment regimen combining NSD2 depletion with MEK1/2 inhibition causes nearly complete regression of LUAD tumors. Our work identifies NSD2 as a bona fide LUAD therapeutic target and suggests a pivotal epigenetic role of the NSD2-H3K36me2 axis in sustaining oncogenic signaling.

TET2 Inhibits PD-L1 Gene Expression in Breast Cancer Cells through Histone Deacetylation.

Activation of PD-1/PD-L1 checkpoint is a critical step for the immune evasion of malignant tumors including breast cancer. However, the epigenetic mechanism underlying the aberrant expression of PD-L1 in breast cancer cells remains poorly understood. To investigate the role of TET2 in the regulation of PD-L1 gene expression, quantitative reverse transcription PCR (RT-qPCR), Western blotting, chromatin immunoprecipitation (ChIP) assay and MeDIP/hMeDIP-qPCR were performed on MCF7 and MDA-MB-231 human breast cancer cells. Here, we reported that TET2 depletion upregulated PD-L1 gene expression in MCF7 cells. Conversely, ectopic expression of TET2 inhibited PD-L1 gene expression in MDA-MB-231 cells. Mechanistically, TET2 protein recruits histone deacetylases (HDACs) to PD-L1 gene promoter and orchestrates a repressive chromatin structure to suppress PD-L1 gene transcription, which is likely independent of DNA demethylation. Consistently, treatment with HDAC inhibitors upregulated PD-L1 gene expression in wild-type (WT) but not TET2 KO MCF7 cells. Furthermore, analysis of the CCLE and TCGA data showed a negative correlation between TET2 and PD-L1 expression in breast cancer. Taken together, our results identify a new epigenetic regulatory mechanism of PD-L1 gene transcription, linking the catalytic activity-independent role of TET2 to the anti-tumor immunity in breast cancer.

Control of Early B Cell Development by the RNA N6-Methyladenosine Methylation.

The RNA N6-methyladenosine (m6A) methylation is installed by the METTL3-METTL14 methyltransferase complex. This modification has critical regulatory roles in various biological processes. Here, we report that deletion of Mettl14 dramatically reduces mRNA m6A methylation in developing B cells and severely blocks B cell development in mice. Deletion of Mettl14 impairs interleukin-7 (IL-7)-induced pro-B cell proliferation and the large-pre-B-to-small-pre-B transition and causes dramatic abnormalities in gene expression programs important for B cell development. Suppression of a group of transcripts by cytoplasmic m6A reader YTHDF2 is critical to the IL-7-induced pro-B cell proliferation. In contrast, the block in the large-pre-B-to-small-pre-B transition is independent of YTHDF1 or YTHDF2 but is associated with a failure to properly upregulate key transcription factors regulating this transition. Our data highlight the important regulatory roles of the RNA m6A methylation and its reader proteins in early B cell development.

Author Correction: Kethoxal-assisted single-stranded DNA sequencing captures global transcription dynamics and enhancer activity in situ.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SETD5-Coordinated Chromatin Reprogramming Regulates Adaptive Resistance to Targeted Pancreatic Cancer Therapy.

Molecular mechanisms underlying adaptive targeted therapy resistance in pancreatic ductal adenocarcinoma (PDAC) are poorly understood. Here, we identify SETD5 as a major driver of PDAC resistance to MEK1/2 inhibition (MEKi). SETD5 is induced by MEKi resistance and its deletion restores refractory PDAC vulnerability to MEKi therapy in mouse models and patient-derived xenografts. SETD5 lacks histone methyltransferase activity but scaffolds a co-repressor complex, including HDAC3 and G9a. Gene silencing by the SETD5 complex regulates known drug resistance pathways to reprogram cellular responses to MEKi. Pharmacological co-targeting of MEK1/2, HDAC3, and G9a sustains PDAC tumor growth inhibition in vivo. Our work uncovers SETD5 as a key mediator of acquired MEKi therapy resistance in PDAC and suggests a context for advancing MEKi use in the clinic.

Kethoxal-assisted single-stranded DNA sequencing captures global transcription dynamics and enhancer activity in situ.

Transcription is a highly dynamic process that generates single-stranded DNA (ssDNA) in the genome as 'transcription bubbles'. Here we describe a kethoxal-assisted single-stranded DNA sequencing (KAS-seq) approach, based on the fast and specific reaction between N3-kethoxal and guanines in ssDNA. KAS-seq allows rapid (within 5 min), sensitive and genome-wide capture and mapping of ssDNA produced by transcriptionally active RNA polymerases or other processes in situ using as few as 1,000 cells. KAS-seq enables definition of a group of enhancers that are single-stranded and enrich unique sequence motifs. These enhancers are associated with specific transcription-factor binding and exhibit more enhancer-promoter interactions than typical enhancers do. Under conditions that inhibit protein condensation, KAS-seq uncovers a rapid release of RNA polymerase II (Pol II) from a group of promoters. KAS-seq thus facilitates fast and accurate analysis of transcription dynamics and enhancer activities simultaneously in both low-input and high-throughput manner.

Refined spatial temporal epigenomic profiling reveals intrinsic connection between PRDM9-mediated H3K4me3 and the fate of double-stranded breaks.

Meiotic recombination is initiated by the formation of double-strand breaks (DSBs), which are repaired as either crossovers (COs) or noncrossovers (NCOs). In most mammals, PRDM9-mediated H3K4me3 controls the nonrandom distribution of DSBs; however, both the timing and mechanism of DSB fate control remain largely undetermined. Here, we generated comprehensive epigenomic profiles of synchronized mouse spermatogenic cells during meiotic prophase I, revealing spatiotemporal and functional relationships between epigenetic factors and meiotic recombination. We find that PRDM9-mediated H3K4me3 at DSB hotspots, coinciding with H3K27ac and H3K36me3, is intimately connected with the fate of the DSB. Our data suggest that the fate decision is likely made at the time of DSB formation: earlier formed DSBs occupy more open chromatins and are much more competent to proceed to a CO fate. Our work highlights an intrinsic connection between PRDM9-mediated H3K4me3 and the fate decision of DSBs, and provides new insight into the control of CO homeostasis.

Regulation of Co-transcriptional Pre-mRNA Splicing by m6A through the Low-Complexity Protein hnRNPG.

N6-methyladenosine (m6A) modification occurs co-transcriptionally and impacts pre-mRNA processing; however, the mechanism of co-transcriptional m6A-dependent alternative splicing regulation is still poorly understood. Heterogeneous nuclear ribonucleoprotein G (hnRNPG) is an m6A reader protein that binds RNA through RRM and Arg-Gly-Gly (RGG) motifs. Here, we show that hnRNPG directly binds to the phosphorylated carboxy-terminal domain (CTD) of RNA polymerase II (RNAPII) using RGG motifs in its low-complexity region. Through interactions with the phosphorylated CTD and nascent RNA, hnRNPG associates co-transcriptionally with RNAPII and regulates alternative splicing transcriptome-wide. m6A near splice sites in nascent pre-mRNA modulates hnRNPG binding, which influences RNAPII occupancy patterns and promotes exon inclusion. Our results reveal an integrated mechanism of co-transcriptional m6A-mediated splicing regulation, in which an m6A reader protein uses RGG motifs to co-transcriptionally interact with both RNAPII and m6A-modified nascent pre-mRNA to modulate RNAPII occupancy and alternative splicing.