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

Oral cancers preceded by proliferative verrucous leukoplakia exhibit distinctive molecular features

OBJECTIVE

Proliferative verrucous leukoplakia (PVL) has high rates of malignant transformation into oral squamous cell carcinoma (OSCC), but the clinical and evolutionary pattern of OSCC from PVL (PVL-OSCC) is more favorable than that of OSCC not preceded by PVL (OSCC). Here, we aimed to explore the pathophysiologic differences between PVL-OSCC and OSCC through transcriptomic and DNA methylation analyses.

MATERIALS AND METHODS

In this case-control study, oral biopsies from 8 PVL-OSCC and 10 OSCC patients were obtained for global sequencing using RNAseq and a genome-wide DNA methylation analysis via the Infinium EPIC Platform (graphical abstract).

RESULTS

One hundred and thirty-three differentially expressed genes (DEGs) were detected, 94 of them upregulated in OSCC. Most of these genes were previously described in cancer and associated with prognosis. The integrative analysis revealed 26 DEGs, corresponding to 37 CpGs, whose promoters were regulated by DNA methylation. Twenty-nine of the CpGs were found as hypermethylated in PVL-OSCC. Only 5 of the genes that were aberrantly methylated and differentially expressed were upregulated in PVL-OSCC patients, whereas 21 were underexpressed.

CONCLUSIONS

PVL-OSCC patients presented lower expression of cancer-related genes. Hypermethylation of the promoter region of many genes was also noticed, indicating that DNA methylation could be a regulatory mechanism.

Network Pharmacology and Molecular Docking Analysis Explores the Mechanisms of Cordyceps sinensis in the Treatment of Oral Lichen Planus

Objective

Oral lichen planus (OLP) is the most common potentially malignant disorder of the oral cavity. This study aimed to investigate the mechanism of action of Cordyceps sinensis in the treatment of OLP and provides a theoretical support for improving current treatment regimens for OLP.

Methods

The active components and therapeutic targets of Cordyceps sinensis were predicted and screened using the TCMSP, SymMap, PubMed, HIT 2.0, and PharmMapper databases, while the relevant OLP targets were predicted and screened using the DisGeNET and GeneCards databases. Protein-protein interactions (PPI) were examined using the String database, and Cytoscape was used to combine and illustrate the findings. GO and KEGG pathway enrichment analyses were carried out using RStudio, and AutoDock Vina and Pymol were used for molecular docking and visualization, respectively.

Results

A total of 404 potential target genes were discovered after evaluating 21 active compounds from Cordyceps sinensis. Potential therapeutic targets included 67 targets that matched and overlapped with OLP, including TNF, IL-6, CD4, EGFR, and IL1B. Key targets were predominantly engaged in the PI3K-Akt signaling pathway and the MAPK signaling pathway, according to the GO and KEGG analyses. These targets have a connection to biological processes including apoptosis signaling pathway regulation, T cell activation, and oxidative stress response. The molecular docking results showed that TNF, IL-6, CD4, EGFR, and IL1B could bind to their corresponding active components.

Conclusions

Cordyceps sinensis contains multiple components and acts on multiple targets and multiple pathways. Particularly, Cordyceps sinensis targets TNF, IL-6, CD4, EGFR, and IL1B, regulates PI3K-Akt and MAPK signaling pathways, as well as takes part in biological processes including apoptosis, T cell activation, and oxidative stress. Cordyceps sinensis could be a crucial choice in the therapy of OLP.

Evolutionary Relationships Between Dysregulated Genes in Oral Squamous Cell Carcinoma and Oral Microbiota

Oral squamous cell carcinoma (OSCC) is one of the most prevalent cancers in the world. Changes in the composition and abundance of oral microbiota are associated with the development and metastasis of OSCC. To elucidate the exact roles of the oral microbiota in OSCC, it is essential to reveal the evolutionary relationships between the dysregulated genes in OSCC progression and the oral microbiota. Thus, we interrogated the microarray and high-throughput sequencing datasets to obtain the transcriptional landscape of OSCC. After identifying differentially expressed genes (DEGs) with three different methods, pathway and functional analyses were also performed. A total of 127 genes were identified as common DEGs, which were enriched in extracellular matrix organization and cytokine related pathways. Furthermore, we established a predictive pipeline for detecting the coevolutionary of dysregulated host genes and microbial proteomes based on the homology method, and this pipeline was employed to analyze the evolutionary relations between the seven most dysregulated genes (MMP13, MMP7, MMP1, CXCL13, CRISPO3, CYP3A4, and CRNN) and microbiota obtained from the eHOMD database. We found that cytochrome P450 3A4 (CYP3A4), a member of the cytochrome P450 family of oxidizing enzymes, was associated with 45 microbes from the eHOMD database and involved in the oral habitat of Comamonas testosteroni and Arachnia rubra. The peptidase M10 family of matrix metalloproteinases (MMP13, MMP7, and MMP1) was associated with Lacticaseibacillus paracasei, Lacticaseibacillus rhamnosus, Streptococcus salivarius, Tannerella sp._HMT_286, and Streptococcus infantis in the oral cavity. Overall, this study revealed the dysregulated genes in OSCC and explored their evolutionary relationship with oral microbiota, which provides new insight for exploring the microbiota–host interactions in diseases.