A Minimal DNA Methylation Signature in Oral Tongue Squamous Cell Carcinoma Links Altered Methylation with Tumor Attributes

Universal penalized regression (Elastic-net) model with differentially methylated promoters for oral cancer prediction

BACKGROUND

DNA methylation showed notable potential to act as a diagnostic marker in many cancers. Many studies proposed DNA methylation biomarker in OSCC detection, while most of these studies are limited to specific cohorts or geographical location. However, the generalizability of DNA methylation as a diagnostic marker in oral cancer across different geographical locations is yet to be investigated.

METHODS

We used genome-wide methylation data from 384 oral cavity cancer and normal tissues from TCGA HNSCC and eastern India. The common differentially methylated CpGs in these two cohorts were used to develop an Elastic-net model that can be used for the diagnosis of OSCC. The model was validated using 812 HNSCC and normal samples from different anatomical sites of oral cavity from seven countries. Droplet Digital PCR of methyl-sensitive restriction enzyme digested DNA (ddMSRE) was used for quantification of methylation and validation of the model with 22 OSCC and 22 contralateral normal samples. Additionally, pyrosequencing was used to validate the model using 46 OSCC and 25 adjacent normal and 21 contralateral normal tissue samples.

RESULTS

With ddMSRE, our model showed 91% sensitivity, 100% specificity, and 95% accuracy in classifying OSCC from the contralateral normal tissues. Validation of the model with pyrosequencing also showed 96% sensitivity, 91% specificity, and 93% accuracy for classifying the OSCC from contralateral normal samples, while in case of adjacent normal samples we found similar sensitivity but with 20% specificity, suggesting the presence of early disease methylation signature at the adjacent normal samples. Methylation array data of HNSCC and normal tissues from different geographical locations and different anatomical sites showed comparable sensitivity, specificity, and accuracy in detecting oral cavity cancer with across. Similar results were also observed for different stages of oral cavity cancer.

CONCLUSIONS

Our model identified crucial genomic regions affected by DNA methylation in OSCC and showed similar accuracy in detecting oral cancer across different geographical locations. The high specificity of this model in classifying contralateral normal samples from the oral cancer compared to the adjacent normal samples suggested applicability of the model in early detection.

Regulation of the Lewis Blood Group Antigen Expression: A Literature Review Supplemented with Computational Analysis

Background

The Lewis (Le) blood group system, unlike most other blood groups, is not defined by antigens produced internally to the erythrocytes and their precursors but rather by glycan antigens adsorbed on to the erythrocyte membrane from the plasma. These oligosaccharides are synthesized by the two fucosyltransferases FUT2 and FUT3 mainly in epithelial cells of the digestive tract and transferred to the plasma. At their place of synthesis, some Lewis blood group carbohydrate antigen variants also seem to be involved in various gastrointestinal malignancies. However, relatively little is known about the transcriptional regulation of FUT2 and FUT3.

Summary

To address this question, we screened existing literature and additionally used in silico prediction tools to identify novel candidate regulators for FUT2 and FUT3 and combine these findings with already known data on their regulation. With this approach, we were able to describe a variety of transcription factors, RNA binding proteins and microRNAs, which increase FUT2 and FUT3 transcription and translation upon interaction.

Key Messages

Understanding the regulation of FUT2 and FUT3 is crucial to fully understand the blood group system Lewis (ISBT 007 LE) phenotypes, to shed light on the role of the different Lewis antigens in various pathologies, and to identify potential new diagnostic targets for these diseases.

Genome-wide DNA methylation profiling in tongue squamous cell carcinoma

OBJECTIVES

To provide a comprehensive characterization of DNA methylome of oral tongue squamous cell carcinoma (OTSCC) and identify novel tumor-specific DNA methylation markers for early detection using saliva.

MATERIAL AND METHODS

Genome-wide DNA methylation analysis including six OTSCC matched adjacent non-tumoral tissue and saliva was performed using Infinium MethylationEPIC array. Differentially methylated levels of selected genes in our OTSCC cohort were further validated using OTSCC methylation data from The Cancer Genome Atlas database (TCGA). The methylation levels of a set of tumor-specific hypermethylated genes associated with a downregulated expression were evaluated in saliva. Receiver operating characteristic (ROC) curves were performed to assess the diagnostic value of DNA methylation markers.

RESULTS

A total of 25,890 CpGs (20,505 hypomethylated and 5385 hypermethylated) were differentially methylated (DMCpGs) between OTSCC and adjacent non-tumoral tissue. Hypermethylation of 11 tumor-specific genes was validated in OTSCC TCGA cohort. Of these 11 genes, A2BP1, ANK1, ALDH1A2, GFRA1, TTYH1, and PDE4B were also hypermethylated in saliva. These six salivary methylated genes showed high diagnostic accuracy (≥0.800) for discriminating patients from controls.

CONCLUSIONS

This is the first largest genome-wide DNA methylation study on OTSCC that identifies a group of novel tumor-specific DNA methylation markers with diagnostic potential in saliva.

HOXA1 3'UTR Methylation Is a Potential Prognostic Biomarker in Oral Squamous cell Carcinoma

BACKGROUND

HOXA1 is a prognostic marker and a potential predictive biomarker for radioresistance in head and neck tumors. Its overexpression has been associated with promoter methylation and a worse prognosis in oral squamous cell carcinoma (OSCC) patients. However, opposite outcomes are also described. The effect of the methylation of this gene on different gene regions, other than the promoter, remains uncertain. We investigated the methylation profile at different genomic regions of HOXA1 in OSCC and correlated differentially methylated CpG sites with clinicopathological data.

METHODS

The HOXA1 DNA methylation status was evaluated by analyzing data from The Cancer Genome Atlas and three Gene Expression Omnibus datasets. Significant differentially methylated CpG sites were considered with a |∆β| ≥ 0.10 and a Bonferroni-corrected p-value < 0.01. Differentially methylated CpGs were validated by pyrosequencing using two independent cohorts of 15 and 47 OSCC patients, respectively.

RESULTS

Compared to normal tissues, we found significantly higher DNA methylation levels in the 3'UTR region of HOXA1 in OSCC. Higher methylation levels in tumor samples were positively correlated with smoking habits and patients' overall survival.

CONCLUSIONS

Our findings suggest that HOXA1 gene body methylation is a promising prognostic biomarker for OSCC with potential clinical applications in patient monitoring.

Regulation of Cellular-Signaling Pathways by Mammalian Proteins Containing Bacterial EPIYA or EPIYA-Like Motifs Predicted to be Phosphorylated

The effector proteins of several pathogenic bacteria contain the Glu-Pro-Ile-Tyr-Ala (EPIYA) motif or other similar motifs. The EPIYA motif is delivered into the host cells by type III and IV secretion systems, through which its tyrosine residue undergoes phosphorylation by host kinases. These motifs atypically interact with a wide range of Src homology 2 (SH2) domain-containing mammalian proteins through tyrosine phosphorylation, which leads to the perturbation of multiple signaling cascades, the spread of infection, and improved bacterial colonization. Interestingly, it has been reported that EPIYA (or EPIYA-like) motifs exist in mammalian proteomes and regulate mammalian cellular-signaling pathways, leading to homeostasis and disease pathophysiology. It is possible that pathogenic bacteria have exploited EPIYA (or EPIYA-like) motifs from mammalian proteins and that the mammalian EPIYA (or EPIYA-like) motifs have evolved to have highly specific interactions with SH2 domain-containing proteins. In this review, we focus on the regulation of mammalian cellular-signaling pathways by mammalian proteins containing these motifs.

Investigation of Roles of SLC38A1 in Proliferation and Differentiation of Mouse Tongue Epithelium and Expression in Human Oral Tongue Squamous Cell Carcinoma

The aerobic glycolytic pathway, boosting lactate formation, and glutamine addiction are two hallmarks of cancer pathophysiology. Consistent with this, several cell membrane glutamine transporters, belonging to different solute carrier (SLC) families, have been shown to be upregulated in a cell-specific manner to furnish the cells with glutamine and glutamine-derived metabolic intermediates. Among them, the system A transporter Slc38a1 has a higher affinity for glutamine compared to other SLC transporters, and it undergoes highly multifaceted regulation at gene and protein levels. The current study aimed to investigate the functional role of Slc38a1 in the proliferation and maturation of the mouse tongue epithelium. Secondly, we aimed to examine the expression of SLC38A1 and its regulation in human tongue oral squamous cell carcinoma (OTSCC). Employing Slc38a1 wild-type and knockout mice, we showed that Slc38a1 was not directly linked to the regulation of the proliferation and differentiation of the mouse tongue epithelium. External transcriptomic datasets and Western blot analyses showed upregulation of SLC38A1 mRNA/protein in human OTSCC and oral cancer cell lines as compared to the corresponding controls. Further, an investigation of external datasets indicated that mechanisms other than the amplification of the SLC38A1 chromosomal locus or hypomethylation of the SLC38A1 promoter region might be important for the upregulation of SLC38A1 in OTSCC.

ATHENA: an independently validated autophagy-related epigenetic prognostic prediction model of head and neck squamous cell carcinoma

The majority of these existing prognostic models of head and neck squamous cell carcinoma (HNSCC) have unsatisfactory prediction accuracy since they solely utilize demographic and clinical information. Leveraged by autophagy-related epigenetic biomarkers, we aim to develop a better prognostic prediction model of HNSCC incorporating CpG probes with either main effects or gene-gene interactions. Based on DNA methylation data from three independent cohorts, we applied a 3-D analysis strategy to develop An independently validated auTophagy-related epigenetic prognostic prediction model of HEad and Neck squamous cell carcinomA (ATHENA). Compared to prediction models with only demographic and clinical information, ATHENA has substantially improved discriminative ability, prediction accuracy and more clinical net benefits, and shows robustness in different subpopulations, as well as external populations. Besides, epigenetic score of ATHENA is significantly associated with tumor immune microenvironment, tumor-infiltrating immune cell abundances, immune checkpoints, somatic mutation and immunity-related drugs. Taken together these results, ATHENA has the demonstrated feasibility and utility of predicting HNSCC survival ( http://bigdata.njmu.edu.cn/ATHENA/ ).

Role of epigenetics in OSCC: an understanding above genetics

Oral cavity cancer is categorized under head and neck cancer that frequently develops from squamous cells hence also known as oral squamous cell carcinoma (OSCC). Although molecular markers for oral cavity cancer are already known, epigenetic signatures for the same haven't been explored much. Epigenetic and genetic alterations were initially thought to be discrete mechanisms driving the tumour but the whole exome sequencing of various cancers has revealed the interdependency of epigenetics and genetic alterations. The reversible nature of these epigenetic changes makes them an alluring target for cancer therapeutics. The primary epigenetic alterations in cancer include DNA methylation and histone modifications. These alterations are useful for patient early detection and prognostication. This review summarizes the epigenetic perspective to understand the etiology, epigenetic biomarkers, and epi-drugs for better predictive diagnosis and treatment of OSCC.

Functional characterization of TSPAN7 as a novel indicator for immunotherapy in glioma

Glioma is the most common primary malignant tumor of the central nervous system in clinical practice. Most adult diffuse gliomas have poor efficacy after standard treatment, especially glioblastoma. With the in-depth understanding of brain immune microenvironment, immunotherapy as a new treatment has attracted much attention. In this study, through analyzing a large number of glioma cohorts, we reported that TSPAN7, a member of the tetraspanin family, decreased in high-grade gliomas and low expression was associated with poor prognosis in glioma patients. Meanwhile, the expression pattern of TSPAN7 was verified in glioma clinical samples and glioma cell lines by qPCR, Western Blotting and immunofluorescence. In addition, functional enrichment analysis showed that cell proliferation, EMT, angiogenesis, DNA repair and MAPK signaling pathways were activated in the TSPAN7 lower expression subgroup. Lentiviral plasmids were used to overexpress TSPAN7 in U87 and LN229 glioma cell lines to explore the anti-tumor role of TSPAN7 in glioma. Moreover, by analyzing the relationship between TSPAN7 expression and immune cell infiltration in multiple datasets, we found that TSPAN7 was significantly negatively correlated with the immune infiltration of tumor-related macrophages, especially M2-type macrophages. Further analysis of immune checkpoints showed that, the expression level of TSPAN7 was negatively correlated with the expression of PD-1, PD-L1 and CTLA-4. Using an independent anti-PD-1 immunotherapy cohorts of GBM, we demonstrated that TSPAN7 expression may had a synergistic effect with PD-L1 on the response to immunotherapy. Based on the above findings, we speculate that TSPAN7 can serve as a biomarker for prognosis and a potential immunotherapy target in glioma patients.

DNA Methylation as a Diagnostic, Prognostic, and Predictive Biomarker in Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC) is a term collectively used to describe all cancers that develop in the oral and nasal cavities, the paranasal sinuses, the salivary glands, the pharynx, and the larynx. The majority (75%) of all newly diagnosed cases are observed in patients with locally advanced and aggressive disease, associated with significant relapse rates (30%) and poor prognostic outcomes, despite advances in multimodal treatment. Consequently, there is an unmet need for the identification and application of tools that would enable diagnosis at the earliest possible stage, accurately predict prognostic outcomes, contribute to the timely detection of relapses, and aid in the decision for therapy selection. Recent evidence suggests that DNA methylation can alter the expression of genes in a way that it favors tumorigenesis and tumor progression in HNSCC, and therefore represents a potential source for biomarker identification. This study summarizes the current knowledge on how abnormally methylated DNA profiles in HNSCC patients may contribute to the pathogenesis of HNSCC and designate the methylation patterns that have the potential to constitute clinically valuable biomarkers for achieving significant advances in the management of the disease and for improving survival outcomes in these patients.

Unraveling molecular mechanisms of Head and neck cancer

Malignancies that develop from mucosal epithelium of the upper aerodigestive tract are known as head and neck squamous cell carcinomas (HNSCC). Heterogeneity, late stage diagnosis and high recurrence rate are big hurdles in head and neck treatment regimen. Presently, the biomarkers available for diagnosis and prognosis of HNSCC are based on smoking as the major risk habit. This review shed light on the differential environment of HNSCC in smokeless tobacco consuming Indian patients. Frequent mutation in genes involved in DNA repair pathway (p53), cell proliferation (PIK3CA, HRAS) and cell death (CASP8, FADD) are common in western population. On the contrary, the genes involved in metastasis (MMPs, YAP1), lymphocyte proliferation (TNFRSF4, CD80), cell-cell adhesion (DCC, EDNRB), miRNA processing (DROSHA) and inflammatory responses (TLR9, IL-9) are mutated in Indian HNSCC patients. Gene ontology enrichment analysis highlighted that responses to chemical stimulus, immune pathways and stress pathways are highly enriched in Indian patients.

Diagnostic classification of cancers using DNA methylation of paracancerous tissues

The potential role of DNA methylation from paracancerous tissues in cancer diagnosis has not been explored until now. In this study, we built classification models using well-known machine learning models based on DNA methylation profiles of paracancerous tissues. We evaluated our methods on nine cancer datasets collected from The Cancer Genome Atlas (TCGA) and utilized fivefold cross-validation to assess the performance of models. Additionally, we performed gene ontology (GO) enrichment analysis on the basis of the significant CpG sites selected by feature importance scores of XGBoost model, aiming to identify biological pathways involved in cancer progression. We also exploited the XGBoost algorithm to classify cancer types using DNA methylation profiles of paracancerous tissues in external validation datasets. Comparative experiments suggested that XGBoost achieved better predictive performance than the other four machine learning methods in predicting cancer stage. GO enrichment analysis revealed key pathways involved, highlighting the importance of paracancerous tissues in cancer progression. Furthermore, XGBoost model can accurately classify nine different cancers from TCGA, and the feature sets selected by XGBoost can also effectively predict seven cancer types on independent GEO datasets. This study provided new insights into cancer diagnosis from an epigenetic perspective and may facilitate the development of personalized diagnosis and treatment strategies.

The Expression Profile, Clinical Application and Potential Tumor Suppressing Mechanism of hsa_circ_0001675 in Head and Neck Carcinoma

Purpose

This study sought to identify circular RNAs (circRNA) that participate in the regulation of head and neck cancer (HNC), analyze their clinical application, and predict their molecular mechanism during HNC.

Materials and Methods

High-throughput sequencing was used to analyze circRNA expression in 18 matched HNC and adjacent normal tissues. Target circRNAs with significantly differential expression were obtained. In 103 HNC and adjacent normal tissues, real-time fluorescent quantitative PCR (qRT-PCR) was used to verify the differential expression of target circRNAs. This data was combined with clinicopathological information to analyze the diagnostic value of target circRNA. Bioinformatics was used to find target circRNAs that acted as competitive endogenous RNA (ceRNA) and construct a circRNA-miRNA-mRNA regulatory network. mRNA expression was verified by immunohistochemistry (IHC).

Results

A total of 714 differentially expressed circRNAs were detected in HNC, and the low expression of hsa_circ_0001675 was particularly significant (fold change [FC] = -4.85, P = 6.305E-05). hsa_circ_0001675 had significantly lower expression in HNC than in normal tissue (P < 0.01). Low hsa_circ_0001675 expression was positively associated with tumor invasion and clinical staging (P < 0.05), and its area under the ROC curve (AUC) was 0.7776. Low hsa_circ_0001675 expression also correlated with the overall survival (OS) rate and the progression-free survival (PFS) rate of HNC patients (P < 0.001). Bioinformatics was used to construct a ceRNA network of hsa_circ_0001675 with six differentially expressed miRNAs (hsa-miR-330-5p, hsa-miR-498, hsa-miR-532-3p, hsa-miR-577, hsa-miR-1248, and hsa-miR-1305) and 411 differentially expressed mRNAs and found that the neuroactive ligand-receptor interaction, and the cAMP and calcium signaling pathways were particularly enriched. Further bioinformatics and IHC analysis showed that miR577/TESC is the likely downstream signaling pathway for hsa_circ_0001675.

Conclusion

This study showed that hsa_circ_0001675 is downregulated in HNC and could be an effective biomarker for HNC diagnosis. In addition, hsa_circ_0001675 may have a potential ceRNA mechanism and suppress HNC disease progression through the hsa_circ_0001675-miRNA-mRNA axis.

Role of fibulin-1 gene promoter methylation in the carcinogenesis and development of tongue squamous cell carcinoma

OBJECTIVE

Recently, fibulin-1 (FBLN1) has been shown to be downregulated in various cancers via promoter hypermethylation. Our study aimed to determine the expression and methylation status of FBLN1 in tongue squamous cell carcinoma (TSCC) tissues and cells.

METHODS

Methylation-specific polymerase chain reaction was implemented to detect the methylation status of the FBLN1 gene in TSCC tissues and Western blot analysis was used to detect the expression of FBLN1 protein. The human TSCC cell lines CAL27 and SCC9 were cultured in vitro and treated with 5-aza-deoxycytidine (5-Aza-dC). CCK-8, colony formation, and Transwell assays were performed to test TSCC cell proliferation, migration, and invasion following 5-Aza-dC treatment or overexpression of FBLN1, which was further verified in in vivo experiments.

RESULTS

FBLN1 was hypermethylated and the protein expression was reduced in TSCC tissues. After human TSCC cell lines (CAL27 and SCC9) were treated with 5-Aza-dC or overexpressed FBLN1, FBLN1 expression was upregulated and the TSCC cell proliferation, migration, and invasion abilities were suppressed. In vivo experiments further showed that demethylation or overexpression of FBLN1 slowed tumor growth in nude mice.

CONCLUSION

This study demonstrated that 5-Aza-dC treatment or overexpression of FBLN1 inhibited the growth of human TSCC.

Immune Checkpoint Gene Expression Profiling Identifies Programmed Cell Death Ligand-1 Centered Immunologic Subtypes of Oral and Squamous Cell Carcinoma With Favorable Survival

Objective

This study aimed to identify the programmed death ligand-1 (PDL1, also termed as CD274) and its positively correlated immune checkpoint genes (ICGs) and to determine the immune subtypes of CD274-centered ICG combinations in oral and squamous cell carcinoma (OSCC).

Materials and Methods

Firstly, the 95 ICGs obtained via literature reviews were identified in the Cancer Genome Atlas (TCGA) database in relation to OSCC, and such 88 ICG expression profiles were extracted. ICGs positively correlated with CD274 were utilized for subsequent analysis. The relationship between ICGs positively correlated with CD274 and immunotherapy biomarkers (tumor mutation burden (TMB), and adaptive immune resistance pathway genes) was investigated, and the relationships of these genes with OSCC clinical features were explored. The prognostic values of CD274 and its positively correlated ICGs and also their associated gene pairs were revealed using the survival analysis.

Results

Eight ICGs, including CTLA4, ICOS, TNFRSF4, CD27, B- and T-lymphocyte attenuator (BTLA), ADORA2A, CD40LG, and CD28, were found to be positively correlated with CD274. Among the eight ICGs, seven ICGs (CTLA4, ICOS, TNFRSF4, CD27, BTLA, CD40LG, and CD28) were significantly negatively correlated with TMB. The majority of the adaptive immune resistance pathway genes were positively correlated with ICGs positively correlated with CD274. The survival analysis utilizing the TCGA-OSCC data showed that, although CD274 was not significantly associated with overall survival (OS), the majority of ICGs positively correlated with CD274 (BTLA, CD27, CTLA4, CD40LG, CD28, ICOS, and TNFRSF4) were significantly correlated with OS, whereby their low-expression predicted a favorable prognosis. The survival analysis based on the gene pair subtypes showed that the combination subtypes of CD274_low/BTLA_low, CD274_low/CD27_low, CD274_low/CTLA4_low, CD8A_high/BTLA_low, CD8A_high/CD27_low, and CD8A_high/CTLA4_low predicted favorable OS.

Conclusion

The results in this study provide a theoretical basis for prognostic immune subtyping of OSCC and highlight the importance of developing future immunotherapeutic strategies for treating oral cancer.

Molecular Subtypes of Oral Squamous Cell Carcinoma Based on Immunosuppression Genes Using a Deep Learning Approach

Background: The mechanisms through which immunosuppressed patients bear increased risk and worse survival in oral squamous cell carcinoma (OSCC) are unclear. Here, we used deep learning to investigate the genetic mechanisms underlying immunosuppression in the survival of OSCC patients, especially from the aspect of various survival-related subtypes. Materials and methods: OSCC samples data were obtained from The Cancer Genome Atlas (TCGA), International Cancer Genome Consortium (ICGC), and OSCC-related genetic datasets with survival data in the National Center for Biotechnology Information (NCBI). Immunosuppression genes (ISGs) were obtained from the HisgAtlas and DisGeNET databases. Survival analyses were performed to identify the ISGs with significant prognostic values in OSCC. A deep learning (DL)-based model was established for robustly differentiating the survival subpopulations of OSCC samples. In order to understand the characteristics of the different survival-risk subtypes of OSCC samples, differential expression analysis and functional enrichment analysis were performed. Results: A total of 317 OSCC samples were divided into one inferring cohort (TCGA) and four confirmation cohorts (ICGC set, GSE41613, GSE42743, and GSE75538). Eleven ISGs (i.e., BGLAP, CALCA, CTLA4, CXCL8, FGFR3, HPRT1, IL22, ORMDL3, TLR3, SPHK1, and INHBB) showed prognostic value in OSCC. The DL-based model provided two optimal subgroups of TCGA-OSCC samples with significant differences (p = 4.91E-22) and good model fitness [concordance index (C-index) = 0.77]. The DL model was validated by using four external confirmation cohorts: ICGC cohort (n = 40, C-index = 0.39), GSE41613 dataset (n = 97, C-index = 0.86), GSE42743 dataset (n = 71, C-index = 0.87), and GSE75538 dataset (n = 14, C-index = 0.48). Importantly, subtype Sub1 demonstrated a lower probability of survival and thus a more aggressive nature compared with subtype Sub2. ISGs in subtype Sub1 were enriched in the tumor-infiltrating immune cells-related pathways and cancer progression-related pathways, while those in subtype Sub2 were enriched in the metabolism-related pathways. Conclusion: The two survival subtypes of OSCC identified by deep learning can benefit clinical practitioners to divide immunocompromised patients with oral cancer into two subpopulations and give them target drugs and thus might be helpful for improving the survival of these patients and providing novel therapeutic strategies in the precision medicine area.

Wnt/β-catenin signaling pathway participates in the effect of miR-626 on oral squamous cell carcinoma by targeting RASSF4

BACKGROUND

The role of miR-626 in oral squamous cell carcinoma (OSCC) was investigated by targeting RASSF4.

METHODS

The miR-626 and RASSF4 expression was detected in normal oral mucosa or OSCC tissues and OSCC or normal cells. The methylation status of RASSF4 was analyzed using methylation-specific polymerase chain reaction (PCR). The cytoplasmic/nuclear ratios (C/N ratios) targeted by miR-626 were examined using microarray, followed by a dual-luciferase reporter assay. The subcellular localization of RASSF4 and miR-626 in OSCC cells was determined using RNA fluorescence in situ hybridization (FISH) and immunocytochemistry (ICC), respectively. Ca9-22 and HSC2 cells were divided into mock, inhibitor NC, miR-626 inhibitor, scramble, RASSF4 and miR-626 mimic + RASSF4 groups, and then CCK-8, Annexin V-FITC/PI, wound healing, Transwell, qRT-PCR and western blotting assays were performed.

RESULTS

OSCC tissues and cells had increased miR-626 expression and decreased RASSF4 expression. Patients with RASSF4 methylation had lower RASSF4 expression than those without methylation. In addition, a negative correlation between miR-626 and RASSF4 was found in OSCC tissues, both of which were correlated with the pathological grade, pathological stage, lymph node metastasis and patient prognosis. MiR-626 targeted RASSF4 in OSCC cells. Overexpressed RASSF4 inhibited the proliferation, invasion, migration and epithelial-mesenchymal transition (EMT) of OSCC cells, promoted cell apoptosis, and blocked the Wnt/β-Catenin pathway, which was reversed by miR-626 overexpression.

CONCLUSIONS

Inhibiting miR-626 can regulate the biological characteristics of OSCC cells, including proliferation, invasion, migration, EMT and apoptosis, by targeting RASSF4, which may be related to the Wnt/β-Catenin pathway.

Contribution of Genomics to the Surgical Management and Study of Oral Cancer

Background

Oral squamous cell carcinoma (OSCC) is the most frequent type of tumor arising from the oral cavity. Surgery is the cornerstone of the treatment of these cancers. Tumor biology has long been overlooked as an important contributor to the outcome of surgical procedures, but recent studies are challenging this concept. Molecular analyses of tumor DNA or RNA provide a rich source of information about the biology of OSCC.

Methods

We searched for relevant articles using PubMed. We examined in particular the prospect of applying molecular methods for minimally invasive exploration of OSCC biology.

Results

We examined five potential applications of genomics to the surgical management and study of OSCC: i) assessing oral potentially malignant lesions; ii) tumor staging prior to surgery; iii) predicting postoperative risk in locally advanced tumors; iv) measuring minimal residual disease and optimizing the longitudinal monitoring of OSCC; and v) predicting the efficacy of medical treatment.

Conclusions

Genomic information can be harnessed in order to identify new biomarkers that could improve the staging, choice of therapy and management of OSCC. The identification of new biomarkers is awaited for better personalization of the surgical treatment of OSCC.

TSPAN7 Exerts Anti-Tumor Effects in Bladder Cancer Through the PTEN/PI3K/AKT Pathway

The tetraspanin protein superfamily participate in the dynamic regulation of cellular membrane compartments expressed in a variety of tumor types, which may alter the biological properties of cancer cells such as cell development, activation, growth and motility. The role of tetraspanin 7 (TSPAN7) has never been investigated in bladder cancer (BCa). In this study, we aimed to investigate the biological function of TSPAN7 and its therapeutic potential in human BCa. First, via reverse transcription and quantitative real-time PCR (qRT-PCR), we observed downregulation of TSPAN7 in BCa tissues samples and cell lines and found that this downregulation was associated with a relatively high tumor stage and tumor grade. Low expression of TSPAN7 was significantly correlated with a much poorer prognosis for BCa patients than was high expression. Immunohistochemistry (IHC) showed that low TSPAN7 expression was a high-risk predictor of BCa patient overall survival. Furthermore, the inhibitory effects of TSPAN7 on the proliferation and migration of BCa cell lines were detected by CCK-8, wound-healing, colony formation and transwell assays in vitro. Flow cytometry analysis revealed that TSPAN7 induced BCa cell lines apoptosis and cell cycle arrest. In vivo, tumor growth in nude mice bearing tumor xenografts could be obviously affected by overexpression of TSPAN7. Western blotting showed that overexpression of TSPAN7 activated Bax, cleaved caspase-3 and PTEN but inactivated Bcl-2, p-PI3K, and p-AKT to inhibit BCa cell growth via the PTEN/PI3K/AKT pathway. Taken together, our study will help identify a potential marker for BCa diagnosis and supply a target molecule for BCa treatment.

Sodium hydrogen exchanger 9 NHE9 (SLC9A9) and its emerging roles in neuropsychiatric comorbidity

Variations in SLC9A9 gene expression and protein function are associated with multiple human diseases, which range from Attention-deficit/hyperactivity disorder (ADHD) to glioblastoma multiforme. In an effort to determine the full spectrum of human disease associations with SLC9A9, we performed a systematic review of the literature. We also review SLC9A9's biochemistry, protein structure, and function, as well as its interacting partners with the goal of identifying mechanisms of disease and druggable targets. We report gaps in the literature regarding the genes function along with consistent trends in disease associations that can be used to further research into treating the respective diseases. We report that SLC9A9 has strong associations with neuropsychiatric diseases and various cancers. Interestingly, we find strong overlap in SLC9A9 disease associations and propose a novel role for SLC9A9 in neuropsychiatric comorbidity. In conclusion, SLC9A9 is a multifunctional protein that, through both its endosome regulatory function and its protein-protein interaction network, has the ability to modulate signaling axes, such as the PI3K pathway, among others.

Epigenomic dysregulation-mediated alterations of key biological pathways and tumor immune evasion are hallmarks of gingivo-buccal oral cancer

BACKGROUND

Gingivo-buccal oral squamous cell carcinoma (OSCC-GB) is the most common cancer among men in India and is associated with high mortality. Although OSCC-GB is known to be quite different from tongue cancer in its genomic presentation and its clinical behavior, it is treated identically as tongue cancer. Predictive markers of prognosis and therapy that are specific to OSCC-GB are, therefore, required. Although genomic drivers of OSCC-GB have been identified by whole exome and whole genome sequencing, no epigenome-wide study has been conducted in OSCC-GB; our study has filled this gap, and has discovered and validated epigenomic hallmarks of gingivobuccal oral cancer.

METHODS

We have carried out integrative analysis of epigenomic (n = 87) and transcriptomic (n = 72) profiles of paired tumor-normal tissues collected from OSCC-GB patients from India. Genome-wide DNA methylation assays and RNA-sequencing were performed on high-throughput platforms (Illumina) using a half-sample of randomly selected patients to discover significantly differentially methylated probes (DMPs), which were validated on the remaining half-sample of patients.

RESULTS

About 200 genes showed significant inverse correlation between promoter methylation and expression, of which the most significant genes included genes that act as transcription factors and genes associated with other cancer types. Novel findings of this study include identification of (a) potential immunosuppressive effect in OSCC-GB due to significant promoter hypomethylation driven upregulation of CD274 and CD80, (b) significant dysregulation by epigenetic modification of DNMT3B (upregulation) and TET1 (downregulation); and (c) known drugs that can reverse the direction of dysregulation of gene expression caused by promoter methylation.

CONCLUSIONS

In OSCC-GB patients, there are significant alterations in expression of key genes that (a) regulate normal cell division by maintenance of balanced DNA methylation and transcription process, (b) maintain normal physiological signaling (PPAR, B cell receptor) and metabolism (arachidonic acid) pathways, and (c) provide immune protection against antigens, including tumor cells. These findings indicate novel therapeutic targets, including immunotherapeutic, for treatment of OSCC-GB.

Hypermethylation of DcR1, DcR2, DR4, DR5 gene promoters and clinical significance in tongue carcinoma

OBJECTIVE

Tongue squamous cell carcinoma (TSCC) is one of the most common malignancies in the oral cavity, and its incidence and mortality have been constantly increasing these years. A large number of tumor suppressor genes are involved in the development of the TSCC and it has been reported that the aberrant hypermethylation of tumor suppressor genes may play a key role in the process of the TSCC. In this study, we sought to analyze the association of methylation of DcR1, DcR2, DR4 and DR5 gene promoters and clinical significance in the TSCC to evaluate association between methylation of DcR1, DcR2, DR4 and DR5 gene and Clinical Significance in tongue squamous cell carcinoma.

METHODS

Methylation-specific PCR(MSP) was used to analyze the methylation of the promoters of TRAIL (tumor necrosis factor-related apoptosis-inducing ligand) receptors in 45 TSCC cases. Real-Time PCR was used to detect the expression of the DcR1, DcR2, DR4 and DR5 gene.

RESULTS

All the four genes (DcR1, DcR2, DR4 and DR5) showed different methylation of promoters in TSCC, while methylation of these promoters in paired adjacent normal tissues were almost undetectable. Patients with high methylation index were diagnosed at younger age when compared with the ones with low methylation index. DcR1 and DR4 hypermethylation was correlated significantly with patients' TNM stage.

CONCLUSIONS

Methylation of DcR1, DcR2,DR4 and DR5 promoters are found in TSCC and may associate with its occurrence and development. Taking the reversibility of methylation into account，methylation is a potential targeted therapy of TSCC.

Five extracellular matrix-associated genes upregulated in oral tongue squamous cell carcinoma: An integrated bioinformatics analysis

Despite advancements in treatment regimens, the mortality rate of patients with oral tongue squamous cell carcinoma (OTSCC) is high. In addition, the signaling pathways and oncoproteins involved in OTSCC progression remain largely unknown. Therefore, the aim of the present study was to identify specific prognostic marker for patients at a high risk of developing OTSCC. The present study used four original microarray datasets to identify the key candidate genes involved in OTSCC pathogenesis. Expression profiles of 93 OTSCC tissues and 76 normal tissues from GSE9844, GSE13601, GSE31056 and GSE75538 datasets were investigated. Differentially expressed genes (DEGs) were determined, and gene ontology enrichment and gene interactions were analyzed. The four GSE datasets reported five upregulated and six downregulated DEGs. Five upregulated genes (matrix metalloproteinase 1, 3, 10 and 12 and laminin subunit gamma 2) were localized in the extracellular region of cells and were associated with extracellular matrix disassembly. Furthermore, analysis for The Cancer Genome Atlas database revealed that the aforementioned five upregulated genes were also highly expressed in OTSCC and head and neck squamous cell carcinoma tissues. These results demonstrated that the five upregulated genes may be considered as potential prognostic biomarkers of OTSCC and may serve at understanding OTSCC progression. Upregulated DEGs may therefore represent valuable therapeutic targets to prevent or control OTSCC pathogenesis.

Characterization of a prognostic four‑gene methylation signature associated with radiotherapy for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains one of the most common malignancies associated with poor prognosis. DNA methylation has emerged as an important mechanism underlying the radio-resistance of tumors. Prognostic biomarkers based on radiotherapy-related aberrant DNA methylation are limited. Methylation profiles of 388 patients with HNSCC were acquired from The Cancer Genome Atlas (TCGA) portal. Genes with differentially methylated CpG sites (DMGs) were screened between patients with a favorable and poor prognosis with or without radiotherapy. A weight gene co-methylation network was constructed using a Weighted Gene Co-expression Network Analysis (WGCNA) package. A lasso Cox-PH model was used to identify the optimal panel of genes with the ability to predict survival in these patients. Prognostic performance of the multi-gene methylation signature was assessed in a training set and confirmed in a validation set. A total of 976 DMGs were observed between favorable and poor prognostic samples. Four DMG-enriched co-methylation modules were identified. A four-gene methylation signature was determined by the lasso Cox-PH model that consisted of ZNF10, TMPRSS12, ERGIC2, and RNF215. The risk score based on the four-gene signature was able to divide the training or validation set into two risk groups with significantly different overall survival. Thus, the present study revealed a radiotherapy-related four-gene methylation signature to predict survival outcomes of patients with HNSCC, providing candidate therapeutic targets for novel therapy against HNSCC. However, substantial validation experiments are required.

Detection of High-Risk Human Papillomavirus in Oral Cavity Squamous Cell Carcinoma Using Multiple Analytes and Their Role in Patient Survival

PURPOSE

Accurate detection of human papillomavirus (HPV) in oral cavity squamous cell carcinoma (OSCC) is essential to understanding the role of HPV in disease prognosis and management of patients. We used different analytes and methods to understand the true prevalence of HPV in a cohort of patients with OSCC with different molecular backgrounds, and we correlated HPV data with patient survival.

METHODS

We integrated data from multiple analytes (HPV DNA, HPV RNA, and p16), assays (immunohistochemistry, polymerase chain reaction [PCR], quantitative PCR [qPCR], and digital PCR), and molecular changes (somatic mutations and DNA methylation) from 153 patients with OSCC to correlate p16 expression, HPV DNA, and HPV RNA with HPV incidence and patient survival.

RESULTS

High prevalence (33% to 58%) of HPV16/18 DNA did not correlate with the presence of transcriptionally active viral genomes (15%) in tumors. Eighteen percent of the tumors were p16 positive and only 6% were both HPV DNA and HPV RNA positive. Most tumors with relatively high copy number HPV DNA and/or HPV RNA, but not with HPV DNA alone (irrespective of copy number), were wild-type for TP53 and CASP8 genes. In our study, p16 protein, HPV DNA, and HPV RNA, either alone or in combination, did not correlate with patient survival. Nine HPV-associated genes stratified the virus-positive from the virus-negative tumor group with high confidence ( P < .008) when HPV DNA copy number and/or HPV RNA were considered to define HPV positivity, and not HPV DNA alone, irrespective of copy number ( P < .2).

CONCLUSION

In OSCC, the presence of both HPV RNA and p16 is rare. HPV DNA alone is not an accurate measure of HPV positivity and therefore may not be informative. HPV DNA, HPV RNA, and p16 do not correlate with patients' outcome.

CAFE MOCHA: An Integrated Platform for Discovering Clinically Relevant Molecular Changes in Cancer-An Example of Distant Metastasis- and Recurrence-Linked Classifiers in Head and Neck Squamous Cell Carcinoma

PURPOSE

With large amounts of multidimensional molecular data on cancers generated and deposited into public repositories such as The Cancer Genome Atlas and International Cancer Genome Consortium, a cancer type agnostic and integrative platform will help to identify signatures with clinical relevance. We devised such a platform and showcase it by identifying a molecular signature for patients with metastatic and recurrent (MR) head and neck squamous cell carcinoma (HNSCC).

METHODS

We devised a statistical framework accompanied by a graphical user interface-driven application, Clinical Association of Functionally Established MOlecular CHAnges ( CAFE MOCHA; https://github.com/binaypanda/CAFEMOCHA), to discover molecular signatures linked to a specific clinical attribute in a cancer type. The platform integrates mutations and indels, gene expression, DNA methylation, and copy number variations to discover a classifier first and then to predict an incoming tumor for the same by pulling defined class variables into a single framework that incorporates a coordinate geometry-based algorithm called complete specificity margin-based clustering, which ensures maximum specificity. CAFE MOCHA classifies an incoming tumor sample using either its matched normal or a built-in database of normal tissues. The application is packed and deployed using the install4j multiplatform installer. We tested CAFE MOCHA in HNSCC tumors (n = 513) followed by validation in tumors from an independent cohort (n = 18) for discovering a signature linked to distant MR.

RESULTS

CAFE MOCHA identified an integrated signature, MR44, associated with distant MR HNSCC, with 80% sensitivity and 100% specificity in the discovery stage and 100% sensitivity and 100% specificity in the validation stage.

CONCLUSION

CAFE MOCHA is a cancer type and clinical attribute agnostic statistical framework to discover integrated molecular signatures.

A minimal set of internal control genes for gene expression studies in head and neck squamous cell carcinoma

Selection of the right reference gene(s) is crucial in the analysis and interpretation of gene expression data. The aim of the present study was to discover and validate a minimal set of internal control genes in head and neck tumor studies. We analyzed data from multiple sources (in house whole-genome gene expression microarrays, previously published quantitative real-time PCR (qPCR) data and RNA-seq data from TCGA) to come up with a list of 18 genes (discovery set) that had the lowest variance, a high level of expression across tumors, and their matched normal samples. The genes in the discovery set were ranked using four different algorithms (BestKeeper, geNorm, NormFinder, and comparative delta Ct) and a web-based comparative tool, RefFinder, for their stability and variance in expression across tissues. Finally, we validated their expression using qPCR in an additional set of tumor:matched normal samples that resulted in five genes (RPL30, RPL27, PSMC5, MTCH1, and OAZ1), out of which RPL30 and RPL27 were most stable and were abundantly expressed across the tissues. Our data suggest that RPL30 or RPL27 in combination with either PSMC5 or MTCH1 or OAZ1 can be used as a minimal set of control genes in head and neck tumor gene expression studies.

A 16-gene signature predicting prognosis of patients with oral tongue squamous cell carcinoma

Background

Oral tongue squamous cell carcinoma (OTSCC) is the most common subtype of oral cancer. A predictive gene signature is necessary for prognosis of OTSCC.

Methods

Five microarray data sets of OTSCC from the Gene Expression Omnibus (GEO) and one data set from The Cancer Genome Atlas (TCGA) were obtained. Differentially expressed genes (DEGs) of GEO data sets were identified by integrated analysis. The DEGs associated with prognosis were screened in the TCGA data set by univariate survival analysis to obtain a gene signature. A risk score was calculated as the summation of weighted expression levels with coefficients by Cox analysis. The signature was used to distinguish carcinoma, estimated by receiver operator characteristic curves and the area under the curve (AUC). All were validated in the GEO and TCGA data sets.

Results

Integrated analysis of GEO data sets revealed 300 DEGs. A 16-gene signature and a risk score were developed after survival analysis. The risk score was effective to stratify patients into high-risk and low-risk groups in the TCGA data set (P < 0.001). The 16-gene signature was valid to distinguish the carcinoma from normal samples (AUC 0.872, P < 0.001).

Discussion

We identified a useful 16-gene signature for prognosis of OTSCC patients, which could be applied to clinical practice. Further studies were needed to prove the findings.

Gene expression profiles and protein‑protein interaction networks during tongue carcinogenesis in the tumor microenvironment

Oral tongue squamous cell carcinoma (OTSCC) has a high incidence and is associated with a high mortality rate. Studies regarding the potential molecular mechanism of OTSCC in the tumor microenvironment (TME) are required. The present study aimed to perform bioinformatic analysis to identify important nodes, clusters and functional pathways during tongue carcinogenesis in the TME. After downloading the gene expression data of GSE42780, differentially expressed genes (DEGs) among carcinoma, dysplastic and normal samples in epithelia and fibroblasts were identified using the affy and limma packages with R version 3.3. Subsequently, the Database for Annotation, Visualization and Integrated Discovery was employed to perform Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. Furthermore, a protein-protein interaction (PPI) network was constructed by using the Search Tool for the Retrieval of Interacting Genes/Proteins and analyzed by Cytoscape software. In total, 85 DEGs were identified for tongue epithelia and 46 DEGs were identified for fibroblasts. Neutrophil chemotaxis and inflammatory response from GO, and cytokine-cytokine receptor interaction from KEGG were enriched for epithelia and fibroblasts. The PPI network revealed that C-X-C motif chemokine ligand (Cxcl)1, Cxcl10, Cxcl13, Cxcl2 and pro-platelet basic protein were a key cluster for epithelia, and interleukin (Il)1β, Il1 receptor 2, Il1a and Il1 receptor antagonist were a key cluster for fibroblasts. Therefore, the results indicate that fibroblasts and cytokines associated with an inflammatory immune response contributed substantially to tongue carcinogenesis in the TME, which is useful for the development of OTSCC targeted therapy. However, further investigation is required to elucidate the molecular and cellular mechanisms underlying the inflammatory immune network in the TME.

Seven-CpG-based prognostic signature coupled with gene expression predicts survival of oral squamous cell carcinoma

Background

DNA methylation has started a recent revolution in genomics biology by identifying key biomarkers for multiple cancers, including oral squamous cell carcinoma (OSCC), the most common head and neck squamous cell carcinoma.

Methods

A multi-stage screening strategy was used to identify DNA-methylation-based signatures for OSCC prognosis. We used The Cancer Genome Atlas (TCGA) data as training set which were validated in two independent datasets from Gene Expression Omnibus (GEO). The correlation between DNA methylation and corresponding gene expression and the prognostic value of the gene expression were explored as well.

Results

The seven DNA methylation CpG sites were identified which were significantly associated with OSCC overall survival. Prognostic signature, a weighted linear combination of the seven CpG sites, successfully distinguished the overall survival of OSCC patients and had a moderate predictive ability for survival [training set: hazard ratio (HR) = 3.23, P = 5.52 × 10⁻¹⁰, area under the curve (AUC) = 0.76; validation set 1: HR = 2.79, P = 0.010, AUC = 0.67; validation set 2: HR = 3.69, P = 0.011, AUC = 0.66]. Stratification analysis by human papillomavirus status, clinical stage, age, gender, smoking status, and grade retained statistical significance. Expression of genes corresponding to candidate CpG sites (AJAP1, SHANK2, FOXA2, MT1A, ZNF570, HOXC4, and HOXB4) was also significantly associated with patient’s survival. Signature integrating of DNA methylation, gene expression, and clinical information showed a superior ability for prognostic prediction (AUC = 0.78).

Conclusion

Prognostic signature integrated of DNA methylation, gene expression, and clinical information provides a better prognostic prediction value for OSCC patients than that with clinical information only.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0392-9) contains supplementary material, which is available to authorized users.

Epigenetic Bases of Aberrant Glycosylation in Cancer

In this review, the sugar portions of glycoproteins, glycolipids, and glycosaminoglycans constitute the glycome, and the genes involved in their biosynthesis, degradation, transport and recognition are referred to as “glycogenes“. The extreme complexity of the glycome requires the regulatory layer to be provided by the epigenetic mechanisms. Almost all types of cancers present glycosylation aberrations, giving rise to phenotypic changes and to the expression of tumor markers. In this review, we discuss how cancer-associated alterations of promoter methylation, histone methylation/acetylation, and miRNAs determine glycomic changes associated with the malignant phenotype. Usually, increased promoter methylation and miRNA expression induce glycogene silencing. However, treatment with demethylating agents sometimes results in silencing, rather than in a reactivation of glycogenes, suggesting the involvement of distant methylation-dependent regulatory elements. From a therapeutic perspective aimed at the normalization of the malignant glycome, it appears that miRNA targeting of cancer-deranged glycogenes can be a more specific and promising approach than the use of drugs, which broad target methylation/acetylation. A very specific type of glycosylation, the addition of GlcNAc to serine or threonine (O-GlcNAc), is not only regulated by epigenetic mechanisms, but is an epigenetic modifier of histones and transcription factors. Thus, glycosylation is both under the control of epigenetic mechanisms and is an integral part of the epigenetic code.

Gene promoter-associated CpG island hypermethylation in squamous cell carcinoma of the tongue

The present study was undertaken to explore and validate novel hypermethylated DNA regions in squamous cell carcinoma of the tongue (SCCT). Genome-wide methylation changes were identified by differential methylation hybridization (DMH) microarray and validated by bisulfite genome sequencing (BGS). The results were compared against datasets from The Cancer Genome Atlas head and neck squamous cell carcinoma (TCGA-HNSCC), Gene Expression Omnibus (GSE26549), and ArrayExpress (E-MTAB-1328). DMH identified 116 hypomethylated and 241 hypermethylated regions. Of the latter, 24 were localized to promoter or 5'-UTR regions. By BGS, promoter sequences of DAPK1, LRPPRC, RAB6C, and ZNF471 were significantly hypermethylated in tumors when compared with matched normal tissues (P < 0.0001). A TCGA-HNSCC dataset (516 cases of cancer and 50 normal tissue samples) further confirmed hypermethylation of DAPK1, RAB6C, and ZNF471. Sensitivity and specificity of methylation markers for a diagnosis of cancer were in the range of 70-100% in our study and from TCGA-HNSCC datasets, with an area under curve (AUC) of 0.83 and above. Kaplan-Meier survival analysis of TCGA-HNSCC expression data revealed that patients with low expressions of DAPK1, RAB6C, and ZNF471 showed poorer survival than patients with high expression (P = 0.02). Human papillomavirus (HPV) was found in 55% of cases, HPV16 being the predominant genotype. DAPK1 immunohistochemical staining was lower in SCCT than in normal buccal epithelial cells. This is the first study to report hypermethylation of LRPPRC, RAB6C, and ZNF471 in SCCT and its diagnostic and prognostic potentials in a specific head and neck squamous cell carcinoma.