Identification and validation of neurotrophic factor-related gene signatures in glioblastoma and Parkinson's disease

Machine Learning Integration with Single-Cell Transcriptome Sequencing Datasets Reveals the Impact of Tumor-Associated Neutrophils on the Immune Microenvironment and Immunotherapy Outcomes in Gastric Cancer

The characteristics of neutrophils play a crucial role in defining the tumor inflammatory environment. However, the function of tumor-associated neutrophils (TANs) in tumor immunity and their response to immune checkpoint inhibitors (ICIs) remains incompletely understood. By analyzing single-cell RNA sequencing data from over 600,000 cells in gastric cancer (GSE163558 and GSE183904), colorectal cancer (GSE205506), and lung cancer (GSE207422), we identified neutrophil subsets in primary gastric cancer that are associated with the treatment response to ICIs. Specifically, we focused on neutrophils with high expression of CD44 (CD44_NEU), which are abundant during tumor progression and exert significant influence on the gastric cancer immune microenvironment. Machine learning analysis revealed 22 core genes associated with CD44_NEU, impacting inflammation, proliferation, migration, and oxidative stress. In addition, multiple immunofluorescence staining and gastric cancer spatial transcriptome data (GSE203612) showed a correlation between CD44_NEU and T-cell infiltration in gastric cancer tissues. A risk score model derived from seven essential genes (AQP9, BASP1, BCL2A1, PLEK, PDE4B, PROK2, and ACSL1) showed better predictive capability for patient survival compared to clinical features alone, and integrating these scores with clinical variables resulted in a prognostic nomogram. Overall, this study highlights the heterogeneity of TANs, particularly the CD44_NEU critical influence on immunotherapy outcomes, paving the way for personalized treatment strategies.

Casein Kinase 2 Inhibitor, CX-4945, Induces Apoptosis and Restores Blood-Brain Barrier Homeostasis in In Vitro and In Vivo Models of Glioblastoma

Background: In oncology, casein kinase 2 (CK2), a serine/threonine kinase, has a dual action, regulating cellular processes and acting as an oncogenic promoter. Methods: This study examined the effect of CX-4945, a selective CK2 inhibitor, in a human U-87 glioblastoma (GBM) cell line, treated with CX-4945 (5, 10, and 15 μM) for 24 h. Similarly, the hCMEC/D3 cell line was used to mimic the blood-brain barrier (BBB), examining the ability of CX-4945 to restore BBB homeostasis, after stimulation with lipopolysaccharide (LPS) and then treated with CX-4945 (5, 10, and 15 μM). Results: We reported that CX-4945 reduced the proliferative activity and modulated the main pathways involved in tumor progression including apoptosis. Furthermore, in confirmation of the in vitro study, performing a xenograft model, we demonstrated that CX-4945 exerted promising antiproliferative effects, also restoring the tight junctions' expression. Conclusions: These new insights into the molecular signaling of CK2 in GBM and BBB demonstrate that CX-4945 could be a promising approach for future GBM therapy, not only in the tumor microenvironment but also at the BBB level.

Don't ignore the middle: Distinct mnemonic functions of intermediate hippocampus

The dorsal region of the hippocampus (dHC) mediates many of the mnemonic functions traditionally associated with the hippocampus proper, such as spatial and episodic memory, whereas ventral hippocampus (vHC) has been extensively implicated in emotional memory and motivational processes. By contrast, the functions of the intermediate hippocampus (iHC) are far less understood. In this study, we aimed to investigate the mnemonic functions of iHC by reversibly inactivating iHC prior to testing memory in behavioral tasks dependent on the integrity of dHC, iHC, or vHC, namely, rapid place water maze, inhibitory avoidance, spontaneous alternation, and temporal ordering of odors. Given our previous findings showing that dHC and vHC are involved in mnemonic control of ingestive behavior, we also assessed the effects of iHC inactivation on sucrose intake. The results showed that pharmacological inhibition of iHC impairs rapid place water maze memory, which has been previously shown to be dependent on iHC but not dHC or vHC. iHC inactivation does not impact memory dependent on dHC (spontaneous alternation), vHC (temporal odor memory), or either dHC or vHC (inhibitory avoidance), and only modestly affects sucrose intake. These findings provide support for the involvement of iHC in mnemonic functions that are distinct from dHC and vHC and highlight the need to further advance our understanding of the functions of this hippocampal region that has been relatively understudied.

Discovery and Validation of Ferroptosis-Associated Genes of Ulcerative Colitis

Background

Ulcerative colitis (UC) is a long-lasting idiopathic condition, but its precise mechanisms remain unclear. Meanwhile, evidence has demonstrated that ferroptosis seems to interlock with the progress of UC. This research sought to identify hub genes of UC related to ferroptosis.

Methods

First, the relevant profiles for this article were obtained from GEO database. From the FerrDb, 479 genes linked to ferroptosis were retrieved. Using analysis of the difference and WGCNA on colonic samples from GSE73661, the remaining six hub genes linked to ferroptosis and UC were discovered. Through logistic regression analyses, the diagnostic model was constructed and was then evaluated by external validation using dataset GSE92415. Afterwards, the correlation between immune cell filtration in UC and hub genes was examined. Finally, a mice model of colitis was established, and the results were verified using qRT-PCR.

Results

We acquired six hub genes linked to ferroptosis and UC. In order to create a diagnostic model for UC, we used logistic regression analysis to screen three of the six ferroptosis related genes (HIF1A, SLC7A11, and LPIN1). The ROC curve showed that the three hub genes had outstanding potential for disease diagnosis (AUC = 0.976), which was subsequently validated in samples from GSE92415 (AUC = 0.962) and blood samples from GSE3365 (AUC = 0.847) and GSE94648 (AUC = 0.769). These genes might be crucial for UC immunity based upon the results on the immune system. Furthermore, mouse samples examined using qRT-PCR also verified our findings.

Conclusion

In conclusion, the findings have important implications for ferroptosis and UC, and these hub genes may also offer fresh perspectives on the aetiology and therapeutic approaches of UC.

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.

Identification and validation of a novel Parkinson-Glioma feature gene signature in glioma and Parkinson's disease

Introduction

The prognosis for glioma is generally poor, and the 5-year survival rate for patients with this disease has not shown significant improvement over the past few decades. Parkinson's disease (PD) is a prevalent movement disorder, ranking as the second most common neurodegenerative disease after Alzheimer's disease. Although Parkinson's disease and glioma are distinct diseases, they may share certain underlying biological pathways that contribute to their development.

Objective

This study aims to investigate the involvement of genes associated with Parkinson's disease in the development and prognosis of glioma.

Methods

We obtained datasets from the TCGA, CGGA, and GEO databases, which included RNA sequencing data and clinical information of glioma and Parkinson's patients. Eight machine learning algorithms were used to identify Parkinson-Glioma feature genes (PGFGs). PGFGs associated with glioma prognosis were identified through univariate Cox analysis. A risk signature was constructed based on PGFGs using Cox regression analysis and the Least Absolute Shrinkage and Selection Operator (LASSO) method. We subsequently validated its predictive ability using various methods, including ROC curves, calibration curves, KM survival analysis, C-index, DCA, independent prognostic analysis, and stratified analysis. To validate the reproducibility of the results, similar work was performed on three external test datasets. Additionally, a meta-analysis was employed to observe the heterogeneity and consistency of the signature across different datasets. We also compared the differences in genomic variations, functional enrichment, immune infiltration, and drug sensitivity analysis based on risk scores. This exploration aimed to uncover potential mechanisms of glioma occurrence and prognosis.

Results

We identified 30 PGFGs, of which 25 were found to be significantly associated with glioma survival. The prognostic signature, consisting of 19 genes, demonstrated excellent predictive performance for 1-, 2-, and 3-year overall survival (OS) of glioma. The signature emerged as an independent prognostic factor for glioma overall survival (OS), surpassing the predictive performance of traditional clinical variables. Notably, we observed differences in the tumor microenvironment (TME), levels of immune cell infiltration, immune gene expression, and drug resistance analysis among distinct risk groups. These findings may have significant implications for the clinical treatment of glioma patients.

Conclusion

The expression of genes related to Parkinson's disease is closely associated with the immune status and prognosis of glioma patients, potentially regulating glioma pathogenesis through multiple mechanisms. The interaction between genes associated with Parkinson's disease and the immune system during glioma development provides novel insights into the molecular mechanisms and targeted therapies for glioma.

Integrating single-cell and multi-omic approaches reveals Euphorbiae Humifusae Herba-dependent mitochondrial dysfunction in non-small-cell lung cancer

Euphorbiae Humifusae Herba (EHH) is a pivotal therapeutic agent with diverse pharmacological effects. However, a substantial gap exists in understanding its pharmacological properties and anti-tumour mechanisms. This study aimed to address this gap by exploring EHH's pharmacological properties, identifying NSCLC therapy-associated protein targets, and elucidating how EHH induces mitochondrial disruption in NSCLC cells, offering insights into novel NSCLC treatment strategies. String database was utilized to explore protein-protein interactions. Subsequently, single-cell analysis and multi-omics further unveiled the impact of EHH-targeted genes on the immune microenvironment of NSCLC, as well as their influence on immunotherapeutic responses. Finally, both in vivo and in vitro experiments elucidated the anti-tumour mechanisms of EHH, specifically through the assessment of mitochondrial ROS levels and alterations in mitochondrial membrane potential. EHH exerts its influence through engagement with a cluster of 10 genes, including the apoptotic gene CASP3. This regulatory impact on the immune milieu within NSCLC holds promise as an indicator for predicting responses to immunotherapy. Besides, EHH demonstrated the capability to induce mitochondrial ROS generation and perturbations in mitochondrial membrane potential in NSCLC cells, ultimately leading to mitochondrial dysfunction and consequent apoptosis of tumour cells. EHH induces mitochondrial disruption in NSCLC cells, leading to cell apoptosis to inhibit the progress of NSCLC.

KHDRBS1 as a novel prognostic signaling biomarker influencing hepatocellular carcinoma cell proliferation, migration, immune microenvironment, and drug sensitivity

Background

Human tumors pose significant challenges, with targeted therapy against specific molecular targets or signaling pathways being a mainstay alongside surgical resection. Previous studies have implicated KHDRBS1 in the oncogenesis of certain human tumors such as colorectal and prostate cancers, underscoring its potential as a therapeutic target. However, the comprehensive expression pattern of KHDRBS1 in hepatocellular carcinoma (HCC) warrants further exploration.

Methods

Integrating and analyzing multi-omics, multi-cohort data from public databases, coupled with clinical samples and molecular biology validation, we elucidate the oncogenic role of KHDRBS1 in HCC progression. Additionally, leveraging HCC single-cell sequencing data, we segregate malignant cells into KHDRBS1-positive and negative subsets, uncovering significant differences in their expression profiles and functional roles.

Results

Our study identifies KHDRBS1 as a tumor-promoting factor in HCC, with its positivity correlating with tumor progression. Furthermore, we highlight the clinical significance of KHDRBS1-positive malignant cells, aiming to further propel its clinical utility.

Conclusion

KHDRBS1 plays a key role in HCC development. This study provides crucial insights for further investigation into KHDRBS1 as a therapeutic target in HCC.

Experimentally validated oxidative stress -associated prognostic signatures describe the immune landscape and predict the drug response and prognosis of SKCM

Background

Skin Cutaneous Melanoma (SKCM) incidence is continually increasing, with chemotherapy and immunotherapy being among the most common cancer treatment modalities. This study aims to identify novel biomarkers for chemotherapy and immunotherapy response in SKCM and explore their association with oxidative stress.

Methods

Utilizing TCGA-SKCM RNA-seq data, we employed Weighted Gene Co-expression Network Analysis (WGCNA) and Protein-Protein Interaction (PPI) networks to identify six core genes. Gene co-expression analysis and immune-related analysis were conducted, and specific markers associated with oxidative stress were identified using Gene Set Variation Analysis (GSVA). Single-cell analysis revealed the expression patterns of Oxidative Stress-Associated Genes (OSAG) in the tumor microenvironment. TIDE analysis was employed to explore the association between immune therapy response and OSAG, while CIBERSORT was used to analyze the tumor immune microenvironment. The BEST database demonstrated the impact of the Oxidative Stress signaling pathway on chemotherapy drug resistance. Immunohistochemical staining and ROC curve evaluation were performed to assess the protein expression levels of core genes in SKCM and normal samples, with survival analysis utilized to determine their diagnostic value.

Results

We identified six central genes associated with SKCM metastasis, among which the expression of DSC2 and DSC3 involved in the oxidative stress pathway was closely related to immune cell infiltration. DSC2 influenced drug resistance in SKMC patients. Furthermore, downregulation of DSC2 and DSC3 expression enhanced the response of SKCM patients to immunotherapy.

Conclusion

This study identified two Oxidative Stress-Associated genes as novel biomarkers for SKCM. Additionally, targeting the oxidative stress pathway may serve as a new strategy in clinical practice to enhance SKCM chemotherapy and sensitivity.

The comprehensive landscape of prognosis, immunity, and function of the GLI family by pan-cancer and single-cell analysis

The Hedgehog (Hh) signaling pathway has been implicated in the pathogenesis of various cancers. However, the roles of the downstream GLI family (GLI1, GLI2, and GLI3) in tumorigenesis remain elusive. This study aimed to unravel the genetic alterations of GLI1/2/3 in cancer and their association with the immune microenvironment and related signaling pathways. Firstly, we evaluated the expression profiles of GLI1/2/3 in different cancer types, analyzed their prognostic and predictive values, and assessed their correlation with tumor-infiltrating immune cells. Secondly, we explored the relationships between GLI1/2/3 and genetic mutations, epigenetic modifications, and clinically relevant drugs. Finally, we performed enrichment analysis to decipher the underlying mechanisms of GLI1/2/3 in cancer initiation and progression. Our results revealed that the expression levels of GLI1/2/3 were positively correlated in most cancer tissues, suggesting a cooperative role of these factors in tumorigenesis. We also identified tissue-specific expression patterns of GLI1/2/3, which may reflect the distinct functions of these factors in different cell types. Furthermore, GLI1/2/3 expression displayed significant associations with poor prognosis in several cancers, indicating their potential as prognostic biomarkers and therapeutic targets. Importantly, we found that GLI1/2/3 modulated the immune microenvironment by regulating the recruitment, activation, and polarization of cancer-associated fibroblasts, endothelial cells, and macrophages. Additionally, functional enrichment analyses indicated that GLI1/2/3 are involved in the regulation of epithelial-mesenchymal transition (EMT). Together, our findings shed new light on the roles of GLI1/2/3 in tumorigenesis and provide a potential basis for the development of novel therapeutic strategies targeting GLI-mediated signaling pathways in cancer.

The Interplay Between Epilepsy and Parkinson's Disease: Gene Expression Profiling and Functional Analysis

The results of many epidemiological studies suggest a bidirectional causality may exist between epilepsy and Parkinson's disease (PD). However, the underlying molecular landscape linking these two diseases remains largely unknown. This study aimed to explore this possible bidirectional causality by identifying differentially expressed genes (DEGs) in each disease as well as their intersection based on two respective disease-related datasets. We performed enrichment analyses and explored immune cell infiltration based on an intersection of the DEGs. Identifying a protein-protein interaction (PPI) network between epilepsy and PD, and this network was visualised using Cytoscape software to screen key modules and hub genes. Finally, exploring the diagnostic values of the identified hub genes. NetworkAnalyst 3.0 and Cytoscape software were also used to construct and visualise the transcription factor-micro-RNA regulatory and co-regulatory networks, the gene-microRNA interaction network, as well as gene-disease association. Based on the enrichment results, the intersection of the DEGs mainly revealed enrichment in immunity-, phosphorylation-, metabolism-, and inflammation-related pathways. The boxplots revealed similar trends in infiltration of many immune cells in epilepsy and Parkinson's disease, with greater infiltration in patients than in controls. A complex PPI network comprising 186 nodes and 512 edges were constructed. According to node connection degree, top 15 hub genes were considered the kernel targets of epilepsy and PD. The area under curve values of hub gene expression profiles confirmed their excellent diagnostic values. This study is the first to analyse the molecular landscape underlying the epidemiological link between epilepsy and Parkinson's disease. The two diseases are closely linked through immunity-, inflammation-, and metabolism-related pathways. This information was of great help in understanding the pathogenesis, diagnosis, and treatment of the diseases. The present results may provide guidance for further in-depth analysis about molecular mechanisms of epilepsy and PD and novel potential targets.

Machine learning and single-cell transcriptome profiling reveal regulation of fibroblast activation through THBS2/TGFβ1/P-Smad2/3 signalling pathway in hypertrophic scar

Hypertrophic scar (HS) is a chronic inflammatory skin disorder characterized by excessive deposition of extracellular matrix, and the mechanisms underlying their formation remain poorly understood. We analysed scRNA-seq data from samples of normal skin and HS. Using the hdWGCNA method, key gene modules of fibroblasts in HS were identified. Non-negative matrix factorization was employed to perform subtype analysis of HS patients using these gene modules. Multiple machine learning algorithms were applied to screen and validate accurate gene signatures for identifying and predicting HS, and a convolutional neural network (CNN) based on deep learning was established and validated. Quantitative reverse transcription-polymerase chain reaction and western blotting were performed to measure mRNA and protein expression. Immunofluorescence was used for gene localization analysis, and biological features were assessed through CCK8 and wound healing assay. Single-cell sequencing revealed distinct subpopulations of fibroblasts in HS. HdWGCNA identified key gene characteristics of this population, and pseudotime analysis was conducted to investigate gene variation during fibroblast differentiation. By employing various machine learning algorithms, the gene range was narrowed down to three key genes. A CNN was trained using the expression of these key genes and immune cell infiltration, enabling diagnosis and prediction of HS. Functional experiments demonstrated that THBS2 is associated with fibroblast proliferation and migration in HS and affects the formation and development of HS through the TGFβ1/P-Smad2/3 pathway. Our study identifies unique fibroblast subpopulations closely associated with HS and provides biomarkers for the diagnosis and treatment of HS.

Glycogen metabolism-mediated intercellular communication in the tumor microenvironment influences liver cancer prognosis

Glycogen metabolism plays a key role in the development of hepatocellular carcinoma (HCC), but the function of glycogen metabolism genes in the tumor microenvironment (TME) is still to be elucidated. Single-cell RNA-seq data were obtained from ten HCC tumor samples totaling 64,545 cells, and 65 glycogen metabolism genes were analyzed by a nonnegative matrix factorization (NMF). The prognosis and immune response of new glycogen TME cell clusters were predicted by using HCC and immunotherapy cohorts from public databases. HCC single-cell analysis was divided into fibroblasts, NT T cells, macrophages, endothelial cells, and B cells, which were separately divided into new cell clusters by glycogen metabolism gene annotation. Pseudo-temporal trajectory analysis demonstrated the temporal differentiation trajectory of different glycogen subtype cell clusters. Cellular communication analysis revealed extensive interactions between endothelial cells with glycogen metabolizing TME cell-related subtypes and different glycogen subtype cell clusters. SCENIC analysis of transcription factors upstream of TME cell clusters with different glycogen metabolism. In addition, TME cell clusters of glycogen metabolism were found to be enriched in expression in CAF subtypes, CD8 depleted, M1, and M2 types. Bulk-seq analysis showed the prognostic significance of glycogen metabolism-mediated TME cell clusters in HCC, while a significant immune response was found in the immunotherapy cohort in patients treated with immune checkpoint blockade (ICB), especially for CAFs, T cells, and macrophages. In summary, our study reveals for the first time that glycogen metabolism mediates intercellular communication in the hepatocellular carcinoma microenvironment while elucidating the anti-tumor mechanisms and immune prognostic responses of different subtypes of cell clusters.

Integrative analysis of gene expression profiles of substantia nigra identifies potential diagnosis biomarkers in Parkinson's disease

Parkinson's disease (PD) is a progressive neurodegenerative disease whose etiology is attributed to development of Lewy bodies and degeneration of dopaminergic neurons in the substantia nigra (SN). Currently, there are no definitive diagnostic indicators for PD. In this study, we aimed to identify potential diagnostic biomarkers for PD and analyzed the impact of immune cell infiltrations on disease pathogenesis. The PD expression profile data for human SN tissue, GSE7621, GSE20141, GSE20159, GSE20163 and GSE20164 were downloaded from the Gene Expression Omnibus (GEO) database for use in the training model. After normalization and merging, we identified differentially expressed genes (DEGs) using the Robust rank aggregation (RRA) analysis. Simultaneously, DEGs after batch correction were identified. Gene interactions were determined through venn Diagram analysis. Functional analyses and protein-protein interaction (PPI) networks were used to the identify hub genes, which were visualized through Cytoscape. A Lasso Cox regression model was employed to identify the potential diagnostic genes. The GSE20292 dataset was used for validation. The proportion of infiltrating immune cells in the samples were determined via the CIBERSORT method. Sixty-two DEGs were screened in this study. They were found to be enriched in nerve conduction, dopamine (DA) metabolism, and DA biosynthesis Gene Ontology (GO) terms. The PPI network and Lasso Cox regression analysis revealed seven potential diagnostic genes, namely SLC18A2, TAC1, PCDH8, KIAA0319, PDE6H, AXIN1, and AGTR1, were subsequently validated in peripheral blood samples obtained from healthy control (HC) and PD patients, as well as in the GSE20292 dataset. The results revealed the exceptional sensitivity and specificity of these genes in PD diagnosis and monitoring. Moreover, PD patients exhibited a higher number of plasma cells, compared to HC individuals. The SLC18A2, TAC1, PCDH8, KIAA0319, PDE6H, AXIN1, and AGTR1 are potential diagnostic biomarkers for PD. Our findings also reveal the essential roles of immune cell infiltration in both disease onset and trajectory.

Comprehending the cuproptosis and cancer-immunity cycle network: delving into the immune landscape and its predictive role in breast cancer immunotherapy responses and clinical endpoints

Background

The role of cuproptosis, a phenomenon associated with tumor metabolism and immunological identification, remains underexplored, particularly in relation to the cancer-immunity cycle (CIC) network. This study aims to rigorously examine the impact of the cuproptosis-CIC nexus on immune reactions and prognostic outcomes in patients with breast cancer (BC), striving to establish a comprehensive prognostic model.

Methods

In the study, we segregated data obtained from TCGA, GEO, and ICGC using CICs retrieved from the TIP database. We constructed a genetic prognostic framework using the LASSO-Cox model, followed by its validation through Cox proportional hazards regression. This framework's validity was further confirmed with data from ICGC and GEO. Explorations of the tumor microenvironment were carried out through the application of ESTIMATE and CIBERSORT algorithms, as well as machine learning techniques, to identify potential treatment strategies. Single-cell sequencing methods were utilized to delineate the spatial distribution of key genes within the various cell types in the tumor milieu. To explore the critical role of the identified CICs, experiments were conducted focusing on cell survival and migration abilities.

Results

In our research, we identified a set of 4 crucial cuproptosis-CICs that have a profound impact on patient longevity and their response to immunotherapy. By leveraging these identified CICs, we constructed a predictive model that efficiently estimates patient prognoses. Detailed analyses at the single-cell level showed that the significance of CICs. Experimental approaches, including CCK-8, Transwell, and wound healing assays, revealed that the protein HSPA9 restricts the growth and movement of breast cancer cells. Furthermore, our studies using immunofluorescence techniques demonstrated that suppressing HSPA9 leads to a notable increase in ceramide levels.

Conclusion

This research outlines a network of cuproptosis-CICs and constructs a predictive nomogram. Our model holds great promise for healthcare professionals to personalize treatment approaches for individuals with breast cancer. The work provides insights into the complex relationship between the cuproptosis-CIC network and the cancer immune microenvironment, setting the stage for novel approaches to cancer immunotherapy. By focusing on the essential gene HSPA9 within the cancer-immunity cycle, this strategy has the potential to significantly improve the efficacy of treatments against breast cancer.

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.

Clinical Significance of Disulfidptosis-related Genes and Functional Analysis in Gastric Cancer

Background: Worldwide, gastric cancer (GC) remains intractable due to its poor prognosis and high morbidity and mortality. Disulfidptosis is a novel kind of cell death mediated by abnormal accumulation of intracellular disulphides. The correlation between disulfidptosis and GC is still unknown. Therefore, it is necessary to elucidate the pathogenesis and mechanism of disulfidptosis and GC for clinical diagnosis and intervention. Methods: RNA-sequencing data from several public data portals and clinical samples were collected. We compared the expression levels of four key genes of disulfidptosis, including SLC7A11, SLC3A2, RPN1, and NCKAP1, in GC and selected prognostic genes to build a novel GC prognosis-related nomogram model. The biological functions and immune landscape of the identified prognostic genes were explored. Results: Overexpressed NCKAP1 and SLC7A11 were prognostic disulfidptosis-related genes in GC. We combined these genes and several clinicopathological factors to build a prognostic nomogram model for GC. Meanwhile, the ROC curves showed that NCKAP1 and SLC7A11 were promising biomarkers for GC screening. The biological and cellular functions were focused on actin activities, GTPase and immunoreaction. The tumour immune microenvironment and immune therapy targets were identified. Competing endogenous RNA network was built to explore the downstream regulatory mechanisms. Finally, the elevated NCKAP1 and SLC7A11 expression in GC was validated via qRT-PCR in a cell line and tissue line. Conclusion: In conclusion, NCKAP1 and SLC7A11 are promising prognostic and diagnostic biomarkers for GC that correlate with the activities of actin, energy metabolism of GTPase, immune infiltration and immunotherapy.

Interplay of Sphingolipid Metabolism in Predicting Prognosis of GBM Patients: Towards Precision Immunotherapy

Background: In spite of numerous existing bio-surveillance systems for predicting glioma (GBM) prognosis, enhancing the efficacy of immunotherapy remains an ongoing conundrum. The continual scrutiny of the dynamic interplay between the sphingolipid metabolic pathway and tumor immunophenotypes has unveiled potential implications. However, the intricate orchestration of functional and regulatory mechanisms by long non-coding RNAs (lncRNAs) in GBM, particularly in the context of sphingolipid metabolism, remains cryptic. Methods: We harnessed established R packages to intersect gene expression profiles of GBM patients within the The Cancer Genome Atlas (TCGA) database with the compilation of sphingolipid metabolism genes from GeneCards. This enabled us to discern markedly distinct lncRNAs, which were subsequently deployed to construct a robust prognostic model utilizing Lasso-Cox regression analysis. We then scrutinized the immune microenvironment across various risk strata using the ssGSEA and CIBERSORT algorithms. To evaluate mutation patterns and drug resistance profiles within patient subgroups, we devised the "Prophytic" and "Maftools" packages, respectively. Results: Our investigation scrutinized lncRNAs linked to sphingolipid metabolism, utilizing glioma specimens from TCGA. We meticulously curated 1224 sphingolipid-associated genes gleaned from GeneCards and pinpointed 272 differentially expressed mRNAs via transcriptomic analysis. Enrichment analyses underscored their significance in sphingolipid processes. A prognostic model founded on 17 meticulously selected lncRNAs was systematically constructed and validated. This model adeptly stratified GBM patients into high- and low-risk categories, yielding highly precise prognostic insights. We also discerned correlations between immune cell infiltration and genetic mutation discrepancies, along with distinct therapeutic responses through drug sensitivity analysis. Notably, computational findings were corroborated through experimental validation by RT-PCR. Conclusion: In summation, our exhaustive inquiry underscores the multifaceted utility of the sphingolipid metabolic pathway as an autonomous diagnostic and prognostic indicator for glioma patients. Furthermore, we amalgamate a profusion of substantiated evidence concerning immune infiltration and gene mutations, thereby reinforcing the proposition that sphingolipid metabolism may function as a pivotal determinant in the panorama of immunotherapeutic interventions.

Combined signature of G protein-coupled receptors and tumor microenvironment provides a prognostic and therapeutic biomarker for skin cutaneous melanoma

BACKGROUND

G protein-coupled receptors (GPCRs) have been shown to have an important role in tumor development and metastasis, and abnormal expression of GPCRs is significantly associated with poor prognosis of tumor patients. In this study, we analyzed the GPCRs-related gene (GPRGs) and tumor microenvironment (TME) in skin cutaneous melanoma (SKCM) to construct a prognostic model to help SKCM patients obtain accurate clinical treatment strategies.

METHODS

SKCM expression data and clinical information were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Differential expression analysis, LASSO algorithm, and univariate and multivariate cox regression analysis were used to screen prognosis-related genes (GPR19, GPR146, S1PR2, PTH1R, ADGRE5, CXCR3, GPR143, and OR2I1P) and multiple prognosis-good immune cells; the data set was analyzed according to above results and build up a GPR-TME classifier. The model was further subjected to immune infiltration, functional enrichment, tumor mutational load, immunotherapy prediction, and scRNA-seq data analysis. Finally, cellular experiments were conducted to validate the functionality of the key gene GPR19 in the model.

RESULTS

The findings indicate that high expression of GPRGs is associated with a poor prognosis in patients with SKCM, highlighting the significant role of GPRGs and the tumor microenvironment (TME) in SKCM development. Notably, the group characterized by low GPR expression and a high TME exhibited the most favorable prognosis and immunotherapeutic efficacy. Furthermore, cellular assays demonstrated that knockdown of GPR19 significantly reduced the proliferation, migration, and invasive capabilities of melanoma cells in A375 and A2058 cell lines.

CONCLUSION

This study provides novel insights for the prognosis evaluation and treatment of melanoma, along with the identification of a new biomarker, GPR19.

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.

Interaction, immune infiltration characteristics and prognostic modeling of efferocytosis-related subtypes in glioblastoma

BACKGROUND

Efferocytosis is a biological process in which phagocytes remove apoptotic cells and vesicles from tissues. This process is initiated by the release of inflammatory mediators from apoptotic cells and plays a crucial role in resolving inflammation. The signals associated with efferocytosis have been found to regulate the inflammatory response and the tumor microenvironment (TME), which promotes the immune escape of tumor cells. However, the role of efferocytosis in glioblastoma multiforme (GBM) is not well understood and requires further investigation.

METHODS

In this study, we conducted a comprehensive analysis of 22 efferocytosis-related genes (ERGs) by searching for studies related to efferocytosis. Using bulk RNA-Seq and single-cell sequencing data, we analyzed the expression and mutational characteristics of these ERGs. By using an unsupervised clustering algorithm, we obtained ERG clusters from 549 GBM patients and evaluated the immune infiltration characteristics of each cluster. We then identified differential genes (DEGs) in the two ERG clusters and classified GBM patients into different gene clusters using univariate cox analysis and unsupervised clustering algorithms. Finally, we utilized the Boruta algorithm to screen for prognostic genes and reduce dimensionality, and the PCA algorithm was applied to create a novel efferocytosis-related scoring system.

RESULTS

Differential expression of ERGs in glioma cell lines and normal cells was analyzed by rt-PCR. Cell function experiments, on the other hand, validated TIMD4 as a tumor risk factor in GBM. We found that different ERG clusters and gene clusters have distinct prognostic and immune infiltration profiles. The ERG signature we developed provides insight into the tumor microenvironment of GBM. Patients with lower ERG scores have a better survival rate and a higher likelihood of benefiting from immunotherapy.

CONCLUSIONS

Our novel efferocytosis-related signature has the potential to be used in clinical practice for risk stratification of GBM patients and for selecting individuals who are likely to respond to immunotherapy. This can help clinicians design appropriate targeted therapies before initiating clinical treatment.

Unraveling the role of disulfidptosis-related LncRNAs in colon cancer: a prognostic indicator for immunotherapy response, chemotherapy sensitivity, and insights into cell death mechanisms

Background: Colon cancer, a prevalent and deadly malignancy worldwide, ranks as the third leading cause of cancer-related mortality. Disulfidptosis stress triggers a unique form of programmed cell death known as disulfidoptosis, characterized by excessive intracellular cystine accumulation. This study aimed to establish reliable bioindicators based on long non-coding RNAs (LncRNAs) associated with disulfidptosis-induced cell death, providing novel insights into immunotherapeutic response and prognostic assessment in patients with colon adenocarcinoma (COAD). Methods: Univariate Cox proportional hazard analysis and Lasso regression analysis were performed to identify differentially expressed genes strongly associated with prognosis. Subsequently, a multifactorial model for prognostic risk assessment was developed using multiple Cox proportional hazard regression. Furthermore, we conducted comprehensive evaluations of the characteristics of disulfidptosis response-related LncRNAs, considering clinicopathological features, tumor microenvironment, and chemotherapy sensitivity. The expression levels of prognosis-related genes in COAD patients were validated using quantitative real-time fluorescence PCR (qRT-PCR). Additionally, the role of ZEB1-SA1 in colon cancer was investigated through CCK8 assays, wound healing experiment and transwell experiments. Results: disulfidptosis response-related LncRNAs were identified as robust predictors of COAD prognosis. Multifactorial analysis revealed that the risk score derived from these LncRNAs served as an independent prognostic factor for COAD. Patients in the low-risk group exhibited superior overall survival (OS) compared to those in the high-risk group. Accordingly, our developed Nomogram prediction model, integrating clinical characteristics and risk scores, demonstrated excellent prognostic efficacy. In vitro experiments demonstrated that ZEB1-SA1 promoted the proliferation and migration of COAD cells. Conclusion: Leveraging medical big data and artificial intelligence, we constructed a prediction model for disulfidptosis response-related LncRNAs based on the TCGA-COAD cohort, enabling accurate prognostic prediction in colon cancer patients. The implementation of this model in clinical practice can facilitate precise classification of COAD patients, identification of specific subgroups more likely to respond favorably to immunotherapy and chemotherapy, and inform the development of personalized treatment strategies for COAD patients based on scientific evidence.

Identifying endoplasmic reticulum stress-related molecular subtypes and prognostic model for predicting the immune landscape and therapy response in pancreatic cancer

Endoplasmic reticulum stress (ERS) is caused by the accumulation of intracellular misfolded or unfolded proteins and is associated with cancer development. In this study, pan-cancer analysis revealed complex genetic variations, including copy number variation, methylation, and somatic mutations for ERS-related genes (ERGs) in 33 kinds of cancer. Consensus clustering divided pancreatic cancer (PC) patients from TCGA and GEO databases into two ERS-related subtypes: ERGcluster A and B. Compared with ERGcluster A, ERGcluster B had a more active ERS state and worse prognosis. Subsequently, the ERS-related prognostic model was established to quantify the ERS score for a single sample. The patient with a low ERS score had remarkably longer survival times. ssGSEA and CIBERSORT algorithms revealed that activated B cells and CD8+ T cells had higher infiltration in the low ERS score group, but higher infiltration of activated CD4+ T cells, activated dendritic cells, macrophages, and neutrophils in the high ERS score group. Drug sensitivity analysis indicated the low ERS score group had a better response to gemcitabine, paclitaxel, 5-fluorouracil, oxaliplatin, and irinotecan. RT-qPCR validated that MET, MUC16, and KRT7 in the model had higher expression levels in pancreatic tumour tissues. Single-cell analysis further revealed that MET, MUC16, and KRT7 were mainly expressed in cancer cells in PC tumour microenvironment. In all, we first constructed the ERS-related molecular subtypes and prognostic model in PC. Our research highlighted the vital role of ERS in PC and contributed to further research on molecular mechanisms and novel therapeutic strategies for PC in the future.

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.

Molecular characterization of cancer-intrinsic immune evasion genes indicates prognosis and tumour microenvironment infiltration in osteosarcoma

Cancer-intrinsic immune evasion (IE) to cells is a critical factor in tumour growth and progression, yet the molecular characterization of IE genes (IEGs) in osteosarcoma remains underexplored. In this study, 85 osteosarcoma patients were comprehensively analyzed based on 182 IEGs, leading to the identification of two IE clusters linked to distinct biological processes and clinical outcomes. In addition, two IE clusters demonstrated diverse immune cell infiltration patterns, with IEGcluster A displaying increased levels compared to IEGcluster B. Moreover, an IE score was identified as an independent prognostic factor and nomogram may serve as a practical tool for the individual prognostic evaluation of patients with osteosarcoma. Finally, GBP1, a potential biomarker with high expression in osteosarcoma was identified. The findings of this study highlight the presence of two IE clusters, each associated with differing patient outcomes and immune infiltration properties. The IE score may serve to assess individual patient IE characteristics, enhance comprehension of immune features, and guide more efficacious treatment approaches.

Integrated single-cell and bulk RNA sequencing revealed the molecular characteristics and prognostic roles of neutrophils in pancreatic cancer

Pancreatic cancer, one of the most prevalent tumors of the digestive system, has a dismal prognosis. Cancer of the pancreas is distinguished by an inflammatory tumor microenvironment rich in fibroblasts and different immune cells. Neutrophils are important immune cells that infiltrate the microenvironment of pancreatic cancer tumors. The purpose of this work was to examine the complex mechanism by which neutrophils influence the carcinogenesis and development of pancreatic cancer and to construct a survival prediction model based on neutrophil marker genes. We incorporated the GSE111672 dataset, comprising RNA expression data from 27,000 cells obtained from 3 patients with PC, and conducted single-cell data analysis. Thorough investigation of pancreatic cancer single-cell RNA sequencing data found 350 neutrophil marker genes. Using The Cancer Genome Atlas (TCGA), GSE28735, GSE62452, GSE57495, and GSE85916 datasets to gather pancreatic cancer tissue transcriptome data, and consistent clustering was used to identify two categories for analyzing the influence of neutrophils on pancreatic cancer. Using the Random Forest algorithm and Cox regression analysis, a survival prediction model for pancreatic cancer was developed, the model showed independent performance for survival prognosis, clinic pathological features, immune infiltration, and drug sensitivity. Multivariate Cox analysis findings revealed that the risk scores derived from predictive models is independent prognostic markers for pancreatic patients. In conclusion, based on neutrophil marker genes, this research created a molecular typing and prognostic grading system for pancreatic cancer, this system was very accurate in predicting the prognosis, tumor immune microenvironment status, and pharmacological treatment responsiveness of pancreatic cancer patients.

Unveiling the role of regulatory T cells in the tumor microenvironment of pancreatic cancer through single-cell transcriptomics and in vitro experiments

Background

In order to investigate the impact of Treg cell infiltration on the immune response against pancreatic cancer within the tumor microenvironment (TME), and identify crucial mRNA markers associated with Treg cells in pancreatic cancer, our study aims to delve into the role of Treg cells in the anti-tumor immune response of pancreatic cancer.

Methods

The ordinary transcriptome data for this study was sourced from the GEO and TCGA databases. It was analyzed using single-cell sequencing analysis and machine learning. To assess the infiltration level of Treg cells in pancreatic cancer tissues, we employed the CIBERSORT method. The identification of genes most closely associated with Treg cells was accomplished through the implementation of weighted gene co-expression network analysis (WGCNA). Our analysis of single-cell sequencing data involved various quality control methods, followed by annotation and advanced analyses such as cell trajectory analysis and cell communication analysis to elucidate the role of Treg cells within the pancreatic cancer microenvironment. Additionally, we categorized the Treg cells into two subsets: Treg1 associated with favorable prognosis, and Treg2 associated with poor prognosis, based on the enrichment scores of the key genes. Employing the hdWGCNA method, we analyzed these two subsets to identify the critical signaling pathways governing their mutual transformation. Finally, we conducted PCR and immunofluorescence staining in vitro to validate the identified key genes.

Results

Based on the results of immune infiltration analysis, we observed significant infiltration of Treg cells in the pancreatic cancer microenvironment. Subsequently, utilizing the WGCNA and machine learning algorithms, we ultimately identified four Treg cell-related genes (TRGs), among which four genes exhibited significant correlations with the occurrence and progression of pancreatic cancer. Among them, CASP4, TOB1, and CLEC2B were associated with poorer prognosis in pancreatic cancer patients, while FYN showed a correlation with better prognosis. Notably, significant differences were found in the HIF-1 signaling pathway between Treg1 and Treg2 cells identified by the four genes. These conclusions were further validated through in vitro experiments.

Conclusion

Treg cells played a crucial role in the pancreatic cancer microenvironment, and their presence held a dual significance. Recognizing this characteristic was vital for understanding the limitations of Treg cell-targeted therapies. CASP4, FYN, TOB1, and CLEC2B exhibited close associations with infiltrating Treg cells in pancreatic cancer, suggesting their involvement in Treg cell functions. Further investigation was warranted to uncover the mechanisms underlying these associations. Notably, the HIF-1 signaling pathway emerged as a significant pathway contributing to the duality of Treg cells. Targeting this pathway could potentially revolutionize the existing treatment approaches for pancreatic cancer.

Identification and validation of disulfidptosis-associated molecular clusters in non-alcoholic fatty liver disease

Objective: Non-alcoholic fatty liver disease (NAFLD) is the most prevalent liver disease in the world, and its pathogenesis is not fully understood. Disulfidptosis is the most recently reported form of cell death and may be associated with NAFLD progression. Our study aimed to explore the molecular clusters associated with disulfidptosis in NAFLD and to construct a predictive model. Methods: First, we analyzed the expression profile of the disulfidptosis regulators and immune characteristics in NAFLD. Using 104 NAFLD samples, we investigated molecular clusters based on differentially expressed disulfidptosis-related genes, along with the related immune cell infiltration. Cluster-specific differentially expressed genes were then identified by using the WGCNA method. We also evaluated the performance of four machine learning models before choosing the optimal machine model for diagnosis. Nomogram, calibration curves, decision curve analysis, and external datasets were used to confirm the prediction effectiveness. Finally, the expression levels of the biomarkers were assessed in a mouse model of a high-fat diet. Results: Two differentially expressed DRGs were identified between healthy and NAFLD patients. We revealed the expression profile of DRGs in NAFLD and the correlation with 22 immune cells. In NAFLD, two clusters of molecules connected to disulfidptosis were defined. Significant immunological heterogeneity was shown by immune infiltration analysis among the various clusters. A significant amount of immunological infiltration was seen in Cluster 1. Functional analysis revealed that Cluster 1 differentially expressed genes were strongly linked to energy metabolism and immune control. The highest discriminatory performance was demonstrated by the SVM model, which had a higher area under the curve, relatively small residual and root mean square errors. Nomograms, calibration curves, and decision curve analyses were used to show how accurate the prediction of NAFLD was. Further analysis revealed that the expression of three model-related genes was significantly associated with the level of multiple immune cells. In animal experiments, the expression trends of DDO, FRK and TMEM19 were consistent with the results of bioinformatics analysis. Conclusion: This study systematically elucidated the complex relationship between disulfidptosis and NAFLD and developed a promising predictive model to assess the risk of disease in patients with disulfidptosis subtypes and NAFLD.

Revealing the role of regulatory T cells in the tumor microenvironment of lung adenocarcinoma: a novel prognostic and immunotherapeutic signature

Background

Regulatory T cells (Tregs), are a key class of cell types in the immune system. In the tumor microenvironment (TME), the presence of Tregs has important implications for immune response and tumor development. Relatively little is known about the role of Tregs in lung adenocarcinoma (LUAD).

Methods

Tregs were identified using but single-cell RNA sequencing (scRNA-seq) analysis and interactions between Tregs and other cells in the TME were investigated. Next, we used multiple bulk RNA-seq datasets to construct risk models based on marker genes of Tregs and explored differences in prognosis, mutational landscape, immune cell infiltration and immunotherapy between high- and low-risk groups, and finally, qRT-PCR and cell function experiments were performed to validate the model genes.

Results

The cellchat analysis showed that MIF-(CD74+CXCR4) pairs play a key role in the interaction of Tregs with other cell subpopulations, and the Tregs-associated signatures (TRAS) could well classify multiple LUAD cohorts into high- and low-risk groups. Immunotherapy may offer greater potential benefits to the low-risk group, as indicated by their superior survival, increased infiltration of immune cells, and heightened expression of immune checkpoints. Finally, the experiment verified that the model genes LTB and PTTG1 were relatively highly expressed in cancer tissues, while PTPRC was relatively highly expressed in paracancerous tissues. Colony Formation assay confirmed that knockdown of PTTG1 reduced the proliferation ability of LUAD cells.

Conclusion

TRAS were constructed using scRNA-seq and bulk RNA-seq to distinguish patient risk subgroups, which may provide assistance in the clinical management of LUAD patients.

The combined signatures of telomere and immune cell landscape provide a prognostic and therapeutic biomarker in glioma

Background

Gliomas, the most prevalent primary malignant tumors of the central nervous system in adults, exhibit slow growth in lower-grade gliomas (LGG). However, the majority of LGG cases progress to high-grade gliomas, posing challenges for prognostication. The tumor microenvironment (TME), characterized by telomere-related genes and immune cell infiltration, strongly influences glioma growth and therapeutic response. Therefore, our objective was to develop a Telomere-TME (TM-TME) classifier that integrates telomere-related genes and immune cell landscape to assess prognosis and therapeutic response in glioma.

Methods

This study encompassed LGG patients from the TCGA and CCGA databases. TM score and TME score were derived from the expression signatures of telomere-related genes and the presence of immune cells in LGG, respectively. The TM-TME classifier was established by combining TM and TME scores to effectively predict prognosis. Subsequently, we conducted Kaplan-Meier survival estimation, univariate Cox regression analysis, and receiver operating characteristic curves to validate the prognostic prediction capacity of the TM-TME classifier across multiple cohorts. Gene Ontology (GO) analysis, biological processes, and proteomaps were performed to annotate the functional aspects of each subgroup and visualize the cellular signaling pathways.

Results

The TM_low+TME_high subgroup exhibited superior prognosis and therapeutic response compared to other subgroups (P<0.001). This finding could be attributed to distinct tumor somatic mutations and cancer cellular signaling pathways. GO analysis indicated that the TM_low+TME_high subgroup is associated with the neuronal system and modulation of chemical synaptic transmission. Conversely, the TM_high+TME_low subgroup showed a strong association with cell cycle and DNA metabolic processes. Furthermore, the classifier significantly differentiated overall survival in the TCGA LGG cohort and served as an independent prognostic factor for LGG patients in both the TCGA cohort (P<0.001) and the CGGA cohort (P<0.001).

Conclusion

Overall, our findings underscore the significance of the TM-TME classifier in predicting prognosis and immune therapeutic response in glioma, shedding light on the complex immune landscape within each subgroup. Additionally, our results suggest the potential of integrating risk stratification with precision therapy for LGG.

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.

Unveiling efferocytosis-related signatures through the integration of single-cell analysis and machine learning: a predictive framework for prognosis and immunotherapy response in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) represents a prominent gastrointestinal malignancy with a grim clinical outlook. In this regard, the discovery of novel early biomarkers holds substantial promise for ameliorating HCC-associated mortality. Efferocytosis, a vital immunological process, assumes a central position in the elimination of apoptotic cells. However, comprehensive investigations exploring the role of efferocytosis-related genes (EFRGs) in HCC are sparse, and their regulatory influence on HCC immunotherapy and targeted drug interventions remain poorly understood.

Methods

RNA sequencing data and clinical characteristics of HCC patients were acquired from the TCGA database. To identify prognostically significant genes in HCC, we performed the limma package and conducted univariate Cox regression analysis. Subsequently, machine learning algorithms were employed to identify hub genes. To assess the immunological landscape of different HCC subtypes, we employed the CIBERSORT algorithm. Furthermore, single-cell RNA sequencing (scRNA-seq) was utilized to investigate the expression levels of ERFGs in immune cells and to explore intercellular communication within HCC tissues. The migratory capacity of HCC cells was evaluated using CCK-8 assays, while drug sensitivity prediction reliability was determined through wound-healing assays.

Results

We have successfully identified a set of nine genes, termed EFRGs, that hold significant potential for the establishment of a hepatocellular carcinoma-specific prognostic model. Furthermore, leveraging the individual risk scores derived from this model, we were able to stratify patients into two distinct risk groups, unveiling notable disparities in terms of immune infiltration patterns and response to immunotherapy. Notably, the model's capacity to accurately predict drug responses was substantiated through comprehensive experimental investigations, encompassing wound-healing assay, and CCK8 experiments conducted on the HepG2 and Huh7 cell lines.

Conclusions

We constructed an EFRGs model that serves as valuable tools for prognostic assessment and decision-making support in the context of immunotherapy and chemotherapy.

Prognostic model for hepatocellular carcinoma based on anoikis-related genes: immune landscape analysis and prediction of drug sensitivity

Background

Hepatocellular carcinoma (HCC) represents a complex ailment characterized by an unfavorable prognosis in advanced stages. The involvement of immune cells in HCC progression is of significant importance. Moreover, metastasis poses a substantial impediment to enhanced prognostication for HCC patients, with anoikis playing an indispensable role in facilitating the distant metastasis of tumor cells. Nevertheless, limited investigations have been conducted regarding the utilization of anoikis factors for predicting HCC prognosis and assessing immune infiltration. This present study aims to identify hepatocellular carcinoma-associated anoikis-related genes (ANRGs), establish a robust prognostic model for HCC, and delineate distinct immune characteristics based on the anoikis signature. Cell migration and cytotoxicity experiments were performed to validate the accuracy of the ANRGs model.

Methods

Consensus clustering based on ANRGs was employed in this investigation to categorize HCC samples obtained from both TCGA and Gene Expression Omnibus (GEO) cohorts. To assess the differentially expressed genes, Cox regression analysis was conducted, and subsequently, prognostic gene signatures were constructed using LASSO-Cox methodology. External validation was performed at the International Cancer Genome Conference. The tumor microenvironment (TME) was characterized utilizing ESTIMATE and CIBERSORT algorithms, while machine learning techniques facilitated the identification of potential target drugs. The wound healing assay and CCK-8 assay were employed to evaluate the migratory capacity and drug sensitivity of HCC cell lines, respectively.

Results

Utilizing the TCGA-LIHC dataset, we devised a nomogram integrating a ten-gene signature with diverse clinicopathological features. Furthermore, the discriminative potential and clinical utility of the ten-gene signature and nomogram were substantiated through ROC analysis and DCA. Subsequently, we devised a prognostic framework leveraging gene expression data from distinct risk cohorts to predict the drug responsiveness of HCC subtypes.

Conclusion

In this study, we have established a promising HCC prognostic ANRGs model, which can serve as a valuable tool for clinicians in selecting targeted therapeutic drugs, thereby improving overall patient survival rates. Additionally, this model has also revealed a strong connection between anoikis and immune cells, providing a potential avenue for elucidating the mechanisms underlying immune cell infiltration regulated by anoikis.

Identification of key biomarkers in ischemic stroke: single-cell sequencing and weighted co-expression network analysis

PURPOSE

At present, there is a lack of accurate early diagnostic markers for ischemic stroke.

METHODS

By using dimensionality reduction cluster analysis, differential expression analysis, weighted co-expression network analysis, protein-protein interaction network analysis, cell heterogeneity and key pathogenic genes were identified in ischemic stroke. Immunomicroenvironment analysis was used to explore the immune landscape and immune associations of key genes in ischemic stroke. The analysis platform we use is R software (version 4.0.5). PCR experiments were used to verify the expression of key genes.

RESULTS

Single cell sequencing data in ischemic stroke can be annotated as fibroblast cells, pre-B cell CD34, neutrophils cells, bone marrow (BM), keratinocytes, macrophage, neurons and mesenchymal stem cells (MSC). By the intersection of differential expression analysis and WGCNA analysis, 385 genes were obtained. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these genes were highly correlated with multiple functions and pathways. Protein-protein interaction network analysis revealed that MRPS11 and MRPS12 were key genes, both of which were down-regulated in ischemic stroke. The Pseudo-time series analysis found that the expression of MRPS12 decreased gradually with the differentiation of pre-B cell CD34 cells in ischemic stroke, suggesting that the downregulation of MRPS12 expression may play an important role in ischemic stroke. At last, PCR showed that MRPS11 and MRPS12 were significantly down-regulated in peripheral blood of patients with ischemic stroke.

CONCLUSIONS

Our study provides a reference for the study of pathogenesis and key targets of ischemic stroke.

An immune-related gene prognostic index for predicting prognosis in patients with colorectal cancer

Background

Colorectal cancer (CRC) is one of the most common solid malignant burdens worldwide. Cancer immunology and immunotherapy have become fundamental areas in CRC research and treatment. Currently, the method of generating Immune-Related Gene Prognostic Indices (IRGPIs) has been found to predict patient prognosis as an immune-related prognostic biomarker in a variety of tumors. However, their role in patients with CRC remains mostly unknown. Therefore, we aimed to establish an IRGPI for prognosis evaluation in CRC.

Methods

RNA-sequencing data and clinical information of CRC patients were retrieved from The Cancer Genome Atlas (TCGA) and The Gene Expression Omnibus (GEO) databases as training and validation sets, respectively. Immune-related gene data was obtained from the ImmPort and InnateDB databases. The weighted gene co-expression network analysis (WGCNA) was used to identify hub immune-related genes. An IRGPI was then constructed using Cox regression methods. Based on the median risk score of IRGPI, patients could be divided into high-risk and low-risk groups. To further investigate the immunologic differences, Gene set variation analysis (GSVA) studies were conducted. In addition, immune cell infiltration and related functional analysis were used to identify the differential immune cell subsets and related functional pathways.

Results

We identified 49 immune-related genes associated with the prognosis of CRC, 17 of which were selected for an IRGPI. The IRGPI model significantly differentiates the survival rates of CRC patients in the different groups. The IRGPI as an independent prognostic factor significantly correlates with clinico-pathological factors such as age and tumor stage. Furthermore, we developed a nomogram to improve the clinical utility of the IRGPI score. Immuno-correlation analysis in different IRGPI groups revealed distinct immune cell infiltration (CD4⁺ T cells resting memory) and associated pathways (macrophages, Type I IFNs responses, iDCs.), providing new insights into the tumor microenvironment. At last, drug sensitivity analysis revealed that the high-risk IRGPI group was sensitive to 11 and resistant to 15 drugs.

Conclusion

Our study established a promising immune-related risk model for predicting survival in CRC patients. This could help to better understand the correlation between immunity and the prognosis of CRC providing a new perspective for personalized treatment of CRC.

Identification of peripheral blood immune infiltration signatures and construction of monocyte-associated signatures in ovarian cancer and Alzheimer's disease using single-cell sequencing

Background

Ovarian cancer (OC) is a common tumor of the female reproductive system, while Alzheimer's disease (AD) is a prevalent neurodegenerative disease that primarily affects cognitive function in the elderly. Monocytes are immune cells in the blood that can enter tissues and transform into macrophages, thus participating in immune and inflammatory responses. Overall, monocytes may play an important role in Alzheimer's disease and ovarian cancer.

Methods

The CIBERSORT algorithm results indicate a potential crucial role of monocytes/macrophages in OC and AD. To identify monocyte marker genes, single-cell RNA-seq data of peripheral blood mononuclear cells (PBMCs) from OC and AD patients were analyzed. Enrichment analysis of various cell subpopulations was performed using the "irGSEA" R package. The estimation of cell cycle was conducted with the "tricycle" R package, and intercellular communication networks were analyzed using "CellChat". For 134 monocyte-associated genes (MRGs), bulk RNA-seq data from two diseased tissues were obtained. Cox regression analysis was employed to develop risk models, categorizing patients into high-risk (HR) and low-risk (LR) groups. The model's accuracy was validated using an external GEO cohort. The different risk groups were evaluated in terms of immune cell infiltration, mutational status, signaling pathways, immune checkpoint expression, and immunotherapy. To identify characteristic MRGs in AD, two machine learning algorithms, namely random forest and support vector machine (SVM), were utilized.

Results

Based on Cox regression analysis, a risk model consisting of seven genes was developed in OC, indicating a better prognosis for patients in the LR group. The LR group had a higher tumor mutation burden, immune cell infiltration abundance, and immune checkpoint expression. The results of the TIDE algorithm and the IMvigor210 cohort showed that the LR group was more likely to benefit from immunotherapy. Finally, ZFP36L1 and AP1S2 were identified as characteristic MRGs affecting OC and AD progression.

Conclusion

The risk profile containing seven genes identified in this study may help further guide clinical management and targeted therapy for OC. ZFP36L1 and AP1S2 may serve as biomarkers and new therapeutic targets for patients with OC and AD.

Innovative breakthroughs facilitated by single-cell multi-omics: manipulating natural killer cell functionality correlates with a novel subcategory of melanoma cells

Background

Melanoma is typically regarded as the most dangerous form of skin cancer. Although surgical removal of in situ lesions can be used to effectively treat metastatic disease, this condition is still difficult to cure. Melanoma cells are removed in great part due to the action of natural killer (NK) and T cells on the immune system. Still, not much is known about how the activity of NK cell-related pathways changes in melanoma tissue. Thus, we performed a single-cell multi-omics analysis on human melanoma cells in this study to explore the modulation of NK cell activity.

Materials and methods

Cells in which mitochondrial genes comprised > 20% of the total number of expressed genes were removed. Gene ontology (GO), gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and AUCcell analysis of differentially expressed genes (DEGs) in melanoma subtypes were performed. The CellChat package was used to predict cell-cell contact between NK cell and melanoma cell subtypes. Monocle program analyzed the pseudotime trajectories of melanoma cells. In addition, CytoTRACE was used to determine the recommended time order of melanoma cells. InferCNV was utilized to calculate the CNV level of melanoma cell subtypes. Python package pySCENIC was used to assess the enrichment of transcription factors and the activity of regulons in melanoma cell subtypes. Furthermore, the cell function experiment was used to confirm the function of TBX21 in both A375 and WM-115 melanoma cell lines.

Results

Following batch effect correction, 26,161 cells were separated into 28 clusters and designated as melanoma cells, neural cells, fibroblasts, endothelial cells, NK cells, CD4+ T cells, CD8+ T cells, B cells, plasma cells, monocytes and macrophages, and dendritic cells. A total of 10137 melanoma cells were further grouped into seven subtypes, i.e., C0 Melanoma BIRC7, C1 Melanoma CDH19, C2 Melanoma EDNRB, C3 Melanoma BIRC5, C4 Melanoma CORO1A, C5 Melanoma MAGEA4, and C6 Melanoma GJB2. The results of AUCell, GSEA, and GSVA suggested that C4 Melanoma CORO1A may be more sensitive to NK and T cells through positive regulation of NK and T cell-mediated immunity, while other subtypes of melanoma may be more resistant to NK cells. This suggests that the intratumor heterogeneity (ITH) of melanoma-induced activity and the difference in NK cell-mediated cytotoxicity may have caused NK cell defects. Transcription factor enrichment analysis indicated that TBX21 was the most important TF in C4 Melanoma CORO1A and was also associated with M1 modules. In vitro experiments further showed that TBX21 knockdown dramatically decreases melanoma cells' proliferation, invasion, and migration.

Conclusion

The differences in NK and T cell-mediated immunity and cytotoxicity between C4 Melanoma CORO1A and other melanoma cell subtypes may offer a new perspective on the ITH of melanoma-induced metastatic activity. In addition, the protective factors of skin melanoma, STAT1, IRF1, and FLI1, may modulate melanoma cell responses to NK or T cells.

Efferocytosis signatures as prognostic markers for revealing immune landscape and predicting immunotherapy response in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is a highly lethal liver cancer with late diagnosis; therefore, the identification of new early biomarkers could help reduce mortality. Efferocytosis, a process in which one cell engulfs another cell, including macrophages, dendritic cells, NK cells, etc., plays a complex role in tumorigenesis, sometimes promoting and sometimes inhibiting tumor development. However, the role of efferocytosis-related genes (ERGs) in HCC progression has been poorly studied, and their regulatory effects in HCC immunotherapy and drug targeting have not been reported. Methods: We downloaded efferocytosis-related genes from the Genecards database and screened for ERGs that showed significant expression changes between HCC and normal tissues and were associated with HCC prognosis. Machine learning algorithms were used to study prognostic gene features. CIBERSORT and pRRophetic R packages were used to evaluate the immune environment of HCC subtypes and predict treatment response. CCK-8 experiments conducted on HCC cells were used to assess the reliability of drug sensitivity prediction. Results: We constructed a prognostic prediction model composed of six genes, and the ROC curve showed good predictive accuracy of the risk model. In addition, two ERG-related subgroups in HCC showed significant differences in tumor immune landscape, immune response, and prognostic stratification. The CCK-8 experiment conducted on HCC cells confirmed the reliability of drug sensitivity prediction. Conclusion: Our study emphasizes the importance of efferocytosis in HCC progression. The risk model based on efferocytosis-related genes developed in our study provides a novel precision medicine approach for HCC patients, allowing clinicians to customize treatment plans based on unique patient characteristics. The results of our investigation carry noteworthy implications for the development of individualized treatment approaches involving immunotherapy and chemotherapy, thereby potentially facilitating the realization of personalized and more efficacious therapeutic interventions for HCC.

Multi-omic analysis of glycolytic signatures: exploring the predictive significance of heterogeneity and stemness in immunotherapy response and outcomes in hepatocellular carcinoma

Background: Hepatocellular carcinoma (HCC) is a global health challenge with complex pathophysiology, characterized by high mortality rates and poor early detection due to significant tumor heterogeneity. Stemness significantly contributes to the heterogeneity of HCC tumors, and glycolysis is crucial for maintaining stemness. However, the predictive significance of glycolysis-related metabolic genes (GMGs) in HCC remains unknown. Therefore, this study aimed to identify critical GMGs and establish a reliable model for HCC prognosis. Methods: GMGs associated with prognosis were identified by evaluating genes with notable expression changes between HCC and normal tissues retrieved from the MsigDB database. Prognostic gene characteristics were established using univariate and multivariate Cox regression studies for prognosis prediction and risk stratification. The "CIBERSORT" and "pRRophetic" R packages were respectively used to evaluate the immunological environment and predict treatment response in HCC subtypes. The HCC stemness score was obtained using the OCLR technique. The precision of drug sensitivity prediction was evaluated using CCK-8 experiments performed on HCC cells. The miagration and invasion ability of HCC cell lines with different riskscores were assessed using Transwell and wound healing assays. Results: The risk model based on 10 gene characteristics showed high prediction accuracy as indicated by the receiver operating characteristic (ROC) curves. Moreover, the two GMG-related subgroups showed considerable variation in the risk of HCC with respect to tumor stemness, immune landscape, and prognostic stratification. The in vitro validation of the model's ability to predict medication response further demonstrated its reliability. Conclusion: Our study highlights the importance of stemness variability and inter-individual variation in determining the HCC risk landscape. The risk model we developed provides HCC patients with a novel method for precision medicine that enables clinical doctors to customize treatment plans based on unique patient characteristics. Our findings have significant implications for tailored immunotherapy and chemotherapy methods, and may pave the way for more personalized and effective treatment strategies for HCC.

An amino acid metabolism-based seventeen-gene signature correlates with the clinical outcome and immune features in pancreatic cancer

Background: Pancreatic cancer is an aggressive tumor with a low 5-year survival rate and primary resistance to most therapy. Amino acid (AA) metabolism is highly correlated with tumor growth, crucial to the aggressive biological behavior of pancreatic cancer; nevertheless, the comprehensive predictive significance of genes that regulate AA metabolism in pancreatic cancer remains unknown. Methods: The mRNA expression data downloaded from The Cancer Genome Atlas (TCGA) were derived as the training cohort, and the GSE57495 cohort from Gene Expression Omnibus (GEO) database was applied as the validation cohort. Random survival forest (RSF) and the least absolute shrinkage and selection operator (LASSO) regression analysis were employed to screen genes and construct an AA metabolism-related risk signature (AMRS). Kaplan-Meier analysis and receiver operating characteristic (ROC) curve were performed to assess the prognostic value of AMRS. We performed genomic alteration analysis and explored the difference in tumor microenvironment (TME) landscape associated with KRAS and TP53 mutation in both high- and low-AMRS groups. Subsequently, the relationships between AMRS and immunotherapy and chemotherapy sensitivity were evaluated. Results: A 17-gene AA metabolism-related risk model in the TCGA cohort was constructed according to RSF and LASSO. After stratifying patients into high- and low-AMRS groups based on the optimal cut-off value, we found that high-AMRS patients had worse overall survival (OS) in the training cohort (a median OS: 13.1 months vs. 50.1 months, p < 0.0001) and validation cohort (a median OS: 16.2 vs. 30.5 months, p = 1e-04). Genetic mutation analysis revealed that KRAS and TP53 were significantly more mutated in high-AMRS group, and patients with KRAS and TP53 alterations had significantly higher risk scores than those without. Based on the analysis of TME, low-AMRS group displayed significantly higher immune score and more enrichment of T Cell CD8⁺ cells. In addition, high-AMRS-group exhibited higher TMB and significantly lower tumor immune dysfunction and exclusion (TIDE) score and T Cells dysfunction score, which suggested a higher sensitive to immunotherapy. Moreover, high-AMRS group was also more sensitive to paclitaxel, cisplatin, and docetaxel. Conclusion: Overall, we constructed an AA-metabolism prognostic model, which provided a powerful prognostic predictor for the clinical treatment of pancreatic cancer.

Identification of high-risk factors associated with mortality at 1-, 3-, and 5-year intervals in gastric cancer patients undergoing radical surgery and immunotherapy: an 8-year multicenter retrospective analysis

Background

Combining immunotherapy with surgical intervention is a prevailing and radical therapeutic strategy for individuals afflicted with gastric carcinoma; nonetheless, certain patients exhibit unfavorable prognoses even subsequent to this treatment regimen. This research endeavors to devise a machine learning algorithm to recognize risk factors with a high probability of inducing mortality among patients diagnosed with gastric cancer, both prior to and during their course of treatment.

Methods

Within the purview of this investigation, a cohort of 1015 individuals with gastric cancer were incorporated, and 39 variables encompassing diverse features were recorded. To construct the models, we employed three distinct machine learning algorithms, specifically extreme gradient boosting (XGBoost), random forest (RF), and k-nearest neighbor algorithm (KNN). The models were subjected to internal validation through employment of the k-fold cross-validation technique, and subsequently, an external dataset was utilized to externally validate the models.

Results

In comparison to other machine learning algorithms employed, the XGBoost algorithm demonstrated superior predictive capacity regarding the risk factors that affect mortality after combination therapy in gastric cancer patients for a duration of one year, three years, and five years posttreatment. The common risk factors that significantly impacted patient survival during the aforementioned time intervals were identified as advanced age, tumor invasion, tumor lymph node metastasis, tumor peripheral nerve invasion (PNI), multiple tumors, tumor size, carcinoembryonic antigen (CEA) level, carbohydrate antigen 125 (CA125) level, carbohydrate antigen 72-4 (CA72-4) level, and H. pylori infection.

Conclusion

The XGBoost algorithm can assist clinicians in identifying pivotal prognostic factors that are of clinical significance and can contribute toward individualized patient monitoring and management.

A novel T-cell exhaustion-related feature can accurately predict the prognosis of OC patients

The phenomenon of T Cell exhaustion (TEX) entails a progressive deterioration in the functionality of T cells within the immune system during prolonged conflicts with chronic infections or tumors. In the context of ovarian cancer immunotherapy, the development, and outcome of treatment are closely linked to T-cell exhaustion. Hence, gaining an in-depth understanding of the features of TEX within the immune microenvironment of ovarian cancer is of paramount importance for the management of OC patients. To this end, we leveraged single-cell RNA data from OC to perform clustering and identify T-cell marker genes utilizing the Unified Modal Approximation and Projection (UMAP) approach. Through GSVA and WGCNA in bulk RNA-seq data, we identified 185 TEX-related genes (TEXRGs). Subsequently, we transformed ten machine learning algorithms into 80 combinations and selected the most optimal one to construct TEX-related prognostic features (TEXRPS) based on the mean C-index of the three OC cohorts. In addition, we explored the disparities in clinicopathological features, mutational status, immune cell infiltration, and immunotherapy efficacy between the high-risk (HR) and low-risk (LR) groups. Upon the integration of clinicopathological features, TEXRPS displayed robust predictive power. Notably, patients in the LR group exhibited a superior prognosis, higher tumor mutational load (TMB), greater immune cell infiltration abundance, and enhanced sensitivity to immunotherapy. Lastly, we verified the differential expression of the model gene CD44 using qRT-PCR. In conclusion, our study offers a valuable tool to guide clinical management and targeted therapy of OC.

FAM family gene prediction model reveals heterogeneity, stemness and immune microenvironment of UCEC

Background: Endometrial cancer (UCEC) is a highly heterogeneous gynecologic malignancy that exhibits variable prognostic outcomes and responses to immunotherapy. The Familial sequence similarity (FAM) gene family is known to contribute to the pathogenesis of various malignancies, but the extent of their involvement in UCEC has not been systematically studied. This investigation aimed to develop a robust risk profile based on FAM family genes (FFGs) to predict the prognosis and suitability for immunotherapy in UCEC patients. Methods: Using the TCGA-UCEC cohort from The Cancer Genome Atlas (TCGA) database, we obtained expression profiles of FFGs from 552 UCEC and 35 normal samples, and analyzed the expression patterns and prognostic relevance of 363 FAM family genes. The UCEC samples were randomly divided into training and test sets (1:1), and univariate Cox regression analysis and Lasso Cox regression analysis were conducted to identify the differentially expressed genes (FAM13C, FAM110B, and FAM72A) that were significantly associated with prognosis. A prognostic risk scoring system was constructed based on these three gene characteristics using multivariate Cox proportional risk regression. The clinical potential and immune status of FFGs were analyzed using CiberSort, SSGSEA, and tumor immune dysfunction and rejection (TIDE) algorithms. qRT-PCR and IHC for detecting the expression levels of 3-FFGs. Results: Three FFGs, namely, FAM13C, FAM110B, and FAM72A, were identified as strongly associated with the prognosis of UCEC and effective predictors of UCEC prognosis. Multivariate analysis demonstrated that the developed model was an independent predictor of UCEC, and that patients in the low-risk group had better overall survival than those in the high-risk group. The nomogram constructed from clinical characteristics and risk scores exhibited good prognostic power. Patients in the low-risk group exhibited a higher tumor mutational load (TMB) and were more likely to benefit from immunotherapy. Conclusion: This study successfully developed and validated novel biomarkers based on FFGs for predicting the prognosis and immune status of UCEC patients. The identified FFGs can accurately assess the prognosis of UCEC patients and facilitate the identification of specific subgroups of patients who may benefit from personalized treatment with immunotherapy and chemotherapy.

Identification of fatty acid-related subtypes, the establishment of a prognostic signature, and immune infiltration characteristics in lung adenocarcinoma

Abnormal fatty acid (FA) metabolism can change the inflammatory microenvironment and promote tumor progression and metastasis, however, the potential association between FA-related genes (FARGs) and lung adenocarcinoma (LUAD) is still unclear. In this study, we described the genetic and transcriptomic changes of FARGs in LUAD patients and identified two different FA subtypes, which were significantly correlated with overall survival and tumor microenvironment infiltrating cells in LUAD patients. In addition, the FA score was also constructed through the LASSO Cox to evaluate the FA dysfunction of each patient. Multivariate Cox analysis proved that the FA score was an independent predictor and created the FA score integrated nomogram, which offered a quantitative tool for clinical practice. The performance of the FA score has been substantiated in numerous datasets for its commendable accuracy in estimating overall survival in LUAD patients. The groups with high and low FA scores exhibited different mutation spectrums, copy number variations, enrichment pathways, and immune status. Noteworthy differences between the two groups in terms of immunophenoscore and Tumor Immune Dysfunction and Exclusion were observed, suggesting that the group with a low FA score was more responsive to immunotherapy, and similar results were also confirmed in the immunotherapy cohort. In addition, seven potential chemotherapeutic drugs related to FA score targeting were predicted. Ultimately, we ascertained that the attenuation of KRT6A expression impeded the proliferation, migration, and invasion of LUAD cell lines. In summary, this research offers novel biomarkers to facilitate prognostic forecasting and clinical supervision for individuals afflicted with LUAD.

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.

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.

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.

Identification of cuproptosis-related biomarkers in dilated cardiomyopathy and potential therapeutic prediction of herbal medicines

Background: Dilated cardiomyopathy (DCM) is one of the significant causes of heart failure, and the mechanisms of metabolic ventricular remodelling due to disturbances in energy metabolism are still poorly understood in cardiac pathology. Understanding the biological mechanisms of cuproptosis in DCM is critical for drug development. Methods: The DCM datasets were downloaded from Gene Expression Omnibus, their relationships with cuproptosis-related genes (CRGs) and immune signatures were analyzed. LASSO, RF, and SVM-RFE machine learning algorithms were used to identify signature genes and the eXtreme Gradient Boosting (XGBoost) model was used to assess diagnostic efficacy. Molecular clusters of CRGs were identified, and immune Infiltration analysis was performed. The WGCNA algorithm was used to identify specific genes in different clusters. In addition, AUCell was used to analyse the cuproptosis scores of different cell types in the scRNA-seq dataset. Finally, herbal medicines were predicted from an online database, and molecular docking and molecular dynamics simulations were used to support the confirmation of the potential of the selected compounds. Results: We identified dysregulated cuproptosis genes and activated immune responses between DCM and healthy controls. Two signature genes (FDX1, SLC31A1) were identified and performed well in an external validation dataset (AUC = 0.846). Two molecular clusters associated with cuproptosis were further defined in DCM, and immune infiltration analysis showed B-cell naive, Eosinophils, NK cells activated and T-cell CD4 memory resting is significant immune heterogeneity in the two clusters. AUCell analysis showed that cardiomyocytes had a high cuproposis score. In addition, 19 and 3 herbal species were predicted based on FDX1 and SLC31A1. Based on the molecular docking model, the natural compounds Rutin with FDX1 (-9.3 kcal/mol) and Polydatin with SLC31A1 (-5.5 kcal/mol) has high stability and molecular dynamics simulation studies further validated this structural stability. Conclusion: Our study systematically illustrates the complex relationship between cuproptosis and the pathological features of DCM and identifies two signature genes (FDX1 and SLC31A1) and two natural compounds (Rutin and Polydatin). This may enhance our diagnosis of the disease and facilitate the development of clinical treatment strategies for DCM.

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.