RNA-Seq analysis of peripheral blood mononuclear cells reveals unique transcriptional signatures associated with radiotherapy response of nasopharyngeal carcinoma and prognosis of head and neck cancer

Synaptotagmin-1 attenuates myocardial programmed necrosis and ischemia/reperfusion injury through the mitochondrial pathway

Programmed necrosis/necroptosis greatly contributes to the pathogenesis of cardiac disorders including myocardial infarction, ischemia/reperfusion (I/R) injury and heart failure. However, the fundamental mechanism underlying myocardial necroptosis, especially the mitochondria-dependent death pathway, is poorly understood. Synaptotagmin-1 (Syt1), a Ca2+ sensor, is originally identified in nervous system and mediates synchronous neurotransmitter release. The later findings of Syt1 expressions in many non-neuronal tissues including muscles suggest that Syt1 may exert important functions beyond regulation of neurotransmitter release. Syt1 is highly expressed in cardiomyocytes and has been used as an extracellular molecular probe for SPECT imaging of cardiac cell death in acute myocardial infarction. However, whether Syt1 functions in the pathogenesis of cardiac disorders and what is the molecular etiology have not yet been clarified. We showed here that Syt1 expression was significantly down-regulated in mice I/R injured heart tissues, H2O2-challenged cardiomyocytes and hypoxia/reoxygenation (H/R)-damaged cardiomyocytes. Enforced expression of Syt1 significantly inhibited myocardial necrotic cell death and interstitial fibrosis, and improved cardiac function in mice subjected to I/R operation. In exploring the underlying mechanisms, we found that Syt1 interacted with Parkin and promoted Parkin-catalyzed CypD ubiquitination, thus inhibited mitochondrial membrane permeability transition pore (mPTP) opening and ultimately suppressed cardiomyocyte necrosis. We further found that Syt1 expression was negatively regulated by miR-193b-3p. MiR-193b-3p regulated cardiomyocyte necrosis and mPTP opening by targeting Syt1. Our present work revealed a novel regulatory model of myocardial necrosis composed of miR-193b-3p, Syt1, Parkin, and CypD, which may provide potential therapeutic targets and strategies for heart protection.

Identification of a Novel Signature Based on Ferritinophagy-Related Genes to Predict Prognosis in Lung Adenocarcinoma: Focus on AHNAK2

BACKGROUND

Lung adenocarcinoma (LUAD) accounts for over 40% of all non-small cell lung cancer (NSCLC) cases and continues to be difficult to treat despite advancements in diagnostics and therapies. Ferritinophagy, a newly recognized autophagy process linked to ferroptosis, has been associated with LUAD development. Recent studies have shown a dysregulation of genes related to ferritinophagy in LUAD, indicating its potential as a therapeutic target.

METHODS

We constructed a predictive model using seven genes associated with ferritinophagy. The model's accuracy was evaluated across three independent gene expression datasets. We analyzed the biological functions, immune environment, mutations, and drug sensitivities in groups with high and low risk. Utilizing a single-cell sequencing (scRNA-seq) dataset, we confirmed the expression of the model genes and identified a subtype of epithelial cells expressing AHNAK2. We further investigated the impact of the ferritinophagy-related gene AHNAK2 on LUAD cell proliferation, invasion, migration, and ferroptosis in vitro.

RESULTS

Our prediction model, comprising seven genes (AHNAK2, ARNTL2, CD27, LTB, SLC15A1, SLC2A1, and SYT1), has shown potential in predicting the prognosis of individuals diagnosed with LUAD. Notably, AHNAK2 impedes ferroptosis, promoting LUAD progression in vitro.

CONCLUSIONS

Our research suggests that ferritinophagy-associated genes are promising prognostic markers for LUAD and lay the groundwork for further exploration of ferritinophagy's role in LUAD. Furthermore, we present AHNAK2 as a novel regulator of ferroptosis, which requires further investigation to understand its mechanism.

Prognostic Analysis of Lactic Acid Metabolism Genes in Oral Squamous Cell Carcinoma

OBJECTIVES

Oral squamous cell carcinoma (OSCC) is the most common malignant tumour in the oral and maxillofacial region. Lactic acid accumulation in the tumour microenvironment (TME) has gained attention for its dual role as an energy source for cancer cells and an activator of signalling pathways crucial to tumour progression. This study aims to reveal the impact of lactate-related genes (LRGs) on the prognosis, TME, and immune characteristics of OSCC, with the ultimate goal of developing a novel prognostic model.

METHODS

Unsupervised clustering analysis of LRGs in OSCC patients from The Cancer Genome Atlas database was conducted to evaluate and compare TME, immune features, and clinical characteristics across various lactate subtypes. A refined prognostic model was developed through the application of Cox and Least absolute shrinkage and selection operator (LASSO) regression techniques. External validation sets were then utilised to improve model accuracy, along with a detailed correlation analysis of drug sensitivity.

RESULTS

The Cancer Genome Atlas-OSCC patients were categorised into 4 distinct lactate subtypes based on LRGs. Notably, patients in subtype 1 and subtype 2 exhibited the least and most favourable prognoses, respectively. Subtype 1 patients showed elevated expression levels of immune checkpoint genes. Further analysis identified 1086 genes with significant expression differences between cancer and noncancer tissues, as well as between subtype 1 and subtype 2 patients. Selected genes for the prognostic model included ZNF662, CGNL1, VWCE, and ZFP42. The high-risk group defined by this model had a significantly poorer prognosis (P < .0001) and functioned as an independent prognostic factor (P < .001), accurately predicting 1-, 3-, and 5-year survival rates. Additionally, individuals in the high-risk category exhibited heightened sensitivity to chemotherapy drugs such as AZ6102 and Venetoclax.

CONCLUSIONS

The predictive model based on the genes ZNF662, CGNL1, VWCE, and ZFP42 can serve as a reliable biomarker, providing accurate prognostic predictions for OSCC patients and potential opportunities for pharmaceutical interventions.

PBMC-mediated modulation of macrophage polarization in RAW264.7 cells through STAT1/STAT6 signaling cascades

Peripheral blood mononuclear cells (PBMC), sourced autologously, offer numerous advantages when procured: easier acquisition process, no in vitro amplification needed, decreased intervention and overall increased acceptability make PBMC an attractive candidate for cell therapy treatment. However, the exact mechanism by which PBMC treat diseases remains poorly understood. Immune imbalance is the pathological basis of many diseases, with macrophages playing a crucial role in this process. However, research on the role and mechanisms of PBMC in regulating macrophages remains scarce. This study employed an in vitro co-culture model of PBMC and RAW264.7 macrophages to explore the role and mechanisms of PBMC in regulating macrophages. The results showed that the co-culturing led to decreased expression of inflammatory cytokines and increased expression of anti-inflammatory cytokines in RAW264.7 or in the culture supernatant. Additionally, the pro-inflammatory, tissue matrix-degrading M1 macrophages decreased, while the anti-inflammatory, matrix-synthesizing, regenerative M2 macrophages increased in both RAW264.7 and monocytes within PBMC. Moreover, co-cultured macrophages exhibited a significantly decreased p-STAT1/STAT1 ratio, while the p-STAT6/STAT6 ratio significantly increased. This suggests that PBMC may inhibit M1 macrophage polarization by blocking STAT1 signaling cascades and may promote M2 macrophage polarization through the activation of STAT6 signaling cascades. Overall, this study sheds light on the role and mechanism of PBMC in regulating macrophages. Moreover, it was found that monocytes within co-cultured PBMC differentiated into M2 macrophages in the presence of macrophages. This finding provides experimental evidence for the use of PBMC in treating inflammatory diseases, especially macrophage-depleting inflammatory diseases such as osteoarthritis.

Predicting radiotherapy efficacy and prognosis in tongue squamous cell carcinoma through an in-depth analysis of a radiosensitivity gene signature

Background

Tongue squamous cell carcinoma (TSCC) is a prevalent tumor that affects many people worldwide. Radiotherapy is a common treatment option, but its efficacy varies greatly. This study seeks to validate the identified gene signature associated with radiosensitivity in TSCC, and its potential in predicting radiotherapy response and prognosis.

Methods

We analyzed 122 TSCC patients from TCGA database using the radiosensitivity signature and classified them into radiosensitive (RS) and radioresistant (RR) groups. Immune infiltration analysis methods were applied to investigate the immune status between different subgroups. Immunophenotype Score (IPS) and pRRophetic algorithm were employed to estimate the efficiency of treatment. A radioresistant TSCC cell line was established by gradually increasing radiation doses. Cell radiosensitivity was evaluated using the CCK-8 and colony formation assays. The expression of radiosensitivity-related genes was validated by qRT-PCR.

Results

Our study validated the predictive capacity of a previously identified "31-gene signature" in the TCGA-TSCC cohort, which effectively stratified patients into RS and RR groups. We observed that the RS group exhibited superior overall survival and progression-free survival rates relative to the RR group when treated with radiotherapy. The RS group was significantly enriched in most immune-related hallmark pathways, and may therefore benefit from immune checkpoint inhibitors. However, the RS group displayed lower sensitivity to first-line chemotherapy. A radioresistant TSCC cell line (CAL-27R) exhibited increased clonogenic potential and cell viability following irradiation, accompanied by downregulation of three radiosensitivity-related genes compared to its parental non-resistant cell (CAL-27). In addition, we constructed and validated a radiosensitivity-related prognostic index (PI) using 4 radiosensitivity-related genes associated with TSCC prognosis.

Conclusion

We assessed the ability of the radiosensitivity gene signature to predict outcomes in TSCC patients. our research provided valuable insights into the molecular pathways associated with radiosensitivity in TSCC and offered clinicians a practical tool to predict patient radiotherapy effectiveness and prognosis.

Identification of oxidative stress signatures of lung adenocarcinoma and prediction of patient prognosis or treatment response with single-cell RNA sequencing and bulk RNA sequencing data

Lung adenocarcinoma (LUAD), the most common subtype of lung cancer globally, has seen improved prognosis with advancements in diagnostic, surgical, radiotherapy, and molecular therapy techniques, while its 5-year survival rate remains low. Molecular biomarkers provide prognostic value. Oxidative stress factors, such as reactive nitrogen species and ROS, are crucial in various stages of tumor progression, influencing cell transformation, proliferation, angiogenesis, and metastasis. ROS demonstrate dual roles, affecting tumor cells, hypoxia sensitivity, and the microenvironment. Comprehensive analysis of oxidative stress in LUAD has not been conducted to date. Therefore, we systematically investigated the regulatory patterns of oxidative stress in LUAD based on oxidative stress-related genes and correlated these patterns with cellular infiltration characteristics of the tumor immune microenvironment. The model utilizes single-factor Cox analysis to screen key differential genes with prognostic value and employs least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct a prognostic-related prediction model. Ten candidate genes were selected based on this model. The risk score was constructed using the coefficients and expression levels of these ten genes. Furthermore, the impact of this risk score on overall survival (OS) was determined. Two genes with the most significant differential expression, SFTPB and S100P, were selected through qRT-PCR. Cell experiments including CCK-8, Edu, transwell assays confirmed their effects on lung cancer cells growth, consistent with the results of bioinformatics analysis. These findings suggested that this model held potential clinical value for evaluating the prognosis of lung adenocarcinoma.

Transcriptomics integrated with metabolomics reveals partial molecular mechanisms of nutritional risk and neurodevelopment in children with congenital heart disease

Purpose

To explore molecular mechanisms affecting nutritional risk and neurodevelopment in children with congenital heart disease (CHD) by combining transcriptome and metabolome analysis.

Methods

A total of 26 blood and serum samples from 3 groups of children with CHD low nutritional risk combined with normal neurodevelopment (group A), low nutritional risk combined with neurodevelopmental disorders (group B) and high nutritional risk combined with normal neurodevelopment (group C) were analyzed by transcriptome and metabolomics to search for differentially expressed genes (DEGs) and metabolites (DEMs). Functional analysis was conducted for DEGs and DEMs. Further, the joint pathway analysis and correlation analysis of DEGs and DEMs were performed.

Results

A total of 362 and 1,351 DEGs were detected in group B and C compared to A, respectively. A total of 6 and 7 DEMs were detected in group B and C compared to A in positive mode, respectively. There were 39 and 31 DEMs in group B and C compared to A in negative mode. Transcriptomic analysis indicated that neurodevelopment may be regulated by some genes such as NSUN7, SLC6A8, CXCL1 and LCN8, nutritional risk may be regulated by SLC1A3 and LCN8. Metabolome analysis and joint pathway analysis showed that tryptophan metabolism, linoleic and metabolism and glycerophospholipid metabolism may be related to neurodevelopment, and glycerophospholipid metabolism pathway may be related to nutritional risk.

Conclusion

By integrating transcriptome and metabolome analyses, this study revealed key genes and metabolites associated with nutritional risk and neurodevelopment in children with CHD, as well as significantly altered pathways. It has important clinical translational significance.

Identification of PATL1 as a prognostic and immunotherapeutic predictive factor for nasal-type natural killer/T-cell lymphoma and head and neck squamous cell carcinoma

This research examines the function of protein associated with topoisomerase II homolog 1 (PATL1) in nasal-type natural killer/T-cell lymphoma (NKTCL) and head and neck squamous cell carcinoma (HNSCC). We analyzed bulk RNA-seq data from NKTCL, nasal polyps, and normal nasal mucosa, identifying 439 differentially expressed genes. Machine learning algorithms highlighted PATL1 as a hub gene. PATL1 exhibited significant upregulation in NKTCL and HNSCC tumor samples in comparison to normal tissues, showing high diagnostic accuracy (AUC = 1.000) for NKTCL. Further analysis of local hospital data identified PATL1 as an independent prognostic risk factor for NKTCL. Data analysis of TCGA and GEO datasets revealed that high PATL1 expression correlated with poorer prognosis in HNSCC patients (p < 0.05). We also constructed a PATL1-based nomogram, which emerged as an independent prognostic predictor for HNSCC after addressing missing values. Additionally, we found a strong correlation between PATL1 and various immune cell infiltrates (e.g., activated.CD4 T cell), and a significant association with the expression of 37 immune checkpoints genes (e.g., CTLA4, PDCD1) and 20 N6-methyladenosine-related genes (e.g., ZC3H13, METTL3) (all p < 0.05). Both TCIA and TIDE algorithms suggested that PATL1 could potentially predict immunotherapy efficacy (p < 0.05). Cellular experiments demonstrated that transfection with a silencing plasmid of PATL1 significantly inhibited the malignant behaviors of SNK6 and FaDu cell lines(p < 0.05). In conclusion, our findings suggest that PATL1 may serve as a valuable prognostic and predictive biomarker in NKTCL and HNSCC, highlighting its significant role in these cancers.

Type I conventional dendritic cells and CD8+ T cells predict favorable clinical outcome of head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is one of the most common tumor entities worldwide, with human papillomavirus (HPV) infection contributing to cancer development. Conventional therapies achieve only limited efficiency, especially in recurrent or metastatic HNSCC. As the immune landscape decisively impacts the survival of patients and treatment efficacy, this study comprehensively investigated the immunological tumor microenvironment (TME) and its association with patient outcome, with special focus on several dendritic cell (DC) and T lymphocyte subpopulations. Therefore, formalin-fixed paraffin-embedded tumor samples of 56 HNSCC patients, who have undergone resection and adjuvant radiotherapy, were analyzed by multiplex immunohistochemistry focusing on the detailed phenotypic characterization and spatial distribution of DCs, CD8+ T cells, and T-helper cell subsets in different tumor compartments. Immune cell densities and proportions were correlated with clinical characteristics of the whole HNSCC cohort and different HPV- or hypoxia-associated subcohorts. Tumor stroma was highly infiltrated by plasmacytoid DCs and T lymphocytes. Among the T-helper cells and CD8+ T cells, stromal regulatory T cells and intraepithelial exhausted CD8+ T cells expressing programmed cell death protein-1 (PD-1+) and/or lymphocyte-activation gene-3 (LAG-3+) were the predominant phenotypes, indicating an immunosuppressive TME. HPV-associated tumors showed significantly higher infiltration of type I and type II conventional DCs (cDC1, cDC2) as well as several CD8+ T cell phenotypes including exhausted, activated, and proliferating T cells. On the contrary, tumors with hypoxia-associated gene signatures exhibited reduced infiltration for these immune cells. By multivariate Cox regression, immune-related prognostic factors were identified. Patient clusters defined by high infiltration of DCs and T lymphocytes combined with HPV positivity or low hypoxia showed significantly prolonged survival. Thereby, cDC1 and CD8+ T cells emerged as independent prognostic factors for local and distant recurrence. These results might contribute to the implementation of an immune cell infiltration score predicting HNSCC patients' survival and such patient stratification might improve the design of future individualized radiochemo-(immuno)therapies.

Novel prognostic biomarkers in nasopharyngeal carcinoma unveiled by mega-data bioinformatics analysis

Nasopharyngeal carcinoma (NPC) is commonly diagnosed at an advanced stage with a high incidence rate in Southeast Asia and Southeast China. However, the limited availability of NPC patient survival data in public databases has resulted in less rigorous studies examining the prediction of NPC survival through construction of Kaplan-Meier curves. These studies have primarily relied on small samples of NPC patients with progression-free survival (PFS) information or data from head and neck squamous cell carcinoma (HNSCC) studies almost without NPC patients. Thus, we coanalyzed RNA expression profiles in eleven datasets (46 normal (control) vs 160 tumor (NPC)) downloaded from the Gene Expression Omnibus (GEO) database and survival data provided by Jun Ma from Sun Yat-sen University. Then, differential analysis, gene ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and network analysis were performed using STRING database. After that, 2142 upregulated differentially expressed genes (DEGs) and 3857 downregulated DEGs were screened. Twenty-five of them were identified as hub genes, which were enriched in several pathways (cilium movement, extracellular matrix structural constituent, homologous recombination and cell cycle). Utilizing the comprehensive dataset we amassed from GEO database, we conducted a survival analysis of DEGs and subsequently constructed survival models. Seven DEGs (RASGRP2, MOCOS, TTC9, ARHGAP4, DPM3, CD37, and CD72) were identified and closely related to the survival prognosis of NPC. Finally, qRT-PCR, WB and IHC were performed to confirm the elevated expression of RASGRP2 and the decreased expression of TTC9, CD37, DPM3 and ARHGAP4, consistent with the DEG analysis. Conclusively, our findings provide insights into the novel prognostic biomarkers of NPC by mega-data bioinformatics analysis, which suggests that they may serve special targets in the treatment of NPC.

ADAM22 acts as a novel predictive biomarker for unfavorable prognosis and facilitates metastasis via PI3K/AKT signaling pathway in nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) is a type of epithelial malignancy known for its high likelihood of metastasizing to distant organs, which remains the primary obstacle in the treatment of NPC. The present study aimed to identify potential intervention target for NPC metastasis.

METHODS

The differentially expressed genes (DEGs) were firstly analyzed and intersected across various NPC related datasets in the Gene Expression Omnibus database. Subsequently, various techniques including quantitative polymerase chain reaction (qPCR), western blotting, immunohistochemistry, migration and invasion assays, in conjunction with bioinformatics and prognostic modeling, were utilized to elucidate the role of candidate genes in NPC metastasis.

RESULTS

We discerned the gene a disintegrin and metalloprotease 22 (ADAM22) as a distinct and significant factor in the progression and metastasis of NPC through five datasets. The elevated expression of ADAM22 was observed in clinical tissue and plasma samples with advanced NPC, as well as in high metastatic cells. Furthermore, we highlighted its essential role in a prognostic model that demonstrated strong prediction performance for NPC. Notably, overexpression of ADAM22 was found to enhance the aggressiveness and epithelial-mesenchymal transition (EMT) of low metastatic NPC cells, whereas the downregulation of ADAM22 resulted in suppressed effect in high metastatic cells. Delving into the mechanism, ADAM22 activated the PI3K/Akt signaling pathway through the mediation of Rac Family Small GTPase 2 (RAC2), thereby facilitating EMT and metastasis in NPC.

CONCLUSIONS

The study provided pioneering insights that ADAM22 had the potential to act as an oncogene by promoting EMT and metastasis of NPC through the RAC2-mediated PI3K/Akt signaling pathway. Thus, ADAM22 could serve as a novel prognostic indicator in NPC.

A novel CD8+ T cell-related gene signature as a prognostic biomarker in hepatocellular carcinoma

CD8+ T cells have great roles in tumor suppression and elimination of various tumors including hepatocellular carcinoma (HCC). Nonetheless, potential prognostic roles of CD8+ T cell-related genes (CD8Gs) in HCC remains unknown. In our study, 416 CD8Gs were identified in HCC, which were enriched in inflammatory and immune signaling pathways. Using The Cancer Genome Atlas dataset, a 5-CD8Gs risk model (KLRB1, FYN, IL2RG, FCER1G, and DGKZ) was constructed, which was verified in International Cancer Genome Consortium and gene expression omnibus datasets. Furthermore, we found that overall survival was independently correlated with the CD8Gs signature, and it was associated with immune- and cancer-related signaling pathways and immune cells infiltration. Finally, drug sensitivity data indicated that 10 chemotherapeutic drugs held promise as therapeutics for HCC patients with high-risk. In conclusion, multi-databases analysis showed that 5-CD8Gs and their signature could be an indicator to predict candidate drugs for HCC therapy.

Integrating single-cell and bulk expression data to identify and analyze cancer prognosis-related genes

Compared with traditional evaluation methods of cancer prognosis based on tissue samples, single-cell sequencing technology can provide information on cell type heterogeneity for predicting biomarkers related to cancer prognosis. Therefore, the bulk and single-cell expression profiles of breast cancer and normal cells were comprehensively analyzed to identify malignant and non-malignant markers and construct a reliable prognosis model. We first screened highly reliable differentially expressed genes from bulk expression profiles of multiple breast cancer tissues and normal tissues, and inferred genes related to cell malignancy from single-cell data. Then we identified eight critical genes related to breast cancer to conduct Cox regression analysis, calculate polygenic risk score (PRS), and verify the predictive ability of PRS in two data groups. The results show that PRS can divide breast cancer patients into high-risk group and low-risk group. PRS is related to the overall survival time and relapse-free interval and is a prognosis factor independent of conventional clinicopathological characteristics. Breast cancer is usually regarded as a cancer with a relatively good prognosis. In order to further explore whether this workflow can be applied to cancer with poor prognosis, we selected lung cancer for a comparative study. The results show that this workflow can also build a reasonable prognosis model for lung cancer. This study provides new insight and practical source code for further research on cancer biomarkers and drug targets. It also provides basis for survival prediction, treatment response prediction, and personalized treatment.

Generalized machine learning based on multi-omics data to profile the effect of ferroptosis pathway on prognosis and immunotherapy response in patients with bladder cancer

INTRODUCTION

Bladder cancer (BLCA) affects millions of people worldwide, with high rates of incidence and mortality. Ferroptosis proves to be a novel form of cell death process that is triggered by oxidative stress.

METHODS

We procured a total of 25 single nuclear RNA-seq (snRNA-seq) samples from GSE169379 in GEO database. We obtained different cohorts of BLCA patients from the TCGA and GEO databases for model training and validation. A total of 369 ferroptosis-related genes (FRGs) were selected from the FerrDb database. AUCell analysis was performed to assign ferroptosis scores to all the cell types. Weighted Gene Co-Expression Network Analysis (WGCNA), COX, and LASSO regression analysis were conducted to retain and finalize the genes of prognostic values. Various bioinformatic approaches were utilized to depict immune infiltration profile. We conducted a series of colony formation analysis, flow cytometry and western blot (WB) analysis to determine the role of SKAP1 in BLCA.

RESULTS

We divided the cells into high ferroptosis group and low ferroptosis group according to ferroptosis activity score, and then screened 2150 genes most associated with ferroptosis by differential expression analysis, which are related to UV-induced DNA damage, male hormone response, fatty acid metabolism and hypoxia. Subsequently, WGCNA algorithm further screened 741 ferroptosis related genes from the 2150 genes for the construction of prognostic model. Lasso-Cox regression analysis was used to construct the prognostic model, and the prognostic model consisting of 6 genes was obtained, namely JUN, SYT1, MAP3K8, GALNT14, TCIRG1, and SKAP1. Next, we constructed a nomogram model that integrated clinical factors to improving the accuracy. In addition, we performed drug sensitivity analyses in different subgroups and found that Staurosporine, Rapamycin, Gemcitabine, and BI-2536 may be candidates for the drugs treatment in high-risk populations. The ESTIMATE results showed higher stromal scores, immune scores, and ESTIMATE scores in the low-risk group, indicating a higher overall immunity level and immunogenicity of tumor microenvironment (TME) in this group, and tumor immune dysfunction and exclusion (TIDE) analysis confirmed a better response to immunotherapy in the low-risk group. Finally, we selected the oncogene SKAP1 in the prognostic gene for in vitro validation, and found that SKAP1 directly regulated BLCA cell proliferation and apoptosis.

CONCLUSION

We identified a set of six genes, JUN, SYT1, MAP3K8, GALNT14, TCIRG1, and SKAP1, that exhibited significant potential in stratification of BLCA patients with varying prognosis. In addition, we uncovered the direct regulatory effect of SKAP1 on BLCA cell proliferation and apoptosis, shedding some light on the role of FRGs in pathogenesis of BLCA.

Association between spatial distribution of leukocyte subsets and clinical presentation of head and neck squamous cell carcinoma

Background

Interactions between tumor cells and cells in the microenvironment contribute to tumor development and metastasis. The spatial arrangement of individual cells in relation to each other influences the likelihood of whether and how these cells interact with each other.

Methods

This study investigated the effect of spatial distribution on the function of leukocyte subsets in the microenvironment of human head and neck squamous cell carcinoma (HNSCC) using multiplex immunohistochemistry (IHC). Leukocyte subsets were further classified based on analysis of two previously published HNSCC single-cell RNA datasets and flow cytometry (FC).

Results

IHC revealed distinct distribution patterns of leukocytes differentiated by CD68 and CD163. While CD68hiCD163lo and CD68hiCD163hi cells accumulated near tumor sites, CD68loCD163hi cells were more evenly distributed in the tumor stroma. PD-L1hi and PD-1hi cells accumulated predominantly around tumor sites. High cell density of PD-L1hi CD68hiCD163hi cells or PD-1hi T cells near the tumor site correlated with improved survival. FC and single cell RNA revealed high variability within the CD68/CD163 subsets. CD68hiCD163lo and CD68hiCD163hi cells were predominantly macrophages (MΦ), whereas CD68loCD163hi cells appeared to be predominantly dendritic cells (DCs). Differentiation based on CD64, CD80, CD163, and CD206 revealed that TAM in HNSCC occupy a broad spectrum within the classical M1/M2 polarization. Notably, the MΦ subsets expressed predominantly CD206 and little CD80. The opposite was observed in the DC subsets.

Conclusion

The distribution patterns and their distinct interactions via the PD-L1/PD-1 pathway suggest divergent roles of CD68/CD163 subsets in the HNSCC microenvironment. PD-L1/PD-1 interactions appear to occur primarily between specific cell types close to the tumor site. Whether PD-L1/PD-1 interactions have a positive or negative impact on patient survival appears to depend on both the spatial localization and the entity of the interacting cells. Co-expression of other markers, particularly CD80 and CD206, supports the hypothesis that CD68/CD163 IHC subsets have distinct functions. These results highlight the association between spatial leukocyte distribution patterns and the clinical presentation of HNSCC.

A novel oxidative stress-related gene signature as an indicator of prognosis and immunotherapy responses in HNSCC

PURPOSE

To identify molecular subtypes of oxidative stress-related genes in head and neck squamous cell carcinoma (HNSCC) and to construct a scoring model of oxidative stress-related genes.

METHODS

R language based scRNA-seq and bulk RNA-seq analyses were used to identify molecular isoforms of oxidative stress-related genes in HNSCC. An oxidative stress-related gene scoring (OSRS) model was constructed, which were verified through online data and immunohistochemical staining of clinical samples.

RESULTS

Using TCGA-HNSCC datasets, nine predictive genes for overall patient survival, rarely reported in previous similar studies, were screened. AREG and CES1 were identified as prognostic risk factors. CSTA, FDCSP, JCHAIN, IFFO2, PGLYRP4, SPOCK2 and SPINK6 were identified as prognostic factors. Collectively, all genes formed a prognostic risk signature model for oxidative stress in HNSCC, which were validated in GSE41613, GSE103322 and PRJEB23709 datasets. Immunohistochemical staining of SPINK6 in nasopharyngeal cancer samples validated the gene panel. Subsequent analysis indicated that subgroups of the oxidative stress prognostic signature played important roles during cellular communication, the immune microenvironment, the differential activation of transcription factors, oxidative stress and immunotherapeutic responses.

CONCLUSIONS

The risk model might predict HNSCC prognosis and immunotherapeutic responses.

Elevated FBXL6 activates both wild-type KRAS and mutant KRASG12D and drives HCC tumorigenesis via the ERK/mTOR/PRELID2/ROS axis in mice

BACKGROUND

Kirsten rat sarcoma (KRAS) and mutant KRASG12D have been implicated in human cancers, but it remains unclear whether their activation requires ubiquitination. This study aimed to investigate whether and how F-box and leucine-rich repeat 6 (FBXL6) regulates KRAS and KRASG12D activity in hepatocellular carcinoma (HCC).

METHODS

We constructed transgenic mouse strains LC (LSL-Fbxl6KI/+;Alb-Cre, n = 13), KC (LSL-KrasG12D/+;Alb-Cre, n = 10) and KLC (LSL-KrasG12D/+;LSL-Fbxl6KI/+;Alb-Cre, n = 12) mice, and then monitored HCC for 320 d. Multiomics approaches and pharmacological inhibitors were used to determine oncogenic signaling in the context of elevated FBXL6 and KRAS activation. Co‑immunoprecipitation (Co-IP), Western blotting, ubiquitination assay and RAS activity detection assay were employed to investigate the underlying molecular mechanism by which FBXL6 activates KRAS. The pathological relevance of the FBXL6/KRAS/extracellular signal-regulated kinase (ERK)/mammalian target of rapamycin (mTOR)/proteins of relevant evolutionary and lymphoid interest domain 2 (PRELID2) axis was evaluated in 129 paired samples from HCC patients.

RESULTS

FBXL6 is highly expressed in HCC as well as other human cancers (P < 0.001). Interestingly, FBXL6 drives HCC in transgenic mice. Mechanistically, elevated FBXL6 promotes the polyubiquitination of both wild-type KRAS and KRASG12D at lysine 128, leading to the activation of both KRAS and KRASG12D and promoting their binding to the serine/threonine-protein kinase RAF, which is followed by the activation of mitogen-activated protein kinase kinase (MEK)/ERK/mTOR signaling. The oncogenic activity of the MEK/ERK/mTOR axis relies on PRELID2, which induces reactive oxygen species (ROS) generation. Furthermore, hepatic FBXL6 upregulation facilitates KRASG12D to induce more severe hepatocarcinogenesis and lung metastasis via the MEK/ERK/mTOR/PRELID2/ROS axis. Dual inhibition of MEK and mTOR effectively suppresses tumor growth and metastasis in this subtype of cancer in vivo. In clinical samples, FBXL6 expression positively correlates with p-ERK (χ2 = 85.067, P < 0.001), p-mTOR (χ2 = 66.919, P < 0.001) and PRELID2 (χ2 = 20.891, P < 0.001). The Kaplan-Meier survival analyses suggested that HCC patients with high FBXL6/p-ERK levels predicted worse overall survival (log‑rank P < 0.001).

CONCLUSIONS

FBXL6 activates KRAS or KRASG12D via ubiquitination at the site K128, leading to activation of the ERK/mTOR/PRELID2/ROS axis and tumorigenesis. Dual inhibition of MEK and mTOR effectively protects against FBXL6- and KRASG12D-induced tumorigenesis, providing a potential therapeutic strategy to treat this aggressive subtype of liver cancer.

Synaptotagmin 1 Suppresses Colorectal Cancer Metastasis by Inhibiting ERK/MAPK Signaling-Mediated Tumor Cell Pseudopodial Formation and Migration

Although synaptotagmin 1 (SYT1) has been identified participating in a variety of cancers, its role in colorectal cancer (CRC) remains an enigma. This study aimed to demonstrate the effect of SYT1 on CRC metastasis and the underlying mechanism. We first found that SYT1 expressions in CRC tissues were lower than in normal colorectal tissues from the CRC database and collected CRC patients. In addition to this, SYT1 expression was also lower in CRC cell lines than in the normal colorectal cell line. SYT1 expression was downregulated by TGF-β (an EMT mediator) in CRC cell lines. In vitro, SYT1 overexpression repressed pseudopodial formation and reduced cell migration and invasion of CRC cells. SYT1 overexpression also suppressed CRC metastasis in tumor-bearing nude mice in vivo. Moreover, SYT1 overexpression promoted the dephosphorylation of ERK1/2 and downregulated the expressions of Slug and Vimentin, two proteins tightly associated with EMT in tumor metastasis. In conclusion, SYT1 expression is downregulated in CRC. Overexpression of SYT1 suppresses CRC cell migration, invasion, and metastasis by inhibiting ERK/MAPK signaling-mediated CRC cell pseudopodial formation. The study suggests that SYT1 is a suppressor of CRC and may have the potential to be a therapeutic target for CRC.

Mitochondrial-related genes markers that predict survival in patients with head and neck squamous cell carcinoma affect immunomodulation through hypoxia, glycolysis, and angiogenesis pathways

Mitochondria play a crucial role in the occurrence and development of tumors. We used mitochondria-related genes for consistent clustering to identify three stable molecular subtypes of head and neck squamous cell carcinoma (HNSCC) with different prognoses, mutations, and immune characteristics. Significant differences were observed in clinical characteristics, immune microenvironment, immune cell infiltration, and immune cell scores. TP53 was the most significantly mutated; cell cycle-related pathways and tumorigenesis-related pathways were activated in different subtypes. Risk modeling was conducted using a multifactor stepwise regression method, and nine genes were identified as mitochondria-related genes affecting prognosis (DKK1, EFNB2, ITGA5, AREG, EPHX3, CHGB, P4HA1, CCND1, and JCHAIN). Risk score calculations revealed significant differences in prognosis, immune cell scores, immune cell infiltration, and responses to conventional chemotherapy drugs. Glycolysis, angiogenesis, hypoxia, and tumor-related pathways were positively correlated with the RiskScore. Clinical samples were subjected to qPCR to validate the results. In this work, we constructed a prognostic model based on the mitochondrial correlation score, which well reflects the risk and positive factors for the prognosis of patients with HNSCC. This model can be used to guide individualized adjuvant and immunotherapy in patients with HNSCC.

The prognostic significance of KLRB1 and its further association with immune cells in breast cancer

Background

Killer cell lectin-like receptor B1 (KLRB1) is an important member of the natural killer cell gene family. This study explored the potential value of KLRB1 as a breast cancer (BC) biomarker and its close association with the tumor immune microenvironment during the development of BC.

Methods

We examined the differential expression of KLRB1 in pan-cancer. Clinical and RNA-Seq data from BC samples were evaluated in The Cancer Genome Atlas (TCGA) and validated in Gene Expression Omnibus (GEO) datasets and by immunohistochemistry (IHC) staining. The relationship between KLRB1 and clinical parameters was explored through Chi-square tests. The diagnostic value of KLRB1 was evaluated using a receiver operating characteristic (ROC) curve. Survival analysis was tested by Kaplan-Meier curves to demonstrate the relationship between KLRB1 and survival. Univariable and multivariate cox regression analyses were carried out as well. The analysis of immune infiltration level and gene set enrichment analysis (GSEA) were conducted to examine KLRB1's mechanism during the progression of BC. We used the Tumor Immune Estimation Resource (TIMER), the Cancer Single-cell Expression Map (CancerSCEM) database, the Tumor Immune Single-cell Hub (TISCH) database, and the Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) method to explore KLRB1's association with immune infiltration level and different quantitative distribution of immune cells. The relevant signaling pathways in BC associated with KLRB1 were identified using GSEA.

Results

The expression of KLRB1 was downregulated across the majority of cancers including BC. The lower KLRB1 expression group exhibited shorter relapse free survival (RFS) and overall survival (OS). IHC staining showed that KLRB1 staining was weaker in breast tumor tissues than in paratumors. Additionally, GSEA identified several pathway items distinctly enriched in BC. KLRB1 expression level was also positively related to the infiltrating number of immune cells in BC. Moreover, the CancerSCEM and TISCH databases as well as the CIBERSORT method demonstrated the close relationship between KLRB1 and immune cells, particularly macrophages.

Conclusion

Low KLRB1 expression was considered an independent prognostic biomarker and played an important role in the tumor immune microenvironment of BC patients.

CircTTLL13 Promotes TMZ Resistance in Glioma via Modulating OLR1-Mediated Activation of the Wnt/β-Catenin Pathway

Glioma, originating from neuroglial progenitor cells, is a type of intrinsic brain tumor with poor prognosis. temozolomide (TMZ) is the first-line chemotherapeutic agent for glioma. Exploring the mechanisms of circTTLL13 underlying TMZ resistance in glioma is of great significance to improve glioma treatment. Bioinformatics was adopted to identify target genes. The circular structure of circTTLL13 and its high expression in glioma cells were disclosed by quantitative real time-PCR (qRT-PCR) and PCR-agarose gel electrophoresis. Functional experiments proved that oxidized LDL receptor 1 (OLR1) promotes TMZ resistance of glioma cells. CircTTLL13 enhances TMZ resistance of glioma cells via modulating OLR1. Luciferase reporter, RNA-binding protein immunoprecipitation (RIP), RNA pulldown, mRNA stability, N6-methyladenosine (m6A) dot blot and RNA total m6A quantification assays were implemented, indicating that circTTLL13 stabilizes OLR1 mRNA via recruiting YTH N6-methyladenosine RNA binding protein 1 (YTHDF1) and promotes m6A methylation of OLR1 pre-mRNA through recruiting methyltransferase-like 3 (METTL3). TOP/FOP-flash reporter assay and western blot verified that circTTLL13 activates Wnt/β-catenin signaling pathway by regulating OLR1. CircTTLL13 promotes TMZ resistance in glioma through regulating OLR1-mediated Wnt/β-catenin pathway activation. This study offers an insight into the efficacy improvement of TMZ for glioma treatment.

Construction of diagnostic and prognostic models based on gene signatures of nasopharyngeal carcinoma by machine learning methods

Background

Diagnostic models based on gene signatures of nasopharyngeal carcinoma (NPC) were constructed by random forest (RF) and artificial neural network (ANN) algorithms. Least absolute shrinkage and selection operator (Lasso)-Cox regression was used to select and build prognostic models based on gene signatures. This study contributes to the early diagnosis and treatment, prognosis, and molecular mechanisms associated with NPC.

Methods

Two gene expression datasets were downloaded from the Gene Expression Omnibus (GEO) database, and differentially expressed genes (DEGs) associated with NPC were identified by gene expression differential analysis. Subsequently, significant DEGs were identified by a RF algorithm. ANN were used to construct a diagnostic model for NPC. The performance of the diagnostic model was evaluated by area under the curve (AUC) values using a validation set. Lasso-Cox regression examined gene signatures associated with prognosis. Overall survival (OS) and disease-free survival (DFS) prediction models were constructed and validated from The Cancer Genome Atlas (TCGA) database and the International Cancer Genome Consortium (ICGC) database.

Results

A total of 582 DEGs associated with NPC were identified, and 14 significant genes were identified by the RF algorithm. A diagnostic model for NPC was successfully constructed using ANN, and the validity of the model was confirmed on the training set AUC =0.947 [95% confidence interval (CI): 0.911-0.969] and the validation set AUC =0.864 (95% CI: 0.828-0.901). The 24-gene signatures associated with prognosis were identified by Lasso-Cox regression, and prediction models for OS and DFS of NPC were constructed on the training set. Finally, the ability of the model was validated on the validation set.

Conclusions

Several potential gene signatures associated with NPC were identified, and a high-performance predictive model for early diagnosis of NPC and a prognostic prediction model with robust performance were successfully developed. The results of this study provide valuable references for early diagnosis, screening, treatment and molecular mechanism research of NPC in the future.

A multi-omics-based investigation of the prognostic and immunological impact of necroptosis-related mRNA in patients with cervical squamous carcinoma and adenocarcinoma

Necroptosis is a kind of programmed necrosis mode that plays a double-edged role in tumor progression. However, the role of necroptosis-related Messenger RNA (mRNA) in predicting the prognosis and immune response of cervical squamous carcinoma and adenocarcinoma (CESC) has not been fully studied. Firstly, the incidence of somatic mutation rate and copy number variation for 74 necroptosis-related mRNAs (NRmRNAs) were analyzed. Secondly, CESC patients were divided into four stable clusters based on the consensus clustering results and analyzed for correlations with a series of clinical factors. Subsequently, a total of 291 The Cancer Genome Atlas samples were randomly divided into either training or validation cohorts. A Cox proportional hazard model consisting of three NRmRNAs (CXCL8, CLEC9A, and TAB2) was constructed by univariate, least absolute shrinkage and selection operator and multivariate COX regression analysis to identify the prognosis and immune response. Its performance and stability were further validated in another testing dataset (GSE44001) from Gene Expression Omnibus database. The results of the receiver operating characteristic curve, principal component analysis, t-SNE, and nomogram indicated that the prognostic model we constructed can serve as an independent prognostic factor. The combination of the prognostic model and the classic TNM staging system could improve the performance in predicting the survival of CESC patients. In addition, differentially expressed genes from high and low-risk patients are screened by R software for functional analysis and pathway enrichment analysis. Besides, single-sample gene set enrichment analysis revealed that tumor-killing immune cells were reduced in the high-risk group. Moreover, patients in the low-risk group are more likely to benefit from immune checkpoint inhibitors. The analysis of tumor immune dysfunction and exclusion scores, M6A-related genes, stem cell correlation and Tumor mutational burden data with clinical information has quantified the expression levels of NRmRNAs between the two risk subgroups. According to tumor immune microenvironment scores, Spearman’s correlation analysis, and drug sensitivity, immunotherapy may have a higher response rate and better efficacy in patients of the low-risk subgroup. In conclusion, we have reported the clinical significance of NRmRNAs for the prognosis and immune response in CESC patients for the first time. Screening of accurate and effective prognostic markers is important for designing a multi-combined targeted therapeutic strategy and the development of individualized precision medicine.

Single-Cell RNA Sequencing Reveals the Role of Epithelial Cell Marker Genes in Predicting the Prognosis of Colorectal Cancer Patients

Single-cell RNA sequencing (scRNA-seq) is increasingly used in studies on gastrointestinal cancers. This study investigated the prognostic value of epithelial cell-associated biomarkers in colorectal cancer (CRC) using scRNA-seq data. We downloaded and analysed scRNA-seq data from four CRC samples from the Gene Expression Omnibus (GEO), and we identified marker genes of malignant epithelial cells (MECs) using CRC transcriptome and clinical data downloaded from The Cancer Genome Atlas (TCGA) and GEO as training and validation cohorts, respectively. In the TCGA training cohort, weighted gene correlation network analysis, univariate Cox proportional hazard model (Cox) analysis, and least absolute shrinkage and selection operator regression analysis were performed on the marker genes of MEC subsets to identify a signature of nine prognostic MEC-related genes (MECRGs) and calculate a risk score based on the signature. CRC patients were divided into high- and low-risk groups according to the median risk score. We found that the MECRG risk score was significantly correlated with the clinical features and overall survival of CRC patients, and that CRC patients in the high-risk group showed a significantly shorter survival time. The univariate and multivariate Cox regression analyses showed that the MECRG risk score can serve as an independent prognostic factor for CRC patients. Gene set enrichment analysis revealed that the MECRG signature genes are involved in fatty acid metabolism, p53 signalling, and other pathways. To increase the clinical application value, we constructed a MECRG nomogram by combining the MECRG risk score with other independent prognostic factors. The validity of the nomogram is based on receiver operating characteristics and calibration curves. The MECRG signature and nomogram models were well validated in the GEO dataset. In conclusion, we established an epithelial cell marker gene-based risk assessment model based on scRNA-seq analysis of CRC samples for predicting the prognosis of CRC patients.

Genomic and Immunological Characterization of Pyroptosis in Lung Adenocarcinoma

Pyroptosis is a programmed cell death that may either promote or hinder cancer growth under different circumstances. Pyroptosis-related genes (PRGs) could be a useful target for cancer therapy, and are uncommon in lung adenocarcinoma (LUAD). The expression profiles, mutation data and clinical information of LUAD patients were included in this study. A pyroptosis-related prognostic risk score (PPRS) model was constructed by performing Cox regression, weighted gene co-expression network analysis (WGCNA), and least absolute shrinkage and selection operator (LASSO) analysis to score LUAD patients. Somatic mutation and copy number variation (CNV), tumor immunity, and sensitivity to immunotherapy/chemotherapy were compared between different PPRS groups. Clinical parameters of LUAD were combined with PPRS to construct a decision tree and nomogram. Red module was highly positively correlated with pyroptosis. Seven genes (FCRLB, COTL1, GNG10, CASP4, DOK1, CCR2, and AQP8) were screened from the red module to construct a PPRS model. Significantly lower overall survival (OS), higher incidence of somatic mutation and CNV, elevated infiltration level of the immune cell together with increased probability of immune escape were observed in LUAD patients with higher PPRS, and were more sensitive to Cisplatin, Docetaxel, and Vinorelbine. We constructed a new PPRS model for patients with LUAD. The model might have clinical significance in the prediction of the prognosis of patients with LUAD and in the efficacy of chemotherapy and immunotherapy.

Identification of the circRNA-miRNA-mRNA Prognostic Regulatory Network in Lung Adenocarcinoma

Background: Numerous studies have identified that circular RNAs (circRNAs) can serve as competing endogenous RNAs (ceRNAs) to regulate tumor progression. However, there are still a large number of circRNAs to be deciphered. Objective: The purpose of this study was to reveal novel circRNAs and their potential role in lung adenocarcinoma (LUAD). Methods: To unveil LUAD-related circRNAs, microRNA (miRNAs), and messenger RNA (mRNA) and elucidate their possible molecular mechanisms, we employed a strategy combining extensive data mining and bioinformatics methods. According to the results of bioinformatics workflow analysis, a novel circRNA-miRNA-mRNA network was constructed. Results: Ten circRNAs with different expressions were acquired from four Gene Expression Omnibus (GEO) microarray datasets. Seven Prognostic-related differential miRNAs of LUAD were gained from The Cancer Genome Atlas (TCGA). Simultaneously, the miRNA reaction components corresponding to the ten circRNAs were predicted. Two circRNA–miRNA interactions including two circRNAs (hsa_circ_0008234 and hsa_circ_0002360) and two miRNAs (hsa-miR-490-3p and hsa-miR-1293) were identified above. Then, target genes of the two miRNAs and differently expressed genes (DEGs) from TCGA on LUAD were collected. Three hub-genes (ADCY9, NMUR1, SYT1) were determined according to prognosis in patients with LUAD ulteriorly. Conclusions: hsa_circ_0008234/hsa-miR-490-3p/SYT1 and hsa_circ_0002360/hsa-miR-1293/ (ADCY9, NMUR1) networks were established, and identified molecules may be involved in pathogenesis and prognosis in patients with LUAD.

A web-based calculator for predicting the prognosis of patients with sarcoma on the basis of antioxidant gene signatures

Background: Oxidative stress plays a critical role in tumorigenesis, tumor development, and resistance to therapy. A systematic analysis of the interactions between antioxidant gene expression and the prognosis of patients with sarcoma is lacking but urgently needed.

Methods: Gene expression and clinical data of patients with sarcoma were derived from The Cancer Genome Atlas Sarcoma (training cohort) and Gene Expression Omnibus (validation cohorts) databases. Least absolute shrinkage, selection operator regression, and Cox regression were used to develop prognostic signatures for overall survival (OS) and disease-free survival (DFS). Based on the signatures and clinical features, two nomograms for predicting 2-, 4-, and 6-year OS and DFS were established.

Results: On the basis of the training cohort, we identified five-gene (CHAC2, GPX5, GSTK1, PXDN, and S100A9) and six-gene (GGTLC2, GLO1, GPX7, GSTK1, GSTM5, and IPCEF1) signatures for predicting OS and DFS, respectively, in patients with sarcoma. Kaplan–Meier survival analysis of the training and validation cohorts revealed that patients in the high-risk group had a significantly poorer prognosis than those in the low-risk group. On the basis of the signatures and other independent risk factors, we established two models for predicting OS and DFS that showed excellent calibration and discrimination. For the convenience of clinical application, we built web-based calculators (OS: https://quankun.shinyapps.io/sarcOS/; DFS: https://quankun.shinyapps.io/sarcDFS/).

Conclusions: The antioxidant gene signature models established in this study can be novel prognostic predictors for sarcoma.

Differential gene expression and network analysis in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a prevalent malignancy with a poor prognosis, whose biomarkers have not been studied in great detail. We have collected genomic data of HNSCC patients from The Cancer Genome Atlas (TCGA) and analyzed them to get deeper insights into the gene expression pattern. Initially, 793 differentially expressed genes (DEGs) were categorized, and their enrichment analysis was performed. Later, a protein-protein interaction network for the DEGs was constructed using the STRING plugin in Cytoscape to study their interactions. A set of 10 hub genes was selected based on Maximal Clique Centrality score, and later their survival analysis was studied. The elucidated set of 10 genes, i.e., PRAME, MAGEC2, MAGEA12, LHX1, MAGEA3, CSAG1, MAGEA6, LCE6A, LCE2D, LCE2C, referred to as potential candidates to be explored as HNSCC biomarkers. The Kaplan-Meier overall survival of the selected genes suggested that the alterations in the candidate genes were linked to the decreased survival of the HNSCC patients. Altogether, the results of this study signify that the genomic alterations and differential expression of the selected genes can be explored in therapeutic interpolations of HNSCC, exploiting early diagnosis and target-propelled therapy.

Radiotherapy resistance: identifying universal biomarkers for various human cancers

Radiotherapy (RT) is considered as a standard in the treatment of most solid cancers, including glioblastoma, lung, breast, rectal, prostate, colorectal, cervical, esophageal, and head and neck cancers. The main challenge in RT is tumor cell radioresistance associated with a high risk of locoregional relapse and distant metastasis. Despite significant progress in understanding mechanisms of radioresistance, its prediction and overcoming remain unresolved. This review presents the state-of-the-art for the potential universal biomarkers correlated to the radioresistance and poor outcome in different cancers. We describe radioresistance biomarkers functionally attributed to DNA repair, signal transduction, hypoxia, and angiogenesis. We also focus on high throughput genetic and proteomic studies, which revealed a set of molecular biomarkers related to radioresistance. In conclusion, we discuss biomarkers which are overlapped in most several cancers.

Application of liquid biopsy as multi-functional biomarkers in head and neck cancer

Head and neck squamous cell carcinoma (HNSCC) is a molecularly heterogeneous disease, with a 5-year survival rate that still hovers at ~60% despite recent advancements. The advanced stage upon diagnosis, limited success with effective targeted therapy and lack of reliable biomarkers are among the key factors underlying the marginally improved survival rates over the decades. Prevention, early detection and biomarker-driven treatment adaptation are crucial for timely interventions and improved clinical outcomes. Liquid biopsy, analysis of tumour-specific biomarkers circulating in bodily fluids, is a rapidly evolving field that may play a striking role in optimising patient care. In recent years, significant progress has been made towards advancing liquid biopsies for non-invasive early cancer detection, prognosis, treatment adaptation, monitoring of residual disease and surveillance of recurrence. While these emerging technologies have immense potential to improve patient survival, numerous methodological and biological limitations must be overcome before their implementation into clinical practice. This review outlines the current state of knowledge on various types of liquid biopsies in HNSCC, and their potential applications for diagnosis, prognosis, grading treatment response and post-treatment surveillance. It also discusses challenges associated with the clinical applicability of liquid biopsies and prospects of the optimised approaches in the management of HNSCC.

The Highly Expressed IFIT1 in Nasopharyngeal Carcinoma Enhances Proliferation, Migration, and Invasion of Nasopharyngeal Carcinoma Cells

In this study, we aimed to identify potential targets modulating the progression of nasopharyngeal carcinoma (NPC) using integrated bioinformatics analysis and functional assays. Differentially expressed genes (DEGs) between NPC and normal tissues samples were obtained from publicly availably microarray datasets (GSE68799, GSE34573, and GSE53819) in the Gene Expression Omnibus (GEO) database. The bioinformatics analysis identified 49 common DEGs from three GEO datasets, which were mainly enriched in cytokine/chemokine pathways and extracellular matrix organization pathway. Further protein-protein interaction network analysis identified 11 hub genes from the 49 DEGs. The 11 hub genes were significantly up-regulated in the NPC tissues when compared to normal tissues by analyzing the Oncomine database. The 8 hub genes including COL5A1, COL7A1, COL22A1, CXCL11, IFI44L, IFIT1, RSAD2, and USP18 were significantly up-regulated in the NPC tissues when compared to normal tissues by using the Oncomine database. Further validation studies showed that IFIT1 was up-regulated in the NPC cells. Knockdown of IFI1T1 suppressed the proliferation, migration, and invasion of NPC cells; while IFIT1 overexpression promoted the proliferation, migration, and invasion of NPC cells. In conclusion, a total of 49 DEGs and 11 hub genes in NPC using the integrated bioinformatics analysis. IFIT1 was up-regulated in the NPC cells lines, and IFIT1 may act as an oncogene by promoting NPC cell proliferation, migration, and invasion.

Knockdown of lncRNA CCAT1 Inhibits the Progression of Colorectal Cancer via hsa-miR-4679 Mediating the Downregulation of GNG10

Great concerns have raised crucial roles of long noncoding RNAs (lncRNAs) on colorectal cancer progression due to the increasing number of studies in cancer development. Previous studies reveal that lncRNA CCAT1 plays an important role in the progression of a variety of cancers. However, the role of lncRNA CCAT1 in colorectal cancer is still unclear. In this study, we found that in both colorectal tissues and cell lines the level of lncRNA CCAT1 was increased. Downregulation of lncRNA CCAT1 inhibited the proliferation, migration, and invasion of colorectal cell lines and promoted apoptosis. We then found that hsa-miR-4679 could bind to lncRNA CCAT1 directly, and with further functional analyses, we confirmed that lncRNA CCAT1 sponged hsa-miR-4679 to promote the progression of colorectal cancer. Next, we found that hsa-miR-4679 was directly bound to 3'UTR of GNG10 (guanine nucleotide-binding protein, gamma 10). GNG10 overexpression promoted the progression of colorectal cancer, and this phenotype could be reversed by miR-4679 mimics. At last, we knocked down CCAT1 in vivo and found that sh-CCAT1 reduced the tumor size and the number of proliferating cells. In summary, our findings revealed that lncRNA CCAT1 facilitated colorectal cancer progression via the hsa-miR-4679/GNG10 axis and provided new potential therapeutic targets for colorectal cancer.

Integrative analysis-based identification and validation of a prognostic immune cell infiltration-based model for patients with advanced gastric cancer

BACKGROUNDS

Advanced gastric cancer (GC) remains difficult to conduct individualized prognostic evaluations owing to the highly heterogeneous nature and the low level of immune cell infiltration (ICI) within GC tumors. This study thus sought to develop a model capable of classifying GC patients according to the degree of tumor ICI and gauging prognosis.

METHODS

The degree of ICI in GC patients from the GSE15459, GSE57303, and GSE62254 datasets were estimated, and these values were used to group patients via an unsupervised clustering approach, after which ICI cluster-related genes were identified the association with prognosis through Cox and LASSO regression analyses. The primary risk genes were then verified by immunohistochemical staining of GC tumor tissue samples.

RESULTS

570 patients were clustered into three clusters and 289 ICI cluster-related genes were identified. A prognostic model based on the expression of six crucial ICI risk genes (CXCL11, RBPMS2, LOC400043, JCHAIN, CT83, and ORM1) wa constructed. Patients identified as being high risk based upon the model have poorer clinical features and survival outcomes compared to the other patients. Adjuvant intervention was found to be more beneficial for patients expressing high levels of RBPMS2, JCHAIN, or ORM1. Furthermore, patients expressing low levels of JCHAIN or CT83 in GC tumor tissues were verified to exhibit a significantly better prognosis in a CMU cohort.

CONCLUSION

Advanced GC patients were successfully grouped into clusters based on the degree of intratumoral ICI, and a prognostic evaluation model based on 6 ICI risk genes was developed and validated.

Cross-link between ferroptosis and nasopharyngeal carcinoma: New approach to radiotherapy sensitization

Ferroptosis is a recently discovered special type of regulated cell death that is strongly associated with both homeostasis maintenance and cancer development. Previous studies have indicated that a number of small-molecular agents inducing ferroptosis have great potential in the treatment of different types of cancer, including breast, pancreatic, prostate and head and neck cancer. However, the role of ferroptosis in nasopharyngeal carcinoma (NPC) has remained to be fully determined. To the best of our knowledge, no review of the currently available studies on this subject has been published to date. The metabolism and expression of specific genes that regulate ferroptosis may represent a promising radiosensitization target in cancer treatment. The aim of the present review was to describe the cross-link between ferroptosis and NPC and to discuss the potential value of regulators and the possible mechanism underlying the role of ferroptosis in the radiosensitization of NPC, in the hope that linking the mechanism of ferroptosis with the development of NPC will accelerate the development of novel ferroptosis-based targets and radiotherapy strategies in NPC.

Current Status and Future Perspectives about Molecular Biomarkers of Nasopharyngeal Carcinoma

Simple Summary

Nasopharyngeal carcinoma is a serious major public health problem in its endemic countries. Up to 80% of NPC patients with locally advanced disease or distant metastasis at diagnosis were associated with poor prognosis and with median survival less than 4 months. The mortality rate of NPC metastasis is up to 91%. To date, there is no available curative treatment or reliable early diagnosis or prognosis for NPC. Discovery and development of reliable early diagnosis and prognosis biomarkers for nasopharyngeal carcinoma are urgent needed. Hence, we have here listed the potential early diagnosis and prognosis biomarker candidates for nasopharyngeal carcinoma. This review will give an insight to readers on the progress of NPC biomarker discovery to date, as well as future prospective biomarker development and their translation to clinical use.

Abstract

Nasopharyngeal carcinoma (NPC) is an epithelial malignancy that shows a remarkable ethnic and geographical distribution. It is one of the major public health problems in some countries, especially Southern China and Southeast Asia, but rare in most Western countries. Multifactorial interactions such as Epstein–Barr virus infection, individual’s genetic susceptibility, as well as environmental and dietary factors may facilitate the pathogenesis of this malignancy. Late presentation and the complex nature of the disease have led it to become a major cause of mortality. Therefore, an effective, sensitive, and specific molecular biomarker is urgently needed for early disease diagnosis, prognosis, and prediction of metastasis and recurrence after treatment. In this review, we discuss the recent research status of potential biomarker discovery and the problems that need to be explored further for better NPC management. By studying the aberrant pattern of these candidate biomarkers that promote NPC development and progression, we are able to understand the complexity of this malignancy better, hence positing our stands better towards strategies that may provide a way forward to the discovery of more reliable and specific biomarkers for diagnosis and targeted therapeutic development.

Investigations on the Role of the MicroRNA-338-5p/Wnt Family Member 2B (WNT2B) Axis in Regulating the Pathogenesis of Nasopharyngeal Carcinoma (NPC)

Background

The involvement of microRNA-338-5p in modulating NPC pathogenesis is still largely unknown, and this study aimed to investigate this issue.

Methods

The expressions of cancer associated genes were determined by Real-Time qPCR and Western Blot, and cell apoptosis was determined by flow cytometer (FCM). CCK-8 assay and colony formation assay were respectively used to determine cell proliferation and colony formation abilities. Transwell assay was used to evaluate cell migration. The expression levels of Ki67 protein in mice tissues were measured by Immunohistochemistry (IHC) assay.

Results

The present study found that microRNA-338-5p suppressed NPC progression by degrading its downstream target, Wnt family member 2B (WNT2B). Specifically, microRNA-338-5p tended to be low-expressed in NPC tissues and cell lines, compared to the non-tumor nasopharyngeal mucosa tissues and normal nasopharyngeal cell line (NP69). Upregulation of microRNA-338-5p inhibited proliferation, mobility, and epithelial-mesenchymal transition (EMT) in NPC cells in vitro, while silencing of microRNA-338-5p had opposite effects. Consistently, microRNA-338-5p suppressed tumorigenesis of NPC cells in vivo. In addition, microRNA-338-5p targeted WNT2B for degradation and inhibition, and the inhibiting effects of microRNA-338-5p overexpression on NPC development were reversed by upregulating WNT2B.

Conclusions

Taken together, we concluded that microRNA-338-5p targeted WNT2B to hinder NPC development.

Subtype-dependent regulation of Gβγ signalling

G protein-coupled receptors (GPCRs) transmit information to the cell interior by transducing external signals to heterotrimeric G protein subunits, Gα and Gβγ subunits, localized on the inner leaflet of the plasma membrane. Though the initial focus was mainly on Gα-mediated events, Gβγ subunits were later identified as major contributors to GPCR-G protein signalling. A broad functional array of Gβγ signalling has recently been attributed to Gβ and Gγ subtype diversity, comprising 5 Gβ and 12 Gγ subtypes, respectively. In addition to displaying selectivity towards each other to form the Gβγ dimer, numerous studies have identified preferences of distinct Gβγ combinations for specific GPCRs, Gα subtypes and effector molecules. Importantly, Gβ and Gγ subtype-dependent regulation of downstream effectors, representing a diverse range of signalling pathways and physiological functions have been found. Here, we review the literature on the repercussions of Gβ and Gγ subtype diversity on direct and indirect regulation of GPCR/G protein signalling events and their physiological outcomes. Our discussion additionally provides perspective in understanding the intricacies underlying molecular regulation of subtype-specific roles of Gβγ signalling and associated diseases.

c-Myc-induced circ-NOTCH1 promotes aggressive phenotypes of nasopharyngeal carcinoma cells by regulating the miR-34c-5p/c-Myc axis

Nasopharyngeal carcinoma (NPC) is the subclass of head and neck cancer with the highest incidence among otolaryngology malignancies. A growing amount of evidence has proven that circular RNAs (circRNAs) play key roles in the progression of multiple cancers. It has been reported that circ-NOTCH1 is a novel circRNA and functions as an oncogene in gastric cancer, while the regulatory mechanism of circ-NOTCH1 in NPC remains unknown. In the present research, our findings revealed that circ-NOTCH1 was overexpressed in NPC tissues and cells. Circ-NOTCH1 knockdown suppressed NPC cell proliferation, invasion, and migration. Subsequently, we discovered that c-Myc can activate circ-NOTCH1 by binding to the NOTCH1 promoter. c-Myc functioned as a tumor promoter in NPC cells. Mechanistically, circ-NOTCH1 served as a competitive endogenous RNA to modulate c-Myc expression by sponging miR-34c-5p. Additionally, overexpression of c-Myc reversed the circ-NOTCH1 knockdown-mediated inhibition of NPC cellular progression. Overall, this study suggested that c-Myc-induced circ-NOTCH1 promoted malignant phenotypes of NPC cells by regulating the miR-34c-5p/c-Myc axis.

Research status and prospects of biomarkers for nasopharyngeal carcinoma in the era of high‑throughput omics (Review)

As a malignant tumor type, nasopharyngeal carcinoma (NPC) is characterized by distinct geographical, ethnic and genetic differences; presenting a major threat to human health in many countries, especially in Southern China. At present, no accurate and effective methods are available for the early diagnosis, efficacious evaluation or prognosis prediction for NPC. As such, a large number of patients have locoregionally advanced NPC at the time of initial diagnosis. Many patients show toxic reactions to overtreatment and have risks of cancer recurrence and distant metastasis owing to insufficient treatment. To solve these clinical problems, high‑throughput '‑omics' technologies are being used to screen and identify specific molecular biomarkers for NPC. Because of the lack of comprehensive descriptions regarding NPC biomarkers, the present study summarized the research progress that has been made in recent years to discover NPC biomarkers, highlighting the existing problems that require exploration. In view of the lack of authoritative reports at present, study design factors that affect the screening of biomarkers are also discussed here and prospects for future research are proposed to provide references for follow‑up studies of NPC biomarkers.

A signature of estimate-stromal-immune score-based genes associated with the prognosis of lung adenocarcinoma

Background

Immune and stromal component evaluation is necessary to establish accurate prognostic markers for the prediction of clinical outcomes in lung adenocarcinoma (LUAD). We aimed to develop a gene signature based on the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE)-stromal-immune score in LUAD.

Methods

The transcriptomic profiles of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA), and the immune and stromal scores were derived using the ESTIMATE algorithm. The prognostic signature genes were selected from the differentially expressed genes (DEGs) using the robust partial likelihood-based cox proportional hazards regression method. The negative log-likelihood and the Akaike Information Criterion (AIC) were used to identify the optimal gene signature. The validation was carried out in 2 independent datasets from the Gene Expression Omnibus (GSE68571 and GSE72094).

Results

Patients with high ESTIMATE, stromal, and immune scores had better overall survivals (P=0.0035, P=0.066, and P=0.0077). The expression of thirty-seven genes was related to ESTIMATE-stromal-immune score. A risk stratification model was developed based on a gene signature containing CD74, JCHAIN, and PTGDS. The ESTIMATE-stromal-immune risk score was revealed to be a prognostic factor (P=0.009) after multivariate analysis. Four groups were classified based on this risk stratification model, yielding increasing survival outcomes (log-rank test, P=0.0051). This risk stratification model and other clinicopathological factors were combined to generate a nomogram. The calibration curves showed perfect agreement between the nomogram-predicted outcomes and those actually observed. Similar observations were made in 2 independent cohorts.

Conclusions

The gene signature based on the ESTIMATE-stromal-immune score could predict the prognosis of patients with LUAD.

Lipid metabolism-related proteins of relevant evolutionary and lymphoid interest (PRELI) domain containing family proteins in cancer

Metabolic reprogramming of tumor cells plays a critical role in the tumor microenvironment, including disorder of lipid metabolism. Recently, lipid metabolism has received increasing attention in cancer research. The proteins of relevant evolutionary and lymphoid interest (PRELI) domain containing family contains 6 proteins. Functionally, the PRELI-like family proteins were mainly involved in mitochondrial lipid transport and correlated with several types of diseases and malignant tumors. Here we review current knowledge of the functions, structures, biological functions and underlying mechanisms of the PRELI-like family proteins in cancer progression, which provide insights into the clinical translational application.

HPV infection related immune infiltration gene associated therapeutic strategy and clinical outcome in HNSCC

Background

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common tumor in human. Research has shown that HPV status HNSCC is a unique prognosis factor, which may due to its immune infiltration landscape. But the underlying mechanism is unclear.

Methods

In this study, we used a combination of several bioinformatics tools, including WCGNA, ssGSEA, CIBERSORT, TIDE,etc., to explore significant genes both related to HPV infection status and immune cell infiltration in HNSCC patients.

Results

Combined with several bioinformatics algorithms, eight hub genes were identified, including LTB, CD19, CD3D, SKAP1, KLRB1, CCL19, TBC1D10C and ARHGAP4. In HNSCC population, the hub genes had a stable co-expression, which was related to immune cell infiltration, especially CD8+ T cells, and the infiltrative immune cells were in a dysfunctional status. Samples with high hub genes expression presented with better response to immune check point block (ICB) therapy and sensitivity to bleomycin and methotrexate.

Conclusions

The eight hub genes we found presented with a stable co-expression in immune cell infiltration of HPV + ve HNSCC population. The co-expression of hub genes related to an immune microenvironment featuring an increase in immune cells but high degree of immune dysfunction status. Patients with high hub gene expression had a better response to ICB treatment, bleomycin and methotrexate. The co-expression of hub genes may be related to immune infiltration status in patients. The concrete molecular mechanism of hub genes function demands further exploration.

CLIC4 regulates radioresistance of nasopharyngeal carcinoma by iNOS after γ-rays but not carbon ions irradiation

Nasopharyngeal carcinoma (NPC) is a major health problem in the East and Southeast Asia, and the intensity modulated radiotherapy (IMRT) is the current preferred treatment method of NPC, but radioresistance-induced residual and recurrent tumors are the main cause of treatment failure. Till now, the mechanism of radioresistance and prognostic biomarkers of NPC are still unrevealed. In this study, we collected clinical NPC samples and established radioresistant NPC-R cell lines by irradiating NPC cells with fractionation doses of γ-rays. Using genechip assay between radioresistance and radiosensitive clinical samples and TMT assay between NPC and NPC-R cells, differential expressed genes were examined and the potential biomarker of radioresistance was screened. Immunohistochemical assay of NPC clinical specimens showed that CLIC4 was significantly up-regulated in radioresistance tumor tissues. In vitro studies confirmed that up-regulation of CLIC4 gene enhanced radioresistance in comparison with the alterations of intracellular oxidative metabolism of reactive oxygen species (ROS) and nitric oxide (NO) in an opposite way. Correspondingly, inhibition of CLIC4 sensitized NPC cells to irradiation and decreased nuclear translocation of iNOS and intracellular level of NO in NPC cells. Interestingly, the capacity for DNA repair had no difference between NPC and NPC-R cells. Moreover, because of great interests in using carbon ion irradiation to treat NPC effectively, we demonstrated that, after carbon ion irradiation, NPC-R and NPC cells had similar survival even under the status of up- or down-regulation of CLIC4. Conclusively, CLIC4 contributes to radioresistance of NPC to γ-rays but not carbon ions by regulating intracellular oxidative metabolism of nuclear translocation of iNOS.