Comprehensive Molecular Analyses of a Novel Mutational Signature Classification System with Regard to Prognosis, Genomic Alterations, and Immune Landscape in Glioma

IDH mutation, glioma immunogenicity, and therapeutic challenge of primary mismatch repair deficient IDH-mutant astrocytoma PMMRDIA: a systematic review

In 2021, Suwala et al. described Primary Mismatch Repair Deficient IDH-mutant Astrocytoma (PMMRDIA) as a distinct group of gliomas. In unsupervised clustering, PMMRDIA forms distinct cluster, separate from other IDH-mutant gliomas, including IDH-mutant gliomas with secondary mismatch repair (MMR) deficiency. In the published cohort, three patients received treatment with an immune checkpoint blocker (ICB), yet none exhibited a response, which aligns with existing knowledge about the decreased immunogenicity of IDH-mutant gliomas in comparison to IDH-wildtype. In the case of PMMRDIA, the inherent resistance to the standard-of-care temozolomide caused by MMR deficiency is an additional challenge. It is known that a gain-of-function mutation of IDH1/2 genes produces the oncometabolite R-2-hydroxyglutarate (R-2-HG), which increases DNA and histone methylation contributing to the characteristic glioma-associated CpG island methylator phenotype (G-CIMP). While other factors could be involved in remodeling the tumor microenvironment (TME) of IDH-mutant gliomas, this systematic review emphasizes the role of R-2-HG and the subsequent G-CIMP in immune suppression. This highlights a potential actionable pathway to enhance the response of ICB, which might be relevant for addressing the unmet therapeutic challenge of PMMRDIA.

Experimental validation and pan-cancer analysis identified COL10A1 as a novel oncogene and potential therapeutic target in prostate cancer

BACKGROUND

Type X collagen (COL10) is a homologous trimeric non-fibrillar collagen found in the extracellular matrix of human tissues, and it exhibits a distinctive white appearance. Type X collagen α1 chain (COL10A1) is a specific cleaved fragment of type X collagen. However, the expression, prognostic significance, clinicopathological attributes and immune-related associations of COL10A1 in prostate cancer as well as in pan-cancer contexts remain poorly understood.

METHODS

Using bioinformatic analysis of data from the most recent databases (TCGA, GTEx and GEO databases), we have extensively elucidated the role played by COL10A1 in terms of its expression patterns, prognostic implications, and immune efficacy across a pan-cancer spectrum. Subsequently, the biological functions of COL10A1 in prostate cancer were elucidated by experimental validation.

RESULTS

Our findings have confirmed that COL10A1 was highly expressed in most cancers and was associated with poorer prognosis in cancer patients. Immune correlation analysis of COL10A1 in various cancers showed its significant correlation with Tumor mutational burden (TMB), microsatellite instability (MSI) and immune cell infiltration. In addition, knockdown of COL10A1 in prostate cancer resulted in a substantial reduction in the proliferation, migration, and invasive potential of prostate cancer cells.

CONCLUSION

Our pan-cancer analysis of COL10A1 gene provided novel insights into its pivotal role in cancer initiation, progression, and therapeutic implications, underscoring its potential significance in prognosis and immunotherapeutic interventions for cancer, particularly prostate cancer.

A machine learning-derived gene signature for assessing rupture risk and circulatory immunopathologic landscape in patients with intracranial aneurysms

Background

Intracranial aneurysm (IA) is an uncommon but severe subtype of cerebrovascular disease, with high mortality after aneurysm rupture. Current risk assessments are mainly based on clinical and imaging data. This study aimed to develop a molecular assay tool for optimizing the IA risk monitoring system.

Methods

Peripheral blood gene expression datasets obtained from the Gene Expression Omnibus were integrated into a discovery cohort. Weighted gene co-expression network analysis (WGCNA) and machine learning integrative approaches were utilized to construct a risk signature. QRT-PCR assay was performed to validate the model in an in-house cohort. Immunopathological features were estimated using bioinformatics methods.

Results

A four-gene machine learning-derived gene signature (MLDGS) was constructed for identifying patients with IA rupture. The AUC of MLDGS was 1.00 and 0.88 in discovery and validation cohorts, respectively. Calibration curve and decision curve analysis also confirmed the good performance of the MLDGS model. MLDGS was remarkably correlated with the circulating immunopathologic landscape. Higher MLDGS scores may represent higher abundance of innate immune cells, lower abundance of adaptive immune cells, and worse vascular stability.

Conclusions

The MLDGS provides a promising molecular assay panel for identifying patients with adverse immunopathological features and high risk of aneurysm rupture, contributing to advances in IA precision medicine.

A novel pyroptosis-regulated gene signature for predicting prognosis and immunotherapy response in hepatocellular carcinoma

Background: Pyroptosis, a newly discovered type of programmed cell death, has both anti-tumor and tumor-promoting effects on carcinogenesis. In hepatocellular carcinoma (HCC), however, the associations between pyroptosis-regulated genes and prognosis, immune microenvironment, and immunotherapy response remain unclear. Samples and methods: Sequencing data were collected from The Cancer Genome Atlas database, The International Cancer Genome Consortium (ICGC), and The Integrative Molecular Database of Hepatocellular Carcinoma (HCCDB). First, we investigated the expression levels and copy number variations (CNVs) of 56 pyroptosis genes in HCC and pan-cancer. Next, we identified 614 genes related to 56 pyroptosis-associated genes at the expression, mutation, and CNVs levels. Pathway enrichment analysis of 614 genes in the Hallmark, KEGG, and Reactome databases yielded a total of 253 significant signaling pathways. The pyroptosis-regulated genes (PRGs) comprised 108 genes that were derived from the top 20 signaling pathways, of which 57 genes had prognostic value in HCC. The least absolute shrinkage and selection operator (LASSO) analysis was performed to screen for PRGs with prognostic values. Ultimately, we constructed a risk score model with seven PRGs to predict HCC prognosis and validated its predictive value in three independent HCC cohorts. Risk scores were used to illustrate receiver operating characteristic (ROC) curves predicting 1, 3, and 5-years overall survival (OS). Single-sample gene set enrichment analysis (ssGSEA), was performed to study 28 types of immune cells infiltrated in HCC. The relationship between the risk signature and six immune checkpoint genes and immunotherapy was analyzed. Results: A total of seven PRGs were obtained following multiple screening steps. The risk score model containing seven PRGs was found to correlate significantly with the HCC prognosis of the training group. In addition, we validated the risk score model in two additional HCC cohorts. The risk score significantly correlated with infiltrating immune cells (i. e. CD4+ T cells, etc.), ICB key molecules (i. e. HAVCR2, etc.), and ICB response. Conclusions: This study demonstrated a vital role of PRGs in predicting the prognosis and immunotherapy response of HCC patients. The risk model could pave the way for drugs targeting pyroptosis and immune checkpoints in HCC.

Systematic transcriptome profiling of pyroptosis related signature for predicting prognosis and immune landscape in lower grade glioma

Background

Pyroptosis is a programmed cell death mediated by the gasdermin superfamily, accompanied by inflammatory and immune responses. Exogenously activated pyroptosis is still not well characterized in the tumor microenvironment. Furthermore, whether pyroptosis-related genes (PRGs) in lower-grade glioma (LGG) may be used as a biomarker remains unknown.

Methods

The RNA-Sequencing and clinical data of LGG patients were downloaded from publicly available databases. Bioinformatics approaches were used to analyze the relationship between PRGs and LGG patients’ prognosis, clinicopathological features, and immune status. The NMF algorithm was used to differentiate phenotypes, the LASSO regression model was used to construct prognostic signature, and GSEA was used to analyze biological functions and pathways. The expression of the signature genes was verified using qRT-PCR. In addition, the L1000FWD and CMap tools were utilized to screen potential therapeutic drugs or small molecule compounds and validate their effects in glioma cell lines using CCK-8 and colony formation assays.

Results

Based on PRGs, we defined two phenotypes with different prognoses. Stepwise regression analysis was carried out to identify the 3 signature genes to construct a pyroptosis-related signature. After that, samples from the training and test cohorts were incorporated into the signature and divided by the median RiskScore value (namely, Risk-H and Risk-L). The signature shows excellent predictive LGG prognostic power in the training and validation cohorts. The prognostic signature accurately stratifies patients according to prognostic differences and has predictive value for immune cell infiltration and immune checkpoint expression. Finally, the inhibitory effect of the small molecule inhibitor fedratinib on the viability and proliferation of various glioma cells was verified using cell biology-related experiments.

Conclusion

This study developed and validated a novel pyroptosis-related signature, which may assist instruct clinicians to predict the prognosis and immunological status of LGG patients more precisely. Fedratinib was found to be a small molecule inhibitor that significantly inhibits glioma cell viability and proliferation, which provides a new therapeutic strategy for gliomas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09982-7.

Hallmark-guided subtypes of hepatocellular carcinoma for the identification of immune-related gene classifiers in the prediction of prognosis, treatment efficacy, and drug candidates

Hepatocellular carcinoma (HCC), accounting for ~90% of all primary liver cancer, is a prevalent malignancy worldwide. The intratumor heterogeneity of its causative etiology, histology, molecular landscape, and immune phenotype makes it difficult to precisely recognize individuals with high mortality risk or tumor-intrinsic treatment resistance, especially immunotherapy. Herein, we comprehensively evaluated the activities of cancer hallmark gene sets and their correlations with the prognosis of HCC patients using gene set variation analysis (GSVA) and identified two HCC subtypes with distinct prognostic outcomes. Based on these subtypes, seven immune-related genes (TMPRSS6, SPP1, S100A9, EPO, BIRC5, PLXNA1, and CDK4) were used to construct a novel prognostic gene signature [hallmark-guided subtypes-based immunologic signature (HGSIS)] via multiple statistical approaches. The HGSIS-integrated nomogram suggested an enhanced predictive performance. Interestingly, oncogenic hallmark pathways were significantly enriched in the high-risk group and positively associated with the risk score. Distinct mutational landscapes and immune profiles were observed between different risk groups. Moreover, immunophenoscore (IPS) and tumor immune dysfunction and exclusion (TIDE) analysis showed different sensitivities of HGSIS risk groups for immune therapy efficacy, and the pRRophetic algorithm indicated distinguishable responses for targeted/chemotherapies in different groups. KIF2C was picked out as the key target concerning HGSIS, and the top 10 small molecules were predicted to bind to the active site of KIF2C via molecular docking, which might be further used for candidate drug discovery of HCC. Taken together, our study offers novel insights for clinically significant subtype recognition, and the proposed signature may be a helpful guide for clinicians to improve the treatment regimens.

Construction of a novel signature and prediction of the immune landscape in gastric cancer based on necroptosis-related genes

Necroptosis, a type of programmed cell death, has become a potential therapeutic target for solid tumors. Nevertheless, the potential roles of necroptosis-related genes (NRGs) in gastric cancer (GC) remain unknown. The objective of the present study was to create a necroptosis-related prognostic signature that can provide more accurate assessment of prognosis in GC. Using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) data, we identified differentially expressed NRGs. Univariate analysis and Lasso regression were performed to determine the prognostic signature. Risk scores were calculated and all GC patients were divided into high- and low-risk score group according to the median risk score value. The robustness of this signature was externally validated with data from GSE84437 cohort (n = 431). Survival analysis revealed high-risk score patients had a worse prognosis. Results evidenced that the signature was an independent prognosis factor for survival. Single-sample sequence set enrichment analysis (ssGSEA) exhibited different enrichment of immune cells and immune-related pathways in the two risk groups. Furthermore, a predictive nomogram was generated and showed excellent predictive performance based on discrimination and calibration. In addition, the risk score positively correlated with tumor mutational burden and was associated with sensitivity to multiple anti-cancer drugs. Overall, our work demonstrates a close relationship between necroptosis and the prognosis of GC. The signature we constructed with potential clinical application value, can be used for prognosis prediction and being a potential therapeutic responses indicator in GC patients.

A Novel Artificial Neural Network Prognostic Model Based on a Cancer-Associated Fibroblast Activation Score System in Hepatocellular Carcinoma

Introduction

Hepatocellular carcinoma (HCC) ranks fourth as the most common cause of cancer-related death. It is vital to identify the mechanism of progression and predict the prognosis for patients with HCC. Previous studies have found that cancer-associated fibroblasts (CAFs) promote tumor proliferation and immune exclusion. However, the information about CAF-related genes is still elusive.

Methods

The data were obtained from The Cancer Genome Atlas, International Cancer Genome Consortium, and Gene Expression Omnibus databases. On the basis of single-cell transcriptome and ligand–receptor interaction analysis, CAF-related genes were selected. By performing Cox regression and random forest, we filtered 12 CAF-related prognostic genes for the construction of the ANN model based on the CAF activation score (CAS). Then, functional, immune, mutational, and clinical analyses were performed.

Results

We constructed a novel ANN prognostic model based on 12 CAF-related prognostic genes. Cancer-related pathways were enriched, and higher activated cell crosstalk was identified in high-CAS samples. High immune activity was observed in high-CAS samples. We detected three differentially mutated genes (NBEA, RYR2, and FRAS1) between high- and low-CAS samples. In clinical analyses, we constructed a nomogram to predict the prognosis of patients with HCC. 5-Fluorouracil had higher sensitivity in high-CAS samples than in low-CAS samples. Moreover, some small-molecule drugs and the immune response were predicted.

Conclusion

We constructed a novel ANN model based on CAF-related genes. We revealed information about the ANN model through functional, mutational, immune, and clinical analyses.

The Identification of Necroptosis-Related Subtypes, the Construction of a Prognostic Model, and the Characterization of the Tumor Microenvironment in Gliomas

Necroptosis is a recently discovered form of cell death that plays a vital role in the progression of cancer, the spread of metastases, and the immunologic response to tumors. Due to the dual role of necrotic apoptotic processes in tumor pathogenesis and the heterogeneity of gliomas, the function of necroptosis in the glioma microenvironment is still poorly understood. We characterized the expression of necroptosis-related genes (NRGs) within glioma samples at both the genetic and transcriptional levels, identifying three distinct subtypes. Additionally, we constructed a risk score, which is capable of accurately predicting patient prognosis, correlates with tumor mutation burden (TMB), tumor stem cell index (CSC), immune checkpoints, and predicts tumor drug sensitivity. To facilitate its application in the clinic, we developed a nomogram and demonstrated that it predicts the prognosis of glioma patients with good accuracy and reliability using multiple datasets. We examined the function of necroptosis in the tumor microenvironment (TME) and the prognosis of gliomas, which may be useful for guiding individualized treatment plans for gliomas targeting necroptosis.

Title: Multi-Omics and Immune Landscape of Proliferative LncRNA Signatures: Implications for Risk Stratification and Immunotherapy in Hepatocellular Carcinoma

Background: Long noncoding RNAs (lncRNAs) are significantly implicated in tumor proliferation. Nevertheless, proliferation-derived lncRNAs and their latent clinical significance remain largely unrevealed in hepatocellular carcinoma (HCC).

Methods: This research enrolled 658 HCC patients from five independent cohorts. We retrieved 50 Hallmark gene sets from the MSigDB portal. Consensus clustering was applied to identify heterogeneous proliferative subtypes, and the nearest template prediction (NTP) was utilized to validate the subtypes. We introduced an integrative framework (termed “ProLnc”) to identify proliferation-derived lncRNAs. Moreover, a proliferation-related signature was developed and verified in four independent cohorts.

Results: In 50 Hallmarks, seven proliferation pathways were significantly upregulated and correlated with a worse prognosis. Subsequently, we deciphered two heterogeneous proliferative subtypes in TCGA-LIHC. Subtype 2 displayed enhanced proliferative activities and a worse prognosis, whereas subtype 1 was associated with hyperproliferative HCC and a favorable prognosis. The NTP further verified the robustness and reproducibility of two subtypes in four cohorts derived from different platforms. Combining the differentially expressed lncRNAs from two subtypes with proliferative lncRNA modulators from our ProLnc pipeline, we determined 230 proliferation-associated lncRNAs. Based on the bootstrapping channel and the verification of multiple cohorts, we further identified ten lncRNAs that stably correlated with prognosis. Subsequently, we developed and validated a proliferative lncRNA signature (ProLncS) that could independently and accurately assess the overall survival (OS) and relapse-free survival (RFS) of HCC patients in the four cohorts. Patients with high ProLncS score displayed significantly genomic alterations (e.g., TP53 mutation, 8p23-8p24 copy number variation) and higher abundances of immune cells and immune checkpoint molecules, which suggested immunotherapy was more suitable for patients with high ProLncS score.

Conclusion: Our work provided new insights into the heterogeneity of tumor proliferation, and ProLncS could be a prospective tool for tailoring the clinical decision and management of HCC.

Integrative Analysis Revealed Stemness Features and a Novel Stemness-Related Classification in Colorectal Cancer Patients

Cancer stem cells play crucial roles in colorectal cancer (CRC) tumorigenesis and treatment response. This study aimed to determine the value of the mRNA stemness index (mRNAsi) in CRC and introduce a stemness-related classification to predict the outcome of patients. mRNAsi scores and RNA sequence data of CRC patients were analyzed. We found that high mRNAsi scores were related to early-stage CRC and a better patient prognosis. Two stemness-based subtypes (subtype I and II) were identified. Patients in subtype I presented a significantly better prognosis than those in subtype II. Patients in these two subtype groups presented significantly different tumor immunity scores and immune cell infiltration patterns. Genomic variations revealed that patients in subtype I had a lower tumor mutation burden than those in subtype II. A three-gene stemness subtype predictor was established, showing good diagnostic value in discriminating patients in different subtypes. A prognostic signature based on five stemness-related genes was established and validated in two independent cohorts and clinical samples, showing a better predictive performance than other clinical parameters. We concluded that mRNAsi scores were associated with the clinical outcome in CRC patients. The stemness-related classification was a promising prognostic predictor for CRC patients.

Transcriptome-Based Dissection of Intracranial Aneurysms Unveils an “Immuno-Thermal” Microenvironment and Defines a Pathological Feature-Derived Gene Signature for Risk Estimation

Immune inflammation plays an essential role in the formation and rupture of intracranial aneurysm (IA). However, the current limited knowledge of alterations in the immune microenvironment of IA has hampered the mastery of pathological mechanisms and technological advances, such as molecular diagnostic and coated stent-based molecular therapy. In this study, seven IA datasets were enrolled from the GEO database to decode the immune microenvironment and relevant biometric alterations. The ssGSEA algorithm was employed for immune infiltration assessment. IAs displayed abundant immune cell infiltration, activated immune-related pathways, and high expression of immune-related genes. Several immunosuppression cells and genes were also coordinately upregulated in IAs. Five immune-related hub genes, including CXCL10, IL6, IL10, STAT1, and VEGFA, were identified from the protein-protein interaction network and further detected at the protein level. CeRNA networks and latent drugs targeting the hub genes were predicted for targeted therapy reference. Two gene modules recognized via WCGNA were functionally associated with contractile smooth muscle loss and extracellular matrix metabolism, respectively. In blood datasets, a pathological feature-derived gene signature (PFDGS) for IA diagnosis and rupture risk prediction was established using machine learning. Patients with high PFDGS scores may possess adverse biological alterations and present with a high risk of morbidity or IA rupture, requiring more vigilance or prompt intervention. Overall, we systematically unveiled an “immuno-thermal” microenvironment characterized by co-enhanced immune activation and immunosuppression in IA, which provides a novel insight into molecular pathology. The PFDGS is a promising signature for optimizing risk surveillance and clinical decision-making in IA patients.

Identification and Validation of Immune Molecular Subtypes and Immune Landscape Based on Colon Cancer Cohort

Background

The incidence and mortality rates of colon adenocarcinoma (COAD), which is the fourth most diagnosed cancer worldwide, are high. A subset of patients with COAD has shown promising responses to immunotherapy. However, the percentage of patients with COAD benefiting from immunotherapy is unclear. Therefore, gaining a better understanding of the immune milieu of colon cancer could aid in the development of immunotherapy and suitable combination strategies.

Methods

In this study, gene expression profiles and clinical follow-up data were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and molecular subtypes were identified using the ConsensusClusterPlus package in R. Univariate and multivariate Cox regression analyses were performed to evaluate the prognostic value of immune subtypes. The graph structure learning method was used to reduce the dimension to reveal the internal structure of the immune system. Weighted correlation network analysis (WGCNA) was performed to identify immune-related gene modules. Finally, western blotting was performed to verify the gene expression patterns in COAD samples.

Results

The results showed that 424 COAD samples could be divided into three subtypes based on 1921 immune cell-related genes, with significant differences in prognosis between subtypes. Furthermore, immune-related genes could be divided into five functional modules, each with a different distribution pattern of immune subtypes. Immune subtypes and gene modules were highly reproducible across many data sets. There were significant differences in the distribution of immune checkpoints, molecular markers, and immune characteristics among immune subtypes. Four core genes, namely, CD2, FGL2, LAT2, and SLAMF1, with prognostic significance were identified by WGCNA and univariate Cox analysis.

Conclusion

Overall, this study provides a conceptual framework for understanding the tumor immune microenvironment of colon cancer.

Identification of a Tumor Immunological Phenotype-Related Gene Signature for Predicting Prognosis, Immunotherapy Efficacy, and Drug Candidates in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the predominant subtype of primary liver cancer and represents a highly heterogeneous disease, making it hard to predict the prognosis and therapy efficacy. Here, we established a novel tumor immunological phenotype-related gene index (TIPRGPI) consisting of 11 genes by Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) algorithm to predict HCC prognosis and immunotherapy response. TIPRGPI was validated in multiple datasets and exhibited outstanding performance in predicting the overall survival of HCC. Multivariate analysis verified it as an independent predictor and a TIPRGPI-integrated nomogram was constructed to provide a quantitative tool for clinical practice. Distinct mutation profiles, hallmark pathways, and infiltration of immune cells in tumor microenvironment were shown between the TIPRGPI high and low-risk groups. Notably, significant differences in tumor immunogenicity and tumor immune dysfunction and exclusion (TIDE) were observed between the two risk groups, suggesting a better response to immune checkpoint blockade (ICB) therapy of the low-risk group. Besides, six potential drugs binding to the core target of the TIPRGPI signature were predicted via molecular docking. Taken together, our study shows that the proposed TIPRGPI was a reliable signature to predict the risk classification, immunotherapy response, and drugs candidate with potential application in the clinical decision and treatment of HCC. The novel “TIP genes”-guided strategy for predicting the survival and immunotherapy efficacy, we reported here, might be also applied to more cancers other than HCC.

Construction of a Lactate-Related Prognostic Signature for Predicting Prognosis, Tumor Microenvironment, and Immune Response in Kidney Renal Clear Cell Carcinoma

Kidney renal clear cell carcinoma (KIRC) is one of the most prevalent primary malignancies with high heterogeneity in the urological system. Growing evidence implies that lactate is a significant carbon source for cell metabolism and plays a vital role in tumor development, maintenance, and therapeutic response. However, the global influence of lactate-related genes (LRGs) on prognostic significance, tumor microenvironment characteristics, and therapeutic response has not been comprehensively elucidated in patients with KIRC. In the present study, we collected RNA sequencing and clinical data of KIRC from The Cancer Genome Atlas (TCGA), E-MTAB-1980, and GSE22541 cohorts. Unsupervised clustering of 17 differentially expressed LRG profiles divided the samples into three clusters with distinct immune characteristics. Three genes (FBP1, HADH, and TYMP) were then identified to construct a lactate-related prognostic signature (LRPS) using the least absolute shrinkage and selection operator (LASSO) and Cox regression analyses. The novel signature exhibited excellent robustness and predictive ability for the overall survival of patients. In addition, the constructed nomogram based on the LRPS-based risk scores and clinical factors (age, gender, tumor grade, and stage) showed a robust predictive performance. Furthermore, patients classified by risk scores had distinguishable immune status, tumor mutation burden, response to immunotherapy, and sensitivity to drugs. In conclusion, we developed an LRPS for KIRC that was closely related to the immune landscape and therapeutic response. This LRPS may guide clinicians to make more precise and personalized treatment decisions for KIRC patients.

Integrated Analysis of Necroptosis-Related Genes for Prognosis, Immune Microenvironment Infiltration, and Drug Sensitivity in Colon Cancer

Background

Necroptosis, is intimately linked to tumor development and prognosis and has been considered as a target for anticancer therapy. However, the role of necroptosis-related genes (NRGs) in colon cancer is unclear.

Methods

In the present study, we screened 76 NRGs from previous studies and described the landscape of transcriptomic and genetic variation of NRGs in colon cancer (CC) patient samples. Molecular subtypes of necroptosis in colon cancer were identified by clustering analysis, and these molecular subtypes were linked to patient prognosis and TME cell infiltration characteristics. Then, the NRS-score for predicting overall survival (OS) was built based on the TCGA database and validated in the GSE39582 cohort for its predictive power in CC patients. Besides, the ESTIMATE and CIBERSORT algorithms were applied to explore the relationship between NRS-score and tumor immune microenvironment.

Results

We identified two molecular subtypes associated with necroptosis in CC, which have diverse prognosis and immune microenvironment characteristics. Based on the differentially expressed genes between the two molecular subtypes, we further developed a necroptosis risk score signature, referred to as NRS-score. High NRS-score was associated with poor prognosis in CC through immunosuppressive microenvironment and immune escape mechanisms. The nomogram based on NRS-score showed excellent ability to predict prognosis. In addition, NRS-score presented a positive correlation with tumor mutational burden (TMB) and immune checkpoint blockade (ICB) expression and was closely correlated with multiple anticancer agent susceptibility.

Conclusion

This work revealed a close relationship between necroptosis and the prognosis and immune microenvironment of colon cancer. The NRS-score based on the 8-gene signature may be used to predict the sensitivity of immunotherapy and chemotherapy in colon cancer patients, and provides a foundation for future studies targeting necroptosis and its immune microenvironment.

Extracellular Matrix-Based Gene Expression Signature Defines Two Prognostic Subtypes of Hepatocellular Carcinoma With Different Immune Microenvironment Characteristics

Background: Accumulating evidence has suggested that the extracellular matrix (ECM) plays a vital role in the development and progression of cancer, and could be recognized as a biomarker of the response to immunotherapy. However, the effect of the ECM signature in hepatocellular carcinoma (HCC) is not well understood.

Methods: HCC patients derived from the TCGA-LIHC dataset were clustered according to the ECM signature. The differences in prognosis, functional enrichment, immune infiltration, and mutation characteristics between distinct molecular clusters were examined, and its predictive value on the sensitivities to chemotherapy and immunotherapy was further analyzed. Then, a prognostic model was built based on the ECM-related gene expression pattern.

Results: HCC patients were assigned into two molecular subtypes. Approximately 80% of HCC patients were classified into cluster A with poor prognosis, more frequent TP53 mutation, and lower response rate to immunotherapy. In contrast, patients in cluster B had better survival outcomes and higher infiltration levels of dendritic cells, macrophages, and regulatory T cells. The prognostic risk score model based on the expression profiles of six ECM-related genes (SPP1, ADAMTS5, MMP1, BSG, LAMA2, and CDH1) demonstrated a significant association with higher histologic grade and advanced TNM stage. Moreover, the prognostic risk score showed good performance in both the training dataset and validation dataset, as well as improved prognostic capacity compared with TNM stage.

Conclusions: We characterized two HCC subtypes with distinct clinical outcomes, immune infiltration, and mutation characteristics. A novel prognostic model based on the ECM signature was further developed, which may contribute to individualized prognostic prediction and aid in clinical decision-making.

Construction and Validation of Two Hepatocellular Carcinoma-Progression Prognostic Scores Based on Gene Set Variation Analysis

Background: Liver hepatocellular carcinoma (LIHC) remains a global health challenge with a low early diagnosis rate and high mortality. Therefore, finding new biomarkers for diagnosis and prognosis is still one of the current research priorities.

Methods: Based on the variation of gene expression patterns in different stages, the LIHC-development genes (LDGs) were identified by differential expression analysis. Then, prognosis-related LDGs were screened out to construct the LIHC-unfavorable gene set (LUGs) and LIHC-favorable gene set (LFGs). Gene set variation analysis (GSVA) was conducted to build prognostic scoring models based on the LUGs and LFGs. ROC curve analysis and univariate and multivariate Cox regression analysis were carried out to verify the diagnostic and prognostic utility of the two GSVA scores in two independent datasets. Additionally, the key LCGs were identified by the intersection analysis of the PPI network and univariate Cox regression and further evaluated their performance in expression level and prognosis prediction. Single-sample GSEA (ssGSEA) was performed to understand the correlation between the two GSVA enrichment scores and immune activity.

Result: With the development of LIHC, 83 LDGs were gradually upregulated and 247 LDGs were gradually downregulated. Combining with LIHC survival analysis, 31 LUGs and 32 LFGs were identified and used to establish the LIHC-unfavorable GSVA score (LUG score) and LIHC-favorable GSVA score (LFG score). ROC curve analysis and univariate/multivariate Cox regression analysis suggested the LUG score and LFG score could be great indicators for the early diagnosis and prognosis prediction. Four genes (ESR1, EHHADH, CYP3A4, and ACADL) were considered as the key LCGs and closely related to good prognosis. The frequency of TP53 mutation and copy number variation (CNV) were high in some LCGs. Low-LFG score patients have active metabolic activity and a more robust immune response. The high-LFG score patients characterized immune activation with the higher infiltration abundance of type I T helper cells, DC, eosinophils, and neutrophils, while the high-LUG score patients characterized immunosuppression with the higher infiltration abundance of type II T helper cells, TRegs, and iDC. The high- and low-LFG score groups differed significantly in immunotherapy response scores, immune checkpoints expression, and IC50 values of common drugs.

Conclusion: Overall, the LIHC-progression characteristic genes can be great diagnostic and prognostic signatures and the two GSVA score systems may become promising indices for guiding the tumor treatment of LIHC patients.

A Prognostic Model Using Immune-Related Genes for Colorectal Cancer

There is evidence suggesting that immune genes play pivotal roles in the development and progression of colorectal cancer (CRC). Colorectal carcinoma patient data from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) were randomly classified into a training set, a test set, and an external validation set. Differentially expressed gene (DEG) analyses, univariate Cox regression, and the least absolute shrinkage and selection operator (LASSO) were used to identify survival-associated immune genes and develop a prognosis model. Receiver operating characteristic (ROC) analysis and principal component analysis (PCA) were used to evaluate the discrimination of the risk models. The model genes predicted were verified using the Human Protein Atlas (HPA) databases, colorectal cell lines, and fresh CRC and adjacent tissues. To understand the relationship between IRGs and immune invasion and the TME, we analyzed the content of immune cells and scored the TME using CIBERSORT and ESTIMATE algorithms. Finally, we predicted the potential sensitive chemotherapeutic drugs in different risk score groups by the Genomics of Drug Sensitivity in Cancer (GDSC). A total of 491 IRGs were screened, and 14 IRGs were identified to be significantly related to overall survival (OS) and applied to construct an immune-related gene (IRG) prognostic signature (IRGSig) for CRC patients. Calibration plots showed that nomograms have powerful predictive ability. PCA and ROC analysis further verified the predictive value of this fourteen-gene prognostic model in three independent databases. Furthermore, we discovered that the tumor microenvironment changed significantly during the tumor development process, from early to middle to late stage, which may be an essential factor for tumor deterioration. Finally, we selected six commonly used chemotherapeutic drugs that have the potential to be useful in the treatment of CRC. Altogether, immune genes were used to construct a prognosis model for CRC patients, and a variety of methods were used to test the accuracy of this model. In addition, we explored the immune mechanisms of CRC through immune cell infiltration and TME in CRC. Furthermore, we assessed the therapeutic sensitivity of many commonly used chemotherapeutic medicines in individuals with varying risk factors. Finally, the immune risk model and immune mechanism of CRC were thoroughly investigated in this paper.

Pan-Cancer Analysis Identified CD93 as a Valuable Biomarker for Predicting Patient Prognosis and Immunotherapy Response

Background: The rapid development of immunotherapy has significantly improved patient outcomes in recent years. CD93, a novel biomarker expressed on vascular endothelial cells, is essential for tumor angiogenesis. Recent studies have shown that CD93 is closely related to immune cell infiltration and immunotherapy. However, its role in pan-cancer has not been reported.

Methods: The Cancer Genome Atlas (TCGA), Human Protein Atlas (HPA), cbioportal, Gene Expression Omnibus (GEO), Tumor Immune Estimation Resource (TIMER2.0), and the Tumor–Immune System Interactions and Drug Bank (TISIDB) databases were used to analyze CD93 in pan-cancers. R software was used for statistical analysis and mapping.

Results: There were significant differences in the expression of CD93 between tumor tissues and adjacent normal tissues in pan-cancer. The high expression of CD93 was associated with poor prognosis and high TNM stage in multiple tumor types. However, a high expression of CD93 was a protective factor in kidney renal clear cell carcinoma (KIRC). In addition, CD93 was closely related to immune cell infiltration in tumor tissues. Moreover, CD93 presented a robust correlation with immune modulators and immunotherapeutic markers [e.g., tumor mutation burden (TMB) and microsatellite instability (MSI)]. The results of gene set enrichment analysis (GSEA) showed that CD93 was correlated with tumor angiogenesis. Importantly, patients with a low expression of CD93 were more sensitive to immunotherapy in urothelial cancer.

Conclusion: CD93, which is involved in various immune responses, controls immune cell infiltration and impacts on the malignant properties of various cancer types. Therefore, CD93 has potential value to be biomarker for determining the prognosis and immune infiltration in multiple cancers.

Ferroptosis Mediation Patterns Reveal Novel Tool to Implicate Immunotherapy and Multi-Omics Characteristics in Bladder Cancer

Background: The regulatory role of ferroptosis in malignant tumours has been recently demonstrated. However, the potential roles of ferroptosis mediation patterns in bladder cancer remain elusive. Materials and Methods: The ferroptosis mediation patterns of 889 bladder cancer samples were comprehensively evaluated based on ferroptosis-related genes. The underlying correlations between these mediation patterns and multi-omic characteristics of bladder cancer were systematically analysed. The ferroptosis mediation patterns of individual samples were quantified by ferropscore using the principal component analysis algorithm. The typical ferroptosis-related genes with prognostic roles were further randomly validated using immunohistochemical staining, real-time polymerase chain reaction and western blotting. Results: Three different ferroptosis mediation patterns were identified. The abundance of infiltration of 23 immune cells was different among the three mediation patterns. The quantification of ferroptosis mediation patterns in individual samples served as a promising tool for predicting patient survival outcomes; immune cell infiltration abundance; tumour mutation burden; oncogenic mutation status and tumour grade, stage and molecular subtypes. Low ferropscore combined with high tumour mutation burden was associated with the best survival prognosis. Expressions of PD-L1 (p < 0.001), PD-1 (p = 0.002) and CTLA-4 (p = 0.003) were all significantly upregulated in the high ferropscore group. Low ferropscores also predicted good immunotherapy response for anti-CTLA4 strategy. The mRNA and protein levels of FADS2, a typical ferroptosis-related gene used in the study, were higher in bladder cancer cell lines than in controlled SV-HUC-1 cells. In addition, immunohistochemical staining revealed significantly higher expression levels of FADS2 in human bladder cancer tumour tissues than in normal tissues. Conclusion: This study identified three distinct ferroptosis mediation patterns in bladder cancer. Quantification of ferroptosis mediation patterns in individual samples may help to improve the understanding of multiomic characteristics and guide future immunotherapy responses to bladder cancer.

Immune Landscape Refines the Classification of Colorectal Cancer With Heterogeneous Prognosis, Tumor Microenvironment and Distinct Sensitivity to Frontline Therapies

The immune microenvironment has profound impacts on the initiation and progression of colorectal cancer (CRC). Therefore, the goal of this article is to identify two robust immune subtypes in CRC, further provide novel insights for the underlying mechanisms and clinical management. In this study, two CRC immune subtypes were identified using the consensus clustering of immune-related gene expression profiles in the meta-GEO dataset (n = 1,198), and their reproducibility was further verified in the TCGA-CRC dataset (n = 638). Subsequently, we characterized the immune escape mechanisms, gene alterations, and clinical features of two immune subtypes. Cluster 1 (C1) was defined as the “immune cold subtype” with immune cell depletion and deficiency, while cluster 2 (C2) was designed as the “immune hot subtype”, with abundant immune cell infiltration and matrix activation. We also underlined the potential immune escape mechanisms: lack of MHC molecules and defective tumor antigen presentation capacity in C1, increased immunosuppressive molecules in C2. The prognosis and sensitivity to 5-FU, Cisplatin and immunotherapy differed between two subtypes. According to the two immune subtypes, we developed a prognosis associated risk score (PARS) with the accurate performance for predicting the prognosis. Additionally, two nomograms for overall survival (OS) and disease-free survival (DFS) were further constructed to facilitate clinical management. Overall, our research provides new references and insights for understanding and refining the CRC.

A Novel Lactate Metabolism-Related Gene Signature for Predicting Clinical Outcome and Tumor Microenvironment in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is the main subtype of primary liver cancer with high malignancy and poor prognosis. Metabolic reprogramming is a hallmark of cancer and has great importance on the tumor microenvironment (TME). As an abundant metabolite, lactate plays a crucial role in cancer progression and the immunosuppressive TME. Nonetheless, the potential roles of lactate in HCC remain unclear. In this study, we downloaded transcriptomic data of HCC patients with corresponding clinical information from the TCGA and ICGC portals. The TCGA-HCC dataset used as the training cohort, while the ICGC-LIRI-JP dataset was served as an external validation cohort. Cox regression analysis and the LASSO regression model were combined to construct the lactate metabolism-related gene signature (LMRGS). Then, we assessed the clinical significance of LMRGS in HCC. Besides, enriched molecular functions, tumor mutation burden (TMB), infiltrating immune cells, and immune checkpoint were comprehensively analyzed in different LMRGS subgroups. In total, 66 differentially expressed lactate metabolism-related genes (LMRGs) were screened. The functions of LMRGs were mainly enriched in mitochondrial activity and metabolic processes. The LMRGS comprised of six key LMRGs (FKTN, PDSS1, PET117, PUS1, RARS1, and RNASEH1) had significant clinical value for independently predicting the prognosis of HCC patients. The overall survival and median survival of patients in the LMRGS-high group were significantly shorter than in the LMRGS-low group. In addition, there were differences in TMB between the two LMRGS subgroups. The probability of genetic mutations was higher in the LMRGS-high group. Most importantly, the LMRGS reflected the TME characteristics. In the LMRGS-high group, the immune microenvironment presented a suppressed state, accompanied by more inhibitory immune cell infiltration, including follicular helper T cells and regulatory T cells. Additionally, the expression of inhibitory checkpoint molecules was much higher in the LMRGS-high group. Our study suggested that the LMRGS was a robust biomarker to predict the clinical outcomes and evaluate the TME of patients with HCC.

LCK and CD3E Orchestrate the Tumor Microenvironment and Promote Immunotherapy Response and Survival of Muscle-Invasive Bladder Cancer Patients

Background: Studies have demonstrated the significance of multiple biomarkers for bladder cancer. Here, we attempt to present biomarkers potentially predictive of the prognosis and immunotherapy response of muscle-invasive bladder cancer (MIBC). Method: Immune and stromal scores were calculated for MIBC patients from The Cancer Genome Atlas (TCGA). Core differential expression genes (DEGs) with prognostic value were identified and validated using an independent dataset GSE31684. The clinical implications of prognostic genes and the inter-gene correlation were presented. The distribution of tumor-infiltrating immune cells (TICs), the correlation with tumor mutation burden (TMB), and the expression of eight immune checkpoint-relevant genes and CD39 were accordingly compared. Two bladder cancer cohorts (GSE176307 and IMvigor210) receiving immunotherapy were recruited to validate the prognostic value of LCK and CD3E for immunotherapy. Results: 361 MIBC samples from TCGA revealed a worse overall survival for higher stromal infiltration (p = 0.009) but a better overall survival for higher immune infiltration (p = 0.042). CD3E and LCK were independently validated by TCGA and GSE31684 to be prognostic for MIBC. CD3E was the most correlative gene of LCK, with a coefficient of r = 0.86 (p < 0.001). CD8⁺ T cells and macrophage M1 are more abundant in favor of a higher expression of CD3E and LCK in MIBC and across pan-cancers. Immune checkpoints like CTLA4, CD274 (PD-1), and PDCD1 (PD-L1) were highly expressed in high-CD3E and high-LCK groups for MIBC and also for pan-cancers, except for thymoma. LCK and CD3E had a moderate positive correlation with CD39 expression. Importantly, high-LCK and high-CD3E groups had a higher percentage of responders than the low-expression groups both in GSE176307 (LCK: 22.73vs. 13.64%, CD3E: 22.00 vs. 13.16%) and IMvigor210 cohorts (LCK: 28.19 vs. 17.45%, CD3E: 25.50 vs. 20.13%). Conclusion: CD3E and LCK were potential biomarkers of MIBC. CD3E and LCK were positively correlated with several regular immunotherapy biomarkers, which is supported by real-world outcomes from two immunotherapy cohorts.

CD8+ T Cell-Based Molecular Classification With Heterogeneous Immunogenomic Landscapes and Clinical Significance of Clear Cell Renal Cell Carcinoma

The tumor microenvironment (TME) exerts a high impact on tumor biology and immunotherapy. The heterogeneous phenotypes and the clinical significance of CD8⁺ T cells in TME have not been fully elucidated. Here, a comprehensive immunogenomic analysis based on multi-omics data was performed to investigate the clinical significance and tumor heterogeneity between CD8⁺ T cell-related molecular clusters. We identified two distinct molecular clusters of ccRCC (C1 and C2) in TCGA and validated in E-MTAB-1980 cohorts. The C1 cluster was characterized by unfavorable prognosis, increased expression levels of CD8⁺ T cell exhaustion markers, high immune infiltration levels as well as more immune escape mechanisms. The C2 cluster was featured by favorable prognosis, elevated expression levels of CD8⁺ T cell effector markers, low load of copy number loss and low frequency of 9p21.3 deletion. Moreover, the effect of molecular classifications on Nivolumab therapeutic efficacy in the CheckMate 025 cohort was examined, and the C2 cluster exhibited a better prognosis. Taken together, we determine two CD8⁺ T cell-related molecular clusters in ccRCC, and provide new insights for evaluating the functions of CD8⁺ T cells. Our molecular classification is a potential strategy for prognostic prediction and immunotherapeutic guidance for ccRCC patients.

Computational Recognition and Clinical Verification of TGF-β-Derived miRNA Signature With Potential Implications in Prognosis and Immunotherapy of Intrahepatic Cholangiocarcinoma

MicroRNAs (miRNAs) were recently implicated in modifying the transforming growth factor β (TGF-β) signaling in multiple cancers. However, TGF-β-derived miRNAs and their potential clinical significance remain largely unexplored in intrahepatic cholangiocarcinoma (ICC). In this study, we proposed an integrated framework that enables the identification of TGF-β-derived miRNAs in ICC (termed “TGFmitor”). A total of 36 TGF-β-derived miRNAs were identified, of which nine significantly correlated with overall survival (OS) and aberrantly expressed in ICC. According to these miRNAs, we discovered and validated a TGF-β associated miRNA signature (TAMIS) in GSE53870 (n =63) and TCGA-CHOL (n =32). To further confirm the clinical interpretation of TAMIS, another validation based on qRT-PCR results from 181 ICC tissues was performed. TAMIS was proven to be an independent risk indicator for both OS and relapse-free survival (RFS). TAMIS also displayed robust performance in three cohorts, with satisfactory AUCs and C-index. Besides, patients with low TAMIS were characterized by superior levels of CD8+ T cells infiltration and PD-L1 expression, while patients with high TAMIS possessed enhanced CMTM6 expression. Kaplan-Meier analysis suggested CMTM6 could further stratify TAMIS. The TAMIS^highCMTM6^high subtype had the worst prognosis and lowest levels of CD8A and PD-L1 expression relative to the other subtypes, indicating this subtype might behave as “super-cold” tumors. Notably, the improved discrimination was observed when CMTM6 was combined with TAMIS. Overall, our signature could serve as a powerful tool to help improve prognostic management and immunotherapies of ICC patients.

Immunogenomic Analyses of the Prognostic Predictive Model for Patients With Renal Cancer

Background

Renal cell carcinoma (RCC) is associated with poor prognostic outcomes. The current stratifying system does not predict prognostic outcomes and therapeutic benefits precisely for RCC patients. Here, we aim to construct an immune prognostic predictive model to assist clinician to predict RCC prognosis.

Methods

Herein, an immune prognostic signature was developed, and its predictive ability was confirmed in the kidney renal clear cell carcinoma (KIRC) cohorts based on The Cancer Genome Atlas (TCGA) dataset. Several immunogenomic analyses were conducted to investigate the correlations between immune risk scores and immune cell infiltrations, immune checkpoints, cancer genotypes, tumor mutational burden, and responses to chemotherapy and immunotherapy.

Results

The immune prognostic signature contained 14 immune-associated genes and was found to be an independent prognostic factor for KIRC. Furthermore, the immune risk score was established as a novel marker for predicting the overall survival outcomes for RCC. The risk score was correlated with some significant immunophenotypic factors, including T cell infiltration, antitumor immunity, antitumor response, oncogenic pathways, and immunotherapeutic and chemotherapeutic response.

Conclusions

The immune prognostic, predictive model can be effectively and efficiently used in the prediction of survival outcomes and immunotherapeutic responses of RCC patients.