- Open Access
Aberrantly hypermethylated tumor suppressor genes were identified in oral squamous cell carcinoma (OSCC)
Clinical Epigeneticsvolume 11, Article number: 116 (2019)
Oral squamous cell carcinoma (OSCC) is a genetic and epigenetic disease. There is growing evidence to suggest that environmental factors due to epigenetic changes can be involved in the OSCC pathogenesis. Although tumor suppressor genes (TSGs) are commonly inactivated by promoter hypermethylation in human cancers, the epigenetic changes and the mechanism of TSGs in human OSCC remain unclear. We therefore assessed the methylation status of the TSGs, which are associated with epigenetic silencing in human cancers, OSCC cell lines, primary tumors, and normal oral mucosa.
We used 14 TSGs that were originally identified in colon cancer to investigate the aberrant hypermethylation of these genes associated with transcriptional silencing in 10 OSCC cell lines. We found three TSGs, TFPI2, SOX17, and GATA4, that are robustly hypermethylated and are associated with transcriptional silencing in OSCC cell lines. The re-expression of the three genes was induced by 5-aza-2′-deoxycytidine (5-aza-dC) in cells in which these genes were not expressed or had a lack of expression. In 33 cases of primary OSCC tumors, promoter hypermethylation was detected for the TFPI2, SOX17, and GATA4 genes at (32/33) 97%, (22/33) 67%, and (11/33) 33%, respectively. Eleven normal oral mucosa samples showed no promoter hypermethylation for all three genes, which suggests that this promoter hypermethylation is cancer-specific. Bisulfite sequencing analysis confirmed the cancer-specific methylation of the TFPI2, SOX17, and GATA4 promoters in the OSCC cell lines and tumors but not in the normal oral mucosa samples. More importantly, the methylation status of TFPI2, GATA4, and SOX17 was significantly associated with OSCC patients’ overall survival through TCGA DNA methylation database.
We identified that TFPI2, SOX17, and GATA4 are frequently hypermethylated in human OSCC cells in a cancer-specific manner and that the transcriptional expression of these genes is regulated by promoter hypermethylation in OSCC. Our results highlight the great potential used as a synergistic biomarker set to improve the prognosis and therapeutic treatment for patients with OSCC.
Epigenetic modification, including DNA methylation and many types of histone modifications, is responsible for the altered gene expression patterns that allow for specific phenotypes . DNA methylation is the primary and most studied epigenetic modification  and plays an important role in normal mammalian development, but aberrant methylation patterns are correlated with several differentiation-related diseases, including many types of human cancers. Gene promoter hypermethylation in cancer can silence gene expression and regulate biological processes, especially tumor suppressor genes (TSGs), which have an important role in cancer initiation, progression, and metastasis . Aberrant DNA methylation is a common phenomenon in malignancies, and the methylation profiles are altered in various tumors, which might be associated with clinical outcomes .
We and others have reported that the transcriptional silencing of tumor suppressor candidate genes is regulated by promoter hypermethylation in colorectal cancer [3, 4]. High-throughput genome-wide methylation studies offer a sophisticated strategy to understand the significance of DNA methylation and its impact on gene regulation [3, 5].
Numerous studies have used these well-known tumor suppressor candidate genes in other cancer types, such as lung, breast, and pancreatic cancer types [4, 6]. However, there is a lack of studies on the DNA methylation profiles of these TSGs in oral squamous carcinoma cancer (OSCC). Very recent findings suggest that the oral microbiota could contribute to CRC development, implying that the profile of epigenetic changes may be shared between two different cancer types [7, 8].
Oral squamous cell carcinoma (OSCC), the most prevalent type of head and neck SCC (HNSCC), typically behaves in an aggressive manner, frequently leading to local invasion and early lymph node metastasis . The overall 5-year survival rate for OSCC is approximately 60%  and shows only modest improvement over the past two decades, despite considerable improvements in the treatment of OSCC [11, 12]. In the last decade, there has been a lack of understanding of OSCC carcinogenesis, progression, and metastasis in terms of epigenetic regulation and its clinical application. Therefore, the identification of molecular changes for the oncogenes or TSGs associated with OSCC will help improve diagnosis predictions and early treatment [13, 14].
In this study, we investigated transcriptional silencing by the promoter hypermethylation of 18 tumor suppressor candidate genes (MGMT, sFRP2, HIC1, sFRP4, Timp3, sFRP5, TFPI2, p16, sFRP1, E-cad, SOX17, GATA4, GATA5, p14, FBN2, p15, and TCERG1L) in 10 OSCC cell lines and a small set of OSCC primary tumor samples (n = 33). Remarkably, we found that most of the genes are silenced by promoter hypermethylation in OSCC cell lines. Among these genes, TFPI2, SOX17, and GATA4 were frequently hypermethylated in both OSCC and oral cancer patient samples in a cancer-specific manner, suggesting that these genes may play a role as tumor suppressors in OSCC.
Identification of aberrantly hypermethylated tumor suppressor candidates in oral squamous cell carcinoma (OSCC)
In this study, we investigated transcriptional silencing by the promoter hypermethylation of 14 well-known tumor suppressor candidate genes (MGMT, sFRP2, sFRP1, HIC1, sFRP4, Timp3, sFRP5, TFPI2, p16, SOX17, GATA4, GATA5, FBN2, and TCERG1L) in OSCC cell lines. These genes were originally identified from comprehensive studies of colon cancer as tumor suppressor candidate genes, and their transcriptional silencing is regulated by promoter hypermethylation in many other cancer types . First, we examined whether the expression of these genes is regulated by promoter methylation in 10 OSCC cell lines (Ca9-22, HSC4, OSC20, SAS, SCC25, HN22, YD10B, YD38, YD9, and HSC3) by reverse transcriptase polymerase chain reaction (RT-PCR). We treated the OSCC cell lines with the demethylation agent 5-aza-2′-deoxycytidine (5-aza-dC) to determine whether these TSGs are re-expressed after 5-aza-dC treatment. To correlate the transcriptional expression of TSGs with promoter hypermethylation, we assessed the methylation analysis of these TSGs in OSCC cell lines by methylation-specific PCR (MSP). We performed massive RT-PCR and MSP to identify whether these TSGs are silenced by promoter hypermethylation in OSCC. The results for transcriptional expression and promoter hypermethylation by the RT-PCR and MSP analyses of TSGs are summarized in Fig. 1a (Additional file 1: Figure S1) with the following three categories: (1) genes are silenced by promoter hypermethylation with re-expression by 5-aza-dC treatment, (2) genes are silenced by promoter hypermethylation without re-expression by 5-aza-dC treatment, and (3) gene expression is not correlated with promoter methylation. The methylation frequencies of 14 TSGs in 10 OSCC cell lines ranged from 50 to 100% (Fig. 1b), which suggests that all of the genes we tested are frequently hypermethylated in OSCC cell lines.
The silencing of most TSGs is regulated by promoter hypermethylation, and the genes are re-expressed by 5-aza-dC treatment. We found that the transcriptional silencing of 9 genes (MGMT, HIC1, sFRP4, Timp3, TFPI2, SOX17, GATA4, FBN2, and TCERG1L) out of 14 genes by promoter hypermethylation was correlated with gene re-expression by 5-aza-dC treatment in 2 or more OSCC cell lines. Nonetheless, we decided to focus on 4 TSGs (TFPI2, SOX17, GATA4, and FBN2) because we later confirmed that these genes had a cancer-specific methylation pattern in primary OSCC tumor samples (Fig. 3). We determined whether these 4 genes are re-expressed after 5-aza-dC treatment in 10 OSCC cell lines by real-time RT-PCR (Fig. 2). We postulate that these genes were re-expressed in OSCC cells after 5-aza-dC treatment, and we should detect their promoter hypermethylation in those cells. To verify that transcriptional silencing was due to promoter hypermethylation, we assessed the methylation levels of the TFPI2, SOX17, GATA4, and FBN2 genes by MSP. Overall, we identified four TSGs (TFPI2, SOX17, GATA4, and FBN2) that are regulated by promoter hypermethylation associated with transcriptional silencing in OSCC cell lines.
TFPI2, SOX17, and GATA4 are frequently hypermethylated in OSCC cells and primary OSCC tumors in a cancer-specific manner
To validate that four TSGs were indeed regulated by promoter hypermethylation, we extended our methylation analysis for these genes in a panel of primary OSCC tumor samples (n = 33) and normal oral mucosa samples (n = 11) obtained from cancer-free individuals. The clinical characteristics of the patients are documented in Table 1. Four genes (TFPI2, SOX17, GATA4, and FBN2) satisfied the criteria of “cancer-specific methylation” with high-frequency methylation in cell lines, no/undetectable methylation in normal oral mucosa (Fig. 3a), and frequent methylation in primary OSCC tumor samples (Fig. 3b). Based on our data, we decided to eliminate FBN2 for further experimental analyses because of the low frequency (< 20%) of methylation in primary OSCC tissue samples.
We next assessed changes in DNA methylation status in the promoter region of the three genes by bisulfite genomic sequencing in representative OSCC cell lines (Ca9-22, HSC3, OSC-20, SAS, and YD10B), OSCC primary tumor samples (n = 5), and normal oral mucosa samples (n = 3) (Fig. 4). The location of promoter regions (TFPI2: upstream region from − 286 to − 76, 30 CpG sites; SOX17: region from + 339 to + 603, 34 CpG sites; and GATA4: upstream region from − 229 to 161, 46 CpG sites) for sequencing is relative to the transcription start sites (TSS) of exon 1. The bisulfite sequencing of five individual clones of PCR products from five OSCC cell lines revealed densely methylated CpG sites within the promoter regions of all of the clones. The methylation status of the TFPI2 promoter region in OSCC cells and primary tissues showed that most of the CpG sites (> 99%) we sequenced were methylated compared to those in the normal oral mucosal tissues (< 11%). Both the SOX17 and GATA4 genes in OSCC cells and primary tissues (59–77% for SOX17 and 71–91% for GATA4) showed relatively dense methylation than those in normal tissues (< 32% for SOX17 and< 23% for GATA4). Accordingly, these bisulfite sequencing data confirm the complete methylation of the three gene promoters in OSCC cell lines, primary tumor tissues, and normal oral mucosal samples determined by MSP (Fig. 2), but they also explain the lack of these three gene expression results in OSCC cell lines obtained by quantitative RT-PCR analysis. Notably, all three genes were significantly hypermethylated in primary OSCC tumors compared to those in the normal tissues. Overall, these results strongly suggest that the TFPI2, SOX17, and GATA4 genes showed increasing DNA methylation levels in primary OSCC tumors compared with the DNA methylation levels in normal oral tissues.
TFPI2, SOX17, and GATA4 protein expression is downregulated in primary OSCC tumor tissues
Next, we examined the protein expression of the three genes by immunohistochemical staining. Normal oral mucosal tissue showed strong nuclear positivity against TFPI2, SOX17, and GATA4 as well as cytoplasmic staining for TFPI2 and SOX17 (Fig. 5, Additional file 1: Figure S2). In contrast, those proteins were negatively expressed in both the nucleus and cytoplasm of the primary OSCC tissue.
Methylation of TFPI2, SOX17, and GATA4 is associated with OSCC patient survival
To investigate whether the methylation status of the three potential biomarkers (TFPI2, SOX17, and GATA4) described above had an impact on overall survival of OSCC patients, TCGA DNA methylation data (CpG sites) with clinical information (n = 344) were divided into two groups (hypermethylated group, n = 172; hypomethylated group, n = 172) according to their methylation status on TFPI2, SOX17, or GATA4. Kaplan–Meier analysis was performed on either each CpG on the promoter region or combinations of CpGs. To explore the potential clinical significance of these data, we tested these above three tumor suppressor genes for their prognostic potential in a large TCGA cohort (head and neck squamous cell carcinoma; HNSCC) . We subtracted the patients with OSCC from HNSCC TCGA cohort and well-annotated clinical data that could be correlated with survival and gene methylation status. Strikingly, we observed a statistically significant increased risk for mortality when either individual genes or combination of genes were methylated. For example, Kaplan–Meier survival curves for GATA4 and SOX17, individually, show that DNA methylation of each is associated with decreased survival but no significant statistical results for TFPI2 alone (Fig. 6a, p = 0.016 for GATA4, p = 0.0066 for SOX17), even though TPFI2 gene was frequently hypermethylated in both OSCC cells and OSCC primary tumors. We wondered whether the combination of two- or three-gene methylation could have an impact on overall survival in this cohort. Very interestingly, combination of two-gene (SOX17 and GATA4) methylation is associated with significantly worse survival (Fig. 6b, p = 0.013). In addition, SOX17 combination of TFPI2 gene methylation is associated with worse survival as well (Fig, 6)b, p = 0.019, but there are no statistically significant results for combination of GATA4 and TFPI2 genes. More importantly, Kaplan–Meier survival curves for TFPI2, GATA4, and SOX17 show that DNA methylation of three-gene combination is associated with significantly worse survival (Fig. 6c, p = 0.0041). These findings support our hypothesis that DNA methylation of TSGs may have an important role during OSCC cancer progression as well as clinical application as prognostic biomarkers.
The aberrant promoter hypermethylation in the adjacent of transcription start sites (TSSs) often leads to alterations in gene function and abnormal cellular pathway in human cancer. Epigenetic events linked to TSG inactivation through promoter methylation are as frequent as somatic mutations in cancer and contribute to drive tumor initiation and progression. The promoter methylation of 14 selected TSGs in this study has previously been identified in colon cancer [3, 16] and has been characterized in other cancer types . Recently, it has been demonstrated that OSCC involves genetic alterations and epigenetic mechanisms play crucial roles . Epigenetically regulated genes have become an important tool for better understanding cancer initiation, progression, and development.
Interestingly, recent studies have demonstrated that Fusobacterium nucleatum (Fn), a specific species of the gut microbiota, may be closely associated with colorectal carcinogenesis . A growing body of evidence suggests that Fn promotes tumor progression and inhibits the antitumor immune response in the colorectum by regulating the β-catenin pathway [8, 18]. In addition, Fn DNA is significantly enriched in the microsatellite instability-high (MSI-H) molecular subtypes of CRC samples [19, 20]. These studies encouraged us to investigate whether the epigenome profile of colon cancer might be shared with OSCC. We therefore tested the hypothesis that selected TSGs, which are frequently hypermethylated in colon cancer, predominantly promote methylation in OSCC tumors. We observed interesting relationships among the most methylated genes (TFPI2, SOX17, and GATA4) in OSCC cell lines and primary tumors in a cancer-specific manner.
Tissue factor pathway inhibitor-2 (TFPI2), a Kunitz-type serine proteinase inhibitor, can inhibit a variety of serine proteases, including factor VIIa/tissue factor, factor Xa, plasmin, trypsin, chymotrypsin, and plasma kallikrein. Over the last decade, TFPI2 has been identified as a TSG in several types of cancer, including colorectal cancer (CRC) [21, 22]. TFPI2 methylation frequently existed in CRC patients’ sera  and stool samples . Moreover, hypermethylated TFPI2 was associated with recurrence and early-stage CRC , and TFPI2 was significant in CRC patients’ sera with large, poorly differentiated carcinoma, deep invasion, lymph node metastasis, or distant metastasis . Although the methylation of TFPI2 in OSCC has been identified using a genome-wide methylation array, its promoter hypermethylation associated with transcriptional silencing was not validated in OSCC tumor samples . In this study, we validated that TFPI2 is frequently methylated in a panel of OSCC primary tumor samples in a cancer-specific manner and that TFPI2 is transcriptionally regulated by promoter hypermethylation in OSCC cells. Our study thus highlights the finding that TFPI2 is hypermethylated in most OSCC tumor samples we tested in this study, suggesting that TFPI2 may be a useful biomarker for screening OSCC patients. Therefore, further studies in a larger series of samples will be necessary to confirm that the TFPI2 gene can be a useful methylation biomarker for screening OSCC patients.
SOX17 encodes an HMG box transcription factor and has been implicated in oligodendrocyte development, vascular development, the formation of definitive endoderm, and embryonic hematopoiesis [26,27,28,29]. The SOX17 promoter is hypermethylated in cholangiocarcinoma (CCA) tissues, lung cancer , gastric cancer , liver cancer , and breast cancer . Additionally, epigenetically regulated SOX17 may contribute to the abnormal activation of the Wnt signaling pathway in human cancer . The high frequency of SOX17 methylation in 50% of non-small-cell lung cancers and nearly 90% of esophageal squamous cancers strongly supports this possibility and is consistent with the known role of the aberrant activation of Wnt signaling during tumorigenesis for multiple cancer types .
The GATA4 gene encodes a member of the GATA family of zinc-finger transcription factors, which recognize the GATA motif present in the promoters of many genes. Transcription factors GATA4 play an essential role in the development and differentiation of the gastrointestinal tract and are suggested to be involved in colorectal cancer development [36, 37]. GATA4 is known to be more likely to behave as a tumor suppressor gene since the increased expression levels correlate with the terminal differentiation in the intestinal epithelium and the terminal differentiation induced in CRC cells by sodium butyrate .
The loss of GATA4 expression due to promoter hypermethylation has been reported in primary colorectal, gastric, esophageal, lung, and ovarian cancers [37,38,39,40]. To our knowledge, this is the first study to describe how the promoter hypermethylation of both the SOX17 and GATA4 genes is regulated in OSCC cells and how these promoters are frequently hypermethylated in OSCC primary tumors.
Epigenetic alterations of specific genes have recently emerged as potential candidate biomarkers for the early detection of cancer . Recent advances attest to the great promise of DNA methylation markers as powerful tools for future use in the clinic. It would be of great value to find useful biomarkers and prognostic molecular signatures in OSCC to develop novel therapeutic strategies or chemotherapeutic agents. Numerous potential clinical applications of epigenetics for diagnostic and therapeutic applications have been reported .
The promoter hypermethylation of specific genes has already been studied in oral premalignancy [42, 43], but little is known about the description of early DNA methylation changes during oral carcinogenesis using genome-wide profiles. Recently, several studies have reported that early DNA promoter methylation alterations are associated with histological changes during oral tumorigenesis . They confirmed that these alterations are early events during oral tumorigenesis and offer an opportunity for biomarker development .
According to our KM analysis for patients with OSCC from TCGA data, we found a statistically significant increased risk for mortality when either individual genes or combination of genes were methylated. Unlikely SOX17 and GATA4, we did not see any significant clinical correlation with overall survival of TFPI2 alone. Promoter methylation of TFPI2 is detected in most of tumors (over 90% tested cancer patients samples) with OSCC and other cancer types  as well which explained TFPI2 methylation can be useful for early detection biomarker rather than prognostic biomarkers. Previous studies have reported survival-related OSCC methylation biomarkers [46, 47]. These reports described novel methylation biomarkers associated with prognostic biomarkers in OSCC, showing that the methylation frequency of candidate genes ranged from 10 to 65%. However, TFPI2 and SOX17 in our study showed a significant hypermethylation pattern (95% and 67%, respectively) in OSCC primary tumors compared to those observed in previous studies because we used a well-established validation approach using a specific promoter region to identify methylation biomarker candidates. Based on our data, the dysregulation of the TFPI2, SOX17, and GATA4 genes regulated by epigenetic changes may explain a main mechanism linked to OSCC development, and these three-gene panels could serve as synergistic prognostic biomarkers in OSCC treatment. Further studies are needed to verify TFPI2 and SOX17 as methylation biomarkers in a large cohort of OSCC primary tumors along with normal controls, and these genes should be evaluated by testing in blood or sputum DNA.
In summary, we used 14 putative TSGs originally identified in colon cancer to screen the epigenetic regulation of these genes in 10 OSCC cell lines, a small set of primary OSCC tumors, and normal oral mucosa. We identified the transcriptional expression of 9 genes regulated by promoter hypermethylation in OSCC and finally provided 3 genes (TFPI2, SOX17, and GATA4) that were frequently hypermethylated in primary OSCC tumors in a cancer-specific manner. The methylation level in the promoter region of these genes was confirmed by bisulfite genomic sequencing analysis, which suggested that the methylation level is significantly increased in OSCC cell lines and primary tumors compared to the methylation level in normal tissues. Our study thus highlights the finding that the dysregulation of TFPI2, SOX17, and GATA4 genes regulated by epigenetic changes may explain a main mechanism linked to OSCC development, and this three-gene panel has a potential to be used as a synergistic biomarker set capable of improving the prognosis and treatment for patients with OSCC.
For OSCC samples, paraffin-embedded tissue specimens from patients who were diagnosed and surgically treated at the Department of Oral and Maxillofacial Surgery, Pusan National University Dental Hospital, between 2009 and 2014 were retrieved. For normal oral mucosal tissue for comparison, fresh healthy tissue that had been obtained during dental procedures, such as gingivoplasty, third molar extraction, or implant surgery, was used. Clinicopathologic features, including age, gender, histopathologic grade, and the TNM stage, are summarized in Table 1. Tumor grading was based on the ethical permission for this study and was granted by the Institutional Review Board (IRB) of Pusan National University Dental Hospital (IRB No. PNUDH-2019-006).
Cell culture and 5-aza-dC treatment
Ten human OSCC cell lines (Ca9-22, HSC4, OSC20, SAS, SCC25, HN22, YD10B, YD38, YD9, and HSC3) were used in this study. The Ca9-22 and HSC4 cell lines were cultured in MEM/EBSS (HyClone, Logan, UT) medium. The OSC20, SAS, and SCC25 cell lines were cultured in a 1:1 mixture of Dulbecco’s modified Eagle’s medium and Ham’s F-12 Nutrient Mixture (DMEM/F12; HyClone, Logan, UT). The HN22, YD10B, YD38, YD9, and HSC3 cell lines were cultured in a 1:1 mixture of Dulbecco’s modified Eagle’s medium (HyClone, Logan, UT). All cell culture medium was supplemented with 10% fetal bovine serum (HyClone, Logan, UT, USA) and 1% antibiotic–antimycotic (Gibco, Grand Island, NY, USA). All cell lines were incubated at 37 °C in atmospheric conditions of 20% O2 and 5% CO2. To investigate the effects of 5-aza-dC treatment, the cells were treated with 5 μM 5-aza-dC (Sigma, St. Louis, MO, USA) for 72 h.
DNA methylation analysis
For methylation analyses, genomic DNA was isolated from 10 OSCC cell lines, 33 primary OSCC tissues, and 11 normal oral mucosa samples using phenol/chloroform. The bisulfite modification of 2 μg genomic DNA was performed using the EZ DNA Methylation Kit (Zymo Research, CA, USA). For positive and negative controls, in vitro methylated DNA (IVD) and H2O were used. Primer pairs for methylation analyses were preferentially designed for CpG islands of the target genes [21, 34, 37, 48]. The primer sequences are listed in Additional file 2: Table S1. Methylation analyses, including MSP and bisulfite sequencing analysis, were performed as previously described .
One microgram of genomic DNA from each sample was bisulfite converted using the EZ DNA Methylation Kit (Zymo Research, CA, USA) following the manufacturer’s protocol. The PCR conditions and primer sequences are provided in Additional file 2: Table S1. The PCR amplicons were gel-purified and subcloned into the pCRII-TOPO vector (Invitrogen, Carlsbad, CA). At least five to seven clones were randomly selected and sequenced on an ABI3730xl DNA analyzer to ascertain the methylation patterns of each locus.
Quantitative real-time RT-PCR
Total RNA was isolated from the human normal colon and UC patient tissues using TRI-Solution (Bio Science Technology) following the manufacturer’s protocol. RNA quantity was measured using a NanoDrop 2000/2000c instrument (Thermo Scientific, MA, USA), and 1 μg of total RNA was reverse-transcribed into cDNA using the iScriptTMcDNA Synthesis Kit (Bio-Rad). For expression studies, the primers are listed in Additional file 2: Table S1. Quantitative RT-PCR was performed on a CFX96TM Real-Time PCR Detection System (Bio-Rad) using SYBR Green Master Mix (Thermo Scientific, MA, USA). The expression levels of the target genes were normalized to the GAPDH levels, and all the relative quantifications of expression were calculated using the 2−ΔΔCt method.
Formalin-fixed, paraffin-embedded specimens of OSCC and normal oral mucosal tissues were sectioned to a 5-μm thickness. The sections were deparaffinized in xylene and rehydrated through graded alcohol solutions. Heat-induced antigen retrieval was performed in a microwave oven for 10 min in Tris/EDTA buffer (pH 9.0). Endogenous peroxidase activity was inactivated by treating with 3% H2O2 in PBS for 15 min. After blocking nonspecific binding (0.75% BSA in PBS), the sections were incubated with the anti-TFPI2 (Abcam, dilution 1:100), anti-SOX17 (Abcam, dilution 1:100), and anti-GATA4 antibodies (Abcam, dilution 1:100) overnight at 4 °C. The staining was visualized by a peroxidase-conjugated secondary antibody and diaminobenzidine (Vector labs, Burlingame, CA, USA). Finally, the sections were counterstained by Mayer’s hematoxylin and were mounted and photographed with an Axioplan microscope (Carl Zeiss, Germany). Primary antibodies were omitted for the negative controls.
Kaplan–Meier survival analysis using TCGA DNA methylation data
Head and neck squamous cell carcinoma (HNSCC) DNA methylation dataset (HumanMethylation 450 k Illumina; n = 528) in The Cancer Genome Atlas (TCGA) database  was initially downloaded from the firebrowse website (http://firebrowse.org/). Then, samples related to OSCC (n = 344) were selected for survival analysis from HNSCC TCGA cohort (patients (n = 184) from tonsil, larynx, oropharynx, hypopharynx, and lip were excluded based on anatomic subdivision). Patients were divided into two groups (high—hypermethylated, n = 174; low—hypomethylated, n = 174) according to the methylation level (β value) on each CpG site located in the promoter regions of TFPI2 and SOX17, but in gene body region of GATA4. The R package survminer (https://github.com/kassambara/survminer) was used for the analysis. P values were calculated using the log-rank test.
Quantified data are expressed as the mean ± standard deviation (SD). Significance testing was conducted using Student’s t test.
Availability of data and materials
- KM survival curve:
Kaplan–Meier survival curves
Oral squamous cell carcinoma
Reverse transcriptase polymerase chain reaction
Tumor suppressor genes
Transcription start sites
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128(4):683–92.
Herman JG, Baylin SB. Gene silencing in cancer in association with promoter hypermethylation. N Engl J Med. 2003;349(21):2042–54.
Schuebel KE, Chen W, Cope L, Glöckner SC, Suzuki H, Yi JM, et al. Comparing the DNA hypermethylome with gene mutations in human colorectal cancer. PLoS Genetics. 2007;3(9):1709–23.
Esteller M. Epigenetic gene silencing in cancer: the DNA hypermethylome. Hum Mol Genet. 2007;16:R50–9.
Shen L, Kondo Y, Guo Y, Zhang J, Zhang L, Ahmed S, et al. Genome-wide profiling of DNA methylation reveals a class of normally methylated CpG island promoters. PLOS Genetics. 2007;3(10):e181.
Heyn H, Esteller M. DNA methylation profiling in the clinic: applications and challenges. Nat Rev Genet. 2012;13(10):679–92.
Castellarin M, Warren RL, Freeman JD, Dreolini L, Krzywinski M, Strauss J, et al. Fusobacterium nucleatum infection is prevalent in human colorectal carcinoma. Genome Res. 2012;22(2):299–306.
Kostic Aleksandar D, Chun E, Robertson L, Glickman Jonathan N, Gallini Carey A, et al. Fusobacterium nucleatum potentiates intestinal tumorigenesis and modulates the tumor-immune microenvironment. Cell Host Microbe. 2013;14(2):207–15.
Marur S, Forastiere AA. Head and neck squamous cell carcinoma: update on epidemiology, diagnosis, and treatment. Mayo Clin Proc. 2016;91(3):386–96.
Argiris A, Karamouzis MV, Raben D, Ferris RL. Head and neck cancer. Lancet. 2008;371(9625):1695–709.
Scully C, Bagan J. Oral squamous cell carcinoma overview. Oral Oncol. 2009;45(4):301–8.
Scully C, Bagan JV. Recent advances in Oral Oncology 2007: Imaging, treatment and treatment outcomes. Oral Oncol. 2008;44(3):211–5.
Pitiyage G, Tilakaratne WM, Tavassoli M, Warnakulasuriya S. Molecular markers in oral epithelial dysplasia: review. J Oral Pathol Med. 2009;38(10):737–52.
Choi S, Myers JN. Molecular pathogenesis of oral squamous cell carcinoma: implications for therapy. J Dent Res. 2008;87(1):14–32.
The Cancer Genome Atlas. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature. 2015;517:576.
Suzuki H, Watkins DN, Jair KW, Schuebel KE, Markowitz SD, Chen WD, et al. Epigenetic inactivation of SFRP genes allows constitutive WNT signaling in colorectal cancer. Nat Genet. 2004;36(4):417–22.
Jithesh PV, Risk JM, Schache AG, Dhanda J, Lane B, Liloglou T, Shaw RJ. The epigenetic landscape of oral squamous cell carcinoma. Br J Cancer. 2013;108:370.
Flemer B, Warren RD, Barrett MP, et al. The oral microbiota in colorectal cancer is distinctive and predictive. Gut. 2018;67:1454–63.
Tahara T, Yamamoto E, Suzuki H, Maruyama R, Chung W, Garriga J, et al. Fusobacterium in colonic flora and molecular features of colorectal carcinoma. Cancer Res. 2014;74(5):1311.
Mima K, Nishihara R, Qian ZR, Cao Y, Sukawa Y, Nowak JA, et al. Fusobacterium nucleatum in colorectal carcinoma tissue and patient prognosis. Gut. 2016;65(12):1973.
Glockner SC, Dhir M, Yi JM, McGarvey KE, Van Neste L, Louwagie J, et al. Methylation of TFPI2 in stool DNA: a potential novel biomarker for the detection of colorectal cancer. Cancer Res. 2009;69(11):4691–9.
Liu Z, Zhang J, Gao Y, Pei L, Zhou J, Gu L, et al. Large-scale characterization of DNA methylation changes in human gastric carcinomas with and without metastasis. Clin Cancer Res. 2014;20(17):4598.
Hibi K, Goto T, Shirahata A, Saito M, Kigawa G, Nemoto H, et al. Detection of TFPI2 methylation in the serum of colorectal cancer patients. Cancer Lett. 2011;311(1):96–100.
Rasmussen SL, Krarup HB, Sunesen KG, Pedersen IS, Madsen PH, Thorlacius-Ussing O. Hypermethylated DNA as a biomarker for colorectal cancer: a systematic review. Color Dis. 2016;18(6):549–61.
Li Y-F, Hsiao Y-H, Lai Y-H, Chen Y-C, Chen Y-J, Chou J-L, et al. DNA methylation profiles and biomarkers of oral squamous cell carcinoma. Epigenetics. 2015;10(3):229–36.
Sohn J, Natale J, Chew L-J, Belachew S, Cheng Y, Aguirre A, et al. Identification of SOX17 as a transcription factor that regulates oligodendrocyte development. J Neurosci. 2006;26(38):9722.
Matsui T, Kanai-Azuma M, Hara K, Matoba S, Hiramatsu R, Kawakami H, et al. Redundant roles of SOX17 and Sox18 in postnatal angiogenesis in mice. J Cell Sci. 2006;119(17):3513.
Park K-S, Wells JM, Zorn AM, Wert SE, Whitsett JA. SOX17 influences the differentiation of respiratory epithelial cells. Dev Biol. 2006;294(1):192–202.
Kim I, Saunders TL, Morrison SJ. SOX17 dependence distinguishes the transcriptional regulation of fetal from adult hematopoietic stem cells. Cell. 2007;130(3):470–83.
Hulbert A, Jusue-Torres I, Stark A, Chen C, Rodgers K, Lee B, et al. Early detection of lung cancer using DNA promoter hypermethylation in plasma and sputum. Clin Cancer Res. 2017;23(8):1998.
Oishi Y, Watanabe Y, Yoshida Y, Sato Y, Hiraishi T, Oikawa R, et al. Hypermethylation of SOX17 gene is useful as a molecular diagnostic application in early gastric cancer. Tumor Biol. 2012;33(2):383–93.
Jia Y, Yang Y, Liu S, Liu S, Herman JG, Lu F, et al. SOX17 antagonizes WNT/β-catenin signaling pathway in hepatocellular carcinoma. Epigenetics. 2010;5(8):743–9.
Fu D-Y, Wang Z-M, Li C, Wang B-L, Shen Z-Z, Huang W, et al. SOX17, the canonical Wnt antagonist, is epigenetically inactivated by promoter methylation in human breast cancer. Breast Cancer Res Treat. 2010;119(3):601–12.
Zhang W, Glockner SC, Guo M, Machida M, Ota E, Wang DH, et al. Epigenetic inactivation of the canonical Wnt antagonist SRY-box containing gene 17 in colorectal cancer. Cancer Res. 2008;68(8):2764–72.
Clevers H. Wnt/β-catenin signaling in development and disease. Cell. 2006;127(3):469–80.
Gao X, Sedgwick T, Shi YB, Evans T. Distinct functions are implicated for the GATA-4, -5, and -6 transcription factors in the regulation of intestine epithelial cell differentiation. Mol Cell Biol. 1998;18(5):2901–11.
Akiyama Y, Watkins N, Suzuki H, Jair KW, van Engeland M, Esteller M, et al. GATA-4 and GATA-5 transcription factor genes and potential downstream antitumor target genes are epigenetically silenced in colorectal and gastric cancer. Mol Cell Biol. 2003;23(23):8429–39.
Guo M, Akiyama Y, House MG, Hooker CM, Heath E, Gabrielson E, et al. Hypermethylation of the GATA genes in lung cancer. Clin Cancer Res. 2004;10(23):7917–24.
Guo M, House MG, Akiyama Y, Qi Y, Capagna D, Harmon J, Baylin SB, Brock MV, Herman JG. Hypermethylation of the GATA gene family in esophageal cancer. Int J Cancer. 2006;119(9):2078–83.
Wakana K, Akiyama Y, Aso T, Yuasa Y. Involvement of GATA-4/-5 transcription factors in ovarian carcinogenesis. Cancer Lett. 2006;241(2):281–8.
Laird PW. The power and the promise of DNA methylation markers. Nat Rev Cancer. 2003;3(4):253–66.
Pattani KM, Zhang Z, Demokan S, Glazer C, Loyo M, Goodman S, et al. Independently associated with risk of oral cavity cancer and premalignancy. Cancer Prev Res. 2010;3(9):1093.
Schussel J, Zhou XC, Zhang Z, Pattani K, Bermudez F, Jean-Charles G, et al. EDNRB and DCC salivary rinse hypermethylation has a similar performance as expert clinical examination in discrimination of oral cancer/dysplasia versus benign lesions. Clin Cancer Res. 2013;19(12):3268.
Towle R, Truong D, Hogg K, Robinson WP, Poh CF, Garnis C. Global analysis of DNA methylation changes during progression of oral cancer. Oral Oncol. 2013;49(11):1033–42.
Arantes LMRB, de Carvalho AC, Melendez ME, Carvalho AL, Goloni-Bertollo EM. Methylation as a biomarker for head and neck cancer. Oral Oncol. 2014;50(6):587–92.
Foy J-P, Pickering CR, Papadimitrakopoulou VA, Jelinek J, Lin SH, William WN, et al. New DNA methylation markers and global DNA hypomethylation are associated with oral cancer development. Cancer Prev Res. 2015;8(11):1027.
Langevin SM, Butler RA, Eliot M, Pawlita M, Maccani JZJ, McClean MD, et al. Novel DNA methylation targets in oral rinse samples predict survival of patients with oral squamous cell carcinoma. Oral Oncol. 2014;50(11):1072–80.
Yi JM, Dhir M, Guzzetta AA, Iacobuzio-Donahue CA, Heo K, Yang KM, et al. DNA methylation biomarker candidates for early detection of colon cancer. Tumor Biol. 2012;33(2):363–72.
Kim J-G, Kim T-O, Bae J-H, Shim J-W, Kang MJ, Yang K, et al. Epigenetically regulated MIR941 and MIR1247 target gastric cancer cell growth and migration. Epigenetics. 2014;9(7):1018–30.
We would like to thank Prof. Khadijah A. Mitchell of the Lafayette College for the critical reading of the manuscript and providing language editing.
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean Government (MSIT) (No. NRF-2018R1A5A2023879 and No. NRF-2017R1A2B4005588).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board of Pusan National University Dental Hospital (IRB No. PNUDH-2019-006), and written, informed consent was obtained for all study participants prior to data collection.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Correlation between transcriptional expression by qRT-PCR and the promoter hypermethylation of tumor suppressor genes by MSP analysis in OSCC cell lines. Related to Fig. 1a. Figure S2. Protein expression levels of TFPI2, SOX17, and GATA4 in additional set of primary OSCC tumors and normal oral mucosa. Related to Fig. 5. (PPTM 3917 kb)
Table S1. Selected gene primers for MSP, bisulfite sequencing, and RT-PCR analyses. (DOCX 13 kb)