Consensus Clustering Analysis Research Articles

Subgrouping testicular germ cell tumors based on immunotherapy and chemotherapy associated lncRNAs

Testicular germ cell tumors (TGCT) are the most common reproductive system malignancies in men aged 15–44 years, accounting for 95 % of all testicular tumors. Our previous studies have been shown that long non-coding RNAs (lncRNAs), such as LINC00313, TTTY14 and RFPL3S, were associated with development of TGCT. Subgrouping TGCT according to differential expressed lncRNAs and immunological characteristics is helpful to comprehensively describe the characteristics of TGCT and implement precise treatment. In this study, the TGCT transcriptome data in The Cancer Genome Atlas Program (TCGA) database was used to perform consensus clustering analysis to construct a prognostic model for TGCT. TGCT was divided into 3 subtypes C1, C2, and C3 based on the differentially expressed lncRNAs. C1 subtype was sensitive to chemotherapy drugs, while the C2 subtype was not sensitive to chemotherapy drugs, and C3 subtype may benefit from immunotherapy. We defined the C1 subtype as epidermal progression subtype, the C2 subtype as mesenchymal progression subtype, and the C3 subtype as T cell activation subtype. Subgrouping based on differentially expressed genes (DEGs) and immunological characteristics is helpful for the precise treatment of TGCT.

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.

Comprehensive analysis of disulfidptosis-related lncRNA features for prognosis and immune landscape prediction in colorectal cancer.

Disulfidptosis is a novel mechanism underlying actin-cytoskeleton-associated cell death, but its function in colorectal cancer (CRC) is still elusive. In this study, we investigated the potential role of Disulfidptosis-Related Long Non-Coding RNAs (DRLs) as prognostic indicators in CRC. Through transcriptome data from TCGA CRC cases, we identified 44 prognosis-correlated DRLs by Univariate Cox Regression Analysis and observed a differential expression pattern of these DRLs between CRC and normal tissues. Consensus clustering analysis based on DRL expression led to subgroup classification of CRC patients with distinct molecular fingerprints, accompanied by a superior survival outcome in cluster 2. We are encouraged to develop a score model incorporating 12 key DRLs to predict patient outcomes. Notably, this model displayed more reliable accuracy than other predictive indicators since DRLs are intimately related to tumor immune cell infiltration, suggesting a considerable potential of our DRL-score model for tumor therapy. Our data offered a valuable insight into the prognostic significance of DRLs in CRC and broke a new avenue for tumor prognosis prediction.

Neutrophil infiltration associated genes on the prognosis and tumor immune microenvironment of lung adenocarcinoma.

The neutrophils exhibit both anti-tumor and pro-tumor effects in cancers. The correlation between neutrophils and tumor development in lung adenocarcinoma (LUAD) is still uncertain, possibly due to a lack of specific neutrophil infiltration evaluation methods. In this study, we identified 30 hub genes that were significantly associated with neutrophil infiltration in LUAD through data mining, survival analysis, and multiple tumor-infiltrating immune cells (TICs) analysis, including TIMER, CIBERSORT, QUANTISEQ, XCELL, and MCPCOUNTER. Consensus clustering analysis showed that these 30 hub genes were correlated with clinical features in LUAD. We further developed a neutrophil scoring system based on these hub genes. The neutrophil score was significantly correlated with prognosis and tumor immune microenvironment (TIME) in LUAD. It was also positively associated with PD-L1 expression and negatively associated with tumor mutational burden (TMB). When combined with the neutrophil score, the predictive capacity of PD-L1 and TMB for prognosis was significantly improved. Thus, the 30 hub genes might play an essential role in the interaction of neutrophils and LUAD, and the neutrophil scoring system might effectually assess the infiltration of neutrophils. Furthermore, we verified the expression of these 30 genes in the LUAD tumor tissues collected from our department. We further found that overexpressed TNFAIP6 and TLR6 and downregulated P2RY13, SCARF1, DPEP2, PRAM1, CYP27A1, CFP, GPX3, and NCF1 in LUAD tissue might be potentially associated with neutrophils pro-tumor effects. The following in vitro experiments demonstrated that TNFAIP6 and TLR6 were significantly overexpressed, and P2RY13 and CYP27A1 were significantly downregulated in LUAD cell lines, compared to BEAS-2B cells. Knocking down TNFAIP6 in A549 and PC9 resulted in the upregulation of FAS, CCL3, and ICAM-1, and the downregulation of CCL2, CXCR4, and VEGF-A in neutrophils when co-culturing with the conditioned medium (CM) from LUAD cells. Knocking down TNFAIP6 in LUAD also led to an elevated early apoptosis rate of neutrophils. Therefore, overexpressed TNFAIP6 in LUAD cancer cells might lead to neutrophils "N2" polarization, which exhibited pro-tumor effects. Further research based on the genes identified in this pilot study might shed light on neutrophils' effects on LUAD in the future.

Advanced hemodynamic and cluster analysis for identifying novel RV function subphenotypes in patients with pulmonary hypertension

BackgroundQuantifying right ventricular (RV) function is important to describe the pathophysiology of in pulmonary hypertension (PH). Current phenotyping strategies in PH rely on few invasive hemodynamic parameters to quantify RV dysfunction severity. The aim of this study was to identify novel RV phenotypes using unsupervised clustering methods on advanced hemodynamic features of RV function. MethodsParticipants were identified from the University of Arizona Pulmonary Hypertension Registry (n=190). RV-pulmonary artery coupling (Ees/Ea), RV systolic (Ees) and diastolic function (Eed) was quantified from stored RV pressure waveforms. Consensus clustering analysis with bootstrapping was used to identify the optimal clustering method. Pearson correlation analysis was used to reduce collinearity between variables. RV cluster subphenotypes were characterized using clinical data and compared to pulmonary vascular resistance (PVR) quintiles. ResultsFive distinct RV clusters (C1-C5) with distinct RV subphenotypes were identified using k-medoids with a Pearson distance matrix. Clusters 1 and 2 both have low diastolic stiffness (Eed) and afterload (Ea) but RV-PA coupling (Ees/Ea) is decreased in C2. Intermediate cluster (C3) has a similar Ees/Ea as C2 but with higher PA pressure and afterload. Clusters C4 and C5 have increased Eed and Ea but C5 has a significant decrease in Ees/Ea. Cardiac output was high in C3 distinct from the other clusters. In the PVR quintiles, contractility increased and stroke volume decreased as a function of increased afterload. World Symposium PH classifications were distributed across clusters and PVR quintiles. ConclusionsRV-centric phenotyping offers an opportunity for a more precise-medicine based management approach.

Development of a ferroptosis-based molecular markers for predicting RFS in prostate cancer patients

The goal of this study was to develop a ferroptosis-based molecular signature that can predict recurrence-free survival (RFS) in patients with prostate cancer (PCa). In this study, we obtained ferroptosis-related genes (FRGs) in FerrDb database and clinical transcriptome data in TCGA database and GEO database. Consensus cluster analysis was used to identify three molecular markers of ferroptosis in PCa with differential expression of 40 FRGs, including PD-L1 expression levels. We conducted a new ferroptosis-related signature for PCa RFS using four FRGs identified through univariate and multivariate Cox regression analyses. The signature was validated in the training, testing, and validation cohorts, and it demonstrated remarkable results in the area under the time-dependent receiver operating characteristic (ROC) curve of 0.757, 0.715, and 0.732, respectively. Additionally, we observed that younger patients, those with stage T III and stage T IV, stage N0, cluster 1, and cluster 2 PCa were more accurately predicted by the signature as independent predictors of RFS. DU-145 and RWPE-1 cells were successfully analyzed by qRT-PCR and Western blot for ASNS, GPT2, RRM2, and NFE2L2. In summary, we developed a novel ferroptosis-based signature for RFS in PC, utilizing four FRGs identified through univariate and multivariate Cox regression analyses. This signature was rigorously validated across training, testing, and validation cohorts, demonstrating exceptional performance as evidenced by its ROC curves. Notably, our findings indicate that this signature is particularly effective as an independent predictor of RFS in younger patients or those with stage T III and T IV, stage N0, and in clusters 1 and 2. Finally, we confirmed the expression of these four FRGs in DU-145 and RWPE-1 cell lines.

Comprehensive Analysis of N6-Methyladenosine RNA Methylation Regulators in the Diagnosis and Subtype Classification of Rheumatoid Arthritis.

m6A modification is the most abundant mRNA modifications and plays an integral role in various biological processes in eukaryotes. However, the role of m6A regulators in rheumatoid arthritis remains unknown. To determine the expression of m6A RNA methylation regulators in rheumatoid arthritis and their possible functional and prognostic value. In this study, we performed differential analysis in the comprehensive gene expression database GSE93272 dataset between non-rheumatoid arthritis patients and rheumatoid arthritis patients to obtain 15 important m6A regulators. A random forest model and lasso regression were used to screen the five most important m6A regulators to predict the risk of developing rheumatoid arthritis. After further validation using in vitro qPCR experiments, a nomogram model was developed based on the four most important m6A regulators (ELAVL1, WTAP, YTHDF1, and ALKBH5). Immuno-infiltration analysis and consensus clustering analysis were then performed. An analysis of the decision curve showed that the nomogram model could be beneficial to patients. According to selected important m6A regulators, patients with rheumatoid arthritis were classified into two m6A models (ClusterA and ClusterB) via consensus approach. Activated B cells, CD56dim natural killer cells, immature B cells, monocytes, natural killer T cells, and T lymphocytes were associated with ClusterA in immune infiltration analysis. Importantly, immune infiltration in patients with high ELAVL1 expression was strikingly similar to ClusterA. m6A regulators play a non-negligible role in the development of rheumatoid arthritis. A study of m6A patterns may provide future therapeutic options for rheumatoid arthritis.

Integrated analysis of necroptosis related gene signature to predict clinical outcomes, immune status and drug sensitivity in lower grade Glioma

BackgroundThe treatment of lower grade gliomas (LGG) is currently the most challenging dilemma in the management of intracranial tumors. Necroptosis is a type of programmed cell death that is closely associated with tumor progression, However, the role of necroptosis related genes in LGG is not yet well elucidated. MethodsOnline databases were used to obtain gene expression and clinical information. After gene differential expression analysis, a risk score model based on prognostic differentially expressed necroptosis-related genes (DENGs) were constructed to predict prognosis for LGG patients. The validity of the risk score model was then assessed with Kaplan-Meier survival curve. The prognostic DENGs included in the risk score model were then subjected to gene expression analysis, functional enrichment analysis, consensus clustering analysis, and single cell sequencing analysis. Finally, we investigated the correlation of the risk score and immune infiltration in LGG tumor microenvironment and drug sensitivity for LGG patients in different risk groups. ResultsA survival risk score model was constructed based on seven prognostic DENGs, which demonstrated satisfactory performance in predicting the prognosis of LGG patients. According to functional enrichment analyses, these seven DENGs may play a regulatory role in LGG tumorigenesis through several immune and metabolic pathways. LGG patients could be categorized into two clusters with distinct prognosis and clinicopathologic characteristics based on the expression of seven DENGs. Single-cell sequencing analysis demonstrated that the DENG signature was differentially expressed in various types of cells in LGG and may play a vital role in oncogenesis. Additionally, drug sensitivity analysis suggested that the seven-gene signature could guide clinical medication for LGG patients. ConclusionOur study developed a reliable necroptosis-related signature to predict the prognosis of LGG patients. This gene signature may also help estimate immune status and anti-cancer drug sensitivity in LGG patients. Our findings may pave the way to enhance our understanding of necroptosis in LGG.

Intrinsic circuitry of the rhombicbrain (central nervous system's intermediate sector) in a mammal.

The rhombicbrain (rhombencephalon or intermediate sector) is the vertebrate central nervous system part between the forebrain-midbrain (rostral sector) and spinal cord (caudal sector), and it has three main divisions: pons, cerebellum, and medulla. Using a data-driven approach, here we examine intrinsic rhombicbrain (intrarhombicbrain) network architecture that in rat consists of 52,670 possible axonal connections between 230 gray matter regions (115 bilaterally symmetrical pairs). Our analysis indicates that only 8,089 (15.4%) of these connections exist. Multiresolution consensus cluster analysis yields a nested hierarchy model of rhombicbrain subsystems that at the top level are associated with 1) the cerebellum and vestibular nuclei, 2) orofacial-pharyngeal-visceral integration, and 3) auditory connections; the bottom level has 68 clusters, ranging in size from 2 to 11 regions. The model provides a basis for functional hypothesis development and interrogation. More granular network analyses performed on the intrinsic connectivity of individual and combined main rhombicbrain divisions (pons, cerebellum, medulla, pons + cerebellum, and pons + medulla) demonstrate the mutability of network architecture in response to the addition or subtraction of connections. Clear differences between the structure-function network architecture of the rhombicbrain and forebrain-midbrain are discussed, with a stark comparison provided by the subsystem and small-world organization of the cerebellar cortex and cerebral cortex. Future analysis of the connections within and between the forebrain-midbrain and rhombicbrain will provide a model of brain neural network architecture in a mammal.

Construction and validation of a prognostic signature based on seven endoplasmic reticulum stress-related lncRNAs for patients with head and neck squamous cell carcinoma

Endoplasmic reticulum stress (ERS) occurs when misfolded or unfolded proteins accumulate in the endoplasmic reticulum (ER), and it is often observed in tumors, including head and neck squamous cell carcinoma (HNSCC). Relevant studies have demonstrated the prognostic significance of ERS-related long non-coding RNAs (lncRNAs) in various cancers. However, the relationship between ERS and lncRNAs in HNSCC has received limited attention in previous studies. In this study, we aimed to develop an ERS-related lncRNAs prognostic model using correlation analysis, Cox regression analysis, least absolute shrinkage, and selection operator (LASSO) regression analysis based on data from The Cancer Genome Atlas (TCGA) database. The survival and predictive ability of this model were evaluated using Kaplan–Meier analysis and time-dependent receiver operating characteristics (ROC), while nomograms and calibration curves were constructed. Then, functional enrichment analyses, tumor mutation burden (TMB), tumor infiltration of immune cells, single sample Gene Set Enrichment Analysis (ssGSEA), and drug sensitivity analysis were performed. Additionally, we conducted a consensus cluster analysis to compare differences between subtypes of tumors. Finally, we validated the expression of the ERS-related lncRNAs that constructed prognostic risk score model in HNSCC tissues through quantitative real-time PCR (qRT-PCR). We developed a prognostic signature based on seven ERS-related lncRNAs, which showed better predictive performance than other clinicopathological features. The high-risk poor prognosis group had a poorer prognosis in comparison to the low-risk good prognosis. The area under the ROC curve (AUC) predicted by this model for 3-year survival rates of HNSCC patients was 0.805. Enrichment analysis revealed that the differentially expressed genes were primarily enriched in pathways related to immune responses and signal transduction. Low-risk patients had lower TMB, more immune cell infiltrations, and enhanced anti-tumor immunity. Cluster analysis indicated that cluster 3 may have a better prognosis and immunotherapy effect. In addition, the result of qRT-PCR was consistent with our analysis. This prognostic model based on seven ERS-related lncRNAs is a promising tool for risk stratification, survival prediction, and immune cell infiltration status assessment.

The identification of N6-methyladenosine-related miRNAs predictive of hepatocellular carcinoma prognosis and immunotherapy efficacy.

Hepatocellular carcinoma (HCC) has a high degree of malignancy and poor prognosis. N6-methyladenosine (m6A) modifications and microRNAs (miRNAs) play pivotal roles in tumorigenesis and development. However, the role of m6A-related miRNAs in HCC has not been clarified yet. This study aimed to identify the role of m6A-miRNAs in HCC prognosis through bioinformatics analysis. The clinicopathological information and RNA sequencing data of 369 HCC tumor tissues and 49 tumor-adjacent tissues were downloaded from the TCGA database. A total of 23 m6A regulators were extracted to evaluated the m6A-related miRNAs using Pearson's correlation analysis. Then, we selected prognosis-related m6A-miRNAs using a univariate Cox regression model and used the consensus cluster analysis to explore the characteristics of the m6A-miRNAs. The coefficient of the least absolute shrinkage and selection operator (LASSO) Cox regression was applied to construct a prognostic risk score model. The receiver operated characteristic (ROC) analysis was applied to evaluate the prognostic value of the signature. The biological functions of targeted genes were predicted by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses. Then, to validate the potential predictive value for prognosis, the miRNA expression profiles from the GSE76903 and GSE6857 were used. Single sample Gene Set Enrichment Analysis (ssGSEA) and Estimation of Stromal and Immune cells in Malignant Tumor tissues using Expression data (ESTIMATE) were applied to assess the immune microenvironment of HCC. Additionally, a meta-analysis was used to verify the prognostic value of the m6A-microRNAs. RT-PCR was applied to validated the expression of miRNAs in HCC tissues. Cell viability, transwell assay and RNA m6A dot blot assays of HCC cells was applied to access the function of miR-17-5p. The expression of 48 m6A-related miRNAs was identified and 17 prognostic m6A-miRNAs was discovered. The expression profile of those 17 miRNAs was divided into three clusters, and these clusters were associated with the tumor microenvironment (TME) and prognosis. The nine m6A-related miRNA signature was associated with the prognosis of HCC, the AUC of the ROC was 0.771(TCGA dataset), 0.788(GSE76903) and 0.646(GSE6857). The TME and the expression of immune checkpoint molecules were associated with the risk score. The meta-analysis also validated the prognostic value of the m6A-related miRNAs (miR182-5p (HR:1.58, 95%CI:1.04-2.40) and miR-17-5p (HR:1.58, 95%CI: 1.04-2.40)). The expression of miR-17-5p was upregulated in HCC tissues and miR-17-5p showed an oncogenic role in HCC cells. The clinical innovation is the use of m6A-miRNAs as biomarkers for predicting prognosis regarding immunotherapy response in HCC patients.

Dry and wet experiments reveal diagnostic clustering and immune landscapes of cuproptosis patterns in patients with ankylosing spondylitis

Cuproptosis is a new manner of mitochondrial cell death induced by copper. There is evidence that serum copper has a crucial impact on ankylosing spondylitis (AS) by copper-induced inflammatory response. However, the molecular mechanisms of cuproptosis modulators in AS remain unknown. We aimed to use a bioinformatics-based method to comprehensively investigate cuproptosis-related subtype identification and immune microenvironment infiltration of AS. Additionally, we further verified the results by in vitro experiments, in which peripheral blood and fibroblast cells from AS patients were used to evaluate the functions of significant cuproptosis modulators on AS. Finally, eight significant cuproptosis modulators were identified by analysis of differences between controls and AS cases from GSE73754 dataset. Eight prognostic cuproptosis modulators (LIPT1, DLD, PDHA1, PDHB, SLC31A1, ATP7A, MTF1, CDKN2A) were identified using a random forest model for prediction of AS risk. A nomogram model of the 8 prognostic cuproptosis modulators was then constructed; the model could be beneficial in clinical settings, as indicated by decision curve analysis. Consensus clustering analysis was used to divide AS patients into two cuproptosis subtypes (clusterA & B) according to significant cuproptosis modulators. The cuproptosis score of each sample was calculated by principal component analysis to quantify cuproptosis subtypes. The cuproptosis scores were higher in clusterB than in clusterA. Additionally, cases in clusterA were closely associated with the immunity of activated B cells, Activated CD4 T cell, Type17 T helper cell and Type2 T helper cell, while cases in clusterB were linked to Mast cell, Neutrophil, Plasmacytoid dendritic cell immunity, indicating that clusterB may be more correlated with AS. Notably, key cuproptosis genes including ATP7A, MTF1, SLC31A1 detected by RT-qPCR with peripheral blood exhibited significantly higher expression levels in AS cases than controls; LIPT1 showed the opposite results; High MTF1 expression is correlated with increased osteogenic capacity. In general, this study of cuproptosis patterns may provide promising biomarkers and immunotherapeutic strategies for future AS treatment.

Multiple programmed cell death patterns and immune landscapes in bladder cancer: Evidence based on machine learning and multi-cohorts.

Bladder cancer (BLCA) is the most prevalent malignant neoplasm of the urinary tract, and ranks seventh as the most frequent systemic neoplasm in males. Dysregulation of programmed cell death (PCD) has been implicated in various stages of cancer progression, including tumorigenesis, invasion, and metastasis. However, the correlation between multiple PCD modes and BLCA is lacking. Thus, a risk prediction model was built based on 12 models of PCD to predict prognosis and immunotherapy response in patients with BLCA. The RNA sequencing transcriptome data of BLCA were collected from the Cancer Genome Atlas Program (TCGA) and GEO datasets. Univariate Cox and LASSO regression analyzes were performed to identify PCD-related genes (PCDRGs) significant for prognosis. Multivariate Cox regression analysis was used to develop a prognostic model for PCD. Survival analysis and chi-squared test were employed to analyze the survival variations between different risk groups. Univariate and multivariate Cox analyses were performed to evaluate the model as an independent prognostic predictor. A nomogram was formulated using both clinical data and the model to predict the survival rates of BLCA patients. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were performed to analyze and elucidate the molecular mechanisms and pathways operating within different risk score groups. Furthermore, the immune landscape was investigated and the efficacy of various anti-tumor drugs was evaluated for BLCA. Finally, consensus clustering analysis was adopted to explore the association between different PCD clusters and clinical characteristics. Assessment of the public datasets and multivariate Cox analysis yielded 1254 PCDRGs, of which 10 PCDRGs for BLCA were identified. Based on the PCDRGs, a prognostic model was built for BLCA patient prognosis. Compared with the low-risk group, the high-risk group had a poorer prognosis. The model predicted area under the curve (AUC) values of 0.751, 0.753, and 0.763, respectively, for 1-, 3-, and 5-year survival of BLCA patients. The nomogram further demonstrated the credibility of the prognosis model. The low-risk group patients exhibited lower TIDE scores and higher TMB scores, implying better response of the low-risk group to immunotherapy. The consensus clustering analysis indicated that compared with PCD cluster A, PCD cluster B was significantly more expressed in PCDRGs, suggesting a closer relation of PCD cluster B to PCDRGs. Patients in PCD cluster B had lower risk scores. To summarize, the effects of 12 PCD patterns on BLCA were synthesized and the correlation between PCD and BLCA was explored. These findings provide new and convincing evidence for individualized treatment of BLCA, and help guide the treatment strategy and improve the prognosis of BLCA patients.

Multi-dimension metabolic prognostic model for gastric cancer.

Metabolic reprogramming is involved in different stages of tumorigenesis. There are six widely recognized tumor-associated metabolic pathways, including cholesterol catabolism process, fatty acid metabolism, glutamine metabolic process, glycolysis, one carbon metabolic process, and pentose phosphate process. This study aimed to classify gastric cancer patients into different metabolic bio-similar clusters. We analyzed six tumor-associated metabolic pathways and calculated the metabolic pathway score through RNA-seq data using single sample gene set enrichment analysis. The consensus clustering analysis was performed to classify patients into different bio-similar clusters by multi-dimensional scaling. Kaplan-Meier curves were presented between different metabolic bio-similar groups for OS analysis. A training set of 370 patients from the Cancer Genome Atlas database with primary gastric cancer was chosen. Patients were classified into four metabolic bio-similar clusters, which were identified as metabolic non-specificity, metabolic-active, cholesterol-silence, and metabolic-silence clusters. Survival analysis showed that patients in metabolic-active cluster and metabolic-silence cluster have significantly poor prognosis than other patients (p=0.031). Patients in metabolic-active cluster and metabolic-silence cluster had significantly higher intra-tumor heterogeneity than other patients (p=0.032). Further analysis was performed in metabolic-active cluster and cholesterol-silence cluster. Three cell-cycle-related pathways, including G2M checkpoints, E2F targets, and MYC targets, were significantly upregulated in metabolic-active cluster than in cholesterol-silence cluster. A validation set of 192 gastric cancer patients from the Gene Expression Omnibus data portal verified that metabolic bio-similar cluster can predict prognosis in gastric cancer. Our study established a multi-dimension metabolic prognostic model in gastric cancer, which may be feasible for predicting clinical outcome.

Leveraging a disulfidptosis‑related lncRNAs signature for predicting the prognosis and immunotherapy of glioma

BackgroundGliomas, a prevalent form of primary brain tumors, are linked with a high mortality rate and unfavorable prognoses. Disulfidptosis, an innovative form of programmed cell death, has received scant attention concerning disulfidptosis-related lncRNAs (DRLs). The objective of this investigation was to ascertain a prognostic signature utilizing DRLs to forecast the prognosis and treatment targets of glioma patients.MethodsRNA-seq data were procured from The Cancer Genome Atlas database. Disulfidptosis-related genes were compiled from prior research. An analysis of multivariate Cox regression and the least absolute selection operator was used to construct a risk model using six DRLs. The risk signature’s performance was evaluated via Kaplan-Meier survival curves and receiver operating characteristic curves. Additionally, functional analysis was carried out using GO, KEGG, and single-sample GSEA to investigate the biological functions and immune infiltration. The research also evaluated tumor mutational burden, therapeutic drug sensitivity, and consensus cluster analysis. Reverse transcription quantitative PCR was conducted to validate the expression level of DRLs.ResultsA prognostic signature comprising six DRLs was developed to predict the prognosis of glioma patients. High-risk patients had significantly shorter overall survival than low-risk patients. The robustness of the risk model was validated by receiver operating characteristic curves and subgroup survival analysis. Risk model was used independently as a prognostic indicator for the glioma patients. Notably, the low-risk patients displayed a substantial decrease in the immune checkpoints, the proportion of immune cells, ESTIMATE and immune score. IC50 values from the different risk groups allowed us to discern three drugs for the treatment of glioma patients. Lastly, the potential clinical significance of six DRLs was determined.ConclusionsA novel six DRLs signature was developed to predict prognosis and may provide valuable insights for patients with glioma seeking novel immunotherapy and targeted therapy.

Disulfidptosis-related signatures for prognostic and immunotherapy reactivity evaluation in hepatocellular carcinoma

BackgroundHepatocellular carcinoma (HCC) is one of the most common cancers in the world and a nonnegligible health concern on a worldwide scale. Disulfidptosis is a novel mode of cell death, which is mainly caused by the collapse of the actin skeleton. Although many studies have demonstrated that various types of cell death are associated with cancer treatment, the relationship between disulfidptosis and HCC has not been elucidated.MethodsHere, we mainly applied bioinformatics methods to construct a disulfidptosis related risk model in HCC patients. Specifically, transcriptome data and clinical information were downloaded from the Gene Expression Omnibus (GEO), International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA) database. A total of 45 co-expressed genes were extracted between the disulfidptosis-related genes (DRGs) and the differential expression genes (DEGs) of liver hepatocellular carcinoma (LIHC) in the TCGA database. The LIHC cohort was divided into two subgroups with different prognosis by k-mean consensus clustering and functional enrichment analysis was performed. Subsequently, three hub genes (CDCA8, SPP2 and RDH16) were screened by Cox regression and LASSO regression analysis. In addition, a risk signature was constructed and the HCC cohort was divided into high risk score and low risk score subgroups to compare the prognosis, clinical features and immune landscape between the two subgroups. Finally, the prognostic model of independent risk factors was constructed and verified.ConclusionsHigh DRGs-related risk score in HCC individuals predict poor prognosis and are associated with poor immunotherapy response, which indicates that risk score assessment model can be utilized to guide clinical treatment strategy.

Integrated analysis identities Rho GTPases related molecular map in patients with gastric carcinoma

The intricate involvement of Rho GTPases in a multitude of human malignancies and their diverse array of biological functions has garnered substantial attention within the scientific community. However, their expression pattern and potential role in gastric cancer (GC) remain unclear. In this study, we successfully identified two distinct subtypes associated with Rho GTPase-related gene (RGG) through consensus clustering analysis, which exhibited significant disparities in overall survival and the tumor microenvironment. Subsequently, an extensively validated risk model termed RGGscore was meticulously constructed to prognosticate the outcomes of GC patients. This model was further assessed and validated using an external cohort. Notably, the high RGGscore group was indicative of a poorer prognosis. Univariate and multivariate Cox regression analyses unveiled the RGGscore as an autonomous prognostic indicator for GC patients. Subsequent external validation, utilizing two cohorts of patients who underwent immunotherapy, demonstrated a significant correlation between a low RGGscore and improved response to immunotherapy. Additionally, the expression levels of three genes associated with RGGscore were examined using qRT-PCR. Taken together, a pioneering RGGscore model has been successfully established, showcasing its potential efficacy in offering valuable therapeutic guidance for GC.

Re-assessing prolonged cold ischemia time in kidney transplantation through machine learning consensus clustering.

We aimed to cluster deceased donor kidney transplant recipients with prolonged cold ischemia time (CIT) using an unsupervised machine learning approach. We performed consensus cluster analysis on 11615 deceased donor kidney transplant patients with CIT exceeding 24h using OPTN/UNOS data from 2015 to 2019. Cluster characteristics of clinical significance were identified, and post-transplant outcomes were compared. Consensus cluster analysis identified two clinically distinct clusters. Cluster 1 was characterized by young, non-diabetic patients who received kidney transplants from young, non-hypertensive, non-ECD deceased donors with lower KDPI scores. In contrast, the patients in cluster 2 were older and more likely to have diabetes. Cluster 2 recipients were more likely to receive transplants from older donors with a higher KDPI. There was lower use of machine perfusion in Cluster 1 and incrementally longer CIT in Cluster 2. Cluster 2 had a higher incidence of delayed graft function (42%vs. 29%), and lower 1-year patient (95%vs. 98%) and death-censored (95%vs. 97%) graft survival compared to Cluster 1. Unsupervised machine learning characterized deceased donor kidney transplant recipients with prolonged CIT into two clusters with differing outcomes. Although Cluster 1 had more favorable recipient and donor characteristics and better survival, the outcomes observed in Cluster 2 were also satisfactory. Overall, both clusters demonstrated good survival suggesting opportunities for transplant centers to incrementally increase CIT.

Machine learning and molecular subtype analyses provide insights into PANoptosis-associated genes in rheumatoid arthritis

BackgroundPANoptosis represents a newly identified form of programmed cell death that plays a significant role in the autoimmune diseases. Rheumatoid arthritis (RA) is characterized by the presence of autoantibodies. Nevertheless, the specific biomarkers and molecular mechanisms responsible for the apoptotic characteristics of RA remain largely uninvestigated.MethodsWe utilized 8 synovial tissue RA datasets. We selected genes associated with PANoptosis from the GeneCard database. By employing the limma, WGCNA, and machine learning algorithms we identified core genes. We utilized consensus clustering analysis to identify distinct PANoptosis subtypes of RA. Boruta algorithm was employed to construct a PANoptosis signature score. The sensitivity of distinct subtypes to drug treatment was verified using an independent dataset.ResultsThe SPP1 emerged as the significant gene, with its elevated expression in RA patients. We identified two PANoptosis RA subtypes. Cluster 1 showed high expression of Tregs, resting dendritic cells, and resting mast cells. Cluster 2 exhibited high expression of CD4 memory T cells and follicular helper T cells. Cluster 2 exhibited a higher degree of sensitivity towards immune checkpoint therapy. Employing the Boruta algorithm, a subtype score was devised for 37 PANoptosis genes, successfully discerning the subtypes (AUC = 0.794), wherein patients with elevated scores demonstrated enhanced responsiveness to Rituximab treatment.ConclusionOur analysis revealed that SPP1 holds potential biomarker for the diagnosis of RA. Cluster 2 exhibited enhanced sensitivity to immune checkpoint therapy, higher PANoptosis scores, and improved responsiveness to drug treatment. This study offers potential implications in the realm of diagnosis and treatment.

TGF-β signaling pathway-related genes in predicting the prognosis of colon cancer and guiding immunotherapy

BackgroundColon cancer is a malignant tumor with high malignancy and a low survival rate whose heterogeneity limits systemic immunotherapy. Transforming growth factor-β (TGF-β) signaling pathway-related genes are associated with multiple tumors, but their role in prognosis prediction and tumor microenvironment (TME) regulation in colon cancer is poorly understood. Using bioinformatics, this study aimed to construct a risk prediction signature for colon cancer, which may provide a means for developing new effective treatment strategies. MethodsUsing consensus clustering, patients in The Cancer Genome Atlas (TCGA) with colon adenocarcinoma were classified into several subtypes based on the expression of TGF-β signaling pathway-related genes, and differences in survival, molecular, and immunological TME characteristics and drug sensitivity were examined in each subtype. Ten genes that make up a TGF-β-related predictive signature were found by least absolute shrinkage and selector operation (LASSO) regression using colon cancer data from the TCGA database and confirmed using a Gene Expression Omnibus (GEO) dataset. A nomogram incorporating risk scores and clinicopathologic factors was developed to stratify the prognosis of patients with colon cancer for accurate clinical diagnosis and therapy. ResultsTwo TGF-β subtypes were identified, with the TGF-β-high subtype being associated with a poorer prognosis and superior sensitivity to immunotherapy. Mutation analyses showed a high incidence of gene mutations in the TGF-β-high subtype. After completing signature construction, patients with colon cancer were categorized into high- and low-risk subgroups based on the median risk score of the TGF-β-related predictive signature. The risk score exhibited superior predictive performance relative to age, gender, and stage, as evidenced by its AUC of 0.686. Patients in the high-risk subgroup had higher levels of immunosuppressive cell infiltration and immune checkpoints in the TME, suggesting that these patients had better responses to immunotherapy. ConclusionsPatients with colon cancer were divided into two subtypes with different survival and immune characteristics using consensus clustering analysis based on TGF-β signaling pathway-related genes. The constructed risk prediction signature may show promise as a biomarker for evaluating the prognosis of colon cancer, with potential utility for screening individuals for immunotherapy.