Pseudouridine synthase 1 promotes hepatocellular carcinoma through mRNA pseudouridylation to enhance the translation of oncogenic mRNAs.

BACKGROUND AND AIMS: Pseudouridine is a prevalent RNA modification and is highly present in the serum and urine of patients with HCC. However, the role of pseudouridylation and its modifiers in HCC remains unknown. We investigated the function and underlying mechanism of pseudouridine synthase 1 (PUS1) in HCC. APPROACH AND RESULTS: By analyzing the TCGA data set, PUS1 was found to be significantly upregulated in human HCC specimens and positively correlated with tumor grade and poor prognosis of HCC. Knockdown of PUS1 inhibited cell proliferation and the growth of tumors in a subcutaneous xenograft mouse model. Accordingly, increased cell proliferation and tumor growth were observed in PUS1-overexpressing cells. Furthermore, overexpression of PUS1 significantly accelerates tumor formation in a mouse HCC model established by hydrodynamic tail vein injection, while knockout of PUS1 decreases it. Additionally, PUS1 catalytic activity is required for HCC tumorigenesis. Mechanistically, we profiled the mRNA targets of PUS1 by utilizing surveying targets by apolipoprotein B mRNA-editing enzyme 1 (APOBEC1)-mediated profiling and found that PUS1 incorporated pseudouridine into mRNAs of a set of oncogenes, thereby endowing them with greater translation capacity. CONCLUSIONS: Our study highlights the critical role of PUS1 and pseudouridylation in HCC development, and provides new insight that PUS1 enhances the protein levels of a set of oncogenes, including insulin receptor substrate 1 (IRS1) and c-MYC, by means of pseudouridylation-mediated mRNA translation.

RNA-binding protein CCDC137 activates AKT signaling and promotes hepatocellular carcinoma through a novel non-canonical role of DGCR8 in mRNA localization.

BACKGROUND: RNA binding proteins (RBPs)-regulated gene expression play a vital role in various pathological processes, including the progression of cancer. However, the role of RBP in hepatocellular carcinoma (HCC) remains much unknown. In this study, we aimed to explore the contribution of RBP CCDC137 in HCC development. METHODS: We analyzed the altered expression level and clinical significance of CCDC137 in database and HCC specimens. In vitro cell assays and in vivo spontaneous mouse models were used to assess the function of CCDC137. Finally, the molecular mechanisms of how CCDC137 regulates gene expression and promotes HCC was explored. RESULTS: CCDC137 is aberrantly upregulated in HCC and correlates with poor clinical outcomes in HCC patients. CCDC137 markedly promoted HCC proliferation and progression in vitro and in vivo. Mechanistically, CCDC137 binds with FOXM1, JTV1, LASP1 and FLOT2 mRNAs, which was revealed by APOBEC1-mediated profiling, to increase their cytoplasmic localization and thus enhance their protein expressions. Upregulation of FOXM1, JTV1, LASP1 and FLOT2 subsequently synergistically activate AKT signaling and promote HCC. Interestingly, we found that CCDC137 binds with the microprocessor protein DGCR8 and DGCR8 has a novel non-canonical function in mRNA subcellular localization, which mediates the cytoplasmic distribution of mRNAs regulated by CCDC137. CONCLUSIONS: Our results identify a critical proliferation-related role of CCDC137 and reveal a novel CCDC137/DGCR8/mRNA localization/AKT axis in HCC progression, which provide a potential target for HCC therapy.

The Conserved LncRNA DIO3OS Restricts Hepatocellular Carcinoma Stemness by Interfering with NONO-Mediated Nuclear Export of ZEB1 mRNA.

Hepatocellular carcinoma (HCC) is an aggressive and fatal disease caused by a subset of cancer stem cells (CSCs). It is estimated that there are approximately 100 000 long noncoding RNAs (lncRNAs) in humans. However, the mechanisms by which lncRNAs affect tumor stemness remain poorly understood. In the present study, it is found that DIO3OS is a conserved lncRNA that is generally downregulated in multiple cancers, including HCC, and its low expression correlates with poor clinical outcomes in HCC. In in vitro cancer cell lines and an in vivo spontaneous HCC mouse model, DIO3OS markedly represses tumor development via its suppressive role in CSCs through downregulation of zinc finger E-box binding homeobox 1 (ZEB1). Interestingly, DIO3OS represses ZEB1 post-transcriptionally without affecting its mRNA levels. Subsequent experiments show that DIO3OS interacts with the NONO protein and restricts NONO-mediated nuclear export of ZEB1 mRNA. Overall, these findings demonstrate that the DIO3OS-NONO-ZEB1 axis restricts HCC development and offers a valuable candidate for CSC-targeted therapeutics for HCC.

Characterization of Long Non-coding RNAs Modified by m6A RNA Methylation in Skeletal Myogenesis.

Proper development of mammalian skeletal muscle relies on precise gene expression regulation. Our previous studies revealed that muscle development is regulated by both mRNA and long non-coding RNAs (lncRNAs). Accumulating evidence has demonstrated that N6-methyladenosine (m6A) plays important roles in various biological processes, making it essential to profile m6A modification on a transcriptome-wide scale in developing muscle. Patterns of m6A methylation in lncRNAs in developing muscle have not been uncovered. Here, we reveal differentially expressed lncRNAs and report temporal m6A methylation patterns in lncRNAs expressed in mouse myoblasts and myotubes by RNA-seq and methylated RNA immunoprecipitation (MeRIP) sequencing. Many lncRNAs exhibit temporal differential expression, and m6A-lncRNAs harbor the consensus m6A motif "DRACH" along lncRNA transcripts. Interestingly, we found that m6A methylation levels of lncRNAs are positively correlated with the transcript abundance of lncRNAs. Overexpression or knockdown of m6A methyltransferase METTL3 alters the expression levels of these lncRNAs. Furthermore, we highlight that the function of m6A genic lncRNAs might correlate to their nearby mRNAs. Our work reveals a fundamental expression reference of m6A-mediated epitranscriptomic modifications in lncRNAs that are temporally expressed in developing muscle.

Dynamic m6A mRNA Methylation Reveals the Role of METTL3/14-m6A-MNK2-ERK Signaling Axis in Skeletal Muscle Differentiation and Regeneration.

N6-methyladenosine (m6A) RNA methylation has emerged as an important factor in various biological processes by regulating gene expression. However, the dynamic profile, function and underlying molecular mechanism of m6A modification during skeletal myogenesis remain elusive. Here, we report that members of the m6A core methyltransferase complex, METTL3 and METTL14, are downregulated during skeletal muscle development. Overexpression of either METTL3 or METTL14 dramatically blocks myotubes formation. Correspondingly, knockdown of METTL3 or METTL14 accelerates the differentiation of skeletal muscle cells. Genome-wide transcriptome analysis suggests ERK/MAPK is the downstream signaling pathway that is regulated to the greatest extent by METTL3/METTL14. Indeed, METTL3/METTL14 expression facilitates ERK/MAPK signaling. Via MeRIP-seq, we found that MNK2, a critical regulator of ERK/MAPK signaling, is m6A modified and is a direct target of METTL3/METTL14. We further revealed that YTHDF1 is a potential reader of m6A on MNK2, regulating MNK2 protein levels without affecting mRNA levels. Furthermore, we discovered that METTL3/14-MNK2 axis was up-regulated notably after acute skeletal muscle injury. Collectively, our studies revealed that the m6A writers METTL3/METTL14 and the m6A reader YTHDF1 orchestrate MNK2 expression posttranscriptionally and thus control ERK signaling, which is required for the maintenance of muscle myogenesis and may contribute to regeneration.

mascRNA and its parent lncRNA MALAT1 promote proliferation and metastasis of hepatocellular carcinoma cells by activating ERK/MAPK signaling pathway.

MALAT1-associated small cytoplasmic RNA (mascRNA) is a cytoplasmic tRNA-like small RNA derived from nucleus-located long noncoding RNA (lncRNA) metastasis-associated lung adenocarcinoma transcript 1 (MALAT1). While MALAT1 was extensively studied and was found to function in multiple cellular processes, including tumorigenesis and tumor progression, the role of mascRNA was largely unknown. Here we show that mascRNA is upregulated in multiple cancer cell lines and hepatocellular carcinoma (HCC) clinical samples. Using HCC cells as model, we found that mascRNA and its parent lncRNA MALAT1 can both promote cell proliferation, migration, and invasion in vitro. Correspondingly, both of them can enhance the tumor growth in mice subcutaneous tumor model and can promote metastasis by tail intravenous injection of HCC cells. Furthermore, we revealed that mascRNA and MALAT1 can both activate ERK/MAPK signaling pathway, which regulates metastasis-related genes and may contribute to the aggressive phenotype of HCC cells. Our results indicate a coordination in function and mechanism of mascRNA and MALAT1 during development and progress of HCC, and provide a paradigm for deciphering tRNA-like structures and their parent transcripts in mammalian cells.

METTL3 regulates skeletal muscle specific miRNAs at both transcriptional and post-transcriptional levels.

METTL3 increasing the mature miRNA levels via N6-Methyladenosine (m6A) modification of primary miRNA (pri-miRNA) transcripts has emerged as an important post-transcriptional regulation of miRNA biogenesis. Our previous studies and others have showed that muscle specific miRNAs are essential for skeletal muscle differentiation. Whether these miRNAs are also regulated by METTL3 is still unclear. Here, we found that m6A motifs were present around most of these miRNAs, which were indeed m6A modified as confirmed by m6A-modified RNA immunoprecipitation (m6A RIP). However, we surprisingly found that these muscle specific miRNAs were repressed instead of increased by METTL3 in C2C12 in vitro differentiation and mouse skeletal muscle regeneration after injury in vivo model. To elucidate the underlined mechanism, we performed reporter assays in 293T cells and validated METTL3 increasing these miRNAs at post-transcriptional level as expected. Furthermore, in myogenic C2C12 cells, we found that METTL3 not only repressed the expression of myogenic transcription factors (TFs) which can enhance the muscle specific miRNAs, but also increased the expression of epigenetic regulators which can repress these miRNAs. Thus, METTL3 could repress the muscle specific miRNAs at transcriptional level indirectly. Taken together, our results demonstrated that skeletal muscle specific miRNAs were repressed by METTL3 and such repression is likely synthesized transcriptional and post-transcriptional regulations.

N6-methyladenine demethylase ALKBH1 inhibits the differentiation of skeletal muscle.

DNA N6-methyladenine (N6-mA) was recently recognized as a new epigenetic modification in mammalian genome, and ALKBH1 was discovered as its demethylase. Knock-out mice studies revealed that ALKBH1 was indispensable for normal embryonic development. However, the function of ALKBH1 in myogenesis is largely unknown. In this study, we found that N6-mA showed a steady increase, going along with a strong decrease of ALKBH1 during skeletal muscle development. Our results also showed that ALKBH1 enhanced proliferation and inhibited differentiation of C2C12 cells. Genome-wide transcriptome analysis and reporter assays further revealed that ALKBH1 accomplished the differentiation inhibiting function by regulating a core set of genes and multiple signaling pathways, including increasing chemokine (C-X-C motif) ligand 14 (CXCL14) and activating ERK signaling. Taken together, our results demonstrated that ALKBH1 is critical for the myogenic differentiation of C2C12 cells, and suggested that N6-mA might be a new epigenetic mechanism for the regulation of myogenesis.

The landscape of long noncoding RNA-involved and tumor-specific fusions across various cancers.

The majority of the human genome encodes long noncoding RNA (lncRNA) genes, critical regulators of various cellular processes, which largely outnumber protein-coding genes. However, lncRNA-involved fusions have not been surveyed and characterized yet. Here, we present a systematic study of the lncRNA fusion landscape across cancer types and identify >30 000 high-confidence tumor-specific lncRNA fusions (using 8284 tumor and 6946 normal samples). Fusions positively correlated with DNA damage and cancer stemness and were specifically low in microsatellite instable (MSI)-High or virus-infected tumors. Moreover, fusions distribute differently among cancer molecular subtypes, but with shared enrichment in tumors that are microsatellite stable (MSS), with high somatic copy number alterations (SCNA), and with poor survival. Importantly, we find a potentially new mechanism, mediated by enhancer RNAs (eRNA), which generates secondary fusions that form densely connected fusion networks with many fusion hubs targeted by FDA-approved drugs. Finally, we experimentally validate functions of two tumor-promoting chimeric proteins derived from mRNA-lncRNA fusions, KDM4B-G039927 and EPS15L1-lncOR7C2-1. The EPS15L1 fusion protein may regulate (Gasdermin E) GSDME, critical in pyroptosis and anti-tumor immunity. Our study completes the fusion landscape in cancers, sheds light on fusion mechanisms, and enriches lncRNA functions in tumorigenesis and cancer progression.

RNA sequencing: new technologies and applications in cancer research.

Over the past few decades, RNA sequencing has significantly progressed, becoming a paramount approach for transcriptome profiling. The revolution from bulk RNA sequencing to single-molecular, single-cell and spatial transcriptome approaches has enabled increasingly accurate, individual cell resolution incorporated with spatial information. Cancer, a major malignant and heterogeneous lethal disease, remains an enormous challenge in medical research and clinical treatment. As a vital tool, RNA sequencing has been utilized in many aspects of cancer research and therapy, including biomarker discovery and characterization of cancer heterogeneity and evolution, drug resistance, cancer immune microenvironment and immunotherapy, cancer neoantigens and so on. In this review, the latest studies on RNA sequencing technology and their applications in cancer are summarized, and future challenges and opportunities for RNA sequencing technology in cancer applications are discussed.

Schisandrin B attenuates renal fibrosis via miR-30e-mediated inhibition of EMT.

Tubulointerstitial fibrosis (TIF) is the main pathologic feature of end-stage renal disease. Epithelial-mesenchymal transition (EMT) of proximal tubular cells (PTCs) is one of the most significant features of TIF. MicroRNAs play critical roles during EMT in TIF. However, whether miRNAs can be used as therapeutic targets in TIF therapy remains undetermined. We found that miR-30e, a member of the miR-30 family, is deregulated in TGF-ß1-induced PTCs, TIF mice and human fibrotic kidney tissues. Moreover, transcription factors that induce EMT, such as snail, slug, and Zeb2, were direct targets of miR-30e. Using a cell-based miR-30e promoter luciferase reporter system, Schisandrin B (Sch B) was selected for the enhancement of miR-30e transcriptional activity. Our results indicate that Sch B can decrease the expression of snail, slug, and Zeb2, thereby attenuating the EMT of PTCs during TIF by upregulating miR-30e, both in vivo and in vitro. This study shows that miR-30e can serve as a therapeutic target in the treatment of patients with TIF and that Sch B may potentially be used in therapy against renal fibrosis.

Delineation of the role of chromatin assembly and the Rtt101Mms1 E3 ubiquitin ligase in DNA damage checkpoint recovery in budding yeast.

The DNA damage checkpoint is activated in response to DNA double-strand breaks (DSBs). We had previously shown that chromatin assembly mediated by the histone chaperone Asf1 triggers inactivation of the DNA damage checkpoint in yeast after DSB repair, also called checkpoint recovery. Here we show that chromatin assembly factor 1 (CAF-1) also contributes to chromatin reassembly after DSB repair, explaining its role in checkpoint recovery. Towards understanding how chromatin assembly promotes checkpoint recovery, we find persistent presence of the damage sensors Ddc1 and Ddc2 after DSB repair in asf1 mutants. The genes encoding the E3 ubiquitin ligase complex Rtt101Mms1 are epistatic to ASF1 for survival following induction of a DSB, and Rtt101Mms1 are required for checkpoint recovery after DSB repair but not for chromatin assembly. By contrast, the Mms22 substrate adaptor that is degraded by Rtt101Mms1 is required for DSB repair per se. Deletion of MMS22 blocks loading of Rad51 at the DSB, while deletion of ASF1 or RTT101 leads to persistent Rad51 loading. We propose that checkpoint recovery is promoted by Rtt101Mms1-mediated ubiquitylation of Mms22 in order to halt Mms22-dependent loading of Rad51 onto double-stranded DNA after DSB repair, in concert with the chromatin assembly-mediated displacement of Rad51 and checkpoint sensors from the site of repair.

Conservation and divergence of transcriptional coregulations between box C/D snoRNA and ribosomal protein genes in Ascomycota.

Coordinated assembly of the ribosome is essential for proper translational activity in eukaryotic cells. It is therefore critical to coordinate the expression of components of ribosomal programs with the cell's nutritional status. However, coordinating expression of these components is poorly understood. Here, by combining experimental and computational approaches, we systematically identified box C/D snoRNAs in four fission yeasts and found that the expression of box C/D snoRNA and ribosomal protein (RP) genes were orchestrated by a common Homol-D box, thereby ensuring a constant balance of these two genetic components. Interestingly, such transcriptional coregulations could be observed in most Ascomycota species and were mediated by different cis-regulatory elements. Via the reservation of cis elements, changes in spatial configuration, the substitution of cis elements, and gain or loss of cis elements, the regulatory networks of box C/D snoRNAs evolved to correspond with those of the RP genes, maintaining transcriptional coregulation between box C/D snoRNAs and RP genes. Our results indicate that coregulation via common cis elements is an important mechanism to coordinate expression of the RP and snoRNA genes, which ensures a constant balance of these two components.

The ribosomal protein rpl26 promoter is required for its 3' sense terminus ncRNA transcription in Schizosaccharomyces pombe, implicating a new transcriptional mechanism for ncRNAs.

Transcriptome studies have revealed that many non-coding RNAs (ncRNAs) are located near the 3' sense terminus of protein-coding genes. However, the transcription and function of these RNAs remain elusive. Here, we identify a 3' sense termini-associated sRNA (TASR) downstream of rpl26 in Schizosaccharomyces pombe (S. pombe). Structure and function assays indicate that the TASR is an H/ACA box snoRNA required for 18S rRNA pseudouridylation at U121 and U305 sites and is therefore a cognate of snR49 from the budding yeast. Transcriptional studies show that pre-snR49 overlaps most of the coding sequence (CDS) of rpl26. Using scanning deletion analysis within promoter region, we show that the rpl26 promoter is required for the 3' TASR transcription. Interestingly, chromosomal synteny of rpl26-snR49 is found in the Schizosaccharomyces groups. Taken together, we have revealed a new transcriptional mechanism for 3' sense TASRs, which are transcribed by the same promoter as their upstream protein genes. These results further suggest that the origin and function of 3' sense ncRNAs are associated with upstream genes in higher eukaryotes.

Deciphering the transcriptional regulation of microRNA genes in humans with ACTLocater.

Understanding the transcriptional regulation of microRNAs (miRNAs) is extremely important for determining the specific roles they play in signaling cascades. However, precise identification of transcription factor binding sites (TFBSs) orchestrating the expressions of miRNAs remains a challenge. By combining accessible chromatin sequences of 12 cell types released by the ENCODE Project, we found that a significant fraction (~80%) of such integrated sequences, evolutionary conserved and in regions upstream of human miRNA genes that are independently transcribed, were preserved across cell types. Accordingly, we developed a computational method, Accessible and Conserved TFBSs Locater (ACTLocater), incorporating this chromatin feature and evolutionary conservation to identify the TFBSs associated with human miRNA genes. ACTLocater achieved high positive predictive values, as revealed by the experimental validation of FOXA1 predictions and by the comparison of its predictions of some other transcription factors (TFs) to empirical ChIP-seq data. Most notably, ACTLocater was widely applicable as indicated by the successful prediction of TF → miRNA interactions in cell types whose chromatin accessibility profiles were not incorporated. By applying ACTLocater to TFs with characterized binding specificities, we compiled a novel repository of putative TF → miRNA interactions and displayed it in ACTViewer, providing a promising foundation for future investigations to elucidate the regulatory mechanisms of miRNA transcription in humans.

RTL-P: a sensitive approach for detecting sites of 2'-O-methylation in RNA molecules.

2'-O-methylation is present within various cellular RNAs and is essential to RNA biogenesis and functionality. Several methods have been developed for the identification and localization of 2'-O-methylated sites in RNAs; however, the detection of RNA modifications, especially in low-abundance RNAs and small non-coding RNAs with a 2'-O-methylation at the 3'-end, remains a difficult task. Here, we introduce a new method to detect 2'-O-methylated sites in diverse RNA species, referred to as RTL-P [Reverse Transcription at Low deoxy-ribonucleoside triphosphate (dNTP) concentrations followed by polymerase chain reaction (PCR)] that demonstrates precise mapping and superior sensitivity compared with previous techniques. The main procedures of RTL-P include a site-specific primer extension by reverse transcriptase at a low dNTP concentration and a semi-quantitative PCR amplification step. No radiolabeled or fluorescent primers are required. By designing specific RT primers, we used RTL-P to detect both previously identified and novel 2'-O-methylated sites in human and yeast ribosomal RNAs (rRNAs), as well as mouse piwi-interacting RNAs (piRNAs). These results demonstrate the powerful application of RTL-P for the systematic analysis of fully or partially methylated residues in diverse RNA species, including low-abundance RNAs or small non-coding RNAs such as piRNAs and microRNAs (miRNAs).