Expression of hub genes of endothelial cells in glioblastoma-A prognostic model for GBM patients integrating single-cell RNA sequencing and bulk RNA sequencing

Comprehensive machine learning-based integration develops a novel prognostic model for glioblastoma

In this study, we developed a new prognostic model for glioblastoma (GBM) based on an integrated machine learning algorithm. We used univariate Cox regression analysis to identify prognostic genes by combining six GBM cohorts. Based on the prognostic genes, 10 machine learning algorithms were integrated into 117 algorithm combinations, and the artificial intelligence prognostic signature (AIPS) with the greatest average C-index was chosen. The AIPS was compared with 10 previously published models by univariate Cox analysis and the C-index. We compared the differences in prognosis, tumor immune microenvironment (TIME), and immunotherapy sensitivity between the high and low AIPS score groups. The AIPS based on the random survival forest algorithm with the highest average C-index (0.868) was selected. Compared with the previous 10 prognostic models, our AIPS has the highest C-index. The AIPS was closely linked to the clinical features of GBM. We discovered that patients in the low score group had improved prognoses, a more active TIME, and were more sensitive to immunotherapy. Finally, we verified the expression of several key genes by western blotting and immunohistochemistry. We identified an ideal prognostic signature for GBM, which might provide new insights into stratified treatment approaches for GBM patients.

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

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

Exploring the prognostic value of BRMS1 + microglia based on single-cell anoikis regulator patterns in the immunologic microenvironment of GBM

BACKGROUND

Anoikis is a specialized form of programmed cell death induced by the loss of cell adhesion to the extracellular matrix (ECM). Acquisition of anoikis resistance is a significant marker for cancer cell invasion, metastasis, therapy resistance, and recurrence. Although current research has identified multiple factors that regulate anoikis resistance, the pathological mechanisms of anoikis-mediated tumor microenvironment (TME) in glioblastoma (GBM) remain largely unexplored.

METHODS

Utilizing single-cell RNA sequencing (scRNA-seq) data and employing non-negative matrix factorization (NMF), we identified and characterized TME cell clusters with distinct anoikis-associated gene signatures. Prognostic and therapeutic response analyses were conducted using TCGA and CGGA datasets to assess the clinical significance of different TME cell clusters. The spatial relationship between BRMS1 + microglia and tumor cells was inferred from spatial transcriptome RNA sequencing (stRNA-seq) data. To simulate the tumor immune microenvironment, co-culture experiments were performed with microglia (HMC3) and GBM cells (U118/U251), and microglia were transfected with a BRMS1 overexpression lentivirus. Western blot or ELISA were used to detect BRMS1, M2 macrophage-specific markers, PI3K/AKT signaling proteins, and apoptosis-related proteins. The proliferation and apoptosis capabilities of tumor cells were evaluated using CCK-8, colony formation, and apoptosis assays, while the invasive and migratory abilities of tumor cells were assessed using Transwell assays.

RESULTS

NMF-based analysis successfully identified CD8 + T cell and microglia cell clusters with distinct gene signature characteristics. Trajectory analysis, cell communication, and gene regulatory network analyses collectively indicated that anoikis-mediated TME cell clusters can influence tumor cell development through various mechanisms. Notably, BRMS1 + AP-Mic exhibited an M2 macrophage phenotype and had significant cell communication with malignant cells. Moreover, high expression of BRMS1 + AP-Mic in TCGA and CGGA datasets was associated with poorer survival outcomes, indicating its detrimental impact on immunotherapy. Upregulation of BRMS1 in microglia may lead to M2 macrophage polarization, activate the PI3K/AKT signaling pathway through SPP1/CD44-mediated cell interactions, inhibit tumor cell apoptosis, and promote tumor proliferation and invasion.

CONCLUSION

This pioneering study used NMF-based analysis to reveal the important predictive value of anoikis-regulated TME in GBM for prognosis and immunotherapeutic response. BRMS1 + microglial cells provide a new perspective for a deeper understanding of the immunosuppressive microenvironment of GBM and could serve as a potential therapeutic target in the future.

Spinal cord and brain atrophy patterns in neuromyelitis optica spectrum disorder and multiple sclerosis

BACKGROUND

Spinal cord and brain atrophy are common in neuromyelitis optica spectrum disorder (NMOSD) and relapsing-remitting multiple sclerosis (RRMS) but harbor distinct patterns accounting for disability and cognitive impairment.

METHODS

This study included 209 NMOSD and 304 RRMS patients and 436 healthy controls. Non-negative matrix factorization was used to parse differences in spinal cord and brain atrophy at subject level into distinct patterns based on structural MRI. The weights of patterns were obtained using a linear regression model and associated with Expanded Disability Status Scale (EDSS) and cognitive scores. Additionally, patients were divided into cognitive impairment (CI) and cognitive preservation (CP) groups.

RESULTS

Three patterns were observed in NMOSD: (1) Spinal Cord-Deep Grey Matter (SC-DGM) pattern was associated with high EDSS scores and decline of visuospatial memory function; (2) Frontal-Temporal pattern was associated with decline of language learning function; and (3) Cerebellum-Brainstem pattern had no observed association. Patients with CI had higher weights of SC-DGM pattern than CP group. Three patterns were observed in RRMS: (1) DGM pattern was associated with high EDSS scores, decreased information processing speed, and decreased language learning and visuospatial memory functions; (2) Frontal-Temporal pattern was associated with overall cognitive decline; and (3) Occipital pattern had no observed association. Patients with CI trended to have higher weights of DGM and Frontal-Temporal patterns than CP group.

CONCLUSION

This study estimated the heterogeneity of spinal cord and brain atrophy patterns in NMOSD and RRMS patients at individual level, and evaluated the clinical relevance of these patterns, which may contribute to stratifying participants for targeted therapy.

mRNA markers for survival prediction in glioblastoma multiforme patients: a systematic review with bioinformatic analyses

BACKGROUND

Glioblastoma multiforme (GBM) is a type of fast-growing brain glioma associated with a very poor prognosis. This study aims to identify key genes whose expression is associated with the overall survival (OS) in patients with GBM.

METHODS

A systematic review was performed using PubMed, Scopus, Cochrane, and Web of Science up to Journey 2024. Two researchers independently extracted the data and assessed the study quality according to the New Castle Ottawa scale (NOS). The genes whose expression was found to be associated with survival were identified and considered in a subsequent bioinformatic study. The products of these genes were also analyzed considering protein-protein interaction (PPI) relationship analysis using STRING. Additionally, the most important genes associated with GBM patients' survival were also identified using the Cytoscape 3.9.0 software. For final validation, GEPIA and CGGA (mRNAseq_325 and mRNAseq_693) databases were used to conduct OS analyses. Gene set enrichment analysis was performed with GO Biological Process 2023.

RESULTS

From an initial search of 4104 articles, 255 studies were included from 24 countries. Studies described 613 unique genes whose mRNAs were significantly associated with OS in GBM patients, of which 107 were described in 2 or more studies. Based on the NOS, 131 studies were of high quality, while 124 were considered as low-quality studies. According to the PPI network, 31 key target genes were identified. Pathway analysis revealed five hub genes (IL6, NOTCH1, TGFB1, EGFR, and KDR). However, in the validation study, only, the FN1 gene was significant in three cohorts.

CONCLUSION

We successfully identified the most important 31 genes whose products may be considered as potential prognosis biomarkers as well as candidate target genes for innovative therapy of GBM tumors.

Identification Of Endothelial Cell Immune-related Gene Signature for Lung Adenocarcinoma by Integrated Analysis of Single-cell and Bulk RNA Sequencing Data

Background: The role of endothelial cells in tumor progression is considerable, yet the effect of endothelial cell immune-related genes (EIRGs) is still unclear. This research aimed to scrutinize the prognostic value of EIRGs in lung adenocarcinoma (LUAD) and provide further insights into the abovementioned uncertainties. Methods: After single-cell RNA sequencing (scRNA-seq) samples were obtained from the Gene Expression Omnibus (GEO) database, they were integrated with bulk RNA sequencing data from The Cancer Genome Atlas (TCGA). Prognostic markers were determined and a prognostic model was developed. From this model, a nomogram was constructed. We analyzed the biological mechanism of the EIRGs in LUAD, including functional enrichment, tumor mutational burden (TMB), tumor microenvironment (TME) analyses and drug sensitivity. We validated the signature by validating the external cohort GSE31210 and RT-qPCR. Results: After analyzing the model constructed from eight EIRGs, we observed that high-risk group (HG) LUAD patients (a risk score exceeding 4.65) exhibited unfavorable outcomes according to Kaplan‒Meier survival curves. This outcome was confirmed by GSE31210. The nomogram based on the model demonstrated significant predictive value. HG was influenced primarily by steroid hormone biosynthesis and ECM receptor interactions. The TMB in HGs was greater than that in the LG. Analysis of drug sensitivity revealed the direction for individualized treatment for both risk cohorts. Variations in the expression of EIRGs have been confirmed via RT-qPCR in several LUAD cell lines. Conclusions: The prognostic model and nomogram above are valuable for determining the survival rate and treatment options for LUAD patients.

Exploration of functional relations among differentially co-expressed genes identifies regulators in glioblastoma

The conventional computational approaches to investigating a disease confront inherent constraints as they often need to improve in delving beyond protein functional associations and grasping their deeper contextual significance within the disease framework. Such context-specificity can be explored using clinical data by evaluating the change in interaction between the biological entities in different conditions by investigating the differential co-expression relationships. We believe that the integration and analysis of differential co-expression and the functional relationships, primarily focusing on the source nodes, will open novel insights about disease progression as the source proteins could trigger signaling cascades, mostly because they are transcription factors, cell surface receptors, or enzymes that respond instantly to a particular stimulus. A thorough contextual investigation of these nodes could lead to a helpful beginning point for identifying potential causal linkages and guiding subsequent scientific investigations to uncover mechanisms underlying observed associations. Our methodology includes functional protein-protein Interaction (PPI) data and co-expression information and filters functional linkages through a series of critical steps, culminating in the identification of a robust set of regulators. Our analysis identified eleven key regulators-AKT1, BRCA1, CAMK2G, CUL1, FGFR3, KIF3A, NUP210, PRKACB, RAB8A, RPS6KA2 and TGFB3-in glioblastoma. These regulators play a pivotal role in disease classification, cell growth control, and patient survivability and exhibit associations with immune infiltrations and disease hallmarks. This underscores the importance of assessing correlation towards causality in unraveling complex biological insights.

XRCC1: a potential prognostic and immunological biomarker in LGG based on systematic pan-cancer analysis

X-ray repair cross-complementation group 1 (XRCC1) is a pivotal contributor to base excision repair, and its dysregulation has been implicated in the oncogenicity of various human malignancies. However, a comprehensive pan-cancer analysis investigating the prognostic value, immunological functions, and epigenetic associations of XRCC1 remains lacking. To address this knowledge gap, we conducted a systematic investigation employing bioinformatics techniques across 33 cancer types. Our analysis encompassed XRCC1 expression levels, prognostic and diagnostic implications, epigenetic profiles, immune and molecular subtypes, Tumor Mutation Burden (TMB), Microsatellite Instability (MSI), immune checkpoints, and immune infiltration, leveraging data from TCGA, GTEx, CELL, Human Protein Atlas, Ualcan, and cBioPortal databases. Notably, XRCC1 displayed both positive and negative correlations with prognosis across different tumors. Epigenetic analysis revealed associations between XRCC1 expression and DNA methylation patterns in 10 cancer types, as well as enhanced phosphorylation. Furthermore, XRCC1 expression demonstrated associations with TMB and MSI in the majority of tumors. Interestingly, XRCC1 gene expression exhibited a negative correlation with immune cell infiltration levels, except for a positive correlation with M1 and M2 macrophages and monocytes in most cancers. Additionally, we observed significant correlations between XRCC1 and immune checkpoint gene expression levels. Lastly, our findings implicated XRCC1 in DNA replication and repair processes, shedding light on the precise mechanisms underlying its oncogenic effects. Overall, our study highlights the potential of XRCC1 as a prognostic and immunological pan-cancer biomarker, thereby offering a novel target for tumor 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.

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.

Combining single-cell sequencing and spatial transcriptome sequencing to identify exosome-related features of glioblastoma and constructing a prognostic model to identify BARD1 as a potential therapeutic target for GBM patients

Background

Glioblastoma (GBM) is a malignant primary brain tumor. This study focused on exploring the exosome-related features of glioblastoma to better understand its cellular composition and molecular characteristics.

Methods

Single-cell RNA sequencing (scRNA-seq) and spatial transcriptome RNA sequencing (stRNA-seq) were used to analyze the heterogeneity of glioblastomas. After data integration, cell clustering, and annotation, five algorithms were used to calculate scores for exosome-related genes(ERGs). Cell trajectory analysis and intercellular communication analysis were performed to explore exosome-related communication patterns. Spatial transcriptome sequencing data were analyzed to validate the findings. To further utilize exosome-related features to aid in clinical decision-making, a prognostic model was constructed using GBM's bulk RNA-seq.

Results

Different cell subpopulations were observed in GBM, with Monocytes/macrophages and malignant cells in tumor samples showing higher exosome-related scores. After identifying differentially expressed ERGs in malignant cells, pseudotime analysis revealed the cellular status of malignant cells during development. Intercellular communication analysis highlighted signaling pathways and ligand-receptor interactions. Spatial transcriptome sequencing confirmed the high expression of exosome-related gene features in the tumor core region. A prognostic model based on six ERGs was shown to be predictive of overall survival and immunotherapy outcome in GBM patients. Finally, based on the results of scRNA-seq and prognostic modeling as well as a series of cell function experiments, BARD1 was identified as a novel target for the treatment of GBM.

Conclusion

This study provides a comprehensive understanding of the exosome-related features of GBM in both scRNA-seq and stRNA-seq, with malignant cells with higher exosome-related scores exhibiting stronger communication with Monocytes/macrophages. In terms of spatial data, highly scored malignant cells were also concentrated in the tumor core region. In bulk RNA-seq, patients with a high exosome-related index exhibited an immunosuppressive microenvironment, which was accompanied by a worse prognosis as well as immunotherapy outcomes. Prognostic models constructed using ERGs are expected to be independent prognostic indicators for GBM patients, with potential implications for personalized treatment strategies for GBM. Knockdown of BARD1 in GBM cell lines reduces the invasive and value-added capacity of tumor cells, and thus BARD1-positively expressing malignant cells are a risk factor for GBM patients.

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.

Genomic analysis of immunogenic cell death-related subtypes for predicting prognosis and immunotherapy outcomes in glioblastoma multiforme

Immunogenic cell death (ICD), a unique form of cancer cell death, has therapeutic potential in anti-tumour immunotherapy. The aim of this study is to explore the predictive potential of ICD in the prognosis and immunotherapy outcomes of glioblastoma multiforme (GBM). RNA sequencing data and clinical information were downloaded from three databases. Unsupervised consistency clustering analysis was used to identify ICD-related clusters and gene clusters. Additionally, the ICD scores were determined using principal component analysis and the Boruta algorithm via dimensionality reduction techniques. Subsequently, three ICD-related clusters and three gene clusters with different prognoses were identified, with differences in specific tumour immune infiltration-related lymphocytes in these clusters. Moreover, the ICD score was well differentiated among patients with GBM, and the ICD score was considered an independent prognostic factor for patients with GBM. Furthermore, two datasets were used for the external validation of ICD scores as predictors of prognosis and immunotherapy outcomes. The validation analysis suggested that patients with high ICD scores had a worse prognosis. Additionally, a higher proportion of patients with high ICD scores were non-responsive to immunotherapy. Thus, the ICD score has the potential as a biomarker to predict the prognosis and immunotherapy outcomes of patients with GBM.

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.

Integrating multiple machine learning methods to construct glutamine metabolism-related signatures in lung adenocarcinoma

Background

Glutamine metabolism (GM) is known to play a critical role in cancer development, including in lung adenocarcinoma (LUAD), although the exact contribution of GM to LUAD remains incompletely understood. In this study, we aimed to discover new targets for the treatment of LUAD patients by using machine learning algorithms to establish prognostic models based on GM-related genes (GMRGs).

Methods

We used the AUCell and WGCNA algorithms, along with single-cell and bulk RNA-seq data, to identify the most prominent GMRGs associated with LUAD. Multiple machine learning algorithms were employed to develop risk models with optimal predictive performance. We validated our models using multiple external datasets and investigated disparities in the tumor microenvironment (TME), mutation landscape, enriched pathways, and response to immunotherapy across various risk groups. Additionally, we conducted in vitro and in vivo experiments to confirm the role of LGALS3 in LUAD.

Results

We identified 173 GMRGs strongly associated with GM activity and selected the Random Survival Forest (RSF) and Supervised Principal Components (SuperPC) methods to develop a prognostic model. Our model's performance was validated using multiple external datasets. Our analysis revealed that the low-risk group had higher immune cell infiltration and increased expression of immune checkpoints, indicating that this group may be more receptive to immunotherapy. Moreover, our experimental results confirmed that LGALS3 promoted the proliferation, invasion, and migration of LUAD cells.

Conclusion

Our study established a prognostic model based on GMRGs that can predict the effectiveness of immunotherapy and provide novel approaches for the treatment of LUAD. Our findings also suggest that LGALS3 may be a potential therapeutic target for LUAD.

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

Background

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

Methods

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

Results

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

Conclusion

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

Systematic review registration

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

T-cell exhaustion signatures characterize the immune landscape and predict HCC prognosis via integrating single-cell RNA-seq and bulk RNA-sequencing

Background

Hepatocellular carcinoma (HCC), the third most prevalent cause of cancer-related death, is a frequent primary liver cancer with a high rate of morbidity and mortality. T-cell depletion (TEX) is a progressive decline in T-cell function due to continuous stimulation of the TCR in the presence of sustained antigen exposure. Numerous studies have shown that TEX plays an essential role in the antitumor immune process and is significantly associated with patient prognosis. Hence, it is important to gain insight into the potential role of T cell depletion in the tumor microenvironment. The purpose of this study was to develop a trustworthy TEX-based signature using single-cell RNA-seq (scRNA-seq) and high-throughput RNA sequencing, opening up new avenues for evaluating the prognosis and immunotherapeutic response of HCC patients.

Methods

The International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) databases were used to download RNA-seq information for HCC patients. The 10x scRNA-seq. data of HCC were downloaded from GSE166635, and UMAP was used for clustering descending, and subgroup identification. TEX-related genes were identified by gene set variance analysis (GSVA) and weighted gene correlation network analysis (WGCNA). Afterward, we established a prognostic TEX signature using LASSO-Cox analysis. External validation was performed in the ICGC cohort. Immunotherapy response was assessed by the IMvigor210, GSE78220, GSE79671, and GSE91061cohorts. In addition, differences in mutational landscape and chemotherapy sensitivity between different risk groups were investigated. Finally, the differential expression of TEX genes was verified by qRT-PCR.

Result

11 TEX genes were thought to be highly predictive of the prognosis of HCC and substantially related to HCC prognosis. Patients in the low-risk group had a greater overall survival rate than those in the high-risk group, according to multivariate analysis, which also revealed that the model was an independent predictor of HCC. The predictive efficacy of columnar maps created from clinical features and risk scores was strong.

Conclusion

TEX signature and column line plots showed good predictive performance, providing a new perspective for assessing pre-immune efficacy, which will be useful for future precision immuno-oncology studies.

Pro- vs. Anti-Inflammatory Features of Monocyte Subsets in Glioma Patients

Monocytes constitute a heterogenous group of antigen-presenting cells that can be subdivided based on CD14, CD16 and SLAN expression. This division reflects the functional diversity of cells that may play different roles in a variety of pathologies including gliomas. In the current study, the three monocyte subpopulations: classical (CD14⁺ CD16⁺ SLAN^-), intermediate (CD14^dim CD16⁺ SLAN^-) and non-classical (CD14^low/- CD16⁺ SLAN⁺) in glioma patients' peripheral blood were analysed with flow cytometry. The immune checkpoint molecule (PD-1, PD-L1, SIRPalpha, TIM-3) expression along with pro- and anti-inflammatory cytokines (TNF, IL-12, TGF-beta, IL-10) were assessed. The significant overproduction of anti-inflammatory cytokines by intermediate monocytes was observed. Additionally, SLAN-positive cells overexpressed IL-12 and TNF when compared to the other two groups of monocytes. In conclusion, these results show the presence of different profiles of glioma patient monocytes depending on CD14, CD16 and SLAN expression. The bifold function of monocyte subpopulations might be an additional obstacle to the effectiveness of possible immunotherapies.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

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.

Meta-Analysis of RNA-Seq Datasets Identifies Novel Players in Glioblastoma

Glioblastoma is a devastating grade IV glioma with poor prognosis. Identification of predictive molecular biomarkers of disease progression would substantially contribute to better disease management. In the current study, we performed a meta-analysis of different RNA-seq datasets to identify differentially expressed protein-coding genes (PCGs) and long non-coding RNAs (lncRNAs). This meta-analysis aimed to improve power and reproducibility of the individual studies while identifying overlapping disease-relevant pathways. We supplemented the meta-analysis with small RNA-seq on glioblastoma tissue samples to provide an overall transcriptomic view of glioblastoma. Co-expression correlation of filtered differentially expressed PCGs and lncRNAs identified a functionally relevant sub-cluster containing DANCR and SNHG6, with two novel lncRNAs and two novel PCGs. Small RNA-seq of glioblastoma tissues identified five differentially expressed microRNAs of which three interacted with the functionally relevant sub-cluster. Pathway analysis of this sub-cluster identified several glioblastoma-linked pathways, which were also previously associated with the novel cell death pathway, ferroptosis. In conclusion, the current meta-analysis strengthens evidence of an overarching involvement of ferroptosis in glioblastoma pathogenesis and also suggests some candidates for further analyses.