Identification of the key genes implicated in the transformation of OLP to OSCC using RNA-sequencing

Oral lichen planus (OLP) is a chronic inflammatory disease that may transform to oral squamous cell carcinoma (OSCC), while its carcinogenesis mechanisms are not entirely clear. This study was designed to identify the important genes involved in the malignant transformation of OLP to OSCC. After RNA-sequencing, the differently expressed genes (DEGs) in OLP vs. normal and OSCC vs. normal groups, respectively, were identified by limma package in R language, and then clustering analysis were conducted by Pheatmap package in R language. Weighed gene co-expression network analysis (WGCNA) was performed for the DEGs to screen disease-associated modules. Using Cytoscape software, co-expression networks were constructed for the genes involved in the modules. Enrichment analysis was conducted for the genes involved in the co-expression networks using GOstat package in R language. Finally, quantitative real-time PCR (qRT-PCR) experiments were conducted to validate the key genes. There were, respectively, 223 and 548 DEGs in OLP vs. normal and OSCC vs. normal groups. WGCNA identified the blue modules for the DEGs in the two groups as disease-associated modules. Moreover, 19 common DEGs (including upregulated BCL9L, PER2 and TSPAN33, and downregulated GMPS and HES1) associated with both OLP and OSCC were identified. In the co-expression networks, BCL9L, HES1, PER2 and TSPAN33 might function in OLP via interactions (such as BCL9L-TSPAN33 and HES1-PER2). qRT-PCR analysis showed that BCL9L, PER2 and TSPAN33 were significantly upregulated, and GMP and HES1 were downregulated. These findings indicated that BCL9L, GMPS, HES1, PER2 and TSPAN33 affected the transformation of OLP to OSCC.