Pseudouridine Detection and Quantification Using Bisulfite Incorporation Hindered Ligation.

Pseudouridine (Ψ) is a widespread RNA modification found in various RNA species, including rRNA, tRNA, snRNA, mRNA, and long noncoding RNA (lncRNA). Understanding the function of Ψ in these RNA types requires a robust method for the detection and quantification of the Ψ level at single-nucleotide resolution. A previously used method utilizes Ψ labeling by N-cyclohexyl-N'-ß-(4-methylmorpholinium)ethylcarbodiimide (CMC). The quantification of Ψ is based on the stop ratio after reverse transcription. However, the use of CMC followed by strong alkaline treatment causes severe RNA degradation, often requiring a large amount of RNA. The removal of CMC and recovery of RNA by ethanol precipitation are also time-consuming. Here, we introduce a Bisulfite Incorporation Hindered ligation-based method (BIHIND), which can detect and quantify Ψ sites on rRNA, mRNA, and noncoding RNA. BIHIND can be coupled with quantitative PCR (BIHIND-qPCR) for quantitative detection of Ψ fraction at individual modification sites, as well as with next-generation sequencing (BIHIND-seq) for high-throughput sequencing of Ψ without requiring reverse transcription. We validated the robustness of BIHIND with the elucidation of Ψ dynamics following pseudouridine synthase depletion.

Small-molecule inhibition of the METTL3/METTL14 complex suppresses neuroblastoma tumor growth and promotes differentiation.

The N6-methyladenosine (m6A) RNA modification is an important regulator of gene expression. m6A is deposited by a methyltransferase complex that includes methyltransferase-like 3 (METTL3) and methyltransferase-like 14 (METTL14). High levels of METTL3/METTL14 drive the growth of many types of adult cancer, and METTL3/METTL14 inhibitors are emerging as new anticancer agents. However, little is known about the m6A epitranscriptome or the role of the METTL3/METTL14 complex in neuroblastoma, a common pediatric cancer. Here, we show that METTL3 knockdown or pharmacologic inhibition with the small molecule STM2457 leads to reduced neuroblastoma cell proliferation and increased differentiation. These changes in neuroblastoma phenotype are associated with decreased m6A deposition on transcripts involved in nervous system development and neuronal differentiation, with increased stability of target mRNAs. In preclinical studies, STM2457 treatment suppresses the growth of neuroblastoma tumors in vivo. Together, these results support the potential of METTL3/METTL14 complex inhibition as a therapeutic strategy against neuroblastoma.

Ultrafast bisulfite sequencing detection of 5-methylcytosine in DNA and RNA.

Bisulfite sequencing (BS-seq) to detect 5-methylcytosine (5mC) is limited by lengthy reaction times, severe DNA damage, overestimation of 5mC level and incomplete C-to-U conversion of certain DNA sequences. We present ultrafast BS-seq (UBS-seq), which uses highly concentrated bisulfite reagents and high reaction temperatures to accelerate the bisulfite reaction by ~13-fold, resulting in reduced DNA damage and lower background noise. UBS-seq allows library construction from small amounts of purified genomic DNA, such as from cell-free DNA or directly from 1 to 100 mouse embryonic stem cells, with less overestimation of 5mC level and higher genome coverage than conventional BS-seq. Additionally, UBS-seq quantitatively maps RNA 5-methylcytosine (m5C) from low inputs of mRNA and allows the detection of m5C stoichiometry in highly structured RNA sequences. Our UBS-seq results identify NSUN2 as the major 'writer' protein responsible for the deposition of ~90% of m5C sites in HeLa mRNA and reveal enriched m5C sites in 5'-regions of mammalian mRNA, which may have functional roles in mRNA translation regulation.

BID-seq for transcriptome-wide quantitative sequencing of mRNA pseudouridine at base resolution.

Pseudouridine (Ψ) is an abundant RNA modification that is present in and affects the functions of diverse non-coding RNA species, including rRNA, tRNA and small nuclear RNA. Ψ also exists in mammalian mRNA and probably exhibits functional roles; however, functional investigations of mRNA Ψ modifications in mammals have been hampered by the lack of a quantitative method that detects Ψ at base precision. We have recently developed bisulfite-induced deletion sequencing (BID-seq), which provides the community with a quantitative method to map RNA Ψ distribution transcriptome-wide at single-base resolution. Here, we describe an optimized BID-seq protocol for mapping Ψ distribution across cellular mRNAs, which includes fast steps in both library preparation and data analysis. This protocol generates highly reproducible results by inducing high deletion ratios at Ψ modification within diverse sequence contexts, and meanwhile displayed almost zero background deletions at unmodified uridines. When used for transcriptome-wide Ψ profiling in mouse embryonic stem cells, the current protocol uncovered 8,407 Ψ sites from as little as 10 ng of polyA+ RNA input. This optimized BID-seq workflow takes 5 days to complete and includes four main sections: RNA preparation, library construction, next-generation sequencing (NGS) and data analysis. Library construction can be completed by researchers who have basic knowledge and skills in molecular biology and genetics. In addition to the experimental protocol, we provide BID-pipe ( https://github.com/y9c/pseudoU-BIDseq ), a user-friendly data analysis pipeline for Ψ site detection and modification stoichiometry quantification, requiring only basic bioinformatic and computational skills to uncover Ψ signatures from BID-seq data.

Quantitative sequencing using BID-seq uncovers abundant pseudouridines in mammalian mRNA at base resolution.

Functional characterization of pseudouridine (Ψ) in mammalian mRNA has been hampered by the lack of a quantitative method that maps Ψ in the whole transcriptome. We report bisulfite-induced deletion sequencing (BID-seq), which uses a bisulfite-mediated reaction to convert pseudouridine stoichiometrically into deletion upon reverse transcription without cytosine deamination. BID-seq enables detection of abundant Ψ sites with stoichiometry information in several human cell lines and 12 different mouse tissues using 10-20 ng input RNA. We uncover consensus sequences for Ψ in mammalian mRNA and assign different 'writer' proteins to individual Ψ deposition. Our results reveal a transcript stabilization role of Ψ sites installed by TRUB1 in human cancer cells. We also detect the presence of Ψ within stop codons of mammalian mRNA and confirm the role of Ψ in promoting stop codon readthrough in vivo. BID-seq will enable future investigations of the roles of Ψ in diverse biological processes.

The m6A methyltransferase METTL3 regulates muscle maintenance and growth in mice.

Skeletal muscle serves fundamental roles in organismal health. Gene expression fluctuations are critical for muscle homeostasis and the response to environmental insults. Yet, little is known about post-transcriptional mechanisms regulating such fluctuations while impacting muscle proteome. Here we report genome-wide analysis of mRNA methyladenosine (m6A) dynamics of skeletal muscle hypertrophic growth following overload-induced stress. We show that increases in METTL3 (the m6A enzyme), and concomitantly m6A, control skeletal muscle size during hypertrophy; exogenous delivery of METTL3 induces skeletal muscle growth, even without external triggers. We also show that METTL3 represses activin type 2 A receptors (ACVR2A) synthesis, blunting activation of anti-hypertrophic signaling. Notably, myofiber-specific conditional genetic deletion of METTL3 caused spontaneous muscle wasting over time and abrogated overload-induced hypertrophy; a phenotype reverted by co-administration of a myostatin inhibitor. These studies identify a previously unrecognized post-transcriptional mechanism promoting the hypertrophic response of skeletal muscle via control of myostatin signaling.

METTL3 Regulates Liver Homeostasis, Hepatocyte Ploidy, and Circadian Rhythm-Controlled Gene Expression in Mice.

N6-methyladenosine (m6A), the most abundant internal modifier of mRNAs installed by the methyltransferase 13 (METTL3) at the (G/A)(m6A)C motif, plays a critical role in the regulation of gene expression. METTL3 is essential for embryonic development, and its dysregulation is linked to various diseases. However, the role of METTL3 in liver biology is largely unknown. In this study, METTL3 function was unraveled in mice depleted of Mettl3 in neonatal livers (Mettl3fl/fl; Alb-Cre). Liver-specific Mettl3 knockout (M3LKO) mice exhibited global decrease in m6A on polyadenylated RNAs and pathologic features associated with nonalcoholic fatty liver disease (eg, hepatocyte ballooning, ductular reaction, microsteatosis, pleomorphic nuclei, DNA damage, foci of altered hepatocytes, focal lobular and portal inflammation, and elevated serum alanine transaminase/alkaline phosphatase levels). Mettl3-depleted hepatocytes were highly proliferative, with decreased numbers of binucleate hepatocytes and increased nuclear polyploidy. M3LKO livers were characterized by reduced m6A and expression of several key metabolic transcripts regulated by circadian rhythm and decreased nuclear protein levels of the core clock transcription factors BMAL1 and CLOCK. A significant decrease in total Bmal1 and Clock mRNAs but an increase in their nuclear levels were observed in M3LKO livers, suggesting impaired nuclear export. Consistent with the phenotype, methylated (m6A) RNA immunoprecipitation coupled with sequencing and RNA sequencing revealed transcriptome-wide loss of m6A markers and alterations in abundance of mRNAs involved in metabolism in M3LKO. Collectively, METTL3 and m6A modifications are critical regulators of liver homeostasis and function.

Corrigendum: Loss of TARBP2 Drives the Progression of Hepatocellular Carcinoma via miR-145-SERPINE1 Axis.

[This corrects the article DOI: 10.3389/fonc.2021.620912.].

Loss of TARBP2 Drives the Progression of Hepatocellular Carcinoma via miR-145-SERPINE1 Axis.

The clinical outcomes of hepatocellular carcinoma (HCC) remain dismal. Elucidating the molecular mechanisms for the progression of aggressive HCC holds the promise for developing novel intervention strategies. The transactivation response element RNA-binding protein (TRBP/TARBP2), a key component of microRNA (miRNA) processing and maturation machinery has been shown to play conflicting roles in tumor development and progression. We sought to investigate the expression of TARBP2 in HCC using well-characterized HCC cell lines, patient-derived tissues and blood samples. Additionally, the potential prognostic and diagnostic value of TARBP2 in HCC were analyzed using Kaplan-Meier plots and ROC curve. Cell counting kit-8 (CCK-8), wound healing and transwell assays examined the ability of TARBP2 to induce cell proliferation, migration, and invasion in HCC cell lines. RNA sequencing was applied to identify the downstream elements of TARBP2. The interaction of potential targets of TARBP2, miR-145 and serpin family E member 1 (SERPINE1), was assessed using luciferase reporter assay. TARBP2 expression was down-regulated in HCC cell lines relative to normal hepatocyte cells, with a similar pattern further confirmed in tissue and blood samples. Notably, the loss of TARBP2 was demonstrated to promote proliferation, migration, and invasion in HCC cell lines. Interestingly, the reduction of TARBP2 was shown to result in the upregulation of SERPINE1, also known as plasminogen activator inhibitor (PAI-1), which is a vital gene of the HIF-1 signaling pathway. Knockdown of SERPINE1 rescued the TARBP2-lost phenotype. Moreover, TARBP2 depletion induced the upregulation of SERPINE1 through reducing the processing of miR-145, which directly targets SERPINE1. Finally, overexpression of miR-145 repressed SERPINE1 and rescued the functions in sh-TARBP2 HCC cells. Our findings underscore a linear TARBP2-miR-145-SERPINE1 pathway that drives HCC progression, with the potential as a novel intervention target for aggressive HCC.

E-cadherin Interacts With Posttranslationally-Modified AGO2 to Enhance miRISC Activity.

MicroRNAs (miRNAs) are small non-coding RNAs which post-transcriptionally suppress target mRNAs expression and/or translation to modulate pathophyological processes. Expression and function of miRNAs are fine-tuned by a conserved biogenesis machinery involves two RNase-dependent processing steps of miRNA maturation and the final step of miRNA-induced silencing complex (miRISC)-mediated target silencing. A functional miRISC requires Argonaute 2 (AGO2) as an essential catalytic component which plays central roles in miRISC function. We uncovered a post-translational regulatory mechanism of AGO2 by E-cadherin. Mechanistically, E-cadherin activates ERK to phosphorylate AGO2, along with enhanced protein glycosylation. Consequently, the phosphorylated AGO2 was stabilized and ultimately resulted in induced miRISC activity on gene silencing. This study revealed a novel pathway for miRNA regulation through an E-cadherin-mediated miRISC activation.

Stabilization of ERK-Phosphorylated METTL3 by USP5 Increases m6A Methylation.

N6-methyladenosine (m6A) is the most abundant mRNA modification and is installed by the METTL3-METTL14-WTAP methyltransferase complex. Although the importance of m6A methylation in mRNA metabolism has been well documented recently, regulation of the m6A machinery remains obscure. Through a genome-wide CRISPR screen, we identify the ERK pathway and USP5 as positive regulators of the m6A deposition. We find that ERK phosphorylates METTL3 at S43/S50/S525 and WTAP at S306/S341, followed by deubiquitination by USP5, resulting in stabilization of the m6A methyltransferase complex. Lack of METTL3/WTAP phosphorylation reduces decay of m6A-labeled pluripotent factor transcripts and traps mouse embryonic stem cells in the pluripotent state. The same phosphorylation can also be found in ERK-activated human cancer cells and contribute to tumorigenesis. Our study reveals an unrecognized function of ERK in regulating m6A methylation.

A New Model of Spontaneous Colitis in Mice Induced by Deletion of an RNA m6A Methyltransferase Component METTL14 in T Cells.

BACKGROUND AND AIMS: Mouse models of colitis have been used to study the pathogenesis of inflammatory bowel disease (IBD) and for pre-clinical development of therapeutic agents. Various epigenetic pathways have been shown to play important regulatory roles in IBD. Reversible N6-methyladenosine (m6A) methylation represents a new layer of post-transcriptional gene regulation that affects a variety of biological processes. We aim to study how deletion of a critical component of m6A writer complex, METTL14, in T cells affects the development of colitis. METHODS: Conditional Mettl14 was lineage specifically deleted with CD4-regulated Cre in T cells. Colitis phenotype was determined by H&E staining, colon weight-to-length ratio and cytokine expression. We additionally utilized T cell transfer model of colitis and adoptive transfer of regulatory T cells. Mice were treated with antibiotics to determine if the colitis could be attenuated. RESULTS: METTL14 deficiency in T cells induced spontaneous colitis in mice. This was characterized by increased inflammatory cell infiltration, increased colonic weight-to-length ratio and increased Th1 and Th17 cytokines. The colitis development was due to dysfunctional regulatory T (Treg) cells, as adoptive transfer of WT Treg cells attenuated the colitis phenotype. The METTL14-deficient Treg cells have decreased RORγt expression compared with WT controls. METTL14 deficiency caused impaired induction of naïve T cells into induced Treg cells. Antibiotic treatment notably attenuated the colitis development. CONCLUSION: Here we report a new mouse model of spontaneous colitis based on perturbation of RNA methylation in T cells. The colitis is T cell-mediated and dependent on the microbiome. This model represents a new tool for elucidating pathogenic pathways, studying the contribution of intestinal microbiome and preclinical testing of therapeutic agents for inflammatory bowel disease.

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.

YTHDF2 promotes mitotic entry and is regulated by cell cycle mediators.

The N6-methyladenosine (m6A) modification regulates mRNA stability and translation. Here, we show that transcriptomic m6A modification can be dynamic and the m6A reader protein YTH N6-methyladenosine RNA binding protein 2 (YTHDF2) promotes mRNA decay during cell cycle. Depletion of YTHDF2 in HeLa cells leads to the delay of mitotic entry due to overaccumulation of negative regulators of cell cycle such as Wee1-like protein kinase (WEE1). We demonstrate that WEE1 transcripts contain m6A modification, which promotes their decay through YTHDF2. Moreover, we found that YTHDF2 protein stability is dependent on cyclin-dependent kinase 1 (CDK1) activity. Thus, CDK1, YTHDF2, and WEE1 form a feedforward regulatory loop to promote mitotic entry. We further identified Cullin 1 (CUL1), Cullin 4A (CUL4A), damaged DNA-binding protein 1 (DDB1), and S-phase kinase-associated protein 2 (SKP2) as components of E3 ubiquitin ligase complexes that mediate YTHDF2 proteolysis. Our study provides insights into how cell cycle mediators modulate transcriptomic m6A modification, which in turn regulates the cell cycle.

N6-Deoxyadenosine Methylation in Mammalian Mitochondrial DNA.

N6-Methyldeoxyadenosine (6mA) has recently been shown to exist and play regulatory roles in eukaryotic genomic DNA (gDNA). However, the biological functions of 6mA in mammals have yet to be adequately explored, largely due to its low abundance in most mammalian genomes. Here, we report that mammalian mitochondrial DNA (mtDNA) is enriched for 6mA. The level of 6mA in HepG2 mtDNA is at least 1,300-fold higher than that in gDNA under normal growth conditions, corresponding to approximately four 6mA modifications on each mtDNA molecule. METTL4, a putative mammalian methyltransferase, can mediate mtDNA 6mA methylation, which contributes to attenuated mtDNA transcription and a reduced mtDNA copy number. Mechanistically, the presence of 6mA could repress DNA binding and bending by mitochondrial transcription factor (TFAM). Under hypoxia, the 6mA level in mtDNA could be further elevated, suggesting regulatory roles for 6mA in mitochondrial stress response. Our study reveals DNA 6mA as a regulatory mark in mammalian mtDNA.

Transcriptome-wide Mapping of Internal N7-Methylguanosine Methylome in Mammalian mRNA.

N7-methylguanosine (m7G) is a positively charged, essential modification at the 5' cap of eukaryotic mRNA, regulating mRNA export, translation, and splicing. m7G also occurs internally within tRNA and rRNA, but its existence and distribution within eukaryotic mRNA remain to be investigated. Here, we show the presence of internal m7G sites within mammalian mRNA. We then performed transcriptome-wide profiling of internal m7G methylome using m7G-MeRIP sequencing (MeRIP-seq). To map this modification at base resolution, we developed a chemical-assisted sequencing approach that selectively converts internal m7G sites into abasic sites, inducing misincorporation at these sites during reverse transcription. This base-resolution m7G-seq enabled transcriptome-wide mapping of m7G in human tRNA and mRNA, revealing distribution features of the internal m7G methylome in human cells. We also identified METTL1 as a methyltransferase that installs a subset of m7G within mRNA and showed that internal m7G methylation could affect mRNA translation.

Link Between m6A Modification and Cancers.

N6-methyladenosine (m6A) epitranscriptional modification has recently gained much attention. Through the development of m6A sequencing, the molecular mechanism and importance of m6A have been revealed. m6A is the most abundant internal modification in higher eukaryotic mRNAs, which plays crucial roles in mRNA metabolism and multiple biological processes. In this review, we introduce the characteristics of m6A regulators, including "writers" that create m6A mark, "erasers" that show demethylation activity and "readers" that decode m6A modification to govern the fate of modified transcripts. Moreover, we highlight the roles of m6A modification in several common cancers, including solid and non-solid tumors. The regulators of m6A exert enormous functions in cancer development, such as proliferation, migration and invasion. Especially, with the underlying mechanisms being uncovered, m6A and its regulators are expected to be the targets for the diagnosis and treatment of cancers.

Hepatitis B Virus Covalently Closed Circular DNA-Selective Droplet Digital PCR: A Sensitive and Noninvasive Method for Hepatocellular Carcinoma Diagnosis?

This commentary highlights the article by Huang et al that reports a highly sensitive assay for detection of closed circular DNA of hepatitis B virus.

Pin1 impairs microRNA biogenesis by mediating conformation change of XPO5 in hepatocellular carcinoma.

MicroRNA (miRNA) dysregulation is associated with the tumorigenesis and development of numerous human cancers. The defect in miRNA biogenesis is the main cause of miRNA dysregulation. We previously demonstrated that ERK-induced phosphorylation of XPO5 followed by peptidyl-prolyl cis/trans isomerase Pin1-mediated isomerization downregulates miRNA expression and contributes to hepatocellular carcinoma (HCC) development. However, how Pin1 precisely regulates miRNA biogenesis in HCC remains elusive. Here we reveal that Pin1 has a pivotal role in the miRNA maturation process by modulating phosphorylated Serine-Proline (pS-P) motif of XPO5 in a phosphorylation-dependent manner. By recognizing and binding to XPO5 via its WW domain, Pin1 catalyzes the conformation change of XPO5 and diminishes XPO5 ability to export pre-miRNAs from the nucleus to the cytoplasm, resulting in the reduced mature miRNA levels and promoted HCC development. Furthermore, downregulation of Pin1 by shRNA restores XPO5-dependent pre-miRNA export and effective biogenesis of mature miRNAs, leading to both in vitro and in vivo HCC inhibition. Therefore, our research discloses a new posttranscriptional regulatory mechanism of miRNA biosynthesis and provides the experimental basis for a novel HCC therapy by targeting Pin1.

HIF-1α promotes autophagic proteolysis of Dicer and enhances tumor metastasis.

HIF-1α, one of the most extensively studied oncogenes, is activated by a variety of microenvironmental factors. The resulting biological effects are thought to depend on its transcriptional activity. The RNAse enzyme Dicer is frequently downregulated in human cancers, which has been functionally linked to enhanced metastatic properties; however, current knowledge of the upstream mechanisms regulating Dicer is limited. In the present study, we identified Dicer as a HIF-1α-interacting protein in multiple types of cancer cell lines and different human tumors. HIF-1α downregulated Dicer expression by facilitating its ubiquitination by the E3 ligase Parkin, thereby enhancing autophagy-mediated degradation of Dicer, which further suppressed the maturation of known tumor suppressors, such as the microRNA let-7 and microRNA-200b. Consequently, expression of HIF-1α facilitated epithelial-mesenchymal transition (EMT) and metastasis in tumor-bearing mice. Thus, this study uncovered a connection between oncogenic HIF-1α and the tumor-suppressive Dicer. This function of HIF-1α is transcription independent and occurs through previously unrecognized protein interaction-mediated ubiquitination and autophagic proteolysis.