Quantitative analysis of KLF4 and SOX2 expression in oral carcinomas reveals independent association with oral tongue subsite location and histological grade

Abstract

BACKGROUND:

Stemness factors associated with tumorigenesis in different types of cancers have not been specifically studied in oral tongue SCC (OTSSC). Here, we aimed to quantify expression levels and distribution of KLF4 and SOX2, two relevant stemness factors, in oral SCC including OTSCC samples from different subsites.

METHODS AND RESULTS:

We determined KLF4 and SOX2 expression levels by immunostaining 35 biopsies of OSCC. Stained wholeslide images were digitized and subjected to automatic cell detection and unbiased quantification using Qupath software. We found statistically significant reduction in KLF4 positive cells density ( $p=$ 0.024), and fraction ( $p=$ 0.022) in OTSCC from tongue borders compared with other tongue subsites. Instead, quantitative SOX2 analysis did not show differences in expression levels between OTSCC from the borders versus OTSCC developed in others subsites. Notably SOX2 expression was revealed increased in moderately and poorly differentiated OSCC compared with well differentiated ones (positive cells density $p=$ 0.025, fraction $p=$ 0.006). No significant correlation between KLF4 and SOX2 expression was observed, neither in OSCC nor in OTSCC.

CONCLUSIONS:

KLF4 and SOX2 exhibit opposite expression profiles regarding subsite localization and differentiation level in OSCC. Our study prompts future OTSCC prospective studies looking for clinical prognosis to incorporate detailed subsite information in the analysis.

Keywords

KLF4 SOX2 oral squamous cell carcinoma oral tongue cancer whole slide scanning quantitative image analysis

1. Introduction

Head and neck cancers (HNC), including cancers of the oral cavity, are highly prevalent worldwide, characterized by high morbidity and limited therapeutic options [1, 2]. Oral cancers show variations in aetiology, incidence rates, clinical progression and prognosis based on the anatomical subsite affected [3, 4, 5]. The most common histologic type is squamous cell carcinoma (SCC) comprising over 90% of oral cancers [6, 7]. In this regard, oral (mobile) tongue squamous cell carcinoma (OTSCC) is one of the most common malignant tumors of the oral cavity [8, 9]. Moreover, current epidemiological records show that oral tongue cancers are becoming more common amongst the younger patients than in previous decades, although the exact reason of this trend is yet to be elucidated [8, 9]. OTSCC is generally associated with poor prognosis [9, 10] and behaves aggressively with higher rates of occult metastasis and neck nodal metastasis than any other cancer of the oral cavity [10, 11]. Patients with OTSCC have a significantly worse prognosis than those with similar lesions of the oropharynx, larynx, hypopharynx and others oral cavity sites [11]. These clinical and epidemiological information highlight the importance of developing more effective diagnostic modalities associated with different oral cancer subsites that allow for targeted therapeutic approaches.

In the last years, new research on molecular mechanism underlying HNC development has revolutionized the treatment for patients with metastatic disease. In particular, immune modulation with checkpoint inhibitors has been established and applied to HNC patients [12, 13]. However, the overall survival of treated patients after development of metastatic disease has barely improved [14], in part due to treatment resistance. The latter is thought to be related to intratumoral heterogeneity as well as the presence of specific population of resistant cancer stem cells (CSC) [15, 16]. In this regard, previous studies have shown that the expression of the Yamanaka factors (Oct3/4, Sox2, Klf4, c-Myc) is sufficient to reprogram differentiated cells to pluripotent stem cells [17]. One of these factors, the transcription factor Krppel-like factor 4 (KLF4), in addition to its role in cell stemness, suppresses cell proliferation, promotes differentiation, and thus helps to maintain homeostasis of adult epithelial cells [18, 19]. Furthermore, KLF4 has a tumor suppressor function in many human malignancies including gastric, colorectal and lung cancers [19, 20]. Nevertheless KLF4 was also reported to have an oncogenic role in breast and other cancers in a context dependent-manner [18, 20]. However, the role of KLF4 in head and neck cancer remains controversial. Increased KLF4 expression has been reported in head and neck SCC by Shyh-Kuan Tai et al. [21]. Conversely, KLF4 expression was shown decreased in oral SCC (OSCC) [22]. Furthermore, sex-determining region Y-box 2 (SOX2) is other Yamanaka factor that has been suggested as a marker of CSC [23, 24]. SOX2 is a transcription factor expressed in a wide range of stem and progenitor cells in diverse adult tissues, including squamous epithelia of the tongue and genetically amplified in different human cancers [24, 25]. However, to date, the expression pattern of KLF4 and SOX2 has not been specifically studied in OTSSC. We hypostatize that OTSCC may show a differential expression of these Yamanaka factors connected with tumor localization and histologic grade. In the present study we quantified KLF4 and SOX2 expression levels by immunohistochemical staining in OSCC biopsies including OTSCC samples from different oral tongue subsites.

2. Materials and methods

Thirty five excisional specimens with histological diagnosis of squamous cell carcinoma were identified retrospectively at the Surgical Pathology Laboratory, Oral Pathology Department of the School of Dentistry of the University of Buenos Aires, Argentina. Clinicopathological information, including age, gender and location was obtained from patient records (Table 1). Tumors were graded using World Health Organization criteria [26]. The protocol for this study was approved by the Ethics Committee of the School of Dentistry of the University of Buenos Aires.

Table 1
Clinicopathologic characteristic for 35 selected cases of OSCC

Case	Age of diagnosis (years)	Sex	Site	Histological grade
1	82	F	TONGUE (B)	WD
2	67	M	TONGUE (B)	WD
3	84	M	TONGUE	MD
4	61	F	TONGUE	MD
5	46	M	TONGUE (B)	MD
6	56	M	TONGUE	MD
7	63	M	TONGUE (B)	MD
8	73	F	TONGUE (B)	MD
9	68	F	TONGUE	MD
10	61	F	TONGUE (B)	MD
11	69	M	TONGUE	MD
12	78	M	LIP	MD
13	45	M	TONGUE (B)	MD
14	65	M	TONGUE (B)	MD
15	80	F	CHEEK	WD
16	78	F	PALATE	WD
17	76	F	ALVEOLAR RIDGE	WD
18	79	F	ALVEOLAR RIDGE	WD
19	76	M	CHEEK	WD
20	78	M	GINGIVA	WD
21	72	M	LIP	WD
22	80	F	GINGIVA	WD
23	69	M	TONGUE	PD
24	73	M	ALVEOLAR RIDGE	PD
25	60	M	GINGIVA	PD
26	70	F	ALVEOLAR RIDGE	PD
27	55	F	GINGIVA	PD
28	50	F	TONGUE (B)	PD
29	61	M	TONGUE (B)	MD
30	65	F	TONGUE (B)	WD
31	78	F	TONGUE	MD
32	77	F	TONGUE	MD
33	59	M	TONGUE (B)	MD
34	56	M	TONGUE	MD
35	78	F	GINGIVA	MD

B, lateral border of the tongue; F, female; M, male; WD, well differentiated; MD, moderately differentiated; PD, poorly differentiated.

2.1 Histopathology and immunohistochemistry

Before immunohistochemical studies, H&E-stained section from each specimen was evaluated and the suitability of inclusion for the study was determined. Five micrometer serial sections were used for all subsequent analysis. The tissues slides were dewaxed and hydrated through graded alcohols and incubated in 3% hydrogen peroxide in PBS for 30 min to quench the endogenous peroxidase. Antigen retrieval was done using 10 mmol/L citric acid in a microwave and the following steps were performed as described previously [27]. Primary antibodies used were: KLF4 (H180) at dilution 1:400 and Sox-2 (E4) at dilution 1:50 from Santa Cruz Biotechnology. The slides were washed and developed in 3, 3-diaminobenzidine (DAB) (FASTDAB tablet; Sigma) under microscopic control and after counterstained with Mayer’s hematoxylin.

2.2 Qualitative staining evaluation

Initially, tissue sections immunohistochemically stained for KLF4 and SOX2 were independently classified by two pathologists blinded to the clinicopathological information. A case was marked negative when both observers considered that nuclear DAB labeling was absent or only present within a few limited tumor cells (less than $\sim$ 1%).

2.3 Image acquisition and analysis

Stained histological slides were scanned and digitized for unbiased cell detection and automatic quantification using Qupath software (https://qupath.githubio) [28]. Stained histological slides were scanned using a fully motorized AxioImager.M2 microscope (Zeiss) controlled by Stereoinvestigator software (MBF Bioscience). Brightfield images were automatically acquired under 20x magnification and stitched together into a whole-slide composite image. Images were subsequently exported as high-resolution TIFF images (final resolution 0.322 $\mu$ m/pixel) and off-line analyzed using Qupath 0.1.2. Regions of interest (ROIs) were manually defined and annotated (by a trained pathologist) in each whole-slide image including all tumor areas present in the image. Areas with cutting or staining artifacts were excluded from the ROIs (Fig. 1A).

Figure 1.

Automatic cell detection with mean DAB nuclear intensity measurement in ROIs using Qupath software. A. Low magnification whole slide composite image depicting the ROIs drawn by a trained pathologist in a tongue case immunostained for KLF4. B. Same image with the superimposed color-coded cell classification obtained after Qupath analysis (red positive, blue negative). C. Digital zoom within a ROI of the whole slide composite image depicted in A showing cells with different degrees of nuclear DAB staining intensity. D. Frequency histogram of DAB nuclear intensity values obtained for all detected cells in the ROIs of the case presented in A. The arrow indicates positive cell threshold in this case. Cells with DAB nuclear signal above threshold are considered positive (black). Inset: schematic representation of threshold value determination for case A, frequency histogram of DAB nuclear signal of case A lymphocytes. Red bar indicates mean $+$ / $-$ SD values for lymphocyte population. E. Same image depicted in C with cell and nuclear boundaries delineated by Qupath automatic analysis. F. Same image with color coded overlay indicating positive and negative cells classified using threshold calculated in D. Note that cells were classified according to DAB signal intensity without interference of Hematoxylin signal.

Stain separation (Hematoxylin/DAB) was achieved with Qupath color deconvolution tool using each image RGB modal values for automatic classification. Next, automatic cell detection was achieved using optical density sum. Detection parameters initially optimized by visual inspection and kept constant across images with few exceptions. Qupath output parameters were collected for all individual cells included in the ROIs. No filtering or exclusion criteria were adopted (Fig. 1B).

To establish an observer-independent threshold for positive cell detection (allowing automatic cell classification as positive/negative) we quantified DAB nuclear staining in a subset of 50–300 lymphocytes per slide. Nuclear lymphocytes DAB intensity followed a normal/quasi-normal distribution (Fig. 1D) reflecting a homogenous population (Fig. 1D). Threshold was independently set for each case as the mean nuclear DAB signal $+$ one standard deviation of their lymphocyte population (Fig. 1D). Cells were considered positive when nuclear DAB signal exceeded the case threshold value. When frequency distribution of nuclear DAB signal from tumor cells exhibited a bimodal distribution threshold values coincides with the histogram valley (Fig. 1D). Lymphocytes were selected for stain normalization because they are easy to identify, they are present in all cases (inside or outside the ROIs), display a low level of KLF4 and SOX2 immunoreactivity comparable with background staining levels in non-cellular regions of the sample [29, 30].

2.4 Quantification

Positive cells fraction was defined as the number of positive tumor cells divided by total number of detected cells, positive cells density was calculated as the number of positive tumor cells divided by ROIs combined area (positive cells/mm ${}^{2}$ ). Nuclear labeling intensity corresponds to the mean value of the DAB stain for the nuclear region of the positive tumor cel. A single value was obtained for each case by averaging all individual measured cells from all ROIs annotated in the case. Population data is presented as mean $\pm$ SD of individual cases. KLF4 and SOX2 staining were independently analyzed.

2.5 The Cancer Genome Atlas (TCGA) data

TCGA mRNA expression analysis was conducted using a publicly available dataset. Patient data from TCGA (Head and Neck Squamous Cell Carcinoma. Source data from GDAC Firehose. Previously known as TCGA Provisional), including site location, were retrieved using cBioportal web site (https://www.cbioportal.org). All samples from selected locations were included irrespectively of their HPV status.

2.6 Statistical analysis

Data are presented as mean $\pm$ SD. Student’s $t$ -test was used for comparison of two independent groups, comparison of three or more groups was done by one way ANOVA followed by Newman-Keuls post hoc test when required. CV differences between groups were assessed by Shapiro-Wilk test using Real Statistics Resource Pack software (Release 6.8), other statistical analyses were performed using Statistica 8.0. Association between two variables was analyzed by Pearson correlation. All tests were conducted two-side and considered statistically significant at $p<$ 0.05.

Figure 2.

Assessment of KLF4 immunohistochemical expression in oral squamous cell carcinomas. A. Normal representative oral epithelium with positive cells with KLF4 expression in the intermediate layer and strongest staining at upper layers. B. Representative image of SCC developed in the cheek with clear nuclear staining for KLF4 throughout the malignant epithelial cells. C. A representative carcinoma from gingiva immunostained for KLF4 detection displays positive cells with a wide range of intensity labeling. D. OTSCC developed in the dorsal tongue shows KLF4 positive staining with variable intensity. E. OTSCC from the lateral border of the tongue positive for KLF4 staining. F, G. Detail at higher magnification of the boxed areas in E. Images F and G depict two different areas with positive KLF4 expression exemplifying regional differences in fraction and density of KLF4 positive cells.

Figure 3.

Quantitative analysis of KLF4 expression in OTSCC versus OSCC. A. Fraction, density and nuclear labeling intensity of KLF4 positive cells for OTSCC (T) and the rest of the SCC developed in the oral cavity (NT, non-tongue). Bars represent mean $+$ SD. Non-significant differences were observed by Student $t$ -tests (fraction: t ${}_{\text{33}}=$ 0.83 $p=$ 0.41; density: t ${}_{\text{33}}=$ 1.58 $p=$ 0.1; intensity: t ${}_{\text{33}}=$ 1.01 $p=$ 0.32). B. Individual values of fraction, density and nuclear labeling intensity of KLF4 positive cells for each case according to the oral cavity subsite (the subsites are color coded as shown in graph legend). Coefficient of variation (CV, in blue) was calculated for each group to assess differences in data dispersion. No significant differences were observed between groups for fraction ( $p=$ 0.15), density ( $p=$ 0.42), and intensity ( $p=$ 0.39) after Shapiro-Wilk test. A.R., alveolar ridge. Each dot represents one individual biopsy from the data set.

Figure 4.

Quantitative analysis of KLF4 expression in OTSCC developed in different tongue subsites. Fraction (A), density (B) and intensity (C) of KLF4 positive cells for OTSCC from the borders of the tongue (tongue border) and the rest of the oral tongue subsites (other tongue site). Bars represent mean $+$ SD. Dots represent values of individual cases. ${}^{*}p<$ 0.05, Student $t$ -test, A: t ${}_{{19}}=$ 2.50 $p=$ 0.02; B: t ${}_{{19}}=$ 2.45 $p=$ 0.024; C: t ${}_{{19}}=$ 0.032 $p=$ 0.90.

Figure 5.

Histologic degree of OSCC and KLF4 expression. Fraction (A), density (B) and intensity (C) of KLF4 positive cells for OSCC grouped according to histologic degree WD, well differentiated; MD, moderately differentiated; PD, poorly differentiated. Bars represent mean $+$ SD. Dots represent values of individual cases. Non-significant one-way ANOVAs, A: F ${}_{2,32}=$ 0.39 $p=$ 0.6; B: F ${}_{2,32}=$ 3.27 $p=$ 0.051; C: F ${}_{2,32}=$ 0.0004 $p=$ 0.97.

3. Results

With the aim to characterize KLF4 protein expression in tongue carcinomas and its relation to oral cancer we performed immunohistochemistry in OSCC biopsies from different oral locations (Table 1). Nuclear staining for KLF4 in the intermediate and upper layers of the normal oral mucosa epithelium was observed (Fig. 2A). KLF4 protein was detected at variable levels and clearly stained in the cell nucleus throughout the epithelium in 94.3% (33 out of 35) of the OSCC analyzed by qualitative evaluation (Fig. 2). However the pattern and intensity of expression for KLF4 in OSCCs was very diverse. We described positive cases displaying regions with near full cell labelling intermingled with areas almost lacking of KLF4 staining. Therefore we performed whole slide scanning in all stained biopsies in order to do quantitative analysis on digital images using Qupath software. On average, a ROI area equivalent to 4.72 mm ${}^{2}$ was analyzed in each case with 4062 to 18403 tumor cells detected and analyzed. No statistically significant differences were observed regarding KLF4 expression when comparing OTSCC versus SCCs from others oral cavity anatomical locations neither in fraction or density of KLF4 positive cells nor in KLF4 nuclear labeling intensity (Fig. 3A). Surprisingly, the variability of positive cells fraction and density within tongue samples was similar to the one observed for all others oral anatomical locations combined together (Fig. 3B). These results suggest the presence of regional heterogeneity within oral tongue subsites. To extend our results we examined KLF4 expression at mRNA level by using a public database containing over 500 HNSCC cases. Supplementary Fig. 1A shows that the regional heterogeneity of KLF4 protein expression observed immunohistochemically is also present at mRNA level. Since OTSCC behaves aggressively [11] and considering the lack of available information regarding oncogene or tumor suppressor genes expression in OTSCC from different subsites we next evaluated if there is a differential KLF4 expression in different areas of the oral tongue.

We analyzed KLF4 immunostaining in OTSCC developed in the borders or in the rest of the tongue (Fig. 4). We found statistically significant differences in KLF4 positive cells fraction in OTSCC from the borders of the tongue when compared with OTSCC from the rest of the tongue (Fig. 4A). The same results were found when we calculated the density of KLF4 positive cells (Fig. 4B). Although we did not observed changes in the KLF4 nuclear labeling intensity (Fig. 4C). The lack of subsite categorization in the available TCGA databases precluded the comparative assessment of KLF4 mRNA levels whiting the tongue. Considering that KLF4 has been reported as a differentiation marker for skin and oral epithelia as well as OSCC [31, 32] we analyzed a possible relation between the histologic degree and KLF4 expression in our group of cases. We did not find any association between the grade of differentiation of the carcinomas and KLF4 expression (Fig. 5). Furthermore no association was found between differentiation degree and KLF4 expression pattern when the analysis was restricted to oral tongue carcinomas (one-way ANOVAs, cell fraction: $F_{2,18}=$ 0.23 $p=$ 0.79, cell density: $F_{2,18}=$ 0.83 $p=$ 0.45). These results reinforce the notion that in oral cavity carcinomas KLF4 expression varies among location rather than differentiation.

Figure 6.

SOX2 immunohistochemical expression in oral squamous cell carcinomas. Representative sections of WD, MD and PD oral carcinomas with SOX2 expression. A, B. Image from a SCC developed in the gingiva and with low expression of SOX2. SOX2 immunostaining was localized in the nucleus of the SCC cells (B). C, D. Image from a carcinoma from the lateral border of the tongue presenting a moderate to strong expression of SOX2 (D). E, F. Representative image from a SCC developed in the gingiva showing strong SOX2 expression throughout the parenchyma of the carcinoma (F). Insets in A, C and E show high power view from each carcinoma with different histologic degree studied.

Figure 7.

Histologic degree of OSCC and SOX2 expression. Fraction (A), density (B) and intensity (C) of SOX2 positive cells for OSCC grouped according to the histologic degree WD, well differentiated; MD, moderately differentiated; PD, poorly differentiated. Bars represent mean $+$ SD. Dots represent values of individual cases. ${}^{*}p<$ 0.05, Newman-Keuls post hoc test after significant one-way ANOV; A: F ${}_{2,32}=$ 4.14 $p=$ 0.025, B: F ${}_{2,32}=$ 5.92 $p=$ 0.006. C: non-significant one-way ANOVA F ${}_{2,32}=$ 0.55 $p=$ 0.58.

Figure 8.

Quantitative analysis of SOX2 expression in OTSCC developed in different tongue subsites. Fraction (A), density (B) and intensity (C) of SOX2 positive cells for OTSCC from the borders of the tongue (tongue border) and the rest of the oral tongue subsites (other tongue site). Bars represent mean $+$ SD. Dots represent values of individual cases. Non-significant differences were observed by Student $t$ -tests (fraction: t ${}_{{19}}=$ 0.57 $p=$ 0.57; density: t ${}_{19}=$ 0.91 $p=$ 0.3; intensity: t ${}_{19}=$ 1.41 $p=$ 0.17).

Figure 9.

Correlation analysis between SOX2 and KLF4 expression in OSCC and OTSCC. Scatter plots depicting KLF4 and SOX2 positive cells fractions (A) and nuclear labeling intensities (B) for all analyzed cases. Each dot represents an individual case. A. No significant correlation was observed between KLF4 and SOX2 positive cells fraction in OTSCC (T) R ${}^{2}=$ 0.0001, F ${}_{1,19}=$ 0.0018 $p=$ 0.96, nor in the rest of the SCC developed in the oral cavity (NT) R ${}^{2}=$ 0.067, F ${}_{1,12}=$ 1.94 $p=$ 0.19, neither in all regions combined together R ${}^{2}=$ 0.051, F ${}_{1,33}=$ 1.77 $p=$ 0.19. B. No significant correlation was observed between KLF4 and SOX2 nuclear labeling intensity in OTSCC (T) R ${}^{2}=$ 0.0018, F ${}_{1,19}=$ 0.033 $p=$ 0.86, nor in the rest of the SCC developed in the oral cavity (NT) R ${}^{2}=$ 0.074, F ${}_{1,12}=$ 0.96 $p=$ 0.35, neither in all regions combined together R ${}^{2}=$ 0.0087, F ${}_{1,33}=$ 0.29 $p=$ 0.59.

To extend our analysis we analyzed SOX2 expression by immunohistochemistry in the OSCC cases already studied for KLF4 expression. SOX2 was localized in the nucleus of the OSCC cells throughout the parenchyma of the carcinomas at different levels of expression (Fig. 6). The quantitative analysis of the SOX2 immunostaining revealed statistically significant differences between OSCC of different histological grade (Fig. 7). We found increased fraction and density of SOX2 positive cells in moderate and poorly differentiated carcinomas when comparing with well differentiated OSCC. Our data is in concordance with previous ones already reported for OSCC [33]. Therefore we studied if there is a differential expression for SOX2 in OTSCC from different areas of the oral tongue. No statistically significant differences were found regarding SOX2 expression parameters between OTSCC from the border of the tongue versus OTSCC developed in others areas of the tongue (Fig. 8). These results suggest that KLF4 differential expression within tongue subsites is not related to SOX2 differences between these areas. To further evaluate if KLF4 expression is associated with SOX2 labeling we calculated the correlation between these markers using mRNA data from TGCA database and immunolabeling results from our cases. Although KLF4 mRNA expression displayed a significant correlation with SOX2 mRNA levels (Supplementary Fig. 1B) the extremely low R2 obtained indicates that KLF4 expression is not a good predictor of SOX2 mRNA levels for tongue or non-tongue locations. Furthermore, when assessing our immunolabeled samples we found that none of the parameters studied show a significant correlation between KLF4 and SOX2 expression in OSCC (Fig. 9). Limiting the analysis to oral tongue samples produced the same non-significant results supporting the idea that KLF4 and SOX2 can be considered as independent markers in OSCC in general and in tongue SCC in particular.

4. Discussion

We describe here for the first time a differential KLF4 pattern of expression in OTSCC from different tongue subsites. We performed a quantitative analysis on digital images and we showed a reduction in KLF4 positives cells in OTSCC from the lateral borders of the tongue when compared with carcinomas developed in the rest of the oral tongue. However OTSCC did not differ in KLF4 expression versus SCC developed in others oral cavity anatomical locations. Previous reports assessing KLF4 expression levels in human cancer support the idea that KLF4 inhibit carcinogenesis in several tissues and organs as mentioned in the introduction. However, when considering squamous epithelium carcinogenesis the evidence is not compelling. While a qualitative assessment of KLF4 expression in skin SCC showed a decrease in its expression [34] other report supports the notion that KLF4 expression was increased and significantly associated with moderately and poorly differentiated skin SCC [35]. Regarding to OSCC the data is also contradictory. For instance, it has been reported that KLF4 expression was related with tumor grade and its expression was much lower in poorly differentiated oral cancers than in well differentiated cancer cells [22]. Conversely, no significant differences between KLF4 expression and the degree of tumor differentiation have been reported elsewhere [21, 36]. These discrepancies may arise from the methodology used since these studies had assessed KLF4 expression in a non-quantitative manner. In most cases a categorical classification based on visual inspection or semi-quantitative methods was used and the analysis was performed on a rather limited number of sampled areas. In our study, using a quantitative large-scale assessment, we could not find any differential expression of KLF4 across tumors with dissimilar differentiation levels confirming that KLF4 is not associated with tumor differentiation at least in OSCC. The lack of association between KLF4 levels and differentiation degree could be attributed to the diverse OSCC localization analyzed, however, when our analysis was restricted to the tongue, the results were similar to the ones obtained combining tumors from all oral locations, supporting our interpretation that KLF4 is not linked to histological degree in OSCC.

In this regard, inspection of mRNA expression data from a public oral TCGA database (42 samples) that includes case differentiation status, suggests an association between mRNA KLF4 levels and differentiation status of the tumor [37]. However, recent studies have shown that discrepancies between mRNA and protein abundance are common [38, 39]. A comprehensive tissue- and gene-specific analysis revealed that across the 14 tissues studied (including normal and cancer samples) correlation between mRNA and protein expression was positive in only 6.1% of the samples [39]. Furthermore, neither KLF4 nor SOX2 were reported in the group of genes and proteins with highly correlated expression, probably as consequence of differences in mRNA/protein in vivo half-life or varied post-transcriptional mechanisms involved in turning mRNA into protein.

In particular, our results in OTSCC show a significant reduction in the fraction and density of KLF4 positives cells in carcinomas developed in the border of the tongue. In this regard, epidemiological evidence suggests that biological differences and distinct clinical behaviors do exist between tumors arising from different oral cavity subsites and should be examined as different entities [40]. Furthermore carcinomas arising from different subsites within mobile tongue present distinct prognosis [41]. For instance, different series of cases reported the lateral side of the tongue as the more affected area of the mobile tongue as well as a high risk oral site distinctive from the rest of the oral cavity [42, 43].

Animal studies also support the oral tongue subsites differences reported in humans regarding the biological behavior from tumors originated in the lateral border or other tongue subsite. Evidence from genetic manipulation of KLF4 expression when it was restricted to a cellular compartment that includes the oral epithelial stem cells, showed oral premalignant lesions development in the tongue margins [44, 45]. Furthermore, the concomitant Klf4 deletion and Ras activation was sufficient to promote the development of SCC in the margins of the tongue [44]. The results presented here constitute a first attempt to uncover the biological basis aimed to explain clinical behavioral differences reported in tumors arising from diverse locations within the tongue, and it might be revealing an anatomical-dependent function during oral tongue carcinogenesis as proposed for KLF4 in others cancer models [20, 46].

The combined evidence obtained from our previous studies using mouse models of KLF4 ablation in a tongue cell compartment containing epithelial stem cells lead us to explore the relationship in human cancer between KLF4 and SOX2, a stem cell marker involved in carcinogenesis. The transcription factor SOX2 it is extensively known for its roles during mammalian development and its critical roles in the self-renewal and maintenance of embryonic and adult normal stem cells as well as in CSC [24]. In particular, SOX2 expression has been reported increased at protein level for glioblastoma, esophageal, breast and prostate cancer and the SOX2 increased expression correlates with increased histologic grade [24, 47]. Regarding to SOX2 role in oral cancer the reported data is controversial. Bayo et al. showed that tumors with high SOX2 levels of expression have a highest progression-free survival when compared to SOX2 expressing tumors at lower levels [48]. However this seems contradictory with the reported results demonstrating that SOX2 amplifications are common in OSCC and SOX2 amplification is significantly associated with the pathological grade [33, 49]. Our SOX2 data reported here is aligned with the later proposal since we found increased fraction and density of SOX2 positive cells in moderate and poorly differentiated carcinomas when comparing with well differentiated OSCC. This discrepancy could be apparent when considering a recent proposal where histological degree does not correlate with survival outcome [50, 51, 52].

Considering that SOX2 expression was linked to tumor differentiation status and KLF4 expression differ across tongue subsites we did not anticipated a strong association between these two markers of stemness. However, previous studies qualitatively assessing their expression have described some degree of association between them in OSCCs [36]. Thus, to further clarify this point we took advantage of our quantitative data and formally tested the relationship between the expression levels of these two markers by calculating the correlation index in a case by case manner. Our result could not identify any significant correlation, neither in OSCC nor in tongue samples, between KLF4 and SOX2 expression. These results reinforce the notion that they might be involved in different stages of the oral carcinogenesis process.

5. Conclusion

We show here for the first time that OTSCC developed in the borders of the tongue contain fewer tumor cells positive for KLF4 expression than their counterpart from the rest of the tongue. Interestingly this reduction takes place without changes in KLF4 expression level per cell. KLF4 may be useful in understanding the biological behavior of carcinomas originated in different oral tongue subsites. In addition the present work highlights the potential benefits of unbiased quantitative whole slide scanning analysis for diagnostic. Our study prompts future OTSCC prospective studies looking for clinical prognosis to incorporate detailed subsite information in the analysis.

Footnotes

Acknowledgments

We thank German La Iacona from the Microscopy Facility IFIBIO-Houssay and Juan Belforte who kindly provided image analysis technical support. The valuable technical advice for immunostaining of Valeria Sanchez is also acknowledged.

This work was supported by grants from the National Agency for the Promotion of Science and Technology, Argentina (ANPCyT, PICT 2017-0361) and the Argentinean National Research Council (CONICET, PIP11220150100251CO).

Conflict of interest

The authors state that they have no conflicts of interest.

Author’s contributions

Conception: Ana Rosa Raimondi.

Interpretation or analysis of data: Maria Luisa Paparella, Dario Martin Ferri, Krissya Maria Villegas.

Preparation of the manuscript: Maria Luisa Paparella, Ana Rosa Raimondi.

Revision for important intelectual content: Ana Rosa Raimondi.

Supervision: Ana Rosa Raimondi.

Supplementary data

Supplementary Fig. 1.

A. Left. KLF4 mRNA expression levels in OTSCC (T) and the rest of the SCC developed in the oral cavity (NT). Bars represent mean $+$ SD of 130 (T) and 186 (NT) cases analyzed from the TCGA database. Non-significant differences were observed by Student $t$ -tests (t ${}_{314}=$ 0.82, $p=$ 0.42). Right. Individual values of KLF4 mRNA expression for each case according to the oral cavity subsite (the subsites are color coded as shown in graph legend). No significant differences were observed between groups for CV ( $p=$ 0.85), after Shapiro-Wilk test. Each dot represents one individual biopsy from the TCGA data set. B. Scatter plots depicting KLF4 and SOX2 mRNA expression levels. Each dot represents an individual case from tongue (left) and the rest of the SCCs developed in the oral cavity (right). Results from Pearson correlation analyses are presented inside the graphs (blue texts). The analyzed data set included 130 oral tongue SCC and 186 oral cavity SCCs excluding the tongue (18 from the alveolar ridge, 22 from the buccal mucosa, 63 from the floor of mouth, 7 from the hard palate, 3 from the lip, and 73 form the oral cavity).

References

Ferlay

Soerjomataram

Dikshit

Eser

Mathers

Rebelo

Parkin

D.M.

Forman

and Bray

, Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012, International Journal of Cancer 136 (2015), E359–E386.

Siegel

R.L.

Miller

K.D.

and Jemal

, Cancer statistics, 2018, CA: A Cancer Journal for Clinicians 68 (2018), 7–30.

Sathyan

K.M.

Sailasree

Jayasurya

Lakshminarayanan

Abraham

Nalinakumari

K.R.

Abraham

E.K.

and Kannan

, Carcinoma of tongue and the buccal mucosa represent different biological subentities of the oral carcinoma, Journal of Cancer Research and Clinical Oncology 132 (2006), 601–609.

Rautava

Luukkaa

Heikinheimo

Alin

Grenman

and Happonen

R.P.

, Squamous cell carcinomas arising from different types of oral epithelia differ in their tumor and patient characteristics and survival, Oral Oncology 43 (2007), 911–919.

Hussein

A.A.

Forouzanfar

Bloemena

de Visscher

Brakenhoff

R.H.

Leemans

C.R.

and Helder

M.N.

, A review of the most promising biomarkers for early diagnosis and prognosis prediction of tongue squamous cell carcinoma, British Journal of Cancer 119 (2018), 724–736.

Warnakulasuriya

, Global epidemiology of oral and oropharyngeal cancer, Oral Oncology 45 (2009), 309–316.

Slootweg

and Eveson

, Tumours of the Oral Cavity and Oropharynx, in: Barnes

Eveson

Reichart

and Sidransky

(Eds.), Pathology and Genetics of Head and Neck Tumours, World Health Organization Classification of Tumours, 3er ed., IARC Press, 2005: pp. 166–167.

Tota

J.E.

Anderson

W.F.

Coffey

Califano

Cozen

Ferris

R.L.

, John

St.

Cohen

E.E.W.

and Chaturvedi

A.K.

, Rising incidence of oral tongue cancer among white men and women in the United States, 1973–2012, Oral Oncology 67 (2017), 146–152.

Hussein

A.A.

Helder

M.N.

de Visscher

J.G.

Leemans

C.R.

Braakhuis

B.J.

de Vet

H.C.W.

and Forouzanfar

, Global incidence of oral and oropharynx cancer in patients younger than 45 years versus older patients: a systematic review, European Journal of Cancer (Oxford, England: 1990) 82 (2017), 115–127.

10.

Bello

I.O.

Soini

and Salo

, Prognostic evaluation of oral tongue cancer: means, markers and perspectives (II), Oral Oncology 46 (2010), 636–643.

11.

Rusthoven

Ballonoff

Raben

and Chen

, Poor prognosis in patients with stage I and II oral tongue squamous cell carcinoma, Cancer 112 (2008), 345–351.

12.

Demeng

and Cun-Yu

, Targeting cancer stem cells in squamous cell carcinoma, Precision Clinical Medicine 2 (2019), 152–165.

13.

Nör

J.E.

and Gutkind

J.S.

, Head and neck cancer in the new era of precision medicine, Journal of Dental Research 97 (2018), 601–602.

14.

Ferris

R.L.

Blumenschein

Fayette

Guigay

Colevas

A.D.

Licitra

Harrington

Kasper

Vokes

E.E.

Even

Worden

Saba

N.F.

Docampo

L.C.I.

Haddad

Rordorf

Kiyota

Tahara

Monga

Lynch

Geese

W.J.

et al., Nivolumab for recurrent squamous-cell carcinoma of the head and neck, New England Journal of Medicine 375 (2016), 1856–1867.

15.

Dagogo-Jack

and Shaw

A.T.

, Tumour heterogeneity and resistance to cancer therapies, Nature Reviews Clinical Oncology 15 (2018), 81–94.

16.

S.S.

and Cirillo

, The molecular markers of cancer stem cells in head and neck tumors, Journal of Cellular Physiology 235 (2020), 65–73.

17.

Takahashi

Tanabe

Ohnuki

Narita

Ichisaka

Tomoda

and Yamanaka

, Induction of pluripotent stem cells from adult human fibroblasts by defined factors, Cell 131 (2007), 861–872.

18.

Rowland

B.D.

and Peeper

D.S.

, KLF4, p21 and context-dependent opposing forces in cancer, Nat.Rev.Cancer 6 (2006), 11–23.

19.

Ghaleb

A.M.

and Yang

V.W.

, Krüppel-like factor 4 (KLF4): what we currently know, Gene 611 (2017), 27–37.

20.

Tetreault

M.P.

Yang

and Katz

J.P.

, Krüppel-like factors in cancer, Nature Reviews Cancer 13 (2013), 701–713.

21.

Tai

S.K.

Yang

M.H.

Chang

S.Y.

Chang

Y.C.

W.Y.

Tsai

T.L.

Wang

Y.F.

Chu

P.Y.

and Hsieh

S.L.

, Persistent Krüppel-like factor 4 expression predicts progression and poor prognosis of head and neck squamous cell carcinoma, Cancer Science 102 (2011), 895–902.

22.

Liu

Zhou

Liu

and Zhang

, The janus-faced roles of Krüppel-like factor 4 in oral squamous cell carcinoma cells, Oncotarget 6 (2015), 44480–44494.

23.

Sarkar

and Hochedlinger

, The sox family of transcription factors: versatile regulators of stem and progenitor cell fate, Cell Stem Cell 12 (2013), 15–30.

24.

Novak

Hüser

Elton

J.J.

Umansky

Altevogt

and Utikal

, SOX2 in development and cancer biology, Seminars in Cancer Biology 67 (2020), 74–82.

25.

Okubo

Clark

and Hogan

B.L.M.

, Cell lineage mapping of taste bud cells and keratinocytes in the mouse tongue and soft palate, Stem Cells 27 (2009), 442–450.

26.

Takata

and Slootweg

, Tumors of the oral cavity and mobile tongue, in: L-Naggar

Chan

Grandis

Takata

(Eds.), World Health Organization Classification of Head and Neck Tumours, 4th ed., IARC: Lyon, 2017: pp. 105–131.

27.

Raimondi

A.R.

Vitale-Cross

Amornphimoltham

Gutkind

J.S.

and Molinolo

, Rapid development of salivary gland carcinomas upon conditional expression of K-ras driven by the cytokeratin 5 promoter, Am.J.Pathol. 168 (2006), 1654–1665.

28.

Bankhead

Loughrey

M.B.

Fernández

J.A.

Dombrowski

McArt

D.G.

Dunne

P.D.

McQuaid

Gray

R.T.

Murray

L.J.

Coleman

H.G.

James

J.A.

Salto-Tellez

and Hamilton

P.W.

, QuPath: open source software for digital pathology image analysis, Scientific Reports 7 (2017), 16878.

29.

Gelebart

Hegazy

S.A.

Wang

Bone

K.M.

Anand

Sharon

Hitt

Pearson

J.D.

Ingham

R.J.

and Lai

, Aberrant expression and biological significance of Sox2, an embryonic stem cell transcriptional factor, in ALK-positive anaplastic large cell lymphoma, Blood Cancer Journal 2 (2012), e82.

30.

Hart

G.T.

Hogquist

K.A.

and Jameson

S.C.

, Krüppel-like factors in lymphocyte biology, The Journal of Immunology 188 (2012), 521–526.

31.

Segre

J.A.

, Epidermal barrier formation and recovery in skin disorders, J.Clin.Invest 116 (2006), 1150–1158.

32.

Shibata Masaki

I.K.

, Chiba Tadashige, Matsuoka Takanori, Mihara Nozomi, Kawashiri Shuichi, Krüppel-like factors 4 and 5 expression and their involvement in differentiation of oral carcinomas, Int J Clin Exp Pathol 8 (2015), 3701–3709.

33.

and ZJ

, CD163+ tumor-associated macrophages correlated with poor prognosis and cancer stem cells in oral squamous cell carcinoma, BioMed Research International 2014 (2014), 838632.

34.

Choi

W.J.

Youn

S.H.

Back

J.H.

Park

E.J.

Kim

K.J.

Park

H.R.

Kim

A.L.

and Kim

K.H.

, The role of KLF4 in UVB-induced murine skin tumor development and its correlation with cyclin D1, p53, and p21(Waf1/Cip1) in epithelial tumors of the human skin, Arch.Dermatol.Res. 303 (2011), 191–200.

35.

Chen

Y.J.

C.Y.

Chang

C.C.

C.J.

M.C.

and Chen

C.M.

, Nuclear Krüppel-like factor 4 expression is associated with human skin squamous cell carcinoma progression and metastasis, Cancer Biology and Therapy 7 (2008), 777–782.

36.

Yoshihama

Yamaguchi

Imajyo

Mine

Hiyake

Akimoto

Kobayashi

Chigita

Kumamaru

Kiyoshima

Mori

and Sugiura

, Expression levels of SOX2, KLF4 and brachyury transcription factors are associated with metastasis and poor prognosis in oral squamous cell carcinoma, Oncology Letters 11 (2016), 1435–1446.

37.

Pickering

C.R.

Zhang

Yoo

S.Y.

Bengtsson

Moorthy

Neskey

D.M.

Zhao

Ortega Alves

M.V.

Chang

Drummond

Cortez

Xie

T.X.

Zhang

Chung

Issa

J.P.J.

Zweidler-McKay

P.A.

El-Naggar

A.K.

Weinstein

J.N.

Wang

et al., Integrative genomic characterization of oral squamous cell carcinoma identifies frequent somatic drivers, Cancer Discovery 3 (2013), 770–781.

38.

Greenbaum

Colangelo

Williams

and Gerstein

, Comparing protein abundance and mRNA expression levels on a genomic scale, Genome Biology 4 (2003), 117.

39.

Kosti

Jain

Aran

Butte

A.J.

and Sirota

, Cross-tissue analysis of gene and protein expression in normal and cancer tissues, Scientific Reports 6 (2016), 24799.

40.

Montero

P.H.

Patel

P.D.

Palmer

F.L.

Patel

S.G.

Shah

J.P.

Hayes

R.B.

and Ganly

, Changing trends in smoking and alcohol consumption in patients with oral cancer treated at Memorial Sloan-Kettering Cancer Center from 1985 to 2009, Archives of Otolaryngology – Head and Neck Surgery 138 (2012), 817–822.

41.

Sano

and Myers

J.N.

, Metastasis of squamous cell carcinoma of the oral tongue, Cancer and Metastasis Reviews 26 (2007), 645–662.

42.

Zini

Czerninski

and Sgan-Cohen

H.D.

, Oral cancer over four decades: epidemiology, trends, histology, and survival by anatomical sites, Journal of Oral Pathology and Medicine 39 (2010), 299–305.

43.

Ion Ciucă Mărăşescu

F.I.

Marasescu

P.C.

Matei

Florescu

A.M.

Margaritescu

Petrescu

S.M.S.

and Dumitrescu

C.I.

, Epidemiological and histopathological aspects of tongue squamous cell carcinomas-retrospective study, Current Health Sciences Journal 44 (2018), 211–224.

44.

Abrigo

Alvarez

Paparella

M.L.

Calb

D.E.

Bal de Kier Joffe

Gutkind

J.S.

and Raimondi

A.R.

, Impairing squamous differentiation by Klf4 deletion is sufficient to initiate tongue carcinoma development upon K-Ras activation in mice, Carcinogenesis 35 (2014), 662–669.

45.

Paparella

M.L.

Abrigo

Bal de Kier Joffe

and Raimondi

A.R.

, Oral-specific ablation of Klf4 disrupts epithelial terminal differentiation and increases premalignant lesions and carcinomas upon chemical carcinogenesis, Journal of Oral Pathology & Medicine 44 (2015), 801–809.

46.

Rowland

B.D.

Bernards

and Peeper

D.S.

, The KLF4 tumour suppressor is a transcriptional repressor of p53 that acts as a context-dependent oncogene, Nat.Cell Biol. 7 (2005), 1074–1082.

47.

Wang

Z.S.Q.

Wang

Giercksky

K.E.

and Nesland

J.M.

, Oct3/4 and Sox2 are significantly associated with an unfavorable clinical outcome in human esophageal squamous cell carcinoma – PubMed, Anticancer Res (2009), 1233–1241.

48.

Bayo

Jou

Stenzinger

Shao

Gross

Jensen

Grabe

Mende

C.H.

Rados

P.V.

Debus

Weichert

Plinkert

P.K.

Lichter

Freier

and Hess

, Loss of SOX2 expression induces cell motility via vimentin up-regulation and is an unfavorable risk factor for survival of head and neck squamous cell carcinoma, Molecular Oncology 9 (2015), 1704–1719.

49.

Kokalj Vokač

Čizmarević

Zagorac

Zagradišnik

and Lanišnik

, An evaluation of SOX2 and hTERC gene amplifications as screening markers in oral and oropharyngeal squamous cell carcinomas, Molecular Cytogenetics 7 (2014), 5.

50.

O-charoenrat

Pillai

Patel

Fisher

Archer

Eccles

and Rhys-Evans

, Tumour thickness predicts cervical nodal metastases and survival in early oral tongue cancer, Oral Oncology 39 (2003), 386–390.

51.

Okamoto

Nishimine

Kishi

Kirita

Sugimura

Nakamura

and Konishi

, Prediction of delayed neck metastasis in patients with stage I/II squamous cell carcinoma of the tongue, Journal of Oral Pathology and Medicine 31 (2002), 227–233.

52.

Al-Rajhi

Khafaga

El-Husseiny

Saleem

Mourad

Al-Otieschan

and Al-Amro

, Early stage carcinoma of oral tongue: prognostic factors for local control and survival, Oral Oncology 36 (2000), 508–514.

Quantitative analysis of KLF4 and SOX2 expression in oral carcinomas reveals independent association with oral tongue subsite location and histological grade

Abstract

BACKGROUND:

METHODS AND RESULTS:

CONCLUSIONS:

Keywords

1. Introduction

2. Materials and methods

Table 1 Clinicopathologic characteristic for 35 selected cases of OSCC

2.2 Qualitative staining evaluation

2.3 Image acquisition and analysis

2.5 The Cancer Genome Atlas (TCGA) data

2.6 Statistical analysis

5. Conclusion

Footnotes

Acknowledgments

Conflict of interest

Author’s contributions

Supplementary data

References

Table 1
Clinicopathologic characteristic for 35 selected cases of OSCC