oncoPredict: an R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data

Potential crosstalk between Naïve CD4+ T cells and SPP1+ Macrophages is associated with clinical outcome and therapeutic response in hepatocellular carcinoma

BACKGROUND

The highly heterogeneity of the tumor microenvironment (TME) in hepatocellular carcinoma (HCC) results in diverse clinical outcomes and therapeutic responses. This study aimed to investigate potential intercellular crosstalk and its impact on clinical outcomes and therapeutic responses.

METHODS

Single-cell RNA sequencing (scRNA-seq), spatial transcriptomics (ST) and bulk RNA sequencing (RNA-seq) datasets were integrated to comprehensively analyze the intercellular interactions within the TME. Multiplex immunohistochemistry was conducted to validate the intercellular interactions. A machine learning-based integrative procedure was used in bulk RNA-seq datasets to generate a risk model to predict prognosis and therapeutic responses.

RESULTS

Survival analyses based on the bulk RNA-seq datasets revealed the negative impact of the naïve Cluster of Differentiation 4+ (CD4) T cells and Secreted Phosphoprotein 1+ (SPP1) macrophages on prognosis. Furthermore, their intricate intercellular crosstalk and spatial colocalization were also observed by scRNA-seq and ST analyses. Based on this crosstalk, a machine learning model, termed the naïve CD4+ T cell and SPP1+ macrophage prognostic score (TMPS), was established in the bulk-RNA seq datasets for prognostic prediction. The TMPS achieved C-index values of 0.785, 0.715, 0.692 and 0.857, respectively, across 4 independent cohorts. A low TMPS was associated with a significantly increased survival rates, improved response to immunotherapy and reduced infiltration of immunosuppressive cells, such as. regulatory T cells. Finally, 8 potential sensitive drugs and 6 potential targets were predicted for patients based on their TMPS.

CONCLUSION

The crosstalk between naïve CD4+ T cells and SPP1+ macrophages play a crucial role in the TME. TMPS can reflect this crosstalk and serve as a valuable tool for prognostic stratification and guiding clinical decision-making.

The immunological landscape and silico analysis of key paraptosis regulator LPAR1 in gastric cancer patients

This study aims to identify key regulators of paraptosis in gastric cancer (GC) and explore their potential in guiding therapeutic strategies, especially in stomach adenocarcinoma (STAD). Genes associated with paraptosis were identified from the references and subjected to Cox regression analysis in the TCGA-STAD cohort. Using machine learning models, LPAR1 consistently ranked highest in feature importance. Multiple sequencing data showed that LPAR1 was significantly overexpressed in cancer-associated fibroblasts (CAFs). LPAR1 expression was significantly higher in normal tissues, and ROC analysis demonstrated its discriminative ability. Copy number alterations and microsatellite instability were significantly associated with LPAR1 expression. High LPAR1 expression correlated with advanced tumor grades and specific cancer immune subtypes, and multivariate analysis confirmed LPAR1 as an independent predictor of poor prognosis. LPAR1 expression was associated with different immune response metrics, including immune effector activation and upregulated chemokine secretion. High LPAR1 expression also correlated with increased sensitivity to compounds, such as BET bromodomain inhibitors I-BET151 and RITA, suggesting LPAR1 as a biomarker for predicting drug activity. FOXP2 showed a strong positive correlation with LPAR1 transcriptional regulation, while increased methylation of LPAR1 promoter regions was negatively correlated with gene expression. Knockdown of LPAR1 affected cell growth in most tumor cell lines, and in vitro experiments demonstrated that LPAR1 influenced extracellular matrix (ECM) contraction and cell viability in the paraptosis of CAFs. These findings suggest that LPAR1 is a critical regulator of paraptosis in GC and a potential biomarker for drug sensitivity and immunotherapy response. This underscores the role of CAFs in mediating tumorigenic effects and suggests that targeting LPAR1 could be a promising strategy for precision medicine in GC.

Genomic and transcriptomic profiling of inflammatory breast cancer reveals distinct molecular characteristics to non-inflammatory breast cancers

PURPOSE

Inflammatory breast cancer (IBC), a rare and highly aggressive form of breast cancer, accounts for 10% of breast cancer-related deaths. Previous omics studies of IBC have focused solely on one of genomics or transcriptomics and did not discover common differences that could distinguish IBC from non-IBC.

METHODS

Seventeen IBC patients and five non-IBC patients as well as additional thirty-three Asian breast cancer samples from TCGA-BRCA were included for the study. We performed whole-exon sequencing (WES) to investigate different somatic genomic alterations, copy number variants, and large structural variants between IBC and non-IBC. Bulk RNA sequencing (RNA-seq) was performed to examine the differentially expressed genes, pathway enrichment, and gene fusions. WES and RNA-seq data were further investigated in combination to discover genes that were dysregulated in both genomics and transcriptomics.

RESULTS

Copy number variation analysis identified 10 cytobands that showed higher frequency in IBC. Structural variation analysis showed more frequent deletions in IBC. Pathway enrichment and immune infiltration analysis indicated increased immune activation in IBC samples. Gene fusions including CTSC-RAB38 were found to be more common in IBC. We demonstrated more commonly dysregulated RAS pathway in IBC according to both WES and RNA-seq. Inhibitors targeting RAS signaling and its downstream pathways were predicted to possess promising effects in IBC treatment.

CONCLUSION

We discovered differences unique in Asian women that could potentially explain IBC etiology and presented RAS signaling pathway as a potential therapeutic target in IBC treatment.

Patient-derived rhabdomyosarcoma cells recapitulate the genetic and transcriptomic landscapes of primary tumors

Rhabdomyosarcoma (RMS) is the most common soft tissue sarcoma in childhood and adolescence. The availability of appropriate and well-characterized preclinical models for RMS is limited, posing a challenge for investigating the molecular mechanisms and evaluating new targeted compounds in preclinical settings. Here, we collected 51 RMS specimens (referred to as ZJUCH-RMS cohort) and established 9 patient-derived cells (PDCs) and validated the identity of these cells by the expression of RMS-specific markers. Whole-transcriptome analysis identified high-confidence mutations in ZJUCH-RMS cohort including RAS, TP53, ARID1A, MYOD1, and MYCN. Further studies showed that RMS PDCs retained the genetic alterations and the expression of RMS hallmark and dependency genes in matched primary tumors and acted as valuable tools to assess drug responses and pharmacogenomic interactions. Our study provides unique PDCs that are available for preclinical studies of RMS and further advances the feasibility of RMS PDCs as valuable tools for developing personalized treatments for patients.

Integrative analysis of multiple cell death model for precise prognosis and drug response prediction in gastric cancer

BACKGROUND

Gastric cancer (GC) is a common upper gastrointestinal tumor. However, the evaluation of prognosis and treatment response in patients with gastric cancer remains a challenge. Programmed cell death (PCD) is one of the important terminal paths for the cells of metazoans, and is involved in a variety of biological events that include morphogenesis, maintenance of tissue homeostasis, and elimination of harmful cells. The objective of this project is to investigate the predictive significance of cell death pathways and create prognostic signatures associated to cell death, with the purpose of forecasting prognosis and providing guidance for the treatment of gastric cancer.

METHODS

Gene transcription profiles and corresponding clinical data of gastric cancer patients were collected from The Cancer Genome Atlas (TCGA-STAD, n = 448) and the Gene Expression Comprehensive Database (GSE84437, n = 483). Thirteen types of cell death-related genes, including apoptosis, necroptosis, pyroptosis, ferroptosis, autophagy, cuprotosis, parthanatos, entotic cell death, netotic cell death, lysosome-dependent cell death, alkaliptosis, oxeiptosis, and disulfidptosis, were analysed. Cell death-related genes associated with prognosis were identified in the TCGA-STAD training cohort using Lasso-Cox regression to generate a risk score. Patients were categorized into high and low-risk groups based on the median risk score for survival difference analysis. Cell death-related genes associated with prognosis were identified in the TCGA-STAD training cohort using Lasso-Cox regression to generate a risk score. Additionally, the response to immunotherapy in the high-risk and low-risk groups was calculated using the oncoPredict algorithm. Futhermore, the model genes were validated in the GEO validation set.

RESULTS

A total of 324 differential programmed cell death (PCD)-related genes were identified, and 65 were selected through single-factor Cox analysis. Six PCD-related genes were ultimately identified by Lasso regression to construct a prognostic risk score model. The log-rank test revealed that patients in the high-risk group had inferior survival time compared with those in the low-risk group. The area under the ROC curve (AUC) for the training group at years 1, 3, and 5 were 0.684, 0.713, 0.743, respectively, while the AUC for the validation cohort at years 1, 3, and 5 were 0.695, 0.704, and 0.707, respectively. Unsupervised clustering identified potential subtypes included in the model, and a survival difference was also observed between the two subgroups. Multifactor Cox results, combined with clinical information, demonstrated that the prognostic risk score can serve as an independent prognostic factor, irrespective of other clinical features.

CONCLUSION

By comprehensively analyzing multiple cell death patterns, we have established a novel model that accurately forecasts the clinical prognosis and drug sensitivity of gastric cancer. It was found that all 12 representative drugs may not be suitable for patients in high-risk groups.

Immune-related glycosylation genes based classification predicts prognosis and therapy options of osteosarcoma

Osteosarcoma is the most common primary bone malignancy, with a very poor prognosis. Aberrant glycosylation is close involvement in osteosarcoma. Accordingly, this study aimed at investigating the role of glycosylation genes in the prognosis and therapy options of osteosarcoma. The microenvironment of osteosarcoma was assessed using estimate algorithm. A total of 20 immune-related glycosylation genes (IRGGs) was identified using Pearson correlation analysis. Accordingly, osteosarcoma patients were divided into C1 and C2 type using consensus clustering. Multiple algorithms (Xcell, MCP-counter, ssGSEA, epic, quantiseq), cancer immune cycle analysis, and GSVA were applied to estimate the immune, molecule and metabolism characteristics of osteosarcoma, indicating that C1 type was featured with high immune infiltration, high glycosylation, enriched MEK signaling, and good prognosis, while C2 type was characterized by more metastasis, enriched immunotherapy-positive gene signatures, high tumor mutation burden, and poor prognosis. Results from TIDE algorithm and immunotherapy datasets suggested the C2 type's preference of immune checkpoint inhibitors (ICIs), while data of GDSC, CMap analysis and cell experiments indicated that C1 type was sensitivity to MEK inhibitor PD0325901. In addition, univariate Cox and Lasso analysis was combined to establish an IRGGs' risk score containing 6 genes (B3GNT8, FUT7, GAL3ST4, GALNT14, HS3ST2, and MFNG). The data of DCA and ROC indicated its well prediction of prognosis in osteosarcoma. Finally, cellular location analysis showed that the 6 genes not only distributed in tumor cells but also in immune cells. In summary, the classification and risk score based on IRGGs effectively predicted the prognosis and therapy options of osteosarcoma. Further studies on IRGGs may contribute to the understanding of cancer immunity in osteosarcoma.

A Novel Mitochondrial-Related Gene Signature for the Prediction of Prognosis and Therapeutic Efficacy in Lower-Grade Glioma

Lower-grade glioma (LGG) is a common primary brain tumor with a highly heterogeneous clinical presentation, and its prognosis cannot be accurately predicted by current histopathology. It has been found that mitochondria play an important role in hypoxia, angiogenesis, and energy metabolism in glioma, and mitochondrial function may have an important impact on LGG prognosis. The goal of this study was to develop a novel prognostic model based on Mitochondrial-related genes (MRGs). We first analyzed the somatic alterations profiles of MRGs in patients with LGG and found that somatic alterations were common in LGG and correlated with prognosis. Using RNA-seq data from TCGA and CGGA, 12 prognosis-related MRGs were identified to construct a mitochondrial activation score (MiAS) model by combining univariate regression and LASSO regression analysis. The model and nomogram were evaluated using the area under the ROC curve with AUC = 0.910. The model was closely correlated with the clinical characteristics of LGG patients and performed well in predicting the prognosis of LGG patients with significantly shorter overall survival (OS) time in the high-MiAS group. GSVA and GSEA results showed that oxidative stress, pro-cancer, and immune-related pathways were significantly enriched in the high-MiAS group. CIBERSORT results showed that MiAS was significantly associated with immune cell infiltration in LGG. Macrophage M1 and follicular helper T cells had increased infiltration in the high-MiAS group. TIDE predicted a better immunotherapy outcome in patients in the low-MiAS group. Finally, using data from the CTRPv2 and GDSC2 datasets to assess chemotherapy response in LGG, it was predicted that the chemotherapeutic agents AZD6482, MG-132, and PLX-4720 might be potential agents for patients in the high-MiAS group of LGG. In addition, we performed in vitro experiments and found that knockdown of OCIAD2 expression reduced the abilities of glioma cells to proliferate, migrate, and invade. In contrast, overexpression of OCIAD2 enhanced these abilities of glioma cells. This study found that MRGs were correlated with LGG patient prognosis, which is expected to provide new treatment strategies for LGG patients with different MiAS.

Pan-cancer characterization of cellular senescence reveals its inter-tumor heterogeneity associated with the tumor microenvironment and prognosis

Cellular senescence (CS) is characterized by the irreversible cell cycle arrest and plays a key role in aging and diseases, such as cancer. Recent years have witnessed the burgeoning exploration of the intricate relationship between CS and cancer, with CS recognized as either a suppressing or promoting factor and officially acknowledged as one of the 14 cancer hallmarks. However, a comprehensive characterization remains absent from elucidating the divergences of this relationship across different cancer types and its involvement in the multi-facets of tumor development. Here we systematically assessed the cellular senescence of over 10,000 tumor samples from 33 cancer types, starting by defining a set of cancer-associated CS signatures and deriving a quantitative metric representing the CS status, called CS score. We then investigated the CS heterogeneity and its intricate relationship with the prognosis, immune infiltration, and therapeutic responses across different cancers. As a result, cellular senescence demonstrated two distinct prognostic groups: the protective group with eleven cancers, such as LIHC, and the risky group with four cancers, including STAD. Subsequent in-depth investigations between these two groups unveiled the potential molecular and cellular mechanisms underlying the distinct effects of cellular senescence, involving the divergent activation of specific pathways and variances in immune cell infiltrations. These results were further supported by the disparate associations of CS status with the responses to immuno- and chemo-therapies observed between the two groups. Overall, our study offers a deeper understanding of inter-tumor heterogeneity of cellular senescence associated with the tumor microenvironment and cancer prognosis.

IGFBP7 regulates cell proliferation and migration through JAK/STAT pathway in gastric cancer and is regulated by DNA and RNA methylation

New biomarkers for early diagnosis of gastric cancer (GC), the second leading cause of cancer-related death, are urgently needed. IGFBP7, known to play various roles in multiple tumours, is complexly regulated across diverse cancer types, as evidenced by our pancancer analysis. Bioinformatics analysis revealed that IGFBP7 expression was related to patient prognosis, tumour clinicopathological characteristics, tumour stemness, microsatellite instability and immune cell infiltration, as well as the expression of oncogenes and immune checkpoints. GSEA links IGFBP7 to several cancer-related pathways. IGFBP7 deficiency inhibited GC cell proliferation and migration in vitro. Furthermore, an in vivo nude mouse model revealed that IGFBP7 downregulation suppressed the tumorigenesis of GC cells. Western blotting analysis showed that the JAK1/2-specific inhibitor ruxolitinib could rescue alterations induced by IGFBP7 overexpression in GC cells. Additionally, our bioinformatics analysis and in vitro assays suggested that IGFBP7 is regulated by DNA methylation at the genetic level and that the RNA m6A demethylase FTO modulates it at the posttranscriptional level. This study emphasizes the clinical relevance of IGFBP7 in GC and its influence on cell proliferation and migration via the JAK/STAT signalling pathway. This study also highlights the regulation of IGFBP7 in GC by DNA and m6A RNA methylation.

Comprehensive evaluation of immunological attributes and immunotherapy responses of positive T cell function regulators in colorectal cancer

BACKGROUND

Positive regulators of T-cell function (PTFRs), integral to T-cell proliferation and activation, have been identified as potential prognostic markers in colorectal cancer (CRC). Despite this, their role within the tumor microenvironment (TME) and their response to immunotherapy are not yet fully understood.

METHODS

This study delved into PTFR-related CRC subtypes by analyzing four independent transcriptome datasets, emphasizing the most significant prognostic PTFRs. We identified differentially expressed genes (DEGs) between two subtypes and developed a PTFR risk model using LASSO and Cox regression methods. The model's associations with survival time, clinical features, TME characteristics, tumor mutation profiles, microsatellite instability (MSI), cancer stem cell (CSC) index, and responses to chemotherapy, targeted therapy, and immunotherapy were subsequently explored.

RESULTS

The PTFR risk model demonstrated a strong predictive capacity for CRC. It facilitated the estimation of immune cell composition, HLA expression levels, immune checkpoint expression, mutation burden, CSC index features, and the effectiveness of immunotherapy.

CONCLUSIONS

This study enhances our understanding of the role of PTFRs in CRC progression and introduces an innovative assessment framework for CRC immunotherapy. This framework improves the prediction of treatment outcomes and aids in the customization of therapeutic strategies.

Single-cell RNA-seq reveals heterogeneity in metastatic renal cell carcinoma and effect of anti-angiogenesis therapy in the pancreas metastatic lesion

Metastatic clear cell renal cell carcinoma has heterogenous tumor microenvironment (TME). Among the metastatic lesions, pancreas metastasis is rare and controversy in treatment approaches. Here, extensive primary and metastatic lesion samples were included by single-cell RNA-seq to decipher the distinct metastasis TME. The hypoxic and inflammatory TME of pancreas metastasis was decoded in this study, and the activation of PAX8-myc signaling, and metabolic reprogramming were observed. The active components including endothelial cells, fibroblasts and T cells were profiled. Meanwhile, we also evaluated the effect of anti-angiogenesis treatment in the pancreas metastasis patient. The potential mechanisms of pancreatic tropism, instability of genome, and the response of immunotherapy were also discussed in this work. Taken together, our findings suggest a clue to the heterogeneity in metastasis TME and provide evidence for the treatment of pancreas metastasis in renal cell carcinoma patients.

Advanced machine learning unveils CD8 + T cell genetic markers enhancing prognosis and immunotherapy efficacy in breast cancer

BACKGROUND

Breast cancer (BC) is the most common cancer in women and poses a significant health burden, especially in China. Despite advances in diagnosis and treatment, patient variability and limited early detection contribute to poor outcomes. This study examines the role of CD8 + T cells in the tumor microenvironment to identify new biomarkers that improve prognosis and guide treatment strategies.

METHODS

CD8 + T-cell marker genes were identified using single-cell RNA sequencing (scRNA-seq), and a CD8 + T cell-related gene prognostic signature (CTRGPS) was developed using 10 machine-learning algorithms. The model was validated across seven independent public datasets from the GEO database. Clinical features and previously published signatures were also analyzed for comparison. The clinical applications of CTRGPS in biological function, immune microenvironment, and drug selection were explored, and the role of hub genes in BC progression was further investigated.

RESULTS

We identified 71 CD8 + T cell-related genes and developed the CTRGPS, which demonstrated significant prognostic value, with higher risk scores linked to poorer overall survival (OS). The model's accuracy and robustness were confirmed through Kaplan-Meier and ROC curve analyses across multiple datasets. CTRGPS outperformed existing prognostic signatures and served as an independent prognostic factor. The role of the hub gene TTK in promoting malignant proliferation and migration of BC cells was validated.

CONCLUSION

The CTRGPS enhances early diagnosis and treatment precision in BC, improving clinical outcomes. TTK, a key gene in the signature, shows promise as a therapeutic target, supporting the CTRGPS's potential clinical utility.

Matrix stiffness-related extracellular matrix signatures and the DYNLL1 protein promote hepatocellular carcinoma progression through the Wnt/β-catenin pathway

BACKGROUND

In hepatocellular carcinoma (HCC) treatment, first-line targeted therapy in combination with immune checkpoint inhibitors (ICIs) has improved patient prognosis, but the 5-year survival rate is far from satisfactory. Studies have shown that the extracellular matrix (ECM) is an essential part of the tumour microenvironment (TME) and participates in the progression of malignant tumours. ECM remodelling can enhance matrix stiffness in cirrhosis patients, induce an immunosuppressive microenvironment network, and affect the efficacy of targeted therapies and ICIs for treating HCC. However, the exact mechanism is still unclear.

METHODS

We downloaded data from public databases, selected differentially expressed ECM proteins associated with matrix stiffness, constructed and validated a prognostic model of HCC using Lasso Cox regression, and investigated the roles and mechanism of one of the ECM proteins, dynein light chain LC8-type 1 (DYNLL1), in HCC proliferation, migration, and apoptosis via in vitro experiments.

RESULTS

In this study, the risk score of the matrix stiffness-related ECM protein model effectively predicted the prognosis of HCC patients. The high- and low-risk subgroups of the model also showed differences in immune cells, immune functions, and drug sensitivity. DYNLL1 promoted HCC cell progression and migration and inhibited HCC cell apoptosis through the Wnt/β-catenin pathway in vitro.

CONCLUSION

The expression of matrix stiffness-related ECM proteins could be an independent predictor of HCC prognosis. DYNLL1, an oncogenic gene in HCC, has the potential to be a new target for HCC treatment.

Copper metabolism-related signature for prognosis prediction and MMP13 served as malignant factor for breast cancer

Objectives

To comprehensively analyze the copper metabolism in Breast cancer, we established a prognostic signature for breast cancer (BC) related to copper metabolism.

Methods

Copper metabolism-related genes were sourced from previous literatures and were selected by the Univariate Cox regression. Cu-enrichment scores were calculated via ssGSEA. Differentially expressed genes were identified with limma between high and low Cu-enrichment scores group, then we used the Random Survival Forest and LASSO to build the CuScore for BC. Kaplan-Meier analysis, ROC curves, and Cox regression were used to evaluate CuScore. Genomic mutations were analyzed with GISTIC. Immune cells were examined using ESTIMATE, ssGSEA and TIMER. Enrichment analysis used clusterProfiler and GSVA. The GDSC database and oncoPredict package analyzed chemotherapeutic sensitivity. MMP13 was selected for in vitro assays.

Results

Four copper metabolism-related genes (UBE2D2, SLC31A1, ATP7A, and MAPK1) with prognostic value were identified. Higher expression levels of these genes were associated with higher Cu-enrichment scores, a factor of malignancy in breast cancer. Among 115 differentially expressed genes, 19 prognostic genes were identified, with three (CEACAM5, MMP13, and CRISP3) highlighted by Random Survival Forest and LASSO. Higher CuScores correlated with worse prognoses and were effective in predicting breast cancer outcomes. CuScore and metastasis were independent prognostic factors. Tumor-infiltrating immune cells were associated with lower CuScores. GO-GSEA analysis indicated six immune-related pathways might be regulated by CuScore. Patients with higher CuScores had lower TMB and were more sensitive to Sapitinib and LCL161, while those with lower CuScores might respond better to anti-PD1 therapy. High MMP13 expression in breast cancer was linked to malignancy, affecting cell proliferation and migration.

Conclusion

The identified copper metabolism-related gene signature has the potential to predict prognosis and guide clinical treatment for BC. Among these genes, MMP13 may act as a malignant factor in BC.

Identification of target genes co-regulated by four key histone modifications of five key regions in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a cancer with high morbidity and mortality. Studies have shown that histone modification plays an important regulatory role in the occurrence and development of HCC. However, the specific regulatory effects of histone modifications on gene expression in HCC are still unclear. This study focuses on HepG2 cell lines and hepatocyte cell lines. First, the distribution of histone modification signals in the two cell lines was calculated and analyzed. Then, using the random forest algorithm, we analyzed the effects of different histone modifications and their modified regions on gene expression in the two cell lines, four key histone modifications (H3K36me3, H3K4me3, H3K79me2, and H3K9ac) and five key regions that co-regulate gene expression were obtained. Subsequently, target genes regulated by key histone modifications in key regions were screened. Combined with clinical data, Cox regression analysis and Kaplan-Meier survival analysis were performed on the target genes, and four key target genes (CBX2, CEBPZOS, LDHA, and UMPS) related to prognosis were identified. Finally, through immune infiltration analysis and drug sensitivity analysis of key target genes, the potential role of key target genes in HCC was confirmed. Our results provide a theoretical basis for exploring the occurrence of HCC and propose potential biomarkers associated with histone modifications, which may be potential drug targets for the clinical treatment of HCC.

ENO2 in progression and treatment of colon adenocarcinoma: integrative bioinformatics analysis on non-apoptotic cell death

Colon adenocarcinoma (COAD) is one of the most common types of cancer. The interconnection between non-apoptotic cell death and COAD has not been adequately addressed. In our study, an integrative bioinformatics analysis was performed to explore non-apoptotic cell death-related biomarkers in COAD. ENO2 was determined as a potent biomarker for prognosis, drug response, immunity, and immunotherapy prediction. We used EdU and RT-qPCR assays to test our hypothesis and investigate how the ENO2 gene may influence or regulate cancer-related processes. ENO2 was expected to be a potential target in COAD.

Identification of key ferroptosis genes and subtypes in kidney renal clear cell carcinoma

Tumour immunity is highly important for the occurrence and development of tumours, and many cancers are resistant to ferroptosis. This study aims to explore the relationship between ferroptosis-related genes (FRGs) and the immunological characteristics of kidney renal clear cell carcinoma (KIRC). We obtained RNA-seq profiles and clinical data of KIRC patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases and identified CD44 and GLRX5 as the key FRGs involved in KIRC immune infiltration through Spearman's correlation analysis. Based on the expression of CD44 and GLRX5, the consensus clustering algorithm was used to classify the TCGA-KIRC samples into two clusters. A nomogram was constructed to evaluate the prognosis of KIRC patients. ESTIMATE, CIBERSORT, and single-sample gene set enrichment analysis (ssGSEA) were performed to evaluate immune infiltration between the two clusters. A weighted gene co-expression network analysis (WGCNA) was used to identify the most relevant genes to the clusters and immunity. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed. The external dataset GSE53757 was used to validate the immunological features between the two clusters. Cluster 2 patients had more active immune infiltration and might be more sensitive to immunotherapy; Cluster 2 patients also had a worse prognosis and might be at a more advanced stage of KIRC. We identified key ferroptosis-related genes and subgroups involved in the immune infiltration of KIRC, which is highly important for exploring the molecular mechanisms and treatments of KIRC.

Unraveling the immune landscape of lung adenocarcinoma: insights for tailoring therapeutic approaches

Lung adenocarcinoma (LUAD), a prevalent type of non-small cell lung cancer (NSCLC), was known for its diversity and intricate tumor microenvironment (TME). Comprehending the interaction among human immune-related genes (IRGs) and the TME is vital in the creation of accurate predictive models and specific treatments. We created a risk score based on IRGs and designed a nomogram to predict the prognosis of LUAD accurately. This involved a thorough examination of TME and the infiltration of immune cells in both high-risk and low-risk LUAD groups. Furthermore, the examination of the association between characteristic genes (BIRC5 and BMP5) and immune cells, along with immune checkpoints in the TME, was also conducted. The findings of our research unveiled unique immune profiles and interactions among individuals in the high- and low-risk categories, which contribute to variations in prognosis. LUAD demonstrated significant associations between BIRC5, BMP5, immune cells, and checkpoints, suggesting their involvement in disease advancement and resistance to medication. Furthermore, by correlating our findings with a multidrug database, we identified specific LUAD patient subsets that might benefit from tailored treatments. Our study establishes a groundbreaking prognostic model for LUAD, which not only underscores the importance of the immune context in LUAD but also paves the way for advancing precision medicine strategies in this complex malignancy.

Integrating machine learning and multi-omics analysis to develop an asparagine metabolism immunity index for improving clinical outcome and drug sensitivity in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a malignancy affecting the respiratory system. Most patients are diagnosed with advanced or metastatic lung cancer due to the fact that most of their clinical symptoms are insidious, resulting in a bleak prognosis. Given that abnormal reprogramming of asparagine metabolism (AM) has emerged as an emerging therapeutic target for anti-tumor therapy. However, the clinical significance of abnormal reprogramming of AM in LUAD patients is unclear. In this study, we collected 864 asparagine metabolism-related genes (AMGs) and used a machine-learning computational framework to develop an asparagine metabolism immunity index (AMII) for LUAD patients. Through the utilization of median AMII scores, LUAD patients were segregated into either a low-AMII group or a high-AMII group. We observed outstanding performance of AMII in predicting survival prognosis in LUAD patients in the TCGA-LUAD cohort and in three externally independently validated GEO cohorts (GSE72094, GSE37745, and GSE30219), and poorer prognosis for LUAD patients in the high-AMII group. The results of univariate and multivariate analyses showed that AMII can be used as an independent risk factor for LUAD patients. In addition, the results of C-index analysis and decision analysis showed that AMII-based nomograms had a robust performance in terms of accuracy of prognostic prediction and net clinical benefit in patients with LUAD. Excitingly, LUAD patients in the low-AMII group were more sensitive to commonly used chemotherapeutic drugs. Consequently, AMII is expected to be a novel diagnostic tool for clinical classification, providing valuable insights for clinical decision-making and personalized management of LUAD patients.

Novel prognostic signature for hepatocellular carcinoma using a comprehensive machine learning framework to predict prognosis and guide treatment

Background

Hepatocellular carcinoma (HCC) is highly aggressive, with delayed diagnosis, poor prognosis, and a lack of comprehensive and accurate prognostic models to assist clinicians. This study aimed to construct an HCC prognosis-related gene signature (HPRGS) and explore its clinical application value.

Methods

TCGA-LIHC cohort was used for training, and the LIRI-JP cohort and HCC cDNA microarray were used for validation. Machine learning algorithms constructed a prognostic gene label for HCC. Kaplan-Meier (K-M), ROC curve, multiple analyses, algorithms, and online databases were used to analyze differences between high- and low-risk populations. A nomogram was constructed to facilitate clinical application.

Results

We identified 119 differential genes based on transcriptome sequencing data from five independent HCC cohorts, and 53 of these genes were associated with overall survival (OS). Using 101 machine learning algorithms, the 10 most prognostic genes were selected. We constructed an HCC HPRGS with four genes (SOCS2, LCAT, ECT2, and TMEM106C). Good predictive performance of the HPRGS was confirmed by ROC, C-index, and K-M curves. Mutation analysis showed significant differences between the low- and high-risk patients. The low-risk group had a higher response to transcatheter arterial chemoembolization (TACE) and immunotherapy. Treatment response of high- and low-risk groups to small-molecule drugs was predicted. Linifanib was a potential drug for high-risk populations. Multivariate analysis confirmed that HPRGS were independent prognostic factors in TCGA-LIHC. A nomogram provided a clinical practice reference.

Conclusion

We constructed an HPRGS for HCC, which can accurately predict OS and guide the treatment decisions for patients with HCC.

Machine learning-based discovery of UPP1 as a key oncogene in tumorigenesis and immune escape in gliomas

Introduction

Gliomas are the most common and aggressive type of primary brain tumor, with a poor prognosis despite current treatment approaches. Understanding the molecular mechanisms underlying glioma development and progression is critical for improving therapies and patient outcomes.

Methods

The current study comprehensively analyzed large-scale single-cell RNA sequencing and bulk RNA sequencing of glioma samples. By utilizing a series of advanced computational methods, this integrative approach identified the gene UPP1 (Uridine Phosphorylase 1) as a novel driver of glioma tumorigenesis and immune evasion.

Results

High levels of UPP1 were linked to poor survival rates in patients. Functional experiments demonstrated that UPP1 promotes tumor cell proliferation and invasion and suppresses anti-tumor immune responses. Moreover, UPP1 was found to be an effective predictor of mutation patterns, drug response, immunotherapy effectiveness, and immune characteristics.

Conclusions

These findings highlight the power of combining diverse machine learning methods to identify valuable clinical markers involved in glioma pathogenesis. Identifying UPP1 as a tumor growth and immune escape driver may be a promising therapeutic target for this devastating disease.

Leveraging programmed cell death patterns to predict prognosis and therapeutic sensitivity in OSCC

OBJECTIVES

Intricate associations between programmed cell death (PCD) and cancer development and treatment outcomes have been increasingly appreciated. Here, we integrated 12 PCD patterns to construct a novel biomarker, cell death index (CDI), for oral squamous cell carcinoma (OSCC) prognostication and therapeutic prediction.

MATERIALS AND METHODS

Univariate Cox regression, Kaplan-Meier survival, and LASSO analyses were performed to construct the CDI. A nomogram combining CDI and selected clinicopathological parameters was established by multivariate Cox regression. The associations between CDI and immune landscape and therapeutic sensitivity were estimated. Single-cell RNA-seq data of OSCC was used to infer CDI genes in selected cell types and determine their expression along cell differentiation trajectory.

RESULTS

Ten selected PCD genes derived a novel prognostic signature for OSCC. The predictive prognostic performance of CDI and nomogram was robust and superior across multiple independent patient cohorts. CDI was negatively associated with tumor-infiltrating immune cell abundance and immunotherapeutic outcomes. Moreover, scRNA-seq data reanalysis revealed that GSDMB, IL-1A, PRKAA2, and SFRP1 from this signature were primarily expressed in cancer cells and involved in cell differentiation.

CONCLUSIONS

Our findings established CDI as a novel powerful predictor for prognosis and therapeutic response for OSCC and suggested its potential involvement in cancer cell differentiation.

RNA sequencing identifies lung cancer lineage and facilitates drug repositioning

Recent breakthrough therapies have improved survival rates in non-small cell lung cancer (NSCLC), but a paradigm for prospective confirmation is still lacking. Patientdatasets were mainly downloaded from TCGA, CPTAC and GEO. We conducted downstream analysis by collecting metagenes and generated 42-gene subtype classifiers to elucidate biological pathways. Subsequently, scRNA, eRNA, methylation, mutation, and copy number variation were depicted from a phenotype perspective. Enhancing the clinical translatability of molecular subtypes, preclinical models including CMAP, CCLE, and GDSC were utilized for drug repositioning. Importantly, we verified the presence of previously described three phenotypes including bronchioid, neuroendocrine, and squamoid. Poor prognosis was seen in squamoid and neuroendocrine clusters for treatment-naive and immunotherapy populations. The neuroendocrine cluster was dominated by STK11 mutations and 14q13.3 amplifications, whose related methylated loci are predictive of immunotherapy. And the greatest therapeutic potential lies in the bronchioid cluster. We further estimated the relative cell abundance of the tumor microenvironment (TME), specific cell types could be reflected among three clusters. Meanwhile, the higher portion of immune cell infiltration belonged to bronchioid and squamoid, not the neuroendocrine cluster. In drug repositioning, MEK inhibitors resisted bronchioid but were squamoid-sensitive. To conceptually validate compounds/targets, we employed RNA-seq and CCK-8/western blot assays. Our results indicated that dinaciclib and alvocidib exhibited similar activity and sensitivity in the neuroendocrine cluster. Also, a lineage factor named KLF5 recognized by inferred transcriptional factors activity could be suppressed by verteporfin.

Identification and validation of diagnostic and prognostic biomarkers in prostate cancer based on WGCNA

BACKGROUND

Prostate cancer (PCa) represents a significant health challenge for men, and the advancement of the disease often results in a grave prognosis for patients. Therefore, the identification of biomarkers associated with the diagnosis and prognosis of PCa holds paramount importance in patient health management.

METHODS

The datasets pertaining to PCa were retrieved from the Gene Expression Omnibus (GEO) database. Weighted gene co-expression network analysis (WGCNA) was conducted to investigate the modules specifically associated with the diagnosis of PCa. The hub genes were identified using the LASSO regression analysis. The expression levels of these hub genes were further validated by qRT-PCR experiments. Receiver operating characteristic (ROC) curves and nomograms were employed as evaluative measures for assessing the diagnostic value.

RESULTS

The blue module identified by WGCNA exhibited a strong association with PCa. Six hub genes (SLC14A1, COL4A6, MYOF, FLRT3, KRT15, and LAMB3) were identified by LASSO regression analysis. Further verification confirmed that these six genes were significantly downregulated in tumor tissues and cells. The six hub genes and the nomogram demonstrated substantial diagnostic value, with area under the curve (AUC) values ranging from 0.754 to 0.961. Moreover, patients with low expression levels of these six genes exhibited elevated T/N pathological stage and Gleason score, implying a more advanced disease state. Meanwhile, their progression-free survival (PFS) was observed to be potentially poorer. Finally, a significant association could be observed between the expression of these genes and the dysregulation of immune cells, along with drug sensitivity.

CONCLUSIONS

In summary, our study identified six hub genes, namely SLC14A1, COL4A6, MYOF, FLRT3, KRT15, and LAMB3, which can be utilized to establish a diagnostic model for PCa. The discovery may offer potential molecular targets for clinical diagnosis and treatment of PCa.

Protocol to develop a chemotherapy drug screening process by constructing a cancer prognostic model using public databases

Drug resistance is currently the biggest challenge in cancer chemotherapy. Here, we present a protocol to develop a chemotherapy drug screening process by constructing a cancer prognostic model (PM) using public databases. We describe steps for downloading code and data, preparing the expression matrix and metadata for analysis, screening modeling genes, and constructing a PM. We then detail procedures for constructing predictive websites for cancer patients' survival based on their age, tumor stage, gene expression levels, and risk scores. For complete details on the use and execution of this protocol, please refer to Bai et al.1.

Development and validation of a mitotic catastrophe-related genes prognostic model for breast cancer

Background

Breast cancer has become the most common malignant tumor in women worldwide. Mitotic catastrophe (MC) is a way of cell death that plays an important role in the development of tumors. However, the exact relationship between MC-related genes (MCRGs) and the development of breast cancer is still unclear, and further research is needed to elucidate this complexity.

Methods

Transcriptome data and clinical data of breast cancer were downloaded from the Cancer Genome Atlas (TCGA) database and the Gene Expression Omnibus (GEO) database. We identified differential expression of MCRGs by comparing tumor tissue with normal tissue. Subsequently, we used COX regression analysis and LASSO regression analysis to construct the prognosis risk model of MCRGs. Kaplan-Meier survival curve and receiver operating characteristic (ROC) curve were used to evaluate the predictive ability of prognostic model. Moreover, the clinical relevance, gene set enrichment analysis (GSEA), immune landscape, tumor mutation burden (TMB), and immunotherapy and drug sensitivity analysis between high-risk and low-risk groups were systematically investigated. Finally, we validated the expression levels of genes involved in constructing the prognostic model through real-time quantitative polymerase chain reaction (RT-qPCR) at the cellular and tissue levels.

Results

We identified 12 prognostic associated MCRGs, four of which were selected to construct prognostic model. The Kaplan-Meier analysis suggested that patients in the high-risk group had a shorter overall survival (OS). The Cox regression analysis and ROC analysis indicated that risk model had independent and excellent ability in predicting prognosis of breast cancer patients. Mechanistically, a remarkable difference was observed in clinical relevance, GSEA, immune landscape, TMB, immunotherapy response, and drug sensitivity analysis. RT-qPCR results showed that genes involved in constructing the prognostic model showed significant abnormal expressions and the expression change trends were consistent with the bioinformatics results.

Conclusions

We established a prognosis risk model based on four MCRGs that had the ability to predict clinical prognosis and immune landscape, proposing potential therapeutic targets for breast cancer.

Identification of prognostic biomarkers in neuroblastoma using WGCNA and multi-omics analysis

BACKGROUND

Neuroblastoma (NB) is one of the most frequent parenchymal tumors among children, with a high degree of heterogeneity and wide variation in clinical presentation. Despite significant therapeutic advances in recent years, long-term survival in high-risk patients remains low, emphasizing the urgent need to find new biomarkers and construct reliable prognostic models.

METHODS

In this study, data from neuroblastoma samples in the ArrayExpress database were utilized to identify key gene modules and pivotal genes associated with NB prognosis by weighted gene co-expression network analysis (WGCNA). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis was performed using the DAVID database. Based on these hub genes, survival prognosis models were constructed and validated on an independent validation set in the Gene Expression Omnibus (GEO) database. Differences in biological functions and immune microenvironments and the sensitivity to pharmacological and immunotherapeutic treatments of patients in the high- and low-risk groups were examined by gene set enrichment analysis (GSEA) and immune infiltration analysis.

RESULTS

WGCNA identified 14 gene modules and screened the module with the highest relevance to the International Neuroblastoma Staging System (INSS), containing 60 pivotal genes. GO and KEGG analyses demonstrated that these pivotal genes were mainly implicated in biological processes and signaling pathways including DNA replication, cell division, mitotic cell cycle, and cell cycle. Based on Lasso regression and COX regression analysis, a prognostic model containing DHFR, GMPS and E2F3 was constructed, and the RiskScore was significantly correlated with the 1-, 3- and 5-year survival of the patients. GSEA and immune infiltration analyses revealed significant differences in the levels of cell cycle-related pathways and immune cell infiltration between the high and low RiskScore groups. In particular, patients in the high-risk group are less likely to benefit from immunotherapy and may be better suited for treatment with drugs such as Oxaliplatin and Alpelisib.

CONCLUSION

This research systematically identified biomarkers related to NB prognosis and developed a reliable prognostic model applying WGCNA and multiple bioinformatics methods. The model has important application value in predicting patients' prognosis, evaluating drug sensitivity and immunotherapy effect, and provides new ideas and directions for precise treatment of neuroblastoma.

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.

Characterizing m6A modification factors and their interactions in colorectal cancer: implications for tumor subtypes and clinical outcomes

BACKGROUND

The study aims to comprehensively combine colorectal cancer data cohorts in order to analyze the effects of various DNA methylation-coding genes on colorectal cancer patients. The annual incidence and mortality of colorectal cancer are very high, and there are no effective treatments for advanced colorectal cancer. DNA methylation is a method widely used to regulate epigenetics in the molecular mechanism study of tumors.

METHOD

Three single-cell cohorts GSE166555, GSE146771, and EMTAB8107, and five transcriptome cohorts GSE17536, GSE39582, GSE72970, and TCGA-CRC (TCGA-COAD and TCGA-READ) were applied in this study. 2 erasers (ALKBH5 and FTO), There are 7 writers (METTL3, METTL14, WTAP, VIRMA, RBM15, RBM15B, and ZC3H13) and 11 readers (YTHDC1, IGF2BP1, IGF2BP2, IGF2BP3, YTHDF1, YTHDF3, YTHDC2, and HNRNPA2B1, YTHDF2, HNRNPC and RBMX), a total of 20 M6A regulators, were used as the basis of the dataset in this study and were applied to the construction of molecular typing and prognostic models. Drugs that are differentially sensitive in methylation-regulated gene-related prognostic models were identified using the ConsensusClusterPlus package, which was also used to identify distinct methylation regulatory expression patterns in colorectal cancer and to model the relationship between tissue gene expression profiles and drug IC50 values. Finally, TISCH2 assessed which immune cells were significantly expressed with M6A scores. The immunosuppression of M6A methylation is spatially explained.

RESULTS

This study used data from 583 CRC patients in the TCGA-CRC cohort. Firstly, the mutation frequency and CNV variation frequency of 20 m6A modification-related factors were analyzed, and the corresponding histogram and heat map were drawn. The study next analyzed the expression variations between mutant and wild forms of the VIRMA gene and explored differences in the expression of these variables in tumor and normal tissues. In addition, the samples were divided into different subgroups by molecular clustering method based on m6A modification, and each subgroup's expression and clinicopathological characteristics were analyzed. Finally, we compared prognostic differences, tumor microenvironment (TME) characteristics, immune cell infiltration, and gene function enrichment among different subpopulations. We also developed a colorectal cancer m6A-associated gene signature and validated its prognostic effects across multiple cohorts. Finally, using single-cell RNA sequencing data, we confirmed that tumor cells show elevated expression of m6A-related gene signatures.

DISCUSSION

This study explored the mutation frequency, expression differences, interactions, molecular clustering, prognostic effect, and association with tumor characteristics of m6A modification-related factors in CRC and validated them at the single-cell level. These results clarify the association between m6A alteration and colorectal cancer (CRC) and offer important insights into the molecular recognition and management of cancer.

Machine learning-based identification of an immunotherapy-related signature to enhance outcomes and immunotherapy responses in melanoma

Background

Immunotherapy has revolutionized skin cutaneous melanoma treatment, but response variability due to tumor heterogeneity necessitates robust biomarkers for predicting immunotherapy response.

Methods

We used weighted gene co-expression network analysis (WGCNA), consensus clustering, and 10 machine learning algorithms to develop the immunotherapy-related gene model (ITRGM) signature. Multi-omics analyses included bulk and single-cell RNA sequencing of melanoma patients, mouse bulk RNA sequencing, and pathology sections of melanoma patients.

Results

We identified 66 consensus immunotherapy prognostic genes (CITPGs) using WGCNA and differentially expressed genes (DEGs) from two melanoma cohorts. The CITPG-high group showed better prognosis and enriched immune activities. DEGs between CITPG-high and CITPG-low groups in the TCGA-SKCM cohort were analyzed in three additional melanoma cohorts using univariate Cox regression, resulting in 44 consensus genes. Using 101 machine learning algorithm combinations, we constructed the ITRGM signature based on seven model genes. The ITRGM outperformed 37 published signatures in predicting immunotherapy prognosis across the training cohort, three testing cohorts, and a meta-cohort. It effectively stratified patients into high-risk or low-risk groups for immunotherapy response. The low-risk group, with high levels of model genes, correlated with increased immune characteristics such as tumor mutation burden and immune cell infiltration, indicating immune-hot tumors with a better prognosis. The ITRGM's relationship with the tumor immune microenvironment was further validated in our experiments using pathology sections with GBP5, an important model gene, and CD8 IHC analysis. The ITRGM also predicted better immunotherapy response in eight cohorts, including urothelial carcinoma and stomach adenocarcinoma, indicating broad applicability.

Conclusions

The ITRGM signature is a stable and robust predictor for stratifying melanoma patients into 'immune-hot' and 'immune-cold' tumors, enhancing prognosis and response to immunotherapy.

Prognostic aging gene-based score for colorectal cancer: unveiling links to drug resistance, mutation burden, and personalized treatment strategies

OBJECTIVE

Colorectal cancer (CRC) is characterized by high incidence and mortality rates worldwide. In this study, we present a novel aging-related gene-based risk scoring system (Aging score) as a predictive tool for CRC prognosis.

METHOD

We identified prognostic aging-related genes using univariate Cox regression analysis, revealing key biological processes in CRC progression. We then constructed a robust prognostic model using LASSO and multivariate Cox regression analyses, including four critical genes: CAV1, FOXM1, MAD2L1, and WT1.

RESULT

The Aging score demonstrated high prognostic performance across the training, testing, and entire TCGA-CRC datasets, proving its reliability. High-risk patients identified by the Aging score had significantly shorter overall survival times than low-risk patients, indicating its potential for patient stratification and personalized treatment. The Aging score remained an independent prognostic factor compared to age, gender, and tumor stage. Additionally, the score was linked to tumor mutation burden and microsatellite instability, indicators of immune checkpoint inhibitor response. High-risk patients also showed higher estimated IC50 values for common chemotherapeutic drugs, suggesting possible treatment resistance.

CONCLUSION

Our findings highlight the Aging score's potential to enhance clinical decision-making and pave the way for personalized CRC management.

The synergistic effects of anoikis-related genes and EMT-related genes in the prognostic prediction of Wilms tumor

Wilms tumor (WT) is the most common type of malignant abdominal tumor in children; it exhibits a high degree of malignancy, grow rapidly, and is prone to metastasis. This study aimed to construct a prognosis model based on anoikis-related genes (ARGs) and epithelial-mesenchymal transition (EMT)-related genes (ERGs) for WT patients; we assessed the characteristics of the tumor microenvironment and treatment efficacy, as well as identifying potential therapeutic targets. To this end, we downloaded transcriptome sequencing data and clinical data for WT and normal renal cortices and used R to construct and validate the prognostic model based on ARGs and ERGs. Additionally, we performed clinical feature analysis, nomogram construction, mutation analysis, drug sensitivity analysis, Connectivity Map (cMAP) analysis, functional enrichment analysis, and immune infiltration analysis. Finally, we screened the hub gene using the STRING database and validated it via experiments. In this way, we constructed a model with good accuracy and robustness, which was composed of seven anoikis- and EMT-related genes. Paclitaxel and mesna were selected as potential chemotherapeutic drugs and adjuvant chemotherapeutic drugs for the WT high-risk group by using the Genomics of Drug Sensitivity in Cancer (GDSC) and cMAP compound libraries, respectively. We proved the existence of a strong correlation between invasive immune cells and prognostic genes and risk scores. Next, we selected NTRK2 as the hub gene, and in vitro experiments confirmed that its inhibition can significantly inhibit the proliferation and migration of tumor cells and promote late apoptosis. In summary, we screened out the potential biomarkers and chemotherapeutic drugs that can improve the prognosis of patients with WT.

Reliable estrogen-related prognostic signature for uterine corpus endometrial carcinoma

BACKGROUND

Uterine corpus endometrial carcinoma (UCEC) is a predominant gynecological malignancy worldwide. Overdosed estrogen exposure has been widely known as a crucial risk factor for UCEC patients. The purpose of this work is to explore crucial estrogen-related genes (ERGs) in UCEC.

METHODS

UCEC scRNA-seq data, bulk RNA data, and ERGs were obtained from GEO, TCGA, and Molecular Signature Database, respectively. Differential expression analysis and cross analysis determined the candidate genes, and optimal genes in risk score were obtained after univariate Cox regression analysis, LASSO Cox regression analysis, and multivariate Cox regression analysis. The functional information was revealed by GO, KEGG, and GSVA enrichment analyses. CCK8 assay was used to detect the drug sensitivity.

RESULTS

After cross analysis of the differentially expressed genes and the 8734 ERGs, 86 differentially expressed ERGs were identified in UCEC, which were significantly enriched in some immune related pathways and microbiota related pathways. Of them, the most optimal 8 ERGs were obtained to build prognostic risk score, including GAL, PHGDH, SLC7A2, HNMT, CLU, AREG, MACC1, and HMGA1. The risk score could reliably predict patient prognosis, and high-risk patients had worse prognosis. Higher HMGA1 gene expression exhibited higher sensitivity to Osimertinib.

CONCLUSIONS

Predictive risk score based on 8 ERGs exhibited excellent prognostic value in UCEC patients, and high-risk patients had inferior survival. UCEC patients with distinct prognoses showed different tumor immune microenvironment.

A novel prognostic model based on pyroptosis signature in AML

Acute myeloid leukemia (AML), a highly heterogeneous myeloid malignancy, remains a challenge in terms of proper risk stratification. In this study, we developed a novel pyroptosis prognostic model based on pyroptosis-related gene pairs, which exhibited excellent prognostic performance across multiple cohorts (N = 1506) and accurately predicted both adult and pediatric AML prognosis. Additionally, we integrated the pyroptosis risk model with other clinical risk factors to construct a highly operational nomogram. Moreover, our findings indicate a significant correlation between elevated pyroptosis risk scores and increased stemness of AML. Using CIBERSORT immune analysis, we found a decreased proportion of resting NK cells and activated mast cells in the high-risk group. Through analyzing the correlation between chemotherapy drug response sensitivity and risk scores, we found that AZD1332 and BPD-0008900 were extremely sensitive in the high-risk group.

Prognostic value of anoikis-related genes revealed using multi-omics analysis and machine learning based on lower-grade glioma features and tumor immune microenvironment

Background

The investigation explores the involvement of anoikis-related genes (ARGs) in lower-grade glioma (LGG), seeking to provide fresh insights into the disease's underlying mechanisms and to identify potential targets for therapy.

Methods

We applied unsupervised clustering techniques to categorize LGG patients into distinct molecular subtypes based on ARGs with prognostic significance. Additionally, various machine learning algorithms were employed to pinpoint genes most strongly correlated with patient outcomes, which were then used to develop and assess risk profiles.

Results

Our analysis identified two distinct molecular subtypes of LGG, each with significantly different prognoses. Patients in Cluster 2 had a median survival of 2.036 years, markedly shorter than the 7.994 years observed in Cluster 1 (P < 0.001). We also constructed a six-gene ARG signature that efficiently classified patients into two risk categories, showing median survival durations of 4.084 years for the high-risk group and 10.304 years for the low-risk group (P < 0.001). Significantly, the immune profiles, tumor mutation characteristics, and drug sensitivity varied greatly among these risk groups. The high-risk group was characterized by a "cold" tumor microenvironment (TME), a lower IDH1 mutation rate (61.7 % vs. 91.4 %), a higher TP53 mutation rate (53.7 % vs. 38.9 %), and greater sensitivity to targeted therapies such as QS11 and PF-562271. Furthermore, our nomogram, integrating risk scores with clinicopathological features, demonstrated strong predictive accuracy for clinical outcomes in LGG patients, with an AUC of 0.903 for the first year. The robustness of this prognostic model was further validated through internal cross-validation and across three external cohorts.

Conclusions

The evidence from our research suggests that ARGs could potentially serve as reliable indicators for evaluating immunotherapy effectiveness and forecasting clinical results in patients with LGG.

Energy metabolism-related GLUD1 contributes to favorable clinical outcomes of IDH-mutant glioma

BACKGROUND

Glioma is the most common brain tumor. IDH mutations occur frequently in glioma, indicating a more favorable prognosis. We aimed to explore energy metabolism-related genes in glioma to promote the research and treatment.

METHODS

Datasets were obtained from TCGA and GEO databases. Candidate genes were screened by differential gene expression analysis, then functional enrichment analysis was conducted on the candidate genes. PPI was also carried out to help determine the target gene. GSEA and DO analysis were conducted in the different expression level groups of the target gene. Survival analysis and immune cell infiltrating analysis were performed as well.

RESULTS

We screened 34 candidate genes and selected GLUD1 as the target gene. All candidate genes were significantly enriched in 10 KEGG pathways and 330 GO terms. GLUD1 expression was higher in IDH-mutant samples than IDH-wildtype samples, and higher in normal samples than tumor samples. Low GLUD1 expression was related to poor prognosis according to survival analysis. Most types of immune cells were negatively related to GLUD1 expression, but monocytes and activated mast cells exhibited significantly positive correlation with GLUD1 expression. GLUD1 expression was significantly related to 119 drugs and 6 immune checkpoint genes. GLUD1 was able to serve as an independent prognostic indicator of IDH-mutant glioma.

CONCLUSION

In this study, we identified an energy metabolism-related gene GLUD1 potentially contributing to favorable clinical outcomes of IDH-mutant glioma. In glioma, GLUD1 related clinical outcomes and immune landscape were clearer, and more valuable information was provided for immunotherapy.

Integrated machine learning algorithms identify KIF15 as a potential prognostic biomarker and correlated with stemness in triple-negative breast cancer

Cancer stem cells (CSCs) have the potential to self-renew and induce cancer, which may contribute to a poor prognosis by enabling metastasis, recurrence, and therapy resistance. Hence, this study was performed to identify the association between CSC-related genes and triple-negative breast cancer (TNBC) development. Stemness gene sets were downloaded from StemChecker. Based on the online databases, a consensus clustering algorithm was conducted for unsupervised classification of TNBC samples. The variations between subtypes were assessed with regard to prognosis, tumor immune microenvironment (TIME), and chemotherapeutic sensitivity. The stemness-related gene signature was established and random survival forest analysis was employed to identify the core gene for validation experiments and tumor sphere formation assays. 499 patients with TNBC were classified into three subgroups and the Cluster 1 had a better OS than others. After that, WGCNA study was performed to identify genes important for Cluster 1 subtype. Out of all 8 modules, the subtype of Cluster 1 and the yellow module with 103 genes demonstrated the largest positive association. After that, a four-gene stemness-related signature was established. Based on the yellow module, the 39 potential pivotal genes were subjected to the random forest survival analysis to find out the gene that was relatively important for OS. KIF15 was confirmed as the targeted gene by LASSO and random survival forest analyses. In vitro experiments, the downregulation of KIF15 promoted the stemness of TNBC cells. The expression levels of stem cell markers Nanog, SOX2, and OCT4 were found to be elevated in TNBC cell lines after KIF15 inhibition. A stemness-associated risk model was constructed to forecast the clinical outcomes of TNBC patients. The downregulation of KIF15 expression in a subpopulation of TNBC stem cells may promote stemness and possibly TNBC progression.

Comprehensive identification of a disulfidptosis-associated long non-coding RNA signature to predict the prognosis and treatment options in ovarian cancer

Purpose

Distinguished from cuproptosis and ferroptosis, disulfidptosis has been described as a newly discovered form of non-programmed cell death tightly associated with glucose metabolism. However, the prognostic profile of disulfidptosis-related lncRNAs (DRLRs) in ovarian cancer (OC) and their biological mechanisms need to be further elucidated.

Materials and methods

First, we downloaded the profiles of RNA transcriptome, clinical information for OC patients from the TCGA database. Generated from Cox regression analysis, prognostic lncRNAs were utilized to identify the risk signature by least absolute shrinkage and selection operator analysis. Then, we explored the intimate correlations between disulfidptosis and lncRNAs. What's more, we performed a series of systemic analyses to assess the robustness of the model and unravel its relationship with the immune microenvironment comprehensively.

Results

We identified two DRLR clusters, in which OC patients with low-risk scores exhibited a favorable prognosis, up-regulated immune cell infiltrations and enhanced sensitivity to immunotherapy. Furthermore, validation of the signature by clinical features and Cox analysis demonstrated remarkable consistency, suggesting the universal applicability of our model. It's worth noting that high-risk patients showed more positive responses to immune checkpoint inhibitors and potential chemotherapeutic drugs.

Conclusion

Our findings provided valuable insights into DRLRs in OC for the first time, which indicated an excellent clinical value in the selection of management strategies, spreading brilliant horizons into individualized therapy.

PANoptosis-related molecular clustering and prognostic signature associated with the immune landscape and therapy response in breast cancer

Breast cancer (BC) remains one of the most pervasive and complex malignancies. PANoptosis represents a recently identified cellular mechanism leading to programmed cell death. However, the prognostic implications and influence on the immune microenvironment of BC pertaining to PANoptosis-related genes (PRGs) remain significantly understudied. We conducted differential expression analysis to identify prognostic-Related PRGs by the Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. Next, we identified the PANoptosis-related molecular subtype using the consensus clustering analysis, and constructed and validated the PANoptosis-related prognostic signature using LASSO and Cox regression analyses. ROC curves were employed to assess the performance of the signatures. Furthermore, drug sensitivity between low- and high-risk group were analysis. Finally, we conducted RT-qPCR to assess the gene expression levels involved in this signature. We categorized BC patients into 2 distinct molecular clusters based on PRGs and identified differentially expressed genes associated with prognosis. Subsequently, BC patients were then divided into 2 gene clusters. The identified PRGs molecular clusters and gene clusters demonstrated association with patient survival, immune system functions, and biological processes and pathways of BC. A prognostic signature comprising 5 genes was established, and BC patients were classified into low- and high-risk groups based on the risk scores. The ROC curves demonstrated that those in the low-risk category exhibited notably extended survival compared to the high-risk group. A nomogram model for patient survival was constructed based on the risk score in conjunction with other clinical features. High-risk group had higher tumor burden mutation, CSC index and lower StomalScore, ImmuneScore, and ESTIMATEScore. Subsequently, we established a correlation between the risk score and drug sensitivity among BC patients. Finally, qRT-PCR results showed that the expression of CXCL1, PIGR, and TNFRSF14 significantly decreased, while CXCL13 and NKAIN were significantly increased in BC tissues. We have developed a molecular clustering and prognostic signature based on PANoptosis to improve the prediction of BC prognosis. This discovery has the potential to not only assist in assessing overall patient prognosis but also to deepen our understanding of the underlying mechanisms of PANoptosis in BC pathogenesis.

A scoring model for the expression of genes related to programmed cell death predicts immunotherapy response and prognosis in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) continues to be the leading cause of cancer death worldwide, driven by environmental factors like smoking and genetic predispositions. LUAD has a high mortality rate, and new biomarkers are urgently needed to improve treatment strategies and patient management. Programmed cell death (PCD) is involved in tumor progression and response to treatment. Therefore, there is a need for an extensive study of the role and functions of PCD-related genes (PCDRGs) in lung adenocarcinoma so as to understand the pathophysiologic features of lung adenocarcinoma.

METHODS

Based on TCGA and GEO databases, this research is aimed at screening differentially expressed PCD-related genes in lung adenocarcinoma. We conducted GO, and KEGG analysis to establish the link between these genes and biological processes. By applying various machine learning algorithms such as CoxBoost analysis, we developed PCD-related indices (PCDI) that were used to verify their ability to predict prognosis with the use of other datasets. This was done in addition to exploring the biological functions of PCD genes associated with lung adenocarcinoma by assessing the relationship between immune cell components of tumor microenvironment and PCD genes together with examining how they affect drug sensitivity.

RESULTS

The research presented in this article offers significant insights into LUAD. The authors identified 113 PCDRGs that were differentially expressed in LUAD. These genes are implicated in various biological functions, including High risk ing apoptosis, ferroptosis, and pathways specific to non-small cell lung cancer. Notably, the PCDI proved effective in distinguishing between High risk and Low risk LUAD patients, demonstrating a higher accuracy in prognosis prediction compared to traditional clinical indicators such as age and gender. This high prediction accuracy was validated in both test and validation cohorts. Additionally, these genes showed significant correlations with immune cell infiltration and drug sensitivity in LUAD patients.

CONCLUSION

We analysed the expression and function of PCDRGs in LUAD and revealed their correlation with patient survival, the immune microenvironment and drug sensitivity. The constructed PCDI model provides a scientific basis for the personalised treatment of lung adenocarcinoma, and future optimisation of treatment strategies based on these genes may improve patient clinical outcomes.

Integrative genomics unveils basement membrane-related diagnostic markers and therapeutic targets in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is often diagnosed at advanced stages due to the inherent limitations of current screening methodologies. Central to evaluating tumor invasion and prognostic assessment in ESCC is the integrity of the basement membrane (BM). However, current research on the implications of BM-related genes (BMRGs) in diagnosing ESCC remains sparse.

METHODS

We performed a comprehensive analysis using single-cell RNA-sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database, alongside gene expression profiles acquired from GEO and The Cancer Genome Atlas (TCGA) databases. This identified differentially expressed BMRGs in ESCC. Employing LASSO, RF, and SVM-RFE, we selected potential BM biomarkers and crafted a diagnostic nomogram for ESCC, validated by ROC curves and AUC values. We also explored immune infiltration and biological mechanisms through consensus clustering and GSVA, and utilized single cell trajectory analysis and GSCALite to study gene distributions and pathways. In vitro experiments further elucidated the role of these genes in ESCC carcinogenesis.

RESULTS

Here, we discovered that ESCC cell types exhibited markedly elevated BM-related scores. Our analysis pinpointed seven BM genes upregulated and linked to immune infiltration, showcasing unique gene expression profiles and varying immune cell densities across the BM-related subtypes. Furthermore, a robust positive correlation was observed between these genes expression and EMT activity. The knockdown of BGN significantly suppressed cell proliferation, migration, invasion, while also augmenting cell viability following chemotherapy drug treatment.

CONCLUSION

Our study identified seven key BMRGs (BGN, LAMB3, SPARC, MMP1, LUM, COL4A1, and NELL2) and established a diagnostic nomogram for ESCC. Of noteworthy significance is the discovery of BGN as a promising drug target, indicating a novel strategy for future clinical combination therapies in ESCC.

Regulation of Age-Related Lipid Metabolism in Ovarian Cancer

Although a lot of effort has been dedicated to ovarian cancer (OC) research, the mortality rate is still among the highest in female gynecologic malignancies. The effects of the aged tumor microenvironment are still being undermined despite age being the highest risk factor in ovarian cancer development and progression. In this study, we have conducted RNA sequencing and lipidomics analysis of gonadal adipose tissues from young and aged rat xenografts before and after ovarian cancer formation. We have found significantly higher tumor formation rates and volumes in aged OC xenograft rat models compared to their young counterparts (p<0.05), suggesting the aged adipose microenvironment (AME) is more susceptible to OC outgrowth. We have revealed significant shifts in the gene expression enrichment from groups of young vs. aged rats before tumor formation, groups of young vs. aged rats when the tumor formed, and groups of aged rats before and after tumor formation. We also observed shifts in the lipid components of the gonadal adipose tissues between young and aged rat xenografts when tumors were generated. Additionally, we found that the aged AME was associated with age-related changes in the immune cell composition, especially inflammation-related cells. The top hits showing the most differences between aged and young adipose tissues were eight genes including S100a8, S100a9, Il1rl1, Lcn2, C3, Hba-a1, Fcna, and Pnpla3, 22 lipids including multiple isoforms of free fatty acids (FFA) and triglyceride (TG), as well as four immune cells including neutrophil, myeloid dendritic cell, T cell CD4+ (non-regulatory), and mast cell activation. The functional correlation among S100a8, S100a9, neutrophil, and FFA (18:3) was also determined. Furthermore, FFA (18:3), which was shown to be downregulated in aged xenograft rats, was capable of inhibiting OC cell proliferation. In conclusion, our study suggested that aging promoted OC proliferation through changes in genes/pathways, lipid metabolism, and immune cells. Targeting the aging adipose microenvironment, particularly lipid metabolism reprogramming, holds promise as a therapeutic strategy for OC, which warrants further investigation.

Significance

Aging microenvironment of OC may be regulated by S100a8 and S100a9 secreted by adipocytes, preadipocytes, or neutrophils through affecting the lipid metabolism, such as FFA (18:3).

Machine learning-based autophagy-related prognostic signature for personalized risk stratification and therapeutic approaches in bladder cancer

OBJECTIVE

Bladder cancer (BCa) is a highly lethal urological malignancy characterized by its notable histological heterogeneity. Autophagy has swiftly emerged as a diagnostic and prognostic biomarker in diverse cancer types. Nonetheless, the currently accessible autophagy-related signature specific to BCa remains limited.

METHODS

A refined autophagy-related signature was developed through a 10-fold cross-validation framework, incorporating 101 combinations of machine learning algorithms. The performance of this signature in predicting prognosis and response to immunotherapy was thoroughly evaluated, along with an exploration of potential drug targets and compounds. In vitro and in vivo experiments were conducted to verify the regulatory mechanism of hub gene.

RESULTS

The autophagy-related prognostic signature (ARPS) has exhibited superior performance in predicting the prognosis of BCa compared to the majority of clinical features and other developed markers. Higher ARPS is associated with poorer prognosis and reduced sensitivity to immunotherapy. Four potential targets and five therapeutic agents were screened for patients in the high-ARPS group. In vitro and vivo experiments have confirmed that FKBP9 promotes the proliferation, invasion, and metastasis of BCa.

CONCLUSIONS

Overall, our study developed a valuable tool to optimize risk stratification and decision-making for BCa patients.

Immune-related cell death index and its application for hepatocellular carcinoma

Regulated cell death (RCD) plays a crucial role in the immune microenvironment, development, and progression of hepatocellular carcinoma (HCC). However, reliable immune-related cell death signatures have not been explored. In this study, we collected 12 RCD modes (e.g., apoptosis, ferroptosis, and cuproptosis), including 1078 regulators, to identify immune-related cell death genes based on HCC immune subgroups. Using a developed competitive machine learning framework, nine genes were screened to construct the immune-related cell death index (IRCDI), which is available for online application. Multi-omics data, along with clinical features, were analyzed to explore the HCC malignant heterogeneity. To validate the efficacy of this model, more than 18 independent cohorts, including survival and diverse treatment cohorts and datasets, were utilized. These findings were further validated using in-house samples and molecular biological experiments. Overall, the IRCDI may have a wide application in individual therapeutic decision-making and improving outcomes for HCC patients.

Identification of HDAC10 as a candidate oncogene in clear cell renal carcinoma that facilitates tumor proliferation and metastasis

BACKGROUND

Clear cell renal cell carcinoma (ccRCC) remains one of the most lethal urological malignancies even though a great number of improvements in diagnosis and management have achieved over the past few decades. Accumulated evidence revealed that histone deacetylases (HDACs) play vital role in cell proliferation, differentiation and apoptosis. Nevertheless, the biological functions of histone deacetylation modification related genes in ccRCC remains poorly understood.

METHOD

Bulk transcriptomic data and clinical information of ccRCC patients were obtained from the TCGA database and collected from the Chinese PLA General Hospital. A total of 36 histone deacetylation genes were selected and studied in our research. Univariate cox regression analysis, least absolute shrinkage and selection operator (LASSO) regression, random forest (RF) analysis, and protein-protein interaction (PPI) network analysis were applied to identify key genes affecting the prognosis of ccRCC. The 'oncoPredict' algorithm was utilized for drug-sensitive analysis. Gene Set Enrichment Analysis (GSEA) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was used to explore the potential biological function. The ssGSEA algorithm was used for tumor immune microenvironment analysis. The expression levels of HDAC10 were validated by RT-PCR and immunohistochemistry (IHC). 5-ethynyl-2'-deoxyuridine (EdU assay), CCK-8 assay, cell transwell migration and invasion assay and colony formation assay were performed to detect the proliferation and invasion ability of ccRCC cells. A nomogram incorporating HDAC10 and clinicopathological characteristics was established to predict the prognosis of ccRCC patients.

RESULT

Two machine learning algorithms and PPI analysis identified four histone deacetylation genes that have a significant association with the prognosis of ccRCC, with HDAC10 being the key gene among them. HDAC10 is highly expressed in ccRCC and its high expression is associated with poor prognosis for ccRCC patients. Pathway enrichment and the experiments of EdU staining, CCK-8 assay, cell transwell migration and invasion assay and colony formation assay demonstrated that HDAC10 mediated the proliferation and metastasis of ccRCC cells and involved in reshaping the tumor microenvironment (TME) of ccRCC. A clinically reliable prognostic predictive model was established by incorporating HDAC10 and other clinicopathological characteristics ( https://nomogramhdac10.shinyapps.io/HDAC10_Nomogram/ ).

CONCLUSION

Our study found the increased expression of HDAC10 was closely associated with poor prognosis of ccRCC patients. HDAC10 showed a pro-tumorigenic effect on ccRCC and promote the proliferation and metastasis of ccRCC, which may provide new light on targeted therapy for ccRCC.

mtPCDI: a machine learning-based prognostic model for prostate cancer recurrence

Background

This research seeks to formulate a prognostic model for forecasting prostate cancer recurrence by examining the interaction between mitochondrial function and programmed cell death (PCD).

Methods

The research involved analyzing four gene expression datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) using univariate Cox regression. These analyses identified genes linked with mitochondrial function and PCD that correlate with recurrence prognosis. Various machine learning algorithms were then employed to construct an optimal predictive model.

Results

A key outcome was the creation of a mitochondrial-related programmed cell death index (mtPCDI), which effectively predicts the prognosis of prostate cancer patients. It was observed that individuals with lower mtPCDI exhibited higher immune activity, correlating with better recurrence outcomes.

Conclusion

The study demonstrates that mtPCDI can be used for personalized risk assessment and therapeutic decision-making, highlighting its clinical significance and providing insights into the biological processes affecting prostate cancer recurrence.

Integrative analysis of bulk and single-cell RNA sequencing reveals sphingolipid metabolism and immune landscape in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is characterized by its aggressive behavior and complex molecular heterogeneity, posing significant challenges for treatment and prognostication. This study offers a comprehensive analysis of ccRCC by leveraging both bulk and single-cell RNA sequencing data, with a specific aim to unravel the complexities of sphingolipid metabolism and the intricate dynamics within the tumor microenvironment (TME). By examining ccRCC samples sourced from public databases, our investigation delves deep into the genetic and transcriptomic landscape of this cancer type. Employing advanced analytical techniques, we have identified pivotal patterns in gene expression and cellular heterogeneity, with a special focus on the roles and interactions of various immune cells within the TME. Significantly, our research has unearthed insights into the dynamics of sphingolipid metabolism in ccRCC, shedding light on its potential implications for tumor progression and strategies for immune evasion. A novel aspect of this study is the development of a risk score model designed to enhance prognostic predictions for ccRCC patients, which is currently pending external validation to ascertain its clinical utility. Despite its contributions, the study is mindful of its limitations, including a reliance on observational data from public sources and a primary focus on RNA sequencing data, which may constrain the depth and generalizability of the findings. The study does not encompass critical aspects, such as protein expression, posttranslational modifications, and comprehensive metabolic profiles. Moreover, its retrospective design underscores the necessity for future prospective studies to solidify these preliminary conclusions. Our findings illuminate the intricate interplay between genetic alterations, sphingolipid metabolism, and immune responses in ccRCC. This research not only enhances our understanding of the molecular foundations of ccRCC but also paves the way for the development of targeted therapies and personalized treatment modalities. The study underlines the importance of cautious interpretation of results and champions ongoing research using diverse methodologies to thoroughly comprehend and effectively combat this formidable cancer type.

Characterization of Mitoribosomal Small Subunit unit genes related immune and pharmacogenomic landscapes in renal cell carcinoma

Mitoribosomes are essential for the production of biological energy. The Human Mitoribosomal Small Subunit unit (MRPS) family, responsible for encoding mitochondrial ribosomal small subunits, is actively engaged in protein synthesis within the mitochondria. Intriguingly, MRPS family genes appear to play a role in cancer. A multistep process was employed to establish a risk model associated with MRPS genes, aiming to delineate the immune and pharmacogenomic landscapes in clear cell renal cell carcinoma (ccRCC). MRPScores were computed for individual patients to assess their responsiveness to various treatment modalities and their susceptibility to different therapeutic targets and drugs. While MRPS family genes have been implicated in various cancers as oncogenes, our findings reveal a contrasting tumor suppressor role for MRPS genes in ccRCC. Utilizing an MRPS-related risk model, we observed its excellent prognostic capability in predicting survival outcomes for ccRCC patients. Remarkably, the subgroup with high MRPS-related scores (MRPScore) displayed poorer prognosis but exhibited a more robust response to immunotherapy. Through in silico screening of 2183 drug targets and 1646 compounds, we identified two targets (RRM2 and OPRD1) and eight agents (AZ960, carmustine, lasalocid, SGI-1776, AZD8055_1059, BPD.00008900_1998, MK.8776_2046, and XAV939_1268) with potential therapeutic implications for high-MRPScore patients. Our study represents the pioneering effort in proposing that molecular classification, diagnosis, and treatment strategies can be formulated based on MRPScores. Indeed, a high MRPScore profile appears to elevate the risk of tumor progression and mortality, potentially through its influence on immune regulation. This suggests that the MRPS-related risk model holds promise as a prognostic predictor and may offer novel insights into personalized therapeutic strategies.

Identification and validation of an immune-derived multiple programmed cell death index for predicting clinical outcomes, molecular subtyping, and drug sensitivity in lung adenocarcinoma

OBJECTIVES

Comprehensive cross-interaction of multiple programmed cell death (PCD) patterns in the patients with lung adenocarcinoma (LUAD) have not yet been thoroughly investigated.

METHODS

Here, we collected 19 different PCD patterns, including 1911 PCD-related genes, and developed an immune-derived multiple programmed cell death index (MPCDI) based on machine learning methods.

RESULTS

Using the median MPCDI scores, we categorized the LUAD patients into two groups: low-MPCDI and high-MPCDI. Our analysis of the TCGA-LUAD training cohort and three external GEO cohorts (GSE37745, GSE30219, and GSE68465) revealed that patients with high-MPCDI experienced a more unfavorable prognosis, whereas those with low-MPCDI had a better prognosis. Furthermore, the results of both univariate and multivariate Cox regression analyses further confirmed that MPCDI serves as a novel independent risk factor. By combining clinical characteristics with the MPCDI, we constructed a nomogram that provides an accurate and reliable quantitative tool for personalized clinical management of LUAD patients. The findings obtained from the analysis of C-index and the decision curve revealed that the nomogram outperformed various clinical variables in terms of net clinical benefit. Encouragingly, the low-MPCDI patients are more sensitive to commonly used chemotherapy drugs, which suggests that MPCDI scores have a guiding role in chemotherapy for LUAD patients.

CONCLUSION

Therefore, MPCDI can be used as a novel clinical diagnostic classifier, providing valuable insights into the clinical management and clinical decision-making for LUAD patients.

Mesenchymal-like immune-altered is the fourth robust triple-negative breast cancer molecular subtype

BACKGROUND

Robust molecular subtyping of triple-negative breast cancer (TNBC) is a prerequisite for the success of precision medicine. Today, there is a clear consensus on three TNBC molecular subtypes: luminal androgen receptor (LAR), basal-like immune-activated (BLIA), and basal-like immune-suppressed (BLIS). However, the debate about the robustness of other subtypes is still open.

METHODS

An unprecedented number (n = 1942) of TNBC patient data was collected. Microarray- and RNAseq-based cohorts were independently investigated. Unsupervised analyses were conducted using k-means consensus clustering. Clusters of patients were then functionally annotated using different approaches. Prediction of response to chemotherapy and targeted therapies, immune checkpoint blockade, and radiotherapy were also screened for each TNBC subtype.

RESULTS

Four TNBC subtypes were identified in the cohort: LAR (19.36%); mesenchymal stem-like (MSL/MES) (17.35%); BLIA (31.06%); and BLIS (32.23%). Regarding the MSL/MES subtype, we suggest renaming it to mesenchymal-like immune-altered (MLIA) to emphasize its specific histological background and nature of immune response. Treatment response prediction results show, among other things, that despite immune activation, immune checkpoint blockade is probably less or completely ineffective in MLIA, possibly caused by mesenchymal background and/or an enrichment in dysfunctional cytotoxic T lymphocytes. TNBC subtyping results were included in the bc-GenExMiner v5.0 webtool ( http://bcgenex.ico.unicancer.fr ).

CONCLUSION

The mesenchymal TNBC subtype is characterized by an exhausted and altered immune response, and resistance to immune checkpoint inhibitors. Consensus for molecular classification of TNBC subtyping and prediction of cancer treatment responses helps usher in the era of precision medicine for TNBC patients.