Next Article in Journal
Systematic Review on the Association of Radiomics with Tumor Biological Endpoints
Next Article in Special Issue
Somatic Hypomethylation of Pericentromeric SST1 Repeats and Tetraploidization in Human Colorectal Cancer Cells
Previous Article in Journal
Impact of Planning Method (Conventional versus Virtual) on Time to Therapy Initiation and Resection Margins: A Retrospective Analysis of 104 Immediate Jaw Reconstructions
Previous Article in Special Issue
Follow-Up Study Confirms the Presence of Gastric Cancer DNA Methylation Hallmarks in High-Risk Precursor Lesions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption

by
Sheila Coelho Soares-Lima
1,
Hisham Mehanna
2,
Diego Camuzi
1,
Paulo Thiago de Souza-Santos
3,
Tatiana de Almeida Simão
4,
Pedro Nicolau-Neto
1,
Monique de Souza Almeida Lopes
1,
Cyrille Cuenin
5,
Fazlur Rahman Talukdar
5,
Nikolaos Batis
2,
Izabella Costa
6,
Fernando Dias
6,
Davide Degli Esposti
5,
Mariana Boroni
7,
Zdenko Herceg
5 and
Luis Felipe Ribeiro Pinto
1,4,*
1
Molecular Carcinogenesis Program, Brazilian National Cancer Institute, Rua André Cavalcanti, 37–6° Andar, Bairro de Fátima, Rio de Janeiro 20231-050, Brazil
2
Institute of Head and Neck Studies and Education (InHANSE), Institute of Cancer and Genomic Sciences, University of Birmingham, Birmingham B15 2TT, UK
3
Laboratório de Hanseníase, Instituto Oswaldo Cruz, Fiocruz, Av. Brasil, 4365, Rio de Janeiro 21040-360, Brazil
4
Departamento de Bioquímica, Instituto de Biologia Roberto Alcantara Gomes, Universidade do Estado do Rio de Janeiro, Av. 28 de Setembro 87 fundos, Vila Isabel, Rio de Janeiro 20551-013, Brazil
5
Epigenetics Group, International Agency for Research on Cancer, 150 Cours Albert Thomas, CEDEX 08, 69372 Lyon, France
6
Seção de Cirurgia de Cabeça e Pescoço, Instituto Nacional de Câncer—INCA, Praça da Cruz Vermelha, Rio de Janeiro 20230-130, Brazil
7
Bioinformatics and Computational Biology Lab, Brazilian National Cancer Institute, Rua André Cavalcanti, 37–1° Andar, Bairro de Fátima, Rio de Janeiro 20231-050, Brazil
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(12), 3014; https://doi.org/10.3390/cancers13123014
Submission received: 9 April 2021 / Revised: 27 May 2021 / Accepted: 8 June 2021 / Published: 16 June 2021
(This article belongs to the Special Issue DNA Methylation in Cancer)

Abstract

:

Simple Summary

Squamous cell carcinomas of the upper aerodigestive tract are highly incident, lethal, and share the same epithelial lining of origin, risk factors and genetic alterations. However, their biological and clinical behaviors differ, having an impact on patient survival. This study aimed at identifying the main DNA methylation differences between these tumors, giving an overview of the main genomic regions affected, whether DNA methylation gains or losses are more common, the impact on gene expression and the signaling pathways affected. This knowledge will help identifying potential site-specific biomarkers as well as shedding light on whether epigenetic mechanisms explain, at least in part, the diverse behavior of upper aerodigestive tract tumors.

Abstract

Upper aerodigestive tract (UADT) tumors present different biological behavior and prognosis, suggesting specific molecular mechanisms underlying their development. However, they are rarely considered as single entities (particularly head and neck subsites) and share the most common genetic alterations. Therefore, there is a need for a better understanding of the global DNA methylation differences among UADT tumors. We performed a genome-wide DNA methylation analysis of esophageal (ESCC), laryngeal (LSCC), oral (OSCC) and oropharyngeal (OPSCC) squamous cell carcinomas, and their non-tumor counterparts. The unsupervised analysis showed that non-tumor tissues present markedly distinct DNA methylation profiles, while tumors are highly heterogeneous. Hypomethylation was more frequent in LSCC and OPSCC, while ESCC and OSCC presented mostly hypermethylation, with the latter showing a CpG island overrepresentation. Differentially methylated regions affected genes in 127 signaling pathways, with only 3.1% of these being common among different tumor subsites, but with different genes affected. The WNT signaling pathway, known to be dysregulated in different epithelial tumors, is a frequent hit for DNA methylation and gene expression alterations in ESCC and OPSCC, but mostly for genetic alterations in LSCC and OSCC. UADT tumor subsites present differences in genome-wide methylation regarding their profile, intensity, genomic regions and signaling pathways affected.

1. Introduction

Upper aerodigestive tract (UADT) tumors are highly incident worldwide and their prognosis vary according to tumor subsite and affected population [1,2]. Oral cavity, laryngeal and oropharyngeal tumors, usually grouped together as Head and Neck (HN) cancer, as well as esophageal neoplasms are highly associated with alcohol and tobacco consumption, and usually develop into a squamous cell carcinoma [3]. Furthermore, a significant proportion of patients affected by these cancers develop a second primary tumor, synchronously or metachronously, either in the same subsite or at related anatomic subsites [4], usually resulting in a worse prognosis [5,6]. These cancers are not metastases, but tumor masses that develop independently of the first primary tumor and carry different driver molecular alterations [7]. Despite these similarities, UADT tumors show distinct clinical behavior according to specific subsites, and the identification of their specific characteristics is important to improve patients’ survival.
Oropharyngeal squamous cell carcinoma (OPSCC) is probably the best example of how tumor segregation according to carcinogenesis mechanisms may impact overall survival. Human Papillomavirus (HPV)-positive cases show better clinical outcome compared with HPV-negative cases, and molecular markers, such as p16, are used for discriminating groups and staging OPSCC [8]. Although not yet applied in clinical practice, recent studies have shown that disrupting mutations in NSD1 and NSD2 (Nuclear Receptor Binding SET Domain Proteins 1 and 2) define a group of good prognoses within laryngeal squamous cell carcinoma (LSCC) cases, but not in other HN subsites [9,10]. Therefore, gaining insight into specific molecular alterations present in UADT tumor subsites is of critical importance.
Genome-wide studies have shown that the most common mutational signatures observed in UADT squamous cell carcinomas are those associated with AID/APOBEC (activation-induced cytidine deaminase/apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like) activity [11,12,13], and the most common genetic alteration is TP53 mutations, leading to the inactivation of this tumor suppressor gene [14,15]. Therefore, genetic alterations are in general shared and might not be sufficient to distinguish these tumor subsites. Conversely, DNA methylation alterations are also involved in tumor initiation and progression [16], and aberrant DNA methylation profiles have been shown to be tissue-specific and less heterogeneous than genetic alterations, underscoring their potential as subsite-specific attractive biomarkers [17]. These characteristics together with the reversibility of epigenetic modifications have resulted in an increasing interest in the field. Nevertheless, few studies have compared global methylation profile of UADT tumor subsites, but did not investigate thoroughly subsite-specific alterations, particularly those affecting signaling pathways disruption such as the WNT pathway [14,18,19]. The WNT pathway plays a central role in development and stemness [20,21,22], and its dysregulation in epithelial tumors is recurrent [23,24,25]. Furthermore, WNT signaling pathway disruption was previously shown to impact on cancer patient prognosis and presents the potential of anti-cancer therapeutic approaches targeting this pathway [20,21,22,25].
The present study aimed to compare UADT squamous cell carcinomas subsite DNA methylome changes, pointing out to their main differences, and to identify potential differences among subsites regarding the WNT pathway.

2. Materials and Methods

2.1. Patients

In total, 24 esophageal squamous cell carcinoma (ESCC) patients, 21 LSCC patients, 16 oral squamous cell carcinoma (OSCC) patients and 15 OPSCC patients diagnosed at the Brazilian National Cancer Institute (INCA, Rio de Janeiro, Brazil) were included in the study. Additionally, eight OPSCC patients from the PET-Neck trial (Institute of Head and Neck Studies and Education (InHANSE), University of Birmingham) were also included. Esophageal samples were collected as biopsies through endoscopy procedures, with non-tumor adjacent tissue collected 5 cm from the tumor border. HN tumors and adjacent tissue were collected by the Head and Neck Surgical Division from INCA or from Birmingham University Hospital, from patients who had not undergone chemo- or radiotherapy treatment. For oral cavity and laryngeal sites non-tumor tissue was collected from tumor border free margin sites, selected by a pathologist after patient surgery. For oropharyngeal, the non-tumor tissue consisted of samples collected from tonsillectomies of non-cancer patients. All samples were immediately snap-frozen at liquid nitrogen just after collection (INCA), or formalin-fixed and paraffin embedded (FFPE, PET-Neck). Histopathological profiling of all samples was evaluated by the Pathology Department of INCA, and only tumor samples with >70% of tumor cells were included. Patients’ main characteristics are disclosed in Table S1, but briefly the median age varied between 56–63, most were male (ranging from 87.5% in ESCC to 95.7% in OPSCC), tobacco smokers (smoking at least a cigarette a day for at least one year), alcohol drinkers (drinking alcoholic beverages twice a week for at least one year) and most tumors were diagnosed in late stages (III and IV) in all tumor groups. Patient overall survival (OS) varied according to tumor site (Figure S1), with ESCC showing the lowest median OS (9.1 months), followed by LSCC (20.9 months), and OPSCC (22.1 months). For OSCC, it was not possible to calculate de median OS.
All patients signed informed consents for using biological samples as well as clinical and pathology data from patient records. This study was approved by the institutions Ethics Committees and was conducted according to the Declaration of Helsinki.

2.2. Methylome Analysis

Genomic DNA was extracted from 16 non-tumor adjacent tissues and 24 tumors from ESCC patients; 12 non-tumor adjacent tissues and 20 tumors from LSCC patients; seven non-tumor adjacent tissues and 15 tumors from OSCC patients; as well as 15 OPSCC samples, all snap-frozen, with DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany). All non-tumor adjacent tissues were collected from patients who also donated tumor specimens (matched pairs). Genomic DNA was also extracted from FFPE tonsils from nine non-cancer patients, for whom tobacco and alcohol habits information was not available, and FFPE tumor tissue from eight OPSCC patients with QIAamp DNA FFPE Tissue Kit (Qiagen). Sodium bisulfite-treated DNA (EZ DNA Methylation Kit, Zymo Research, Irvine, CA, USA) was used to assess global DNA methylation profiles by microarray using the Infinium Human Methylation 450K BeadChip (Illumina, San Diego, CA, USA), according to manufacturer’s instructions. In the case of FFPE samples, DNA quality was assessed by the Infinium FFPE DNA QC Kit (Illumina, CA, USA) and only samples that passed the established threshold were used in the following experiments. After sodium bisulfite treatment, restoration was performed using the Infinium FFPE DNA Restoration Kit (Illumina, CA, USA). Microarray experiments were performed in three batches, all containing all sample types, both considering frozen and FFPE samples and tumor types. After checking the built-in-controls with GenomeStudio Software (Illumina, CA, USA), Bioconductor packages were used to perform all analyses in R environment. The idat files were used to obtain average beta-values (AVG_betas) with methylumi [26]. AVG_betas represent the ratio between the intensity of the methylated allele and the intensity of all alleles and vary between 0.0–1.0. Cross-reactive probes and polymorphic probes as described by Chen et al. (2013) [27], samples with 1% of probes with a detection p-value > 0.05, probes with beadcount <3 in 5% of samples as well as probes having 1% of samples with a detection p-value > 0.05 were removed. After applying these filters, one FFPE sample and one adjacent non-tumor samples from a LSCC patient were removed from the analysis (had more than 1% of probes with a detection p-value > 0.05). Finally, Singular Value Decomposition was run to estimate the impact of batch effects and ComBat was run to correct for batch effects, both as functions of the Bioconductor package ‘ChAMP’. Furthermore, color bias adjustment (lumi [28]) and probe bias correction with BMIQ method (watermelon [29]) were performed. In all comparisons, moderated t-statistic was used, and p-values were adjusted for multiple testing by Benjamini and Hochberg’s method. Probes with an adjusted p-value < 0.001 were considered differentially methylated between groups (limma [30]).
DMRs were identified with the methyAnalysis package [31] using the following criteria: false discovery rate (FDR) <0.05, a minimum of 5 probes within the DMR, a maximum gap of 300 bp between two nearby probes to be considered within a same DMR, and DMRs with a gap ≤100 bp were merged.
The 1000 most variable CpG positions were used to calculate the Euclidean distances between samples (separately for non-tumor and tumor samples) that were projected into a 2-d plane using classical multidimensional scaling transformation (minfi [32]).
Genome-wide mean delta beta from each tumor type were used as input to construct heatmaps with the RCircos package.
Over-representation analyses were performed with WEB-based GEne SeT AnaLysis Toolkit [33] using Kyoto Encyclopedia of Genes and Genomes (KEGG) database. For this analysis, genes containing DMRs in each tumor type were used separately and pathways were considered enriched when FDR was less than 0.05.
Information regarding WNT differentially methylated genes by position (promoter and outside promoter) were plotted with PathVisio [34]. When more than one DMR was differentially methylated in a given region, the mean delta beta of all DMRs was represented.
The datasets used during the current study are deposited in the Gene Expression Omnibus (GEO) database, under the accession numbers: GSE178212, GSE178216, GSE178218, GSE178219.

2.3. Cellular Component Analysis

R package MethylCIBERSORT (v0.2.0) [35] was used to construct the stromal and immune cells signature matrix. Methylation beta-values from stromal and immune cells built-in the package were used. Methylation beta-values normalization was performed as already mentioned through the BMIQ method. 100 DMPs with median beta-value differences > 0.25 and FDR < 0.01 were considered in the cellular component signature. Deconvolution was performed in CIBERSORTx [36] server with quantile normalization, in relative mode and 1000 permutations.

2.4. RNA-Seq Analysis

Gene expression data from ESCC and LSCC cases was retrieved from previous work by the group [23]. Briefly, total RNA was isolated with the RNeasy Kit (Qiagen) and integrity was assessed with RNA 6000 Nano chip (Agilent, Santa Clara, CA, USA). RNA samples with RNA integrity number (RIN) ≥ 8 were used for constructing cDNA libraries (TruSeq RNA, Illumina), which were sequenced in an Illumina HiSeq 2500 platform. Data processing and read mapping were performed as previously described [23]. Differential expression analysis was performed with DESeq2 package [37].

2.5. TCGA Data

TCGA gene expression data were retrieved from cBioportal for cancer genomics [38] and DNA methylation data from TCGA Wanderer [39], and used to calculate correlations for selected genes. The DMR beta-value was defined for each sample by calculating the mean of all probes within the DMR. A total of 81 ESCC, 94 LSCC, 233 OSCC and 69 OPSCC samples were included in this analysis.
Information regarding copy number alterations (deletions) and the presence of putative driver mutations in TAP1, AJUBA and NOTCH1 for each tumor type was retrieved from cBioportal for cancer genomics [38]. Data from a total of 96 ESCC, 116 LSCC, 248 OSCC and 79 OPSCC samples were included in this analysis.
The datasets analyzed during the current study are available in the cBioPortal for Cancer Genomics repository [cbioportal.org] [Esophageal Carcinoma (TCGA, Firehose Legacy) and Head and Neck Squamous Cell Carcinoma (TCGA, Firehose Legacy)].

3. Results

3.1. Global DNA Methylation Profile Is Quite Distinct among UADT Squamous Cell Carcinomas

The study workflow is shown in Figure S2. Unsupervised analyses revealed differences in the global DNA methylation profile of non-tumor tissue subsites, with oropharynx presenting the most distinct profile (Figure 1A). Among tumors, a more heterogeneous profile was observed, with no clear stratification by subsite in the unsupervised analysis (Figure 1B). This observation led us to perform a stromal and immune cell component analysis by tumor type, based on DNA methylation levels (Figure S3). Fibroblasts and endothelial cells contributed to ~50% of this non-tumor component in most of the samples. The contribution of the immune components was more variable, with ESCC showing a higher percentage of neutrophils, for example.
The supervised comparison revealed a total of 76,210 differentially methylated probes (DMPs) in ESCC, 42,907 in OPSCC, 18,389 in LSCC, and 9978 in OSCC, when compared to their respective non-tumor tissues (p-value adjusted by Benjamini and Hochberg’s method—BH < 0.001). Most DMPs were hypermethylated in ESCC and in OSCC (67.6% and 71.83%, respectively), but hypomethylated in LSCC and OPSCC (71.73% and 64.5%, respectively), profile also observed when probes were stratified by gene region (Figure 1C). This stratification revealed that most DMPs encompassed gene bodies in all tumor sites (Figure 1C), as expected by the higher coverage of this region in the beadchip. Nevertheless, they were also more affected than expected by chance (p < 0.001) (Figure 1C, Table S2). In ESCC, LSCC and OPSCC, 3′UTR regions were also more frequently affected than by chance (p < 0.001) (Figure 1C, Table S2). Promoter regions (TSS200 and TSS1500) were less affected than expected by chance in general. DMPs mapping to the first exon or 5′UTR region of genes were less affected than would be expected by the representation in the beadchip in all tumor subsites except for OSCC (Figure 1C, Table S2).
CpG islands were more commonly hypermethylated in all subsites (Figure 1D), but DMPs mapped to these regions less frequently than expected by chance (p < 0.001) in ESCC, LSCC and OPSCC, while in OSCC they were more affected than expected (Table S3, Figure 1D). On the opposite direction, DMPs mapped to open sea regions more frequently than expected in ESCC, LSCC and OPSCC, while in OSCC they were less frequently affected (Table S3, Figure 1D).
Figure 1E shows the genome-wide differential methylation profile for each tumor subsite, evidencing that OSCC shows the highest delta betas (median of –0.28 and 0.29 for hypo and hypermethylated probes, respectively), followed by LSCC (hypo: −0.26 and hyper: 0.27), OPSCC (hypo: −0.21 and hyper: 0.26), and ESCC, which showed the lowest differences in comparison with its non-tumor tissue (hypo: −0.15 and hyper: 0.20) (Figure S4A). Since this could be a consequence of the different number of samples (non-tumor and tumor tissues) used for each tumor subsite, we analyzed this profile only considering the common DMPs, and ESCC indeed showed the lowest median delta betas (Figure S4B). Next, we evaluated whether the mean delta betas could reflect differences within non-tumor tissues or tumors. Figure S4C shows the methylation profile of 1114 common DMPs in these tissues. Taking ESCC as an example, for hypomethylated DMPs, it showed the lowest methylation levels among non-tumor tissues and the highest among tumors, explaining the low delta betas. Similarly, for hypermethylated DMPs, both the non-tumor tissues and tumors’ methylation profiles explain the delta beta differences. ESCC showed the highest methylation levels in NTST and intermediate levels among tumors, again corroborating the low delta betas.
Tumors also showed a different profile of over or underrepresentation of affected chromosomes, with no specific chromosome being commonly affected in the four subsites (Table S4). However, head and neck tumor subsites (LSCC, OSCC and OPSCC) shared an underrepresentation of DNA methylation alterations in chromosomes 9, 17 and 22 and an overrepresentation of chromosome 14.

3.2. Differentially Methylated Regions Affect Different Cell Pathways and Different Genes in the Same Cell Signaling Pathways, According to Tumor Subsite

The analysis of differentially methylated regions (DMRs, FDR < 0.05 and CpG sites ≥ 5) corroborated the global methylation profile assessed by DMPs. Hypermethylation was more common in ESCC and OSCC (61.8% of 5833 DMRs and 72.4% of 2145 DMRs, respectively), while LSCC and OPSCC presented mostly hypomethylation (63.2% of 2643 DMRs and 78.5% of 7676 DMRs, respectively). For LSCC and OSCC, the methylation levels of the same DMRs were assessed in TCGA dataset, showing significant differences between tumor and NTST for 98% and 96% of the DMRs, respectively (BH adjusted p-value < 0.05). Only 12 out of a total 18,261 DMRs considered (0.066%) were common among the four subsites.
Next, over-representation analysis (ORA) was performed to identify the KEGG (Kyoto Encyclopedia of Genes and Genomes) signaling pathways enriched for genes affected by DMRs in each tumor subsite (FDR < 0.05, Table S5). Figure 2 shows that ESCC presented the highest number of affected pathways (57), followed by OPSCC (49), OSCC (22), and LSCC (11). The overlap of enriched pathways evidenced that only four out of 127 affect pathways (3.1%) were common among the four subsites. These were cAMP Signaling Pathway, Calcium Signaling Pathway, Glutamatergic Synapse, and Neuroactive Ligand-Receptor Interaction (Figure 2). Nevertheless, a proportion of the genes affected in these four pathways were different among the subsites (Figure 2, Table S5).

3.3. Methylation Levels of DMRs in WNT Signaling Genes Are Potentially Associated with Gene Expression

WNT signaling pathway was found enriched by genes carrying DMRs in ESCC and OPSCC, but not in LSCC and OSCC. Abnormal activation of this pathway has been constantly associated with epithelial tumor development and prognosis [22,23,24,25]. Therefore, we performed a deeper analysis of WNT-associated genes in these tumors.
From the 150 genes included in the WNT signaling pathway map from KEGG database, 51 (34%, with 64 DMRs) were affected by DMRs in OPSCC, 42 (28%, with 57 DMRs) in ESCC, but only 11 (7.3%, with 14 DMRs) in OSCC, and 9 (6%, with 12 DMRs) in LSCC (Figure 3A, Table S6). DMRs in genes involved in the WNT pathway were classified according to their location (promoter and outside promoter) and a representation of the pathway containing this information was designed (Figure 4). To understand whether these DNA methylation alterations would have a functional impact, we assessed gene expression in ESCC and LSCC in our dataset, representing tumor subsites with a high and low number of affected genes, respectively. In ESCC, 24 out of 42 genes carrying DMRs (57.1%) were differentially expressed (Figure S5). In LSCC, only three out of nine DMR-affected genes showed differential expression (30%), (Figure S6).
To validate the potential impact of DMRs identified here on gene expression, we used the TCGA datasets for each tumor subsite included in the present study. In this analysis, the correlation between the mean methylation value of all CpG sites included in the DMR and the mRNA expression of the respective gene was assessed. WNT signaling genes carrying more than one DMR in a specific tumor subsite had the correlation calculated separately for each DMR. The methylation levels of at least one DMR was significantly correlated with mRNA expression (BH adjusted p-value < 0.05) in 81.8% of the genes in OSCC, 40.5% in ESCC, 39.2% in OPSCC, and only 12.5% in LSCC (Figure 3B, Table S6).
Although a relatively low number of genes was affected by DMRs in OSCC, most of these DMRs were significantly correlated with expression in TCGA dataset (Figure 3B). The pathway inhibitors DKK2, SFRP2, SFRP4 and SOX17 showed hypermethylated DMRs in the promoter region that were inversely correlated with mRNA expression (Figure 3B, Table S6), suggesting their downregulation in OSCC. SOX17 also carried a hypermethylated DMR outside the promoter region inversely correlated with expression.
In ESCC, only DMRs outside promoters were significantly correlated with the expression of WNT ligands (Figure 3B). Of note, WNT7A was found hypomethylated and overexpressed in our dataset, and TCGA data evidenced an inverse correlation between methylation and expression (Figure 3B, Table S6). Considering WNT receptors, FZD7 showed a hypomethylated DMR outside the promoter region and mRNA overexpression in ESCC respective to NTST in our dataset, and its methylation levels were inversely correlated with expression in TCGA. The downregulation and hypermethylation (outside promoter) of CTNNBIP1, a pathway inhibitor, was detected in our dataset, followed by an inverse correlation between DMR methylation and mRNA expression in TCGA (Figure 3B, Table S6).
OPSCC showed DMRs both in the promoters and outside the promoter region of WNT signaling genes. However, the latter were more commonly correlated with gene expression (Figure 3B). Regarding WNT ligands, DMRs in WNT5A and WNT7A were hypomethylated and inversely correlated with gene expression (Figure 3B, Table S6). WNT3A was hypermethylated and its DMR methylation levels were positively correlated with expression. In all these cases, an overexpression of the ligands would be expected to take place. Conversely, WNT11 DMR was hypomethylated and positively correlated with expression, suggesting its downregulation.
In LSCC, the methylation levels of only one DMR located outside the promoter region of FZD10 was significantly and positively correlated with mRNA expression in TCGA dataset (Figure 3B). Since this DMR was hypomethylated in tumors in our dataset (Table S6), a gene upregulation would be expected. However, our RNA-seq data showed no significant differences (Figure S6).

3.4. WNT Signaling Pathway May Be Affected by Different Molecular Mechanisms in UADT Subsites

WNT signaling activation in HN tumor subsites has been associated with genetic alterations in AJUBA, FAT1 and NOTCH1 [24]. We evaluated whether the frequency of genetic alterations would be higher in the subsites that presented a lower proportion of WNT genes affected by DMRs. With this purpose, we evaluated the proportion of cases of each tumor subsite carrying potential driver mutations and/or deep deletions in AJUBA, FAT1 and NOTCH1 in TCGA datasets. The frequency of genetic alterations in either of the three genes was 43.1% in LSCC, 35.9% in OSCC, 15.6% in ESCC, and 10.1% in OPSCC (Figure 3C).
Finally, we checked whether the WNT signaling pathway was activated in ESCC, LSCC and OSCC, although the associated molecular mechanisms seem to differ. From a list of 49 direct transcriptional targets of the pathway in mammals (the WNT homepage [40]), 32 (65.3%) were differentially expressed (BH adjusted p-value < 0.05) in ESCC (being 28 up- and 4 downregulated), 29 (59.2%) in OSCC (being 22 up- and 7 downregulated) and 23 (46.9%) in LSCC (being 22 up- and 1 downregulated) (Figure S7), suggesting an activation of the pathway in the evaluated tumor subsites.

4. Discussion

There is limited information regarding the comprehension of the specific epigenetic alterations involved in UADT tumor subsites [41,42,43,44]. In this study we showed that four subsites of UADT tumors analyzed in this study, ESCC, OSCC, OPSCC and LSCC, present marked differences among them regarding general profile, intensity, and genomic region of methylation changes, resulting in different genes and pathways affected, exemplified by the WNT pathway alteration analysis in the four subsites.
We showed that the general DNA methylation profile differs among the non-tumor epithelium of the different UADT subsites, evidencing that, although they are lined by a common squamous mucosa, tissue-specific phenotypes can be noticed even before neoplastic transformation. However, we cannot rule out that at least some of the observed differences can be a consequence of different sample source. For OPSCC, opposite to the other tumor subsites, for which fresh-frozen NTST was collected, the non-tumor tissue consisted of FFPE tonsils from individuals going through uvulopalatopharyngoplasty for conditions not associated with cancer. Tonsils were used based on previous studies and on the high prevalence of OPSCC arising from these organs [45,46,47,48,49,50,51], but future studies on OPSCC NTST should be performed to confirm our findings.
When considering all DMPs, ESCC and OSCC presented a clear hypermethylation, whereas LSCC and OPSCC showed a hypomethylation profile. Poage et al. (2011) [41] compared the global methylation profile of HN tumors and showed that LINE-1 and 27k Illumina Infinium Methylation Beadchip median methylation levels were higher in oral, followed by pharynx and larynx squamous cell carcinomas. In a multivariable regression, AluYB8 methylation levels were also lower in LSCC relative to OSCC [41]. These data are in accordance with our findings, but the authors directly compared tumors, while we compared the differences relative to non-tumor counterparts. In a similar experimental design to ours, LSCC and OSCC were shown to be mostly hypomethylated, while OPSCC was mostly hypermethylated relative to normal adjacent tissue [42]. Apart from LSCC, our results show the opposite profile, what could be associated with the different platforms used in the two studies, 27 K by Lleras and colleagues and 450 K by us. Nevertheless, from the 25 genes carrying hypermethylated DMPs in all tumor subsites in the cited study, 17 (68%) were also commonly hypermethylated in our dataset, including ESCC.
TCGA network has proposed a classification for HNSCC based on DNA methylation, in which tumors could be divided in hypermethylated, hypomethylated, normal-like and CpG island methylated [14]. More specifically, the hypomethylated subgroup was enriched for NSD1-mutated tumors. Since NSD1 mutations have been shown to be enriched in LSCC [52], this genetic alteration could contribute to the establishment of a hypomethylated phenotype in these tumors [14]. Methylation profiles have also been used before to evaluate and subdivide tumors from the same histology within TCGA datasets [43]. Squamous cell carcinomas, including lung, cervical, bladder, head and neck and esophagus, were classified in 5 clusters, based on the methylation profile of 905 differentially methylated and differentially expressed genes in at least one tumor type. Cluster 4, hypermethylated, and cluster 5, hypomethylated, were mostly composed by HNSCC. In cluster 5, a higher proportion of NSD1 mutations was found, what could explain the hypomethylated profile, and suggests a higher representation of LSCC. The hypermethylation profile observed in cluster 4 was associated with CASP8 and HRAS mutations. A higher frequency of these mutations has been described before in OSCC relative to other HNSCC subsites [52]. This suggests that cluster 4 might be enriched for OSCC, corroborating the hypermethylation profile described here.
OSCC clearly presented a unique hypermethylation profile, particularly enriched in CpG islands. A CpG island methylator phenotype (CIMP) was previously reported in HNSCC and was indeed more frequent among OSCC patients [53]. However, this group of patients, called CIMP-Atypical, was enriched for females and never or former smokers (>15 years without smoking), differing from our dataset. The patients included in our study reflect the general OSCC patient profile, being mostly males (93.3%) and current smokers (>80%) [54]. Furthermore, OSCC patient profile was not different from patients with tumors at other subsites in our study. Therefore, etiology might not be the cause of CIMP. Brennan K. and colleagues also associated the hypermethylation phenotype with mutations in CASP8, HRAS and NOTCH1 [53], the latter being more commonly, but not exclusively, mutated in OSCC in our analyses.
Promoter hypermethylation of 11 tumor suppressor genes has been previously correlated with EZH2 expression in HNSCC [55]. This histone methyltransferase is part of the Polycomb complex and catalyzes H3K27 trimethylation, a mark of gene silencing [56]. The association of EZH2 overexpression with both gene-specific and global hypermethylation has been described before for prostate cancer, small cell lung carcinoma and cell sarcoma of the kidney, suggesting the role of this histone modifier in directing hypermethylation [57,58,59]. Although such association has not been described in OSCC, mostly hypermethylated according to our results, EZH2 higher expression was shown to be an independent factor for OSCC development within oral leukoplakia patients [60] and associated with radioresistance [61]. Interestingly, global hypermethylation has been observed in radioresistant HNSCC before [62]. ESCC, the other hypermethylated UADT subsite described here, also presents EZH2 overexpressed [63] and associated with a poor prognosis [64]. These data suggest a role for EZH2 in OSCC and ESCC and, although speculative, this might be associated, at least in part with the hypermethylation profile observed here, and future studies should test this hypothesis.
Another surprising difference among tumor subsites evaluated was the variable median delta betas, both considering hypo and hyper DMPs, with ESCC showing the lowest differences. This could indicate that esophageal mucosa shows the most distinct methylation profile when considering the other UADT linings. This may be a consequence of inherent molecular characteristics of the esophageal squamous epithelium but is also clearly caused by subsite specific neoplastic transformation. The methylation profile of normal tissues from cancer patients have been shown to differ and this might be associated with cancer development [65].
We observed a difference in DMRs similar to that seen for DMPs. Genes affected by DMRs resulted in a total of 127 pathways potentially altered, but only 3.1% (four) of these were common to all subsites. Even among the low proportion of similar target pathways, the genes affected were in most cases different. These data indicate that epigenetic alterations affect distinct pathways and genes among different UADT tumor subsites. One of the differentially affected pathways was the WNT signaling pathway. The WNT signaling pathway activation has been recurrently reported in UADT tumors from preneoplastic lesions to more advanced tumors, being associated with patients’ outcome [66,67,68,69,70,71,72,73,74,75]. This pathway is crucial for the development of the epithelial linings, but in adult epithelial cells it becomes silent, and its reactivation seems to play a driver role in neoplastic transformation [71,73,74]. Studies using in vitro and in vivo models have shown that tobacco-mimicking compounds and HPV infection can induce the activation of WNT pathway [76,77,78,79]. However, the molecular mechanisms that lead to its activation are not so clear. Genetic alterations in WNT-activating genes are rare in UADT tumors, while inactivating mutations in non-canonical genes (FAT1, AJUBA and NOTCH1) might be associated with the disruption of the WNT pathway in a subset of UADTs, leading to beta-catenin stabilization [24].
However, here we show that not only the molecular mechanisms vary depending on subsite, but also the genes and the affected regulatory regions can differ. LSCC and OSCC showed a higher frequency of disrupting genetic alterations in the so-called non-canonical WNT genes (43.1% and 35.9%, respectively), AJUBA, FAT1 or NOTCH1, and a relatively small number of WNT genes affected by DMRs. LSCC had only one gene (FZD10) with a correlation between methylation and expression, suggesting a small contribution of DNA methylation dysregulation of WNT genes to pathway disruption. OSCC also showed a relatively small number of WNT genes affected by DMRs, but their methylation levels were significantly correlated with mRNA expression in most cases. This profile was mainly observed for WNT pathway inhibitors and corroborated previous studies that showed SFRP2, SFRP4, SFRP5, WIF1, DKK3 and DAB2 (all WNT antagonists) promoter hypermethylation [69,80,81]. In addition, for SOX17, recently shown to be epigenetically silenced in OSCC, with its hypermethylation associated with a poor prognosis [66], we detected DMRs both in the promoter and outside the promoter region, both associated with gene expression. This suggests the contribution of both epigenetic and genetic alterations activating WNT pathway in OSCC.
Genetic alterations in the non-canonical WNT genes are present in only 15.6% of ESCC and 10.1% of OPSCC, and epigenetic mechanisms seem to be at least one of the causes of WNT signaling disruption in these tumor subsites. Here we identified DMRs in all classes of WNT genes, including ligands, receptors, activators, inhibitors, targets, and genes involved in WNT/Ca2+ signaling and planar cell polarity (PCP) pathway. While few studies have shown such association for OPSCC before [69,82,83], most data available for ESCC point to a role of aberrant DNA methylation in the silencing of WNT negative regulators [84,85,86,87,88]. Therefore, our study suggests a broader range of genes affected by this epigenetic mechanism that may go beyond promoters.
The data presented here show a quite complex tumor specific subsite methylome landscape. Although we applied stringent filters and quality control approaches, the different tissue source for OPSCC (fresh-frozen and FFPE) was a limitation of the present work and should be addressed in future studies. Nevertheless, overall DNA methylation profile in tumors was highly heterogeneous, corroborating previous studies [13,38,39]. This might be associated both with tumor cell epigenetic clonality and variations of the stromal and immune cell component, as shown in Figure S3 and previously assessed by other authors [89,90]. Indeed, the immune infiltrate, usually assessed by gene expression levels, was already shown to vary among ESCC and HNSCC patients [91,92,93] and might explain, at least in part, tumor heterogeneity in terms of DNA methylation profile. The number of samples analyzed was limited and we were not able to evaluate the impact of etiological factors on the DNA methylation profiles identified due to the quite homogeneous characteristics of patients, mostly heavy smokers, heavy drinkers and HPV-negative. However, it is important to mention that this is the actual profile of UADT patients in general [51,54,94,95]. Therefore, the molecular heterogeneity within these tumor types might be even more pronounced. Defining these profiles may indicate new potential diagnostic and prognostic biomarkers, and point to new therapeutic routes for tumors with very poor outcomes.

5. Conclusions

The present study gives a comprehensive overview of the differences in the DNA methylation profile of ESCC, LSCC, OSCC and OPSCC, showing that these UADT subsite tumors present differences in profile, intensity and regions affected by methylation alterations. The WNT signaling pathway, constantly activated in epithelial tumors, was differentially affected among UADT subsites, presenting methylation alterations particularly in ESCC and OPSCC, whereas genetic alterations capable of activating it are more pronounced in LSCC and OSCC.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13123014/s1, Table S1: Characteristics of the patients included in the study, Figure S1: Patient overall survival according to tumor site, Figure S2: Study workflow, Figure S3: Stromal and immune cell contribution in ESCC and HNSCC, Table S2: Differentially methylated probes stratified by their gene position in ESCC, LSCC, OSCC and OPSCC, Table S3: Differentially methylated probes stratified by their relation to CpG island in ESCC, LSCC, OSCC and OPSCC, Figure S4: DMP delta beta differences reflect non-tumor and tumor differences, Table S4: Differentially methylated probes stratified by chromosomes in ESCC, LSCC, OSCC and OPSCC, Table S5: Over-representation (ORA) analysis of genes carrying DMRs in ESCC, LSCC, OSCC and OPSCC, Table S6: DMR description and correlation with the expression of the associated gene in ESCC, LSCC, OSCC and OPSCC, Figure S5: Gene expression profile of WNT genes carrying DMRs in ESCC, Figure S6: Gene expression profile of WNT genes carrying DMRs in LSCC, Figure S7: WNT target genes are upregulated in ESCC, LSCC and OSCC.

Author Contributions

Conceptualization, S.C.S.-L., H.M., Z.H. and L.F.R.P.; data curation, S.C.S.-L., D.C., P.N.-N., P.T.d.S.-S., M.B.; formal analysis, S.C.S.-L., D.C., P.T.d.S.-S., M.B.; investigation, S.C.S.-L., D.C., P.N.-N., P.T.d.S.-S.; resources, S.C.S.-L., D.C., T.d.A.S., P.N.-N., M.d.S.A.L., C.C., I.C., F.D., F.R.T., N.B., D.D.E., M.B.; software, S.C.S.-L., D.C., P.T.d.S.-S., M.B.; methodology, S.C.S.-L., D.C., T.d.A.S., P.N.-N., M.d.S.A.L., C.C., F.R.T., N.B., D.D.E.; validation, S.C.S.-L., D.C., T.d.A.S., P.N.-N., P.T.d.S.-S.; visualization, S.C.S.-L., D.C.; writing—original draft preparation, S.C.S.-L.; writing—review and editing, H.M., D.C., T.d.A.S., P.T.d.S.-S., F.R.T., M.B., Z.H., L.F.R.P.; supervision, Z.H., L.F.R.P.; project administration, S.C.S.-L., H.M., Z.H., L.F.R.P.; funding acquisition, Z.H., L.F.R.P. Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy or views of the International Agency for Research on Cancer/World Health Organization. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Swiss Bridge Award, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro (E-26/010.00185/2015, E-26/202.563/2019, E-26/202.918/2019), Conselho Nacional de Desenvolvimento Científico e Tecnológico (407992/2016-2) and Ministério da Saúde.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki. Ethical approval for the PET-Neck clinical trial was provided by the Oxfordshire Multi-Research Ethics Committee in May 2007 (Ref No: 07/Q1604/35). Ethical approval was also obtained from the Research Ethics Committee from the Brazilian National Cancer Institute in December 2011 (Ref No: CAAE–0086.0.007.000-11).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets used during the current study are deposited in the Gene Expression Omnibus (GEO) database, under the accession numbers: GSE178212, GSE178216, GSE178218, GSE178219.

Acknowledgments

The authors would like to thank the National Tumor Bank, the Endoscopy Section and Head and Neck Section, all from the Brazilian National Cancer Institute (INCA), for all the support with patient enrollment and sample collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [Green Version]
  2. Rettig, E.M.; D’Souza, G. Epidemiology of head and neck cancer. Surg. Oncol. Clin. N. Am. 2015, 24, 379–396. [Google Scholar] [CrossRef]
  3. Prabhu, A.; Obi, K.O.; Rubenstein, J.H. The synergistic effects of alcohol and tobacco consumption on the risk of esophageal squamous cell carcinoma: A meta-analysis. Am. J. Gastroenterol. 2014, 109, 822–827. [Google Scholar] [CrossRef] [PubMed]
  4. Priante, A.V.; Castilho, E.C.; Kowalski, L.P. Second primary tumors in patients with head and neck cancer. Curr. Oncol. Rep. 2011, 13, 132–137. [Google Scholar] [CrossRef] [PubMed]
  5. Lampri, E.S.; Chondrogiannis, G.; Ioachim, E.; Varouktsi, A.; Mitselou, A.; Galani, A.; Briassoulis, E.; Kanavaros, P.; Galani, V. Biomarkers of head and neck cancer, tools or a gordian knot? Int. J. Clin. Exp. Med. 2015, 8, 10340–10357. [Google Scholar] [PubMed]
  6. Lim, H.; Kim, D.H.; Jung, H.Y.; Gong, E.J.; Na, H.K.; Ahn, J.Y.; Kim, M.Y.; Lee, J.H.; Choi, K.S.; Choi, K.D.; et al. Clinical significance of early detection of esophageal cancer in patients with head and neck cancer. Gut Liver 2015, 9, 159–165. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Leemans, C.R.; Braakhuis, B.J.; Brakenhoff, R.H. The molecular biology of head and neck cancer. Nat. Rev. Cancer 2011, 11, 9–22. [Google Scholar] [CrossRef]
  8. Amin, M.B.; Greene, F.L.; Edge, S.B.; Compton, C.C.; Gershenwald, J.E.; Brookland, R.K.; Meyer, L.; Gress, D.M.; Byrd, D.R.; Winchester, D.P. The Eighth Edition AJCC Cancer Staging Manual: Continuing to build a bridge from a population-based to a more “personalized” approach to cancer staging. CA Cancer J. Clin. 2017, 67, 93–99. [Google Scholar] [CrossRef]
  9. Peri, S.; Izumchenko, E.; Schubert, A.D.; Slifker, M.J.; Ruth, K.; Serebriiskii, I.G.; Guo, T.; Burtness, B.A.; Mehra, R.; Ross, E.A.; et al. NSD1- and NSD2-damaging mutations define a subset of laryngeal tumors with favorable prognosis. Nat. Commun. 2017, 8, 1772. [Google Scholar] [CrossRef] [PubMed]
  10. Pan, C.; Izreig, S.; Yarbrough, W.G.; Issaeva, N. NSD1 mutations by HPV status in head and neck cancer: Differences in survival and response to DNA-damaging agents. Cancers Head Neck 2019, 4, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Roberts, S.A.; Lawrence, M.S.; Klimczak, L.J.; Grimm, S.A.; Fargo, D.; Stojanov, P.; Kiezun, A.; Kryukov, G.V.; Carter, S.L.; Saksena, G.; et al. An APOBEC cytidine deaminase mutagenesis pattern is widespread in human cancers. Nat. Genet. 2013, 45, 970–976. [Google Scholar] [CrossRef]
  12. Alexandrov, L.B.; Nik-Zainal, S.; Wedge, D.C.; Aparicio, S.A.; Behjati, S.; Biankin, A.V.; Bignell, G.R.; Bolli, N.; Borg, A.; Børresen-Dale, A.L.; et al. Signatures of mutational processes in human cancer. Nature 2013, 500, 415–421. [Google Scholar] [CrossRef] [Green Version]
  13. Sawada, G.; Niida, A.; Uchi, R.; Hirata, H.; Shimamura, T.; Suzuki, Y.; Shiraishi, Y.; Chiba, K.; Imoto, S.; Takahashi, Y.; et al. Genomic Landscape of Esophageal Squamous Cell Carcinoma in a Japanese Population. Gastroenterology 2016, 150, 1171–1182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Network, C.G.A. Comprehensive genomic characterization of head and neck squamous cell carcinomas. Nature 2015, 517, 576–582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Network, C.G.A.R. Integrated genomic characterization of oesophageal carcinoma. Nature 2017, 541, 169–175. [Google Scholar]
  16. Jones, P.A.; Baylin, S.B. The fundamental role of epigenetic events in cancer. Nat. Rev. Genet. 2002, 3, 415–428. [Google Scholar] [CrossRef]
  17. Koch, A.; Joosten, S.C.; Feng, Z.; de Ruijter, T.C.; Draht, M.X.; Melotte, V.; Smits, K.M.; Veeck, J.; Herman, J.G.; Van Neste, L.; et al. Analysis of DNA methylation in cancer: Location revisited. Nat. Rev. Clin. Oncol. 2018, 15, 459–466. [Google Scholar] [CrossRef]
  18. Zhou, C.; Ye, M.; Ni, S.; Li, Q.; Ye, D.; Li, J.; Shen, Z.; Deng, H. DNA methylation biomarkers for head and neck squamous cell carcinoma. Epigenetics 2018, 13, 398–409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Bai, G.; Song, J.; Yuan, Y.; Chen, Z.; Tian, Y.; Yin, X.; Niu, Y.; Liu, J. Systematic analysis of differentially methylated expressed genes and site-specific methylation as potential prognostic markers in head and neck cancer. J. Cell Physiol. 2019, 234, 22687–22702. [Google Scholar] [CrossRef] [Green Version]
  20. Taciak, B.; Pruszynska, I.; Kiraga, L.; Bialasek, M.; Krol, M. Wnt signaling pathway in development and cancer. J. Physiol. Pharmacol. 2018, 69. [Google Scholar] [CrossRef]
  21. Zhan, T.; Rindtorff, N.; Boutros, M. Wnt signaling in cancer. Oncogene 2017, 36, 1461–1473. [Google Scholar] [CrossRef]
  22. Tan, S.H.; Barker, N. Wnt Signaling in Adult Epithelial Stem Cells and Cancer. Prog. Mol. Biol. Transl. Sci. 2018, 153, 21–79. [Google Scholar] [CrossRef] [PubMed]
  23. Souza-Santos, P.T.; Soares Lima, S.C.; Nicolau-Neto, P.; Boroni, M.; Meireles Da Costa, N.; Brewer, L.; Menezes, A.N.; Furtado, C.; Moreira, M.A.M.; Seuanez, H.N.; et al. Mutations, Differential Gene Expression, and Chimeric Transcripts in Esophageal Squamous Cell Carcinoma Show High Heterogeneity. Transl. Oncol. 2018, 11, 1283–1291. [Google Scholar] [CrossRef]
  24. Leemans, C.R.; Snijders, P.J.; Brakenhoff, R.H. The molecular landscape of head and neck cancer. Nat. Rev. Cancer 2018, 18, 269. [Google Scholar] [CrossRef]
  25. Han, Q.; Wang, X.; Liao, X.; Han, C.; Yu, T.; Yang, C.; Li, G.; Han, B.; Huang, K.; Zhu, G.; et al. Diagnostic and prognostic value of WNT family gene expression in hepatitis B virus-related hepatocellular carcinoma. Oncol. Rep. 2019, 42, 895–910. [Google Scholar] [CrossRef]
  26. Davis, S.; Du, P.; Bilke, S.; Triche, T.; Bootwalla, M. methylumi: Handle Illumina Methylation Data. R Package Version. 2015. Available online: https://www.bioconductor.org/packages/release/bioc/html/methylumi.html (accessed on 9 April 2021).
  27. Chen, Y.A.; Lemire, M.; Choufani, S.; Butcher, D.T.; Grafodatskaya, D.; Zanke, B.W.; Gallinger, S.; Hudson, T.J.; Weksberg, R. Discovery of cross-reactive probes and polymorphic CpGs in the Illumina Infinium HumanMethylation450 microarray. Epigenetics 2013, 8, 203–209. [Google Scholar] [CrossRef] [Green Version]
  28. Du, P.; Kibbe, W.A.; Lin, S.M. lumi: A pipeline for processing Illumina microarray. Bioinformatics 2008, 24, 1547–1548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Pidsley, R.; Y Wong, C.C.; Volta, M.; Lunnon, K.; Mill, J.; Schalkwyk, L.C. A data-driven approach to preprocessing Illumina 450K methylation array data. BMC Genom. 2013, 14, 293. [Google Scholar] [CrossRef] [Green Version]
  30. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  31. Du, P.; Bourgon, R. methyAnalysis: DNA Methylation Data Analysis and Visualization. R Package Version. 2014, pp. 1–33. Available online: http://bioconductor.org/packages/release/bioc/html/methyAnalysis.html (accessed on 9 April 2021).
  32. Aryee, M.J.; Jaffe, A.E.; Corrada-Bravo, H.; Ladd-Acosta, C.; Feinberg, A.P.; Hansen, K.D.; Irizarry, R.A. Minfi: A flexible and comprehensive Bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 2014, 30, 1363–1369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Liao, Y.; Wang, J.; Jaehnig, E.J.; Shi, Z.; Zhang, B. WebGestalt 2019: Gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res. 2019, 47, W199–W205. [Google Scholar] [CrossRef] [Green Version]
  34. Kutmon, M.; van Iersel, M.P.; Bohler, A.; Kelder, T.; Nunes, N.; Pico, A.R.; Evelo, C.T. PathVisio 3: An extendable pathway analysis toolbox. PLoS Comput. Biol. 2015, 11, e1004085. [Google Scholar] [CrossRef] [Green Version]
  35. Chakravarthy, A.; Furness, A.; Joshi, K.; Ghorani, E.; Ford, K.; Ward, M.J.; King, E.V.; Lechner, M.; Marafioti, T.; Quezada, S.A.; et al. Pan-cancer deconvolution of tumour composition using DNA methylation. Nat. Commun. 2018, 9, 3220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Newman, A.M.; Steen, C.B.; Liu, C.L.; Gentles, A.J.; Chaudhuri, A.A.; Scherer, F.; Khodadoust, M.S.; Esfahani, M.S.; Luca, B.A.; Steiner, D.; et al. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat. Biotechnol. 2019, 37, 773–782. [Google Scholar] [CrossRef]
  37. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [Green Version]
  38. Cerami, E.; Gao, J.; Dogrusoz, U.; Gross, B.E.; Sumer, S.O.; Aksoy, B.A.; Jacobsen, A.; Byrne, C.J.; Heuer, M.L.; Larsson, E.; et al. The cBio cancer genomics portal: An open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012, 2, 401–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Díez-Villanueva, A.; Mallona, I.; Peinado, M.A. Wanderer, an interactive viewer to explore DNA methylation and gene expression data in human cancer. Epigenetics Chromatin 2015, 8, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Nusse, R. The WNT Homepage. Available online: https://web.stanford.edu/group/nusselab/cgi-bin/wnt/target_genes (accessed on 24 July 2020).
  41. Poage, G.M.; Houseman, E.A.; Christensen, B.C.; Butler, R.A.; Avissar-Whiting, M.; McClean, M.D.; Waterboer, T.; Pawlita, M.; Marsit, C.J.; Kelsey, K.T. Global hypomethylation identifies Loci targeted for hypermethylation in head and neck cancer. Clin. Cancer Res. 2011, 17, 3579–3589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Lleras, R.A.; Smith, R.V.; Adrien, L.R.; Schlecht, N.F.; Burk, R.D.; Harris, T.M.; Childs, G.; Prystowsky, M.B.; Belbin, T.J. Unique DNA methylation loci distinguish anatomic site and HPV status in head and neck squamous cell carcinoma. Clin. Cancer Res. 2013, 19, 5444–5455. [Google Scholar] [CrossRef] [Green Version]
  43. Campbell, J.D.; Yau, C.; Bowlby, R.; Liu, Y.; Brennan, K.; Fan, H.; Taylor, A.M.; Wang, C.; Walter, V.; Akbani, R.; et al. Genomic, Pathway Network, and Immunologic Features Distinguishing Squamous Carcinomas. Cell Rep. 2018, 23, 194–212.e196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Talukdar, F.R.; Soares Lima, S.C.; Khoueiry, R.; Laskar, R.S.; Cuenin, C.; Sorroche, B.P.; Boisson, A.C.; Abedi-Ardekani, B.; Carreira, C.; Menya, D.; et al. Genome-wide DNA methylation profiling of esophageal squamous cell carcinoma from global high-incidence regions identifies crucial genes and potential cancer markers. Cancer Res. 2021, 81, 2612–2624. [Google Scholar] [CrossRef]
  45. Pyeon, D.; Newton, M.A.; Lambert, P.F.; den Boon, J.A.; Sengupta, S.; Marsit, C.J.; Woodworth, C.D.; Connor, J.P.; Haugen, T.H.; Smith, E.M.; et al. Fundamental differences in cell cycle deregulation in human papillomavirus-positive and human papillomavirus-negative head/neck and cervical cancers. Cancer Res. 2007, 67, 4605–4619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Guo, T.; Sakai, A.; Afsari, B.; Considine, M.; Danilova, L.; Favorov, A.V.; Yegnasubramanian, S.; Kelley, D.Z.; Flam, E.; Ha, P.K.; et al. A Novel Functional Splice Variant of. Cancer Res. 2017, 77, 5248–5258. [Google Scholar] [CrossRef] [Green Version]
  47. Guo, T.; Gaykalova, D.A.; Considine, M.; Wheelan, S.; Pallavajjala, A.; Bishop, J.A.; Westra, W.H.; Ideker, T.; Koch, W.M.; Khan, Z.; et al. Characterization of functionally active gene fusions in human papillomavirus related oropharyngeal squamous cell carcinoma. Int. J. Cancer 2016, 139, 373–382. [Google Scholar] [CrossRef] [Green Version]
  48. Ren, S.; Gaykalova, D.; Wang, J.; Guo, T.; Danilova, L.; Favorov, A.; Fertig, E.; Bishop, J.; Khan, Z.; Flam, E.; et al. Discovery and development of differentially methylated regions in human papillomavirus-related oropharyngeal squamous cell carcinoma. Int. J. Cancer 2018, 143, 2425–2436. [Google Scholar] [CrossRef] [Green Version]
  49. Elhalawani, H.; Mohamed, A.S.R.; Elgohari, B.; Lin, T.A.; Sikora, A.G.; Lai, S.Y.; Abusaif, A.; Phan, J.; Morrison, W.H.; Gunn, G.B.; et al. Tobacco exposure as a major modifier of oncologic outcomes in human papillomavirus (HPV) associated oropharyngeal squamous cell carcinoma. BMC Cancer 2020, 20, 912. [Google Scholar] [CrossRef]
  50. Baliga, S.; Klamer, B.; Jhawar, S.; Gamez, M.; Mitchell, D.; Blakaj, A.; Grecula, J.; Gardner, U.; Dibs, K.; Old, M.; et al. Identification of Clinical and Socioeconomic Predictors of Adjuvant Therapy after Trans-Oral Robotic Surgery in Patients with Oropharyngeal Squamous Cell Carcinoma. Cancers 2020, 12, 2474. [Google Scholar] [CrossRef] [PubMed]
  51. Buexm, L.A.; Soares-Lima, S.C.; Brennan, P.; Fernandes, P.V.; de Souza Almeida Lopes, M.; Nascimento de Carvalho, F.; Santos, I.C.; Dias, L.F.; de Queiroz Chaves Lourenço, S.; Ribeiro Pinto, L.F. Hpv impact on oropharyngeal cancer patients treated at the largest cancer center from Brazil. Cancer Lett. 2020, 477, 70–75. [Google Scholar] [CrossRef] [PubMed]
  52. Vossen, D.M.; Verhagen, C.V.M.; Verheij, M.; Wessels, L.F.A.; Vens, C.; van den Brekel, M.W.M. Comparative genomic analysis of oral versus laryngeal and pharyngeal cancer. Oral Oncol. 2018, 81, 35–44. [Google Scholar] [CrossRef] [PubMed]
  53. Brennan, K.; Koenig, J.L.; Gentles, A.J.; Sunwoo, J.B.; Gevaert, O. Identification of an atypical etiological head and neck squamous carcinoma subtype featuring the CpG island methylator phenotype. EBioMedicine 2017, 17, 223–236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Giraldi, L.; Leoncini, E.; Pastorino, R.; Wünsch-Filho, V.; de Carvalho, M.; Lopez, R.; Cadoni, G.; Arzani, D.; Petrelli, L.; Matsuo, K.; et al. Alcohol and cigarette consumption predict mortality in patients with head and neck cancer: A pooled analysis within the International Head and Neck Cancer Epidemiology (INHANCE) Consortium. Ann. Oncol. 2017, 28, 2843–2851. [Google Scholar] [CrossRef]
  55. Mochizuki, D.; Misawa, Y.; Kawasaki, H.; Imai, A.; Endo, S.; Mima, M.; Yamada, S.; Nakagawa, T.; Kanazawa, T.; Misawa, K. Aberrant Epigenetic Regulation in Head and Neck Cancer Due to Distinct EZH2 Overexpression and DNA Hypermethylation. Int. J. Mol. Sci. 2018, 19, 3707. [Google Scholar] [CrossRef] [Green Version]
  56. Hussein, Y.R.; Sood, A.K.; Bandyopadhyay, S.; Albashiti, B.; Semaan, A.; Nahleh, Z.; Roh, J.; Han, H.D.; Lopez-Berestein, G.; Ali-Fehmi, R. Clinical and biological relevance of enhancer of zeste homolog 2 in triple-negative breast cancer. Hum. Pathol. 2012, 43, 1638–1644. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Hoffmann, M.J.; Engers, R.; Florl, A.R.; Otte, A.P.; Muller, M.; Schulz, W.A. Expression changes in EZH2, but not in BMI-1, SIRT1, DNMT1 or DNMT3B are associated with DNA methylation changes in prostate cancer. Cancer Biol. Ther. 2007, 6, 1403–1412. [Google Scholar] [CrossRef] [Green Version]
  58. Karlsson, J.; Valind, A.; Jansson, C.; O’Sullivan, M.J.; Holmquist Mengelbier, L.; Gisselsson, D. Aberrant epigenetic regulation in clear cell sarcoma of the kidney featuring distinct DNA hypermethylation and EZH2 overexpression. Oncotarget 2016, 7, 11127–11136. [Google Scholar] [CrossRef] [PubMed]
  59. Poirier, J.T.; Gardner, E.E.; Connis, N.; Moreira, A.L.; de Stanchina, E.; Hann, C.L.; Rudin, C.M. DNA methylation in small cell lung cancer defines distinct disease subtypes and correlates with high expression of EZH2. Oncogene 2015, 34, 5869–5878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Cao, W.; Younis, R.H.; Li, J.; Chen, H.; Xia, R.; Mao, L.; Chen, W.; Ren, H. EZH2 promotes malignant phenotypes and is a predictor of oral cancer development in patients with oral leukoplakia. Cancer Prev. Res. 2011, 4, 1816–1824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Lin, H.Y.; Hung, S.K.; Lee, M.S.; Chiou, W.Y.; Huang, T.T.; Tseng, C.E.; Shih, L.Y.; Lin, R.I.; Lin, J.M.; Lai, Y.H.; et al. DNA methylome analysis identifies epigenetic silencing of FHIT as a determining factor for radiosensitivity in oral cancer: An outcome-predicting and treatment-implicating study. Oncotarget 2015, 6, 915–934. [Google Scholar] [CrossRef] [Green Version]
  62. Chen, X.; Liu, L.; Mims, J.; Punska, E.C.; Williams, K.E.; Zhao, W.; Arcaro, K.F.; Tsang, A.W.; Zhou, X.; Furdui, C.M. Analysis of DNA methylation and gene expression in radiation-resistant head and neck tumors. Epigenetics 2015, 10, 545–561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Koumangoye, R.B.; Andl, T.; Taubenslag, K.J.; Zilberman, S.T.; Taylor, C.J.; Loomans, H.A.; Andl, C.D. SOX4 interacts with EZH2 and HDAC3 to suppress microRNA-31 in invasive esophageal cancer cells. Mol. Cancer 2015, 14, 24. [Google Scholar] [CrossRef] [Green Version]
  64. Liu, F.; Gu, L.; Cao, Y.; Fan, X.; Zhang, F.; Sang, M. Aberrant overexpression of EZH2 and H3K27me3 serves as poor prognostic biomarker for esophageal squamous cell carcinoma patients. Biomarkers 2016, 21, 80–90. [Google Scholar] [CrossRef] [PubMed]
  65. Ghosh, J.; Schultz, B.; Coutifaris, C.; Sapienza, C. Highly variant DNA methylation in normal tissues identifies a distinct subclass of cancer patients. Adv. Cancer Res. 2019, 142, 1–22. [Google Scholar] [CrossRef]
  66. Kim, S.Y.; Han, Y.K.; Song, J.M.; Lee, C.H.; Kang, K.; Yi, J.M.; Park, H.R. Aberrantly hypermethylated tumor suppressor genes were identified in oral squamous cell carcinoma (OSCC). Clin. Epigenetics 2019, 11, 116. [Google Scholar] [CrossRef] [Green Version]
  67. Ma, J.; Li, T.F.; Han, X.W.; Yuan, H.F. Downregulated MEG3 contributes to tumour progression and poor prognosis in oesophagal squamous cell carcinoma by interacting with miR-4261, downregulating DKK2 and activating the Wnt/β-catenin signalling. Artif. Cells Nanomed. Biotechnol. 2019, 47, 1513–1523. [Google Scholar] [CrossRef]
  68. Singh, V.; Singh, A.P.; Sharma, I.; Singh, L.C.; Sharma, J.; Borthakar, B.B.; Rai, A.K.; Kataki, A.C.; Kapur, S.; Saxena, S. Epigenetic deregulations of Wnt/β-catenin and transforming growth factor beta-Smad pathways in esophageal cancer: Outcome of DNA methylation. J. Cancer Res. Ther. 2019, 15, 192–203. [Google Scholar] [CrossRef]
  69. Paluszczak, J.; Kiwerska, K.; Mielcarek-Kuchta, D. Frequent methylation of DAB2, a Wnt pathway antagonist, in oral and oropharyngeal squamous cell carcinomas. Pathol. Res. Pract. 2018, 214, 314–317. [Google Scholar] [CrossRef]
  70. Paluszczak, J. The Significance of the Dysregulation of Canonical Wnt Signaling in Head and Neck Squamous Cell Carcinomas. Cells 2020, 9, 723. [Google Scholar] [CrossRef] [Green Version]
  71. Papagerakis, P.; Pannone, G.; Shabana, A.H.; Depondt, J.; Santoro, A.; Ghirtis, K.; Berdal, A.; Papagerakis, S. Aberrant beta-catenin and LEF1 expression may predict the clinical outcome for patients with oropharyngeal cancer. Int. J. Immunopathol. Pharmacol. 2012, 25, 135–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Reyes, M.; Peña-Oyarzun, D.; Maturana, A.; Torres, V.A. Nuclear localization of β-catenin and expression of target genes are associated with increased Wnt secretion in oral dysplasia. Oral Oncol. 2019, 94, 58–67. [Google Scholar] [CrossRef]
  73. Silva, B.S.; Castro, C.A.; Von Zeidler, S.L.; Sousa, S.C.; Batista, A.C.; Yamamoto-Silva, F.P. Altered β-catenin expression in oral mucosal dysplasia: A comparative study. J. Appl. Oral Sci. 2015, 23, 472–478. [Google Scholar] [CrossRef] [PubMed]
  74. Ishida, K.; Ito, S.; Wada, N.; Deguchi, H.; Hata, T.; Hosoda, M.; Nohno, T. Nuclear localization of beta-catenin involved in precancerous change in oral leukoplakia. Mol. Cancer 2007, 6, 62. [Google Scholar] [CrossRef] [Green Version]
  75. Kaur, J.; Sawhney, M.; DattaGupta, S.; Shukla, N.K.; Srivastava, A.; Walfish, P.G.; Ralhan, R. Clinical significance of altered expression of β-catenin and E-cadherin in oral dysplasia and cancer: Potential link with ALCAM expression. PLoS ONE 2013, 8, e67361. [Google Scholar] [CrossRef] [PubMed]
  76. Rampias, T.; Boutati, E.; Pectasides, E.; Sasaki, C.; Kountourakis, P.; Weinberger, P.; Psyrri, A. Activation of Wnt signaling pathway by human papillomavirus E6 and E7 oncogenes in HPV16-positive oropharyngeal squamous carcinoma cells. Mol. Cancer Res. 2010, 8, 433–443. [Google Scholar] [CrossRef] [Green Version]
  77. Al Moustafa, A.E.; Foulkes, W.D.; Benlimame, N.; Wong, A.; Yen, L.; Bergeron, J.; Batist, G.; Alpert, L.; Alaoui-Jamali, M.A. E6/E7 proteins of HPV type 16 and ErbB-2 cooperate to induce neoplastic transformation of primary normal oral epithelial cells. Oncogene 2004, 23, 350–358. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Osei-Sarfo, K.; Tang, X.H.; Urvalek, A.M.; Scognamiglio, T.; Gudas, L.J. The molecular features of tongue epithelium treated with the carcinogen 4-nitroquinoline-1-oxide and alcohol as a model for HNSCC. Carcinogenesis 2013, 34, 2673–2681. [Google Scholar] [CrossRef] [PubMed]
  79. Kanojia, D.; Vaidya, M.M. 4-nitroquinoline-1-oxide induced experimental oral carcinogenesis. Oral Oncol. 2006, 42, 655–667. [Google Scholar] [CrossRef] [PubMed]
  80. Pannone, G.; Bufo, P.; Santoro, A.; Franco, R.; Aquino, G.; Longo, F.; Botti, G.; Serpico, R.; Cafarelli, B.; Abbruzzese, A.; et al. WNT pathway in oral cancer: Epigenetic inactivation of WNT-inhibitors. Oncol. Rep. 2010, 24, 1035–1041. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Xiao, C.; Wang, L.; Zhu, L.; Zhang, C.; Zhou, J. Secreted frizzled-related protein 2 is epigenetically silenced and functions as a tumor suppressor in oral squamous cell carcinoma. Mol. Med. Rep. 2014, 10, 2293–2298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Paluszczak, J.; Sarbak, J.; Kostrzewska-Poczekaj, M.; Kiwerska, K.; Jarmuż-Szymczak, M.; Grenman, R.; Mielcarek-Kuchta, D.; Baer-Dubowska, W. The negative regulators of Wnt pathway-DACH1, DKK1, and WIF1 are methylated in oral and oropharyngeal cancer and WIF1 methylation predicts shorter survival. Tumour Biol. 2015, 36, 2855–2861. [Google Scholar] [CrossRef] [Green Version]
  83. Paluszczak, J.; Wiśniewska, D.; Kostrzewska-Poczekaj, M.; Kiwerska, K.; Grénman, R.; Mielcarek-Kuchta, D.; Jarmuż-Szymczak, M. Prognostic significance of the methylation of Wnt pathway antagonists-CXXC4, DACT2, and the inhibitors of sonic hedgehog signaling-ZIC1, ZIC4, and HHIP in head and neck squamous cell carcinomas. Clin. Oral Investig. 2017, 21, 1777–1788. [Google Scholar] [CrossRef] [Green Version]
  84. Zhang, M.; Linghu, E.; Zhan, Q.; He, T.; Cao, B.; Brock, M.V.; Herman, J.G.; Xiang, R.; Guo, M. Methylation of DACT2 accelerates esophageal cancer development by activating Wnt signaling. Oncotarget 2016, 7, 17957–17969. [Google Scholar] [CrossRef]
  85. Cao, B.; Yang, W.; Jin, Y.; Zhang, M.; He, T.; Zhan, Q.; Herman, J.G.; Zhong, G.; Guo, M. Silencing NKD2 by Promoter Region Hypermethylation Promotes Esophageal Cancer Progression by Activating Wnt Signaling. J. Thorac. Oncol. 2016, 11, 1912–1926. [Google Scholar] [CrossRef] [Green Version]
  86. Guo, Q.; Wang, H.B.; Li, Y.H.; Li, H.F.; Li, T.T.; Zhang, W.X.; Xiang, S.S.; Sun, Z.Q. Correlations of Promoter Methylation in WIF-1, RASSF1A, and CDH13 Genes with the Risk and Prognosis of Esophageal Cancer. Med. Sci. Monit. 2016, 22, 2816–2824. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  87. Meng, Y.; Wang, Q.G.; Wang, J.X.; Zhu, S.T.; Jiao, Y.; Li, P.; Zhang, S.T. Epigenetic inactivation of the SFRP1 gene in esophageal squamous cell carcinoma. Dig. Dis. Sci. 2011, 56, 3195–3203. [Google Scholar] [CrossRef] [PubMed]
  88. Kishino, T.; Niwa, T.; Yamashita, S.; Takahashi, T.; Nakazato, H.; Nakajima, T.; Igaki, H.; Tachimori, Y.; Suzuki, Y.; Ushijima, T. Integrated analysis of DNA methylation and mutations in esophageal squamous cell carcinoma. Mol. Carcinog. 2016, 55, 2077–2088. [Google Scholar] [CrossRef]
  89. Houseman, E.A.; Accomando, W.P.; Koestler, D.C.; Christensen, B.C.; Marsit, C.J.; Nelson, H.H.; Wiencke, J.K.; Kelsey, K.T. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012, 13, 86. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Lutsik, P.; Slawski, M.; Gasparoni, G.; Vedeneev, N.; Hein, M.; Walter, J. MeDeCom: Discovery and quantification of latent components of heterogeneous methylomes. Genome Biol. 2017, 18, 55. [Google Scholar] [CrossRef] [Green Version]
  91. Zhang, X.; Shi, M.; Chen, T.; Zhang, B. Characterization of the Immune Cell Infiltration Landscape in Head and Neck Squamous Cell Carcinoma to Aid Immunotherapy. Mol. Ther. Nucleic Acids 2020, 22, 298–309. [Google Scholar] [CrossRef]
  92. Wood, O.; Clarke, J.; Woo, J.; Mirza, A.H.; Woelk, C.H.; Thomas, G.J.; Vijayanand, P.; King, E.; Ottensmeier, C.H. Head and Neck Squamous Cell Carcinomas Are Characterized by a Stable Immune Signature Within the Primary Tumor Over Time and Space. Clin. Cancer Res. 2017, 23, 7641–7649. [Google Scholar] [CrossRef] [Green Version]
  93. Barros, L.R.C.; Souza-Santos, P.T.; Pretti, M.A.M.; Vieira, G.F.; Bragatte, M.A.S.; Mendes, M.F.A.; De Freitas, M.V.; Scherer, N.M.; De Oliveira, I.M.; Rapozo, D.C.M.; et al. High infiltration of B cells in tertiary lymphoid structures, TCR oligoclonality, and neoantigens are part of esophageal squamous cell carcinoma microenvironment. J. Leukoc. Biol. 2020, 108, 1307–1318. [Google Scholar] [CrossRef]
  94. Seijas-Tamayo, R.; Fernández-Mateos, J.; Adansa Klain, J.C.; Mesía, R.; Pastor Borgoñón, M.; Pérez-Ruiz, E.; Vázquez Fernández, S.; Salvador Coloma, C.; Rueda Domínguez, A.; Taberna, M.; et al. Epidemiological characteristics of a Spanish cohort of patients diagnosed with squamous cell carcinoma of head and neck: Distribution of risk factors by tumor location. Clin. Transl. Oncol. 2016, 18, 1114–1122. [Google Scholar] [CrossRef] [PubMed]
  95. Arnal, M.J.D.; Arenas, Á.F.; Arbeloa, Á.L. Esophageal cancer: Risk factors, screening and endoscopic treatment in Western and Eastern countries. World J. Gastroenterol. 2015, 21, 7933. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Overall DNA methylation profile of head and neck and esophageal squamous cell carcinomas. (A) Principal component analysis of the 1000 most variable CpG sites in non-tumor tissues of esophagus, larynx, oral cavity and oropharynx. (B) Principal component analysis of the 1000 most variable CpG sites in squamous cell carcinomas of the esophagus, larynx, oral cavity and oropharynx. (C) Bar graphs showing the proportion of hypermethylated and hypomethylated probes, according to their genic location, for each tumor type relative to their non-tumor counterparts. In red: hypermethylated probes; in blue: hypomethylated probes. On the bottom, a table with the proportion of 450 K beadchip probes in each gene region as well as the proportion of DMPs in each gene region by tumor type is shown. p-values of chi-squared tests calculated based on the observed and expected absolute numbers of probes are also shown. (D) Bar graphs showing the proportion of hypermethylated and hypomethylated probes, according to their relation to CpG islands, for each tumor type relative to their non-tumor counterparts. In red: hypermethylated probes; in blue: hypomethylated probes. On the bottom, a table with the proportion of 450 K beadchip probes according to their relation to CpG island as well as the proportion of DMPs according to their relation to CpG island by tumor type is shown. p-values of chi-squared tests calculated based on the observed and expected absolute numbers of probes are also shown. (E) Circos plot showing the delta-betas (differences between the mean methylation in tumors and non-tumor tissues) for each probe in the beadchip according to tumor type. From the outer circle to the inner circle: esophageal squamous cell carcinoma (1), laryngeal squamous cell carcinoma (2), oral cavity squamous cell carcinoma (3) and oropharyngeal squamous cell carcinoma (4). The chromosomes are included in the center of the circus plot and genes carrying the 20 DMPs with the largest variances across tumor subsites are shown. Regions in red: hypermethylated sites in tumors relative to their non-tumor counterparts; regions in blue: hypomethylated sites in tumors relative to their non-tumor counterparts, regions in white: probes removed from the analysis for quality reasons or filtering purposes. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; MDS, multidimensional scaling; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Figure 1. Overall DNA methylation profile of head and neck and esophageal squamous cell carcinomas. (A) Principal component analysis of the 1000 most variable CpG sites in non-tumor tissues of esophagus, larynx, oral cavity and oropharynx. (B) Principal component analysis of the 1000 most variable CpG sites in squamous cell carcinomas of the esophagus, larynx, oral cavity and oropharynx. (C) Bar graphs showing the proportion of hypermethylated and hypomethylated probes, according to their genic location, for each tumor type relative to their non-tumor counterparts. In red: hypermethylated probes; in blue: hypomethylated probes. On the bottom, a table with the proportion of 450 K beadchip probes in each gene region as well as the proportion of DMPs in each gene region by tumor type is shown. p-values of chi-squared tests calculated based on the observed and expected absolute numbers of probes are also shown. (D) Bar graphs showing the proportion of hypermethylated and hypomethylated probes, according to their relation to CpG islands, for each tumor type relative to their non-tumor counterparts. In red: hypermethylated probes; in blue: hypomethylated probes. On the bottom, a table with the proportion of 450 K beadchip probes according to their relation to CpG island as well as the proportion of DMPs according to their relation to CpG island by tumor type is shown. p-values of chi-squared tests calculated based on the observed and expected absolute numbers of probes are also shown. (E) Circos plot showing the delta-betas (differences between the mean methylation in tumors and non-tumor tissues) for each probe in the beadchip according to tumor type. From the outer circle to the inner circle: esophageal squamous cell carcinoma (1), laryngeal squamous cell carcinoma (2), oral cavity squamous cell carcinoma (3) and oropharyngeal squamous cell carcinoma (4). The chromosomes are included in the center of the circus plot and genes carrying the 20 DMPs with the largest variances across tumor subsites are shown. Regions in red: hypermethylated sites in tumors relative to their non-tumor counterparts; regions in blue: hypomethylated sites in tumors relative to their non-tumor counterparts, regions in white: probes removed from the analysis for quality reasons or filtering purposes. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; MDS, multidimensional scaling; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Cancers 13 03014 g001
Figure 2. Signaling pathway enrichment analysis with genes affected by differentially methylated regions. DMRs were identified by comparing each tumor type with its corresponding non-tumor tissue. In all comparisons, only DMRs containing at least 5 differentially methylated CpG sites, with FDR < 0.05 were considered. Genes with at least one DMR were included in the Over-representation analysis (ORA) using KEGG database to identify enriched signaling pathways. In the highlighted rectangle in the Figure, enriched pathways for each tumor type are shown (FDR < 0.05). The overlapping of enriched pathways (color filled Venn diagram) showed that 4 signaling pathways were over-represented in all tumor types, cAMP Signaling Pathway, Calcium Signaling Pathway, Glutamatergic Synapse, and Neuroactive Ligand-Receptor Interaction. However, the overlapping of DMR-carrying genes showed that (non-colored filled Venn diagrams) the genes affected varied in the different tumor types. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Figure 2. Signaling pathway enrichment analysis with genes affected by differentially methylated regions. DMRs were identified by comparing each tumor type with its corresponding non-tumor tissue. In all comparisons, only DMRs containing at least 5 differentially methylated CpG sites, with FDR < 0.05 were considered. Genes with at least one DMR were included in the Over-representation analysis (ORA) using KEGG database to identify enriched signaling pathways. In the highlighted rectangle in the Figure, enriched pathways for each tumor type are shown (FDR < 0.05). The overlapping of enriched pathways (color filled Venn diagram) showed that 4 signaling pathways were over-represented in all tumor types, cAMP Signaling Pathway, Calcium Signaling Pathway, Glutamatergic Synapse, and Neuroactive Ligand-Receptor Interaction. However, the overlapping of DMR-carrying genes showed that (non-colored filled Venn diagrams) the genes affected varied in the different tumor types. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Cancers 13 03014 g002
Figure 3. The molecular mechanisms affecting the WNT signaling pathway are different in ESCC, LSCC, OSCC and OPSCC. (A) Bar graph showing the percentage of genes of the WNT signaling pathway carrying at least one DMR in each tumor type within our dataset. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the frequency of altered and non-altered genes, respectively. (B) Bar graphs showing the correlation between the methylation levels of each affected WNT DMR and the mRNA expression of its associated gene in each tumor type. For this analysis, TCGA datasets were used. Spearman’s correlation rho is shown in X-axis. Genes were grouped according to their function in WNT signaling pathway and DMRs were grouped according to their location, within promoters or outside promoters. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the significant (BH adjusted p-value < 0.05) and not significant correlations, respectively. (C) Bar graph showing the percentage of cases carrying driver mutations (DM) and/or deep deletions in AJUBA, FAT1 and NOTCH1 in each tumor type within TCGA dataset. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the frequency of altered and non-altered genes, respectively. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; NTST, non-tumor surrounding tissue; NTT, non-tumor tissue; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Figure 3. The molecular mechanisms affecting the WNT signaling pathway are different in ESCC, LSCC, OSCC and OPSCC. (A) Bar graph showing the percentage of genes of the WNT signaling pathway carrying at least one DMR in each tumor type within our dataset. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the frequency of altered and non-altered genes, respectively. (B) Bar graphs showing the correlation between the methylation levels of each affected WNT DMR and the mRNA expression of its associated gene in each tumor type. For this analysis, TCGA datasets were used. Spearman’s correlation rho is shown in X-axis. Genes were grouped according to their function in WNT signaling pathway and DMRs were grouped according to their location, within promoters or outside promoters. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the significant (BH adjusted p-value < 0.05) and not significant correlations, respectively. (C) Bar graph showing the percentage of cases carrying driver mutations (DM) and/or deep deletions in AJUBA, FAT1 and NOTCH1 in each tumor type within TCGA dataset. ESCC, blue; LSCC, pink; OSCC, green; OPSCC, yellow. Dark and light colors represent the frequency of altered and non-altered genes, respectively. ESCC, esophageal squamous cell carcinoma; LSCC, laryngeal squamous cell carcinoma; NTST, non-tumor surrounding tissue; NTT, non-tumor tissue; OSCC, oral cavity squamous cell carcinoma; OPSCC, oropharyngeal squamous cell carcinoma.
Cancers 13 03014 g003
Figure 4. WNT pathway disruption by DNA methylation alterations in ESCC, LSCC, OSCC and OPSCC. Schematic representation of the WNT pathway showing the genes involved and their DMR alterations in: (A) esophageal squamous cell carcinoma (ESCC); (B) laryngeal squamous cell carcinoma (LSCC); (C) oral cavity squamous cell carcinoma (OSCC); (D) oropharyngeal squamous cell carcinoma (OPSCC). Each rectangle represents a gene involved in the pathway and is divided into two halves: the left halves represent the DMRs annotated to promoters; and the right halves represent the DMRs annotated outside promoters. The colors indicate whether the region was found to be hypermethylated (red) or hypomethylated (blue) in each tumor relative to their respective non-tumor tissue. When more than one DMR was differentially methylated in a given region, the mean delta beta of all DMRs was represented.
Figure 4. WNT pathway disruption by DNA methylation alterations in ESCC, LSCC, OSCC and OPSCC. Schematic representation of the WNT pathway showing the genes involved and their DMR alterations in: (A) esophageal squamous cell carcinoma (ESCC); (B) laryngeal squamous cell carcinoma (LSCC); (C) oral cavity squamous cell carcinoma (OSCC); (D) oropharyngeal squamous cell carcinoma (OPSCC). Each rectangle represents a gene involved in the pathway and is divided into two halves: the left halves represent the DMRs annotated to promoters; and the right halves represent the DMRs annotated outside promoters. The colors indicate whether the region was found to be hypermethylated (red) or hypomethylated (blue) in each tumor relative to their respective non-tumor tissue. When more than one DMR was differentially methylated in a given region, the mean delta beta of all DMRs was represented.
Cancers 13 03014 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Soares-Lima, S.C.; Mehanna, H.; Camuzi, D.; de Souza-Santos, P.T.; Simão, T.d.A.; Nicolau-Neto, P.; Almeida Lopes, M.d.S.; Cuenin, C.; Talukdar, F.R.; Batis, N.; et al. Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption. Cancers 2021, 13, 3014. https://doi.org/10.3390/cancers13123014

AMA Style

Soares-Lima SC, Mehanna H, Camuzi D, de Souza-Santos PT, Simão TdA, Nicolau-Neto P, Almeida Lopes MdS, Cuenin C, Talukdar FR, Batis N, et al. Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption. Cancers. 2021; 13(12):3014. https://doi.org/10.3390/cancers13123014

Chicago/Turabian Style

Soares-Lima, Sheila Coelho, Hisham Mehanna, Diego Camuzi, Paulo Thiago de Souza-Santos, Tatiana de Almeida Simão, Pedro Nicolau-Neto, Monique de Souza Almeida Lopes, Cyrille Cuenin, Fazlur Rahman Talukdar, Nikolaos Batis, and et al. 2021. "Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption" Cancers 13, no. 12: 3014. https://doi.org/10.3390/cancers13123014

APA Style

Soares-Lima, S. C., Mehanna, H., Camuzi, D., de Souza-Santos, P. T., Simão, T. d. A., Nicolau-Neto, P., Almeida Lopes, M. d. S., Cuenin, C., Talukdar, F. R., Batis, N., Costa, I., Dias, F., Degli Esposti, D., Boroni, M., Herceg, Z., & Ribeiro Pinto, L. F. (2021). Upper Aerodigestive Tract Squamous Cell Carcinomas Show Distinct Overall DNA Methylation Profiles and Different Molecular Mechanisms behind WNT Signaling Disruption. Cancers, 13(12), 3014. https://doi.org/10.3390/cancers13123014

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop