An integrated methylation and gene expression microarray analysis reveals significant prognostic biomarkers in oral squamous cell carcinoma

Inhibition of the Tumor Suppressor Gene SPINK5 via EHMT2 Induces the Oral Squamous Cell Carcinoma Development

Serine protease inhibitor Kazal-type 5 (SPINK5) has been revealed as a significant prognostic biomarker in oral squamous cell carcinoma (OSCC). However, there is little information regarding the detailed epigenetics mechanism underlying its dysregulation in OSCC. Using the Gene Expression Omnibus database, we identified SPINK5 as a significantly downregulated gene in OSCC tissues. Moreover, SPINK5 inhibited the malignant aggressiveness of HSC3 and squamous cell carcinomas (SCC)9 cells, whereas depletion of SPINK5 using shRNAs led to the opposite trend. The euchromatic histone lysine methyltransferase 2 (EHMT2) was found to bind to the SPINK5 promoter, and EHMT2 repressed the SPINK5 expression. SPINK5 reversed the stimulating effects of EHMT2 on the aggressiveness of HSC3 and SCC9 cells by impairing the Wnt/β-catenin pathway. Wnt/β-catenin inhibitor IWR-1 treatment reverted the malignant phenotype of OSCC cells in the presence of short hairpin RNA (sh)-SPINK5. Silencing of EHMT2 inhibited tumor growth and blocked the Wnt/β-catenin signaling in OSCC, which was reversed by SPINK5 knockdown. Our study shows that SPINK5, mediated by the loss of EHMT2, can inhibit the development of OSCC by inhibiting Wnt/β-catenin signaling and may serve as a treatment target for OSCC.

The repertoire of mutational signatures in tobacco- and non-tobacco-induced oral cancer

The use of tobacco products is one of the established contributors toward the development and spread of oral cancer. Additionally, recent research has indicated oral microbiome, infections with Human papilloma virus (HPV), Epstein-Barr virus (EBV), Candida as significant contributing factors to this disease along with lifestyle habits. Deregulation of cellular pathways envisaging metabolism, transcription, translation, and epigenetics caused by these risk factors either individually or in unison is manifold, resulting in the increased risk of oral cancer. Globally, this cancer continues to exist as one of the major causes of cancer-related mortalities; the numbers in the developing South Asian countries clearly indicate yearly escalation. This review encompasses the variety of genetic modifications, including adduct formation, mutation (duplication, deletion, and translocation), and epigenetic changes evident in oral squamous cell carcinoma (OSCC). In addition, it highlights the interference caused by tobacco products in Wnt signaling, PI3K/Akt/mTOR, JAK-STAT, and other important pathways. The information provided also ensures a comprehensive and critical revisit to non-tobacco-induced OSCC. Extensive literature survey and analysis has been conducted to generate the chromosome maps specifically highlighting OSCC-related mutations with the potential to act as spectacles for the early diagnosis and targeted treatment of this disease cancer.

Epigenetic methylation changes: implication as biomarkers in oral and maxillofacial area cancers

Background/Aim

Squamous cell carcinoma (SCC) is the most frequent cancer of the head and neck area in the oral cavity. Epigenetic alterations in oral and maxillofacial area cancers are urgently needed to be investigated, as the observed changes might have crucial diagnostic value for personalized medicine.

Methods

Our study aimed to identify the most frequently hypermethylated tumor suppressor gene promoters in OSCC, followed by correlation analysis with the patients' survival. We evaluated the methylation status of the promoters in a panel of 22 tumor suppressor genes in Romanian (n=9) and Bulgarian (n=12) patient groups suffering from oral and maxillofacial area cancers. The extracted DNA was further digested through EpiTect Methyl II PCR Array System containing methylation-sensitive and methylation-dependent restriction enzymes, followed by specific amplification of the products obtained by qPCR and data analysis using the online platform provided by the producer.

Results

Different methylation patterns were observed in the tumor suppressor genes' promoters. Among them, the methylation profile of Cccnd2, Chd1, Cdh13, Cdkn1c, Neurog1, Gstp1, and Runx3 genes further correlated with overall survival rates.

Conclusions

Our data emphasize that epigenetic alterations are responsible for the clinical heterogeneity of oral and maxillofacial area cancers and significantly impact on patient survival. Additional investigation on a larger patient cohort should validate these potential biomarkers.

Construction and analysis of a conjunctive diagnostic model of HNSCC with random forest and artificial neural network

Head and neck squamous cell carcinoma (HNSCC) is a heterogeneous tumor that is highly aggressive and ranks fifth among the most common cancers worldwide. Although, the researches that attempted to construct a diagnostic model were deficient in HNSCC. Currently, the gold standard for diagnosing head and neck tumors is pathology, but this requires a traumatic biopsy. There is still a lack of a noninvasive test for such a high-incidence tumor. In order to screen genetic markers and construct diagnostic model, the methods of random forest (RF) and artificial neural network (ANN) were utilized. The data of HNSCC gene expression was accessed from Gene Expression Omnibus (GEO) database; we selected three datasets totally, and we combined 2 datasets (GSE6631 and GSE55547) for screening differentially expressed genes (DEGs) and chose another dataset (GSE13399) for validation. Firstly, the 6 DEGs (CRISP3, SPINK5, KRT4, MMP1, MAL, SPP1) were screened by RF. Subsequently, ANN was applied to calculate the weights of 6 genes. Besides, we created a diagnostic model and nominated it as neuralHNSCC, and the performance of neuralHNSCC by area under curve (AUC) was verified using another dataset. Our model achieved an AUC of 0.998 in the training cohort, and 0.734 in the validation cohort. Furthermore, we used the Cell-type Identification using Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm to investigate the difference in immune cell infiltration between HNSCC and normal tissues initially. The selected 6 DEGs and the constructed novel diagnostic model of HNSCC would make contributions to the diagnosis.

A prognostic model for oral squamous cell carcinoma using 7 genes related to tumor mutational burden

Background

Oral squamous cell carcinoma (OSCC) is a rising problem in global public health. The traditional physical and imageological examinations are invasive and radioactive. There is a need for less harmful new biomarkers. Tumor mutational burden (TMB) is a novel prognostic biomarker for various cancers. We intended to explore the relationship between TMB-related genes and the prognosis of OSCC and to construct a prognostic model.

Methods

TMB-related differential expressed genes (DEGs) were screened by differential analysis and optimized via the univariate Cox and LASSO Cox analyses. Risk Score model was constructed by expression values of screened genes multiplying coefficient of LASSO Cox.

Results

Seven TMB-related DEGs (CTSG, COL6A5, GRIA3, CCL21, ZNF662, TDRD5 and GSDMB) were screened. Patients in high-risk group (Risk Score >  − 0.684511507) had worse prognosis compared to the low-risk group (Risk Score <  − 0.684511507). Survival rates of patients in the high-risk group were lower in the gender, age and degrees of differentiation subgroups compared to the low-risk group.

Conclusions

The Risk Score model constructed by 7 TMB-related genes may be a reliable biomarker for predicting the prognosis of OSCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12903-022-02193-3.

The Role of Zinc Finger Proteins in Various Oral Conditions

The zinc finger proteins (ZNFs) are essential transcription factors, and the genes encoding them constitute about 3% of the entire human genome. They are involved in the development of several tissues, and any alterations in their structure may promote chronic conditions like diabetes and tumorigenesis. Lately, their role in the development, progression, and metastasis of Oral Squamous Cell Carcinoma (OSCC), Epithelial Dysplasia, Oral Lichen Planus, and Periodontitis has been found. The present review aims to describe their role in various oral conditions. Electronic databases like Medline (PubMed) and Scopus were searched for original studies related to the role of ZNFs in various oral conditions. It yielded 48 studies included in the review. It was found that the ZNFs influenced chronic conditions like Oral Cancer and Periodontitis. They act both as tumor suppressors and oncogenes and have an anti-inflammatory effect. The knowledge from the present review may be utilized in designing drugs that prevent unusual expression of specific ZNFs. Besides, they may be applied as prognostic markers due to their high expression specificity in some tumors.

SPINKs in Tumors: Potential Therapeutic Targets

The serine protease inhibitor Kazal type (SPINK) family includes SPINK1-14 and is the largest branch in the serine protease inhibitor family. SPINKs play an important role in pancreatic physiology and disease, sperm maturation and capacitation, Nager syndrome, inflammation and the skin barrier. Evidence shows that the unregulated expression of SPINK1, 2, 4, 5, 6, 7, and 13 is closely related to human tumors. Different SPINKs exhibit various regulatory modes in different tumors and can be used as tumor prognostic markers. This article reviews the role of SPINK1, 2, 4, 5, 6, 7, and 13 in different human cancer processes and helps to identify new cancer treatment targets.

Integrative Multi-Omics Analysis Reveals Candidate Biomarkers for Oral Squamous Cell Carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of cancer worldwide. Due to the lack of early detection and treatment, the survival rate of OSCC remains poor and the incidence of OSCC has not decreased during the past decades. To explore potential biomarkers and therapeutic targets for OSCC, we analyzed differentially expressed genes (DEGs) associated with OSCC using RNA sequencing technology. Methylation-regulated and differentially expressed genes (MeDEGs) of OSCC were further identified via an integrative approach by examining publicly available methylomic datasets together with our transcriptomic data. Protein-protein interaction (PPI) networks of MeDEGs were constructed and highly connected hub MeDEGs were identified from these PPI networks. Subsequently, expression and survival analyses of hub genes were performed using The Cancer Genome Atlas (TCGA) database and the Gene Expression Profiling Interactive Analysis (GEPIA) online tool. A total of 56 upregulated MeDEGs and 170 downregulated MeDEGs were identified in OSCC. Eleven hub genes with high degree of connectivity were picked out from the PPI networks constructed by those MeDEGs. Among them, the expression level of four hub genes (CTLA4, CDSN, ACTN2, and MYH11) were found to be significantly changed in the head and neck squamous carcinoma (HNSC) patients. Three hypomethylated hub genes (CTLA4, GPR29, and TNFSF11) and one hypermethylated hub gene (ISL1) were found to be significantly associated with overall survival (OS) of HNSC patients. Therefore, these hub genes may serve as potential DNA methylation biomarkers and therapeutic targets of OSCC.

Resection Margins in Head and Neck Cancer Surgery: An Update of Residual Disease and Field Cancerization

Surgery is one of the mainstays of head and neck cancer treatment, and aims at radical resection of the tumor with 1 cm tumor-free margins to obtain locoregional control. Surgical margins are evaluated by histopathological examination of the resection specimen. It has been long an enigma that approximately 10-30% of surgically treated head and neck cancer patients develop locoregional recurrences even though the resection margins were microscopically tumor-free. However, the origins of these recurrences have been elucidated by a variety of molecular studies. Recurrences arise either from minimal residual disease, cancer cells in the surgical margins that escape detection by the pathologist when examining the specimen, or from precancerous mucosal changes that may remain unnoticed. Head and neck tumors develop in mucosal precursor changes that are sometimes visible but mostly not, fueling research into imaging modalities such as autofluorescence, to improve visualization. Mostly unnoticed, these precancerous changes may stay behind when the tumor is resected, and subsequent malignant progression will cause a local relapse. This led to a clinical trial of autofluorescence-guided surgery, of which the results were reported in 2020. This review focuses on the most recent literature of the improved diagnosis of the resection margins of surgically treated head and neck cancer patients, the pathobiological origin of recurrent disease, and relevant biomarkers to predict local relapse. Directions for further research will be discussed, including potential options for improved and personalized treatment, based on the most recently published data.

Bioinformatics analysis reveals that ANXA1 and SPINK5 are novel tumor suppressor genes in patients with oral squamous cell carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is a solid tumor of squamous epithelial origin. Currently, surgery is still the main treatment for OSCC, with radiotherapy and chemotherapy as important adjuvant treatments. However, the problem of poor prognosis of OSCC patients still exists in clinical practice. To explore further potential biomarkers or treatment targets in OSCC patients, this study used a high-throughput gene expression database to study the potential molecular mechanisms of OSCC carcinogenesis.

METHODS

The GEO database related to OSCC was searched and analyzed using GEO2R. Oncomine and the Human Protein Atlas were used to evaluate the expression level of differentially-expressed genes (DEGs). The cBioPortal dataset was used to analyze the mutations of the potential DEGs and patient survival.

RESULTS

Three GEO datasets, GSE146483, GSE138206, and GSE148944, were downloaded and 7 DEGs were found in common in OSCC tissues. Using Oncomine and the Human Protein Atlas, ANXA1, IL1RN, and SPINK5 were decreased in cancer tissues, while protein levels of APOE and IFI35 were increased accordingly. Interestingly, low levels of ANXA1 and SPINKS were associated with the TNM stage of OSCC patients. No mutations in DEGs were found in OSCC patients, based on the cBioPortal dataset. Survival analysis indicated OSCC patients with high MSR1 had poor overall survival (OS), while low expression of CXCR4, ANXA1, IL1RN, and SPINK5 also predicted poor OS in OSCC patients.

CONCLUSIONS

Our findings uncovered 7 potential biomarkers of OSCC patients, with ANXA1 and SPINK5 serving as potential tumor suppressor genes in OSCC.

Integrative Analysis of lncRNA-mRNA Coexpression in Human Lung Epithelial Cells Exposed to Dimethyl Selenide-Derived Secondary Organic Aerosols

Dimethyl selenide (DMSe) is one of the major volatile organoselenium compounds released into the atmosphere through plant metabolism and microbial methylation. DMSe has been recently revealed as a precursor of secondary organic aerosol (SOA), and its resultant SOA possesses strong oxidizing capability toward thiol groups that can perturb several major biological pathways in human airway epithelial cells and is linked to genotoxicity, DNA damage, and p53-mediated stress responses. Mounting evidence has suggested that long noncoding RNAs (lncRNAs) are involved in stress responses to internal and environmental stimuli. However, the underlying molecular interactions remain to be elucidated. In this study, we performed integrative analyses of lncRNA-mRNA coexpression in the transformed human bronchial epithelial BEAS-2B cell line exposed to DMSe-derived SOA. We identified a total of 971 differentially expressed lncRNAs in BEAS-2B cells exposed to SOA derived from O₃ and OH oxidation of DMSe. Gene ontology (GO) network analysis of cis-targeted genes showed significant enrichment of DNA damage, apoptosis, and p53-mediated stress response pathways. trans-Acting lncRNAs, including PINCR, PICART1, DLGAP1-AS2, and LINC01629, known to be associated with human carcinogenesis, also showed altered expression in cell treated with DMSe-SOA. Overall, this study highlights the regulatory role of lncRNAs in altered gene expression induced by DMSe-SOA exposure.

MMP1 3'UTR facilitates the proliferation and migration of human oral squamous cell carcinoma by sponging miR-188-5p to up-regulate SOX4 and CDK4

Growing evidence indicates that the non-coding 3'-untranslated region (3'UTR) of genes acts as competing endogenous RNAs (ceRNAs) to exert their roles in a number of diseases, including cancer. In the present study, MMP1 messenger RNA was identified to be significantly up-regulated in oral squamous cell carcinoma (OSCC) tissues, and both MMP1 and its 3'UTR promoted tumor growth and cell motility. Further mechanism investigations indicated that MMP1 3'UTR was able to antagonize miR-188-5p; in addition, overexpression of MMP1 3'UTR up-regulated the expression level of SOX4 and CDK4, target genes of miR-188-5p, which have also been identified as oncogenic driver genes in OSCC. Therefore, a ceRNA regulatory network among MMP1, SOX4, and CDK4 mediated via competing for binding to miR-188-5p was proved. Taken together, the present study demonstrates for the first time that MMP1 mRNA participates in the development of OSCC via ceRNA regulatory mechanism and genes involved in the ceRNA network may provide a novel avenue for target therapy.

An Integrated Microarray Analysis Reveals Significant Diagnostic and Prognostic Biomarkers in Pancreatic Cancer

Background

Pancreatic cancer (PAC) is a lethal cancer and it is essential to develop accurate diagnostic and prognostic biomarkers for PAC.

Material/Methods

An integrated microarray analysis of PAC was conducted to identify differentially expressed genes (DEGs) between PAC and non-tumor controls. Expression of DEGs were further confirmed by The Cancer Genome Atlas and the Genotype-Tissue Expression. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis, and protein–protein integration network construction were performed to further research the biological functions of DEGs. Receiver-operating characteristic analysis and survival analysis were used to evaluate the diagnostic and prognostic value of DEGs for PAC.

Results

Seventeen microarray datasets were downloaded from Gene Expression Omnibus to conduct the integrated microarray analysis. A total of 1136 DEGs (596 upregulated and 540 downregulated DEGs) in PAC tissues compared with non-tumor controls were identified. Pancreatic secretion (Kegg: 04972), insulin signaling pathway (Kegg: 04910), and several cancer-related pathways including pathways in cancer (Kegg: 05200), MAPK signaling pathway (Kegg: 04010), and pancreatic cancer (Kegg: 05212) were enriched for DEGs in PAC. Seven DEGs (AHNAK2, CDH3, IFI27, ITGA2, LAMB3, SLC6A14, and TMPRSS4) were found to have both great diagnostic and prognostic value for PAC. High expression of these 7 DEGs were significantly associated with poor prognosis of patients with PAC.

Conclusions

These 7 DEGs might be potential diagnostic and prognostic biomarkers for PAC and help uncovering the mechanism of PAC.

A Novel Tumor Suppressor SPINK5 Serves as an Independent Prognostic Predictor for Patients with Head and Neck Squamous Cell Carcinoma

BACKGROUND

In our previous study, serine protease inhibitor Kazal-type 5 (SPINK5), which encodes the product of serine protease inhibitor lymphoepithelial Kazal-type-related inhibitor (LEKTI) was found to be down-regulated in head and neck squamous cell carcinoma (HNSCC) using oligonucleotide microarrays. However, the function and clinical implications of SPINK5/LEKTI remain obscure in HNSCC.

METHODS

The endogenous expression level of SPINK5/LEKTI was further verified in 9 HNSCC cell lines and HNSCCs by means of reverse transcription-polymerase chain reaction, real-time PCR, Western blotting and immunohistochemistry. The biological function of SPINK5/LEKTI was investigated in vitro and in vivo experiments. Kaplan-Meier survival analysis and Cox proportional hazards regression model were used to determine the correlation between SPINK5/LEKTI expression and clinical outcome.

RESULTS

Down-regulation expression of SPINK5/LEKTI was found in six out of nine HNSCC cell lines and in 85.7% HNSCC specimens (P<0.0001). Upon silencing of SPINK5/LEKTI, the cell proliferation, plate colony formation and cell invasion of WU-HN6 cells were significantly increased, while exogenous overexpression of SPINK5/LEKTI, the proliferation, plate colony and invasion of WU-HN13 and HN30 cells were remarkably inhibited with the arrest of G1 cell cycle (P=0.0001, P=0.003, respectively). HNSCC patients with lower LEKTI levels had significantly inferior overall survival compared to those patients with higher LEKTI (P=0.0017) by Kaplan-Meier survival analysis. Univariate and multivariate Cox proportional hazards regression model analysis revealed that LEKTI expression was an independent prognostic predictor for HNSCC patients (HR=0.114, 95% CI:0.044-0.292, P<0.001).

CONCLUSION

Our results demonstrate that SPINK5/LEKTI might be a tumor suppressor in HNSCCs and serve as an independent prognostic predictor for HNSCC patients.

A panel of Transcription factors identified by data mining can predict the prognosis of head and neck squamous cell carcinoma

Background

Transcription factors (TFs) are responsible for the regulation of various activities related to cancer like cell proliferation, invasion, and migration. It is thought that, the measurement of TFs levels could assist in developing strategies for diagnosis and prognosis of cancer detection. However, due to lack of effective genome-wide tests, this cannot be carried out in clinical settings.

Methods

A complete assessment of RNA-seq data in samples of a head and neck squamous cell carcinoma (HNSCC) cohort in The Cancer Genome Atlas (TCGA) database was carried out. From the expression data of six TFs, a risk score model was developed and further validated in the GSE41613 and GSE65858 series. Potential functional roles were identified for the six TFs via gene set enrichment analysis.

Results

Based on our multi-TF signature, patients are stratified into high- and low-risk groups with significant variations in overall survival (OS) (median survival 2.416 vs. 5.934 years, log-rank test P < 0.001). The sensitivity and specificity evaluation of our multi-TF for 3-year OS in TCGA, GSE41613 and GSE65858 was 0.707, 0.679 and 0.605, respectively, demonstrating good reproducibility and robustness for predicting overall survival of HNSCC patients. Through multivariate Cox regression analyses (MCRA) and stratified analyses, we confirmed that the predictive capability of this risk score (RS) was not dependent on any of other factors like clinicopathological parameters.

Conclusions

With the help of a RS obtained from a panel of TFs expression signatures, effective OS prediction and stratification of HNSCC patients can be carried out.

Transcriptional Analysis of Immunohistochemically Defined Subgroups of Non-Muscle-Invasive Papillary High-Grade Upper Tract Urothelial Carcinoma

Immunohistochemical (IHC) staining for CK5/6 and CK20 was reported to be correlated with the prognosis of early urothelial carcinoma in a way contrary to that of advanced tumors for unknown reasons. We aimed to characterize the gene expression profiles of subgroups of non-muscle-invasive papillary high-grade upper tract urothelial carcinoma (UTUC) classified by CK5/6 and CK20 expression levels: group 1 (CK5/6-high/CK20-low), group 2 (CK5/6-high/CK20-high), and group 3 (CK5/6-low/CK20-high). Expression of group 3 was predictive of worse prognosis of non-muscle-invasive papillary high-grade UTUC. Transcriptional analysis revealed 308 differentially expressed genes across the subgroups. Functional analyses of the genes identified cell adhesion as a common process differentially enriched in group 3 compared to the other groups, which could explain its high-risk phenotype. Late cell cycle/proliferation signatures were also enriched in group 3 and in some of the other groups, which may be used as a prognostic biomarker complementary to CK5/6 and CK20. Group 2, characterized by low levels of genes associated with mitogen-activated protein kinase and tumor necrosis factor signaling pathways, was hypothesized to represent the least cancerous subtype considering its normal urothelium-like IHC pattern. This study would facilitate the application of easily accessible prognostic biomarkers in practice.

Investigation of the Prognostic Role of Carbonic Anhydrase 9 (CAIX) of the Cellular mRNA/Protein Level or Soluble CAIX Protein in Patients with Oral Squamous Cell Carcinoma

Carbonic anhydrase 9 (CAIX) is an important protein that stabilizes the extracellular pH value and is transcriptionally regulated by hypoxia-inducible factor 1 (HIF1), but more stable than HIF1α. Here we show a comparative study that examines the prognostic value of CA9 mRNA, CAIX protein of tumor cells and secreted CAIX protein for oral squamous cell carcinoma (OSCC) patients. Tumor samples from 72 OSCC patients and 24 samples of normal tissue were analyzed for CA9 mRNA levels. A total of 158 OSCC samples were stained for CAIX by immunohistochemistry and 89 blood serum samples were analyzed by ELISA for soluble CAIX protein content. Survival analyses were performed by Kaplan–Meier and Cox’s regression analysis to estimate the prognostic effect of CA9/CAIX in OSCC patients. The CA9 mRNA and CAIX protein levels of tumor cells correlated with each other, but not with those of the secreted CAIX protein level of the blood of patients. ROC curves showed a significant (p < 0.001) higher mRNA-level of CA9 in OSCC samples than in adjacent normal tissue. Cox’s regression analysis revealed an increased risk (i) of death for patients with a high CA9 mRNA level (RR = 2.2; p = 0.02), (ii) of locoregional recurrence (RR = 3.2; p = 0.036) at higher CA9 mRNA levels and (iii) of death at high CAIX protein level in their tumors (RR = 1.7; p = 0.066) and especially for patients with advanced T4-tumors (RR = 2.0; p = 0.04). However, the secreted CAIX protein level was only as a trend associated with prognosis in OSCC (RR = 2.2; p = 0.066). CA9/CAIX is an independent prognostic factor for OSCC patients and therefore a potential therapeutic target.