Development of a murine model of oral carcinogenesis: an accelerated tool for biomarker and anti-tumour drug discovery

Sofia Ali Syed¹, Muhammad Asif Qureshi^2a, Saeed Khan², Rajesh Kumar³, Yusra Shafique⁴, Bilal Ahmed Khan² and Jawad Safdar⁵

¹Department of Oral Pathology, Dr. Ishrat-ul-Ibad Khan Institute of Oral Health Sciences, Dow University of Health Sciences, Karachi 74200, Pakistan

²Department of Pathology, Dow International Medical College, Dow University of Health Sciences, Karachi 74200, Pakistan

³Department of ENT, Dr. Ruth K.M. Pfau, Civil Hospital, Dow University of Health Sciences, Karachi 74200, Pakistan

⁴Department of Pathology, Dow Medical College, Dow University of Health Sciences, Karachi 74200, Pakistan

⁵Department of Oral & Maxillofacial Surgery, Dow Dental College, Dow University of Health Sciences, Karachi 74200, Pakistan

^ahttps://orcid.org/0000-0002-8370-1069

Abstract

Oral squamous cell carcinoma (OSCC) is the most common cancer in Pakistani men and the second most common cancer in women. The objective of our study was to devise a novel accelerated murine model of oral carcinogenesis that can be exploited as a tool to investigate the cancer circuitry involved in OSCC and to identify molecules of diagnostic, therapeutic and prognostic significance. A total of 40 healthy male, 6–8 weeks old, 22 ± 2 gram, Naval Medical Research Institute (NMRI) outbred strain mice were recruited in the experiment. NMRI mice are commonly used for animal experiments in various fields of biology and for drug toxicity. Of these, 25 mice underwent the oral carcinogenesis regimen via topical application of 0.5% 9,10-dimethyl-1,2-benzanthracene (DMBA) on the lower left lip for a maximum of 20 weeks and 15 mice were used as controls (without the carcinogenic regimen). Exophytic tissue masses were harvested, fixed in 10% formalin and stained with haematoxylin and eosin (H&E) for microscopic diagnosis. Additionally, the expression levels of CK 5/6, p53 and Ki-67 were investigated using immunohistochemistry. Of the 25 mice which underwent the carcinogenic regimen, 21 developed moderately differentiated squamous cell carcinoma and 1 showed dysplastic features with foci of invasion. Three mice were found dead with lesion(s). CK 5/6 showed strong positivity (100%) and p53 and Ki-67 showed patchy (<30%) strong positivity in OSCC, suggesting the similarity of our model to human OSCC. We present an accelerated, close-to-human carcinogenesis, model of oral carcinogenesis using DMBA in NMRI mice that can be exploited to study the pathogenesis of oral squamous cell carcinoma and subsequently devise immunotherapy or targeted therapy.

Keywords: carcinogenesis, 9,10-dimethyl-1,2-benzanthracene, mice, squamous cell carcinoma, tumour burden

Correspondence to: Muhammad Asif Qureshi
Email: email2asifqureshi@gmail.com

Published: 15/06/2022
Received: 23/01/2022

Publication costs for this article were supported by ecancer (UK Charity number 1176307).

Copyright: © the authors; licensee ecancermedicalscience. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Introduction

Lip and oral cavity cancer (C00-C06) is the most common cancer in Pakistani men and the second most common cancer in women [1–3]. Moreover, treatment options for these tumours are largely limited to conventional surgical, chemo/radiotherapeutic (cisplatin, fluorouracil) regimens. Despite a major tumour burden in Pakistan, this tumour type is under-researched in the country [4].

Oral carcinogenesis involves the accumulation of various genetic aberrations, including DNA mutation, epigenetic changes (DNA methylation) and chromosomal aberrations. More recently, NOTCH1, PIK3CA, AKT, PTEN and TRIM21 genomic changes have been identified as actionable genetic changes in oral cancers [5–7]. While chemotherapeutic regimens (cisplatin with or without fluorouracil and docetaxel) are most commonly used for the treatment of oral squamous cell carcinoma (OSCC), patients with prolonged chemotherapy develop a relapse due to chemoresistance [8]. Therefore, in order to identify relevant markers for early tumour detection and their response to chemotherapy and prognostic outcomes, it is highly relevant to undertake extensive research on this tumour type.

Diagnostic and therapeutic biomolecules can be investigated in vivo as well as in vitro. In-vivo models resemble more closely the complex biological environment and therefore are the mainstay of preclinical trials. For oral carcinogenesis, various murine models have been described and reviewed earlier, such as the 4-nitroquinolone 1-oxide (4NQO) model, the 9,10-dimethyl-1,2-benzanthracene (DMBA)-induced hamster cheek-pouch carcinogenesis model and the genetically modified mouse models of oral squamous cell carcinoma [9]. However, oral carcinogenesis models described in the literature did not objectively look into the number of days taken for mice to develop tumour(s). Moreover, the majority of the studies were performed using hamsters and rats. There is no anatomic counterpart of the buccal pouch in humans. Moreover, the buccal pouch epithelium is thinner in hamsters compared to mice and humans [10]. It, therefore, becomes highly relevant to modulate/tweak the existing carcinogenic regimens to develop/establish a novel close-to-human cancer murine model of accelerated oral carcinogenesis with high penetrance and replicability.

In the study presented herein, we describe an accelerated close-to-human cancer murine model of oral carcinogenesis that can be exploited not only to understand the mechanisms involved in oral carcinogenesis but also to investigate various molecular players of diagnostic, prognostic and therapeutic significance.

Methods

A total of 40 healthy immunocompetent males, 6–8 weeks old, 22 ± 2 gram, NMRI outbred strain of laboratory mice were purchased from the Department of Animal Sciences, Dow University of Health Sciences (DUHS). All mice (5 mice/cage) were transferred to the animal testing facility and were given a week’s acclimatisation to the facility with 26 ± 2°C temperature, 55% humidity and 12-hour light/dark cycles with food and water ad libitum. Mice were housed in standard murine cages and with standardised bedding, food and water supply. Bedding was changed either weekly or as required. Of the 40 mice, a total of 25 were subjected to the accelerated carcinogenic regimen, while 15 mice were used as controls. The study was approved by the Institutional Ethical Review Board for Animal Research and Ethics (Ref# AR.IRB-013/DUHS/Approval/2018/020). All animal experiments were performed at the animal testing facility of DUHS under the supervision of internationally licensed personnel trained and licensed as per the Animals (Scientific Procedures) Act 1986.

The induction of oral carcinogenesis was performed as previously published [11]. Briefly, oral carcinogenesis was induced in the lower left lip of the mice via topical application of 0.5% 9,10-dimethyl-1,2-benzanthracene (DMBA; Sigma-Aldrich). The 0.5% DMBA solution was prepared by diluting 100 mg of DMBA in 20 ml of acetone at room temperature. The DMBA application was performed twice every alternate day on the lower left side of the lip with a sable brush (no. 0) for a maximum 20 weeks. The mice were inspected daily. Additionally, the weights of the mice were recorded weekly. The mice were culled if any evidence of sickness was found or if the study end point was reached, whichever was earlier.

The collected tissues were fixed in 10% formalin, followed by tissue processing in an automated tissue processor, tissue embedding in paraffin wax moulds/blocks and sectioning at 4 µm thickness. Tissue sections were subsequently stained with haematoxylin and eosin (H&E). All slides were examined and analysed under a light microscope (Motic, BA310E; USA) by a consultant histopathologist (blinded for the samples) for identification of tumour cells. Images were captured using a camera (Motic Images Plus 3.0v; USA).

Expressions of cytokeratin (CK) 5/6, p53 and Ki67 were investigated using horseradish peroxidase (HRP, Dako) immunohistochemistry. Briefly, paraffin-embedded tissue sections were cut at 4 µm thickness using microtome. The sections were placed on pre-coated slides and incubated in a hot air oven at 80°C for 40 minutes. For antigen retrieval, slides were placed in a heat-induced epitope retrieval Dako PT Link containing target retrieval solution (pH 9) with low (65°C) to high temperature (97°C) for 20 minutes, followed by cooling at 65°C and rinsing with phosphate buffered saline (PBS) thrice at pH 7.4. Slides were then placed in a humidity chamber and incubated with mouse monoclonal primary antibodies against CK 5/6 (D5/16 Dako, Denmark), p53 (DO-7 Dako, Denmark) and Ki-67 (MIB-1 Dako, Denmark) for 30 minutes at room temperature, followed by rinsing with PBS thrice. Subsequently, slides were incubated with HRP for 30 minutes at room temperature. After washing with PBS thrice, slides were treated with ready to use (RTU) chromogen/substrate (1:20) for 5 minutes, followed by washing with distilled water. Finally, the slides were counterstained with haematoxylin for 1 minute, washed with distilled water, dried and mounted. Positive control was performed with human OSCC, appendix and ovarian carcinoma for CK 5/6, Ki 67 and p53, respectively. The slides were examined under a light microscope and evaluated by a histopathologist. Images were captured using a camera (Nikon; eclipse-80i; Japan). The expression of CK 5/6 was analysed in the cytoplasm and p53 and Ki-67 were analysed in the nucleus [12–14]. The expressions were evaluated by determining the percentages of positively stained cells (brown in colour) in 1,000 cells (magnification, x400) with staining intensity graded as 0 (negative), 1 (weak), 2 (moderate) and 3 (strong). Ki-67 and p53 positive cells were scored in percentages as 1 (<30%), 2 (30–60%), 3 (>60%) and CK 5/6 as 1%–10%, 11%–50% and more than 50%.

Data were entered and analysed using GraphPad Prism version 9.3 (GraphPad software, San Diego, CA, USA) by log-rank test. p-value ≤ 0.05 was considered statistically significant.

Results

Macroscopic changes

In all mice undergoing the carcinogenic regimen, the earliest macroscopic changes were observed after 2 weeks of DMBA administration. These included hair loss in the cutaneous part of the lower lip, followed by the appearance of ulcerative or exophytic lesions at different weeks post-DMBA administration. All the 25 mice (100%) which underwent the accelerated carcinogenic regimen developed tumours/lesions. The earliest lesion was observed on day 33, while the highest number of days taken for development of a lesion was 126. The average number of days taken by our mice to develop lesion was 84.56 days. Controls that were not administered DMBA did not show any morphological changes throughout the course of study (Figure 1).

Tumour burden

The overall average tumour burden was 4.64 lesions/mouse. The minimum number of lesions recorded in a mouse was 1 and the maximum number of lesions recorded was 13. For further analyses, lesions were divided into small, medium and large categories with ≤2 mm, 2.1–4 mm and >4 mm in diameter, respectively (Table 1).

At all times, a digital Vernier Caliper was used to measure the size of the lesions in two dimensions (length × width) in millimetres. Overall, the small category was the most commonly recorded lesion type. The smallest lesion was 0.87 mm in size, while the largest recorded lesion was 10 mm in diameter.

Histopathological analysis

Of the 25 mice, a total of 3 were found dead on day 120, day 135 and day 137, respectively. Since these mice were found dead, their tissues were not histopathologically analysed. Nevertheless, tissues from these dead mice were morphologically evident as carcinomas. Tissues collected from the remaining 22 mice were investigated for microscopic analysis using haematoxylin and eosin (H&E) staining. Of the 22 tissues investigated, 21 were moderately differentiated OSCC and 1 showed only dysplastic features, which was characterised by full dysplastic epithelium with foci of invasion, as shown in Figure 2.

Figure 1. Graphical representation of DMBA-treated and control mice. (A): Number of days taken to develop lesions in the control and DMBA-treated mice. (B): Quantification of grossly visible tumours in DMBA-treated mice. *p-value ≤ 0.05.

Immunohistochemical analysis

In order to confirm the epithelial origin of tumours, we investigated the expression of CK5/6 in murine tumours. Our data showed a significant expression of CK 5/6 in tumour tissues indicating epithelial origin of tumours developed due to DMBA treatment. The proliferative status of tumours was investigated using antibodies against Ki-67; a strong positive Ki-67 expression in tumours suggests high proliferative potential. We also investigated the p53 expression in cancer tissues and controls. Our data suggest high p53 expression in tumour cells compared to the controls, as shown in Figure 3.

Discussion

Lip and oral cavity cancer is the second most common malignancy, with a 5-year prevalence of 14.78% and 9.1% mortality in Pakistan [1]. Recently, regional cancer registries reported a high burden of lip and oral cavity cancer in Karachi [2,3]. In Pakistan (South Central Asia), the majority of patients present in the late stages of cancer and despite treatment modalities, the overall 5-year survival rate is below 50% [15].

It is therefore highly relevant to develop a cancer model that can be exploited for diagnostic and therapeutic purposes. Histopathologically, the majority of lip and oral cavity cancers are OSCC [16]. Tobacco, whether smoked or chewed, is the major risk factor for OSCC [17]. Therefore, it is of prime significance to induce tumours using a carcinogen which is found in tobacco smoke. The compound DMBA is a polycyclic aromatic hydrocarbon present in smoked/chewed tobacco and is widely used to induce tumours [18,19].

One of the major reasons for slow progression in oral cancer-related research could be the unavailability of appropriate tools to test various drugs in vivo. The available rodent models of oral carcinogenesis are time-consuming (32 weeks) and in vitro studies, on the other hand, have their own limitations in terms of cost, lack of histological nature and lack of heterogeneity of tumour [20,21]. It is therefore highly relevant to devise a novel in vivo model that does not only exhibit accelerated carcinogenesis but is also close to human oral carcinogenesis.

Table 1. Tumour burden in DMBA-treated mice.

Figure 2. Gross (left) and microscopic (right) images of DMBA-treated versus control mice. Representative images of haematoxylin and eosin (H&E)-stained OSCC and control tissue sections at ×40, ×100 and ×400 magnification.

With this context, we developed a murine model of accelerated oral carcinogenesis. We induced lip carcinogenesis via the topical application of DMBA. Several studies have reported the stages of chemically induced oral squamous cell carcinoma development [10,11]; however, the number of days taken to develop the first lesion and tumour size have not yet been reported. In our model, similar to human cancer, some mice developed lesion(s) earlier than others, as early as 5 weeks in 1 mouse while 15 (60%) mice developed exophytic lesions during 12 weeks and the 25 (100%) mice developed lesions by 18 weeks. Our findings are in contrast with previous study where red velvety rough appearance was observed after 4 weeks of DMBA but no macroscopic lesions (papilloma, ulcer and exophytic mass) were evident after 12 weeks. These differences could possibly be due to the differences in the animal strain and induction site, i.e., buccal pouch in hamsters, saliva and lymphocytic infiltrate that inhibited tumour development in a previous study model [22]. Another study revealed the first macroscopic change in tongue 19 weeks post-experiment and squamous cell carcinoma in 32 weeks in 4NQO-treated mice [20]. This could be explained based on the differences in application: 1 stroke of 4NQO in CBA mice, 3 times weekly for 16 weeks, as opposed to our regimen (detailed in the methodology).

In the carcinogenesis model described herein, the average time for onset of lesions was 12.08 weeks, which is in contrast to a previous study where the lesions developed in 21.1 ± 2.13 weeks [23]. The plausible reason for this difference could be the weight and age of mice. In our model, a young NMRI strain of mice, 6–8 weeks old, weighing 20 ± 2 g, were recruited compared to the previous study where Swiss mice, 14.2 weeks old, weighing 150 ± 30 g, were used.

Figure 3. Representative images of immunohistochemical and haematoxylin and eosin staining in DMBA-treated and control mice. High expression of CK 5/6, Ki-67 and p53 in OSCC compared to control tissue at ×40, ×100 and ×400 magnification.

We used Broder’s classification for grading OSCC [24]. Of 25 macroscopic lesions, 21 were diagnosed as moderately differentiated squamous cell carcinoma which was characterised by malignant neoplastic lesion arranged in nests and cords. Individual tumour cells showed individual cell keratinisation, vesicular nuclei, prominent nucleoli, abundant eosinophilic cytoplasm and keratin pearls. Stroma exhibited dense acute and chronic inflammation. These histopathological features mimic human OSCC [25] and one lesion showed dysplastic features. Three (12%) mice with tumours (macroscopic examination) were found dead; therefore, microscopic examination/diagnosis of their lesions was not performed.

In addition, immunohistochemistry was performed to further evaluate the similarity between mice and human OSCC. In our model, DMBA-treated mice demonstrated a strong positive (100%) expression of CK 5/6 in OSCC, indicating that the tumour was of epithelial origin. These findings are parallel with human OSCC [26]. TP53 is the most commonly mutated gene in human oral squamous cell carcinoma [27]. The mutated protein p53 accumulates in tumour cells and can be evaluated by immunohistochemistry. In our model, the expressions of p53 and Ki-67 were found to be strongly positive in DMBA-treated group compared to the control group; however, the positivity score was 1 (<30%). These findings are in accordance with other studies where the expressions of p53 and Ki-67 varied in both human and mice OSCC [28-31]. In summary, histopathological features and immunohistochemical expressions of CK 5/6, p53 and Ki-67 as described reveal that DMBA-treated NMRI mice developed tumours that closely mimic human OSCC. Nevertheless, the limitation of our study was that we did not evaluate the premalignant stages and lymph node metastasis.

Conclusion

Taken together, we provide evidence for the development of an accelerated immunocompetent murine model of oral carcinogenesis which is close-to-human oral carcinogenesis. This model can be exploited as an important tool to understand the molecular circuitry of oral carcinogenesis and to investigate various biomolecules of diagnostic, therapeutic and prognostic significance.

Conflicts of interest

The authors declare that they have no conflicts of interest.

Funding

The authors received no financial support for the research, authorship and/or publication of this article.

Authors’ contribution

SAS performed the experiment, acquired data, analysed and interpreted data and drafted the manuscript. MAQ conceived, designed and supervised the study, analysed and interpreted data and drafted the manuscript. SK, RK, BAK and JS facilitated sample collection and data acquisition. YS performed the histopathological analysis of the slides. All authors critically reviewed and approved the final draft and are responsible for the content and similarity index of the manuscript.

References

1. International Agency for Research on Cancer, World Health Organization [https://gco.iarc.fr/today/data/factsheets/populations/586-pakistan-fact-sheets.pdf] Date accessed: 08/08/21

2. Qureshi MA, Khan S, and Sharafat S, et al (2020) Common cancers in Karachi, Pakistan: 2010-2019 cancer data from the Dow cancer registry Pak J Med Sci 36 1572–1578 https://doi.org/10.12669/pjms.36.7.3056 PMID: 33235577 PMCID: 7674861

3. Pervez S, Jabbar AA, and Haider G, et al (2020) Karachi cancer registry (KCR): age-standardized incidence rate by age-group and gender in a mega city of Pakistan Asian Pac J Cancer Prev 21 3251–3258 https://doi.org/10.31557/APJCP.2020.21.11.3251 PMID: 33247682 PMCID: 8033128

4. Khan Z, Muller S, and Ahmed S, et al (2015) Quantitative review of oral cancer research output from Pakistan Asian Pac J Cancer Prev 16, 4733–4739 https://doi.org/10.7314/APJCP.2015.16.11.4733

5. Wu-Chou YH, Hsieh CH, and Liao CT, et al (2021) NOTCH1 mutations as prognostic marker in oral squamous cell carcinoma Pathol Res Pract 223, 153474 https://doi.org/10.1016/j.prp.2021.153474

6. Chuang CY, Chien YC, and Lin CW, et al (2021) TRIM21 polymorphisms are associated with susceptibility and clinical status of oral squamous cell carcinoma patients Int J Med Sci 18, 2997–3003 https://doi.org/10.7150/ijms.56614

7. Starzyńska A, Sejda A, and Adamska P, et al (2020) Prognostic value of the PIK3CA, AKT, and PTEN mutations in oral squamous cell carcinoma: literature review Arch Med Sci 17 207–217 https://doi.org/10.5114/aoms.2020.100780

8. Maji S, Shriwas O, and Saamal SK, et al (2019) STAT3- and GSK3β-mediated Mcl-1 regulation modulates TPF resistance in oral squamous cell carcinoma Carcinogenesis 40 173–183 https://doi.org/10.1093/carcin/bgy135

9. Ishida K, Tomita H, and Nakashima T, et al (2017) Current mouse models of oral squamous cell carcinoma: genetic and chemically induced models Oral Oncol 73 16–20 https://doi.org/10.1016/j.oraloncology.2017.07.028 PMID: 28939069

10. MacDonald DG (1981) Comparison of epithelial dysplasia in hamster cheek pouch carcinogenesis and human oral mucosa J Oral Pathol 10 186–191 https://doi.org/10.1111/j.1600-0714.1981.tb01264.x PMID: 6798180

11. Ribeiro DR, Alves ÂV, and dos Santos EP, et al (2015) Inhibition of DMBA-induced oral squamous cells carcinoma growth by Brazilian red propolis in Rodent Model Basic Clin Pharmacol Toxicol 117 85–95 https://doi.org/10.1111/bcpt.12374 PMID: 25556639

12. He W, Xiao Y, and Chen WM (2015) Expression of Ki-67 and P53 protein in oral squamous cell carcinoma and its clinical significance Shanghai Kou Qiang Yi Xue 24 228–231 PMID: 25938156

13. Kaufmann O, Fietze E, and Mengs J, et al (2001) Value of p63 and cytokeratin 5/6 as immunohistochemical markers for the differential diagnosis of poorly differentiated and undifferentiated carcinomas Am J Clin Pathol 116 823–830 https://doi.org/10.1309/21TW-2NDG-JRK4-PFJX

14. Prowse DM, Ktori EN, and Chandrasekaran D, et al (2008) Human papillomavirus-associated increase in p16INK4A expression in penile lichen sclerosus and squamous cell carcinoma Br J Dermatol 158 261–265 https://doi.org/10.1111/j.1365-2133.2007.08305.x PMCID: 2268980

15. Taghavi N and Yazdi I (2015) Prognostic factors of survival rate in oral squamous cell carcinoma: clinical, histologic, genetic and molecular concepts Arch Iran Med 18 314–319 PMID: 25959914

16. Dhanuthai K, Rojanawatsirivej S, and Thosaporn W, et al (2018) Oral cancer: a multicenter study Med Oral Patol Oral Cir Bucal 23 e23–e29 PMCID: 5822535

17. Jiang X, Wu J, and Wang J, et al (2019) Tobacco and oral squamous cell carcinoma: a review of carcinogenic pathways Tob Induc Dis 17 29 https://doi.org/10.18332/tid/111652 PMID: 31582940 PMCID: 6752112

18. Stepanov I, Villalta PW, and Knezevich A, et al (2010) Analysis of 23 polycyclic aromatic hydrocarbons in smokeless tobacco by gas chromatography-mass spectrometry Chem Res Toxicol 23 66–73 https://doi.org/10.1021/tx900281u PMCID: 2807893

19. Deng Y, Ma M, and Guo G, et al (2021) Kirenol regulates the cell proliferative and inflammatory markers in DMBA-induced oral squamous cell carcinogenesis in hamster Environ Toxicol 36 328–338 https://doi.org/10.1002/tox.23039

20. Schoop RA, Noteborn MH, and Baatenburg de Jong RJ (2009) A mouse model for oral squamous cell carcinoma J Mol Histol 40 177–181 https://doi.org/10.1007/s10735-009-9228-z PMID: 19685146 PMCID: 2770130

21. Tinhofer I, Braunholz D, and Klinghammer K (2020) Preclinical models of head and neck squamous cell carcinoma for a basic understanding of cancer biology and its translation into efficient therapies Cancers Head Neck 5 9 https://doi.org/10.1186/s41199-020-00056-4 PMID: 32714605 PMCID: 7376675

22. Martínez BDA, Barato Gómez PA, and Iregui Castro CA, et al (2020) DMBA-induced oral carcinoma in Syrian hamster: increased carcinogenic effect by dexamethasone coexposition Biomed Res Int 2020 1470868 https://doi.org/10.1155/2020/1470868

23. Alves ÂV, Ribeiro DR, and Lima SO, et al (2016) Expression of Ki-67 and P16 INK4a in chemically-induced perioral squamous cell carcinomas in mice Rev Col Bras Cir 43 72–79 https://doi.org/10.1590/0100-69912016002002 PMID: 27275587

24. Broders AC (1920) Squamous-cell epithelioma of the lip: a study of five hundred and thirty-seven cases J Am Med Assoc 74 656–664 https://doi.org/10.1001/jama.1920.02620100016007

25. Thompson L (2006) World Health Organization classification of tumours: pathology and genetics of head and neck tumours Ear Nose Throat J 85 74 https://doi.org/10.1177/014556130608500201 PMID: 16579185

26. Frohwitter G, Buerger H, and Van Diest PJ, et al (2016) Cytokeratin and protein expression patterns in squamous cell carcinoma of the oral cavity provide evidence for two distinct pathogenetic pathways Oncol Lett 12 107–113 https://doi.org/10.3892/ol.2016.4588 PMID: 27347109 PMCID: 4906805

27. Nakagaki T, Tamura M, and Kobashi K, et al (2017) Profiling cancer-related gene mutations in oral squamous cell carcinoma from Japanese patients by targeted amplicon sequencing Oncotarget 8 59113–59122 https://doi.org/10.18632/oncotarget.19262 PMID: 28938622 PMCID: 5601718

28. Agarwal VK, Sharma R, and Gahlot G, et al (2021) Clinical and histopathological correlation of p16 and p53 expression in oral cancer Indian J Surg Oncol 12 164–168 https://doi.org/10.1007/s13193-020-01145-z PMID: 33994742 PMCID: 8119508

29. Takkem A, Barakat C, and Zakaraia S, et al (2018) Ki-67 prognostic value in different histological grades of oral epithelial dysplasia and oral squamous cell carcinoma Asian Pac J Cancer Prev 19 3279–3286 https://doi.org/10.31557/APJCP.2018.19.11.3279 PMID: 30486632 PMCID: 6318382

30. Angiero F, Berenzi A, and Benetti A, et al (2008) Expression of p16, p53 and Ki-67 proteins in the progression of epithelial dysplasia of the oral cavity Anticancer Res 28 2535–2539 PMID: 19035275

31. Guo M, Mu Y, and Yu D, et al (2018) Comparison of the expression of TGF-β1, E-cadherin, N-cadherin, TP53, RB1CC1 and HIF-1α in oral squamous cell carcinoma and lymph node metastases of humans and mice Oncol Lett 15 639–645