Integrating Clinical and Genetic Analysis of Perineural Invasion in Head and Neck Squamous Cell Carcinoma

Head and Neck Osteosarcoma: Perineural Invasion is Associated With Disease-Free Survival and Tumor Metastasis

BACKGROUND

Head and neck osteosarcoma (HNOS) is the most common bone malignancy in the head and neck region, accounting for 10% of all osteosarcoma cases. Perineural invasion (PNI) is a notable indication of aggressive tumor behavior, which includes the phenomenon of tumor cells invading any of the 3 layers of the nerve sheath or tumor cells gathering, encircling one-third of the nerve circumference, and infiltrating and metastasizing along the nerve. PNI has been reported in various malignant tumors and is considered to be linked to poor prognosis.

PURPOSE

The study's purpose is to measure the association between PNI and survival outcomes in patients with HNOS.

STUDY DESIGN, SETTING, SAMPLE

This retrospective cohort study focused on HNOS patients who underwent surgery at the Department of Oral and Maxillofacial Head and Neck Oncology, Shanghai Ninth People's Hospital School of Medicine, Shanghai Jiao Tong University, from January 1, 2019 to December 31, 2021. Patients who did not undergo complete surgical resection of the tumor, did not receive a conventional osteosarcoma diagnosis, and had positive surgical margins were eliminated.

PREDICTOR VARIABLE

The predictor variable is PNI status. The pathological section of the tumor was consistent with any of the PNI features, which was considered PNI-positive.

MAIN OUTCOME VARIABLE(S)

The primary outcome variables were 3-year disease-free survival (DFS) and 3-year overall survival. Secondary outcomes were 3-year tumor local recurrence and 3-year metastasis (MT).

COVARIATES

Covariates were categorized into the following categories: demographic variables (age, sex), clinical variables (tumor region, primary tumor), and treatment variables (chemotherapy, radiotherapy).

ANALYSES

Analytic statistical methods were used for the data analysis. Pearson χ2 or Fisher's exact test was used to describe the baseline data. Kaplan-Meier is used to calculate survival rates. The Cox regression model was adapted for univariate and multivariate analysis. A P value less than .05 indicated statistical significance.

RESULTS

The study sample comprised 70 patients; 33 (47.1%) were male, and the mean age was 42.2 (standard deviation: 16.7) years. There were 15 (21.4%) cases of PNI. The 3-year DSF rate and OS rate were 67.3% and 82.0%, respectively. PNI-positive resulted in higher risk for MT (P ＜ .01, hazard ratio: 5.95, 95% confidence interval: 1.62-21.86) and negative impact on DFS (P < .01, hazard ratio: 6.35, 95% confidence interval: 2.11-19.17) for HNOS patients.

CONCLUSION AND RELEVANCE

Positive PNI status was associated with decreased DFS and increased risk of MT.

Neuronal CFL1 upregulation in head and neck squamous cell carcinoma enhances tumor-nerve crosstalk and promotes tumor growth

Head and neck squamous cell carcinoma (HNSCC) is a prevalent cancer type, marked by a pronounced nerve density within the tumor microenvironment and a high rate of perineural invasion (PNI). Growing evidence suggests that the nervous system plays a vital role in HNSCC progression. Yet, the mechanisms governing cancer-nerve interactions remain largely elusive. Our research revealed that cofilin-1 (CFL1) is significantly overexpressed in HNSCC and correlates with both PNI and unfavorable prognosis. Utilizing multiplex fluorescent immunohistochemistry, we have localized CFL1 chiefly to the nerves adjacent to tumor sites. Significantly, it is the elevated expression of CFL1 in neuronal structures, rather than in the tumor cells, that aligns with diminished patient survival rates. We observed that HNSCC cells induced the expression of neuronal CFL1 and that the conditional knockout of neuronal CFL1 impedes tumor-nerve interactions. Both Gene Ontology functional enrichment analyses and Gene Set Enrichment Analysis demonstrate that CFL1 expression in HNSCC is associated with specific biological processes, including "RIBOSOME," "PROTEASOME," and "cadherin binding." In summary, HNSCC promotes the expression of CFL1 in nerves, which is essential for cancer-nerve interactions. The neuronal CFL1 is associated with PNI and may be a potential molecular prognostic marker of poor survival in HNSCC.

Sensory nerve release of CGRP increases tumor growth in HNSCC by suppressing TILs

BACKGROUND

Perineural invasion (PNI) and nerve density within the tumor microenvironment (TME) have long been associated with worse outcomes in head and neck squamous cell carcinoma (HNSCC). This prompted an investigation into how nerves within the tumor microenvironment affect the adaptive immune system and tumor growth.

METHODS

We used RNA sequencing analysis of human tumor tissue from a recent HNSCC clinical trial, proteomics of human nerves from HNSCC patients, and syngeneic orthotopic murine models of HPV-unrelated HNSCC to investigate how sensory nerves modulate the adaptive immune system.

FINDINGS

Calcitonin gene-related peptide (CGRP) directly inhibited CD8 T cell activity in vitro, and blocking sensory nerve function surgically, pharmacologically, or genetically increased CD8 and CD4 T cell activity in vivo.

CONCLUSIONS

Our data support sensory nerves playing a role in accelerating tumor growth by directly acting on the adaptive immune system to decrease Th1 CD4 T cells and activated CD8 T cells in the TME. These data support further investigation into the role of sensory nerves in the TME of HNSCC and points toward the possible treatment efficacy of blocking sensory nerve function or specifically inhibiting CGRP release or activity within the TME to improve outcomes.

FUNDING

1R01DE028282-01, 1R01DE028529-01, 1P50CA261605-01 (to S.D.K.), 1R01CA284651-01 (to S.D.K.), and F31 DE029997 (to L.B.D.).

Tumor Neurobiology in the Pathogenesis and Therapy of Head and Neck Cancer

The neurobiology of tumors has attracted considerable interest from clinicians and scientists and has become a multidisciplinary area of research. Neural components not only interact with tumor cells but also influence other elements within the TME, such as immune cells and vascular components, forming a polygonal relationship to synergistically facilitate tumor growth and progression. This review comprehensively summarizes the current state of the knowledge on nerve-tumor crosstalk in head and neck cancer and discusses the potential underlying mechanisms. Several mechanisms facilitating nerve-tumor crosstalk are covered, such as perineural invasion, axonogenesis, neurogenesis, neural reprogramming, and transdifferentiation, and the reciprocal interactions between the nervous and immune systems in the TME are also discussed in this review. Further understanding of the nerve-tumor crosstalk in the TME of head and neck cancer may provide new nerve-targeted treatment options and help improve clinical outcomes for patients.

Prognostic DNA mutation and mRNA expression analysis of perineural invasion in oral squamous cell carcinoma

This study analyzed oral squamous cell carcinoma (OSCC) genomes and transcriptomes in relation to perineural invasion (PNI) and prognosis using Cancer Genome Atlas data and validated these results with GSE41613 data. Gene set enrichment analysis (GSEA), gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes were conducted. We identified 22 DNA mutations associated with both overall survival (OS) and PNI. Among them, TGFBR1 and RPS6KA4 mRNAs were overexpressed, while TYRO3 and GPR137 mRNAs were underexpressed in PNI patients. Among the 141 mRNA genes associated with both OS and PNI, we found overlap with PNI-related DNA mutations, including ZNF43, TEX10, TPSD1, and PSD3. In GSE41613 data, TGFBR1, RPS6KA4, TYRO3, GPR137, TEX10 and TPSD1 mRNAs were expressed differently according to the prognosis. The 22 DNA-mutated genes clustered into nervous system development, regulation of DNA-templated transcription, and transforming growth factor beta binding. GSEA analysis of mRNAs revealed upregulation of hallmarks epithelial mesenchymal transition (EMT), TNFα signaling via NF-κB, and IL2 STAT5 signaling. EMT upregulation aligned with the TGFBR1 DNA mutation, supporting its significance in PNI. These findings suggest a potential role of PNI genes in the prognosis of OSCC, providing insights for diagnosis and treatment of OSCC.

Association analysis of polymorphisms at GLRB, GRIA2, and GASK1B genes with reproductive traits in Dazu Black Goats

Reproductive traits are essential economic traits in goats. This study aimed to analyze the relationship between single nucleotide polymorphisms (SNPs) within the genes of GLRB, GRIA2, and GASK1B, and reproductive traits (kidding traits and placental traits) in goats. We used the resequencing data of 150 Dazu Black Goats to perform correlation analysis with the average litter size. We screened thirteen SNPs loci in introns and then used the Sanger method to genotype the remaining 150 Dazu Black Goats. The results showed that a total of six SNPs were screened. Three SNPs related to litter size and live litter size (g.28985790T > G, g.28986352A > G, and g.28987976A > G); one SNP related to total cotyledon area (g.29203243G > A); two SNPs related to placental efficiency (g.30189055G > A and g.30193974C > T); one SNP associated with cotyledon support efficiency (g.30193974C > T). The qPCR results showed that GLRB, GRIA2, and GASK1B were all highly expressed in the udder, kidney, uterus, and ovary. It indicated that these three candidate genes might affect the reproductive traits, which could be used as candidate markers for reproductive traits in Dazu Black Goats. Moreover, association studies on a large scale are still needed to figure out what effect these SNPs have on reproductive traits.

DKK3 expression is correlated with poorer prognosis in head and neck squamous cell carcinoma: A bioinformatics study based on the TCGA database

OBJECTIVE

We previously reported that dickkopf WNT signaling pathway inhibitor 3 (DKK3) expression is correlated with poorer prognosis in head and neck squamous cell carcinoma (HNSCC). Here we investigated DKK3 expression by using The Cancer Genome Atlas (TCGA) public database and bioinformatic analyses.

METHODS

We used the RNA sequence data and divided the tumor samples into "DKK3-high" and "DKK3-low" groups according to median DKK3 expression. The correlations between DKK3 expression and the clinical data were investigated. Differentially expressed genes (DEGs) were detected using DESEq2 and analyzed by ShinyGO 0.77. A gene set enrichment analysis (GSEA) was also performed using GSEA software. The DEGs were also analyzed with TargetMine to establish the protein-protein interaction (PPI) network.

RESULTS

DKK3 expression was significantly increased in cancer samples, and a high DKK3 expression was significantly associated with shorter overall survival. We identified 854 DEGs, including 284 up-regulated and 570 down-regulated. Functional enrichment analyses revealed several Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with extracellular matrix remodeling. The PPI network identified COL8A1, AGTR1, FN1, P4HA3, PDGFRB, and CEP126 as the key genes.

CONCLUSIONS

These results suggested the cancer-promoting ability of DKK3, the expression of which is a promising prognostic marker and therapeutic target for HNSCC.

Overexpression of POFUT1 promotes malignant phenotype and mediates perineural invasion in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is one of the most aggressive neoplasms, which requires more effective prevention and treatment modalities. Previous studies found that protein O-fucosyltransferase 1 (POFUT1) upregulation promotes carcinogenesis, although the potential roles, underlying molecular mechanisms, and biological implications of POFUT1 in HNSCC were not investigated. In this study, in silico analyses referred POFUT1 as a potential oncogene in HNSCC. Further analysis of tumor and normal tissue samples as well as HNSCC cells with quantitative real-time polymerase chain reaction, Western blot analysis, and immunohistochemistry showed significant overexpression of POFUT1 in HNSCC clinical tumor tissue specimens and cell lines compared to corresponding controls. In vitro investigations revealed that overexpression of POFUT1 promoted phenotypes associated with cancer aggressiveness and its knockdown in HNSCC cells suppressed those phenotypes. Further xenograft experiments demonstrated that POFUT1 is an oncogene in vivo for HNSCC. Immunohistochemical analysis with human clinical samples and cancer cell-dorsal root ganglion ex-vivo coculture model showed that deregulation of POFUT1 is involved in the perineural invasion of HNSCC cells. These results suggest POFUT1 expression as a potential prognostic marker for patients with head and neck cancer and highlight its potential as a target for HNSCC therapy, although more molecular clues are needed to better define the functions of POFUT1 related to HNSCC carcinogenesis.

New perspectives on biology, disease progression, and therapy response of head and neck cancer gained from single cell RNA sequencing and spatial transcriptomics

Head and neck squamous cell carcinoma (HNSCC) is one of the most frequent cancers worldwide. The main risk factors are consumption of tobacco products and alcohol, as well as infection with human papilloma virus. Approved therapeutic options comprise surgery, radiation, chemotherapy, targeted therapy through epidermal growth factor receptor inhibition, and immunotherapy, but outcome has remained unsatisfactory due to recurrence rates of ~50% and the frequent occurrence of second primaries. The availability of the human genome sequence at the beginning of the millennium heralded the omics era, in which rapid technological progress has advanced our knowledge of the molecular biology of malignant diseases, including HNSCC, at an unprecedented pace. Initially, microarray-based methods, followed by approaches based on next-generation sequencing, were applied to study the genetics, epigenetics, and gene expression patterns of bulk tumors. More recently, the advent of single-cell RNA sequencing (scRNAseq) and spatial transcriptomics methods has facilitated the investigation of the heterogeneity between and within different cell populations in the tumor microenvironment (e.g., cancer cells, fibroblasts, immune cells, endothelial cells), led to the discovery of novel cell types, and advanced the discovery of cell-cell communication within tumors. This review provides an overview of scRNAseq, spatial transcriptomics, and the associated bioinformatics methods, and summarizes how their application has promoted our understanding of the emergence, composition, progression, and therapy responsiveness of, and intercellular signaling within, HNSCC.

Pretreatment pain predicts perineural invasion in patients with head and neck squamous cell carcinoma

OBJECTIVES

Perineural invasion (PNI) in head and neck cancer (HNC) is a distinct pathological feature used to indicate aggressive tumor behavior and drive treatment strategies. Our study examined the prevalence and predictors of PNI in HNC patients stratified by tumor site.

STUDY DESIGN AND METHODS

A retrospective analysis of head and neck squamous cell carcinoma (HNSCC) patients who underwent surgical resection at the University of Pittsburgh Medical Center between 2015 and 2018 was performed. Pretreatment pain was assessed at least 1 week before surgery using the Functional Assessment of Cancer Therapy-Head and Neck (FACT-H&N). Demographics, clinical characteristics, and concomitant medications were obtained from medical records. Patients with cancers at the oropharynx and non-oropharynx (i.e., cancer at oral cavity, mandible, larynx) sites were separately analyzed. Tumor blocks were obtained from 10 patients for histological evaluation of intertumoral nerve presence.

RESULTS

A total of 292 patients (202 males, median age = 60.94 ± 11.06) were assessed. Pain and PNI were significantly associated with higher T stage (p < 0.001) and tumor site (p < 0.001); patients with non-oropharynx tumors reported more pain and had a higher incidence of PNI compared to oropharynx tumors. However, multivariable analysis identified pain as a significant variable uniquely associated with PNI for both tumor sites. Evaluation of nerve presence in tumor tissue showed 5-fold higher nerve density in T2 oral cavity tumors compared to oropharyngeal tumors.

CONCLUSIONS

Our study finds that PNI is associated with pretreatment pain and tumor stage. These data support the need for additional research into the impact of tumor location when investigating targeted therapies of tumor regression.

Establishment of a Machine Learning Model for the Risk Assessment of Perineural Invasion in Head and Neck Squamous Cell Carcinoma

Perineural invasion is a prevalent pathological finding in head and neck squamous cell carcinoma and a risk factor for unfavorable survival. An adequate diagnosis of perineural invasion by pathologic examination is limited due to the availability of tumor samples from surgical resection, which can arise in cases of definitive nonsurgical treatment. To address this medical need, we established a random forest prediction model for the risk assessment of perineural invasion, including occult perineural invasion, and characterized distinct cellular and molecular features based on our new and extended classification. RNA sequencing data of head and neck squamous cell carcinoma from The Cancer Genome Atlas were used as a training cohort to identify differentially expressed genes that are associated with perineural invasion. A random forest classification model was established based on these differentially expressed genes and was validated by inspection of H&E-stained whole image slides. Differences in epigenetic regulation and the mutational landscape were detected by an integrative analysis of multiomics data and single-cell RNA-sequencing data were analyzed. We identified a 44-gene expression signature related to perineural invasion and enriched for genes mainly expressed in cancer cells according to single-cell RNA-sequencing data. A machine learning model was trained based on the expression pattern of the 44-gene set with the unique feature to predict occult perineural invasion. This extended classification model enabled a more accurate analysis of alterations in the mutational landscape and epigenetic regulation by DNA methylation as well as quantitative and qualitative differences in the cellular composition in the tumor microenvironment between head and neck squamous cell carcinoma with or without perineural invasion. In conclusion, the newly established model could not only complement histopathologic examination as an additional diagnostic tool but also guide the identification of new drug targets for therapeutic intervention in future clinical trials with head and neck squamous cell carcinoma patients at a higher risk for treatment failure due to perineural invasion.

LINC00173 blocks GATA6-mediated transcription of COL5A1 to affect malignant development of oral squamous cell carcinoma

BACKGROUND

Aberrant expression of collagen type V alpha 1 chain (COL5A1) has been linked to several forms of human cancers. In this work, we focused on the interaction of the LINC00173/GATA binding protein 6 (GATA6)/COL5A1 axis in the malignant property of oral squamous cell carcinoma (OSCC) cells.

METHODS

We analyzed six publicly accessible datasets GSE160042, GSE74530, GSE138206, GSE23558, GSE31853 and GSE146483 to identify aberrantly expressed genes in OSCC. The expression of COL5A1 in OSCC tissues and cell lines was examined by reverse transcription-quantitative polymerase chain reaction and/or immunohistochemistry. The regulatory mechanism responsible for COL5A1 transcription was predicted via bioinformatics systems, and the interactions of LINC00173, GATA6, and COL5A1 were identified by immunoprecipitation and luciferase assays. Overexpression or downregulation of COL5A1, GATA6, and LINC00173 were induced in OSCC cell lines to determine their roles in the malignant phenotype of the OSCC cells in vitro and in vivo.

RESULTS

COL5A1 showed elevated expression in OSCC tissues and cells. The COLA51 knockdown suppressed proliferation, migration and invasiveness, apoptosis resistance, and pro-angiogenic ability of OSCC cells, and it suppressed the growth and dissemination of xenograft tumors in vivo. GATA6 bound to COL5A1 promoter to activate its transcription, whereas LINC00173 bound to GATA6 to block this transcriptional activation. Overexpression of GATA6 or COL5A1 promoted the malignant phenotype of the OSCC cells, which were blocked upon LINC00173 upregulation.

CONCLUSION

This work demonstrates that LINC00173 blocks GATA6-mediated transcription of COL5A1 to affect malignant development of OSCC.

Important Cells and Factors from Tumor Microenvironment Participated in Perineural Invasion

Perineural invasion (PNI) as the fourth way for solid tumors metastasis and invasion has attracted a lot of attention, recent research reported a new point that PNI starts to include axon growth and possible nerve "invasion" to tumors as the component. More and more tumor-nerve crosstalk has been explored to explain the internal mechanism for tumor microenvironment (TME) of some types of tumors tends to observe nerve infiltration. As is well known, the interaction of tumor cells, peripheral blood vessels, extracellular matrix, other non-malignant cells, and signal molecules in TME plays a key role in the occurrence, development, and metastasis of cancer, as to the occurrence and development of PNI. We aim to summarize the current theories on the molecular mediators and pathogenesis of PNI, add the latest scientific research progress, and explore the use of single-cell spatial transcriptomics in this invasion way. A better understanding of PNI may help to understand tumor metastasis and recurrence and will be beneficial for improving staging strategies, new treatment methods, and even paradigm shifts in our treatment of patients.

Pro-metastatic and mesenchymal gene expression signatures characterize circulating tumor cells of neuroblastoma patients with bone marrow metastases and relapse

Existing marker-based methods of minimal residual disease (MRD) determination in neuroblastoma do not effectively enrich for the circulating disease cell population. Given the relative size differential of neuroblastoma tumor cells over normal hematogenous cells, we hypothesized that cell size-based separation could enrich circulating tumor cells (CTCs) from blood samples and disseminated tumor cells (DTCs) from bone marrow aspirates (BMA) of neuroblastoma patients, and that their gene expression profiles could vary dynamically with various disease states over the course of treatment. Using a spiral microfluidic chip, peripheral blood of 17 neuroblastoma patients at 3 serial treatment timepoints (diagnosis, n=17; post-chemotherapy, n=11; and relapse, n=3), and bone marrow samples at diagnosis were enriched for large intact circulating cells. Profiling the resulting enriched samples with immunohistochemistry and mRNA expression of 1490 cancer-related genes via NanoString, 13 of 17 samples contained CTCs displaying cytologic atypia, TH and PHOX2B expression and/or upregulation of cancer-associated genes. Gene signatures reflecting pro-metastatic processes and the neuroblastoma mesenchymal super-enhancer state were consistently upregulated in 7 of 13 samples, 6 of which also had metastatic high-risk disease. Expression of 8 genes associated with PI3K and GCPR signaling were significantly upregulated in CTCs of patients with bone marrow metastases versus patients without. Correspondingly, in patients with marrow metastases, differentially-expressed gene signatures reflected upregulation of immune regulation in bone marrow DTCs versus paired CTCs samples. In patients who later developed disease relapse, 5 genes involved in immune cell regulation, JAK/STAT signaling and the neuroblastoma mesenchymal super-enhancer state (OLFML2B, STAT1, ARHGDIB, STAB1, TLR2) were upregulated in serial CTC samples over their disease course, despite urinary catecholamines and bone marrow aspirates not indicating the disease recurrences. In summary, using a label-free cell size-based separation method, we enriched and characterized intact circulating cells in peripheral blood indicative of neuroblastoma CTCs, as well as their DTC counterparts in the bone marrow. Expression profiles of pro-metastatic genes in CTCs correlated with the presence of bone marrow metastases at diagnosis, while longitudinal profiling identified persistently elevated expression of genes in CTCs that may serve as novel predictive markers of hematogenous MRD in neuroblastoma patients that subsequently relapse.

Spatial and Transcriptomic Analysis of Perineural Invasion in Oral Cancer

PURPOSE

Perineural invasion (PNI), a common occurrence in oral squamous cell carcinomas, is associated with poor survival. Consequently, these tumors are treated aggressively. However, diagnostic criteria of PNI vary and its role as an independent predictor of prognosis has not been established. To address these knowledge gaps, we investigated spatial and transcriptomic profiles of PNI-positive and PNI-negative nerves.

EXPERIMENTAL DESIGN

Tissue sections from 142 patients were stained with S100 and cytokeratin antibodies. Nerves were identified in two distinct areas: tumor bulk and margin. Nerve diameter and nerve-to-tumor distance were assessed; survival analyses were performed. Spatial transcriptomic analysis of nerves at varying distances from tumor was performed with NanoString GeoMx Digital Spatial Profiler Transcriptomic Atlas.

RESULTS

PNI is an independent predictor of poor prognosis among patients with metastasis-free lymph nodes. Patients with close nerve-tumor distance have poor outcomes even if diagnosed as PNI negative using current criteria. Patients with large nerve(s) in the tumor bulk survive poorly, suggesting that even PNI-negative nerves facilitate tumor progression. Diagnostic criteria were supported by spatial transcriptomic analyses of >18,000 genes; nerves in proximity to cancer exhibit stress and growth response changes that diminish with increasing nerve-tumor distance. These findings were validated in vitro and in human tissue.

CONCLUSIONS

This is the first study in human cancer with high-throughput gene expression analysis in nerves with striking correlations between transcriptomic profile and clinical outcomes. Our work illuminates nerve-cancer interactions suggesting that cancer-induced injury modulates neuritogenesis, and supports reclassification of PNI based on nerve-tumor distance rather than current subjective criteria.

Brain Neural Progenitors are New Predictive Biomarkers for Breast Cancer Hormonotherapy

Heterogeneity of the tumor microenvironment (TME) is one of the major causes of treatment resistance in breast cancer. Among TME components, nervous system role in clinical outcome has been underestimated. Identifying neuronal signatures associated with treatment response will help to characterize neuronal influence on tumor progression and identify new treatment targets. The search for hormonotherapy-predictive biomarkers was implemented by supervised machine learning (ML) analysis on merged transcriptomics datasets from public databases. ML-derived genes were investigated by pathway enrichment analysis, and potential gene signatures were curated by removing the variables that were not strictly nervous system specific. The predictive and prognostic abilities of the generated signatures were examined by Cox models, in the initial cohort and seven external cohorts. Generated signature performances were compared with 14 other published signatures, in both the initial and external cohorts. Underlying biological mechanisms were explored using deconvolution tools (CIBERSORTx and xCell). Our pipeline generated two nervous system-related signatures of 24 genes and 97 genes (NervSign24 and NervSign97). These signatures were prognostic and hormonotherapy-predictive, but not chemotherapy-predictive. When comparing their predictive performance with 14 published risk signatures in six hormonotherapy-treated cohorts, NervSign97 and NervSign24 were the two best performers. Pathway enrichment score and deconvolution analysis identified brain neural progenitor presence and perineural invasion as nervous system-related mechanisms positively associated with NervSign97 and poor clinical prognosis in hormonotherapy-treated patients. Transcriptomic profiling has identified two nervous system-related signatures that were validated in clinical samples as hormonotherapy-predictive signatures, meriting further exploration of neuronal component involvement in tumor progression.

Significance

The development of personalized and precision medicine is the future of cancer therapy. With only two gene expression signatures approved by FDA for breast cancer, we are in need of new ones that can reliably stratify patients for optimal treatment. This study provides two hormonotherapy-predictive and prognostic signatures that are related to nervous system in TME. It highlights tumor neuronal components as potential new targets for breast cancer therapy.

Pan-Cancer Analysis of OLFML2B Expression and Its Association With Prognosis and Immune Infiltration

Background: The function of olfactomedin-like 2B (OLFML2B), as a member of the olfactomedin domain-containing protein family, remains ambiguous, especially in tumors. The current study explores the possible correlation between OLFML2B, prognosis, and immune infiltration in pan-cancer. Methods: We applied a number of bioinformatics techniques to probe the prospective function of OLFML2B, consisting of its association with prognosis, clinicopathology, alteration, GSEA, tumor microenvironment (TME), immune-associated genes, immune infiltration, tumor mutational burden (TMB), microsatellite instability (MSI), and drug sensitivity in several cancer types. qPCR and immunohistochemistry were used to identify OLFML2B expression in LIHC cell lines and liver cancer tissues. Results: We discovered that OLFML2B was overexpressed in 14 cancers and positively related to several cancer type prognoses. The expression of OLFML2B was further validated in the LIHC cell lines. OLFML2B expression was bound up with TMB in 13 cancers, MSI in 10 cancers, and TME in almost all cancers. Furthermore, OLFML2B was highly co-expressed with genes encoding immune activators and immune suppressors. We further found that OLFML2B played a role in infiltrating different types of immune cells, such as macrophages and cancer-associated fibroblasts. OLFML2B may influence various cancer and immune-related pathways, such as the PI3K-Akt signaling pathway, ECM-receptor interaction, focal adhesion, and leukocyte transendothelial migration. In addition, OLFML2B may increase drug resistance of binimetinib, cobimentinib, and trametinib. Conclusion: Our outcomes reveal that OLFML2B may act as a prognostic marker and a potential target in immunotherapy for diverse tumors due to its oncogenesis function and immune infiltration.

The role of matrixmetalloproteinase-2 expression by fibroblasts in perineural invasion by oral cavity squamous cell carcinoma

OBJECTIVES

The mechanisms of perineural invasion (PNI) in oral cavity squamous cell carcinoma are only partially understood, and no studies have specifically investigated the role of perineural fibroblasts in PNI. Here, we identified fibroblasts within the microenvironment of perineural invasion and assessed their expression of matrix metalloproteinase-2 (MMP-2).

MATERIALS AND METHODS

Tumor specimens from 12 patients with oral cavity squamous cell carcinoma and pathologically-confirmed perineural invasion were stained by immunohistochemistry (IHC) for vimentin (positive control) and MMP-2. Scoring was quantified and compared at nerves involved with PNI and nerves uninvolved with PNI.

RESULTS

All 12 patients had perineural fibroblasts around involved and uninvolved nerves as marked by vimentin IHC staining. Perineural fibroblasts had detectable MMP-2 expression at areas of perineural invasion in all 12 patients, but no patients had MMP-2 expression by fibroblasts at nerves without PNI.

CONCLUSION

MMP-2 is expressed by fibroblasts within the microenvironment of perineural invasion, and MMP-2 expression by fibroblasts is a possible mechanism of perineural invasion by oral cavity squamous cell carcinoma. MMP-2 may be an anti-cancer target among oral cavity squamous cell carcinoma patients with PNI.

Development and external validation of a novel nomogram for predicting cancer-specific survival in patients with ascending colon adenocarcinoma after surgery: a population-based study

Background

This study aimed to develop and validate a novel nomogram to predict the cancer-specific survival (CSS) of patients with ascending colon adenocarcinoma after surgery.

Methods

Patients with ascending colon adenocarcinoma were enrolled from the Surveillance, Epidemiology, and End Results (SEER) database from 1973 to 2015 and randomly divided into a training set (5930) and a validation set (2540). The cut-off values for age, tumour size and lymph node ratio (LNR) were calculated via X-tile software. In the training set, independent prognostic factors were identified using univariate and multivariate Cox analyses, and a nomogram incorporating these factors was subsequently built. Data from the validation set were used to assess the reliability and accuracy of the nomogram and then compared with the 8th edition of the American Joint Committee on Cancer (AJCC) tumour-node-metastasis (TNM) staging system. Furthermore, external validation was performed from a single institution in China.

Results

A total of 8470 patients were enrolled from the SEER database, 5930 patients were allocated to the training set, 2540 were allocated to the internal validation set and a separate set of 473 patients was allocated to the external validation set. The optimal cut-off values of age, tumour size and lymph node ratio were 73 and 85, 33 and 75 and 4.9 and 32.8, respectively. Univariate and multivariate Cox multivariate regression revealed that age, AJCC 8th edition T, N and M stage, carcinoembryonic antigen (CEA), tumour differentiation, chemotherapy, perineural invasion and LNR were independent risk factors for patient CSS. The nomogram showed good predictive ability, as indicated by discriminative ability and calibration, with C statistics of 0.835 (95% CI, 0.823–0.847) and 0.848 (95% CI, 0.830–0.866) in the training and validation sets and 0.732 (95% CI, 0.664–0.799) in the external validation set. The nomogram showed favourable discrimination and calibration abilities and performed better than the AJCC TNM staging system.

Conclusions

A novel validated nomogram could effectively predict patients with ascending colon adenocarcinoma after surgery, and this predictive power may guide clinicians in accurate prognostic judgement.

Advances in Head and Neck Cancer Pain

Head and neck cancer (HNC) affects over 890,000 people annually worldwide and has a mortality rate of 50%. Aside from poor survival, HNC pain impairs eating, drinking, and talking in patients, severely reducing quality of life. Different pain phenotype in patients (allodynia, hyperalgesia, and spontaneous pain) results from a combination of anatomical, histopathological, and molecular differences between cancers. Poor pathologic features (e.g., perineural invasion, lymph node metastasis) are associated with increased pain. The use of syngeneic/immunocompetent animal models, as well as a new mouse model of perineural invasion, provides novel insights into the pathobiology of HNC pain. Glial and immune modulation of the tumor microenvironment affect not only cancer progression but also pain signaling. For example, Schwann cells promote cancer cell proliferation, migration, and secretion of nociceptive mediators, whereas neutrophils are implicated in sex differences in pain in animal models of HNC. Emerging evidence supports the existence of a functional loop of cross-activation between the tumor microenvironment and peripheral nerves, mediated by a molecular exchange of bioactive contents (pronociceptive and protumorigenic) via paracrine and autocrine signaling. Brain-derived neurotrophic factor, tumor necrosis factor α, legumain, cathepsin S, and A disintegrin and metalloprotease 17 expressed in the HNC microenvironment have recently been shown to promote HNC pain, further highlighting the importance of proinflammatory cytokines, neurotrophic factors, and proteases in mediating HNC-associated pain. Pronociceptive mediators, together with nerve injury, cause nociceptor hypersensitivity. Oncogenic, pronociceptive mediators packaged in cancer cell-derived exosomes also induce nociception in mice. In addition to increased production of pronociceptive mediators, HNC is accompanied by a dampened endogenous antinociception system (e.g., downregulation of resolvins and µ-opioid receptor expression). Resolvin treatment or gene delivery of µ-opioid receptors provides pain relief in preclinical HNC models. Collectively, recent studies suggest that pain and HNC progression share converging mechanisms that can be targeted for cancer treatment and pain management.

Phosphorylation of MICAL2 by ARG promotes head and neck cancer tumorigenesis by regulating skeletal rearrangement

The actin cytoskeletal architecture provides the structural underpinnings for crucial cellular behaviors. In cancer cells, changes in the actin cytoskeleton may serve as prerequisites for proliferation, invasion, and metastatic dissemination. However, the underlying mechanisms remain largely unknown. Here, we show that MICAL2, which is increased in head and neck squamous cell carcinoma (HNSCC) and inversely associated with patient survival, promotes HNSCC growth, invasion, and migration. MICAL2 serves as a flavoprotein monooxygenase and directly induces actin filament depolymerization by specifically oxidizing the methionine 44 and 47 residues of F-actin. The kinase ARG interacts with MICAL2 and augments MICAL2-mediated actin disassembly. Direct phosphorylation assay and mass spectrometry confirmed that ARG phosphorylates MICAL2 at Tyr445, Tyr463, and Tyr488. Substitution of the Tyr445 or Tyr463 residue of purified recombinant MICAL2-redox with phenylalanine (generating a non-phosphorylatable mutant) abolishes the enhanced MICAL2-mediated F-actin disassembly induced by ARG. Consistently, ectopic expression of non-phosphorylatable MICAL2 mutants (MICAL2^Y445F and MICAL2^Y463F, not MICAL2^Y488F) failed to ameliorate HNSCC cell growth, whereas expression of wild-type MICAL2 or MICAL2^Y488F rescued the impaired proliferation induced by MICAL2 knockdown. Moreover, CCG-1423, an inhibitor of MICAL2, was shown to inhibit HNSCC cell proliferation, invasion, and migration. Taken together, our findings indicate that phosphorylation of MICAL2 at Tyr445 and Tyr463 by ARG mediates F-actin disassembly and promotes HNSCC progression.

Preoperative Evaluation of Perineural Invasion in Cervical Cancer: Development and Independent Validation of a Novel Predictive Nomogram

BACKGROUND

Perineural invasion (PNI) is associated with a poor prognosis for cervical cancer and influences surgical strategies. However, a preoperative evaluation that can determine PNI in cervical cancer patients is lacking.

METHODS

After 1:1 propensity score matching, 162 cervical cancer patients with PNI and 162 cervical cancer patients without PNI were included in the training set. Forty-nine eligible patients were enrolled in the validation set. The PNI-positive and PNI-negative groups were compared. Multivariate logistic regression was performed to build the PNI prediction nomogram.

RESULTS

Age [odds ratio (OR), 1.028; 95% confidence interval (CI), 0.999-1.058], adenocarcinoma (OR, 1.169; 95% CI, 0.675-2.028), tumor size (OR, 1.216; 95% CI, 0.927-1.607), neoadjuvant chemotherapy (OR, 0.544; 95% CI, 0.269-1.083), lymph node enlargement (OR, 1.953; 95% CI, 1.086-3.550), deep stromal invasion (OR, 1.639; 95% CI, 0.977-2.742), and full-layer invasion (OR, 5.119; 95% CI, 2.788-9.799) were integrated in the PNI prediction nomogram based on multivariate logistic regression. The PNI prediction nomogram exhibited satisfactory performance, with areas under the curve of 0.763 (95% CI, 0.712-0.815) for the training set and 0.860 (95% CI, 0.758-0.961) for the validation set. Moreover, after reviewing the pathological slides of patients in the validation set, four patients initially diagnosed as PNI-negative were recognized as PNI-positive. All these four patients with false-negative PNI were correctly predicted to be PNI-positive (predicted p > 0.5) by the nomogram, which improved the PNI detection rate.

CONCLUSION

The nomogram has potential to assist clinicians when evaluating the PNI status, reduce misdiagnosis, and optimize surgical strategies for patients with cervical cancer.

OLFML2B Is a Robust Prognostic Biomarker in Bladder Cancer Through Genome-Wide Screening: A Study Based on Seven Cohorts

Background

Bladder cancer lacks useful and robust prognostic markers to stratify patients at risk. Our study is to identify a robust prognostic marker for bladder cancer.

Methods

The transcriptome and clinical data of bladder cancer were downloaded from multiple databases. We searched for genes with robust prognosis by Kaplan-Meier analysis of the whole genome. CIBERSORT and TIMER algorithm was used to calculate the degree of immune cell infiltration.

Results

We identified OLFML2B as a robust prognostic marker for bladder cancer in five cohorts. Kaplan-Meier analysis showed that patients with a high level of OLFML2B expression had a poor prognosis. The expression of OLFML2B increased with the increase of stage and grade. We found that patients with high expression of OLFML2B still had a poor prognosis in two small bladder cancer cohorts. OLFML2B also has the prognostic ability in ten other tumors, and the prognosis is poor in high expression. The correlation analysis between OLFML2B and immune cells showed that it was positively correlated with the degree of macrophage infiltration and highly co-expressed with tumor-associated macrophage markers. Finally, the Wound-healing assay and Colony formation assay results showed that the migration and proliferation ability of bladder cancer cell lines decreased after the knockdown of OLFML2B.

Conclusions

In summary, OLFML2B is a robust risk prognostic marker, and it can help patients with bladder cancer improve individualized treatment.

BCAR3 promotes head and neck cancer growth and is associated with poor prognosis

Breast cancer anti-estrogen resistance protein 3 (BCAR3) is involved in anti-estrogen resistance and other important aspects of breast cancer. However, the role of BCAR3 in other solid tumors remains unclear. The relationship between the clinicopathologic characteristics of head and neck squamous cell carcinoma (HNSCC) patients and BCAR3 was analyzed using the Wilcoxon’s signed-rank test and logistic regression. The association between BCAR3 expression and clinicopathologic features and survival was analyzed using Cox regression and the Kaplan–Meier method. In vivo and in vitro assays were performed to validate the effect of BCAR3 on HNSCC growth. BCAR3-related mRNAs were determined by calculating the Pearson’s correlation coefficient based on The Cancer Genome Atlas (TCGA). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses, and gene set enrichment analysis (GSEA) were used to predict the potential functions of BCAR3. BCAR3 expression is overexpressed in HNSCC and was shown to be associated with perineural invasion (PNI) and poor survival. BCAR3 silencing significantly attenuated the proliferation of HNSCC cells, whereas BCAR3 depletion inhibited tumor growth in vitro. GO and KEGG functional enrichment analyses, and GSEA showed that BCAR3 expression in HNSCC was associated with biological processes, such as cell adhesion, actin binding, cadherin binding, and angiogenesis. BCAR3, which promotes HNSCC growth, is associated with perineural invasion and may be a potential molecular prognostic marker of poor survival in HNSCC.

Systematic review on the current knowledge and use of single-cell RNA sequencing in head and neck cancer

Single-cell RNA sequencing (scRNA-seq) is a novel method enabling genetic characterization of tumor tissue at a single-cell level. This study systematically reviewed the literature on studies using scRNA-seq to characterize head and neck squamous cell carcinomas (HNSCCs). Seven studies were included, of which two studies performed scRNA-seq on 20 patients in total, and five studies used scRNA-seq data in a subsequent clinical study. The former mentioned two studies found intra-tumoral genetic heterogeneity among malignant cells but genetic uniformity among non-malignant cells. The five latter studies used scRNA-seq data in various ways. Three studies identified biomarkers related to predicting whether a patient would benefit from immunotherapeutic treatment. One study characterized genes related to the perineural invasion. One study identified genes to be used in diagnostics. Further studies performing scRNA-seq on HNSCC are required to continue the ongoing development and use of scRNA-seq.

Identification of Key Genes and Pathways associated with Endometriosis by Weighted Gene Co-expression Network Analysis

Background: Endometriosis is a common gynecological disorder with high rates of infertility and pelvic pain. However, its pathogenesis and diagnostic biomarkers remain unclear. This study aimed to elucidate potential hub genes and key pathways associated with endometriosis in ectopic endometrium (EC) and eutopic endometrium (EU). Material and Method: EC and EU-associated microarray datasets were obtained from the gene expression omnibus (GEO) database. Gene set enrichment analysis was performed to obtain further biological insight into the EU and EC-associated genes. Weighted gene co-expression network analysis (WGCNA) was performed to find clinically significant modules of highly-correlated genes. The hub genes that belong to both the weighted gene co-expression network and protein-protein interaction (PPI) network were identified using a Venn diagram. Results: We obtained EC and EU-associated microarray datasets GSE7305 and GSE120103. Genes in the EC were mainly enriched in the immune response and immune cell trafficking, and genes in the EU were mainly enriched in stress response and steroid hormone biosynthesis. PPI networks and weighted gene co-expression networks were constructed. An EC-associated blue module and an EU-associated magenta module were identified, and their function annotations revealed that hormone receptor signaling or inflammatory microenvironments may promote EU passing through the oviducts and migrating to the ovarian surfaces, and adhesion and immune correlated genes may induce the successful ectopic implantation of the endometrium (EC). Twelve hub genes in the EC and sixteen hub genes in the EU were recognized and further validated in independent datasets. Conclusion: Our study identified, for the first time, the hub genes and enrichment pathways in the EC and EU using WGCNA, which may provide a comprehensive understanding of the pathogenesis of endometriosis and have important clinical implications for the treatment and diagnosis of endometriosis.

The Landscape of Genetic Alterations Stratified Prognosis in Oriental Pancreatic Cancer Patients

Background

Pancreatic cancer is a life-threatening malignant disease with significant diversity among geographic regions and races leading to distinct carcinogenesis and prognosis. Previous studies mainly focused on Western patients, while the genomic landscape of Oriental patients, especially Chinese, remained less investigated.

Methods

A total of 408 pancreatic cancer patients were enrolled. A panel containing 436 cancer-related genes was used to detect genetic alterations in tumor samples.

Results

We profiled the genomic alteration landscape of pancreatic duct adenocarcinoma (PDAC), intraductal papillary mucinous neoplasm (IPMN), periampullary carcinoma (PVC), and solid-pseudopapillary tumor (SPT). Comparison with a public database revealed specific gene mutations in Oriental PDAC patients including higher mutation rates of DNA damage repair-related genes. Analysis of mutational signatures showed potential heterogenous carcinogenic factors caused by diabetes mellitus. KRAS mutation, especially KRAS G12D mutation, was associated with poor survival, while patients not harboring the 17 significant copy number variations (CNVs) had a better prognosis. We further identified multiple correlations between clinicopathologic variables and genetic mutations, as well as CNVs. Finally, by network-based stratification, three classes of PDAC patients were robustly clustered. Among these, class 1 (characterized by the Fanconi anemia pathway) achieved the best outcome, while class 2 (involved in the platinum drug resistance pathway) suffered from the worst prognosis.

Conclusions

In this study, we reported for the first time the genetic alteration landscape of Oriental PDAC patients identifying many Oriental-specific alterations. The relationship between genetic alterations and clinicopathological factors as well as prognosis demonstrated important genomic impact on tumor biology. This study will help to optimize clinical treatment of Oriental PDAC patients and improve their survival.

Network analysis with biological data of cancer patients: A scoping review

BACKGROUND & OBJECTIVE

Network Analysis (NA) is a mathematical method that allows exploring relations between units and representing them as a graph. Although NA was initially related to social sciences, the past two decades was introduced in Bioinformatics. The recent growth of the networks' use in biological data analysis reveals the need to further investigate this area. In this work, we attempt to identify the use of NA with biological data, and specifically: (a) what types of data are used and whether they are integrated or not, (b) what is the purpose of this analysis, predictive or descriptive, and (c) the outcome of such analyses, specifically in cancer diseases.

METHODS & MATERIALS

The literature review was conducted on two databases, PubMed & IEEE, and was restricted to journal articles of the last decade (January 2010 - December 2019). At a first level, all articles were screened by title and abstract, and at a second level the screening was conducted by reading the full text article, following the predefined inclusion & exclusion criteria leading to 131 articles of interest. A table was created with the information of interest and was used for the classification of the articles. The articles were initially classified to analysis studies and studies that propose a new algorithm or methodology. Each one of these categories was further screened by the following clustering criteria: (a) data used, (b) study purpose, (c) study outcome. Specifically for the studies proposing a new algorithm, the novelty presented in each one was detected.

RESULTS

& Conclusions: In the past five years researchers are focusing on creating new algorithms and methodologies to enhance this field. The articles' classification revealed that only 25% of the analyses are integrating multi-omics data, although 50% of the new algorithms developed follow this integrative direction. Moreover, only 20% of the analyses and 10% of the newly developed methodologies have a predictive purpose. Regarding the result of the works reviewed, 75% of the studies focus on identifying, prognostic or not, gene signatures. Concluding, this review revealed the need for deploying predictive and multi-omics integrative algorithms and methodologies that can be used to enhance cancer diagnosis, prognosis and treatment.

Peripheral nerve injury and sensitization underlie pain associated with oral cancer perineural invasion

Cancer invading into nerves, termed perineural invasion (PNI), is associated with pain. Here, we show that oral cancer patients with PNI report greater spontaneous pain and mechanical allodynia compared with patients without PNI, suggesting that unique mechanisms drive PNI-induced pain. We studied the impact of PNI on peripheral nerve physiology and anatomy using a murine sciatic nerve PNI model. Mice with PNI exhibited spontaneous nociception and mechanical allodynia. Perineural invasion induced afterdischarge in A high-threshold mechanoreceptors (HTMRs), mechanical sensitization (ie, decreased mechanical thresholds) in both A and C HTMRs, and mechanical desensitization in low-threshold mechanoreceptors. Perineural invasion resulted in nerve damage, including axon loss, myelin damage, and axon degeneration. Electrophysiological evidence of nerve injury included decreased conduction velocity, and increased percentage of both mechanically insensitive and electrically unexcitable neurons. We conclude that PNI-induced pain is driven by nerve injury and peripheral sensitization in HTMRs.

CARM1 Inhibition Enables Immunotherapy of Resistant Tumors by Dual Action on Tumor Cells and T Cells

A number of cancer drugs activate innate immune pathways in tumor cells but unfortunately also compromise antitumor immune function. We discovered that inhibition of CARM1, an epigenetic enzyme and cotranscriptional activator, elicited beneficial antitumor activity in both cytotoxic T cells and tumor cells. In T cells, Carm1 inactivation substantially enhanced their antitumor function and preserved memory-like populations required for sustained antitumor immunity. In tumor cells, Carm1 inactivation induced a potent type 1 interferon response that sensitized resistant tumors to cytotoxic T cells. Substantially increased numbers of dendritic cells, CD8 T cells, and natural killer cells were present in Carm1-deficient tumors, and infiltrating CD8 T cells expressed low levels of exhaustion markers. Targeting of CARM1 with a small molecule elicited potent antitumor immunity and sensitized resistant tumors to checkpoint blockade. Targeting of this cotranscriptional regulator thus offers an opportunity to enhance immune function while simultaneously sensitizing resistant tumor cells to immune attack. SIGNIFICANCE: Resistance to cancer immunotherapy remains a major challenge. Targeting of CARM1 enables immunotherapy of resistant tumors by enhancing T-cell functionality and preserving memory-like T-cell populations within tumors. CARM1 inhibition also sensitizes resistant tumor cells to immune attack by inducing a tumor cell-intrinsic type 1 interferon response.This article is highlighted in the In This Issue feature, p. 1861.

PSMC2, ORC5 and KRTDAP are specific biomarkers for HPV-negative head and neck squamous cell carcinoma

The prognosis of patients with human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC) is poorer than those with HPV-positive HNSCC. The present study aimed to identify novel and specific biomarkers of HPV-negative HNSCC using bioinformatics analysis and associated experiments. The gene expression profiles of HPV-negative HNSCC tissues and corresponding clinical data were downloaded from The Cancer Genome Atlas database and used in a weighted gene co-expression network analysis. Genes in clinically significant co-expression modules were used to construct a protein-protein interaction (PPI) network. The genes demonstrating a high degree score in the PPI network and a high correlation with tumor grade were considered hub genes. The diagnostic value of the hub genes associated with HPV-negative and HPV-positive HNSCC was analyzed using differential expression gene (DEG) analysis, immunohistochemical (IHC) staining and a receiver operating characteristic (ROC) curve analysis. Seven genes [Serrate RNA effector molecule (SRRT), checkpoint kinase 2 (CHEK2), small nuclear ribonucleoprotein polypeptide E (SNRPE), proteasome 26S subunit ATPase 2 (PSMC2), origin recognition complex subunit 5 (ORC5), S100 calcium binding protein A7 and keratinocyte differentiation associated protein (KRTDAP)] were demonstrated to be hub genes in clinically significant co-expression modules. DEG, IHC and ROC curve analyses revealed that SRRT, CHEK2 and SNRPE were significantly upregulated in HPV-negative and HPV-positive HNSCC tissues compared with in adjacent tissues, and these genes demonstrated a high diagnostic value for distinguishing HNSCC tissues. However, PSMC2, ORC5 and KRTDAP were the only differentially expressed genes identified in HPV-negative HNSCC tissues, and these genes demonstrated a high diagnostic value for HPV-negative HNSCC. PSMC2, ORC5 and KRTDAP may therefore serve as novel and specific biomarkers for HPV-negative HNSCC, potentially improving the diagnosis and treatment of patients with HPV-negative HNSCC.

Integrating Clinical and Genomic Analyses of Hippocampal-Prefrontal Circuit Disorder in Depression

Major depressive disorder (MDD) is a prevalent, devastating and recurrent mental disease. Hippocampus (HIP)-prefrontal cortex (PFC) neural circuit abnormalities have been confirmed to exist in MDD; however, the gene-related molecular features of this circuit in the context of depression remain unclear. To clarify this issue, we performed gene set enrichment analysis (GSEA) to comprehensively analyze the genetic characteristics of the two brain regions and used weighted gene correlation network analysis (WGCNA) to determine the main depression-related gene modules in the HIP-PFC network. To clarify the regional differences and consistency for MDD, we also compared the expression patterns and molecular functions of the key modules from the two brain regions. The results showed that candidate modules related to clinical MDD of HIP and PFC, which contained with 363 genes and 225 genes, respectively. Ninety-five differentially expressed genes (DEGs) were identified in the HIP candidate module, and 51 DEGs were identified in the PFC candidate module, with only 11 overlapping DEGs in these two regional modules. Combined with the enrichment results, although there is heterogeneity in the molecular functions in the HIP-PFC network of depression, the regulation of the MAPK cascade, Ras protein signal transduction and Ephrin signaling were significantly enriched in both brain regions, indicating that these biological pathways play important roles in MDD pathogenesis. Additionally, the high coefficient protein–protein interaction (PPI) network was constructed via STRING, and the top-10 coefficient genes were identified as hub genes via the cytoHubba algorithm. In summary, the present study reveals the gene expression characteristics of MDD and identifies common and unique molecular features and patterns in the HIP-PFC network. Our results may provide novel clues from the gene function perspective to explain the pathogenic mechanism of depression and to aid drug development. Further research is needed to confirm these findings and to investigate the genetic regulation mechanisms of different neural networks in depression.

Identification of Infertility-Associated Topologically Important Genes Using Weighted Co-expression Network Analysis

Endometriosis has been associated with a high risk of infertility. However, the underlying molecular mechanism of infertility in endometriosis remains poorly understood. In our study, we aimed to discover topologically important genes related to infertility in endometriosis, based on the structure network mining. We used microarray data from the Gene Expression Omnibus (GEO) database to construct a weighted gene co-expression network for fertile and infertile women with endometriosis and to identify gene modules highly correlated with clinical features of infertility in endometriosis. Additionally, the protein–protein interaction network analysis was used to identify the potential 20 hub messenger RNAs (mRNAs) while the network topological analysis was used to identify nine candidate long non-coding RNAs (lncRNAs). Functional annotations of clinically significant modules and lncRNAs revealed that hub genes might be involved in infertility in endometriosis by regulating G protein-coupled receptor signaling (GPCR) activity. Gene Set Enrichment Analysis showed that the phospholipase C-activating GPCR signaling pathway is correlated with infertility in patients with endometriosis. Taken together, our analysis has identified 29 hub genes which might lead to infertility in endometriosis through the regulation of the GPCR network.

Perineural invasion in head and neck squamous cell carcinoma: background, mechanisms, and prognostic implications

PURPOSE OF REVIEW

Perineural invasion (PNI) is a pathological feature frequently observed in head and neck squamous cell carcinoma (HNSCC). The difficulties of pathological standardization and the lack of a simple validated experimental model to study PNI render its analysis complex. Here, we aim to summarize the recent advances in the understanding of the biology of PNI in HNSCC and their potential clinical implications.

RECENT FINDINGS

PNI is a multistep process leading to a dialogue between cancer cells and nerve fibers in HNSCC. Recent studies have identified some of the active molecular mechanisms involved in PNI. Comprehensive studies addressing the transcriptional regulation of PNI bring interesting perspectives for a standardized molecular diagnosis of PNI and a better assessment of its contribution to the aggressiveness of HNSCC.

SUMMARY

Perineural invasion is a complex process that reflects specific tumor biology. In addition to unveiling new fundamental concepts about the tumor microenvironment, research on PNI promises to identify new biomarkers, enabling progress in therapeutic development against HNSCC.

A HPV16-related prognostic indicator for head and neck squamous cell carcinoma

Background

The human papillomavirus (HPV) is emerging as an important risk factor in head and neck squamous cell carcinoma (HNSCC) patients. This has been observed particularly in the case of HPV16. The HPV16+ HNSCC subtype has distinct pathological, clinical, molecular, and prognostic characteristics. This study aimed to identify potential microRNAs (miRNAs) and their roles in HPV16+ HNSCC progression.

Method

miRNA, mRNA and the clinical data of 519 HNSCC and 44 HNSCC-negative samples were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed miRNAs (DEMs) in HPV16-related HNSCC tissues with prognostic value were selected. DEM levels were assessed based on clinicopathological parameters and overall survival (OS). Target genes were also predicted and functional analysis based on Gene Set Enrichment Analysis (GSEA) were then performed.

Results

In HPV16+ HNSCC tissues, miR-99a-3p and miR-4746-5p were significantly upregulated. In contrast, miR-411-5p was shown to be downregulated. miR-99a-3p^highmiR-411-5p^lowmiR-4746-5p^high expression could estimate improved OS and low frequent perineural invasion (PNI). Predicted target genes were enriched in cell growth, neuroepithelial cell differentiation, MAPK and FoxO signaling pathways. Epithelial mesenchymal transition (EMT) gene set and invasion related genes were downregulated in miR-99a-3p^highmiR-411-5p^lowmiR-4746-5p^high HNSCC patients.

Conclusion

miR-99a-3p, miR-411-5p and miR-4746-5p might participate in HPV16+ HNSCC progression through EMT related pathways and affect prognosis.

Comprehensive in silico analysis for identification of novel candidate target genes, including DHX36, OPA1, and SENP2, located on chromosome 3q in head and neck cancers

BACKGROUND

Major milestones of head and neck carcinogenesis have been associated with various genetic abnormalities; however, a clear picture of the molecular networks deregulated during the carcinogenesis of head and neck squamous cell carcinoma (HNSC) has not yet completely revealed.

METHODS

In this study, we used in silico tools and online data sets to evaluate the underlying reasons for the expressional changes of genes residing within the chromosome 3q and to help understanding their contributions to HNSC carcinogenesis.

RESULTS

We found that 13 of 20 most upregulated genes in HNSC are localized to 3q. Further analysis revealed a gene signature consisting of DHX36, OPA1, and SENP2, which showed significant correlation in HNSC samples and potentially be deregulated through similar mechanisms including DNA amplification, transcriptional, and posttranscriptional regulation.

CONCLUSIONS

Considering our findings, we suggest DHX36, OPA1, and SENP2 genes as overexpressed in HNSC tumors and that might be concurrently involved in HNSC carcinogenesis, tumor progression, and induction of angiogenic pathways.

Neural invasion: a scenic trail for the nervous tumor and hidden therapeutic opportunity

Neural invasion (NI) is the invasion of cancer cells into nerves, influencing the pathological characteristics of malignant tumors. NI promotes metastasis and is associated with reduced survival of affected patients. Although known for decades, its prognostic and therapeutic implications have not been not much appreciated due to the scattered information available on its clinical complications. The use of multiple nomenclatures to describe NI also generated confusions among researchers to understand this pathological process. Here, we discuss the multiple classifications of NI and review its clinical complications. Recent findings of the regulatory roles of nerves on tumor growth have fuelled research in this field, and there has been several attempts to molecularly define the NI interface and the cancer cells involved. Therefore, in this review, we discuss the large datasets available to characterize the cancer cells in NI and also discuss the roles of Schwann cells and macrophages participating in NI.

Integrated bioinformatics analysis to decipher molecular mechanism of compound Kushen injection for esophageal cancer by combining WGCNA with network pharmacology

Compound Kushen injection (CKI), a medicine in widespread clinical use in China, has proven therapeutic effects on cancer. However, few molecular mechanism analyses have been carried out. To address this problem, bioinformatics approaches combining weighted gene co-expression network analysis with network pharmacology methods were undertaken to elucidate the underlying molecular mechanisms of CKI in the treatment of esophageal cancer (ESCA). First, the key gene modules related to the clinical traits of ESCA were analysed by WCGNA. Based on the results, the hub genes related to CKI treatment for ESCA were explored through network pharmacology. Molecular docking simulation was performed to recognize the binding activity of hub genes with CKI compounds. The results showed that the potential hub targets, including EGFR, ErbB2, CCND1 and IGF1R, are therapeutic targets of CKI for the treatment of ESCA. Moreover, these targets were significantly enriched in many pathways related to cancer and signalling pathways, such as the PI3K-Akt signalling pathway and ErbB signalling pathway. In conclusion, this research partially highlighted the molecular mechanism of CKI in the treatment of ESCA, offering great potential in the identification of the effective compounds in CKI and biomarkers for ESCA treatment.

MMP2 as an independent prognostic stratifier in oral cavity cancers

Background

Around 25% of oral cavity squamous cell carcinoma (OCSCC) are not controlled by the standard of care, but there is currently no validated biomarker to identify those patients. Our objective was to determine a robust biomarker for severe OCSCC, using a biology-driven strategy.

Patients and methods

Tumor and juxtatumor secretome were analyzed in a prospective discovery cohort of 37 OCSCC treated by primary surgery. Independent biomarker validation was performed by RTqPCR in a retrospective cohort of 145 patients with similar clinical features. An 18-gene signature (18 G) predictive of the response to PD-1 blockade was evaluated in the same cohort.

Results

Among 29 deregulated molecules identified in a secretome analysis, including chemokines, cytokines, growth factors, and molecules related to tumor growth and tissue remodeling, only soluble MMP2 was a prognostic biomarker. In our validation cohort, high levels of MMP2 and CD276, and low levels of CXCL10 and STAT1 mRNA were associated with poor prognosis in univariate analysis (Kaplan-Meier). MMP2 (p = .001) and extra-nodal extension (ENE) (p = .006) were independent biomarkers of disease-specific survival (DSS) in multivariate analysis and defined prognostic groups with 5-year DSS ranging from 36% (MMP2highENE+) to 88% (MMP2lowENE-). The expression of 18 G was similar in the different prognostic groups, suggesting comparable responsiveness to anti-PD-1.

Conclusion

High levels of MMP2 were an independent and validated prognostic biomarker, surpassing other molecules of a large panel of the tumor and immune-related processes, which may be used to select poor prognosis patients for intensified neoadjuvant or adjuvant regimens.

Weighted gene co-expression network analysis identifies CCNA2 as a treatment target of prostate cancer through inhibiting cell cycle

Prostate cancer is a malignant tumor disease that seriously harms the lives of middle-aged and elderly men. Weighted gene co-expression analysis can be used to construct gene co-expression networks to explore gene sets and genes that are significantly correlated with clinical features. In this study, the transcriptome data of prostate cancer on TCGA was analyzed by weighted gene co-expression network, and the gene with a significant correlation with disease Gleason stage and tumor T stage was identified: CCNA2. CCNA2 was significantly associated with biochemical recurrence, disease-free survival and overall survival rate of prostate cancer. The ability of cancer cell proliferation, invasion and metastasis was decreased after down-regulated expression of CCNA2 in prostate cancer cell lines. Flow cytometry revealed that tumor cells were arrested in the S phase after down-regulated the expression of CCNA2. Taken together, we used WGCNA and obtain a gene CCNA2 which is significantly associated with the prognosis of prostate cancer, which may be an indicator of the prognosis of prostate cancer and a new therapeutic target.

SEMA3C Promotes Cervical Cancer Growth and Is Associated With Poor Prognosis

Introduction: Aberrant activation of Semaphorin3C(SEMA3C) is widespread in human cancers. We aimed to analyze SEMA3C expression in cervical cancer and investigate the role of SEMA3C in cervical cancer and its underlying mechanism, which is important for exploring new therapeutic targets and prognostic factors.

Materials and Methods: The expression of SEMA3C was examined in paraffin-embedded cervical cancer specimens. In vivo and in vitro assays were performed to validate the effect of SEMA3C on cervical cancer cell proliferation and p-ERK pathway activation. Gene Set Enrichment Analysis (GSEA) was performed using The Cancer Genome Atlas (TCGA) data set.

Results: SEMA3C expression was associated with poor survival in both the TCGA cohort and our cohort. Silencing of SEMA3C suppressed cervical cancer cell proliferation, colony formation ability, and the activation of the p-ERK signaling pathway in vitro. SEMA3C depletion inhibited tumor growth in vitro. GSEA also showed that the epithelial mesenchymal transition (EMT), TGFβ signaling pathway, angiogenesis, and extracellular matrix (ECM) receptor interactions are associated with a high SEMA3C expression phenotype.

Conclusion: SEMA3C is correlated with poor prognosis of cervical cancer patients and promotes tumor growth via the activation of the p-ERK pathway.