| Oral squamous cell carcinoma(OSCC)is one of the ten most common cancers in the world.There are about 300,000 new cases each year and the five-year survival rate is about 50%.Oral cancer is a broad term,including malignant tumors of the cheeks,tongue,base of the tongue base,gums,mouth floor,hard palate,etc.About 95%of oral malignant tumors are squamous cell carcinomas.The molecular mechanism of OSCC is largely unexplored.Therefore,it is necessary to explore the key molecules in the progression of oral squamous cell carcinoma,clarify its role in the progression of oral squamous cell carcinoma,and find some personalized molecular markers for the diagnosis of oral squamous cell carcinoma At the same time,these markers can help us to find therapeutic molecular targets,make prognostic judgments and personalize customized treatment plans for the patients with oral squamous cell carcinoma.At present,with the development of high-throughput sequencing technology,the cost of sequencing is lower,faster and more sensitive than traditional technology.Therefore,with the development of high-throughput sequencing data,it has promoted the development of sequencing data analysis methods and also provided a basis for further exploring the relationship between gene expression regulation in life sciencesThe purpose of this study is to explore the relationship between oral squamous cell carcinoma and gene expression through a series of bioinformatics analysis,screen out differentially expressed RNA,construct ceRNA networks,and explore potential biological markers.First,we download RNA-seq,miRNA-seq and clinical data of oral squamous cell carcinoma patients from the TCGA database and standardize,annotate and integrate them.RNA expression information was obtained from tumor samples of 327 oral squamous cell carcinoma patients and normal tissue samples of which 31 patients were paired.The R package "edgeR" was used to analyze RNA expression information to obtain RNA differentially expressed between tumor and normal tissues(|FC|>2,P<0.01).Then GO gene and KEGG pathway enrichment analysis of the differential genes were performed in the DAVID database.Next,we integrated the miRcode,miRDB,miRTarBase and TargetScan databases,and performed ceRNA network analysis on the differentially expressed RNA of 327 tumor samples through the R package "GDCRNATools"(P<0.01).We used Cytoscape network tool to visualize the ceRNA network,and uncoverd whether single factor COX and KM survival analysis is related to survival.Through the series of analysis,we obtained a total of 2203 differential expression mRNAs,170 differential expression miRNAs and 119 differential expression IncRNAs;Gene function and pathway enrichment results show that KEGG pathway is mainly enriched in focal adhesion,ECM-receptor interaction and cell cycle;The biological process of GO is mainly enriched in extracellular matrix organization,collagen catabolic process and cell adhesion.In the ceRNA network analysis,we obtained a ceRNA network consisting of 39 lncRNAs,31 miRNAs and 67 mRNAs.Among the results of single gene COX and KM survival analysis,the high expression of 6 RNAs was related to the overall survival time of tongue squamous cell carcinoma patients:CEP55,CCNA2,EGLN3,LINC00958,miR-130b and SLC16A1-AS1.The low expression of 1 RNA is related to the overall survival time of patients with tongue squamous cell carcinoma:STARD4-AS1.At the same time,LINC00958,miR-130b,CEP55 and EGLN3 form a prognostic ceRNA network.The analysis above show that 1 prognostic ceRNA network and 7 RNAs may play an important role in the development and metastasis of oral squamous cell carcinoma,and may be used as a biomarker that which is interrelated with the prognosis of OSCC patients.These results also provide a better understanding of the mechanism of occurrence,development and metastasis of tongue squamous cell carcinoma,and provide directions for future research. |