Background: M6A is the most widespread RNA internal modification that contributes to the transcription of RNA-mediated genes. Increasing evidence suggest that m6A has vital role in cancer initiation and progression. However, little is known about its roles in oral squamous cell carcinoma(OSCC). Therefore, this study aimed to explore the prognostic relevance of m6A related genes in OSCC through comprehensive bioinformatics analysis. Materials: Transcriptome expression profile were downloaded from The Cancer Genome Atlas (TCGA) following extracting m6A related genes and long non-coding RNA (lncRNA) based on expression level. The differentially expressed genes between OSCC samples and control samples was analyzed using R software. Furthermore, cluster analysis and principal component analysis (PCA) were used to analyze m6A related genes, and a risk model was constructed on basis of m6A related genes. In addition, co-expression network was established according to m6A related genes and differentially expressed lncRNA (DElncRNA) expression profile and crucial m6A genes were selected based on univariate, multivariate cox regression and survival analysis. Then, CCK-8, colony formation assays, transwell assay and scratch wound healing were performed to explore the biological effects of identified m6A-related genes in OSCC. Furthermore, Western blot assay was conducted to detect the EMT-associated genes expression level. Findings: A total of 13 m6A related genes were extracted and 8 genes were differentially expressed in OSCC. Subsequently, 418 DElncRNA were selected. M6A-Based Clustering Shows 2 Subtypes of OSCC with different clinical outcome.Then, a risk model including HNRNPC, METTL14, YTHDF2, ALKBH5 genes was established on basis of m6A expression profile, which can indicate the survival time and differentiation grade and may be an independent prognostic biomarkers. Furthermore, a co-expression network including 7 m6A related genes and 21 DElncRNA was constructed based on DElncRNA and m6A related genes and HNRNPC was identified as an independent prognostic biomarker in OSCC. Functional studies revealed that HNRNPC knockdown inhibited cell proliferation, migration and invasion of OSCC cells in vitro. Conversely, HNRNPC overexpression showed the opposite biological effects. Moreover, we identified that HNRNPC regulated epithelial-mesenchymal transition in OSCC. Interpretation: This study revealed that m6A may play a significant role in OSCC initiation and progression through integrated bioinformatics analysis. And co-expression network may be helpful for revealing potential correlation of lncRNA and m6A in OSCC. In addition, a risk model was established to predict the prognosis of OSCC.At last, HNRNPC was proved to promote OSCC carcinogenisis and be an independent prognostic biomarker of OSCC, suggesting that it may be a new therapeutic target of OSCC. Funding Statement:This study was supported by the National Natural Science Foundation of China (81472536); the Science and Technology Planning Project of Guangdong Province (No. 2017A020215181 and No. 2014A020212440); Project of Educational Commission of Guangdong Province of China (2018KTSCX026); the Scientific Research Planning Project of Southern Medical University (CX2018N016, PY2018N031) and the Presidential Foundation of the Nanfang Hospital (2014027, 2019Z030). Declaration of Interests: All authors declare no conflicts of interest. Ethics Approval Statement: The study protocol was approved by the Ethics Committees of Nanfang Hospital of Guangdong Province (NO: NFEC-2018-027).
Read full abstract