Unsupervised Consensus Clustering Research Articles

A Gluconeogenesis-Related Genes Model for Predicting Prognosis, Tumor Microenvironment Infiltration, and Drug Sensitivity in Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is a prevalent malignancy within the digestive system, known for its poor prognosis. Gluconeogenesis, a critical metabolic pathway, is responsible for the synthesis of glucose in the normal liver. This study aimed to examine the role of gluconeogenesis-related genes (GRGs) in HCC and evaluate their impact on the tumor microenvironment infiltration and drug sensitivity in HCC. We retrieved gene expression and clinical pathological data of HCC from The Cancer Genome Atlas (TCGA) database. This dataset was utilized to develop a prognosis model. The data from The International Cancer Genome Consortium (ICGC) served as an independent validation cohort. A least absolute shrinkage and selection operator (LASSO) regression analysis was applied to a curated panel of GRGs to construct and validate the predictive model. Furthermore, unsupervised consensus clustering, based on the expression levels of GRGs, categorized HCC patients into distinct subgroups. A four-gene prognostic model, referred to as GRGs, has been successfully developed with high accuracy and stability for the prediction of HCC patient prognosis. This model enables the stratification of patients into high or low risk groups based on individual risk scores, revealing significant differences in immune infiltration patterns and anti-tumor drug responses. Unsupervised consensus clustering analysis delineated four distinct subgroups of patients, each characterized by a unique prognosis and tumor immune microenvironment (TIME). This study is the first to develop a prognostic model incorporating 4-GRGs that effectively predicts the prognosis, tumor microenvironment infiltration, and drug sensitivity in HCC patients. The model based on 4 GRGs may contribute to predict the prognosis, immunotherapy and chemotherapy response of HCC patients.

Prognosis of colorectal cancer, prognostic index of immunogenic cell death associated genes in response to immunotherapy, and potential therapeutic effects of ferroptosis inducers.

This study leverages bioinformatics and medical big data to integrate datasets from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA), providing a comprehensive overview of immunogenic cell death (ICD)-related gene expression in colorectal cancer (CRC). The research aims to elucidate the molecular pathways and gene networks associated with ICD in CRC, with a focus on the therapeutic potential of cell death inducers, including ferroptosis agents, and their implications for precision medicine. We conducted differential expression analysis and utilized advanced bioinformatic techniques to analyze ICD-related gene expression in CRC tissues. Unsupervised consensus clustering was applied to categorize CRC patients into distinct ICD-associated subtypes, followed by an in-depth immune microenvironment analysis and single-cell RNA sequencing to investigate immune responses and cell infiltration patterns. Experimental validation was performed to assess the impact of cell death inducers on ICD gene expression and their interaction with ferroptosis inducers in combination with other clinical drugs. Distinct ICD gene expression profiles were identified in CRC tissues, revealing molecular pathways and intricate gene networks. Unsupervised consensus clustering refined the CRC cohort into unique ICD-associated subtypes, each characterized by distinct clinical and immunological features. Immune microenvironment analysis and single-cell RNA sequencing revealed significant variations in immune responses and cell infiltration patterns across these subtypes. Experimental validation confirmed that cell death inducers directly affect ICD gene expression, highlighting their therapeutic potential. Additionally, combinatorial therapies with ferroptosis inducers and clinical drugs were shown to influence drug sensitivity and resistance in CRC. Our findings underscore the importance of ICD-related genes in CRC prognosis and therapeutic targeting. The study provides actionable insights into the efficacy of cell death-inducing therapies, particularly ferroptosis inducers, and their regulatory mechanisms in CRC. These discoveries support the development of precision medicine strategies targeting ICD genes and offer valuable guidance for translating these therapies into clinical practice, with the potential to enhance CRC treatment outcomes and patient survival.

Identification and Validation of a Necroptosis-Related Prognostic Model in Tumor Recurrence and Tumor Immune Microenvironment in Breast Cancer Management.

Breast cancer is the leading cause of cancer-related death in women. Necroptosis, a form of programmed necrotic cell death, occurs in many solid tumors, including breast cancer, and influences anti-tumor immunity. The role of necroptosis in managing breast cancer recurrence remains unclear. Gene expression profiles and clinical data of breast cancer patients were obtained from the GEO (GSE20685, GSE21653, GSE25055) and TCGA databases. Data analysis and visualization were performed using R. Unsupervised Consensus Clustering and LASSO-COX regression stratified breast cancer patients. GO, KEGG, GSVA, ESTIMATE, and ROC analyses were used to investigate necroptotic signatures. In vitro and in vivo experiments validated necroptosis's role in breast cancer immunity. The potential function of necroptotic signature in immunity was first indicated with GO analysis in BRCA cohort. Next, two prognostic models based on the necroptotic profiles both suggested a link between low-risk group with a particular necroptotic immune signature. And a variety of immune cells and immune pathways were shown to be positively associated with a patient's risk score. As an altered immune checkpoint pattern was observed after regulating necroptotic genes, where TIM-3 and LAGLS9 elevated significantly in low-risk group, further validation in vitro and in vivo demonstrated that manipulating a subset of necroptotic gene set could sensitize tumor response to the co-blockade immunotherapy of anti-TIM-3 and anti-PD-1. We demonstrated two strategies to stratify breast cancer patients based on their necroptotic profiles and showed that necroptotic signature could assign patients with different tumor immune microenvironment patterns and different recurrence-related prognosis. A subset of necroptotic gene set, composed of TLR3, RIPK3, NLRP3, CASP1, ALDH2 and EZH2, was identified as a biomarker set for predicting immunotherapy-response and recurrence-related prognosis. Targeting necroptosis could helpfacilitate the development of novel breast cancer treatments and tailor personalized medical treatment.

Comprehensive Analysis of Ligand-receptor Interactions in Colon Adenocarcinoma to Identify of Tumor Microenvironment Oxidative Stress and Prognosis Model.

Single-cell technology enables a deep study on the mechanism of cancers. This work delineated the function of ligand-receptor interaction in colon adenocarcinoma (COAD), and developed a LR pairs-based prognostic model. For identifying important LR pairs, Single-cell RNA sequencing data of COAD was included. Unsupervised consensus clustering constructed molecular subtypes. LASSO established a prognostic model. Infiltration of 22 immune cells was evaluated by Cibersort. Enrichment score of oxidative stress related pathways was determined by SsGSEA in each patient. Forty-seven LR pairs were closely associated with the prognosis of COAD. Three molecular subtypes were differentiated according to 47 LR pairs, which displayed differential clinical features and molecular features. There were significant differences in immune T cell lytic activity among different subtypes. In clust1 with poor prognosis, significantly enriched oncogenic pathways were found, especially epithelial-mesenchymal transition (EMT). Additionally, it has been found that clust3 had significantly higher immune infiltration. A prognostic model containing eight LR pairs (PDGFB-PDGFRA, FLT4-VEGFC, CSF1R-CSF1, DLL1-NOTCH4, PDGFB-LRP1, DLL1- NOTCH3, FLT4-PDGFC, and NRP2-PGF) was established, which could effectively divide samples into low-risk and high-risk groups. Significantly higher oxidative stress was found among high-risk patients. This study integrated expression data and single-cell data for demonstrating the effectiveness of LR pairs in establishing the prognostic model and constructing molecular subtypes. Prognostic LR pairs may contribute to tumorigenesis and progression in COAD. The prognostic model was the potential for predicting prognosis and guiding immunotherapy for COAD patients.

Roles of RNA m5C modification patterns in prognosis and tumor microenvironment infiltration of diffuse large B-cell lymphoma.

Genetic and epigenetic aberrations display an essential role in the initiation and progression of diffuse large B-cell lymphoma (DLBCL). 5-methylcytosine (m5C), a common RNA modification, regulates various cellular processes and contributes to tumorigenesis and cancer progression. However, m5C alterations in DLBCL remain unclear. Our research constructed an m5C prognostic model utilizing GEO data sets, which can efficiently predict the prognosis of patients with DLBCL, and verified the m5C prognostic model genes by immunohistochemistry analysis. This model was constructed using unsupervised consensus clustering analyses, Least Absolute Shrinkage and Selection Operator (LASSO), and multivariate Cox regression analyses. Based on the expression of m5C genes in the model, patients with DLBCL could be effectively divided into groups with significant survival time differences. The m5C risk-score signature demonstrated a highly significant independent prognostic value. Results from tumor microenvironment analyses revealed that m5C genes altered the infiltration of eosinophils, Tregs, and M2 macrophages. Additionally, they regulated T cell activation by modulating the expression of CTLA4, PDL1, B2M, CD8A, ICOS, and other relevant immune checkpoint expressions. In conclusion, our study presents a robust m5C prognostic model that effectively predicts prognosis in DLBCL. This model may offer a new approach for prognostic stratification and potential therapeutic interventions for patients with DLBCL.

Machine learning identifies ferroptosis-related gene ANXA2 as potential diagnostic biomarkers for NAFLD.

Non-alcoholic fatty liver disease (NAFLD), a major cause of chronic liver disease, still lacks effective therapeutic targets today. Ferroptosis, a type of cell death characterized by lipid peroxidation, has been linked to NAFLD in certain preclinical trials, yet the exact molecular mechanism remains unclear. Thus, we analyzed the relationship between ferroptosis genes and NAFLD using high-throughput data. We utilized a total of 282 samples from five datasets, including two mouse ones, one human one, one single nucleus dataset and one single cell dataset from Gene Expression Omnibus (GEO), as the data basis of our study. To filter robust treatment targets, we employed four machine learning methods (LASSO, SVM, RF and Boruta). In addition, we used an unsupervised consensus clustering algorithm to establish a typing scheme for NAFLD based on the expression of ferroptosis related genes (FRGs). Our study is also the first to investigate the dynamics of FRGs throughout the disease process by time series analysis. Finally, we validated the relationship between core gene and ferroptosis by in vitro experiments on HepG2 cells. We discovered ANXA2 as a central focus in NAFLD and indicated its potential to boost ferroptosis in HepG2 cells. Additionally, based on the results obtained from time series analysis, ANXA2 was observed to significantly define the disease course of NAFLD. Our results demonstrate that implementing a ferroptosis-based staging method may hold promise for the diagnosis and treatment of NAFLD. Our findings suggest that ANXA2 may be a useful biomarker for the diagnosis and characterization of NAFLD.

Comprehensive analysis of mitophagy-related genes in NSCLC diagnosis and immune scenery: based on bulk and single-cell RNA sequencing data.

Lung cancer is the main cause of cancer-related deaths, and non-small cell lung cancer (NSCLC) is the most common type. Understanding the potential mechanisms, prognosis, and treatment aspects of NSCLC is essential. This study systematically analyzed the correlation between mitophagy and NSCLC. Six mitophagy-related feature genes (SRC, UBB, PINK1, FUNDC1, MAP1LC3B, and CSNK2A1) were selected through machine learning and used to construct a diagnostic model for NSCLC. These feature genes are closely associated with the occurrence and development of NSCLC. Additionally, NSCLC was divided into two subtypes using unsupervised consensus clustering, and their differences in clinical characteristics, immune infiltration, and immunotherapy were systematically analyzed. Furthermore, the interaction between mitophagy-related genes (MRGs) and immune cells was analyzed using single-cell sequencing data. The findings of this study will contribute to the development of potential diagnostic biomarkers for NSCLC and the advancement of personalized treatment strategies.

A methylation-related signature for predicting prognosis and sensitivity to first-line therapies in gastric cancer.

Methylation modification patterns play a crucial role in human cancer progression, especially in gastrointestinal cancers. We aimed to use methylation regulators to classify patients with gastric adenocarcinoma and build a model to predict prognosis, promoting the application of precision medicine. We obtained RNA sequencing data and clinical data from The Cancer Genome Atlas (TCGA) database (n=335) and Gene Expression Omnibus (GEO) database (n=865). Unsupervised consensus clustering was used to identify subtypes of gastric adenocarcinoma. We performed functional enrichment analysis, immune infiltration analysis, drug sensitivity analysis, and molecular feature analysis to determine the clinical application for different subtypes. The univariate Cox regression analysis and the LASSO regression analysis were subsequently used to identify prognosis-related methylation regulators and construct a risk model. Through unsupervised consensus clustering, patients were divided into two subtypes (cluster A and cluster B) with different clinical outcomes. Cluster B included patients with a better prognosis outcome and who were more likely to respond to immunotherapy. We then successfully built a predictive model and found five methylation-related genes (CHAF1A, CPNE8, PHLDA3, SPARC, and EHF) potentially significant to the prognosis of patients. The 1-, 3-, and 5-year areas under the curve of the risk model were 0.712, 0.696, and 0.759, respectively. The risk score was an independent prognostic factor and had the highest concordance index among common clinical indicators. Meanwhile, the tumor microenvironment, sensitivity of chemotherapeutic drugs, molecular features, and oncogenic dedifferentiation differed significantly across the risk groups and subtypes. We classified patients with gastric adenocarcinoma based on methylation regulators, which has positive implications for first-line clinical treatment. The prognostic model could predict the prognosis of patients and help to promote the development of precision medicine.

Discovery of distinct cancer cachexia phenotypes using an unsupervised machine-learning algorithm

ObjectivesCancer cachexia is a debilitating condition with widespread negative effects. The heterogeneity of clinical features within patients with cancer cachexia is unclear. The identification and prognostic analysis of diverse phenotypes of cancer cachexia may help develop individualized interventions to improve outcomes for vulnerable populations. The aim of this study was to show that the machine learning–based cancer cachexia classification model generalized well on the external validation cohort. MethodsThis was a nationwide multicenter observational study conducted from October 2012 to April 2021 in China. Unsupervised consensus clustering analysis was applied based on demographic, anthropometric, nutritional, oncological, and quality-of-life data. Key characteristics of each cluster were identified using the standardized mean difference. We used logistic and Cox regression analysis to evaluate 1-, 3-, 5-y, and overall mortality. ResultsA consensus clustering algorithm was performed for 4329 patients with cancer cachexia in the discovery cohort, and four clusters with distinct phenotypes were uncovered. From clusters 1 to 4, the clinical characteristics of patients showed a transition from almost unimpaired to mildly, moderately, and severely impaired. Consistently, an increase in mortality from clusters 1 to 4 was observed. The overall mortality rate was 32%, 40%, 54%, and 68%, and the median overall survival time was 21.9, 18, 16.7, and 13.6 mo for patients in clusters 1 to 4, respectively. Our machine learning-based model performed better in predicting mortality than the traditional model. External validation confirmed the above results. ConclusionsMachine learning is valuable for phenotype classifications of patients with cancer cachexia. Detection of clinically distinct clusters among cachexic patients assists in scheduling personalized treatment strategies and in patient selection for clinical trials.

A comprehensive cuproptosis score and associated gene signatures reveal prognostic and immunological features of idiopathic pulmonary fibrosis.

Cuproptosis, the most recently identified and regulated cell death, depends on copper ions in vivo. Copper regulates the pathogenesis of Idiopathic pulmonary fibrosis (IPF), but the mechanism of action underlying cuproptosis in IPF remains unclear. We identified three cuproptosis patterns based on ten cuproptosis-related genes using unsupervised consensus clustering. We quantified these patterns using a PCA algorithm to construct a cuproptosis score. ssGSEA and the Cibersort algorithm assessed the immune profile of IPF patients. GSEA and GSVA were used to analyze the functional differences in different molecular patterns. Drug susceptibility prediction based on cuproptosis scores and meaningful gene markers was eventually screened in combination with external public data sets,in vitro experiments and our cases. Of the three types of cuproptosis-related clusters identified in the study, patients in the clusterA, geneclusterB, and score-high groups showed improved prognoses. Moreover, each cluster exhibited differential immune characteristics, with the subtype showing a poorer prognosis associated with an immune overreaction. Cuproptosis score can be an independent risk factor for predicting the prognosis of IPF patients. GSEA showed a significant functional correlation between the score and cuproptosis. The genes AKAP9, ANK3, C6orf106, LYRM7, and MBNL1, were identified as prognostic-related signatures in IPF patients. The functional role of immune regulation in IPF was further explored by correlating essential genes with immune factors. Also, the nomogram constructed by cumulative information from gene markers and cuproptosis score showed reliable clinical application. Cuproptosis patterns differ significantly in the prognosis and immune characteristics of IPF patients. The cuproptosis score and five gene signatures can provide a reliable reference in the prognosis and diagnosis of IPF.

Establishment of a new molecular subtyping and prognostic signature with m6A/m5C/m1A/m7G regulatory genes for hepatocellular carcinoma

BackgroundRNA modification, including m6A, m5C, m1A, and m7G, participated in tumor progress. Therefore, the purpose of the present study was to explore the role of m6A/m5C/m1A/m7G regulatory genes in the prognosis and tumor microenvironment (TME) for hepatocellular carcinoma (HCC). Methods71 m6A/m5C/m1A/m7G regulatory genes expression for HCC was detected, differentially expressed genes were screened, and molecular forms were classified by unsupervised consensus clustering. Cox regression and the Least Absolute Shrinkage and Selection Operator (LASSO) analysis were applied to establish a prognostic signature. Time-dependent receiver operating characteristic (ROC) curves were evaluated for clinical effectiveness and accuracy of the prognostic hazard model. In cluster subtypes and risk models, the differences in prognosis, immune cell infiltration, immune checkpoint, immunotherapy, and drug sensitivity between different subtypes were evaluated. ResultsHCC patients were classified into two clusters (cluster 1 and cluster 2) according to the expression of 71 m6A/m5C/m1A/m7G regulatory genes. Cluster 1 had a poor prognosis and different immune cell infiltration. Cluster 1 had higher immune checkpoint expression and TIDE score than cluster 2. Subsequently, we construct a five-gene prognostic model of m6A/m5C/m1A/m7G regulatory genes (YTHDF2, YTHDF1,YBX1, TRMT61A, TRMT10C). The Kaplan-Meier and ROC curve analysis showed that the prognostic signature exhibited good predictability. The risk score was considered an independent poor prognostic index. The high-risk group had higher immune checkpoint expression and higher TIDE scores. 5-Fluorouracil, docetaxel, doxorubicin, etoposide, gemcitabine, paclitaxel, sorafenib, and vinblastine were more suitable for high-risk patients. ECM receptor interaction, cell cycle, and Leishmania infection were enriched in the high-risk group. ConclusionThe clustering subgroups and prognostic model of m6A/m5C/m1A/m7G regulatory genes were linked with bad prognosis and TME for HCC, and had the potential to be a novel tool to evaluate the outcomes of HCC patients.

Molecular classification of human papilloma virus-negative head and neck squamous cell carcinomas: Cell cycle-based classifier and prognostic signature.

The molecular classification of human papillomavirus (HPV)-negative head and neck squamous cell carcinomas (HNSCCs) remains questionable. Differentially expressed genes were detected between tumor and normal tissues and GSEA showed they are associated with cell cycle pathways. This study aimed to classify HPV-negative HNSCCs based on cell cycle-related genes. The established gene pattern was correlated with tumor progression, clinical prognosis, and drug treatment efficacy. Biological analysis was performed using HNSCC patient sample data obtained from the Cancer Genome Atlas (TCGA), Clinical Proteomic Tumor Analysis Consortium (CPTAC), and Gene Expression Omnibus (GEO) databases. All samples included in this study contained survival information. RNA sequencing data from 740 samples were used for the analysis. Previously characterized cell cycle-related genes were included for unsupervised consensus clustering. Two subtypes of HPV-negative HNSCCs (C1, C2) were identified. Subtype C1 displayed low cell cycle activity, 'hot' tumor microenvironment (TME), earlier N stage, lower pathological grade, better prognosis, and higher response rate to the immunotherapy and targeted therapy. Subtype C2 was associated with higher cell cycle activity, 'cold' TME, later N stage, higher pathological grade, worse prognosis, and lower response rate to the treatment. According to the nearest template prediction method, classification rules were established and verified. Our work explored the molecular mechanism of HPV-negative HNSCCs in the view of cell cycle and might provide new sights for personalized anti-cancer treatment.

Identification of Molecular Subtypes and Prognostic Characteristics of Adrenocortical Carcinoma Based on Unsupervised Clustering.

Adrenocortical carcinoma (ACC) is a rare endocrine malignancy with a poor prognosis. Increasing evidence highlights the significant role of immune-related genes (IRGs) in ACC progression and immunotherapy, but the research is still limited. Based on the Cancer Genome Atlas (TCGA) database, immune-related molecular subtypes were identified by unsupervised consensus clustering. Univariate Cox analysis and Least Absolute Shrinkage and Selection Operator (LASSO) regression were employed to further establish immune-related gene signatures (IRGS). An evaluation of immune cell infiltration, biological function, tumor mutation burden (TMB), predicted immunotherapy response, and drug sensitivity in ACC patients was conducted to elucidate the applicative efficacy of IRGS in precision therapy. ACC patients were divided into two molecular subtypes through consistent clustering. Furthermore, the 3-gene signature (including PRKCA, LTBP1, and BIRC5) based on two molecular subtypes demonstrated consistent prognostic efficacy across the TCGA and GEO datasets and emerged as an independent prognostic factor. The low-risk group exhibited heightened immune cell infiltration, TMB, and immune checkpoint inhibitors (ICIs), associated with a favorable prognosis. Pathways associated with drug metabolism, hormone regulation, and metabolism were activated in the low-risk group. In conclusion, our findings suggest IRGS can be used as an independent prognostic biomarker, providing a foundation for shaping future ACC immunotherapy strategies.

RNA modification regulator DDC in endometrial cancer affects the tumor microenvironment and patient prognosis

Uterine corpus endometrial carcinoma (UCEC) is infiltrated by immune cells, which are involved in the growth and proliferation of malignant tumors and resistance to immunotherapy. This study suggested that RNA modification regulators played an important role in the development and prognosis of UCEC. Many studies confirmed that RNA modification played an essential role in tumor immune regulation, and abnormal RNA modification contributed to tumorigenesis and cancer progression. Based on the RNA modification regulatory factors, the UCEC samples from TCGA (The Cancer Genome Atlas) were classified into two clusters, namely Cluster A and Cluster B, using unsupervised consensus clustering. We obtained DEG (differentially expressed genes) between the two clusters, and constructed a risk model of RNA modification-related genes using DEGs. Cluster A had lower RNA modification regulatory factors, richer immune cell infiltration, and better prognosis. The differentially expressed genes between the two clusters were obtained, and these genes were used for modeling. This model divided patients with UCEC into two groups. The low-risk group had better immune infiltration, and the ROC (receiver operating characteristic) curve showed that this model had good predictive efficacy. The low-risk group had a better response to immunotherapy by immune checkpoint prediction. We obtained the key gene l-dopa decarboxylase (DDC) through the intersection of LASSO model genes and GEO dataset GSE17025. We evaluated the potential biological functions of DDC. The differences in the expression of DDC were verified by immunohistochemistry. We evaluated the relationship between DDC and immune cell infiltration and verified this difference using immunofluorescence. Cluster A with low expression of RNA modification regulators has better prognosis and richer immune cell infiltration, therefore, we believed that RNA modification regulators in UCEC were closely related to the tumor microenvironment. Also, the risk score could well predict the prognosis of patients and guide immunotherapy, which might benefit patients with UCEC.

Computational identification and clinical validation of a novel risk signature based on coagulation-related lncRNAs for predicting prognosis, immunotherapy response, and chemosensitivity in colorectal cancer patients.

Coagulation is critically involved in the tumor microenvironment, cancer progression, and prognosis assessment. Nevertheless, the roles of coagulation-related long noncoding RNAs (CRLs) in colorectal cancer (CRC) remain unclear. In this study, an integrated computational framework was constructed to develop a novel coagulation-related lncRNA signature (CRLncSig) to stratify the prognosis of CRC patients, predict response to immunotherapy and chemotherapy in CRC, and explore the potential molecular mechanism. CRC samples from The Cancer Genome Atlas (TCGA) were used as the training set, while the substantial bulk or single-cell RNA transcriptomics from Gene Expression Omnibus (GEO) datasets and real-time quantitative PCR (RT-qPCR) data from CRC cell lines and paired frozen tissues were used for validation. We performed unsupervised consensus clustering of CRLs to classify patients into distinct molecular subtypes. We then used stepwise regression to establish the CRLncSig risk model, which stratified patients into high- and low-risk groups. Subsequently, diversified bioinformatics algorithms were used to explore prognosis, biological pathway alteration, immune microenvironment, immunotherapy response, and drug sensitivity across patient subgroups. In addition, weighted gene coexpression network analysis was used to construct an lncRNA-miRNA-mRNA competitive endogenous network. Expression levels of CRLncSig, immune checkpoints, and immunosuppressors were determined using RT-qPCR. We identified two coagulation subclusters and constructed a risk score model using CRLncSig in CRC, where the patients in cluster 2 and the low-risk group had a better prognosis. The cluster and CRLncSig were confirmed as the independent risk factors, and a CRLncSig-based nomogram exhibited a robust prognostic performance. Notably, the cluster and CRLncSig were identified as the indicators of immune cell infiltration, immunoreactivity phenotype, and immunotherapy efficiency. In addition, we identified a new endogenous network of competing CRLs with microRNA/mRNA, which will provide a foundation for future mechanistic studies of CRLs in the malignant progression of CRC. Moreover, CRLncSig strongly correlated with drug susceptibility. We developed a reliable CRLncSig to predict the prognosis, immune landscape, immunotherapy response, and drug sensitivity in patients with CRC, which might facilitate optimizing risk stratification, guiding the applications of immunotherapy, and individualized treatments for CRC.

Diverse regulated cell death modes predict the immune microenvironment and drug sensitivity in lung adenocarcinoma.

Integrated action modes of regulated cell death (RCD)in lung adenocarcinoma (LUAD)have not been comprehensively dissected. Here, we adopted 15 RCD modes, including 1350 related genes, and established RCD signature scores. We found that LUAD patients with high RCD scores had a significantly worse prognosis in all four different cohorts (TCGA, KM-plotter, GSE31210, and GSE30219). Our nomogram established based on the RCD score and clinical characteristics performed well in both the discovery and validation sets. There was a close correlation between the RCD scores and LUAD molecular subtypes identified by unsupervised consensus clustering. Furthermore, we profiled the tumor microenvironment via deconvolution and found significant differences in immune activity, transcription factor activity and molecular pathway enrichment between the RCD-high and RCD-low groups. More importantly, we revealed that the regulation of antigen presentation is the crucial mechanism underlying RCD. In addition, higher RCD scores predict poorer sensitivity to multiple therapeutic drugs, which indicates that RCD scores may serve as a promising predictor of chemotherapy and immunotherapy outcomes. In summary, this work is the first to reveal the internal links between RCD modes, LUAD, and cancer immunity and highlights the necessity of RCD scores in personalizing treatment plans.

Comprehensive Analysis of METTLs (METTL1/13/18/21A/23/25/2A/2B/5/6/9) and Associated mRNA Risk Signature in Hepatocellular Carcinoma.

Currently, 80%-90% of liver cancers are hepatocellular carcinomas (HCC). HCC patients develop insidiously and have an inferior prognosis. The methyltransferase-like (METTL) family principal members are strongly associated with epigenetic and tumor progression. The present study mainly analyzed the value of METTLs (METTL1/13/18/21A/23/25/2A/2B/5/6/9) and associated mRNA risk signature for HCC. METTLs expression is upregulated in HCC and is a poor prognostic factor in HCC. METTLs were upregulated in patients older than 60 and associated with grade. Except for METTL25, the remaining 10 genes were associated with the HCC stage, invasion depth (T). In addition, METTLs showed an overall alteration rate of 50%. Except for METTL13/2A/25/9, the expression of the other seven genes was significantly associated with overall survival, disease-specific survival, and progression-free survival. Multivariate studies have shown that METTL21A/6 can be an independent prognostic marker in HCC. A total of 664 mRNAs were selected based on Pearson correlation coefficient (R > 0.5), unsupervised consensus clustering, weighted coexpression network analysis, and univariate Cox analysis. These mRNAs were significantly associated with METTLs and were poor prognostic factors in HCC patients. The least absolute shrinkage and selection operator (lasso) was used to construct the best METTLs associated with mRNA risk signature. The mRNA risk signature was significantly associated with age, stage, and t grade. The mRNA high-risk group had higher TP53 and RB1 mutations. This study constructed a nomogram with the mRNA risk profile and clinicopathological features, which could better predict the OS of individuals with HCC. We also analyzed associations between METTLs and mRNA risk signatures in epithelial-mesenchymal transition, immune checkpoints, immune cell infiltration, tumor mutational burden, microsatellite instability, cancer stem cells, tumor pathways, and drug sensitivity. In addition, this study constructed a protein interaction network network including METTLs and mRNA risk signature genes related to tumor microenvironment remodeling based on single-cell sequencing. In conclusion, this study provides a theoretical basis for the mechanism, biomarker screening, and treatment of HCC.

Dissecting tumor antigens and immune subtypes for mRNA vaccine development in breast cancer

PurposesCancer mRNA vaccines are a promising strategy and a hot topic in cancer immunotherapy. However, mRNA vaccines for breast cancer (BRCA) remain undeveloped. This study aimed to identify potential tumor antigens for mRNA vaccine development and a population with BRCA suitable for vaccination.MethodsGene expression profiles and the clinical information of the TCGA-BRCA (the Cancer Genome Atlas Breast Cancer) and METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) cohorts were downloaded from the TCGA and cBioPortal databases, respectively. cBioPortal was used to identify mutant genes. DEG (differentially expressed gene) identification and survival analysis were performed with the GEPIA2 tool. ssGSEA (single-sample gene set enrichment analysis) was applied to estimate abundances of 28 immune cells for each sample. An unsupervised consensus clustering algorithm was used to identify ISs (immune subtypes). A graph learning-based dimensionality reduction analysis algorithm was utilized to construct an immune landscape. WGCNA (weighted correlation network analysis) was performed to identify immune gene modules.ResultsFour potential tumor antigens, i.e., SLC7A5, CHPF, CCNE1, and CENPW, associated with poor prognosis and APCs (antigen-presenting cells) among overexpressed and mutated genes were identified in BRCA. Two ISs (IS1-2) characterized by distinct clinical, immune cell infiltration, and molecular features were observed in both the TCGA-BRCA and METABRIC cohorts. BRCA patients with IS2 tumors related to poor prognosis had an immune "hot" phenotype, while those patients with IS1 tumors related to superior prognosis had an immune "cold" phenotype. Distinct IS tumors were observed in different ICD (immunogenic cell death modulator) and ICP (immune checkpoint) expression profiles. The immune landscape showed an immune distribution in BRCA patients. Additionally, we identified 2 immune gene modules with different biological functions.ConclusionsSLC7A5, CHPF, CCNE1, and CENPW are the potential tumor antigens for mRNA vaccine development with BRCA. Patients with IS2 tumors are a suitable population for mRNA vaccination. This study provides a new insight into mRNA vaccine development, population selection for vaccination, and prognosis prediction.

Identification of prognostic immune subtypes of lung squamous cell carcinoma by unsupervised consistent clustering.

We performed UCC on the expression data of lung squamous cell carcinoma tumor samples to identify the classification of lung squamous cell carcinoma (LUSC) tumor samples, and calculated the levels of different classified immune cells by single-sample gene enrichment analysis (ssGSEA) to obtain a set of immune-related subtype gene tags, which can be used for subtype classification of lung squamous cell carcinoma. TCGA-LUSC and GSE30219 data of lung squamous cell carcinoma were obtained from TCGA and GEO databases. Prognostic-associated subtypes were identified by unsupervised consensus clustering (UCC). Using ssGSEA analysis to calculate the level of immune cells of different subtypes, obtain the connection between subtypes and immunity, identify the gene signatures recognized by subtypes, and verify this group of gene signatures through GSE30219. We effectively identified 2 subtypes that were significantly associated with prognostic survival by UCC, and calculated according to ssGSEA, the 2 subtypes were significantly different at the level of immune cells, followed by introducing a This weighted thinking computes a set of gene signatures that are significantly associated with subtype 1. During validation, this set of gene signatures could efficiently and robustly identify distinct prognostic immune subtypes, demonstrated the validity of this set of gene signatures, as well as 2 subtypes of lung squamous cell carcinoma. We used lung squamous cell carcinoma data from public databases and identified 2 prognostic immunosubtypes of lung squamous cell carcinoma and a set of gene tags that can be used to classify immune subtypes of lung squamous cell carcinoma, which may provide effective evidence for accurate clinical treatment of lung squamous cell carcinoma.

Machine learning-identified stemness features and constructed stemness-related subtype with prognosis, chemotherapy, and immunotherapy responses for non-small cell lung cancer patients

AimThis study aimed to explore a novel subtype classification method based on the stemness characteristics of patients with non-small cell lung cancer (NSCLC).MethodsBased on the Cancer Genome Atlas database to calculate the stemness index (mRNAsi) of NSCLC patients, an unsupervised consensus clustering method was used to classify patients into two subtypes and analyze the survival differences, somatic mutational load, copy number variation, and immune characteristics differences between them. Subsequently, four machine learning methods were used to construct and validate a stemness subtype classification model, and cell function experiments were performed to verify the effect of the signature gene ARTN on NSCLC.ResultsPatients with Stemness Subtype I had better PFS and a higher somatic mutational burden and copy number alteration than patients with Stemness Subtype II. In addition, the two stemness subtypes have different patterns of tumor immune microenvironment. The immune score and stromal score and overall score of Stemness Subtype II were higher than those of Stemness Subtype I, suggesting a relatively small benefit to immune checkpoints. Four machine learning methods constructed and validated classification model for stemness subtypes and obtained multiple logistic regression equations for 22 characteristic genes. The results of cell function experiments showed that ARTN can promote the proliferation, invasion, and migration of NSCLC and is closely related to cancer stem cell properties.ConclusionThis new classification method based on stemness characteristics can effectively distinguish patients' characteristics and thus provide possible directions for the selection and optimization of clinical treatment plans.