Abstract
The typical progression of oral cancer is from hyperplastic epithelial lesions through dysplasia to invasive carcinoma. It is important to investigate malignant oral cancer progression and development in order to determine useful approaches of prevention of dysplastic lesions. The present study aimed to gain insights into the underlying molecular mechanism of oral carcinogenesis by establishing a rat model of oral carcinogenesis using 4‑nitroquinoline 1‑oxide. Subsequently, transcription profile analysis using an integrating microarray was performed. The dynamic gene expression changes of the six stages of rat oral carcinogenesis (normal, mild epithelial dysplasia, moderate dysplasia, severe dysplasia, carcinoma in situ and oral squamous cell carcinomas) were analyzed using component plane presentations (CPP)‑self‑organizing map (SOM). Six genes were verified by quantitative polymerase chain reaction, immunohistochemistry and succinate dehydrogenase (SDH) activity assay kit. Numerous differentially expressed genes (DEGs) were identified during rat oral carcinogenesis. CPP‑SOM determined that these DEGs were primarily enriched during cell cycle, apoptosis, inflammatory response and tricarboxylic acid cycle, indicating the coordinated regulation of molecular networks. In addition, the expression of specific DEGs, such as janus kinase 3, cyclin‑dependent kinase A‑1, B‑cell chronic lymphocytic leukaemia/lymphoma 2‑like 2, nuclear factor‑κB, tumor necrosis factor receptor superfamily member 1A, cyclin D1 and SDH were identified to have high concordance with the results from microarray data. The current study demonstrated that oral carcinogenesis is a multi‑step and multi‑gene process, with a distinct pattern alteration along a continuum of malignant transformation. In addition, this comprehensive investigation provided a theoretical basis for the understanding of the molecular alterations associated with oral carcinogenesis.