Artificial intelligence-based epigenomic, transcriptomic and histologic signatures of tobacco use in oral squamous cell carcinoma.

Oral squamous cell carcinoma (OSCC) biomarker studies rarely employ multi-omic biomarker strategies and pertinent clinicopathologic characteristics to predict mortality. In this study we determine for the first time a combined epigenetic, gene expression, and histology signature that differentiates between patients with different tobacco use history (heavy tobacco use with ≥10 pack years vs. no tobacco use). Using The Cancer Genome Atlas (TCGA) cohort (n = 257) and an internal cohort (n = 40), we identify 3 epigenetic markers (GPR15, GNG12, GDNF) and 13 expression markers (IGHA2, SCG5, RPL3L, NTRK1, CD96, BMP6, TFPI2, EFEMP2, RYR3, DMTN, GPD2, BAALC, and FMO3), which are dysregulated in OSCC patients who were never smokers vs. those who have a ≥ 10 pack year history. While mortality risk prediction based on smoking status and clinicopathologic covariates alone is inaccurate (c-statistic = 0.57), the combined epigenetic/expression and histologic signature has a c-statistic = 0.9409 in predicting 5-year mortality in OSCC patients.

Perineural Invasion Is Associated With Function-evoked Pain and Altered Extracellular Matrix in Patients With Head and Neck Squamous Cell Carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is painful, and perineural invasion (PNI) has been associated with the worst pain. Pain due to HNSCC is diverse and may vary based on clinicopathological factors. This study aims to characterize different pain patterns linked with PNI, its influence on daily functioning, and gain insights into molecular changes and pathways associated with PNI-related pain in HNSCC patients. We conducted a cross-sectional study across 3 medical centers (n = 114), assessing pain phenotypes and their impact on daily functioning using 2 self-reported pain questionnaires, given to patients prior to their cancer surgery. Furthermore, we conducted RNA-seq analysis utilizing the The Cancer Genome Atlas dataset of HNSCC tumor from patients (n = 192) to identify genes relevant to both PNI and pain. Upon adjusting for demographic and clinicopathological variables using linear regression models, we found that PNI independently predicted function-evoked pain according to the University of Calfornia San Francisco Oral Cancer Pain Questionnaire, as well as the worst pain intensity reported in the Brief Pain Inventory. Distinct pain patterns were observed to be associated with daily activities in varying manners. Our molecular analyses revealed significant disruptions in pathways associated with the extracellular matrix structure and organization. The top differentially expressed genes linked to the extracellular matrix are implicated in cancer development, pain, and neurodegenerative diseases. Our data underscore the importance of properly categorizing pain phenotypes in future studies aiming to uncover mechanistic underpinnings of pain. Additionally, we have compiled a list of genes of interest that could serve as targets for both cancer and cancer pain management. PERSPECTIVE: PNI independently predicts function-evoked pain. Different pain phenotypes affect daily activities differently. We identified a list of candidate genes involved in the extracellular matrix structure and function that can be targeted for both cancer and cancer pain control.

Multi-Tiered Assessment of Gene Expression Provides Evidence for Mechanisms That Underlie Risk for Type 2 Diabetes.

Introduction: Integrated transcriptome and microRNA differential gene expression (DEG) analyses may help to explain type 2 diabetes (T2D) pathogenesis in at-risk populations. The purpose of this study was to characterize DEG in banked biospecimens from underactive adult participants who responded to a randomized clinical trial measuring the effects of lifestyle interventions on T2D risk factors. DEGs were further examined within the context of annotated biological pathways. Methods: Participants (n = 52) in a previously completed clinical trial that assessed a 12-week behavioural intervention for T2D risk reduction were included. Participants who showed >6mg/dL decrease in fasting blood glucose were identified as responders. Gene expression was measured by RNASeq, and overrepresentation analysis within KEGG pathways and weighted gene correlation network analysis (WGCNA) were performed. Results: No genes remained significantly differentially expressed after correction for multiple comparisons. One module derived by WGCNA related to body mass index was identified, which contained genes located in KEGG pathways related to known mechanisms underlying risk for T2D as well as pathways related to neurodegeneration and protein misfolding. A network analysis showed indirect connections between genes in this module and islet amyloid polypeptide (IAPP), which has previously been hypothesized as a mechanism for T2D. Discussion: We validated prior studies that showed pathways related to metabolism, inflammation/immunity, and endocrine/hormone function are related to risk for T2D. We identified evidence for new potential mechanisms that include protein misfolding. Additional studies are needed to determine whether these are potential therapeutic targets to decrease risk for T2D.

Neurotrophin Pathway Receptors NGFR and TrkA Control Perineural Invasion, Metastasis, and Pain in Oral Cancer.

Oral squamous cell carcinoma (OSCC) patients suffer from poor survival due to metastasis or locoregional recurrence, processes that are both facilitated by perineural invasion (PNI). OSCC has higher rates of PNI than other cancer subtypes, with PNI present in 80% of tumors. Despite the impact of PNI on oral cancer prognosis and pain, little is known about the genes that drive PNI, which in turn drive pain, invasion, and metastasis. In this study, clinical data, preclinical, and in vitro models are leveraged to elucidate the role of neurotrophins in OSCC metastasis, PNI, and pain. The expression data in OSCC patients with metastasis, PNI, or pain demonstrate dysregulation of neurotrophin genes. TrkA and nerve growth factor receptor (NGFR) are focused, two receptors that are activated by NGF, a neurotrophin expressed at high levels in OSCC. It is demonstrated that targeted knockdown of these two receptors inhibits proliferation and invasion in an in vitro and preclinical model of OSCC, and metastasis, PNI, and pain. It is further determined that TrkA knockdown alone inhibits thermal hyperalgesia, whereas NGFR knockdown alone inhibits mechanical allodynia. Collectively the results highlight the ability of OSCC to co-opt different components of the neurotrophin pathway in metastasis, PNI, and pain.

Schwann Cells Induce Phenotypic Changes in Oral Cancer Cells.

Head and neck cancer (HNC) is the seventh most common cancer worldwide, the majority being oral squamous cell carcinoma. Despite advances in cancer diagnosis and treatment, the survival rate of patients with HNC remains stagnant. The cancer-nerve interaction has been recognized as an important driver of cancer progression. Schwann cells, a type of peripheral glia, have been implicated in promoting cancer cell growth, migration, dispersion, and invasion into the nerve in many cancers. Here, it is demonstrated that the presence of Schwann cells makes oral cancer cells more aggressive by promoting their proliferation, extracellular matrix breakdown, and altering cell metabolism. Furthermore, oral cancer cells became larger, more circular, with more projections and nuclei following co-culturing with Schwann cells. RNA-sequencing analysis in oral cancer cells following exposure to Schwann cells shows corresponding changes in genes involved in the hallmarks of cancer and cell metabolism; the enriched KEGG pathways are spliceosome, RNA transport, cell cycle, axon guidance, signaling pathways regulating pluripotency of stem cells, cAMP signaling, WNT signaling, proteoglycans in cancer and PI3K-Akt signaling. Taken together, these results suggest a significant role for Schwann cells in facilitating oral cancer progression, highlighting their potential as a target to treat oral cancer progression.

Brush swab as a noninvasive surrogate for tissue biopsies in epigenomic profiling of oral cancer.

BACKGROUND: Oral squamous cell carcinoma (OSCC) has poor survival rates. There is a pressing need to develop more precise risk assessment methods to tailor clinical treatment. Epigenome-wide association studies in OSCC have not produced a viable biomarker. These studies have relied on methylation array platforms, which are limited in their ability to profile the methylome. In this study, we use MethylCap-Seq (MC-Seq), a comprehensive methylation quantification technique, and brush swab samples, to develop a noninvasive, readily translatable approach to profile the methylome in OSCC patients. METHODS: Three OSCC patients underwent collection of cancer and contralateral normal tissue and brush swab biopsies, totaling 4 samples for each patient. Epigenome-wide DNA methylation quantification was performed using the SureSelectXT Methyl-Seq platform. DNA quality and methylation site resolution were compared between brush swab and tissue samples. Correlation and methylation value difference were determined for brush swabs vs. tissues for each respective patient and site (i.e., cancer or normal). Correlations were calculated between cancer and normal tissues and brush swab samples for each patient to determine the robustness of DNA methylation marks using brush swabs in clinical biomarker studies. RESULTS: There were no significant differences in DNA yield between tissue and brush swab samples. Mapping efficiency exceeded 90% across all samples, with no differences between tissue and brush swabs. The average number of CpG sites with at least 10x depth of coverage was 2,716,674 for brush swabs and 2,903,261 for tissues. Matched tissue and brush swabs had excellent correlation (r = 0.913 for cancer samples and r = 0.951 for normal samples). The methylation profile of the top 1000 CpGs was significantly different between cancer and normal samples (mean p-value = 0.00021) but not different between tissues and brush swabs (mean p-value = 0.11). CONCLUSIONS: Our results demonstrate that MC-Seq is an efficient platform for epigenome profiling in cancer biomarker studies, with broader methylome coverage than array-based platforms. Brush swab biopsy provides adequate DNA yield for MC-Seq, and taken together, our findings set the stage for development of a non-invasive methylome quantification technique for oral cancer with high translational potential.

The REASON score: an epigenetic and clinicopathologic score to predict risk of poor survival in patients with early stage oral squamous cell carcinoma.

BACKGROUND: Oral squamous cell carcinoma (OSCC) is a capricious cancer with poor survival rates, even for early-stage patients. There is a pressing need to develop more precise risk assessment methods to appropriately tailor clinical treatment. Genome-wide association studies have not produced a viable biomarker. However, these studies are limited by using heterogeneous cohorts, not focusing on methylation although OSCC is a heavily epigenetically-regulated cancer, and not combining molecular data with clinicopathologic data for risk prediction. In this study we focused on early-stage (I/II) OSCC and created a risk score called the REASON score, which combines clinicopathologic characteristics with a 12-gene methylation signature, to predict the risk of 5-year mortality. METHODS: We combined data from an internal cohort (n = 515) and The Cancer Genome Atlas (TCGA) cohort (n = 58). We collected clinicopathologic data from both cohorts to derive the non-molecular portion of the REASON score. We then analyzed the TCGA cohort DNA methylation data to derive the molecular portion of the risk score. RESULTS: 5-year disease specific survival was 63% for the internal cohort and 86% for the TCGA cohort. The clinicopathologic features with the highest predictive ability among the two the cohorts were age, race, sex, tobacco use, alcohol use, histologic grade, stage, perineural invasion (PNI), lymphovascular invasion (LVI), and margin status. This panel of 10 non-molecular features predicted 5-year mortality risk with a concordance (c)-index = 0.67. Our molecular panel consisted of a 12-gene methylation signature (i.e., HORMAD2, MYLK, GPR133, SOX8, TRPA1, ABCA2, HGFAC, MCPH1, WDR86, CACNA1H, RNF216, CCNJL), which had the most significant differential methylation between patients who survived vs. died by 5 years. All 12 genes have already been linked to survival in other cancers. Of the genes, only SOX8 was previously associated with OSCC; our study was the first to link the remaining 11 genes to OSCC survival. The combined molecular and non-molecular panel formed the REASON score, which predicted risk of death with a c-index = 0.915. CONCLUSIONS: The REASON score is a promising biomarker to predict risk of mortality in early-stage OSCC patients. Validation of the REASON score in a larger independent cohort is warranted.

Reduced Expression of the PP2A Methylesterase, PME-1, or the PP2A Methyltransferase, LCMT-1, Alters Sensitivity to Beta-Amyloid-Induced Cognitive and Electrophysiological Impairments in Mice.

Beta-amyloid (Aß) is thought to play a critical role in Alzheimer's disease (AD), and application of soluble oligomeric forms of Aß produces AD-like impairments in cognition and synaptic plasticity in experimental systems. We found previously that transgenic overexpression of the PP2A methylesterase, PME-1, or the PP2A methyltransferase, LCMT-1, altered the sensitivity of mice to Aß-induced impairments, suggesting that PME-1 inhibition may be an effective approach for preventing or treating these impairments. To explore this possibility, we examined the behavioral and electrophysiological effects of acutely applied synthetic Aß oligomers in male and female mice heterozygous for either a PME-1 KO or an LCMT-1 gene-trap mutation. We found that heterozygous PME-1 KO mice were resistant to Aß-induced impairments in cognition and synaptic plasticity, whereas LCMT-1 gene-trap mice showed increased sensitivity to Aß-induced impairments. The heterozygous PME-1 KO mice produced normal levels of endogenous Aß and exhibited normal electrophysiological responses to picomolar concentrations of Aß, suggesting that reduced PME-1 expression in these animals protects against Aß-induced impairments without impacting normal physiological Aß functions. Together, these data provide additional support for roles for PME-1 and LCMT-1 in regulating sensitivity to Aß-induced impairments, and suggest that inhibition of PME-1 may constitute a viable therapeutic approach for selectively protecting against the pathologic actions of Aß in AD.SIGNIFICANCE STATEMENT Elevated levels of ß-amyloid (Aß) in the brain are thought to contribute to the cognitive impairments observed in Alzheimer's disease patients. Here we show that genetically reducing endogenous levels of the PP2A methylesterase, PME-1, prevents the cognitive and electrophysiological impairments caused by acute exposure to pathologic concentrations of Aß without impairing normal physiological Aß function or endogenous Aß production. Conversely, reducing endogenous levels of the PP2A methyltransferase, LCMT-1, increases sensitivity to Aß-induced impairments. These data offer additional insights into the molecular factors that control sensitivity to Aß-induced impairments, and suggest that inhibiting PME-1 may constitute a viable therapeutic avenue for preventing Aß-related impairments in Alzheimer's disease.

Molecular patterning of the embryonic cranial mesenchyme revealed by genome-wide transcriptional profiling.

In the head of an embryo, a layer of mesenchyme surrounds the brain underneath the surface ectoderm. This cranial mesenchyme gives rise to the meninges, the calvaria (top part of the skull), and the dermis of the scalp. Abnormal development of these structures, especially the meninges and the calvaria, is linked to significant congenital defects in humans. It has been known that different areas of the cranial mesenchyme have different fates. For example, the calvarial bone develops from the cranial mesenchyme on the baso-lateral side of the head just above the eye (supraorbital mesenchyme, SOM), but not from the mesenchyme apical to SOM (early migrating mesenchyme, EMM). However, the molecular basis of this difference is not fully understood. To answer this question, we compared the transcriptomes of EMM and SOM using high-throughput sequencing (RNA-seq). This experiment identified a large number of genes that were differentially expressed in EMM and SOM, and gene ontology analyses found very different terms enriched in each region. We verified the expression of about 40 genes in the head by RNA in situ hybridization, and the expression patterns were annotated to make a map of molecular markers for 6 subdivisions of the cranial mesenchyme. Our data also provided insights into potential novel regulators of cranial mesenchyme development, including several axon guidance pathways, lectin complement pathway, cyclic-adenosine monophosphate (cAMP) signaling pathway, and ZIC family transcription factors. Together, information in this paper will serve as a unique resource to guide future research on cranial mesenchyme development.