Animal model and bioinformatics analyses suggest the TIMP1/MMP9 axis as a potential biomarker in oral squamous cell carcinoma

Abstract Oral squamous cell carcinoma (OSCC) is a common malignant tumor of the head and neck. However, the molecular mechanism underlying its development and progression is yet unclear. Genes that are differentially expressed, that is, differentially expressed genes (DEGs), between normal and diseased tissues are believed to be involved in disease development and progression. To identify the DEGs in OSCC and explore their role in occurrence and progression, we established a Chinese hamster OSCC model, determined the DEG, screened the identified DEGs, and performed Gene Ontology (GO) and KEGG enrichment analyses. A protein–protein interaction (PPI) network was generated to screen potential candidate genes. We then analyzed the expression, tumor stage and prognosis of candidate genes using the Gene Expression Profiling Interactive Analysis (GEPIA) database. Finally, we verified the candidate DEGs by quantitative real‐time PCR and Gene Expression Omnibus analysis. The results showed 194 significantly DEGs, 140 enriched GO terms, and 8 KEGG pathways, which suggested that OSCC was closely related to the immune system, cell migration, and extracellular matrix. GEPIA and PPI network analysis revealed that SPP1, TNC, and ACTA1 were significantly related to tumor staging; SPP1, tissue inhibitors of matrix metallopeptidases (MMPs) 1 (TIMP1), and ACTA1 were closely related to prognosis. The scores for the top five highest degree genes were close, and the TIMP1/MMP9 axis appeared to be at the center of the PPI network, indicating that expression changes in the TIMP1/MMP9 axis and related genes may be involved in tumor invasion and metastasis. These findings provide novel insights into the mechanism of oral cancer.