Tumor-Infiltrating CD4+ Central Memory T Cells Correlated with Favorable Prognosis in Oral Squamous Cell Carcinoma

M6A-mediated molecular patterns and tumor microenvironment infiltration characterization in nasopharyngeal carcinoma

N6-methyladenosine (m6A) is the most predominant RNA epigenetic regulation in eukaryotic cells. Numerous evidence revealed that m6A modification exerts a crucial role in the regulation of tumor microenvironment (TME) cell infiltration in several tumors. Nevertheless, the potential role and mechanism of m6A modification in nasopharyngeal carcinoma (NPC) remains unknown. mRNA expression data and clinical information from GSE102349, and GSE53819 datasets obtained from Gene Expression Omnibus (GEO) was used for differential gene expression and subsequent analysis. Consensus clustering was used to identify m6A-related molecular patterns of 88 NPC samples based on prognostic m6A regulators using Univariate Cox analysis. The TME cell-infiltrating characteristics of each m6A-related subclass were explored using single-sample gene set enrichment (ssGSEA) algorithm and CIBERSORT algotithm. DEGs between two m6A-related subclasses were screened using edgeR package. The prognostic signature and predicated nomogram were constructed based on the m6A-related DEGs. The cell infiltration and expression of prognostic signature in NPC was determined using immunohistochemistry (IHC) analysis. Chi-square test was used to analysis the significance of difference of the categorical variables. And survival analysis was performed using Kaplan-Meier plots and log-rank tests. The NPC samples were divided into two m6A-related subclasses. The TME cell-infiltrating characteristics analyses indicated that cluster 1 is characterized by immune-related and metabolism pathways activation, better response to anit-PD1 and anti-CTLA4 treatment and chemotherapy. And cluster 2 is characterized by stromal activation, low expression of HLA family and immune checkpoints, and a worse response to anti-PD1 and anti-CTLA4 treatment and chemotherapy. Furthermore, we identified 1558 DEGs between two m6A-related subclasses and constructed prognostic signatures to predicate the progression-free survival (PFS) for NPC patients. Compared to non-tumor samples, REEP2, TMSB15A, DSEL, and ID4 were upregulated in NPC samples. High expression of REEP2 and TMSB15A showed poor survival in NPC patients. The interaction between REEP2, TMSB15A, DSEL, ID4, and m6A regulators was detected. Our finding indicated that m6A modification plays an important role in the regulation of TME heterogeneity and complexity.

Identification of a novel apoptosis-related genes signature to improve gastric cancer prognosis prediction

Dysregulation of apoptosis occurs in different types of malignant tumors and is likely to influence the tumor evolution, as well as clinical prognosis. However, the limited number of studies investigating the predictive power of apoptosis-related genes (ARGs) in gastric cancer indicates a gap in the current research. 174 ARGs who differentially expressed were screened using public databases, including the Gene Expression Omnibus and the Molecular Signatures Database. Univariate and LASSO regression analyses were rigorous approaches to recognize the 12 optimal genes (CTHRC1, PDGFRL, VCAN, GJA1, LOX, UPP1, ANGPT2, CRIM1, HIF1A, APOD, RNase1, and ID1) that make up the prognostic risk model. Molecular mutations, related signaling pathways, and immune system characteristics in different subgroups defined by the risk model were analyzed using different R packages. Moreover, based on the database of Genomics of Drug Sensitivity in Cancer, chemotherapy sensitivity was predicted among the risk subgroups. As a result, there were differences in mutation profiles, signaling pathways, and infiltrated immune cells between patients in various risk groups. Moreover, the low-risk group displayed greater sensitivity to chemotherapy than the high-risk group. Risk model provided a better prognostic value than the T, N, and M stages, according to the receiver operating characteristic curve. Finally, in a nomogram, the risk model and clinical factors were combined to visualize the survival rates of patients with GC. In response to the differential expression of apoptosis-related genes, a novel model for predicting the prognosis of GC patients was developed. This model may be highly valuable for guiding doctors to deliver treatment plans tailored to the need of patients with GC.

FAM122A functions as a tumor suppressor in oral squamous cell carcinoma

Family with sequence similarity 122a (FAM122A) , identified as an endogenous inhibitor of protein phosphatase 2A (PP2A) previously, is involved in multiple important physiological processes, and essential for the growth of acute myeloid leukemia and hepatocellular carcinoma cells. However, the function of FAM122A in oral squamous cell carcinoma (OSCC) is undetermined. In this study, by analyzing TCGA and GEO databases, we found that the expression of FAM122A was significantly down-regulated in head and neck squamous cell carcinoma and OSCC patients, meanwhile this low expression was tightly associated with the poor prognosis and advanced clinical stage during OSCC development. The similar low expression pattern of FAM122A could also been seen in OSCC cell lines compared with normal human oral keratinocytes. Further, we demonstrated that FAM122A knockdown significantly promoted the growth, clonogenic potential as well as migration capabilities of OSCC cells, while these alterations could be rescued by the re-expression of FAM122A. Over-expression of FAM122A suppressed OSCC cell proliferation and migration. FAM122A also inhibited the epithelial-mesenchymal transition (EMT) in OSCC cells by the up-regulation of epithelial marker E-cadherin and down-regulation of mesenchymal markers Fibronectin and Vimentin, which is presumably mediated by transforming growth factor β receptor 3 (TGFBR3), a novel tumor suppressor. In addition, FAM122A could induce T cell infiltration in OSCC, indicating that FAM122A might influence the immune cell activity of tumor environment and further interfere the tumor development. Collectively, our results suggest that FAM122A functions as a tumor suppressor in OSCC and possibly acts as a predictive biomarker for the diagnosis and/or treatment of OSCC.

Comprehensive molecular characterization of long-term glioblastoma survivors

Fewer than 5 % glioblastoma (GBM) patients survive over five years and are termed long-term survivors (LTS), yet their molecular background is unclear. The present cohort included 72 isocitrate dehydrogenase (IDH)-wildtype GBM patients, consisting of 35 LTS and 37 short-term survivors (STS), and we employed whole exome sequencing, RNA-seq and DNA methylation array to delineate this largest LTS cohort to date. Although LTS and STS demonstrated analogous clinical characters and classical GBM biomarkers, CASC5 (P = 0.002) and SPEN (P = 0.013) mutations were enriched in LTS, whereas gene-to-gene fusions were concentrated in STS (P = 0.007). Importantly, LTS exhibited higher tumor mutation burden (P < 0.001) and copy number (CN) increase (P = 0.013), but lower mutant-allele tumor heterogeneity score (P < 0.001) and CN decrease (P = 0.026). Additionally, LTS demonstrated hypermethylated genome (P < 0.001) relative to STS. Differentially expressed and methylated genes both enriched in olfactory transduction. Further, analysis of the tumor microenvironment revealed higher infiltration of M1 macrophages (P = 0.043), B cells (P = 0.016), class-switched memory B cells (P = 0.002), central memory CD4+ T cells (P = 0.031) and CD4+ Th1 cells (P = 0.005) in LTS. We also separately analyzed a subset of patients who were methylation class-defined GBM, contributing 70.8 % of the entire cohort, and obtained similar results relative to prior analyses. Finally, we demonstrated that LTS and STS could be distinguished using a subset of molecular features. Taken together, the present study delineated unique molecular attributes of LTS GBM.

Identification of ferroptosis related genes and pathways in prostate cancer cells under erastin exposure

BACKGROUND

Few studies are focusing on the mechanism of erastin acts on prostate cancer (PCa) cells, and essential ferroptosis-related genes (FRGs) that can be PCa therapeutic targets are rarely known.

METHODS

In this study, in vitro assays were performed and RNA-sequencing was used to measure the expression of differentially expressed genes (DEGs) in erastin-induced PCa cells. A series of bioinformatic analyses were applied to analyze the pathways and DEGs.

RESULTS

Erastin inhibited the expression of SLC7A11 and cell survivability in LNCaP and PC3 cells. After treatment with erastin, the concentrations of malondialdehyde (MDA) and Fe2+ significantly increased, whereas the glutathione (GSH) and the oxidized glutathione (GSSG) significantly decreased in both cells. A total of 295 overlapping DEGs were identified under erastin exposure and significantly enriched in several pathways, including DNA replication and cell cycle. The percentage of LNCaP and PC3 cells in G1 phase was markedly increased in response to erastin treatment. For four hub FRGs, TMEFF2 was higher in PCa tissue and the expression levels of NRXN3, CLU, and UNC5B were lower in PCa tissue. The expression levels of SLC7A11 and cell survivability were inhibited after the knockdown of TMEFF2 in androgen-dependent cell lines (LNCaP and VCaP) but not in androgen-independent cell lines (PC3 and C4-2). The concentration of Fe2+ only significantly increased in TMEFF2 downregulated LNCaP and VCaP cells.

CONCLUSION

TMEFF2 might be likely to develop into a potential ferroptosis target in PCa and this study extends our understanding of the molecular mechanism involved in erastin-affected PCa cells.

Overexpression of CD73 is associated with recurrence and poor prognosis of gingivobuccal oral cancer as revealed by transcriptome and deep immune profiling of paired tumor and margin tissues

BACKGROUND

For various cancers, differences in response to treatment and subsequent survival period have been reported to be associated with variation in immune contextures.

AIM

We sought to identify whether such association exists in respect of gingivobuccal oral cancer.

MATERIALS AND METHODS

We performed deep immune profiling of tumor and margin tissues collected from 46 treatment naïve, Human Papillomavirus (HPV) negative, patients. Each patient was followed for 24 months and prognosis (recurrence/death) noted. Key findings were validated by comparing with TCGA-HNSC cohort data.

RESULTS

About 28% of patients showed poor post-treatment prognosis. These patients exhibited a high probability of recurrence even within 1 year and death within 2 years. There was restricted immune cell infiltration in tumor, but not in margin, among these patients. Reduced expression of eight immune-related genes (IRGs) (NT5E, THRA, RBP1, TLR4, ITGA6, BMPR1B, ITGAV, SSTR1) in tumor strongly predicted better quality of prognosis, both in our patient cohort and in TCGA-HNSC cohort. Tumors of patients with better prognosis were associated with (a) lower CD73+ cells with concomitant lower expression level of NT5E/CD73, (b) higher proportions of CD4+ and CD8+ T cells, B cells, NK cells, M1 macrophages, (c) higher %Granzyme+ cells, (d) higher TCR and BCR repertoire diversities. CD73 expression in tumor was associated with low CD8+ and CD4+ T cells, low immune repertoire diversity, and advanced cancer stage.

DISCUSSION AND CONCLUSION

High infiltration of anti-tumor immune cells in both tumors and margins results in good prognosis, while in patients with minimal infiltration in tumors in spite of high infiltration in margins results in poor prognosis. Targeted CD73 immune-checkpoint inhibition may improve clinical outcome.

Targeting oral tumor microenvironment for effective therapy

Oral cancers are among the common head and neck malignancies. Different anticancer therapy modalities such as chemotherapy, immunotherapy, radiation therapy, and also targeted molecular therapy may be prescribed for targeting oral malignancies. Traditionally, it has been assumed that targeting malignant cells alone by anticancer modalities such as chemotherapy and radiotherapy suppresses tumor growth. In the last decade, a large number of experiments have confirmed the pivotal role of other cells and secreted molecules in the tumor microenvironment (TME) on tumor progression. Extracellular matrix and immunosuppressive cells such as tumor-associated macrophages, myeloid-derived suppressor cells (MDSCs), cancer-associated fibroblasts (CAFs), and regulatory T cells (Tregs) play key roles in the progression of tumors like oral cancers and resistance to therapy. On the other hand, infiltrated CD4 + and CD8 + T lymphocytes, and natural killer (NK) cells are key anti-tumor cells that suppress the proliferation of malignant cells. Modulation of extracellular matrix and immunosuppressive cells, and also stimulation of anticancer immunity have been suggested to treat oral malignancies more effectively. Furthermore, the administration of some adjuvants or combination therapy modalities may suppress oral malignancies more effectively. In this review, we discuss various interactions between oral cancer cells and TME. Furthermore, we also review the basic mechanisms within oral TME that may cause resistance to therapy. Potential targets and approaches for overcoming the resistance of oral cancers to various anticancer modalities will also be reviewed. The findings for targeting cells and potential therapeutic targets in clinical studies will also be reviewed.

Construction of a novel mRNAsi-related risk model for predicting prognosis and immunotherapy response in osteosarcoma

Background

Targeting cancer stem cells (CSC) may represent a future therapeutic direction for osteosarcoma (OS), which mainly relies on the identification of CSC markers. This study aimed to classify OS based on messenger ribonucleic acid (mRNA) stemness indices (mRNAsi) and construct a mRNAsi-related risk model to predict the prognosis of OS.

Methods

The one-class logistic regression (OCLR) algorithm was applied to the RNA- sequencing (seq) data of human embryonic stem cells (hESC) and induced pluripotent stem cell (iPSC) lines to calculate mRNAsi. Weighted gene co-expression network analysis (WGCNA) was performed on data obtained from the TARGET database to screen the mRNAsi-related genes. Univariate Cox regression analysis was implemented to screen mRNAsi-related genes with prognostic significance for consensus clustering of OS. The least absolute shrinkage and selection operator (LASSO) and COX regression analysis were conducted to construct a risk model based on mRNAsi-related genes.

Results

Six gene modules were identified in the TARGET database. The yellow module showed the strongest negative correlation with mRNAsi and the strongest significant positive correlation with the immune score and stromal score. OS was divided into three molecular subtypes with significant survival differences based on 73 mRNAsi-related genes with prognostic value for OS. The survival rate was ranked as C3 < C1 < C2 from low to high. The levels of immune components in C2 was significantly higher than those in C1 and C3. HSD11B2, GBP1, RNF130, APBB1IP, and NPC2 in the yellow module were used as variables for building the mRNAsi-related risk model. The survival rate of the high-risk group (as defined by this model) was significantly higher than that of the low-risk group, and it had significant survival prediction ability in 28 types of cancer. In addition, the mRNAsi-related risk model was superior to the Tumor Immune Dysfunction and Exclusion (TIDE) model in predicting the prognosis and immunotherapy response in all three immunotherapy cohorts.

Conclusions

This study classified OS and constructed a mRNAsi-related risk model based on mRNAsi-related genes, which provides a potential tool for more accurate risk stratification of OS and prediction of immunotherapy response.

Low-density lipoprotein balances T cell metabolism and enhances response to anti-PD-1 blockade in a HCT116 spheroid model

Introduction

The discovery of immune checkpoints and the development of their specific inhibitors was acclaimed as a major breakthrough in cancer therapy. However, only a limited patient cohort shows sufficient response to therapy. Hence, there is a need for identifying new checkpoints and predictive biomarkers with the objective of overcoming immune escape and resistance to treatment. Having been associated with both, treatment response and failure, LDL seems to be a double-edged sword in anti-PD1 immunotherapy. Being embedded into complex metabolic conditions, the impact of LDL on distinct immune cells has not been sufficiently addressed. Revealing the effects of LDL on T cell performance in tumor immunity may enable individual treatment adjustments in order to enhance the response to routinely administered immunotherapies in different patient populations. The object of this work was to investigate the effect of LDL on T cell activation and tumor immunity in-vitro.

Methods

Experiments were performed with different LDL dosages (LDL^low = 50 μg/ml and LDL^high = 200 μg/ml) referring to medium control. T cell phenotype, cytokines and metabolism were analyzed. The functional relevance of our findings was studied in a HCT116 spheroid model in the context of anti-PD-1 blockade.

Results

The key points of our findings showed that LDL^high skewed the CD4⁺ T cell subset into a central memory-like phenotype, enhanced the expression of the co-stimulatory marker CD154 (CD40L) and significantly reduced secretion of IL-10. The exhaustion markers PD-1 and LAG-3 were downregulated on both T cell subsets and phenotypical changes were associated with a balanced T cell metabolism, in particular with a significant decrease of reactive oxygen species (ROS). T cell transfer into a HCT116 spheroid model resulted in a significant reduction of the spheroid viability in presence of an anti-PD-1 antibody combined with LDL^high.

Discussion

Further research needs to be conducted to fully understand the impact of LDL on T cells in tumor immunity and moreover, to also unravel LDL effects on other lymphocytes and myeloid cells for improving anti-PD-1 immunotherapy. The reason for improved response might be a resilient, less exhausted phenotype with balanced ROS levels.

Sintilimab versus docetaxel as second-line treatment in advanced or metastatic squamous non-small-cell lung cancer: an open-label, randomized controlled phase 3 trial (ORIENT-3)

BACKGROUND

Treatment options for Chinese patients with locally advanced or metastatic squamous-cell non-small-cell lung cancer (sqNSCLC) after failure of first-line chemotherapy are limited. This study (ORIENT-3) aimed to evaluate the efficacy and safety of sintilimab versus docetaxel as second-line treatment in patients with locally advanced or metastatic sqNSCLC.

METHODS

ORIENT-3 was an open-label, multicenter, randomized controlled phase 3 trial that recruited patients with stage IIIB/IIIC/IV sqNSCLC after failure with first-line platinum-based chemotherapy. Patients were randomized in a 1:1 ratio to receive either 200 mg of sintilimab or 75 mg/m² of docetaxel intravenously every 3 weeks, stratified by the Eastern Cooperative Oncology Group performance status. The primary endpoint was overall survival (OS) in the full analysis set (FAS). Secondary endpoints included progression-free survival (PFS), objective response rate (ORR), disease control rate (DCR), duration of response (DoR) and safety.

RESULTS

Between August 25, 2017, and November 7, 2018, 290 patients were randomized. For FAS, 10 patients from the docetaxel arm were excluded. The median OS was 11.79 (n = 145; 95% confidence interval [CI], 10.28-15.57) months with sintilimab versus 8.25 (n = 135; 95% CI, 6.47-9.82) months with docetaxel (hazard ratio [HR]: 0.74; 95% CI, 0.56-0.96; P = 0.025). Sintilimab treatment significantly prolonged PFS (median 4.30 vs. 2.79 months; HR: 0.52; 95% CI, 0.39-0.68; P < 0.001) and showed higher ORR (25.50% vs. 2.20%, P < 0.001) and DCR (65.50% vs. 37.80%, P < 0.001) than the docetaxel arm. The median DoR was 12.45 (95% CI, 4.86-25.33) months in the sintilimab arm and 4.14 (95% CI, 1.41-7.23) months in the docetaxel arm (P = 0.045). Treatment-related adverse events of grade ≥ 3 were reported in 26 (18.1%) patients in the sintilimab arm and 47 (36.2%) patients in the docetaxel arm. Exploratory biomarker analysis showed potential predictive values of expression levels of two transcription factors, including OVOL2 (HR: 0.35; P < 0.001) and CTCF (HR: 3.50; P < 0.001)，for sintilimab treatment.

CONCLUSIONS

Compared with docetaxel, sintilimab significantly improved the OS, PFS, and ORR of Chinese patients with previously treated locally advanced or metastatic sqNSCLC.

Prognostic Role of Combined EGFR and Tumor-Infiltrating Lymphocytes in Oral Squamous Cell Carcinoma

Background

Epidermal growth factor receptor (EGFR) is well known as a general prognostic biomarker for head and neck tumors, however the specific prognostic value of EGFR in oral squamous cell carcinoma (OSCC) is controversial. Recently, the presence of tumor-infiltrating T cells has been associated with significant survival advantages in a variety of disease sites. The present study will determine if the inclusion of T cell specific markers (CD3, CD4 and CD8) would enhance the prognostic value of EGFR in OSCCs.

Methods

Tissue microarrays containing 146 OSCC cases were analyzed for EGFR, CD3, CD4 and CD8 expression using immunohistochemical staining. EGFR and T cell expression scores were correlated with clinicopathological parameters and survival outcomes.

Results

Results showed that EGFR expression had no impact on overall survival (OS), but EGFR-positive (EGFR+) OSCC patients demonstrated significantly worse progression free survival (PFS) compared to EGFR-negative (EGFR-) patients. Patients with CD3, CD4 and CD8-positive tumors had significantly better OS compared to CD3, CD4 and CD8-negative patients respectively, but no impact on PFS. Combined EGFR+/CD3+ expression was associated with cases with no nodal involvement and significantly more favorable OS compared to EGFR+/CD3- expression. CD3 expression had no impact on OS or PFS in EGFR- patients. Combinations of EGFR/CD8 and EGFR/CD4 expression showed no significant differences in OS or PFS among the expression groups.

Conclusion

Altogether these results suggest that the expression of CD3+ tumor-infiltrating T cells can enhance the prognostic value of EGFR expression and warrants further investigation as prognostic biomarkers for OSCC.

Characteristics of Fatty Acid Metabolism in Lung Adenocarcinoma to Guide Clinical Treatment

Background

Lung adenocarcinoma (LUAD) has a very high morbidity and mortality rate, and its pathogenesis and treatment are still in the exploratory stage. Fatty acid metabolism plays a significant role in tumorigenesis, progression, and immune regulation. However, the gene expression of fatty acid metabolism in patients with LUAD and its relationship with prognosis remain unclear.

Methods

We collected 309 fatty acid metabolism-related genes, established a LUAD risk model based on The Cancer Genome Atlas (TCGA) using Least Absolute Shrinkage Selection Operator (LASSO) regression analysis, and divided LUAD patients into high-risk and low-risk groups, which were further validated using the Gene Expression Omnibus (GEO) database. The nomogram, principal component analysis (PCA), and receiver operating characteristic (ROC) curves showed that the model had the best predictive performance. The ROC curves and calibration plots confirmed that the nomogram had good predictive power. We further analyzed the differences in clinical characteristics, immune cell infiltration, immune-related functions, chemotherapy drug sensitivity, and immunotherapy efficacy between the high-risk and low-risk groups. We also analyzed the enrichment pathways and protein–protein interaction (PPI) networks of different genes in the high-risk and low-risk groups to screen for target genes and further explored the correlation between target genes and differences in survival prognosis, clinical characteristics, gene mutations, and immune cells.

Results

Risk score and staging are independent prognostic factors for patients with LUAD. The high-risk group had lower immune cell infiltration, was more sensitive to chemotherapeutic agents, and had a poorer survival prognosis. We also obtained three pivotal genes with poor survival prognosis in the high expression group, which were strongly associated with clinical symptoms and immune cells.

Conclusion

Risk score and staging are independent prognostic factors for patients with LUAD. The high-risk group had lower immune cell infiltration, was more sensitive to chemotherapeutic agents, and had a poorer survival prognosis. We also obtained three survival prognosis-associated target genes that are closely associated with clinical symptoms and immune cells and may be potential targets for immune-targeted therapy in LUAD.