Machine learning to construct sphingolipid metabolism genes signature to characterize the immune landscape and prognosis of patients with uveal melanoma

Crosstalk among disulfidptosis-related lncRNAs in lung adenocarcinoma reveals a correlation with immune profile and clinical prognosis

Disulfidptosis refers to a specific programmed cell death process characterized by the accumulation of disulfides. It has recently been reported in several cancers. However, the impact of disulfidptosis-related long non-coding RNAs (lncRNAs) on malignant tumors has remained largely unknown. In the present work, we screened prognostic disulfidptosis-related lncRNAs and studied their effects on lung adenocarcinoma. Relevant clinical data of lung adenocarcinoma cases were retrieved from The Cancer Genome Atlas (TCGA) database. RNA sequencing was used to identify differentially expressed disulfidptosis-related lncRNAs within lung adenocarcinoma. In addition, prognostic disulfidptosis-related lncRNAs were obtained through univariate Cox regression analysis. LASSO-COX was used to construct new disulfidptosis-related lncRNA signatures. Different statistical approaches were used to validate the practicability and accuracy of the disulfidptosis-related lncRNAs signatures. Furthermore, several bioinformatic approaches were used to study relevant heterogeneities in biological processes and pathways of diverse risk groups. Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) was conducted to analyze the expression of disulfidptosis-related lncRNAs. Finally, seven disulfidptosis-related lncRNA signatures were identified in lung adenocarcinoma cells. The prognosis prediction model constructed efficiently predicted patient survival. Subgroup analysis revealed significant differences in immune cell proportion, including T follicular helper cells and M0 macrophages. In addition, in vitro experimental results demonstrated significant differences in disulfidptosis-related lncRNAs. Altogether, the six disulfidptosis-related lncRNA signatures could serve as a potential prognostic biomarker for lung adenocarcinoma. Furthermore, these can be used as a prediction model in individualized immunotherapy for lung adenocarcinoma.

Characterizing mitochondrial features in osteoarthritis through integrative multi-omics and machine learning analysis

Purpose

Osteoarthritis (OA) stands as the most prevalent joint disorder. Mitochondrial dysfunction has been linked to the pathogenesis of OA. The main goal of this study is to uncover the pivotal role of mitochondria in the mechanisms driving OA development.

Materials and methods

We acquired seven bulk RNA-seq datasets from the Gene Expression Omnibus (GEO) database and examined the expression levels of differentially expressed genes related to mitochondria in OA. We utilized single-sample gene set enrichment analysis (ssGSEA), gene set enrichment analysis (GSEA), and weighted gene co-expression network analysis (WGCNA) analyses to explore the functional mechanisms associated with these genes. Seven machine learning algorithms were utilized to identify hub mitochondria-related genes and develop a predictive model. Further analyses included pathway enrichment, immune infiltration, gene-disease relationships, and mRNA-miRNA network construction based on these hub mitochondria-related genes. genome-wide association studies (GWAS) analysis was performed using the Gene Atlas database. GSEA, gene set variation analysis (GSVA), protein pathway analysis, and WGCNA were employed to investigate relevant pathways in subtypes. The Harmonizome database was employed to analyze the expression of hub mitochondria-related genes across various human tissues. Single-cell data analysis was conducted to examine patterns of gene expression distribution and pseudo-temporal changes. Additionally, The real-time polymerase chain reaction (RT-PCR) was used to validate the expression of these hub mitochondria-related genes.

Results

In OA, the mitochondria-related pathway was significantly activated. Nine hub mitochondria-related genes (SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4) were identified. They constructed predictive models with good ability to predict OA. These genes are primarily associated with macrophages. Unsupervised consensus clustering identified two mitochondria-associated isoforms that are primarily associated with metabolism. Single-cell analysis showed that they were all expressed in single cells and varied with cell differentiation. RT-PCR showed that they were all significantly expressed in OA.

Conclusion

SIRT4, DNAJC15, NFS1, FKBP8, SLC25A37, CARS2, MTHFD2, ETFDH, and PDK4 are potential mitochondrial target genes for studying OA. The classification of mitochondria-associated isoforms could help to personalize treatment for OA patients.

Programmed Cell Death-Related Gene Signature Associated with Prognosis and Immune Infiltration and the Roles of HMOX1 in the Proliferation and Apoptosis were Investigated in Uveal Melanoma

BACKGROUND

Uveal melanoma (UVM) is the most common primary ocular malignancy, with a wide range of symptoms and outcomes. The programmed cell death (PCD) plays an important role in tumor development, diagnosis, and prognosis. There is still no research on the relationship between PCD-related genes and UVM. A novel PCD-associated prognostic model is urgently needed to improve treatment strategies.

OBJECTIVE

We aim to screen PCD-related prognostic signature and investigate its proliferation ability and apoptosis in UVM cells.

METHODS

The clinical information and RNA-seq data of the UVM patients were collected from the TCGA cohort. All the patients were classified using consensus clustering by the selected PCD-related genes. After univariate Cox regression and PPI network analysis, the prognostic PCD-related genes were then submitted to the LASSO regression analysis to build a prognostic model. The level of immune infiltration of 8-PCD signature in high- and low-risk patients was analyzed using xCell. The prediction on chemotherapy and immunotherapy response in UVM patients was assessed by GDSC and TIDE algorithm. CCK-8, western blot and Annexin V-FITC/PI staining were used to explore the roles of HMOX1 in UVM cells.

RESULTS

A total of 8-PCD signature was constructed and the risk score of the PCD signature was negatively correlated with the overall survival, indicating strong predictive ability and independent prognostic value. The risk score was positively correlated with CD8 Tcm, CD8 Tem and Th2 cells. Immune cells in high-risk group had poorer overall survival. The drug sensitivity demonstrated that cisplatin might impact the progression of UVM and better immunotherapy responsiveness in the high-risk group. Finally, Overespression HMOX1 (OE-HMOX1) decreased the cell viability and induced apoptosis in UVM cells. Recuse experiment results showed that ferrostatin-1 (fer-1) protected MP65 cells from apoptosis and necrosis caused by OE-HMOX1.

CONCLUSION

The PCD signature may have a significant role in the tumor microenvironment, clinicopathological characteristics, prognosis and drug sensitivity. More importantly, HMOX1 depletion greatly induced tumor cell growth and inhibited cell apoptosis and fer-1 protected UVM cells from apoptosis and necrosis induced by OE-HMOX1. This work provides a foundation for effective therapeutic strategy in tumour treatment.

Integrated Analysis Construct a Tumor-Associated Macrophage Novel Signature with Promising Implications in Predicting the Prognosis and Immunotherapeutic Response of Gastric Cancer Patients

BACKGROUND

Gastric cancer (GC) remains one of the most prevalent malignant tumors worldwide. At present, tumor-associated macrophages (TAMs) are essential in the progression, metastasis, and drug resistance of tumors. Therefore, TAMs can be a crucial target for tumor treatment.

AIMS

We intended to investigate the TAM characteristics in GC and develop a risk signature based on TAM to predict the prognosis of GC patients.

METHODS

The single-cell RNA sequencing (scRNA-seq) and bulk RNA-seq data were acquired from a publicly available database. We utilized the Seurat pipeline to process the scRNA-seq data and determine TAM cell types using marker genes. Univariate Cox regression analysis was utilized to examine TAM-related prognostic genes, and then we employed Lasso-Cox regression analysis, and Multivariate Cox regression analysis established a novel risk profile to forecast the clinical value of the model with a new nomogram combining risk profiles and clinicopathological characteristics.

RESULTS

The current study employed scRNA-seq data to identify five TAM clusters in GC, among which four were significantly associated with GC prognosis. Accordingly, we further developed a TAM-related risk signature utilizing nine genes. After evaluation, our model accurately predicted the prognosis of gastric cancer. Generally, GC patients with low TAMS scores exhibited a more favorable prognosis, greater benefits from immunotherapy, and higher levels of immune cell infiltration.

CONCLUSIONS

The prognosis of GC can be effectively predicted by TAM-based risk signatures, and the signature may provide a new perspective for comprehensively guiding clinical diagnosis, prediction, and immunotherapy for gastric cancer.

Machine learning to establish three sphingolipid metabolism genes signature to characterize the immune landscape and prognosis of patients with gastric cancer

BACKGROUND

Gastric cancer (GC) is one of the most common malignant tumors worldwide. Nevertheless, GC still lacks effective diagnosed and monitoring method and treating targets. This study used multi omics data to explore novel biomarkers and immune therapy targets around sphingolipids metabolism genes (SMGs).

METHOD

LASSO regression analysis was performed to filter prognostic and differently expression SMGs among TCGA and GTEx data. Risk score model and Kaplan-Meier were built to validate the prognostic SMG signature and prognostic nomogram was further constructed. The biological functions of SMG signature were annotated via multi omics. The heterogeneity landscape of immune microenvironment in GC was explored. qRT-PCR was performed to validate the expression level of SMG signature. Competing endogenous RNA regulatory network was established to explore the molecular regulatory mechanisms.

RESULT

3-SMGs prognostic signature (GLA, LAMC1, TRAF2) and related nomogram were constructed combing several clinical characterizes. The expression difference and diagnostic value were validated by PCR data. Multi omics data reveals 3-SMG signature affects cell cycle and death via several signaling pathways to regulate GC progression. Overexpression of 3-SMG signature influenced various immune cell infiltration in GC microenvironment. RBP-SMGs-miRNA-mRNAs/lncRNAs regulatory network was built to annotate regulatory system.

CONCLUSION

Upregulated 3-SMGs signature are excellent predictive diagnosed and prognostic biomarkers, providing a new perspective for future GC immunotherapy.

Multiomics-based classifier to decipher immune landscape of uveal melanoma and predict patient outcomes

Uveal melanoma (UVM) prognosis and the possibilities for targeted therapy depend on a thorough understanding of immune infiltration features and the analysis of genomic and immune signatures. Leveraging multi-omics data from The Cancer Genome Atlas and GEO datasets, we employed an unsupervised clustering algorithm to categorize UVM into immune-related subgroups. Subsequent multi-omics analysis revealed two distinct UVM subtypes, each characterized by unique genomic mutations and immune microenvironment disparities. The aggressive UMCS2 subtype exhibited higher TNM stage and poorer survival, marked by elevated metabolism and increased immune infiltration. However, UMCS2 displayed heightened tumor mutational burden and immune dysfunction, leading to reduced responsiveness to immunotherapy. Importantly, these subtypes demonstrated differential sensitivity to targeted drugs due to significant variances in metabolic and immune environments, with UMCS2 displaying lower sensitivity. We developed a robust, subtype-specific marker-based risk scoring system. This system's diagnostic accuracy was validated through ROC curves, decision curve analysis, and calibration curves, all yielding satisfactory results. Additionally, cell experiments identified the pivotal function of HTR2B, the most crucial factor in this risk model. Knocking down HTR2B significantly reduced the activity, proliferation, and invasion ability of the UVM cell line. These findings underscored the impact of gene and immune microenvironment alterations in driving distinct molecular subtypes, emphasizing the need for precise treatment strategies. The molecular subtyping-based risk assessment system not only aids in predicting patient prognosis but also guides the identification of populations suitable for combined treatment. Molecules represented by HTR2B in the model may serve as effective therapeutic targets for UVM.Communicated by Ramaswamy H. Sarma.

Diagnostic value of sphingolipid metabolism-related genes CD37 and CXCL9 in nonalcoholic fatty liver disease

The development of nonalcoholic fatty liver disease (NAFLD) has been reported to be caused by sphingolipid family inducing insulin resistance, mitochondrial dysfunction, and inflammation, which can be regulated by multiple sphingolipid metabolic pathways. This study aimed to explore the molecular mechanism of crucial sphingolipid metabolism related genes (SMRGs) in NAFLD. Firstly, the datasets (GSE48452, GSE126848, and GSE63067) from the Gene Expression Omnibus database and sphingolipid metabolism genes (SMGs) from previous research were collected for this study. The differentially expressed genes (DEGs) between different NAFLD and controls were acquired through "limma," and the SMRGs were authenticated via weighted gene co-expression network analysis (WGCNA). After overlapping the DEGs and SMRGs, the causality between the intersection genes (DE-SMRGs) and NAFLD was explored to sort out the candidate biomarkers by Mendelian randomization (MR) study. The receiver operating characteristic (ROC) curves of candidate biomarkers in GSE48452 and GSE126848 were yielded to determine the biomarkers, followed by the nomogram construction and enrichment analysis. Finally, the immune infiltration analysis, the prediction of transcription factors (TFs) and drugs targeting biomarkers were put into effect. A total of 23 DE-SMRGs were acquired based on the differential analysis and weighted gene co-expression network analysis (WGCNA), of which 3 DE-SMRGs (CD37, CXCL9 and IL7R) were picked out for follow-up analysis through univariate and multivariate MR analysis. The values of area under ROC curve of CD37 and CXCL9 were >0.7 in GSE48452 and GSE126848, thereby being regarded as biomarkers, which were mainly enriched in amino acid metabolism. With respect to the Spearman analysis between immune cells and biomarkers, CD37 and CXCL9 were significantly positively associated with M1 macrophages (P < .001), whose proportion was observably higher in NAFLD patients compared with controls. At last, TFs (ZNF460 and ZNF384) of CD37 and CXCL9 and a total of 79 chemical drugs targeting CD37 and CXCL9 were predicted. This study mined the pivotal SMRGs, CD37 and CXCL9, and systematically explored the mechanism of action of both biomarkers based on the public databases, which could tender a fresh reference for the clinical diagnosis and therapy of NAFLD.

A novel sphingolipid metabolism-related long noncoding RNA signature predicts the prognosis, immune landscape and therapeutic response in pancreatic adenocarcinoma

Sphingolipid metabolism affects prognosis and resistance to immunotherapy in patients with cancer and is an emerging target in cancer therapy with promising diagnostic and prognostic value. Long noncoding ribonucleic acids (lncRNAs) broadly regulate tumour-associated metabolic reprogramming. However, the potential of sphingolipid metabolism-related lncRNAs in pancreatic adenocarcinoma (PAAD) is poorly understood. In this study, coexpression algorithms were employed to identify sphingolipid metabolism-related lncRNAs. The least absolute shrinkage and selection operator (LASSO) algorithm was used to develop a sphingolipid metabolism-related lncRNA signature (SMLs). The prognostic predictive stability of the SMLs was validated using Kaplan-Meier. Univariate and multivariate Cox, receiver operating characteristic (ROC) and clinical stratification analyses were used to comprehensively assess the SMLs. Gene set variation analysis (GSVE), gene ontology (GO) and tumor mutation burden (TMB) analysis explored the potential mechanisms. Additionally, single sample gene set enrichment analysis (ssGSEA), ESTIMATE, immune checkpoints and drug sensitivity analysis were used to investigate the potential predictive function of the SMLs. Finally, an SMLs-based consensus clustering algorithm was utilized to differentiate patients and determine the suitable population for immunotherapy. The results showed that the SMLs consists of seven sphingolipid metabolism-related lncRNAs, which can well determine the clinical outcome of individuals with PAAD, with high stability and general applicability. In addition, the SMLs-based consensus clustering algorithm divided the TCGA-PAAD cohort into two clusters, with Cluster 1 showing better survival than Cluster 2. Additionally, Cluster 1 had a higher level of immune cell infiltration than Cluster 2, which combined with the higher levels of immune checkpoints in Cluster 1 suggests that Cluster 1 is more consistent with an immune 'hot tumor' profile and may respond better to immune checkpoint inhibitors (ICIs). This study offers new insights regarding the potential role of sphingolipid metabolism-related lncRNAs as biomarkers in PAAD. The constructed SMLs and the SMLs-based clustering are valuable tools for predicting clinical outcomes in PAAD and provide a basis for clinical selection of individualized treatments.

Systematic analysis of the prognostic value and immunological function of LTBR in human cancer

Lymphotoxin beta receptor (LTBR) is a positive T cell proliferation regulator gene. It is closely associated with the tumor immune microenvironment. However, its role in cancer and immunotherapy is unclear. Firstly, the expression level and prognostic value of LTBR were analyzed. Secondly, the expression of LTBR in clinical stages, immune subtypes, and molecular subtypes was analyzed. The correlation between LTBR and immune regulatory genes, immune checkpoint genes, and RNA modification genes was then analyzed. Correlations between LTBR and immune cells, scores, cancer-related functional status, tumor stemness index, mismatch repair (MMR) genes, and DNA methyltransferase were also analyzed. In addition, we analyzed the role of LTBR in DNA methylation, mutational status, tumor mutation burden (TMB), and microsatellite instability (MSI). Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Gene Set Enrichment Analysis (GSEA) were used to explore the role of LTBR in pan-cancer. Finally, the drugs associated with LTBR were analyzed. The expression of LTBR was confirmed using quantitative real-time PCR and Western blot. LTBR is significantly overexpressed in most cancers and is associated with low patient survival. In addition, LTBR expression was strongly correlated with immune cells, score, cancer-related functional status, tumor stemness index, MMR genes, DNA methyltransferase, DNA methylation, mutational status, TMB, and MSI. Enrichment analysis revealed that LTBR was associated with apoptosis, necroptosis, and immune-related pathways. Finally, multiple drugs targeting LTBR were identified. LTBR is overexpressed in several tumors and is associated with a poor prognosis. It is related to immune-related genes and immune cell infiltration.

Construction of a diagnostic model for hepatitis B-related hepatocellular carcinoma using machine learning and artificial neural networks and revealing the correlation by immunoassay

HBV infection profoundly escalates hepatocellular carcinoma (HCC) susceptibility, responsible for a majority of HCC cases. HBV-driven immune-mediated hepatocyte impairment significantly fuels HCC progression. Regrettably, inconspicuous early HCC symptoms often culminate in belated diagnoses. Nevertheless, surgically treated early-stage HCC patients relish augmented five-year survival rates. In contrast, advanced HCC exhibits feeble responses to conventional interventions like radiotherapy, chemotherapy, and surgery, leading to diminished survival rates. This investigation endeavors to unearth diagnostic hallmark genes for HBV-HCC leveraging a bioinformatics framework, thus refining early HBV-HCC detection. Candidate genes were sieved via differential analysis and Weighted Gene Co-Expression Network Analysis (WGCNA). Employing three distinct machine learning algorithms unearthed three feature genes (HHIP, CXCL14, and CDHR2). Melding these genes yielded an innovative Artificial Neural Network (ANN) diagnostic blueprint, portending to alleviate patient encumbrance and elevate life quality. Immunoassay scrutiny unveiled accentuated immune damage in HBV-HCC patients relative to solitary HCC. Through consensus clustering, HBV-HCC was stratified into two subtypes (C1 and C2), the latter potentially indicating milder immune impairment. The diagnostic model grounded in these feature genes showcased robust and transferrable prognostic potentialities, introducing a novel outlook for early HBV-HCC diagnosis. This exhaustive immunological odyssey stands poised to expedite immunotherapeutic curatives' emergence for HBV-HCC.

Construction of a prognostic model for lung adenocarcinoma based on heat shock protein-related genes and immune analysis

Lung adenocarcinoma (LUAD) represents a prevalent form of cancer, with low early diagnosis rates and high mortality rates, posing a global health challenge. Heat shock proteins (HSPs) assume a crucial role within the tumor immune microenvironment (TME) of LUAD. Here, a collection of 97 HSP-related genes (HSPGs) was assembled based on prior literature reports, of which 36 HSPGs were differentially expressed in LUAD. In The Cancer Genome Atlas (TCGA) cohort, we constructed a prognostic model for risk stratification and prognosis prediction by integrating 13 HSPGs. In addition, the prognostic significance and predictive efficacy of the HSP-related riskscore were examined and validated in the Gene Expression Omnibus (GEO) cohort. To facilitate the clinical use of this riskscore, we also established a nomogram scale by verifying its effectiveness through different methods. In light of these outcomes, we concluded a significant correlation between HSPs and TME in LUAD, and the riskscore can be a reliable prognostic indicator. Furthermore, this study evaluated the differences in immunophenoscore, tumor immune dysfunction and exclusion score, and sensitivity to several common chemotherapy drugs among LUAD individuals in different risk groups, which may aid in clinical decision-making for immune therapy and chemotherapy in LUAD individuals.

Machine learning-enhanced insights into sphingolipid-based prognostication: revealing the immunological landscape and predictive proficiency for immunomotherapy and chemotherapy responses in pancreatic carcinoma

Background: With a poor prognosis for affected individuals, pancreatic adenocarcinoma (PAAD) is known as a complicated and diverse illness. Immunocytes have become essential elements in the development of PAAD. Notably, sphingolipid metabolism has a dual function in the development of tumors and the invasion of the immune system. Despite these implications, research on the predictive ability of sphingolipid variables for PAAD prognosis is strikingly lacking, and it is yet unclear how they can affect PAAD immunotherapy and targeted pharmacotherapy. Methods: The investigation process included SPG detection while also being pertinent to the prognosis for PAAD. Both the analytical capability of CIBERSORT and the prognostic capability of the pRRophetic R package were used to evaluate the immunological environments of the various HCC subtypes. In addition, CCK-8 experiments on PAAD cell lines were carried out to confirm the accuracy of drug sensitivity estimates. The results of these trials, which also evaluated cell survival and migratory patterns, confirmed the usefulness of sphingolipid-associated genes (SPGs). Results: As a result of this thorough investigation, 32 SPGs were identified, each of which had a measurable influence on the dynamics of overall survival. This collection of genes served as the conceptual framework for the development of a prognostic model, which was carefully assembled from 10 chosen genes. It should be noted that this grouping of patients into cohorts with high and low risk was a sign of different immune profiles and therapy responses. The increased abundance of SPGs was identified as a possible sign of inadequate responses to immune-based treatment approaches. The careful CCK-8 testing carried out on PAAD cell lines was of the highest importance for providing clear confirmation of drug sensitivity estimates. Conclusion: The significance of Sphingolipid metabolism in the complex web of PAAD development is brought home by this study. The novel risk model, built on the complexity of sphingolipid-associated genes, advances our understanding of PAAD and offers doctors a powerful tool for developing personalised treatment plans that are specifically suited to the unique characteristics of each patient.

Constructing a prognostic model for head and neck squamous cell carcinoma based on glucose metabolism related genes

Background

Glucose metabolism (GM) plays a crucial role in cancer cell proliferation, tumor growth, and survival. However, the identification of glucose metabolism-related genes (GMRGs) for effective prediction of prognosis in head and neck squamous cell carcinoma (HNSC) is still lacking.

Methods

We conducted differential analysis between HNSC and Normal groups to identify differentially expressed genes (DEGs). Key module genes were obtained using weighted gene co-expression network analysis (WGCNA). Intersection analysis of DEGs, GMRGs, and key module genes identified GMRG-DEGs. Univariate and multivariate Cox regression analyses were performed to screen prognostic-associated genes. Independent prognostic analysis of clinical traits and risk scores was implemented using Cox regression. Gene set enrichment analysis (GSEA) was used to explore functional pathways and genes between high- and low-risk groups. Immune infiltration analysis compared immune cells between the two groups in HNSC samples. Drug prediction was performed using the Genomics of Drug Sensitivity in Cancer (GDSC) database. Quantitative real-time fluorescence PCR (qRT-PCR) validated the expression levels of prognosis-related genes in HNSC patients.

Results

We identified 4973 DEGs between HNSC and Normal samples. Key gene modules, represented by black and brown module genes, were identified. Intersection analysis revealed 76 GMRG-DEGs. Five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) were identified. A nomogram incorporating age, lymph node status (N), and risk score was constructed for survival prediction in HNSC patients. Immune infiltration analysis showed significant differences in five immune cell types (Macrophages M0, memory B cells, Monocytes, Macrophages M2, and Dendritic resting cells) between the high- and low-risk groups. GDSC database analysis identified 53 drugs with remarkable differences between the groups, including A.443654 and AG.014699. DNMT1 and MTHFD2 were up-regulated, while MPZ was down-regulated in HNSC.

Conclusion

Our study highlights the significant association of five prognosis-related genes (MTHFD2, CDKN2A, TPM2, MPZ, and DNMT1) with HNSC. These findings provide further evidence of the crucial role of GMRGs in HNSC.

Decoding tumor heterogeneity in uveal melanoma: basement membrane genes as novel biomarkers and therapeutic targets revealed by multi-omics approaches for cancer immunotherapy

Background: Uveal melanoma (UVM) is a primary intraocular malignancy that poses a significant threat to patients' visual function and life. The basement membrane (BM) is critical for establishing and maintaining cell polarity, adult function, embryonic and organ morphogenesis, and many other biological processes. Some basement membrane protein genes have been proven to be prognostic biomarkers for various cancers. This research aimed to develop a novel risk assessment system based on BMRGs that would serve as a theoretical foundation for tailored and accurate treatment. Methods: We used gene expression profiles and clinical data from the TCGA-UVM cohort of 80 UVM patients as a training set. 56 UVM patients from the combined cohort of GSE84976 and GSE22138 were employed as an external validation dataset. Prognostic characteristics of basement membrane protein-related genes (BMRGs) were characterized by Lasso, stepwise multifactorial Cox. Multivariate analysis revealed BMRGs to be independent predictors of UVM. The TISCH database probes the crosstalk of BMEGs in the tumor microenvironment at the single-cell level. Finally, we investigated the function of ITGA5 in UVM using multiple experimental techniques, including CCK8, transwell, wound healing assay, and colony formation assay. Results: There are three genes in the prognostic risk model (ADAMTS10, ADAMTS14, and ITGA5). After validation, we determined that the model is quite reliable and accurately forecasts the prognosis of UVM patients. Immunotherapy is more likely to be beneficial for UVM patients in the high-risk group, whereas the survival advantage may be greater for UVM patients in the low-risk group. Knockdown of ITGA5 expression was shown to inhibit the proliferation, migration, and invasive ability of UVM cells in vitro experiments. Conclusion: The 3-BMRGs feature model we constructed has excellent predictive performance which plays a key role in the prognosis, informing the individualized treatment of UVM patients. It also provides a new perspective for assessing pre-immune efficacy.

Multi-omics analyses based on genes associated with oxidative stress and phospholipid metabolism revealed the intrinsic molecular characteristics of pancreatic cancer

Oxidative stress (OS), which impacts lipid metabolic reprogramming, can affect the biological activities of cancer cells. How oxidative stress and phospholipid metabolism (OSPM) influence the prognosis of pancreatic cancer (PC) needs to be elucidated. The metabolic data of 35 pancreatic tumor samples, 34 para-carcinoma samples, and 31 normal pancreatic tissues were obtained from the previously published literature. Pan-cancer samples were obtained from The Cancer Genome Atlas (TCGA). And the Gene Expression Omnibus (GEO), International Cancer Genome Consortium (ICGC), ArrayExpress, and the Genotype-Tissue Expression (GTEx) databases were searched for more PC and normal pancreatic samples. The metabolites in PC were compared with normal and para-carcinoma tissues. The characteristics of the key OSPM genes were summarized in pan-cancer. The random survival forest analysis and multivariate Cox regression analysis were utilized to construct an OSPM-related signature. Based on this signature, PC samples were divided into high- and low-risk subgroups. The dysregulations of the tumor immune microenvironment were further investigated. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was conducted to investigate the expression of genes in the signature in PC and normal tissues. The protein levels of these genes were further demonstrated. In PC, metabolomic studies revealed the alteration of PM, while transcriptomic studies showed different expressions of OSPM-related genes. Then 930 PC samples were divided into three subtypes with different prognoses, and an OSPM-related signature including eight OSPM-related genes (i.e., SLC2A1, MMP14, TOP2A, MBOAT2, ANLN, ECT2, SLC22A3, and FGD6) was developed. High- and low-risk subgroups divided by the signature showed different prognoses, expression levels of immune checkpoint genes, immune cell infiltration, and tumor microenvironment. The risk score was negatively correlated with the proportion of TIL, pDC, Mast cell, and T cell co-stimulation. The expression levels of genes in the signature were verified in PC and normal samples. The protein levels of SLC2A1, MMP14, TOP2A, MBOAT2, ANLN, and SLC22A3 showed up-regulation in PC samples compared with normal tissues. After integrating metabolomics and transcriptomics data, the alterations in OSPM in PC were investigated, and an OSPM-related signature was developed, which was helpful for the prognostic assessment and individualized treatment for PC.

GADD45B regulates the carcinogenesis process of chronic atrophic gastritis and the metabolic pathways of gastric cancer

Background

Gastric cancer continues to be a significant global healthcare challenge, and its burden remains substantial. The development of gastric cancer (GC) is closely linked to chronic atrophic gastritis (CAG), yet there is a scarcity of research exploring the underlying mechanisms of CAG-induced carcinogenesis.

Methods

In this study, we conducted a comprehensive investigation into the oncogenes involved in CAG using both bulk transcriptome and single-cell transcriptome data. Our approach employed hdWGCNA to identify pathogenic genes specific to CAG, with non-atrophic gastritis (NAG) serving as the control group. Additionally, we compared CAG with GC, using normal gastric tissue as the control group in the single-cell transcriptome analysis. By intersecting the identified pathogenic genes, we pinpointed key network molecules through protein interaction network analysis. To further refine the gene selection, we applied LASSO, SVM-RFE, and RF techniques, which resulted in a set of cancer-related genes (CRGs) associated with CAG. To identify CRGs potentially linked to gastric cancer progression, we performed a univariate COX regression analysis on the gene set. Subsequently, we explored the relationship between CRGs and immune infiltration, drug sensitivity, and clinical characteristics in gastric cancer patients. We employed GSVA to investigate how CRGs regulated signaling pathways in gastric cancer cells, while an analysis of cell communication shed light on the impact of CRGs on signal transmission within the gastric cancer tumor microenvironment. Lastly, we analyzed changes in metabolic pathways throughout the progression of gastric cancer.

Results

Using hdWGCNA, we have identified a total of 143 pathogenic genes that were shared by CAG and GC. To further investigate the underlying mechanisms, we conducted protein interaction network analysis and employed machine learning screening techniques. As a result, we have identified 15 oncogenes that are specifically associated with chronic atrophic gastritis. By performing ROC reanalysis and prognostic analysis, we have determined that GADD45B is the most significant gene involved in the carcinogenesis of CAG. Immunohistochemical staining and differential analysis have revealed that GADD45B expression was low in GC tissues while high in normal gastric tissues. Moreover, based on prognostic analysis, high expression of GADD45B has been correlated with poor prognosis in GC patients. Additionally, an analysis of immune infiltration has shown a relationship between GADD45B and the infiltration of various immune cells. By correlating GADD45B with clinical characteristics, we have found that it primarily affects the depth of invasion in GC. Through cell communication analysis, we have discovered that the CD99 signaling pathway network and the CDH signaling pathway network are the main communication pathways that significantly alter the microenvironment of gastric tissue during the development of chronic atrophic gastritis. Specifically, GADD45B-low GC cells were predominantly involved in the network communication of the CDH signaling pathway, while GADD45B-high GC cells played a crucial role in both signaling pathways. Furthermore, we have identified several metabolic pathways, including D-Glutamine and D-glutamate metabolism and N-Glycan biosynthesis, among others, that played important roles in the occurrence and progression of GC, in addition to the six other metabolic pathways. In summary, our study highlighted the discovery of 143 pathogenic genes shared by CAG and GC, with a specific focus on 15 oncogenes associated with CAG. We have identified GADD45B as the most important gene in the carcinogenesis of CAG, which exhibited differential expression in GC tissues compared to normal gastric tissues. Moreover, GADD45B expression was correlated with patient prognosis and is associated with immune cell infiltration. Our findings also emphasized the impact of the CD99 and CDH signaling pathway networks on the microenvironment of gastric tissue during the development of CAG. Additionally, we have identified key metabolic pathways involved in GC progression.

Conclusion

GADD45B, an oncogene implicated in chronic atrophic gastritis, played a critical role in GC development. Decreased expression of GADD45B was associated with the onset of GC. Moreover, GADD45B expression levels were closely tied to poor prognosis in GC patients, influencing the infiltration patterns of various cells within the tumor microenvironment, as well as impacting the metabolic pathways involved in GC progression.

Proposing new early detection indicators for pancreatic cancer: Combining machine learning and neural networks for serum miRNA-based diagnostic model

Background

Pancreatic cancer (PC) is a lethal malignancy that ranks seventh in terms of global cancer-related mortality. Despite advancements in treatment, the five-year survival rate remains low, emphasizing the urgent need for reliable early detection methods. MicroRNAs (miRNAs), a group of non-coding RNAs involved in critical gene regulatory mechanisms, have garnered significant attention as potential diagnostic and prognostic biomarkers for pancreatic cancer (PC). Their suitability stems from their accessibility and stability in blood, making them particularly appealing for clinical applications.

Methods

In this study, we analyzed serum miRNA expression profiles from three independent PC datasets obtained from the Gene Expression Omnibus (GEO) database. To identify serum miRNAs associated with PC incidence, we employed three machine learning algorithms: Support Vector Machine-Recursive Feature Elimination (SVM-RFE), Least Absolute Shrinkage and Selection Operator (LASSO), and Random Forest. We developed an artificial neural network model to assess the accuracy of the identified PC-related serum miRNAs (PCRSMs) and create a nomogram. These findings were further validated through qPCR experiments. Additionally, patient samples with PC were classified using the consensus clustering method.

Results

Our analysis revealed three PCRSMs, namely hsa-miR-4648, hsa-miR-125b-1-3p, and hsa-miR-3201, using the three machine learning algorithms. The artificial neural network model demonstrated high accuracy in distinguishing between normal and pancreatic cancer samples, with verification and training groups exhibiting AUC values of 0.935 and 0.926, respectively. We also utilized the consensus clustering method to classify PC samples into two optimal subtypes. Furthermore, our investigation into the expression of PCRSMs unveiled a significant negative correlation between the expression of hsa-miR-125b-1-3p and age.

Conclusion

Our study introduces a novel artificial neural network model for early diagnosis of pancreatic cancer, carrying significant clinical implications. Furthermore, our findings provide valuable insights into the pathogenesis of pancreatic cancer and offer potential avenues for drug screening, personalized treatment, and immunotherapy against this lethal disease.

Development of a novel immune infiltration-related diagnostic model for Alzheimer's disease using bioinformatic strategies

Background

The pathogenesis of Alzheimer's disease (AD) is complex and multi-factorial. Increasing evidence has shown the important role of immune infiltration in AD. Thus the current study was designed to identify immune infiltration-related genes and to explore their diagnostic value in AD.

Methods

The expression data of AD patients were downloaded from the GEO database. The limma R package identified differentially expressed genes (DEGs) between AD and controls. The CIBERSORT algorithm identified differentially infiltrated immune cells (DIICs) between AD and controls. DIIC-correlated DEGs were obtained by Pearson correlation analysis. WGCNA was employed to identify DIIC-related modules. Next, LASSO, RFE, and RF machine learning methods were applied to screen robust DIIC-related gene signatures in AD, followed by the construction and validation of a diagnostic nomogram. Detection of the expression of related genes in the peripheral blood of Alzheimer's disease and healthy volunteers by RT-PCR. In addition, the CTD database predicted chemicals targeting DIIC-related gene signatures in the treatment of AD.

Results

NK cells, M0 macrophages, activated myeloid dendritic cells, resting mast cells, CD8+ T cells, resting memory CD4+ T cells, gamma delta T cells, and M2 macrophages were differentially infiltrated between AD and controls. Pearson analysis identified a total of 277 DIIC-correlated DEGs between AD and controls. Thereafter, 177 DIIC-related genes were further obtained by WGCNA analysis. By LASSO, RFE and RF algorithms, CMTM2, DDIT4, LDHB, NDUFA1, NDUFB2, NDUFS5, RPL17, RPL21, RPL26 and NDUFAF2 were identified as robust gene signature in AD. The results of RT-PCR detection of peripheral blood samples from Alzheimer's disease and healthy volunteers showed that the expression trend of ten genes screened was consistent with the detection results; among them, the expression levels of CMTM2, DDIT4, LDHB, NDUFS5, and RPL21 are significantly different among groups. Thus, a diagnostic nomogram based on a DIIC-related signature was constructed and validated. Moreover, candidate chemicals targeting those biomarkers in the treatment of AD, such as 4-hydroxy-2-nonenal, rosiglitazone, and resveratrol, were identified in the CTD database.

Conclusion

For the first time, we identified 10 immune infiltration-related biomarkers in AD, which may be helpful for the diagnosis of AD and provide guidance in the treatment of AD.

Identification of novel prognostic circRNA biomarkers in circRNA-miRNA-mRNA regulatory network in gastric cancer and immune infiltration analysis

BACKGROUND

Gastric cancer (GC) carries significant morbidity and mortality globally. An increasing number of studies have confirmed that circular RNA (circRNA) is tightly associated with the carcinogenesis and development of GC, especially acting as a competing endogenous RNA for miRNAs.

OBJECTIVE

Our study aimed to construct the circRNA-miRNA-mRNA regulatory network and analyze the function and prognostic significance of the network using bioinformatics tools.

METHODS

We first downloaded the GC expression profile from the Gene Expression Omnibus database and identified differentially expressed genes and differentially expressed circRNAs. Then, we predicted the miRNA-mRNA interaction pairs and constructed the circRNA-miRNA-mRNA regulatory network. Next, we established a protein-protein interaction network and analyzed the function of these networks. Finally, we primarily validated our results by comparison with The Cancer Genome Atlas cohort and by performing qRT-PCR.

RESULTS

We screened the top 15 hub genes and 3 core modules. Functional analysis showed that in the upregulated circRNA network, 15 hub genes were correlated with extracellular matrix organization and interaction. The function of downregulated circRNAs converged on physiological functions, such as protein processing, energy metabolism and gastric acid secretion. We ascertained 3 prognostic and immune infiltration-related genes, COL12A1, COL5A2, and THBS1, and built a nomogram for clinical application. We validated the expression level and diagnostic performance of key prognostic differentially expressed genes.

CONCLUSIONS

In conclusion, we constructed two circRNA-miRNA-mRNA regulatory networks and identified 3 prognostic and screening biomarkers, COL12A1, COL5A2, and THBS1. The ceRNA network and these genes could play important roles in GC development, diagnosis and prognosis.

An N6-methyladenosine regulation- and mRNAsi-related prognostic index reveals the distinct immune microenvironment and immunotherapy responses in lower-grade glioma

BACKGROUND

N6-methyladenosine (m6A) modification is involved in tumorigenesis and progression as well as closely correlated with stem cell differentiation and pluripotency. Moreover, tumor progression includes the acquisition of stemness characteristics and accumulating loss of differentiation phenotype. Therefore, we integrated m6A modification and stemness indicator mRNAsi to classify patients and predict prognosis for LGG.

METHODS

We performed consensus clustering, weighted gene co-expression network analysis, and least absolute shrinkage and selection operator Cox regression analysis to identify an m6A regulation- and mRNAsi-related prognostic index (MRMRPI). Based on this prognostic index, we also explored the differences in immune microenvironments between high- and low-risk populations. Next, immunotherapy responses were also predicted. Moreover, single-cell RNA sequencing data was further used to verify the expression of these genes in MRMRPI. At last, the tumor-promoting and tumor-associated macrophage polarization roles of TIMP1 in LGG were validated by in vitro experiments.

RESULTS

Ten genes (DGCR10, CYP2E1, CSMD3, HOXB3, CABP4, AVIL, PTCRA, TIMP1, CLEC18A, and SAMD9) were identified to construct the MRMRPI, which was able to successfully classify patients into high- and low-risk group. Significant differences in prognosis, immune microenvironment, and immunotherapy responses were found between distinct groups. A nomogram integrating the MRMRPI and other prognostic factors were also developed to accurately predict prognosis. Moreover, in vitro experiments illustrated that inhibition of TIMP1 could inhibit the proliferation, migration, and invasion of LGG cells and also inhibit the polarization of tumor-associated macrophages.

CONCLUSION

These findings provide novel insights into understanding the interactions of m6A methylation regulation and tumor stemness on LGG development and contribute to guiding more precise immunotherapy strategies.

A fibroblast-associated signature predicts prognosis and immunotherapy in esophageal squamous cell cancer

Background

Current paradigms of anti-tumor therapies are not qualified to evacuate the malignancy ascribing to cancer stroma's functions in accelerating tumor relapse and therapeutic resistance. Cancer-associated fibroblasts (CAFs) has been identified significantly correlated with tumor progression and therapy resistance. Thus, we aimed to probe into the CAFs characteristics in esophageal squamous cancer (ESCC) and construct a risk signature based on CAFs to predict the prognosis of ESCC patients.

Methods

The GEO database provided the single-cell RNA sequencing (scRNA-seq) data. The GEO and TCGA databases were used to obtain bulk RNA-seq data and microarray data of ESCC, respectively. CAF clusters were identified from the scRNA-seq data using the Seurat R package. CAF-related prognostic genes were subsequently identified using univariate Cox regression analysis. A risk signature based on CAF-related prognostic genes was constructed using Lasso regression. Then, a nomogram model based on clinicopathological characteristics and the risk signature was developed. Consensus clustering was conducted to explore the heterogeneity of ESCC. Finally, PCR was utilized to validate the functions that hub genes play on ESCC.

Results

Six CAF clusters were identified in ESCC based on scRNA-seq data, three of which had prognostic associations. A total of 642 genes were found to be significantly correlated with CAF clusters from a pool of 17080 DEGs, and 9 genes were selected to generate a risk signature, which were mainly involved in 10 pathways such as NRF1, MYC, and TGF-Beta. The risk signature was significantly correlated with stromal and immune scores, as well as some immune cells. Multivariate analysis demonstrated that the risk signature was an independent prognostic factor for ESCC, and its potential in predicting immunotherapeutic outcomes was confirmed. A novel nomogram integrating the CAF-based risk signature and clinical stage was developed, which exhibited favorable predictability and reliability for ESCC prognosis prediction. The consensus clustering analysis further confirmed the heterogeneity of ESCC.

Conclusion

The prognosis of ESCC can be effectively predicted by CAF-based risk signatures, and a comprehensive characterization of the CAF signature of ESCC may aid in interpreting the response of ESCC to immunotherapy and offer new strategies for cancer treatment.

Mast cell marker gene signature: prognosis and immunotherapy response prediction in lung adenocarcinoma through integrated scRNA-seq and bulk RNA-seq

Background

Mast cells, comprising a crucial component of the tumor immune milieu, modulate neoplastic progression by secreting an array of pro- and antitumorigenic factors. Numerous extant studies have produced conflicting conclusions regarding the impact of mast cells on the prognosis of patients afflicted with lung adenocarcinoma (LUAD).

Methods

Employing single-cell RNA sequencing (scRNA-seq) analysis, mast cell-specific marker genes in LUAD were ascertained. Subsequently, a mast cell-related genes (MRGs) signature was devised to stratify LUAD patients into high- and low-risk cohorts based on the median risk value. Further investigations were conducted to assess the influence of distinct risk categories on the tumor microenvironment. The prognostic import and capacity to prognosticate immunotherapy benefits of the MRGs signature were corroborated using four external cohorts. Ultimately, the functional roles of SYAP1 were validated through in vitro experimentation.

Results

After scRNA-seq and bulk RNA-seq data analysis, we established a prognostic signature consisting of nine MRGs. This profile effectively distinguished favorable survival outcomes in both the training and validation cohorts. In addition, we identified the low-risk group as a population more effective for immunotherapy. In cellular experiments, we found that silencing SYAP1 significantly reduced the proliferation, invasion and migratory capacity of LUAD cells while increasing apoptosis.

Conclusion

Our MRGs signature offers valuable insights into the involvement of mast cells in determining the prognosis of LUAD and may prove instrumental as a navigational aid for immunotherapy selection, as well as a predictor of immunotherapy response in LUAD patients.

Machine learning-based signature of necrosis-associated lncRNAs for prognostic and immunotherapy response prediction in cutaneous melanoma and tumor immune landscape characterization

Background

Cutaneous melanoma (CM) is one of the malignant tumors with a relative high lethality. Necroptosis is a novel programmed cell death that participates in anti-tumor immunity and tumor prognosis. Necroptosis has been found to play an important role in tumors like CM. However, the necroptosis-associated lncRNAs' potential prognostic value in CM has not been identified.

Methods

The RNA sequencing data collected from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression Project (GTEx) was utilized to identify differentially expressed genes in CM. By using the univariate Cox regression analysis and machine learning LASSO algorithm, a prognostic risk model had been built depending on 5 necroptosis-associated lncRNAs and was verified by internal validation. The performance of this prognostic model was assessed by the receiver operating characteristic curves. A nomogram was constructed and verified by calibration. Furthermore, we also performed sub-group K-M analysis to explore the 5 lncRNAs' expression in different clinical stages. Function enrichment had been analyzed by GSEA and ssGSEA. In addition, qRT-PCR was performed to verify the five lncRNAs' expression level in CM cell line (A2058 and A375) and normal keratinocyte cell line (HaCaT).

Results

We constructed a prognostic model based on five necroptosis-associated lncRNAs (AC245041.1, LINC00665, AC018553.1, LINC01871, and AC107464.3) and divided patients into high-risk group and low-risk group depending on risk scores. A predictive nomogram had been built to be a prognostic indicator to clinical factors. Functional enrichment analysis showed that immune functions had more relationship and immune checkpoints were more activated in low-risk group than that in high-risk group. Thus, the low-risk group would have a more sensitive response to immunotherapy.

Conclusion

This risk score signature could be used to divide CM patients into low- and high-risk groups, and facilitate treatment strategy decision making that immunotherapy is more suitable for those in low-risk group, providing a new sight for CM prognostic evaluation.

Identification and validation of TME-related signatures to predict prognosis and response to anti-tumor therapies in skin cutaneous melanoma

The tumor microenvironment (TME) dynamically regulates cancer progression and affects clinical outcomes. This study aimed to identify molecular subtypes and construct a prognostic risk model based on TME-related signatures in skin cutaneous melanoma (SKCM) patients. We categorized SKCM patients based on transcriptome data of SKCM from The Cancer Genome Atlas (TCGA) database and 29 TME-related gene signatures. Differentially expressed genes were identified using univariate Cox regression and Lasso regression analysis, which were used for risk model construction. The robustness of this model was validated in independent external cohorts. Genetic landscape alterations, immune characteristics, and responsiveness to immunotherapy/chemotherapy were evaluated. Three TME-related subtypes were identified, and subtype C3 exhibited the most favorable prognosis, had enriched immune-related pathways, and possessed more infiltration of T_cells_CD8, T_cells_CD4_memory_activated, and Macrophages_M1 but a lower TumorPurity, whereas Macrophages_M2 were increased in subtype C1 and subtype C2. Subtype C1 was more sensitive to Cisplatin, subtype C2 was more sensitive to Temozolomide, and subtype C3 was more sensitive to Paclitaxel; 8 TME-related genes (NOTCH3, HEYL, ZNF703, ABCC2, PAEP, CCL8, HAPLN3, and HPDL) were screened for risk model construction. High-risk patients had dismal prognosis with good prediction performance. Moreover, low-risk patients were more sensitive to Paclitaxel and Temozolomide, whereas high-risk patients were more sensitive to Cisplatin. This risk model had robustness in predicting prognosis in SKCM patients. The results facilitate the understanding of TME-related genes in SKCM and provide a TME-related genes-based predictive model in prognosis and direction of personalized options for SKCM patients.

Biomarkers predicting clinical outcomes in nasopharyngeal cancer patients receiving immune checkpoint inhibitors: A systematic review and meta-analysis

Background

Optimal biomarkers to select patients who will benefit most from immunotherapy remain lacking in nasopharyngeal cancer (NPC). This systematic review and meta-analysis aimed to evaluate the association between various biomarkers and clinical outcomes in NPC patients treated with immune checkpoint inhibitors (ICIs).

Methods

Systematic searches of PubMed, Embase, Cochrane Library, and Web of Science databases were performed up to October 2022. Studies evaluating the association between biomarkers and intended outcomes of ICIs were included. The pooled odds ratio (OR) and hazard ratio (HR) with 95% confidence intervals (CIs) were calculated, respectively, for the objective response rate (ORR) and progression-free survival (PFS) under fixed or random-effect models.

Results

A total of 15 studies involving 1,407 patients were included. The pooled analysis indicated that NPC patients with lower plasma Epstein-Barr virus (EBV) DNA level at baseline (OR = 2.14, 95% CI: 1.46-3.14, P &lt; 0.001), decreased EBV DNA load during immunotherapy (OR = 4.57, 95% CI: 2.24-9.34, P = 0.002) and higher programmed cell death-ligand 1 (PD-L1) expression (OR = 2.35, 95% CI: 1.36-4.09, P = 0.002) had superior ORR than the counterparts. No significant differences of ORR were observed between positive PD-L1 expression and negative PD-L1 expression (OR = 1.50, 95% CI: 0.92-2.45, P = 0.104), as well as higher tumor mutation burden (TMB) and lower TMB (OR = 1.62, 95% CI: 0.41-6.44, P = 0.494). Patients with lower plasma EBV DNA level at baseline obtained a significant benefit on PFS than those with higher plasma EBV DNA level (HR = 0.52, 95% CI: 0.42-0.63, P &lt; 0.001). There were no differences in PFS between decreased EBV DNA load and increased EBV DNA load during immunotherapy (HR = 0.51, 95% CI: 0.22-1.17, P = 0.109), higher PD-L1 expression and lower PD-L1 expression (HR = 0.65, 95% CI: 0.42-1.01, P = 0.054), positive PD-L1 expression and negative PD-L1 expression (HR = 0.90, 95% CI: 0.64-1.26, P = 0.531), lower TMB and higher TMB (HR = 0.84, 95% CI: 0.51-1.38, P = 0.684).

Conclusion

Lower baseline plasma EBV DNA level, decreased plasma EBV DNA during immunotherapy, and higher PD-L1 expression are reliable biomarkers predicting better response to ICIs treatment. Lower baseline plasma EBV DNA level was also associated with longer PFS. It is warranted to further explore and better illuminate the utility of these biomarkers in future clinical trials and real-world practice.

Systematic review registration

https://www.crd.york.ac.uk/PROSPERO/, identifier CRD42022324434.

Identification of novel prognostic biomarkers in the TF-enhancer-target regulatory network in hepatocellular carcinoma and immune infiltration analysis

Background: Hepatocellular carcinoma (HCC) remains notorious for its high malignancy, poor prognosis and high mortality. The exploration of novel therapeutic agents for HCC has remained challenging due to its complex aetiology. Therefore, it is necessary to elucidate the pathogenesis and mechanism of HCC for clinical intervention.

Methods: We collected data from several public data portals and systematically analysed the association between transcription factors (TFs), eRNA-associated enhancers and downstream targets. We next filtered the prognostic genes and established a novel prognosis-related nomogram model. Moreover, we explored the potential mechanisms of the identified prognostic genes. The expression level was validated by several ways.

Results: We first constructed a significant TF-enhancer-target regulatory network and identified DAPK1 as a coregulatory differentially expressed prognosis-related gene. We combined common clinicopathological factors and built a prognostic nomogram model for HCC. We found that our regulatory network was correlated with the processes of synthesizing various substances. Moreover, we explored the role of DAPK1 in HCC and found that it was associated with immune cell infiltration and DNA methylation. Several immunostimulators and targeting drugs could be promising immune therapy targets. The tumor immune microenvironment was analyzed. Finally, the lower DAPK1 expression in HCC was validated via the GEO database, UALCAN cohort, and qRT-PCR.

Conclusion: In conclusion, we established a significant TF-enhancer-target regulatory network and identified downregulated DAPK1 as an important prognostic and diagnostic gene in HCC. Its potential biological functions and mechanisms were annotated using bioinformatics tools.

Development of a risk model to predict prognosis in breast cancer based on cGAS-STING-related genes

Background: Breast cancer (BRCA) is regarded as a lethal and aggressive cancer with increasing morbidity and mortality worldwide. cGAS-STING signaling regulates the crosstalk between tumor cells and immune cells in the tumor microenvironment (TME), emerging as an important DNA-damage mechanism. However, cGAS-STING-related genes (CSRGs) have rarely been investigated for their prognostic value in breast cancer patients.

Methods: Our study aimed to construct a risk model to predict the survival and prognosis of breast cancer patients. We obtained 1087 breast cancer samples and 179 normal breast tissue samples from the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEX) database, 35 immune-related differentially expression genes (DEGs) from cGAS-STING-related genes were systematically assessed. The Cox regression was applied for further selection, and 11 prognostic-related DEGs were used to develop a machine learning-based risk assessment and prognostic model.

Results: We successfully developed a risk model to predict the prognostic value of breast cancer patients and its performance acquired effective validation. The results derived from Kaplan-Meier analysis revealed that the low-risk score patients had better overall survival (OS). The nomogram that integrated the risk score and clinical information was established and had good validity in predicting the overall survival of breast cancer patients. Significant correlations were observed between the risk score and tumor-infiltrating immune cells, immune checkpoints and the response to immunotherapy. The cGAS-STING-related genes risk score was also relevant to a series of clinic prognostic indicators such as tumor staging, molecular subtype, tumor recurrence, and drug therapeutic sensibility in breast cancer patients.

Conclusion: cGAS-STING-related genes risk model provides a new credible risk stratification method to improve the clinical prognostic assessment for breast cancer.

Prognostic signatures of sphingolipids: Understanding the immune landscape and predictive role in immunotherapy response and outcomes of hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a complex disease with a poor outlook for patients in advanced stages. Immune cells play an important role in the progression of HCC. The metabolism of sphingolipids functions in both tumor growth and immune infiltration. However, little research has focused on using sphingolipid factors to predict HCC prognosis. This study aimed to identify the key sphingolipids genes (SPGs) in HCC and develop a reliable prognostic model based on these genes.

Methods

The TCGA, GEO, and ICGC datasets were grouped using SPGs obtained from the InnateDB portal. A prognostic gene signature was created by applying LASSO-Cox analysis and evaluating it with Cox regression. The validity of the signature was verified using ICGC and GEO datasets. The tumor microenvironment (TME) was examined using ESTIMATE and CIBERSORT, and potential therapeutic targets were identified through machine learning. Single-cell sequencing was used to examine the distribution of signature genes in cells within the TME. Cell viability and migration were tested to confirm the role of the key SPGs.

Results

We identified 28 SPGs that have an impact on survival. Using clinicopathological features and 6 genes, we developed a nomogram for HCC. The high- and low-risk groups were found to have distinct immune characteristics and response to drugs. Unlike CD8 T cells, M0 and M2 macrophages were found to be highly infiltrated in the TME of the high-risk subgroup. High levels of SPGs were found to be a good indicator of response to immunotherapy. In cell function experiments, SMPD2 and CSTA were found to enhance survival and migration of Huh7 cells, while silencing these genes increased the sensitivity of Huh7 cells to lapatinib.

Conclusion

The study presents a six-gene signature and a nomogram that can aid clinicians in choosing personalized treatments for HCC patients. Furthermore, it uncovers the connection between sphingolipid-related genes and the immune microenvironment, offering a novel approach for immunotherapy. By focusing on crucial sphingolipid genes like SMPD2 and CSTA, the efficacy of anti-tumor therapy can be increased in HCC cells.

Cuproptosis signature and PLCD3 predicts immune infiltration and drug responses in osteosarcoma

Osteosarcoma (OS) is a cancer that is frequently found in children and adolescents and has made little improvement in terms of prognosis in recent years. A recently discovered type of programmed cell death called cuproptosis is mediated by copper ions and the tricarboxylic acid (TCA) cycle. The expression patterns, roles, and prognostic and predictive capabilities of the cuproptosis regulating genes were investigated in this work. TARGET and GEO provided transcriptional profiling of OS. To find different cuproptosis gene expression patterns, consensus clustering was used. To identify hub genes linked to cuproptosis, differential expression (DE) and weighted gene co-expression network analysis (WGCNA) were used. Cox regression and Random Survival Forest were used to build an evaluation model for prognosis. For various clusters/subgroups, GSVA, mRNAsi, and other immune infiltration experiments were carried out. The drug-responsive study was carried out by the Oncopredict algorithm. Cuproptosis genes displayed two unique patterns of expression, and high expression of FDX1 was associated with a poor outcome in OS patients. The TCA cycle and other tumor-promoting pathways were validated by the functional study, and activation of the cuproptosis genes may also be connected with immunosuppressive state. The robust survival prediction ability of a five-gene prognostic model was verified. This rating method also took stemness and immunosuppressive characteristics into account. Additionally, it can be associated with a higher sensitivity to medications that block PI3K/AKT/mTOR signaling as well as numerous chemoresistances. U2OS cell migration and proliferation may be encouraged by PLCD3. The relevance of PLCD3 in immunotherapy prediction was verified. The prognostic significance, expressing patterns, and functions of cuproptosis in OS were revealed in this work on a preliminary basis. The cuproptosis-related scoring model worked well for predicting prognosis and chemoresistance.

Integrated pan-cancer analysis and experimental verification of the roles of tropomyosin 4 in gastric cancer

Objective

To investigate the function of tropomyosin 4 (TPM4) using pan-cancer data, especially in gastric cancer (GC), using comprehensive bioinformatics analysis and molecular experiments.

Methods

We used UCSC Xena, The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression Project (GTEx), TIMER2.0, GEPIA, cBioPortal, Xiantao tool, and UALCAN websites and databases for the extraction of pan-cancer data on TPM4. TPM4 expression was investigated with respect to prognosis, genetic alterations, epigenetic alterations, and immune infiltration. RNA22, miRWalk, miRDB, Starbase 2.0, and Cytoscape were used for identifying and constructing the regulatory networks of lncRNAs, miRNAs, and TPM4 in GC. Data from GSCALite, drug bank databases, and Connectivity Map (CMap) were used to analyze the sensitivity of drugs dependent on TPM4 expression. Gene Ontology (GO), enrichment analyses of the Kyoto Encyclopedia of Genes and Genomes (KEGG), wound healing assays, and (Matrigel) transwell experiments were used to investigate the biological functions of TPM4 in GC.

Result

The findings of the comprehensive pan-cancer analysis revealed that TPM4 has a certain diagnostic and prognosis value in most cancers. Alterations in the expression of TPM4, including duplications and deep mutations, and epigenetic alterations revealed that TPM4 expression is related to the expression of DNA methylation inhibitors and RNA methylation regulators at high concentrations. Besides, TPM4 expression was found to correlate with immune cell infiltration, immune checkpoint (ICP) gene expression, the tumor mutational burden (TMB), and microsatellite instability (MSI). Neoantigens (NEO) were also found to influence its response to immunotherapy. A lncRNA-miRNA -TPM4 network was found to regulate GC development and progression. TPM4 expression was related to docetaxel,5-fluorouracil, and eight small molecular targeted drugs sensitivity. Gene function enrichment analyses revealed that genes that were co-expressed with TPM4 were enriched within the extracellular matrix (ECM)-related pathways. Wound-healing and (Matrigel) transwell assays revealed that TPM4 promotes cell migration and invasion. TPM4, as an oncogene, plays a biological role, perhaps via ECM remodeling in GC.

Conclusions

TPM4 is a prospective marker for the diagnosis, treatment outcome, immunology, chemotherapy, and small molecular drugs targeted for pan-cancer treatment, including GC treatment. The lncRNA-miRNA-TPM4network regulates the mechanism underlying GC progression. TPM4 may facilitate the invasion and migration of GC cells, possibly through ECM remodeling.

Circadian rhythm-related genes index: A predictor for HNSCC prognosis, immunotherapy efficacy, and chemosensitivity

Background

Head and neck squamous cell carcinoma (HNSCC) is the most common head and neck cancer and is highly aggressive and heterogeneous, leading to variable prognosis and immunotherapy outcomes. Circadian rhythm alterations in tumourigenesis are of equal importance to genetic factors and several biologic clock genes are considered to be prognostic biomarkers for various cancers. The aim of this study was to establish reliable markers based on biologic clock genes, thus providing a new perspective for assessing immunotherapy response and prognosis in patients with HNSCC.

Methods

We used 502 HNSCC samples and 44 normal samples from the TCGA-HNSCC dataset as the training set. 97 samples from GSE41613 were used as an external validation set. Prognostic characteristics of circadian rhythm-related genes (CRRGs) were established by Lasso, random forest and stepwise multifactorial Cox. Multivariate analysis revealed that CRRGs characteristics were independent predictors of HNSCC, with patients in the high-risk group having a worse prognosis than those in the low-risk group. The relevance of CRRGs to the immune microenvironment and immunotherapy was assessed by an integrated algorithm.

Results

6-CRRGs were considered to be strongly associated with HNSCC prognosis and a good predictor of HNSCC. The riskscore established by the 6-CRRG was found to be an independent prognostic factor for HNSCC in multifactorial analysis, with patients in the low-risk group having a higher overall survival (OS) than the high-risk group. Nomogram prediction maps constructed from clinical characteristics and riskscore had good prognostic power. Patients in the low-risk group had higher levels of immune infiltration and immune checkpoint expression and were more likely to benefit from immunotherapy.

Conclusion

6-CRRGs play a key predictive role for the prognosis of HNSCC patients and can guide physicians in selecting potential responders to prioritise immunotherapy, which could facilitate further research in precision immuno-oncology.

Identification of copper metabolism-related subtypes and establishment of the prognostic model in ovarian cancer

BACKGROUND

Ovarian cancer (OC) is one of the most common and most malignant gynecological malignancies in gynecology. On the other hand, dysregulation of copper metabolism (CM) is closely associated with tumourigenesis and progression. Here, we investigated the impact of genes associated with copper metabolism (CMRGs) on the prognosis of OC, discovered various CM clusters, and built a risk model to evaluate patient prognosis, immunological features, and therapy response.

METHODS

15 CMRGs affecting the prognosis of OC patients were identified in The Cancer Genome Atlas (TCGA). Consensus Clustering was used to identify two CM clusters. lasso-cox methods were used to establish the copper metabolism-related gene prognostic signature (CMRGPS) based on differentially expressed genes in the two clusters. The GSE63885 cohort was used as an external validation cohort. Expression of CM risk score-associated genes was verified by single-cell sequencing and quantitative real-time PCR (qRT-PCR). Nomograms were used to visually depict the clinical value of CMRGPS. Differences in clinical traits, immune cell infiltration, and tumor mutational load (TMB) between risk groups were also extensively examined. Tumour Immune Dysfunction and Rejection (TIDE) and Immune Phenotype Score (IPS) were used to validate whether CMRGPS could predict response to immunotherapy in OC patients.

RESULTS

In the TCGA and GSE63885 cohorts, we identified two CM clusters that differed significantly in terms of overall survival (OS) and tumor microenvironment. We then created a CMRGPS containing 11 genes to predict overall survival and confirmed its reliable predictive power for OC patients. The expression of CM risk score-related genes was validated by qRT-PCR. Patients with OC were divided into low-risk (LR) and high-risk (HR) groups based on the median CM risk score, with better survival in the LR group. The 5-year AUC value reached 0.74. Enrichment analysis showed that the LR group was associated with tumor immune-related pathways. The results of TIDE and IPS showed a better response to immunotherapy in the LR group.

CONCLUSION

Our study, therefore, provides a valuable tool to further guide clinical management and tailor the treatment of patients with OC, offering new insights into individualized treatment.

Identification of anoikis-related genes classification patterns and immune infiltration characterization in ischemic stroke based on machine learning

Introduction

Ischemic stroke (IS) is a type of stroke that leads to high mortality and disability. Anoikis is a form of programmed cell death. When cells detach from the correct extracellular matrix, anoikis disrupts integrin junctions, thus preventing abnormal proliferating cells from growing or attaching to an inappropriate matrix. Although there is growing evidence that anoikis regulates the immune response, which makes a great contribution to the development of IS, the role of anoikis in the pathogenesis of IS is rarely explored.

Methods

First, we downloaded GSE58294 set and GSE16561 set from the NCBI GEO database. And 35 anoikis-related genes (ARGs) were obtained from GSEA website. The CIBERSORT algorithm was used to estimate the relative proportions of 22 infiltrating immune cell types. Next, consensus clustering method was used to classify ischemic stroke samples. In addition, we used least absolute shrinkage and selection operator (LASSO), support vector machine-recursive feature elimination (SVM-RFE) and random forest (RF) algorithms to screen the key ARGs in ischemic stroke. Next, we performed receiver operating characteristics (ROC) analysis to assess the accuracy of each diagnostic gene. At the same time, the nomogram was constructed to diagnose IS by integrating trait genes. Then, we analyzed the correlation between gene expression and immune cell infiltration of the diagnostic genes in the combined database. And gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) analysis were performed on these genes to explore differential signaling pathways and potential functions, as well as the construction and visualization of regulatory networks using NetworkAnalyst and Cytoscape. Finally, we investigated the expression pattern of ARGs in IS patients across age or gender.

Results

Our study comprehensively analyzed the role of ARGs in IS for the first time. We revealed the expression profile of ARGs in IS and the correlation with infiltrating immune cells. And The results of consensus clustering analysis suggested that we can classify IS patients into two clusters. The machine learning analysis screened five signature genes, including AKT1, BRMS1, PTRH2, TFDP1 and TLE1. We also constructed nomogram models based on the five risk genes and evaluated the immune infiltration correlation, gene-miRNA, gene-TF and drug-gene interaction regulatory networks of these signature genes. The expression of ARGs did not differ by sex or age.

Discussion

This study may provide a beneficial reference for further elucidating the pathogenesis of IS, and render new ideas for drug screening, individualized therapy and immunotherapy of IS.

Crosstalk of disulfidptosis-related subtypes, establishment of a prognostic signature and immune infiltration characteristics in bladder cancer based on a machine learning survival framework

Background

Bladder cancer (BLCA) is the most common malignancy of the urinary tract. On the other hand, disulfidptosis, a mechanism of disulfide stress-induced cell death, is closely associated with tumorigenesis and progression. Here, we investigated the impact of disulfidptosis-related genes (DRGs) on the prognosis of BLCA, identified various DRG clusters, and developed a risk model to assess patient prognosis, immunological profile, and treatment response.

Methods

The expression and mutational characteristics of four DRGs were first analyzed in bulk RNA-Seq and single-cell RNA sequencing data, IHC staining identified the role of DRGs in BLCA progression, and two DRG clusters were identified by consensus clustering. Using the differentially expressed genes (DEGs) from these two clusters, we transformed ten machine learning algorithms into more than 80 combinations and finally selected the best algorithm to construct a disulfidptosis-related prognostic signature (DRPS). We based this selection on the mean C-index of three BLCA cohorts. Furthermore, we explored the differences in clinical characteristics, mutational landscape, immune cell infiltration, and predicted efficacy of immunotherapy between high and low-risk groups. To visually depict the clinical value of DRPS, we employed nomograms. Additionally, we verified whether DRPS predicts response to immunotherapy in BLCA patients by utilizing the Tumour Immune Dysfunction and Rejection (TIDE) and IMvigor 210 cohorts.

Results

In the integrated cohort, we identified several DRG clusters and DRG gene clusters that differed significantly in overall survival (OS) and tumor microenvironment. After the integration of clinicopathological features, DRPS showed robust predictive power. Based on the median risk score associated with disulfidptosis, BLCA patients were divided into low-risk (LR) and high-risk (HR) groups, with patients in the LR group having a better prognosis, a higher tumor mutational load and being more sensitive to immunotherapy and chemotherapy.

Conclusion

Our study, therefore, provides a valuable tool to further guide clinical management and tailor the treatment of BLCA patients, offering new insights into individualized treatment.