Downregulation of keratin 76 expression during oral carcinogenesis of human, hamster and mouse

Identification of Biomarkers Associated with Cancerous Change in Oral Leukoplakia Based on Integrated Transcriptome Analysis

Objective

Oral leukoplakia (OLK) is the most common precancerous lesion in the oral cavity. This study aimed to explore key biomarkers for monitoring OLK for early diagnosis of oral squamous cell carcinoma (OSCC) and screen small-molecule drugs for the prevention of OSCC.

Method

The Gene Expression Omnibus (GEO) database was explored to extract two microarray datasets, namely, GSE85195 and GSE25099. The data of the normal group, OLK group, and OSCC group were analyzed by weighted gene coexpression network analysis (WGCNA) to identify the most significant gene module and differentially expressed genes (DEGs). The intersection genes were extracted as the key genes of OLK carcinogenesis. Subsequently, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed in the module. Connectivity Map and molecular docking were used to screen small-molecule drugs. The diagnostic values of four key genes were identified and verified in the GSE26549 dataset.

Results

WGCNA obtained the red module (r = −0.91, p < 0.05) with the strongest correlation with cancerous phenotype. GO enrichment analysis showed 60 pathways, including 28 biological processes, 11 cell components, and 21 molecular functions, and KEGG enrichment analysis showed 4 pathways (p < 0.05). In the differential expression analysis, there was no intersection between the upregulated genes and the red module genes. However, the intersection of the downregulated genes and the red module genes yielded 4 key genes: dopachrome tautomerase (DCT), keratin 3 (KRT3), keratin 76 (KRT76), and FAM3 metabolic regulation signal molecule B (FAM3B). The area under the curve of the diagnostic model constructed by these four genes was 0.963 (CI = 0.913–1.000). The sensitivity was 0.933, and the specificity was 0.923. The diagnostic model was successfully verified in GSE26549 (AUC = 0.745, CI = 0.638–0.851). Compared with the diagnostic models of the previous studies, the diagnostic efficiency of this model was the highest. The small-molecule drugs, selumetinib and benidipine, were selected according to the gene expression profile and showed binding activity when docking with the above molecules.

Conclusions

This study provides new targets and drugs for OLK. These targets could be used as the key diagnostic molecules for long-term follow-up of OLK. The small-molecule drugs selumetinib and benidipine could be used for the prevention and treatment of OSCC.

Integrated analysis and identification of nine-gene signature associated to oral squamous cell carcinoma pathogenesis

Oral squamous cell carcinoma (OSCC) is one of the leading cancers with poor disease survival rate. Herein, we explored molecular basis, in silico identification and in vitro verification of genes associated with OSCC. Five gene expression series including, GSE30784, GSE13601, GSE9844, GSE23558 and GSE37991 were screened for differentially expressed genes (DEGs). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched by cluster Profiler. Further, protein-protein interaction network was analysed and hub genes were verified. A total of 6476 (up-regulated: 2848; down-regulated: 3628) DEGs were identified among OSCC patients and healthy controls. Gene Ontology analysis indicated DEGs enrichment in cellular motility, invasion and adhesion processes. KEGG analysis revealed enrichment of PI3K-Akt signalling, focal adhesion and regulation of actin cytoskeleton pathways. Subsequently, nine DEGs including APP, EHMT1, ACACB, PCNA, PLAU, FST, HMGA2, LAMC2 and SPP1 were correlated with TCGA expression data along with significant association towards patient's survival, recognized as hub genes. This dysregulated mRNA signature of genes was validated in two OSCC cell lines with an anti-cancer agent, fisetin. Fisetin inhibited the expression of APP, EHMT1, PCNA, PLAU, FST, HMGA2, LAMC2, SPP1 and upregulated the expression of ACACB gene which were associated with growth inhibition of both the OSCC cell lines. The regulatory effect of fisetin supported crucial role of nine hub genes identified in OSCC. This study signified that hub genes and pathways might influence the aggressiveness of OSCC. Thus, the proposed hub genes could be potential diagnostic biomarker and drug targets for OSCC.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02737-4.

The Hamster Model of Sequential Oral Carcinogenesis: An Update

Animal models are valuable tools for studying human cancer as well as for preclinical trials. The hamster model of chemically induced sequential oral carcinogenesis was developed by our group a decade ago in order to study the multistep process of alterations in gene expression during carcinogenesis. The purpose of this review was to discuss the utility of the hamster model of sequential oral carcinogenesis regarding the deciphering of the main pathways altered. An extended search for articles that cited that specific animal models was performed. Many studies have used the hamster model of sequential oral carcinogenesis either for evaluation of the expression of biomarkers alone, or for applying chemopreventive compounds and other therapeutic methods, or combining the use of biomarkers with the anticancer effect of some compounds. It seems that this animal model is indeed a useful tool that enables the study of cell biology, pathology and therapeutics of oral cancer.

Mouse Tumor-Bearing Models as Preclinical Study Platforms for Oral Squamous Cell Carcinoma

Preclinical animal models of oral squamous cell carcinoma (OSCC) have been extensively studied in recent years. Investigating the pathogenesis and potential therapeutic strategies of OSCC is required to further progress in this field, and a suitable research animal model that reflects the intricacies of cancer biology is crucial. Of the animal models established for the study of cancers, mouse tumor-bearing models are among the most popular and widely deployed for their high fertility, low cost, and molecular and physiological similarity to humans, as well as the ease of rearing experimental mice. Currently, the different methods of establishing OSCC mouse models can be divided into three categories: chemical carcinogen-induced, transplanted and genetically engineered mouse models. Each of these methods has unique advantages and limitations, and the appropriate application of these techniques in OSCC research deserves our attention. Therefore, this review comprehensively investigates and summarizes the tumorigenesis mechanisms, characteristics, establishment methods, and current applications of OSCC mouse models in published papers. The objective of this review is to provide foundations and considerations for choosing suitable model establishment methods to study the relevant pathogenesis, early diagnosis, and clinical treatment of OSCC.

Multifaceted role of keratins in epithelial cell differentiation and transformation

Keratins, the epithelial-predominant members of the intermediate filament superfamily, are expressed in a pairwise, tissuespecific and differentiation-dependent manner. There are 28 type I and 26 type II keratins, which share a common structure comprising a central coiled coil α-helical rod domain flanked by two nonhelical head and tail domains. These domains harbor sites for major posttranslational modifications like phosphorylation and glycosylation, which govern keratin function and dynamics. Apart from providing structural support, keratins regulate various signaling machinery involved in cell growth, motility, apoptosis etc. However, tissue-specific functions of keratins in relation to cell proliferation and differentiation are still emerging. Altered keratin expression pattern during and after malignant transformation is reported to modulate different signaling pathways involved in tumor progression in a context-dependent fashion. The current review focuses on the literature related to the role of keratins in the regulation of cell proliferation, differentiation and transformation in different types of epithelia.

Heterogeneity of the Head and Neck Squamous Cell Carcinoma Immune Landscape and Its Impact on Immunotherapy

Head and neck squamous cell carcinomas (HNSCCs) are highly aggressive, multi-factorial tumors in the upper aerodigestive tract affecting more than half a million patients worldwide each year. Alcohol, tobacco, and human papillomavirus (HPV) infection are well known causative factors for HNSCCs. Current treatment options for HNSCCs are surgery, radiotherapy, chemotherapy, or combinatorial remedies. Over the past decade, despite the marked improvement in clinical outcome of many tumor types, the overall 5-year survival rate of HNSCCs remained ∼40–50% largely due to poor availability of effective therapeutic options for HNSCC patients with recurrent disease. Therefore, there is an urgent and unmet need for the identification of specific molecular signatures that better predict the clinical outcomes and markers that serve as better therapeutic targets. With recent technological advances in genomic and epigenetic analyses, our knowledge of HNSCC molecular characteristics and classification has been greatly enriched. Clinical and genomic meta-analysis of multicohort HNSCC gene expression profile has clearly demonstrated that HPV⁺ and HPV^- HNSCCs are not only derived from tissues of different anatomical regions, but also present with different mutation profiles, molecular characteristics, immune landscapes, and clinical prognosis. Here, we briefly review our current understanding of the biology, molecular profile, and immunological landscape of the HPV⁺ and HPV^- HNSCCs with an emphasis on the diversity and heterogeneity of HNSCC clinicopathology and therapeutic responses. After a review of recent advances and specific challenges for effective immunotherapy of HNSCCs, we then conclude with a discussion on the need to further enhance our understanding of the unique characteristics of HNSCC heterogeneity and the plasticity of immune landscape. Increased knowledge regarding the immunological characteristics of HPV⁺ and HPV^- HNSCCs would improve therapeutic targeting and immunotherapy strategies for different subtypes of HNSCCs.

Immunomodulatory role of Keratin 76 in oral and gastric cancer

Keratin 76 (Krt76) is expressed in the differentiated epithelial layers of skin, oral cavity and squamous stomach. Krt76 downregulation in human oral squamous cell carcinomas (OSCC) correlates with poor prognosis. We show that genetic ablation of Krt76 in mice leads to spleen and lymph node enlargement, an increase in regulatory T cells (Tregs) and high levels of pro-inflammatory cytokines. Krt76^-/- Tregs have increased suppressive ability correlated with increased CD39 and CD73 expression, while their effector T cells are less proliferative than controls. Loss of Krt76 increases carcinogen-induced tumours in tongue and squamous stomach. Carcinogenesis is further increased when Treg levels are elevated experimentally. The carcinogenesis response includes upregulation of pro-inflammatory cytokines and enhanced accumulation of Tregs in the tumour microenvironment. Tregs also accumulate in human OSCC exhibiting Krt76 loss. Our study highlights the role of epithelial cells in modulating carcinogenesis via communication with cells of the immune system.

Raman Spectroscopy of Experimental Oral Carcinogenesis: Study on Sequential Cancer Progression in Hamster Buccal Pouch Model

Oral cancers suffer from poor 5-year survival rates, owing to late detection of the disease. Current diagnostic/screening tools need to be upgraded in view of disadvantages like invasiveness, tedious sample preparation, long output times, and interobserver variances. Raman spectroscopy has been shown to identify many disease conditions, including oral cancers, from healthy conditions. Further studies in exploring sequential changes in oral carcinogenesis are warranted. In this Raman spectroscopy study, sequential progression in experimental oral carcinogenesis in Hamster buccal pouch model was investigated using 3 approaches-ex vivo, in vivo sequential, and in vivo follow-up. In all these studies, spectral changes show lipid dominance in early stages while later stages and tumors showed increased protein to lipid ratio and nucleic acids. On similar lines, early weeks of 7,12-dimethylbenz(a)anthracene-treated and control groups showed higher overlap and low classification. The classification efficiency increased progressively, reached a plateau phase and subsequently increased up to 100% by 14 weeks. The misclassifications between treated and control spectra suggested some changes in controls as well, which was confirmed by a careful reexamination of histopathological slides. These findings suggests Raman spectroscopy may be able to identify microheterogeneity, which may often go unnoticed in conventional biochemistry wherein tissue extracts are employed, as well as in histopathology. In vivo findings, quite comparable to gold-standard supported ex vivo findings, give further proof of Raman spectroscopy being a promising label-free, noninvasive diagnostic adjunct for future clinical applications.

Overexpression of lipocalins and pro-inflammatory chemokines and altered methylation of PTGS2 and APC2 in oral squamous cell carcinomas induced in rats by 4-nitroquinoline-1-oxide

Oral squamous cell carcinomas (OSCC) induced in F344 rats by 4-nitroquinoline-1-oxide (4-NQO) demonstrate considerable phenotypic similarity to human oral cancers. Gene expression studies (microarray and PCR) were coupled with methylation analysis of selected genes to identify molecular markers of carcinogenesis in this model and potential biochemical and molecular targets for oral cancer chemoprevention. Microarray analysis of 11 pairs of OSCC and site-matched phenotypically normal oral tissues from 4-NQO-treated rats identified more than 3500 differentially expressed genes; 1735 genes were up-regulated in rat OSCC versus non-malignant tissues, while 1803 genes were down-regulated. In addition to several genes involved in normal digestion, genes demonstrating the largest fold increases in expression in 4-NQO-induced OSCC include three lipocalins (VEGP1, VEGP2, LCN2) and three chemokines (CCL, CXCL2, CXCL3); both classes are potentially druggable targets for oral cancer chemoprevention and/or therapy. Down-regulated genes in 4-NQO-induced OSCC include numerous keratins and keratin-associated proteins, suggesting that alterations in keratin expression profiles may provide a useful biomarker of oral cancer in F344 rats treated with 4-NQO. Confirming and extending our previous results, PTGS2 (cyclooxygenase-2) and several cyclooxygenase-related genes were significantly up-regulated in 4-NQO-induced oral cancers; up-regulation of PTGS2 was associated with promoter hypomethylation. Rat OSCC also demonstrated increased methylation of the first exon of APC2; the increased methylation was correlated with down-regulation of this tumor suppressor gene. Overexpression of pro-inflammatory chemokines, hypomethylation of PTGS2, and hypermethylation of APC2 may be causally linked to the etiology of oral cancer in this model.

Keratin 76 is required for tight junction function and maintenance of the skin barrier

Keratins are cytoskeletal intermediate filament proteins that are increasingly being recognised for their diverse cellular functions. Here we report the consequences of germ line inactivation of Keratin 76 (Krt76) in mice. Homozygous disruption of this epidermally expressed gene causes neonatal skin flaking, hyperpigmentation, inflammation, impaired wound healing, and death prior to 12 weeks of age. We show that this phenotype is associated with functionally defective tight junctions that are characterised by mislocalization of the integral protein CLDN1. We further demonstrate that KRT76 interacts with CLDN1 and propose that this interaction is necessary to correctly position CLDN1 in tight junctions. The mislocalization of CLDN1 has been associated in various dermopathies, including the inflammatory disease, psoriasis. These observations establish a previously unknown connection between the intermediate filament cytoskeleton network and tight junctions and showcase Krt76 null mice as a possible model to study aberrant tight junction driven skin diseases.

Novel skin phenotypes revealed by a genome-wide mouse reverse genetic screen

Permanent stop-and-shop large-scale mouse mutant resources provide an excellent platform to decipher tissue phenogenomics. Here we analyse skin from 538 knockout mouse mutants generated by the Sanger Institute Mouse Genetics Project. We optimize immunolabelling of tail epidermal wholemounts to allow systematic annotation of hair follicle, sebaceous gland and interfollicular epidermal abnormalities using ontology terms from the Mammalian Phenotype Ontology. Of the 50 mutants with an epidermal phenotype, 9 map to human genetic conditions with skin abnormalities. Some mutant genes are expressed in the skin, whereas others are not, indicating systemic effects. One phenotype is affected by diet and several are incompletely penetrant. In-depth analysis of three mutants, Krt76, Myo5a (a model of human Griscelli syndrome) and Mysm1, provides validation of the screen. Our study is the first large-scale genome-wide tissue phenotype screen from the International Knockout Mouse Consortium and provides an open access resource for the scientific community.

Large-scale efforts are put into the generation of knockout mutant mice for many individual genes. Here, the authors systematically screen skin from 538 mutant mice and identify 50 mutants with epidermal phenotypes, of which 9 are also associated with human skin defects.