M6A-related bioinformatics analysis reveals that HNRNPC facilitates progression of OSCC via EMT

The potential mechanism of WT1-associated protein-induced N-6-methyladenosine modification of colony-stimulating factor 2 in the progression of oral squamous cell carcinoma by JAK/STAT3 pathway regulation

Colony-stimulating factor 2 (CSF2) plays a regulatory role in numerous cancers. However, there is needed to investigate the role of CSF2 in oral squamous cell carcinoma (OSCC) malignant phenotype and the specific mechanisms of CSF2 N-6-methyladenosine (m6A) modification. Therefore, we investigated the regulatory mechanism of m6A-modified CSF2 by WT1-associated protein (WTAP) in OSCC via qRT-PCR, western blot, WTAP and CSF2 overexpression in OSCC. In a panel of OSCCs, Kaplan-Meier plot analysis indicated that high expression of CSF2 was associated with poorer prognosis. Cell functional experiments revealed that enrichment of CSF2 promoted the proliferation and migration of OSCC cells by activating the JAK/STAT3 pathway, whereas the reduced expression of CSF2 resulted in the malignant decline of OSCC cells by blocking the JAK/STAT3 pathway. This study also confirmed that WTAP enhanced the m6A level of CSF2 and facilitated the expression of CSF2 and that CSF2 silencing blocked the invasive phenotype of OSCC cells and reversed the malignancy induced by WTAP overexpression. Overall, this study demonstrated that WTAP mediates the m6A modification of CSF2 and the JAK/STAT3 pathway, which plays an oncogenic role in the development of OSCC and can be a target for the treatment of patients with OSCC.

TFAIP6 facilitates hepatocellular carcinoma cell glycolysis through upregulating c-myc/PKM2 axis

Background

Hepatocellular carcinoma (HCC) is the most prevalent liver cancer. Despite of the improvement of therapies, the durable response rate and survival benefit are still limited for HCC patients. It's urgent to clarify the molecular mechanisms and find therapeutic strategies to improve the clinical outcome. TNFα-stimulated gene-6 (TNFAIP6) plays a critical role in the prognosis of various tumors, but its roles in HCC are still unclear.

Methods

Quantitative real-time PCR (qRT-PCR) and immunohistochemistry (IHC) analysis were employed to evaluate the clinical relevance of TNFAIP6 expressions in HCC patients. Cell counting kit-8 (CCK-8), Edu assay, and transwell assay were performed to evaluate the malignancy of HCC cells. Glucose uptake, lactate production, ATP production, extracellular acidification rate (ECAR) by Seahorse XF analyzer were employed to evaluate the role of TNFAIP6 in the regulation of aerobic glycolysis. The expressions of key proteins involved in glycolysis were examined by Western blot. Co-immunoprecipitation (Co-IP) and chromatin immunoprecipitation (ChIP) were used for protein-protein interactions or protein-RNA interactions respectively. Knockdown and overexpression of TNFAIP6 in HCC cells were employed for analyzing the functions of TNFAIP6 in HCC.

Results

TNFAIP6 was significantly upregulated in HCC and predicted a poor clinical prognosis. Knockdown of TNFAIP6 inhibited in vitro cell proliferation, invasion, migration, as well as glycolysis in HCC cells. Mechanistically, we clarified that TNFAIP6 interacted with heterogeneous nuclear ribonucleoprotein C (HNRNPC), stabilized c-Myc mRNA and upregulated pyruvate kinase M2 (PKM2) to promote glycolysis.

Conclusions

Our study reveals a molecular mechanism by which TNFAIP6 promotes aerobic glycolysis, which is beneficial for malignance of HCC and provides a potential clinical therapy for disease management.

N6-Methyladenosine reader HNRNPC-mediated downregulation of circITCH prevents miR-224-3p sequestering and contributes to tumorigenesis in nasopharyngeal carcinoma

BACKGROUND

N6-Methyladenosine (m6A) RNA methylation modulators are implicated in nasopharyngeal carcinoma (NPC). Circular RNAs (circRNAs) stimulate/inhibit the development of NPC by sponging microRNAs (miRNAs). Herein, m6A modifications affecting the circRNA/miRNA axis in NPC were explored.

METHODS

Twenty prognostic m6A RNA methylation regulators were identified from 504 head/neck squamous cell carcinoma and 44 normal samples from The Cancer Genome Atlas (TCGA). Differentially expressed miRNAs were screened from the TCGA and Gene Expression Omnibus (GEO) databases. RNA-binding protein (RBP)-circRNA and circRNA-miRNA interactive pairs were verified using RBPmap and RNAhybrid, respectively. The RBP/circRNA/miRNA network was constructed using Cytoscape. Furthermore, CircITCH (hsa_circ_00059948), HNRNPC, and miR-224-3p expressions were detected by western blotting and quantitative polymerase chain reaction. The role of circITCH in NPC was examined using apoptosis, scratch wound healing, transwell invasion, and cell counting kit-8 assays. Finally, CircITCH-miR-224-3p and circITCH-HNRNPC interactions were assessed by dual-luciferase reporter and RNA-immunoprecipitation (RIP) assays, respectively.

RESULTS

Bioinformatics analysis revealed that high pathological grade, late-stage tumors, and low survival were associated with increased HNRNPC expression. MiR-224-3p was upregulated in NPC and sequestered by circITCH. Construction of the RBP/circRNA/miRNA network highlighted the HNRNPC/circITCH/miR-224-3p axis. In vitro experiments demonstrated decreased circITCH expression and increased HNRNPC and miR-224-3p expressions in NPC. In NPC cells overexpressing circITCH, HNRNPC and miR-224-3p expressions were significantly decreased. Dual-luciferase assays demonstrated a targeting relationship between circITCH and miR-224-3p, and RIP assays demonstrated interaction of HNRNPC targets with circITCH.

CONCLUSION

CircITCH overexpression inhibited NPC progression by sequestering miR-224-3p, and HNRNPC reduced circITCH expression through direct interaction.

Identification of splicing factors signature predicting prognosis risk and the mechanistic roles of novel oncogenes in HNSCC

Head and neck squamous cell carcinoma (HNSCC) is the most frequent subtype of head and neck cancer, generally with a poor prognosis and limited therapeutic options due to its highly heterogeneous malignancy. In this study, we screened functional splicing regulatory RNA binding proteins (RBPs) that were closely related with the prognosis of HNSCC patients and showed significant expression differences between HNSCC tumors and normal tissues. Based on this finding, we chose six candidate genes (HNRNPC, ZCRB1, RBM12B, SF3A2, SF3B3, and SRSF11) to generate a prognostic prediction model and validated the accuracy of the prognostic model for predicting patient survival outcomes. We found that the risk score predicted by our model can serve as an independent prognostic predictor. Notably, HNSCC tumors showing higher expression of SF3B3, HNRNPC, or ZCRB1 possessed higher risk scores in the discovered prediction model. The investigation of the underlying mechanism validated that knockdown of SF3B3, HNRNPC, and ZCRB1 separately induced a substantial impairment of HNSCC cell survival. Conversely, overexpression of each of the three genes promoted tumor cellular proliferation. High throughput RNA sequencing analysis revealed that changes in the expression of SF3B3 and HNRNPC remarkably affected alternative splicing of genes related to cell cycle regulation, whereas the depletion of ZCRB1 contributed to aberrant splicing events involving in DNA damage response. In addition, the prognostic prediction model's risk score was demonstrated to be related with the immune infiltration score. Particularly, SF3B3 has a negative correlation with CD8A expression. Therefore, our findings provide promising prognosis predictors and potential therapeutic targets for better treatment efficacy of HNSCC.

The m6A reader HNRNPC predicts adverse prognosis and promotes the progression of colorectal cancer

BACKGROUND

As a critical m6A RNA methylation regulator, HNRNPC has been revealed to serve as potential biomarkers in various human cancers. The specific expression and significance of HNRNPC in colorectal cancer remain unknown.

OBJECTIVE

This study aimed to confirm HNRNPC expression level and evaluate its function in colorectal cancer progression.

METHODS

101 paired tissue samples were collected from colorectal cancer patients. HNRNPC levels in colorectal cancer were detected using PCR. CCK8 and transwell assays were conducted to estimate the effect of HNRNPC on cell growth and metastasis with the regulation of HNRNPC by cell transfection.

RESULTS

Upregulated HNRNPC was observed in colorectal cancer compared with normal tissues and cells. The higher HNRNPC levels in tumor tissues were associated with the advanced TNM stage and positive lymph node metastasis. Meanwhile, HNRNPC upregulation could indicate adverse outcomes of colorectal cancer patients. In vitro, the knockdown of HNRNPC significantly suppressed the proliferation, migration, and invasion of colorectal cancer cells.

CONCLUSIONS

Upregulated HNRNPC served as a biomarker for the prognosis and development of colorectal cancer, which provides a novel therapeutic target for colorectal cancer.

Identification and Validation of Genes Related to RNA Methylation Modification in Diabetic Retinopathy

PURPOSE

To identify and validate the differentially expressed genes related to RNA methylation modification in diabetic retinopathy.

METHODS

The data sets GSE12610 and GSE111465 related to diabetic retinopathy in the Gene Expression Omnibus were selected. The R software package was used to identify differentially expressed genes related to RNA methylation modification in diabetic retinopathy. Protein-protein interaction network was constructed to explore the interactions between proteins and predict proteins. Then, Gene Ontology annotation analysis and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis were used to analyze the potential enrichment pathways and clarify the biological functions of these genes. In addition, the correlation between them and immune cells was visualized, and receiver operating characteristic curves were drawn to evaluate the diagnostic performance of each one of them for diabetic retinopathy. To verify the differentially expressed genes, the mRNA expression of rat retinal vascular endothelial cells cultured in low and high glucose medium separately were detected by RT-qPCR.

RESULTS

The expression of Lrpprc, Nsun4, Nsun6 and Trdmt1 were significantly up-regulated in diabetic retinopathy samples, while the expression of Cbll1, Hnrnpc, Mettl3 and Wtap were significantly down-regulated. Differentially expressed genes were mainly enriched in the RNA-methylation-medication pathways and biological function. The results of immune infiltration analysis proved that eosinophils aggregated more in diabetic group, while T cells follicular helper aggregated more in normal samples. These genes of Cbll1 (AUC = 0.986), Hnrnpc (AUC = 0.819), Lrpprc (AUC = 0.806), Mettl3 (AUC = 0.917), Nsun4 (AUC = 0.819), Nsun6 (AUC = 0.819), Trdmt1 (AUC = 0.972) and Wtap (AUC = 0.972) were respectively used as the diagnostic basis of diabetic retinopathy. According to the RT-qPCR results, the expression of Mettl3 was significantly down-regulated (p < 0.0005) in cells cultured in high glucose, while Trdmt1 (p < 0.05), Nsun4 (p < 0.05) and Nsun6 (p < 0.05) were significantly up-regulated.

CONCLUSION

Differentially expressed genes such as Mettl3, Nsun4, Nsun6, and Trdmt1 should be conducted to explore, and the role of RNA methylation in the process of diabetic retinopathy would be revealed in-depth.

METTL14 inhibits malignant progression of oral squamous cell carcinoma by targeting the autophagy-related gene RB1CC1 in an m6A-IGF2BP2-dependent manner

N6-methyladenosine (m6A) plays crucial roles in tumorigenesis and autophagy. However, the underlying mechanisms mediated by m6A and autophagy in the malignant progression of oral squamous cell carcinoma (OSCC) remain unclear. In the present study, we revealed that down-regulated expression of METTL14 was correlated with advanced clinicopathological characteristics and poor prognosis in OSCC. METTL14 knockdown significantly inhibited autophagy and facilitated malignant progression in vitro, and promoted tumor growth and metastasis in vivo. A cell model of rapamycin-induced autophagy was established to identify RB1CC1 as a potential target gene involved in m6A-regulated autophagy in OSCC, through RNA sequencing and methylated RNA immunoprecipitation sequencing (meRIP-seq) analysis. Mechanistically, we confirmed that METTL14 posttranscriptionally enhanced RB1CC1 expression in an m6A-IGF2BP2-dependent manner, thereby affecting autophagy and progression in OSCC, through methylated RNA immunoprecipitation qRT-PCR (meRIP-qPCR), RNA stability assays, mutagenesis assays and dual-luciferase reporter. Collectively, our findings demonstrated that METTL14 serves as an OSCC suppressor by regulating the autophagy-related gene RB1CC1 through m6A modification, which may provide a new insight for the diagnosis and therapy of OSCC.

OCT4's role and mechanism underlying oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC), a common malignancy of the head and neck, ranks sixth worldwide in terms of cancers with the most negative impact, owing to tumor relapse rates, cervical lymphnode metastasis, and the lack of an efficacious systemic therapy. Its prognosis is poor, and its mortality rate is high. Octamer-binding transcription factor 4 (OCT4) is a member of the Pit-Oct-Unc (POU) family and is a key reprogramming factor that produces a marked effect in preserving the pluripotency and self-renewal state of embryonic stem cells (ESCs). According to recent studies, OCT4 participates in retaining the survival of OSCC cancer stem cells (CSCs), which has far-reaching implications for the occurrence, recurrence, metastasis, and prognosis of oral carcinogenesis. Therefore, we summarize the structure, subtypes, and function of OCT4 as well as its role in the occurrence, progression, and prognosis of OSCC.

RNA N6-methyladenosine (m6A) modification in HNSCC: molecular mechanism and therapeutic potential

Head and neck squamous cell carcinoma ranks seventh in incidence of malignant tumours in the world. Although there are treatments including surgery, radiotherapy and chemotherapy, targeted therapy and immunotherapy, drug resistance to treatment is caused by various reasons, and the survival rate of patients remains frustrating. To overcome the bottleneck of treatment at this stage, it is urgent to identify possible diagnostic and prognostic markers. N6-methyladenosine is a methylation modification on the sixth N atom of adenine which is the most abundant epitope transcriptome modification in mammalian genes. N6-methyladenosine modification is reversible and results from the interaction among writers, erasers and readers. A large number of studies have proven that N6-methyladenosine modification has important significance in promoting the progression and treatment of tumours and have made great progress in research. In this review, we introduce how N6-methyladenosine modification promotes the occurrence and development of tumours, the mechanism of drug resistance, and new findings of N6-methyladenosine modification in radiotherapy and chemotherapy, immunotherapy, and targeted therapy. N6-methyladenosine modification provides more possibilities for improving the overall survival rate and prognosis of patients.

FOXM1a Isoform of Oncogene FOXM1 Is a Tumor Suppressor Suppressed by hnRNP C in Oral Squamous Cell Carcinoma

FOXM1 is an oncogenic transcriptional factor and includes several isoforms generated by alternative splicing. Inclusion of alternative exon 9 produces FOXM1a, a transcriptionally inactive isoform. However, the role of FOXM1a in tumorigenesis remains unknown. In addition, the regulatory mechanisms of exon 9 splicing are also unclear. In the present study, we found that overexpression of FOXM1a significantly reduced cell proliferation and colony formation of oral squamous cell carcinoma (OSCC) cell proliferation in vitro. Importantly, OSCC cells with FOXM1a overexpression showed significantly slower tumor formation in nude mice. Moreover, we identified a U-rich exonic splicing suppressor (ESS) which is responsible for exon 9 skipping. Splicing factor heterogeneous nuclear ribonucleoprotein C (hnRNP C) can bind to the ESS and suppress exon 9 inclusion and FOXM1a expression. Silence of hnRNP C also significantly suppresses OSCC cell proliferation. HnRNP C is significantly co-expressed with FOXM1 in cancers. Our study uncovered a novel regulatory mechanism of oncogene FOXM1 expression in OSCC.

N6-methyladenosine with immune infiltration and PD-L1 in hepatocellular carcinoma: novel perspective to personalized diagnosis and treatment

Background

Increasing evidence elucidated N6-methyladenosine (m6A) dysregulation participated in regulating RNA maturation, stability, and translation. This study aimed to demystify the crosstalk between m6A regulators and the immune microenvironment, providing a potential therapeutic target for patients with hepatocellular carcinoma (HCC).

Methods

Totals of 371 HCC and 50 normal patients were included in this study. GSE121248 and GSE40367 datasets were used to validate the expression of HNRNPC. The R package "ConsensusClusterPlus" was performed to screen consensus clustering types based on the expression of m6A regulators in HCC. The R package "pheatmap", "immunedeconv", "survival", "survminer" and "RMS" were applied to investigate the expression, immunity, overall survival, and clinical application in different clusters and expression groups. Comprehensive analysis of HNRNPC in pan-cancer was conducted by TIMER2 database. Besides, HNRNPC mRNA and protein expression were verified by qRT-PCR and immunohistochemistry analysis.

Results

Most of m6A regulators were over-expressed excerpt for ZC3H13 in HCC. Three independent clusters were screened based on m6A regulators expression, and the cluster 2 had a favorable prognosis in HCC. Then, the cluster 2 was positively expression in macrophage, hematopoietic stem cell, endothelial cell, and stroma score, while negatively in T cell CD4⁺ memory and mast cell. We identified HNRNPC was an independent prognostic factor in HCC, and nomogram performed superior application value for clinical decision making. Moreover, PD-L1 was significantly up-regulated in HCC tissues, cluster 1, and cluster 3, and we found PD-L1 expression was positively correlated with HNRNPC. Patients with HCC in high-expression groups was associated with tumor-promoting cells. Besides, HNRNPC was correlated with prognosis, TMB, and immune checkpoints in cancers. Particularly, the experiments confirmed that HNRNPC was positively expression in HCC cells and tissues.

Conclusion

The m6A regulators play irreplaceable roles in prognosis and immune infiltration in HCC, and the relationship of HNRNPC and PD-L1 possesses a promising direction for therapeutic targets of immunotherapy response. Exploration of m6A regulators pattern could be build the prognostic stratification of individual patients and move toward to personalized treatment.

Analyses of m6A regulatory genes and subtype classification in atrial fibrillation

Objective

To explore the role of m6A regulatory genes in atrial fibrillation (AF), we classified atrial fibrillation patients into subtypes by two genotyping methods associated with m6A regulatory genes and explored their clinical significance.

Methods

We downloaded datasets from the Gene Expression Omnibus (GEO) database. The m6A regulatory gene expression levels were extracted. We constructed and compared random forest (RF) and support vector machine (SVM) models. Feature genes were selected to develop a nomogram model with the superior model. We identified m6A subtypes based on significantly differentially expressed m6A regulatory genes and identified m6A gene subtypes based on m6A-related differentially expressed genes (DEGs). Comprehensive evaluation of the two m6A modification patterns was performed.

Results

The data of 107 samples from three datasets, GSE115574, GSE14975 and GSE41177, were acquired from the GEO database for training models, comprising 65 AF samples and 42 sinus rhythm (SR) samples. The data of 26 samples from dataset GSE79768 comprising 14 AF samples and 12 SR samples were acquired from the GEO database for external validation. The expression levels of 23 regulatory genes of m6A were extracted. There were correlations among the m6A readers, erasers, and writers. Five feature m6A regulatory genes, ZC3H13, YTHDF1, HNRNPA2B1, IGFBP2, and IGFBP3, were determined (p < 0.05) to establish a nomogram model that can predict the incidence of atrial fibrillation with the RF model. We identified two m6A subtypes based on the five significant m6A regulatory genes (p < 0.05). Cluster B had a lower immune infiltration of immature dendritic cells than cluster A (p < 0.05). On the basis of six m6A-related DEGs between m6A subtypes (p < 0.05), two m6A gene subtypes were identified. Both cluster A and gene cluster A scored higher than the other clusters in terms of m6A score computed by principal component analysis (PCA) algorithms (p < 0.05). The m6A subtypes and m6A gene subtypes were highly consistent.

Conclusion

The m6A regulatory genes play non-negligible roles in atrial fibrillation. A nomogram model developed by five feature m6A regulatory genes could be used to predict the incidence of atrial fibrillation. Two m6A modification patterns were identified and evaluated comprehensively, which may provide insights into the classification of atrial fibrillation patients and guide treatment.

METTL3/14 and IL-17 signaling contribute to CEBPA-DT enhanced oral cancer cisplatin resistance

OBJECTIVES

Oral squamous cell carcinoma (OSCC) is the most common head and neck cancer. Chemotherapy has been recognized as an optional combination treatment, which enhance the overall survival of OSCC patients. However, the majority of patients would suffer therapeutic resistance, which led to the treatment failure and poor prognosis.

MATERIALS AND METHODS

To explore the mechanism of chemoresistance in OSCC, we first constructed two chemoresistant cell lines using Cal27 and HSC4. Then MeRIP sequencing together with bioinformatics analysis and a series of in vitro experiments were used to assess the possible regulation manner of RNA methylation on OSCC chemoresistance. Finally, xenograft models were constructed to confirm the relationship among OSCC chemoresistance.

RESULTS

METTL3/METTL14 upregulation could enhance OSCC chemoresistance. CEBPA-DT overexpression could regulate METTL3/METTL14 expression and further activate downstream BHLHB9. CEBPA-DT overexpression could inhibit the activity of IL-17 signaling, resulting in the homeostasis breakdown of immune infiltration and cytokine release. CEBPA-DT overexpression could significantly enhance chemoresistance through METTL3/METTL14/BHLHB9 in vivo, which accelerated the tumor growth.

CONCLUSIONS

Our results suggest that CEBPA-DT might regulate OSCC chemoresistance through BHLHB9 gene manipulated by METTL3/METTL14 as well as through IL-17 signaling inhibition, which may contribute to the assessment of potential therapeutic targets in OSCC chemoresistance.

Identification and vitro verification of the potential drug targets of active ingredients of Chonglou in the treatment of lung adenocarcinoma based on EMT-related genes

Lung adenocarcinoma (LUAD) is the main histological type of lung cancer with an unfavorable survival rate. Metastasis is the leading LUAD-related death with Epithelial-Mesenchymal Transition (EMT) playing an essential role. The anticancer efficacies of the active ingredients in Chonglou have been widely reported in various cancers. However, the potential therapeutic targets of the Chonglou active ingredients in LUAD patients remain unknown. Here, the network pharmacology and bioinformatics were performed to analyze the associations of the clinical characteristics, immune infiltration factors and m⁶A-related genes with the EMT-related genes associated with LUAD (EMT-LUAD related genes), and the molecular docking, STRING, GO, and KEGG enrichment for the drug targets of Chonglou active ingredients associated with EMT (EMT-LUAD-Chonglou related genes). And, cell viability analysis and cell invasion and infiltration analysis were used to confirm the theoretical basis of this study. A total of 166 EMT-LUAD related genes were identified and a multivariate Cox proportional hazards regression model with a favorable predictive accuracy was constructed. Meanwhile, the immune cell infiltration, immune cell subsets, checkpoint inhibitors and the expression of m⁶A-related genes were significantly associated with the risk scores for EMT-LUAD related genes with independent significant prognostic value of all included LUAD patients. Furthermore, 12 EMT-LUAD-Chonglou related genes with five core drug targets were identified, which participated in LUAD development through extracellular matrix disassembly, collagen metabolic process, collagen catabolic process, extracellular matrix organization, extracellular structure organization and inflammatory response. Moreover, we found that the active ingredients of Chonglou could indeed inhibit the progression of lung adenocarcinoma cells. These results are oriented towards EMT-related genes to achieve a better understanding of the role of Chonglou and its targets in osteosarcoma development and metastasis, thus guiding future preclinical studies and facilitating clinical translation of LUAD treatment.

Glycolysis-Related Gene Analyses Indicate That DEPDC1 Promotes the Malignant Progression of Oral Squamous Cell Carcinoma via the WNT/β-Catenin Signaling Pathway

Increasing evidence suggests that aerobic glycolysis is related to the progression of oral squamous cell carcinoma (OSCC). Hence, we focused on glycolysis-related gene sets to screen for potential therapeutic targets for OSCC. The expression profiles of OSCC samples and normal controls were obtained from The Cancer Genome Atlas (TCGA). Then, the differentially expressed gene sets were selected from the official GSEA website following extraction of the differentially expressed core genes (DECGs). Subsequently, we tried to build a risk model on the basis of DECGs to predict the prognosis of OSCC patients via Cox regression analysis. Furthermore, crucial glycolysis-related genes were selected to explore their biological roles in OSCC. Two active glycolysis-related pathways were acquired and 66 DECGs were identified. Univariate Cox regression analysis showed that six genes, including HMMR, STC2, DDIT4, DEPDC1, SLC16A3, and AURKA, might be potential prognostic factors. Subsequently, a risk formula consisting of DEPDC1, DDIT4, and SLC16A3 was established on basis of the six molecules. Furthermore, DEPDC1 was proven to be related to advanced stage cancer and lymph node metastasis. Moreover, functional experiments suggested that DEPDC1 promoted the aerobic glycolysis, migration, and invasion of OSCC via the WNT/β-catenin pathway. The risk score according to glycolysis-related gene expression might be an independent prognostic factor in OSCC. In addition, DEPDC1 was identified as playing a carcinogenic role in OSCC progression, suggesting that DEPDC1 might be a novel biomarker and therapeutic target for OSCC.

Identification of m6a-related signature genes in esophageal squamous cell carcinoma by machine learning method

Background: We aimed to construct and validate the esophageal squamous cell carcinoma (ESCC)-related m6A regulators by means of machine leaning. Methods: We used ESCC RNA-seq data of 66 pairs of ESCC from West China Hospital of Sichuan University and the transcriptome data extracted from The Cancer Genome Atlas (TCGA)-ESCA database to find out the ESCC-related m6A regulators, during which, two machine learning approaches: RF (Random Forest) and SVM (Support Vector Machine) were employed to construct the model of ESCC-related m6A regulators. Calibration curves, clinical decision curves, and clinical impact curves (CIC) were used to evaluate the predictive ability and best-effort ability of the model. Finally, western blot and immunohistochemistry staining were used to assess the expression of prognostic ESCC-related m6A regulators. Results: 2 m6A regulators (YTHDF1 and HNRNPC) were found to be significantly increased in ESCC tissues after screening out through RF machine learning methods from our RNA-seq data and TCGA-ESCA database, respectively, and overlapping the results of the two clusters. A prognostic signature, consisting of YTHDF1 and HNRNPC, was constructed based on our RNA-seq data and validated on TCGA-ESCA database, which can serve as an independent prognostic predictor. Experimental validation including the western and immunohistochemistry staining were further successfully confirmed the results of bioinformatics analysis. Conclusion: We constructed prognostic ESCC-related m6A regulators and validated the model in clinical ESCC cohort as well as in ESCC tissues, which provides reasonable evidence and valuable resources for prognostic stratification and the study of potential targets for ESCC.

Breast cancer combined prognostic model based on lactate metabolism genes

To investigate the impact of lactate metabolism genes, lactate metabolism-related genes (LMRG), and immune infiltrating cells on the prognosis of breast cancer. LMRG was identified via single-cell sequencing. Immune cell infiltration was obtained by the CIBERSORT method. The prognostic genes were chosen by cox regression and the least absolute selection operator approach. lactate metabolism-associated immune-infiltrating cells was determined by difference analysis. The GSE20685 dataset was used as an external validation cohort. The model's prognostic usefulness was evaluated utilizing survival, immunological microenvironment, and drug sensitivity assessments. NDUFAF6 was most associated with breast cancer prognosis. We obtained a total of 450 LMRG. SUSD3, IL18, MAL2, and CDKN1C comprised the Model2. NK cell activation was most relevant to lactate metabolism. The combined prognostic model outperformed the individual model, with the area under the curve ranging from 0.7 to 0.8 in all three cohorts. The lactate metabolism-related combination model assisted in evaluating breast cancer prognosis, providing new insights for treatment, particularly immunotherapy.

Emerging roles of alternative RNA splicing in oral squamous cell carcinoma

Alternative RNA splicing (ARS) is an essential and tightly regulated cellular process of post-transcriptional regulation of pre-mRNA. It produces multiple isoforms and may encode proteins with different or even opposite functions. The dysregulated ARS of pre-mRNA contributes to the development of many cancer types, including oral squamous cell carcinoma (OSCC), and may serve as a biomarker for the diagnosis and prognosis of OSCC and an attractive therapeutic target. ARS is mainly regulated by splicing factors, whose expression is also often dysregulated in OSCC and involved in tumorigenesis. This review focuses on the expression and roles of splicing factors in OSCC, the alternative RNA splicing events associated with OSCC, and recent advances in therapeutic approaches that target ARS.

Identification of a five m6A-relevant mRNAs signature and risk score for the prognostication of gastric cancer

BACKGROUND

N6-methyladenosine (m6A) is the most abundant form of methylation modification in eukaryotic cell messenger RNA (mRNA). However, the role of m6A in gastric cancer (GC), which is one of the most common gastrointestinal malignancies, is unclear. In this study, m6A-relevant mRNA signatures and risk scores were determined to predict the prognosis of GC.

METHODS

The expression profiles and clinical information of 367 patients were downloaded from The Cancer Genome Atlas (TCGA). Cluster analysis and univariate Cox analysis were performed to identify the regulatory factors of RNA methylation associated with GC prognosis. A co-expression network was constructed using the WGCNA package in R. The correlations between module eigengenes and clinical traits were then calculated to identify the relevant modules. We used univariate Cox analysis to screen for genes that are significantly associated with prognosis in the module. We identified hub genes by least absolute shrinkage and selection operator (LASSO) and multivariate analysis and developed a Cox prognostic model. Finally, the hub gene expression values weighted by the coefficients from the LASSO regression were applied to generate a risk score for each patient, and receiver operating characteristic (ROC) and Kaplan-Meier curves were used to assess the prognostic capacity of the risk scores. The asporin (ASPN) gene in GC cell lines was verified via quantitative polymerase chain reaction (qPCR) and Western blot. Moreover, 5-ethynyl-2'-deoxyuridine (EdU) and transwell assays were applied to evaluate the effects of the proliferation, migration, and invasion abilities in GC cells after ASPN knockdown. Western blot verified the effects of ASPN on the phosphoinositide 3-kinase (PI3K)/serine/threonine kinase (AKT)/mechanistic target of rapamycin kinase (mTOR) pathway and epithelial-mesenchymal transition (EMT) pathway-related gene expression.

RESULTS

Our results indicated that AARD, ASPN, SLAMF9, MIR3117 and DUSP1 were hub genes affecting the prognosis of GC patients. Besides, we found that ASPN expression was upregulated in GC cells. The knockdown of ASPN expression suppressed GC cell proliferation, migration, and invasion by deactivating the PI3K/AKT/mTOR and EMT pathways, respectively.

CONCLUSIONS

Our findings indicated that ASPN participates in the biological process of GC as an oncogene and may be a promising biomarker in GC.

Transcriptome-wide m6A methylome analysis uncovered the changes of m6A modification in oral pre-malignant cells compared with normal oral epithelial cells

As the most common post-transcriptional RNA modification, m⁶A methylation extensively regulates the structure and function of RNA. The dynamic and reversible modification of m⁶A is coordinated by m⁶A writers and erasers. m⁶A reader proteins recognize m⁶A modification on RNA, mediating different downstream biological functions. mRNA m⁶A modification and its corresponding regulators play an important role in cancers, but its characteristics in the precancerous stage are still unclear. In this study, we used oral precancerous DOK cells as a model to explore the characteristics of transcriptome-wide m⁶A modification and major m⁶A regulator expression in the precancerous stage compared with normal oral epithelial cell HOEC and oral cancer cell SCC-9 through MeRIP-seq and RT-PCR. Compared with HOEC cells, we found 1180 hyper-methylated and 1606 hypo-methylated m⁶A peaks and 354 differentially expressed mRNAs with differential m⁶A peaks in DOK cells. Although the change of m⁶A modification in DOK cells was less than that in SCC-9 cells, mRNAs with differential m⁶A in both cell lines were enriched into many identical GO terms and KEGG pathways. Among the 20 known m⁶A regulatory genes, FTO, ALKBH5, METTL3 and VIRMA were upregulated or downregulated in DOK cells, and the expression levels of 10 genes such as METTL14/16, FTO and IGF2BP2/3 were significantly changed in SCC-9 cells. Our data suggest that precancerous cells showed, to some extent, changes of m⁶A modification. Identifying some key m⁶A targets and corresponding regulators in precancerous stage may provide potential intervention targets for the prevention of cancer development through epigenetic modification in the future.

Construction of a prognostic model related to copper dependence in breast cancer by single-cell sequencing analysis

Purpose: To explore the clinical significance of copper-dependent-related genes (CDRG) in female breast cancer (BC).

Methods: CDRG were obtained by single-cell analysis of the GSE168410 dataset in the Gene Expression Omnibus (GEO) database. According to a 1:1 ratio, the Cancer Genome Atlas (TCGA) cohort was separated into a training and a test cohort randomly. Based on the training cohort, the prognostic model was built using COX and Lasso regression. The test cohort was used to validate the model. The GSE20685 dataset and GSE20711 dataset were used as two external validation cohorts to further validate the prognostic model. According to the median risk score, patients were classified as high-risk or low-risk. Survival analysis, immune microenvironment analysis, drug sensitivity analysis, and nomogram analysis were used to evaluate the clinical importance of this prognostic model.

Results: 384 CDRG were obtained by single-cell analysis. According to the prognostic model, patients were classified as high-risk or low-risk in both cohorts. The high-risk group had a significantly worse prognosis. The area under the curve (AUC) of the model was around 0.7 in the four cohorts. The immunological microenvironment was examined for a possible link between risk score and immune cell infiltration. Veliparib, Selumetinib, Entinostat, and Palbociclib were found to be more sensitive medications for the high-risk group after drug sensitivity analysis.

Conclusion: Our CDRG-based prognostic model can aid in the prediction of prognosis and treatment of BC patients.

The prognostic value and multiomic features of m6A-related risk signature in lung adenocarcinoma

OBJECTIVES

N⁶-methyladenosine (m6A), a predominant RNA modification, has been recently linked to messenger RNA splicing, stability and expression, and its dysregulation may be important in the initiation as well as development of human cancers. The current study was proposed to investigate the clinico-pathological value and multiomic characteristics of m6A-linked genes in the diagnosis and prognosis of lung adenocarcinoma (LUAD).

METHODS

The expression levels and mutation types of 21 previously identified m6A regulators were analyzed using the TCGA (The Cancer Genome Atlas) database. The patients were categorized into two groups, a training group (n=392) and a testing group (n=98). Next, the prognostic score of m6A regulators was determined by the Cox survival analysis and a regression model of LASSO to develop a risk profile for patients with LUAD. Moreover, features of risk signature, including chemosensitivity, tumor immune microenvironment and genetic mutation, were also explored.

RESULTS

In total, 18 of 21 m6A regulators showed significantly differential expression in LUAD (P<0.05). Among them, 6 genes were observed to be associated with the Overall Survival (OS) of patients with LUAD. Three genes (IGF2BP1 and 2, and HNRNPC) were further evaluated as a prognostic signature in LUAD. Patients, grouped as high risk based on the median of risk score, had poorer OS in comparison with those in low-risk group (P<0.05). The accuracy of our prognostic signatures was high: the AUC were 0.67, 0.59, 0.64 (training set), and 0.65, 0.69, 0.64 (testing set) at survival of 1- , 3- and 5-year, respectively. The prognostic performance of IGF2BP1, IGF2BP2 and HNRNPC was successfully validated in two independent external cohorts. High-risk score was an indicator of chemoresistance, TP53 mutation and increased infiltration of immune cells, and in vitro assessment of the cellular function of HNRNPC confirmed that the gene is involved in cell proliferation and invasion.

CONCLUSION

The prognostic signature based on m6A regulators might provide novel insights into prognostic assessment and individualized treatment for patients with LUAD.

A Novel Immune-Related Gene Signature to Identify the Tumor Microenvironment and Prognose Disease Among Patients With Oral Squamous Cell Carcinoma Patients Using ssGSEA: A Bioinformatics and Biological Validation Study

Oral squamous cell carcinoma (OSCC) is the most invasive oral malignancy in adults and is associated with a poor prognosis. Accurate prognostic models are urgently needed, however, knowledge of the probable mechanisms behind OSCC tumorigenesis and prognosis remain limited. The clinical importance of the interplay between the immune system and tumor microenvironment has become increasingly evident. This study explored immune-related alterations at the multi-omics level to extract accurate prognostic markers linked to the immune response and presents a more accurate landscape of the immune genomic map during OSCC. The Cancer Genome Atlas (TCGA) OSCC cohort (n = 329) was used to detect the immune infiltration pattern of OSCC and categorize patients into two immunity groups using single-sample gene set enrichment analysis (ssGSEA) and hierarchical clustering analysis. Multiple strategies, including lasso regression (LASSO), Cox proportional hazards regression, and principal component analysis (PCA) were used to screen clinically significant signatures and identify an incorporated prognosis model with robust discriminative power on the survival status of both the training and testing set. We identified two OSCC subtypes based on immunological characteristics: Immunity-high and immunity low, and verified that the categorization was accurate and repeatable. Immunity_ high cluster with a higher immunological and stromal score. 1047 differential genes (DEGs) integrate with immune genes to obtain 319 immue-related DEGs. A robust model with five signatures for OSCC patient prognosis was established. The GEO cohort (n = 97) were used to validate the risk model’s predictive value. The low-risk group had a better overall survival (OS) than the high-risk group. Significant prognostic potential for OSCC patients was found using ROC analysis and immune checkpoint gene expression was lower in the low-risk group. We also investigated at the therapeutic sensitivity of a number of frequently used chemotherapeutic drugs in patients with various risk factors. The underlying biological behavior of the OSCC cell line was preliminarily validated. This study characterizes a reliable marker of OSCC disease progression and provides a new potential target for immunotherapy against this disease.

N6-methyladenosine (m6A) reader IGF2BP2 stabilizes HK2 stability to accelerate the Warburg effect of oral squamous cell carcinoma progression

PURPOSE

The N⁶-methyladenosine (m⁶A) has been involved in the regulation of cell proliferation and metastasis in multiple cancers. However, the biological significance of m⁶A reader IGF2BP2 in oral squamous cell carcinoma (OSCC) and the mechanism of IGF2BP2 itself have not been fully investigated.

METHODS

The cellular phenotypes of OSCC cells were determined by CCK-8 and transwell migration assays. The energy metabolism was detected using glucose uptake/lactate production assay and extracellular acidification rate analysis. The molecular interaction was tested by  RNA immunoprecipitation  assay.

RESULTS

Here, results indicated that IGF2BP2 was up-regulated in OSCC and that it acted as a predictor of poor prognosis. IGF2BP2 promoted the proliferation, migration and Warburg effect of OSCC cells in vitro. Mechanistical assays illustrated that IGF2BP2 directly interacted with HK2 mRNA by binding the 3'-UTR m⁶A site. Moreover, IGF2BP2 positively promoted the stability of HK2 mRNA and thus the protein level of HK2 increased upon IGF2BP2 overexpression.

CONCLUSIONS

In conclusion, the IGF2BP2/m⁶A/HK2 axis accelerated the abnormal energy metabolism of OSCC. Taken together, these findings revealed a novel mechanism by which IGF2BP2 functions in OSCC progression, which may provide new therapy options for OSCC patients.

FCGBP Is a Promising Prognostic Biomarker and Correlates with Immunotherapy Efficacy in Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common malignancies of the head and neck. In OSCC patients, the prognosis was dramatically different. In this research, we aimed to study the expressions and prognostic values of IgG Fc binding protein (FCGBP) in OSCC patients. The expression of FCGBP was analyzed using TCGA datasets and GEO datasets. FCGBP was evaluated for its predictive significance in OSCC patients by the use of a Kaplan-Meier and Cox regression model. Enrichment analysis for the GO and KEGG databases were conducted. CIBERSORT used TCGA datasets to show immune cell infiltration. In addition, researchers looked into the relationships between FCGBP and immune cells. The levels of FCGBP in OSCC cells was examined through the use of RT-PCR. FCGBP overexpression was tested for its effects on OSCC cell proliferation and invasion using CCK-8 and Transwell assays. We observed that FCGBP expressions were distinctly downregulated in OSCC specimens compared with nontumor tissues in both TCGA and GEO datasets, which was further confirmed by RT-PCR. OSCC patients with advanced clinical stages and poor prognoses had lower levels of FCGBP expression. Many immune-related biological activities and signaling pathways were found to be considerably abundant in KEGG tests and GO analysis results. The correlation analysis indicated that FCGBP was associated with a number of immune cells in a positive way. We found that FCGBP expressions were strongly and distinctly linked to the expressions of known immunological checkpoints, and FCGBP expression had significant positive connections with tumor mutational burden. FCGBP upregulation distinctly slowed the growth and invasion of OSCC cells in functional experiments. FCGBP has the potential to be a therapeutic target for OSCC and a biomarker for OSCC patients' prognosis.

Long non-coding RNA DDX11-AS1 promotes the proliferation and migration of glioma cells by combining with HNRNPC

Glioma is a malignant tumor of the central nervous system with complex pathogenesis, difficult operation, and a high postoperative recurrence rate. At present, there is still a lack of effective treatment. Long non-coding RNA DDX11 antisense RNA 1 (DDX11-AS1) has been shown to promote tumor development, such as hepatocellular carcinoma, esophageal cancer, etc. However, its molecular mechanism in glioma is poorly understood. In this study, we found that the expression of DDX11-AS1 was elevated in glioma tissues, and patients with high expression of DDX11-AS1 had poor prognosis. DDX11-AS1 was a potential prognostic marker. Functionally, DDX11-AS1 promoted glioma cell proliferation and migration. Mechanistically, DDX11-AS1 interacted with RNA-binding protein heterogeneous nuclear ribonucleoprotein C (HNRNPC) to promote Wnt/β-catenin and AKT pathways and the epithelial-mesenchymal transition process. In summary, our study manifests that the DDX11-AS1/HNRNPC axis may play a vital part in the occurrence and development of glioma, which provides new ideas and therapeutic targets for the diagnosis, treatment, and prognosis of glioma.

[Comprehensive Analysis of the Relationship between m6A Methylation Patterns and Immune Microenvironment in Lung Adenocarcinoma]

BACKGROUND

m6A RNA methylation modification plays an important role in the occurrence and progression of lung cancer and regulates tumor immunity. Current studies mostly focus on the differential expression of some specific m6A effectors and infiltrating immune cell. m6A methylation modification is the result of mutual adjustment and balance between effectors, and changes in the expression of one or two effectors are far from enough to reflect the panorama of m6A methylation. The role of m6A in the immune microenvironment of lung adenocarcinoma (LUAD) is still poorly understood. The aim of this study is to investigate the effect of different m6A modification patterns in immune microenvironment of LUAD.

METHODS

LUAD data was obtained from The Cancer Genome Atlas (TCGA), University of California Santa Cruz Xena (UCSC Xena) and Gene Expression Omnibus (GEO) databases. Gene mutation, differential expression and survival analysis were performed for 24 m6A effectors. The m6A modification pattern was constructed by unsupervised clustering method, and the m6A clusters survival analysis, gene set variation analysis, immune score and immune cell infiltration analysis were performed. The association between LRPPRC protein expression levels and infiltration of CD8+ cytotoxic T lymphocytes and CD68+ macrophages in the tumor microenvironment was validated by immunohistochemistry in LUAD tissue microarray with 68 cases.

RESULTS

The mutations of m6A effector were found in 150 of 567 LUAD cases with a frequency of 26.46%. 6 readers and 3 writers were significantly up regulated in LUAD tissues compared with normal tissues. IGF2BP1 and HNRNPC are the independent risk factors for prognosis of LUAD. Abundant cross-talks among writers, erasers and readers were demonstrated. Three m6A modification patterns with different immune cell infiltration characteristics and clinical prognosis were established. Among m6A effectors, LRPPRC was found to be inversely associated with the infiltration of CD8+ cytotoxic T lymphocytes and CD68+ macrophages, and was validated in 68 LUAD tissues.

CONCLUSIONS

m6A modification patterns play non-negligible roles in regulating the immune microenvironment. LRPPRC has potential to be a new biomarker for checkpoint inhibitor immunotherapy.

Metabolism-Related Bioinformatics Analysis Reveals That HPRT1 Facilitates the Progression of Oral Squamous Cell Carcinoma In Vitro

Objectives

Many studies have shown that dysregulation of metabolism contributes to oncogenesis. However, the exact roles of metabolism-related genes (MRGs) in oral squamous cell carcinoma (OSCC) remain unclear. Thus, we aimed to identify a prognostic signature related to MRGs in OSCC.

Methods

The gene sequencing data of OSCC samples and the MRG set were downloaded from The Cancer Genome Atlas (TCGA) and the Molecular Signatures Database (MSigDB). The Wilcoxon rank-sum test was used to identify differentially expressed MRGs. Then, a prognostic signature was established by multivariate Cox regression analysis. Finally, prognosis-related MRGs were selected and further validated in OSCC tissues and cell lines.

Results

A prognostic signature that included 8 MRGs was constructed. Multiple survival analysis revealed that only HPRT1 might be an independent biomarker and indicator of poor overall survival in OSCC patients. The expression of HPRT1 was then found to be upregulated in OSCC tissues and cell lines, and suppression of HPRT1 gene expression by siRNA inhibited the proliferation, migration, and invasion of OSCC cells in vitro.

Conclusions

MRGs play an important role in the development of OSCC. Furthermore, HPRT1 might be an independent biomarker of OSCC and enhance OSCC proliferation, migration, and invasion in vitro; these results emphasize the potential utility of HPRT1 in OSCC therapy.

An analysis of the role of HnRNP C dysregulation in cancers

Heterogeneous nuclear ribonucleoproteins C (HnRNP C) is part of the hnRNP family of RNA-binding proteins. The relationship between hnRNP C and cancers has been extensively studied, and dysregulation of hnRNP C has been found in many cancers. According to existing public data, hnRNP C could promote the maturation of new heterogeneous nuclear RNAs (hnRNA s, also referred to as pre-mRNAs) into mRNAs and could stabilize mRNAs, controlling their translation. This paper reviews the regulation and dysregulation of hnRNP C in cancers. It interacts with some cancer genes and other biological molecules, such as microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and double-stranded RNAs (dsRNAs). Even directly binds to them. The effects of hnRNP C on biological processes such as alternative cleavage and polyadenylation (APA) and N6-methyladenosine (m6A) modification differ among cancers. Its main function is regulating stability and level of translation of cancer genes, and the hnRNP C is regarded as a candidate biomarker and might be valuable for prognosis evaluation.

Roles of Major RNA Adenosine Modifications in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common cancer malignancy worldwide and is known to have poor prognosis. The pathogenesis behind the development of HNSCC is not fully understood. Modifications on RNA are involved in many pathophysiological processes, such as tumor development and inflammation. Adenosine-related RNA modifications have shown to be linked to cancer and may play a role in cancer occurrence and development. To date, there are at least 170 different chemical RNA modifications that modify coding and non-coding RNAs (ncRNAs). These modifications affect RNA stability and transcription efficiency. In this review, we focus on the current understanding of the four major RNA adenosine modifications (N⁶-Methyladenosine, N¹-Methyladenosine, Alternative Polyadenylation Modification and A-to-I RNA editing) and their potential molecular mechanisms related to HNSCC development and progression. We also touch on how these RNA modifications affect treatment of HNSCCs.

m6A Regulator-Mediated Methylation Modification Patterns and Tumor Microenvironment Infiltration Characterization in Acute Myeloid Leukemia

Recent studies have demonstrated epigenetic regulation of immune responses. Nevertheless, the underlying effect of RNA N6-methyladenosine (m6A) modifications on tumor microenvironment cell infiltration remains elusive. In this study, we thoroughly assessed m6A modification patterns of 255 myeloid leukemia specimens based on 23 m6A regulators. Consensus clustering of the 23 m6A regulators was performed to determine three distinct m6A modification patterns that were remarkably consistent with three immunophenotypes of tumors: immunorejection, immune activation, and immune inertness. Further evaluation and prognostic analysis of the m6A modification patterns of individual tumors revealed that low m6A score was characterized by increased mutational burden, immune activation, and survival rates, whereas high m6A score was characterized by poorer survival rates and the absence of effective immune infiltration. In addition, this study investigated the association between m6A regulators and antitumor immune responses and discovered higher expression of the immune regulators PD-L1, PD-L2, MRP1, and MRP2 in low m6A scores. Generally, the expression pattern of m6A regulators was remarkably associated with prognostic results and antitumor immune responses in acute myeloid leukemia and may be an underlying target and biological marker for immune therapies.

Identification of Hypoxia-Related Molecular Classification and Associated Gene Signature in Oral Squamous Cell Carcinoma

The high heterogeneity of oral squamous cell carcinoma (OSCC) is the main obstacle for individualized treatment. Recognizing the characteristics of different subtypes and investigating the promising strategies for each subclass are of great significance in precise treatment. In this study, we systematically evaluated hypoxia-mediated patterns together with immune characteristics of 309 OSCC patients in the TCGA training set and 97 patients in the GSE41613 testing set. We further identified two different hypoxia subtypes with distinct immune microenvironment traits and provided treatment programs for the two subclasses. In order to assess hypoxia level individually, we finally constructed a hypoxia-related risk score, which could predict the clinical outcome and immunotherapy response of OSCC patients. In summary, the recognition of different hypoxia patterns and the establishment of hypoxia-related risk score might enhance our understanding of the tumor microenvironment of OSCC and provide more personalized treatment strategies in the future.

Use of 6 m6A-relevant lncRNA genes as prognostic markers of primary liver hepatocellular carcinoma based on TCGA database

BACKGROUND

Hepatocellular carcinoma (HCC) is diagnosed at the middle and advanced stages, negating radical treatment. Identifying specific and effective prognostic HCC biomarkers is important and can facilitate the discovery of potential therapeutic targets. N6-methyladenosine (m6A) and long non-coding RNAs (lncRNAs) are associated with the development of multiple tumors. The role of m6A-relevant lncRNAs in the initiation and progression of HCC is unclear. The aim of the present study was to investigate the expression of m6A-relevant lncRNAs in HCC and to identify new prognostic markers of the disease.

METHODS

Gene expression and clinical data were retrieved from The Cancer Genome Atlas database. m6A-relevant lncRNAs were identified by co-expression analysis and were screened by univariate Cox regression analysis. Different HCC patient clusters were established via consensus clustering. Gene set enrichment analysis (GSEA) was used to determine the cluster enrichment pathways. A risk score model was constructed, and Kaplan-Meier analysis of the overall survival (OS) between cluster 1 (high risk) and cluster 2 (low risk) was performed. Relationships between the clusters, risk scores, and clinicopathological characteristics were clarified.

RESULTS

Of the 1,852 m6A-relevant lncRNAs identified, 68 had prognostic relevance. The pathological grade, American Joint Committee on Cancer stage, and T stage of cluster 1 were significantly more advanced than those of cluster 2. Based on GSEA, mitotic spindle, G2M_CHECKPOINT, glycolysis, the phosphoinositide 3-kinase (PI3K) protein kinase B (AKT) mammalian target of rapamycin (mTOR) pathway, and DNA repair were more enriched in cluster 1. Six key m6A-relevant lncRNAs were selected to build a risk score model predicting the prognosis of HCC. The OS of patients in the high-risk group was shorter than that of patients in the low-risk group. Risk score was an independent prognostic factor of HCC patients.

CONCLUSIONS

The findings indicated that m6A-relevant lncRNAs may be important in the progression of HCC. The risk score model based on the 6 key m6A-relevant lncRNAs can accurately predict the prognosis of patients with HCC.

Research progress concerning m6A methylation and cancer

N6-methyladenosine (m⁶A) methylation is a type of methylation modification on RNA molecules, which was first discovered in 1974, and has become a hot topic in life science in recent years. m⁶A modification is an epigenetic regulation similar to DNA and histone modification and is dynamically reversible in mammalian cells. This chemical marker of RNA is produced by m⁶A 'writers' (methylase) and can be degraded by m⁶A 'erasers' (demethylase). Methylated reading protein is the 'reader', that can recognize the mRNA containing m⁶A and regulate the expression of downstream genes accordingly. m⁶A methylation is involved in all stages of the RNA life cycle, including RNA processing, nuclear export, translation and regulation of RNA degradation, indicating that m⁶A plays a crucial role in RNA metabolism. Recent studies have shown that m⁶A modification is a complicated regulatory network in different cell lines, tissues and spatio-temporal models, and m⁶A methylation is associated with the occurrence and development of tumors. The present review describes the regulatory mechanism and physiological functions of m⁶A methylation, and its research progress in several types of human tumor, to provide novel approaches for early diagnosis and targeted treatment of cancer.

HNRNPC Promotes Proliferation, Metastasis and Predicts Prognosis in Prostate Cancer

INTRODUCTION

The incidence of prostate cancer remains high worldwide, while exploring new therapeutic targets for prostate cancer is essential. Heterogeneous nuclear ribonucleoproteins have been proved to regulate tumorigeneses in various cancers. This study aimed to explore the role of HNRNPC in prostate cancer progression.

METHODS

HNRNPC expression and its correlation with clinical features and immune infiltration were analyzed by bioinformatics analysis. The effects of HNRNPC on prostate cell proliferation, migration, and invasion were accessed by EdU, colony formation, transwell, and wound-healing assays.

RESULTS

The expression level of HNRNPC was significantly increased in prostate cancer tissues and was correlated with the T stage, N stage, Gleason score, PSA level, residual tumors, overall survival, disease-specific survival, and progression-free interval of prostate cancer patients. Silencing HNRNPC inhibited the proliferation and metastasis of prostate cancer cells. The expression of HNRNPC was negatively correlated with the infiltration level of most immune cells in prostate cancer. Mechanistically, HNRNPC may function through regulating gene expression at the posttranscriptional level.

CONCLUSION

HNRNPC could be a potential marker for the treatment and prognosis prediction of prostate cancer.

METTL3 Intensifies the Progress of Oral Squamous Cell Carcinoma via Modulating the m6A Amount of PRMT5 and PD-L1

N6-Methyladenosine (m6A) modification is one of the commonest chemical modifications in eukaryotic mRNAs, which has essential effects on mRNA translation, splicing, and stability. Currently, there is a rising concern on the regulatory role of m6A in tumorigenesis. As a known component in the m6A methyltransferase complex, METTL3 (methyltransferase-like 3) plays an essential role in m6A methylation. Till now, the functions of METTL3 in oral squamous cell carcinoma (OSCC) and its relative mechanism remain to be explored. In this research, through the GEPIA database, we found that high METTL3 expression has a correlation with poor prognosis of squamous cell carcinoma of head and neck. qRT-PCR displayed that METTL3 was highly expressed in OSCC cells. Functionally, METTL3 knockdown reduced the invasion, migration, and proliferation competence of OSCC cells and attenuated the activation of CD8+ T cells. In contrast, METTL3 overexpression resulted in opposite results. GEPIA, UALCAN, and SRAMP databases, PCR, western blot, and m6A RNA methylation assay confirmed the m6A modification of PRMT5 and PD-L1 mediated by METTL3. In conclusion, our results displayed that METTL3 intensified the metastasis and proliferation of OSCC by modulating the m6A amounts of PRMT5 and PD-L1, suggesting that METTL3 may be a therapeutic target for OSCC patients.

ADAR and hnRNPC deficiency synergize in activating endogenous dsRNA-induced type I IFN responses

Cytosolic double-stranded RNA (dsRNA) initiates type I IFN responses. Endogenous retroelements, notably Alu elements, constitute a source of dsRNA. Adenosine-to-inosine (A-to-I) editing by ADAR induces mismatches in dsRNA and prevents recognition by MDA5 and autoinflammation. To identify additional endogenous dsRNA checkpoints, we conducted a candidate screen in THP-1 monocytes and found that hnRNPC and ADAR deficiency resulted in synergistic induction of MDA5-dependent IFN responses. RNA-seq analysis demonstrated dysregulation of Alu-containing introns in hnRNPC-deficient cells via utilization of unmasked cryptic splice sites, including introns containing ADAR-dependent A-to-I editing clusters. These putative MDA5 ligands showed reduced editing in the absence of ADAR, providing a plausible mechanism for the combined effects of hnRNPC and ADAR. This study contributes to our understanding of the control of repetitive element-induced autoinflammation and suggests that patients with hnRNPC-mutated tumors might maximally benefit from ADAR inhibition-based immunotherapy.

Genes involved in the epithelial-mesenchymal transition in oral cancer: A systematic review

OBJECTIVE

Epithelial-mesenchymal transition (EMT) is considered the initial step in the invasion-metastasis cascade. The aim of this systematic review was to study the signature of genes involved in the EMT process in oral cancer (OC) confirmed by protein expression and its possible relationship with oral squamous cell carcinoma (OSCC) prognostic variables.

MATERIALS AND METHODS

A search of the scientific literature was carried out with no start date restriction until 17 September 2020 in the electronic databases Pubmed/MEDLINE, Web of Science, Cochrane Library and Scopus, following specific eligibility criteria. The methodological quality of the included studies was assessed using the Newcastle-Ottawa tool.

RESULTS

A total of 8 retrospective cohort studies were included, all of them performed in China and with low risk of bias. Overexpression of the genes HNRNPC, ITGA5, HMGA2 and SRSF3, and low expression of ALDH3A1 and ARID2 promote EMT in OC. The more advanced clinical stages of the TNM classification were significantly associated with overexpression of HNRNPC, ITGA5, HMGA2 and SRSF3, and low expression of ARID2.

CONCLUSIONS

HNRNPC, ITGA5, HMGA2, SRSF3, ALDH3A1 and ARID2 genes were associated with EMT process. Over- or under-expression of these genes is associated with worse stages of OSCC and/or worse prognosis of the tumor. Further studies on this topic are needed in different countries to be able to confirm these results, since the detection of these genes can help to know which tumors have a worse prognosis.

Potential role of m6A RNA methylation regulators in osteosarcoma and its clinical prognostic value

BACKGROUND

Osteosarcoma is a disease with high mortality in children and adolescents, and metastasis is one of the important clinical features of osteosarcoma. N6-Methyladenosine (m6A) is the most abundant methylation modification in mRNA, which is regulated by m6A regulators. It is reported that it is related to the occurrence and development of tumors. However, the mechanism of its action in osteosarcoma is rarely known. The purpose of this study was to identify the potential role of m6A regulatory factor in osteosarcoma and its clinical prognostic value.

METHODS

Here, we used The Cancer Genome Atlas (TCGA) to comprehensively analyze the relationship between m6A regulatory factors and osteosarcoma (metastasis group and non-metastasis group). We analyzed their survival relationship and analyzed all the m6A regulatory factors in TCGA tumor data set by using the univariate Cox proportional hazard regression model. Finally, we selected two survival-related methylation regulators (FTO and IGF2BP2) as risk gene signature.

RESULTS

According to the median risk, patients were divided into low-risk group and high-risk group. Multivariate Cox regression analysis showed that these two risk genes were considered to be the key factors independently predicting the prognosis of patients with osteosarcoma. In addition, we verified their characteristics with gene expression omnibus (GEO) DataSets and confirmed that they are related to tumor and immune-related signaling pathways through gene set enrichment analysis (GESA) and immune infiltration analysis.

CONCLUSIONS

In conclusion, m6A regulators might play an important role in the metastasis of osteosarcoma and have potential important value for the prognosis and treatment strategy of osteosarcoma patients.

m6A RNA Methylation Regulators Impact Prognosis and Tumor Microenvironment in Renal Papillary Cell Carcinoma

Accumulating evidence has proven that N6-methyladenosine (m⁶A) RNA methylation plays an essential role in tumorigenesis. However, the significance of m⁶A RNA methylation modulators in the malignant progression of papillary renal cell carcinoma (PRCC) and their impact on prognosis has not been fully analyzed. The present research set out to explore the roles of 17 m⁶A RNA methylation regulators in tumor microenvironment (TME) of PRCC and identify the prognostic values of m⁶A RNA methylation regulators in patients afflicted by PRCC. We investigated the different expression patterns of the m⁶A RNA methylation regulators between PRCC tumor samples and normal tissues, and systematically explored the association of the expression patterns of these genes with TME cell-infiltrating characteristics. Additionally, we used LASSO regression to construct a risk signature based upon the m⁶A RNA methylation modulators. Two-gene prognostic risk model including IGF2BP3 and HNRNPC was constructed and could predict overall survival (OS) of PRCC patients from the Cancer Genome Atlas (TCGA) dataset. The prognostic signature-based risk score was identified as an independent prognostic indicator in Cox regression analysis. Moreover, we predicted the three most significant small molecule drugs that potentially inhibit PRCC. Taken together, our study revealed that m⁶A RNA methylation regulators might play a significant role in the initiation and progression of PRCC. The results might provide novel insight into exploration of m⁶A RNA modification in PRCC and provide essential guidance for therapeutic strategies.

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.

Bioinformatics Analyses Indicate That Cathepsin G (CTSG) is a Potential Immune-Related Biomarker in Oral Squamous Cell Carcinoma (OSCC)

Purpose

Plenty of studies showed that the immune system was associated with cancer initiation and progression. This study aimed to explore the prognostic biomarkers from immune-related genes (IRGs) in oral squamous cell carcinoma (OSCC).

Materials and Methods

RNA-seq data were downloaded from The Cancer Genome Atlas (TCGA) and IRGs and transcription factors (TFs) were extracted. Then, the co-expression network between IRGs and TFs was constructed using the "WGCNA" package in R software. Furthermore, a gene expression signature according to IRGs was constructed to predict OSCC prognosis and its accuracy was validated by survival analysis. Subsequently, correlation analyses between risk-score and immune cells level and clinical parameters were performed. Finally, immune-related biomarkers were selected and further investigated using gain-of-function assays in vitro.

Results

A total of 32 normal cases and 317 OSCC cases were selected in our study. Differentially-expressed analysis indicated that there were 381 differentially-expressed IRGs and 62 TFs in OSCC. Among them, 25 TFs and 21 IRGs were enrolled in the co-expression network. Furthermore, we found that gene expression signature on the basis of 10 IRGs could predict the prognosis accurately and a high-risk score based on gene expression signature meant a high T classification, terminal clinical stage, and low immune cells level in OSCC. Finally, cathepsin G (CTSG) was identified as a potential immune-related biomarker and therapeutic target in OSCC.

Conclusion

In conclusion, IRGs were directly involved in the development and progression of OSCC. Furthermore, CTSG was identified as a potential independent biomarker and might be an immunotherapeutic target in OSCC treatment.

Dual effects of N6-methyladenosine on cancer progression and immunotherapy

According to the latest global cancer statistics, cancer has become a major threat to human health, but cancer treatment has encountered many bottlenecks. As an emerging topic in epigenetics, N⁶-methyladenosine (m⁶A) is the most common internal modification on eukaryotic mRNA, which has attracted increasing attention in recent years. Accumulating studies have shown that aberrant m⁶A modifications have profound effects on the characteristics of tumors, which undoubtedly led to a significant breakthrough in cancer treatment. Although m⁶A function as an oncogene or tumor suppressor is not fully revealed, determining its precise function in the development and evolution of malignant tumors is crucial in improving clinical decisions involving targeted therapies. In this review, we briefly introduce the composition of the m⁶A methylation machinery and mainly summarize the biological mechanism of m⁶A in cancer cell death, angiogenesis, epithelial-mesenchymal transition (EMT), and therapeutic resistance. Subsequently, we present the exogenous regulatory factors of m⁶A and highlight the role of m⁶A on immune cells and cancer immunotherapy. The potential therapeutic strategies of m⁶A in human cancer are also discussed, considering research gaps and future applications.

RNA m6A Methylation Regulators Subclassify Luminal Subtype in Breast Cancer

RNA N⁶-methyladenosine (m⁶A) methylation is the most prevalent epitranscriptomic modification in mammals, with a complex and fine-tuning regulatory system. Recent studies have illuminated the potential of m⁶A regulators in clinical applications including diagnosis, therapeutics, and prognosis. Based on six datasets of breast cancer in The Cancer Genome Atlas (TCGA) database and two additional proteomic datasets, we provide a comprehensive view of all the known m⁶A regulators in their gene expression, copy number variations (CNVs), DNA methylation status, and protein levels in breast tumors and their association with prognosis. Among four breast cancer subtypes, basal-like subtype exhibits distinct expression and genomic alteration in m⁶A regulators from other subtypes. Accordingly, four representative regulators (IGF2BP2, IGF2BP3, YTHDC2, and RBM15) are identified as basal-like subtype-featured genes. Notably, luminal A/B samples are subclassified into two clusters based on the methylation status of those four genes. In line with its similarity to basal-like subtype, cluster1 shows upregulation in immune-related genes and cell adhesion molecules, as well as an increased number of tumor-infiltrating lymphocytes. Besides, cluster1 has worse disease-free and progression-free survival, especially among patients diagnosed with stage II and luminal B subtype. Together, this study highlights the potential functions of m⁶A regulators in the occurrence and malignancy progression of breast cancer. Given the heterogeneity within luminal subtype and high risk of recurrence and metastasis in a portion of patients, the prognostic stratification of luminal A/B subtypes utilizing basal-featured m⁶A regulators may help to improve the accuracy of diagnosis and therapeutics of breast cancer.

Identification of therapeutic targets and prognostic biomarkers from the hnRNP family in invasive breast carcinoma

Heterogeneous nuclear ribonucleoproteins (hnRNPs) are RNA-binding proteins that are reported to play a crucial role in the pathogenic process of multiple malignancies. However, their expression patterns, clinical application significance and prognostic values in invasive breast carcinoma (BRCA) remain unknown. In this study, we investigated hnRNP family members in BRCA using accumulated data from Oncomine 4.5, UALCAN Web portal and other available databases. We explored the expression and prognostic value level of hnRNPs in BRCA. We further analyzed their association with the clinicopathological features of BRCA patients. Subsequently, we calculated the alteration frequency of hnRNPs, constructed the interaction network of hnRNPs, and examined the potential coexpression genes of hnRNPs, revealing that HNRNPU and SYNCRIP are the core molecular genes requiring further investigation for BRCA. We validated the immunohistochemistry (IHC) pattern to simulate clinical applications based on pathology. Cell function experiments conducted in vitro indicated that HNRNPU can promote epithelial-mesenchymal transition, functionally stimulating the invasion capacity and inhibiting the viability of invasive BRCA cells. In summary, our systematic analysis demonstrated that HNRNPU was the key molecule that played a fundamental role in BRCA metastasis, which may facilitate the development of new diagnostic and prognostic markers for the analysis of BRCA progression.

Gene signatures of 6-methyladenine regulators in women with lung adenocarcinoma and development of a risk scoring system: a retrospective study using the cancer genome atlas database

Although the emergence of new treatments has improved the prognosis of women with lung adenocarcinoma (LUAD), the emergence of drug resistance limits their clinical efficacy. Therefore, there is an urgent need to identify new targets and develop a risk scoring system to evaluate the prognosis of patients. 6-methyladenine (M6A), as the most common methyl modification in RNA modification, its clinicopathological features, diagnosis and prognostic value in lung cancer, especially in LUAD remain to be discussed. We analyzed the clinical and sequencing data of the female LUAD cohort from The Cancer Genome Atlas (TCGA), evaluated the expression profiles of 16 M6A regulation-related genes in the cohort and the relationships between genetic changes and clinical characteristics, developed an M6A-related risk scoring system using Cox analysis. Finally, the copy number variations (CNVs) of the related genes in the samples were analyzed and verified using the cBioPortal platform. Compared with other clinical factors, this risk scoring system showed a higher predictive sensitivity and specificity. The M6A-related risk scoring system developed in this study may help to improve the screening of female patients at high risk of LUAD and provides important theoretical bioinformatics support for evaluating the prognosis of such patients.

RNA-Binding Proteins as Regulators of Migration, Invasion and Metastasis in Oral Squamous Cell Carcinoma

Nearly 7.5% of all human protein-coding genes have been assigned to the class of RNA-binding proteins (RBPs), and over the past decade, RBPs have been increasingly recognized as important regulators of molecular and cellular homeostasis. RBPs regulate the post-transcriptional processing of their target RNAs, i.e., alternative splicing, polyadenylation, stability and turnover, localization, or translation as well as editing and chemical modification, thereby tuning gene expression programs of diverse cellular processes such as cell survival and malignant spread. Importantly, metastases are the major cause of cancer-associated deaths in general, and particularly in oral cancers, which account for 2% of the global cancer mortality. However, the roles and architecture of RBPs and RBP-controlled expression networks during the diverse steps of the metastatic cascade are only incompletely understood. In this review, we will offer a brief overview about RBPs and their general contribution to post-transcriptional regulation of gene expression. Subsequently, we will highlight selected examples of RBPs that have been shown to play a role in oral cancer cell migration, invasion, and metastasis. Last but not least, we will present targeting strategies that have been developed to interfere with the function of some of these RBPs.