DNA methylation biomarkers for head and neck squamous cell carcinoma

DNA methylation plays an important role in the etiology and pathogenesis of head and neck squamous cell carcinoma (HNSCC). The current study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) by a comprehensive bioinformatics analysis. In addition, we screened for DEGs affected by DNA methylation modification and further investigated their prognostic values for HNSCC. We included microarray data of DNA methylation (GSE25093 and GSE33202) and gene expression (GSE23036 and GSE58911) from Gene Expression Omnibus. Aberrantly methylated-DEGs were analyzed with R software. The Cancer Genome Atlas (TCGA) RNA sequencing and DNA methylation (Illumina HumanMethylation450) databases were utilized for validation. In total, 27 aberrantly methylated genes accompanied by altered expression were identified. After confirmation by The Cancer Genome Atlas (TCGA) database, 2 hypermethylated-low-expression genes (FAM135B and ZNF610) and 2 hypomethylated-high-expression genes (HOXA9 and DCC) were identified. A receiver operating characteristic (ROC) curve confirmed the diagnostic value of these four methylated genes for HNSCC. Multivariate Cox proportional hazards analysis showed that FAM135B methylation was a favorable independent prognostic biomarker for overall survival of HNSCC patients.     

A robust six-miRNA prognostic signature for head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) remains a major health problem worldwide. We aimed to identify a robust microRNA (miRNA)-based signature for predicting HNSCC prognosis. The miRNA expression profiles of HNSCC were obtained from The Cancer Genome Atlas (TCGA) database. The TCGA HNSCC cohort was randomly divided into the discovery and validation cohort. A miRNA-based prognostic signature was built up based on TGCA discovery cohort, and then further validated. The downstream targets of prognostic miRNAs were subjected to functional enrichment analyses. The role of miR-1229-3p, a prognosis-related miRNA, in tumorigenesis of HNSCC was further evaluated. A total of 305 significantly differentially expressed miRNAs were found between HNSCC samples and normal tissues. A six-miRNA prognostic signature was constructed, which exhibited a strong association with overall survival (OS) in the TCGA discovery cohort. In addition, these findings were successfully confirmed in TCGA validation cohort and our own independent cohort. The miRNA-based signature was demonstrated as an independent prognostic indicator for HNSCC. A risk signature-based nomogram model was constructed and showed good performance for predicting the OS for HNSCC. The functional analyses revealed that the downstream targets of these prognostic miRNAs were closely linked to cancer progression. Mechanistically, in vitro analysis revealed that miR-1229-3p played a tumor promoting role in HNSCC. In conclusion, our study has developed a robust miRNA-based signature for predicting the prognosis of HNSCC with high accuracy, which will contribute to improve the therapeutic outcome.     

DNA methylation regulates Microtubule-associated tumor suppressor 1 in human non-small cell lung carcinoma

Microtubule associated tumor suppressor 1 (MTUS1) has been recognized as a tumor suppressor gene in multiple cancers. However, the molecular mechanisms underlying the regulation of MTUS1 are yet to be investigated. This study aimed to clarify the significance of DNA methylation in silencing MTUS1 expression. We report that MTUS1 acts as tumor suppressor in non-small cell lung carcinoma (NSCLC). Analysis of in silico database and subsequent knockdown of DNMT1 suggested an inverse correlation between DNMT1 and MTUS1 function. Interestingly, increased methylation at MTUS1 promoter is associated with low expression of MTUS1. Treatment with DNA methyltransferases (DNMTs) inhibitor, 5-aza-2′-deoxycytidine (AZA) leads to both reduced promoter methylation accompanied with enrichment of H3K9Ac and enhanced MTUS1 expression. Remarkably, knockdown of MTUS1 showed increased proliferation and migration of NSCLC cells in contrast to diminished proliferation and migration, upon treatment with AZA. We concluded that low expression of MTUS1 correlates to DNA methylation and histone deacetylation in human NSCLC.     

Exploring and comparing of the gene expression and methylation differences between lung adenocarcinoma and squamous cell carcinoma

Lung cancer is one of the most frequently diagnosed malignant tumors and the main reason for cancer-related death around the world, whereas nonsmall cell lung cancer that consists two subtypes: lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) is responsible for an estimated 85% of all lung cancers. The current study aimed to explore gene expression and methylation differences between LUAD and LUSC. EdgeR was used to identify differentially regulated genes between normal and cancer in the LUAD and LUSC extracted from The Cancer Genome Atlas (TCGA), respectively, whereas SAM was used to find genes with differential methylation between normal and cancer in the LUAD and LUSC, respectively. Finally, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed to analyze the function which these genes enriched in. A total of 391 genes with opposite methylation patterns in LUAD and LUSC and four functional pathways were obtained (false discovery rate (FDR) < 0.1). These pathways mainly included fat digestion and absorption, phenylalanine metabolism, bile secretion, and so on, which were related to the airframe nutrition metabolic pathway. Moreover, two genes CTSE (cathepsin E) and solute carrier family 5 member 7 (SLC5A7) were also found, among which CTSE was overexpressed and hypomethylated in LUAD corresponding to normal lung tissues, whereas SLC5A7 showed the opposite in LUAD. In conclusion, this study investigated the differences between the gene expression and methylation patterns in LUAD and LUSC, and explored their different biological characteristics. Further understanding of these differences may promote the discovery and development of new, accurate strategies for the prevention, diagnosis, and treatment of lung cancer.     

Gene promoter methylation signature predicts survival of head and neck squamous cell carcinoma patients

Infection with high-risk types of human papilloma virus (HPV) is currently the best-established prognostic marker for head and neck squamous cell carcinoma (HNSCC), one of the most common and lethal human malignancies worldwide. Clinical trials have been launched to address the concept of treatment de-escalation for HPV-positive HNSCC with the final aim to reduce treatment related toxicity and debilitating long-term impacts on the quality of life. However, HPV-related tumors are mainly restricted to oropharyngeal SCC (OPSCC) and there is an urgent need to establish reliable biomarkers for all patients at high risk for treatment failure. A patient cohort (n = 295) with mainly non-OPSCC (72.9%) and a low prevalence of HPV16-related tumors (8.8%) was analyzed by MassARRAY to determine a previously established prognostic methylation score (MS). Kaplan-Meier revealed a highly significant correlation between a high MS and a favorable survival for OPSCC (P = 0.0004) and for non-OPSCC (P<0.0001), which was confirmed for all HNSCC by multivariate Cox regression models (HR: 9.67, 95% CI [4.61–20.30], P<0.0001). Next, we established a minimal methylation signature score (MMSS), which consists of ten most informative of the originally 62 CpG units used for the MS. The prognostic value of the MMSS was confirmed by Kaplan-Meier analysis for all HNSCC (P<0.0001) and non-OPSCC (P = 0.0002), and was supported by multivariate Cox regression models for all HNSCC (HR: 2.15, 95% CI [1.36–3.41], P = 0.001). In summary, the MS and the MMSS exhibit an excellent performance as prognosticators for survival, which is not limited by the anatomical site, and both could be implemented in future clinical trials.     

An 11-lncRNA expression could be potential prognostic biomarkers in head and neck squamous cell carcinoma

The aim of our study is to construct the competing endogenous RNA (ceRNA) network of head and neck squamous cell carcinoma (HNSCC) and identify key long noncoding RNAs (lncRNAs) to predict prognosis. The genes whose expression were differentially in HNSCC and normal tissues were explored by the Cancer Genome Atlas database. The ceRNA network was constructed by the Cytoscape software. The lncRNAs which could estimate the overall survival were explored from Cox proportional hazards regression. There are 1997, 589, and 82 mRNAs, lncRNAs, and miRNAs whose expression were statistically significant different, respectively. Then, the network between miRNA and mRNA or miRNA and lncRNA was constructed by miRcode, miRDB, TargetScan, and miRanda. Five mRNAs, 10 lncRNAs, and 3 miRNAs were associated with overall survival. Then, 11-lncRNAs were found to be prognostic factors. Therefore, our research analyzed the potential signature of novel 11-lncRNA as candidate prognostic biomarker from the ceRNA network for patients with HNSCC.     

Identification of 4-lncRNA prognostic signature in head and neck squamous cell carcinoma

Deregulated long noncoding RNAs (lncRNA) have been critically implicated in tumorigenesis and serve as novel diagnostic and prognostic biomarkers. Here we sought to develop a prognostic lncRNA signature in patients with head and neck squamous cell carcinoma (HNSCC). Original RNA-seq data of 499 HNSCC samples were retrieved from The Cancer Genome Atlas database, which was randomly divided into training and testing set. Univariate Cox regression survival analysis, robust likelihood-based survival model and random sampling iterations were applied to identify prognostic lncRNA candidates in the training cohort. A prognostic risk score was developed based on the Cox coefficient of four individual lncRNA imputed as follows: (0.14546 × expression level of RP11-366H4.1) + (0.27106 × expression level of LINC01123) + (0.54316 × expression level of RP11-110I1.14) + (−0.48794 × expression level of CTD-2506J14.1). Kaplan-Meier analysis revealed that patients with high-risk score had significantly reduced overall survival as compared with those with low-risk score when patients in training, testing, and validation cohorts were stratified into high- or low-risk subgroups. Multivariate survival analysis further revealed that this 4-lncRNA signature was a novel and important prognostic factor independent of multiple clinicopathological parameters. Importantly, ROC analyses indicated that predictive accuracy and sensitivity of this 4-lncRNA signature outperformed those previously well-established prognostic factors. Noticeably, prognostic score based on quantification of these 4-lncRNA via qRT-PCR in another independent HNSCC cohort robustly stratified patients into subgroups with high or low survival. Taken together, we developed a robust 4-lncRNA prognostic signature for HNSCC that might provide a novel powerful prognostic biomarker for precision oncology.     

Promoter methylation regulates cadherin switching in squamous cell carcinoma

Cadherins are cell adhesion molecules that modulate the epithelial phenotype and regulate tumor invasion. To identify the role of promoter methylation in regulating E-cadherin expression and in the "switching" of cadherins in oral squamous cell carcinoma (SCC), we studied 14 cell lines for cadherin expression. Immunoblotting revealed that only two (HOC-313 and HA-376) showed strong up-regulation of N-cadherin, and neither expressed E-cadherin. These results were confirmed by PCR. Furthermore, analysis of genomic DNA showed that the lack of E-cadherin expression in the two cell lines was not due to gene deletion. In both cell lines, methylation-specific PCR indicated extensive methylation of the 5' CpG island in the E-cadherin promoter. After treatment with a DNA methylation inhibitor (5-Aza-2-deoxycytidine), both immunoblotting and immunofluorescence staining showed that HA-376 cells newly expressed E-cadherin with a parallel decrease in their N-cadherin expression. Multiplex RT-PCR demonstrated that the down-regulation of N-cadherin mRNA was coordinately regulated with E-cadherin expression. Thus, methylation of the 5' CpG island in the E-cadherin promoter induces reciprocal expression of E- and N-cadherins in oral SCC by an unknown mechanism that appears to be mediated at the level of N-cadherin gene expression. These events may play an important role in the regulation of tumor cell mobility and invasion.     

A two-CpG-based prognostic signature for oral squamous cell carcinoma overall survival

Oral squamous cell carcinoma (OSCC) represents one of the most common head and neck cancer that with dire prognosis due partly to the lack of reliable prognostic biomarker. Here, we aimed to develop a CpG site–based prognostic signature through which we could accurately predict overall survival (OS) of patients with OSCC. We obtained OSCC-related DNA methylation and gene expression data sets from the public accessible Gene Expression Omnibus. Correlations between methylation level of CpG sites and OS of patients with OSCC were assessed by univariate Cox regression analysis followed by robust likelihood-based survival analysis on those CpG sites with permutation P < 0.05 for further screening the optimal CpG sites for OSCC OS prediction based on the risk score formula that composed of the methylation level of optimal CpG sites weighted by their regression coefficients. Besides, differential expression genes (DEGs) and differential methylation genes (DMGs) in OSCC samples compared with normal samples were obtained and shared genes were considered as vital genes in OSCC tumorgenesis and progression. As a result, two CpG sites including cg17892178 and cg17378966 that located in NID2 and IDO1, respectively, were identified as the optimal prognostic signatures for OSCC OS. In addition, 12 overlapping genes between DEGs and DMGs that closely associated with inflammation or blood and tissue development–related biological processes were obtained. In conclusions, this study should provide valuable signatures for OSCC diagnosis and treatment.     

Prediction of head and neck squamous cell carcinoma survival based on the expression of 15 lncRNAs

Recent evidence suggests that long noncoding RNAs (lncRNAs) are essential regulators of many cancer-related processes, including cancer cell proliferation, invasion, and migration. There is thus a reason to believe that the detection of lncRNAs may be useful as a diagnostic and prognostic strategy for cancer detection, however, at present no effective genome-wide tests are available for clinical use, constraining the use of such a strategy. In this study, we performed a comprehensive assessment of lncRNAs expressed in samples in the head and neck squamous cell carcinoma (HNSCC) cohort available in The Cancer Genome Atlas database. A risk score (RS) model was constructed based on the expression data of these 15 lncRNAs in the validation data set of HNSCC patients and was subsequently validated in validation data set and the entire data set. We were able to stratify patients into high- and low-risk categories, using our lncRNA expression panel to determine an RS, with significant differences in overall survival (OS) between these two groups in our test set (median survival, 1.863 vs. 5.484 years; log-rank test, p < 0.001). We were able to confirm the predictive value of our 15-lncRNA signature using both a validation data set and a full data set, finding our signature to be reproducible and effective as a means of predicting HNSCC patient OS. Through the multivariate Cox regression and stratified analyses, we were further able to confirm that the predictive value of this RS was independent of other predictive factors such as clinicopathological parameters. The Gene set enrichment analysis revealed potential functional roles for these 15 lncRNAs in tumor progression. Our findings indicate that an RS established based on a panel of lncRNA expression signatures can effectively predict OS and facilitate patient stratification in HNSCC.     

The 3D tissue microenvironment modulates DNA methylation and E-cadherin expression in squamous cell carcinoma

The microenvironment plays a significant role in human cancer progression. However, the role of the tumor microenvironment in the epigenetic control of genes critical to cancer progression remains unclear. As transient E-cadherin expression is central to many stages of neoplasia and is sensitive to regulation by the microenvironment, we have studied if microenvironmental control of E-cadherin expression is linked to transient epigenetic regulation of its promoter, contributing to the unstable and reversible expression of E-cadherin seen during tumor progression. We used 3D, bioengineered human tissue constructs that mimic the complexity of their in vivo counterparts, to show that the tumor microenvironment can direct the re-expression of E-cadherin through the reversal of methylation-mediated silencing of its promoter. This loss of DNA methylation results from the induction of homotypic cell-cell interactions as cells undergo tissue organization. E-cadherin re-expression is associated with multiple epigenetic changes including altered methylation of a small number of CpGs, specific histone modifications, and control of miR-148a expression. These epigenetic changes may drive the plasticity of E-cadherin-mediated adhesion in different tissue microenvironments during tumor cell invasion and metastasis. Thus, we suggest that epigenetic regulation is a mechanism through which tumor cell colonization of metastatic sites occurs as E-cadherin-expressing cells arise from E-cadherin-deficient cells.     

Analysis of DNA methylation and gene expression in radiation-resistant head and neck tumors

Resistance to radiation therapy constitutes a significant challenge in the treatment of head and neck squamous cell cancer (HNSCC). Alteration in DNA methylation is thought to play a role in this resistance. Here, we analyzed DNA methylation changes in a matched model of radiation resistance for HNSCC using the Illumina HumanMethylation450 BeadChip. Our results show that compared to radiation-sensitive cells (SCC-61), radiation-resistant cells (rSCC-61) had a significant increase in DNA methylation. After combining these results with microarray gene expression data, we identified 84 differentially methylated and expressed genes between these 2 cell lines. Ingenuity Pathway Analysis revealed ILK signaling, glucocorticoid receptor signaling, fatty acid α-oxidation, and cell cycle regulation as top canonical pathways associated with radiation resistance. Validation studies focused on CCND2, a protein involved in cell cycle regulation, which was identified as hypermethylated in the promoter region and downregulated in rSCC-61 relative to SCC-61 cells. Treatment of rSCC-61 and SCC-61 with the DNA hypomethylating agent 5-aza-2''deoxycitidine increased CCND2 levels only in rSCC-61 cells, while treatment with the control reagent cytosine arabinoside did not influence the expression of this gene. Further analysis of HNSCC data from The Cancer Genome Atlas found increased methylation in radiation-resistant tumors, consistent with the cell culture data. Our findings point to global DNA methylation status as a biomarker of radiation resistance in HNSCC, and suggest a need for targeted manipulation of DNA methylation to increase radiation response in HNSCC.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.     

DNA methylation profiles and biomarkers of oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) constitutes >90% of oral cancers and is the sixth most common malignancy among males worldwide and the fourth leading cause of death due to cancer among males in Taiwan. However, most patients do not receive a diagnosis of OSCC until the late stages, which have a lower survival rate. The use of molecular marker analysis to identify early-stage OSCC would permit optimal timing for treatments and consequently prolong survival. The aim of this study was to identify biomarkers of OSCC using the Illumina GoldenGate Methylation Cancer Panel, which comprised a total of 1,505 CpG sites covering 807 genes. Samples of buccal mucosa resected from 40 OSCC patients and normal tissue samples obtained from 15 patients (normal mucosa from OSCC patients or from patients undergoing surgery unrelated to OSCC) were analyzed. Fms-related tyrosine kinase 4 (FLT4) methylation exhibited a perfect specificity for detecting OSCC, with an area under the receiver operating characteristic curve of 0.91 for both all-stage and early-stage OSCC. Methylation of 7 genes (ASCL1, FGF3, FLT4, GAS7, KDR, TERT, and TFPI2) constitutes the top-20 panels for detecting OSCC. The top-20 panels for detecting early-stage OSCC contain 8 genes: ADCYAP1, EPHA7, FLT4, GSTM2, KDR, MT1A, NPY, and TFPI2. FLT4 RNA expression and methylation level were validated using RT-PCR and a pyrosequencing methylation assay. The median level of FLT4 expression was 2.14-fold for normal relative to OSCC tissue samples (P < 0.0001). Among the 8 pyrosequenced FLT4 CpG sites, methylation level was much higher in the OSCC samples. In conclusion, methylation statuses of selected genes, and especially FLT4, KDR, and TFPI2, might be of great potential as biomarkers for early detection of buccal OSCC.