Crosstalk between RNA m6A Modification and Non-coding RNA Contributes to Cancer Growth and Progression

m6A methyltransferase METTL3-mediated lncRNA FOXD2-AS1 promotes the tumorigenesis of cervical cancer

Recent studies have indicated that long noncoding RNA (lncRNA) and N⁶-methyladenosine (m⁶A) methylation modification play critical roles in human cancers; however, their regulation on cervical cancer is largely unclear. Here, our study tries to investigate the underlying mechanisms by which lncRNA FOXD2 adjacent opposite strand RNA 1 (FOXD2-AS1) modulates cervical cancer tumorigenesis. Results illuminated that FOXD2-AS1 expression was significantly upregulated in cervical cancer cells and tissue, which was closely correlated to the unfavorable prognosis. Functionally, gain and loss-of-function assays showed that FOXD2-AS1 promoted the migration and proliferation of cervical cancer cells. Besides, FOXD2-AS1 silencing repressed the tumor growth in vivo. Mechanistically, m⁶A methyltransferase methyltransferase-like 3 (METTL3) enhanced the stability of FOXD2-AS1 and maintained its expression. Moreover, FOXD2-AS1 recruited lysine-specific demethylase 1 (LSD1) to the promoter region of p21 to silence its transcription abundance. In conclusion, these findings support that METTL3/FOXD2-AS1 accelerates cervical cancer progression via a m⁶A-dependent modality, which may serve as a potential therapeutic target for cervical cancer.

Analysis of m6A-Related lncRNAs for Prognosis Value and Response to Immune Checkpoint Inhibitors Therapy in Hepatocellular Carcinoma

Introduction

N6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in the progression of hepatocellular carcinoma (HCC). However, how their interaction is involved in the prognostic value of HCC and immune checkpoint inhibitors (ICIs) therapy remains unclear.

Methods

The RNA sequencing and clinical data of HCC patients were collected from TCGA database. The prognostic m6A-related lncRNAs were screened out with Pearson correlation test, univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) Cox regression. Patients with HCC were classified into 2 subtypes by consensus clustering. Survival analyses were performed to assess the prognostic value of different clusters and risk models. Potential tumor correlated biological pathways correlated with different clusters were explored through gene set enrichment analysis. We also identified the relationship of the risk model and clusters with response to immune checkpoint inhibitors (ICIs) therapy and tumor microenvironment (TME). Furthermore, the prognostic value of the 9 m6A-related lncRNAs was validated in the external cohort. Finally, the role of SNHG4 was explored by silencing and overexpression of SNHG4 through conducting proliferation, migration and invasion experiments.

Results

Patients from 2 clusters and different risk groups based on m6A-related lncRNAs had significantly different clinicopathological characteristics and overall survival outcomes. Tumor-correlated biological pathways were found to be correlated with Cluster 2 through GSEA. Moreover, we found that patients from different clusters and risk groups expressed higher levels of immune checkpoint genes and had distinct TME and different responses for ICIs therapy. Prognostic value of this risk model was further confirmed in the external cohort. Finally, consistent with the discovery, SNHG4 played an oncogenic role in vitro.

Conclusion

Our study demonstrated that the 9 m6A-related lncRNA signature may serve as a novel predictor in the prognosis of HCC and optimize (ICIs) therapy. SNHG4 plays an oncogenic role in HCC.

The power and the promise of circRNAs for cancer precision medicine with functional diagnostics and prognostic prediction

Circular RNA (circRNA) is a large class of covalently closed circular RNA. As a member of competitive endogenous RNA (ceRNA), it participates in the regulation of circRNA-miRNA-mRNA network and plays an important role in the regulation of physiology and pathology. CircRNA is produced by the reverse splicing of exon, intron or both, forming exon or intron circRNA. Studies have shown that circRNA is a ubiquitous molecule, which exceeds the linear mRNA distributed in human cells. Because of its covalent closed-loop structure, circRNA is resistant to RNase R, which is more stable than linear mRNA; circRNA is highly conserved in different species. It was found that circRNA competitively adsorbs miRNA, as a miRNA sponge, to involve in the expression regulation of a variety of genes and plays an important role in tumor development, invasion, metastasis and other processes. These molecules offer new potential opportunities for therapeutic intervention and serve as biomarkers for diagnosis. In this paper, the origin, characteristics and functions of circRNA and its role in tumor development, invasion and metastasis, diagnosis and prognosis are reviewed.

N6-Methylandenosine-Related lncRNA Signature Is a Novel Biomarkers of Prognosis and Immune Response in Colon Adenocarcinoma Patients

Background

Colon adenocarcinoma (COAD) is the most common type of colon cancer. To date, however, the prognostic values of m6A RNA methylation-related long non-coding RNAs (lncRNAs) in COAD are largely unknown.

Materials and Methods

The m6A-related lncRNAs were identified from The Cancer Genome Atlas (TCGA) data set. Univariate and multivariate Cox regression analyses were performed to explore the prognostic m6A-related lncRNAs. Consistent clustering analysis was performed to classify the COAD patients into different subgroups based on the expression of m6A-related lncRNAs. The potential biological functions as well as differences in the stemness index and tumor immune microenvironment between different subgroups were analyzed. The prognostic m6A-related lncRNAs were used to establish an m6A-related lncRNA risk model to predict prognosis and survival status.

Results

We identified 31 m6A-associated lncRNAs with prognostic values from the TCGA data set. Based on the expression of prognostic m6A-associated lncRNAs, TCGA-COAD patients were classified into three clusters using consistent clustering analysis. There was a low correlation of tumor stemness between the three clusters but a significant correlation with the tumor immune microenvironment as well as the tumor mutational load. Thirty-one prognostic-related m6A-associated lncRNAs were used to construct a risk model, which was further determined by survival analysis, receiver operating characteristic (ROC) curve, and univariate and multifactor Cox analysis. The m6A-related risk model demonstrates good performance in predicting prognosis and survival status. The model-based high-risk group exhibited poorer overall survival (OS) compared with the low-risk group.

Conclusion

In this study, we construct a risk model that consists of 31 m6A-related lncRNAs with independent prognostic values in COAD. Our study shows the critical roles of these 31 m6A-related lncRNAs in the tumor immune microenvironment, indicating the prospect of informing prognostic stratification and the development of immunotherapeutic strategies for COAD patients.

N6-methyladenosine (m6A) methyltransferase KIAA1429 accelerates the gefitinib resistance of non-small-cell lung cancer

N6-methyladenosine (m6A) modification has been convincingly identified to be a critical regulator in human cancer. However, the contribution of m6A to NSCLC gefitinib resistance is still largely unknown. Here, we screened and identified that m6A methyltransferase KIAA1429 was highly expressed in gefitinib-resistant NSCLC cells (PC9-GR), tissues, and closely related to unfavorable survival. Functionally, KIAA1429 accelerated the gefitinib resistance of NSCLC in vitro. Depletion of KIAA1429 repressed the tumor growth of PC9-GR cells in vivo. Mechanistically, KIAA1429 enhanced the mRNA stability of HOXA1 through targeting its 3'-untranslated regions (3'-UTR). Overall, our findings indicate that KIAA1429 plays essential oncogenic roles in NSCLC gefitinib resistance, which may provide a feasible therapeutic target for NSCLC.

Epigenetics: Roles and therapeutic implications of non-coding RNA modifications in human cancers

As next-generation sequencing (NGS) is leaping forward, more than 160 covalent RNA modification processes have been reported, and they are widely present in every organism and overall RNA type. Many modification processes of RNA introduce a new layer to the gene regulation process, resulting in novel RNA epigenetics. The commonest RNA modification includes pseudouridine (Ψ), N ⁷-methylguanosine (m7G), 5-hydroxymethylcytosine (hm5C), 5-methylcytosine (m5C), N ¹-methyladenosine (m1A), N ⁶-methyladenosine (m6A), and others. In this study, we focus on non-coding RNAs (ncRNAs) to summarize the epigenetic consequences of RNA modifications, and the pathogenesis of cancer, as diagnostic markers and therapeutic targets for cancer, as well as the mechanisms affecting the immune environment of cancer. In addition, we summarize the current status of epigenetic drugs for tumor therapy based on ncRNA modifications and the progress of bioinformatics methods in elucidating RNA modifications in recent years.

N6-methyladenosine modification of circCUX1 confers radioresistance of hypopharyngeal squamous cell carcinoma through caspase1 pathway

Hypopharyngeal squamous cell carcinoma (HPSCC) is one of the most common malignant tumors in otolaryngology head and neck surgery and is one of the worst prognostic malignant tumors. Endogenous circular RNA (circRNA) is more stable than mRNA, microRNA (miRNA), and long non-coding RNA (LncRNA) in exosomes, plasma, and urine, and participates in gene expression regulation to perform different functions. Therefore, circRNA is expected to become a biomarker and therapy target for many tumors. However, the expression and function of circRNA regulated by N6-methyladenosine (m6A) are still unclear in HNSCC. In this study, we demonstrated that a specific circRNA, circCUX1, was upregulated in HPSCC patients who are resistant to radiotherapy and predicts poor survival outcome. We further found that methyltransferase like 3 (METTL3) mediated the m6A methylation of circCUX1 and stabilizes its expression. Knockdown circCUX1 promotes the sensitivity of hypopharyngeal cancer cells to radiotherapy. In addition, circCUX1 binds to Caspase1 and inhibits its expression, resulting in a decrease in the release of inflammatory factors, thereby developing tolerance to radiotherapy. Our findings indicate that circCUX1 is a potential therapeutic target for radiotherapy tolerance in HPSCC patients.

Molecular Characteristics, Prognostic Value, and Immune Characteristics of m6A Regulators Identified in Head and Neck Squamous Cell Carcinoma

N6-methyladenosine (m⁶A) plays crucial roles in a diverse range of physiological and pathological processes, and it is believed that it tremendously promotes neoplasia and progression. However, knowledge of the molecular characteristics of m⁶A modification, its prognostic value, and the infiltration of immune cell populations in head and neck squamous cell carcinoma (HNSCC) is still insufficient. Therefore, a pan-cancer genomic analysis was systematically performed here by examining m⁶A regulators at the molecular level within 33 multiple cancer types, and the correlations between the expression of m⁶A molecules were researched using datasets from The Cancer Genome Atlas (TCGA). Based on the above analysis, insulin-like growth factor 2 mRNA-binding protein 2 (IGF2BP2) is upregulated in HNSCC and may serve as an independent prognostic factor of overall survival, thus showing potential as a prognostic biomarker in HNSCC. Genetic alteration analyses elucidated the reasons for the abnormal upregulation of IGF2BP2 in HNSCC. As a result, IGF2BP2 was selected for further univariate and multivariate analyses. The functions of the related genes were annotated through gene set enrichment analysis, and the activation states of multiple biological pathways were shown by gene set variation analysis. We found that LRRC59 and STIP1 may act as IGF2BP2-associated genes to have a regulatory function in the m⁶A modification. In addition, we found that the status of immune cell infiltration was correlated with the level of IGF2BP2 gene expression. Our results provide supplementation at the molecular level for epigenetic regulation in HNSCC and insight into effective immunotherapy targets and strategies.

Multidimensional crosstalk between RNA-binding proteins and noncoding RNAs in cancer biology

RNA-binding proteins (RBPs) are well-known to bind RNA via a set of RNA-binding domains (RBDs) and determine the fate and function of their RNA targets; inversely, some RBPs, in certain cases, may be modulated by the bound RNAs rather than regulate their RNA partners. Current proteome-wide studies reveal that almost half of RBPs have no canonical RBDs, and the discovery of tens of thousands of noncoding RNAs (ncRNAs), especially those with the size larger than 200 nt (namely long noncoding RNAs, lncRNAs), makes the crosstalk between RBPs and RNAs more complicated. It is clear that macromolecular complexes formed by RBP and RNA are not only a form of existence of their RBP and RNA components in cells, but also represent a functional entity through which those RBPs and regulatory ncRNAs participate in the construction of regulatory networks in organism. In this review, we summarize the multidimensional crosstalk between RBPs and ncRNAs in cancer and discuss how RBPs achieve their function via the bound ncRNAs in different aspects of gene expression as well as how RBPs direct modification and processing of ncRNAs, in order to better understand tumor biology and provide new insights into development of strategies for cancer therapy and early detection.

m6A Regulators Is Differently Expressed and Correlated With Immune Response of Esophageal Cancer

N6 methyladenosine (m6A) RNA methylation regulators play an important role in the development of tumors. However, their function in esophageal cancer (EC) has not been fully elucidated. Here, we analyzed the gene expression data of 24 major m6A RNA methylation regulators from 775 patients with EC from TCGA dataset. The present study showed the aberrations of m6A regulators in genome were correlated to prognosis in human ECs. Meanwhile, 17 m6A regulators showed increased expression in EC samples, including YTHDC1, IGF2BP2, FTO, METTL14, YTHDF3, RBM15, WTAP, HNRNPA2B1, HNRNPC, ALKBH5, YTHDF2, METTL16, IGF2BP3, VIRMA, RBM15B, YTHDF1, KIAA1429, HAKAI, and ZC3H13. Among them, we found HNRNPC, YTHDC2, WTAP, VIRMA, IGF2BP3, and HNRNPA2B1 were significantly correlated to worse outcomes and advanced stage in EC. Furthermore, we showed levels of m6A regulators is correlated with the expression of Immuno-regulators (Immunoinhibitors, Immunostimulators, and MHC molecules) and immune infiltration levels in EC. Bioinformatics further confirm m6A regulators were involved in regulating RNA splicing, RNA stability, and cell proliferation. Our study showed m6A regulators are promising targets and biomarkers for cancer immunotherapy in EC.

Circular RNA CircCOL5A1 Sponges the MiR-7-5p/Epac1 Axis to Promote the Progression of Keloids Through Regulating PI3K/Akt Signaling Pathway

Keloids, as a result of abnormal wound healing in susceptible individuals, are characterized by the hyper-proliferation of fibroblasts and exaggerated deposition of extracellular matrix. Current surgical and therapeutic modalities provide limited satisfactory results. Growing evidence has highlighted the roles of circRNAs in acting as miRNA sponges. However, up to date, the regulatory mechanism of circRNAs in the pathological process of keloids has rarely been reported. In this study, cell proliferation, cell migration, flow cytometry, western blotting, fluorescence in situ hybridization, dual-luciferase activity, and immunohistochemistry assays were applied to explore the roles and mechanisms of the circCOL5A1/miR-7-5p/Epac1 axis in the keloid. The therapeutic potential of circCOL5A1 was investigated by establishing keloid implantation models. The RT-qPCR result revealed that circCOL5A1 expression was obviously higher in keloid tissues and keloid fibroblasts. Subsequent cellular experiments demonstrated that circCOL5A1 knockdown repressed the proliferation, migration, extracellular matrix (ECM) deposition, whereas promoted cell apoptosis, through the PI3K/Akt signaling pathway. Furthermore, RNA-fluorescence in situ hybridization (RNA-FISH) illustrated that both circCOL5A1 and miR-7-5p were located in the cytoplasm. The luciferase reporter gene assay confirmed that exact binding sites were present between circCOL5A1 and miR-7-5p, as well as between miR-7-5p and Epac1. Collectively, the present study revealed that circCOL5A1 functioned as competing endogenous RNA (ceRNA) by adsorbing miR-7-5p to release Epac1, which contributed to pathological hyperplasia of keloids through activating the PI3K/Akt signaling pathway. Our data indicated that circCOL5A1 might serve as a novel promising therapeutic target and represent a new avenue to understand underlying pathogenesis for keloids.

Screening and functional prediction of differentially expressed circular RNAs in human glioma of different grades

Circular RNAs (circRNAs) have a critical regulatory function in human glioma. However, novel circRNAs related to different pathological grades of glioma and their crucial potential function are worth screening and prediction. CircRNA expression profiling was performed for 6 paired high- and low-grade glioma tissues and 5 adjacent normal brain tissues through next-generation sequencing. Quantitative real-time PCR (qRT-PCR) was conducted to validate circRNA expression. Bioinformatics analysis was performed, and circRNA-miRNA-mRNA networks were constructed. The expression and survival data of miRNAs and target genes were examined by GEPIA, Chinese Glioma Genome Atlas (CGGA), ONCOMINE, and cBioPortal databases. The RNA binding proteins (RBPs), open reading frames (ORFs) and N6-methyladenosine (m⁶A) modifications of the identified circRNAs were also predicted. Through multilevel research screening, 4 circRNAs (hsa_circ_0000915, hsa_circ_0127664, hsa_circ_0008362, and hsa_circ_0001467) were associated with glioma of different pathological grades and could be preferred candidates for subsequent functional analysis. Therefore, circRNAs are associated with the different pathological grades of glioma and reveal their potential critical regulatory function. CircRNAs might provide vital molecular biomarkers and potential therapeutic targets for glioma.