Gene Expression Omnibus Cohorts Research Articles

Ferroptosis-related gene signature predicts prognosis and immune microenvironment in prostate cancer.

Ferroptosis, an iron-dependent form of programmed cell death, significantly impacts cancer, yet its link to prostate cancer (PCa) prognosis remains underexplored. This study aims to develop and validate a ferroptosis-related gene signature to predict PCa prognosis and immune microenvironment differences, potentially identifying therapeutic targets. RNA-sequencing data of 478 PCa patients and corresponding clinical data were downloaded from The Cancer Genome Atlas (TCGA) database. We investigated the disease-free survival (DFS) rates of the high- and low-risk groups using the Kaplan-Meier method. Functional differences between the high- and low-risk groups were investigated by a gene set enrichment analysis (GSEA), and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. The link between ferroptosis risk score and immune status was examined using CIBERSORT. The expression levels of core prognostic genes in benign prostatic hyperplasia (BPH) and PCa were verified using quantitative real-time polymerase chain reaction (qRT-PCR), Western blot, and immunohistochemistry (IHC). A novel ferroptosis-related prognostic gene signature was established and tested in the Gene Expression Omnibus (GEO) database based on univariate and multivariate Cox regression analyses. Patients with PCa were classified into high- and low-risk groups based on this ferroptosis signature. Patients in the high-risk group had worse outcomes than those in the low-risk group. The predictive accuracy of the model was demonstrated by a receiver operating characteristic (ROC) analysis. An additional enrichment analysis of TCGA cohort revealed the immune-related pathways were significantly upregulated in the high-risk group, with areas under the curve (AUCs) of 0.85 at 1 year, 0.82 at 3 years, and 0.76 at 5 years. In the GEO cohort, the AUCs reached 0.69 at 1 year, 0.74 at 3 years, and 0.75 at 5 years. An additional enrichment analysis indicated a significant upregulation of cytokine-related pathways, immune receptor activity, and other immune-related pathways in the high-risk group. Furthermore, the analysis revealed that the proportions of mast cells and plasma cells were significantly lower in the high-risk group compared to the low-risk group of PCa patients. Conversely, the proportion of regulatory T cells (Tregs) was significantly higher in the high-risk group than in the low-risk group. According to the qRT-PCR, Western blot, and IHC results, DRD4, SRC, AKR1C2, and AIFM2 expression was significantly higher in PCa than BPH. We also showed that the ferrostatin 1-treated LNCaP cells had higher expression levels of DRD4, SRC, and AKR1C2. A prognostic signature of eight ferroptosis-related genes (FRGs) that may accurately predict PCa patient outcomes was constructed and validated. FRGs may contribute to anti-tumor immunity and serve as therapeutic targets in PCa.

Comprehensive characterization of cytokines in patients under extracorporeal membrane oxygenation: Evidence from integrated bulk and single-cell RNA sequencing data using multiple machine learning approaches.

ECMO (extracorporeal membrane oxygenation) is an effective technique for providing short-term mechanical support to the heart, lungs, or both. During ECMO treatment, the inflammatory response, particularly involving cytokines, plays a crucial role in pathophysiology. However, the potential effects of cytokines on patients receiving ECMO are not comprehensively understood. We acquired three ECMO support datasets, namely two bulk and one single-cell RNA sequencing (RNA-seq), from the GEO (Gene Expression Omnibus) combined with hospital cohorts to investigate the expression pattern and potential biological processes of cytokine-related genes (CRGs) in patients under ECMO. Subsequently, machine learning approaches, including support vector machine (SVM), random forest (RF), modified Lasso penalized regression, extreme gradient boosting (XGBoost), and artificial neural network (ANN), were applied to identify hub CRGs, thus developing a prediction model called CRG classifier. The predictive and prognostic performance of the model was comprehensively evaluated in GEO and hospital cohorts. Finally, we mechanistically analyzed the relationship between hub cytokines, immune cells, and pivotal molecular pathways. Analyzing bulk and single-cell RNA-seq data revealed that most CRGs were significantly differentially expressed, the enrichment scores of cytokine and cytokine cytokine receptor (CCR) interaction were significantly higher during ECMO. Based on multiple machine learning algorithms, nine key CRGs (CCL2, CCL4, IFNG, IL1R2, IL20RB, IL31RA, IL4, IL7, and IL7R) were used to develop the CRG classifier. The CRG classifier exhibited excellent prognostic values (AUC > 0.85), serving as an independent risk factor. It performed better in predicting mortality and yielded a larger net benefit than other clinical features in GEO and hospital cohorts. Additionally, IL1R2, CCL4, and IL7R were predominantly expressed in monocytes, NK cells, and T cells, respectively. Their expression was significantly positively correlated with the relative abundance of corresponding immune cells. Gene set variation analysis (GSVA) revealed that parainflammation, complement and coagulation cascades, and IL6/JAK/STAT3 signaling were significantly enriched in the subgroup that died after receiving ECMO. Spearman correlation analyses and Mantel tests revealed that the expression of hub cytokines (IL1R2, CCL4, and IL7R) and pivotal molecular pathways scores (complement and coagulation cascades, IL6/JAK/STAT3 signaling, and parainflammation) were closely related. A predictive model (CRG classifier) comprising nine CRGs based on multiple machine learning algorithms was constructed, potentially assisting clinicians in guiding individualized ECMO treatment. Additionally, elucidating the underlying mechanistic pathways of cytokines during ECMO will provide new insights into its treatment.

Establishment of an ovarian cancer exhausted CD8+T cells-related genes model by integrated analysis of scRNA-seq and bulk RNA-seq

Ovarian cancer (OC) was the fifth leading cause of cancer death and the deadliest gynecological cancer in women. This was largely attributed to its late diagnosis, high therapeutic resistance, and a dearth of effective treatments. Clinical and preclinical studies have revealed that tumor-infiltrating CD8+T cells often lost their effector function, the dysfunctional state of CD8+T cells was known as exhaustion. Our objective was to identify genes associated with exhausted CD8+T cells (CD8TEXGs) and their prognostic significance in OC. We downloaded the RNA-seq and clinical data from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. CD8TEXGs were initially identified from single-cell RNA-seq (scRNA-seq) datasets, then univariate Cox regression, the least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression were utilized to calculate risk score and to develop the CD8TEXGs risk signature. Kaplan–Meier analysis, univariate Cox regression, multivariate Cox regression, time-dependent receiver operating characteristics (ROC), nomogram, and calibration were conducted to verify and evaluate the risk signature. Gene set enrichment analyses (GSEA) in the risk groups were used to figure out the closely correlated pathways with the risk group. The role of risk score has been further explored in the homologous recombination repair deficiency (HRD), BRAC1/2 gene mutations and tumor mutation burden (TMB). A risk signature with 4 CD8TEXGs in OC was finally built in the TCGA database and further validated in large GEO cohorts. The signature also demonstrated broad applicability across various types of cancer in the pan-cancer analysis. The high-risk score was significantly associated with a worse prognosis and the risk score was proven to be an independent prognostic biomarker. The 1-, 3-, and 5-years ROC values, nomogram, calibration, and comparison with the previously published models confirmed the excellent prediction power of this model. The low-risk group patients tended to exhibit a higher HRD score, BRCA1/2 gene mutation ratio and TMB. The low-risk group patients were more sensitive to Poly-ADP-ribose polymerase inhibitors (PARPi). Our findings of the prognostic value of CD8TEXGs in prognosis and drug response provided valuable insights into the molecular mechanisms and clinical management of OC.

Identification of a basement membrane-related gene signature for predicting prognosis, immune infiltration, and drug sensitivity in colorectal cancer.

Colorectal cancer (CRC) is the most common malignancy affecting the gastrointestinal tract. Extensive research indicates that basement membranes (BMs) may play a crucial role in the initiation and progression of the disease. Data on the RNA expression patterns and clinicopathological information of patients with CRC were sourced from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. A BM-linked risk signature for the prediction of overall survival (OS) was formulated using univariate Cox regression and combined machine learning techniques. Survival outcomes, functional pathways, the tumor microenvironment (TME), and responses to both immunotherapy and chemotherapy within varying risk classifications were also investigated. The expression trends of the model genes were evaluated by reverse transcription polymerase chain reaction (RT-PCR) and the Human Protein Atlas (HPA) database. A nine-gene risk signature containing UNC5C, TINAG, TIMP1, SPOCK3, MMP1, AGRN, UNC5A, ADAMTS4, and ITGA7 was constructed for the prediction of outcomes in patients with CRC. The expression profiles of these candidate genes were verified using RT-PCR and the HPA database and were found to be consistent with the findings on differential gene expression in the TCGA dataset. The validity of the signature was confirmed using the GEO cohort. The patients were stratified into different risk groups according to differences in clinicopathological characteristics, TME features, enrichment functions, and drug sensitivities. Lastly, the prognostic nomogram model based on the risk score was found to be effective in identifying high-risk patients and predicting OS. A basement membrane-related risk signature was constructed and found to be effective for predicting the prognosis of patients with CRC.

Identification of immunotherapy-related subtypes, characterization of tumor microenvironment infiltration, and development of a prognostic signature in gastric carcinoma.

Recent advances in immunotherapy have elicited a considerable amount of attention as viable therapeutic options for several cancer types, the present study aimed to explore the immunotherapy-related genes (IRGs) and develop a prognostic risk signature in gastric carcinoma (GC) based on these genes. IRGs were identified by comparing immunotherapy responders and non-responders in GC. Then, GC patients were divided into distinct subtypes by unsupervised clustering method based on IRGs, and the differences in immune characteristics and prognostic stratification between these subtypes were analyzed. An immunotherapy-related risk score (IRRS) signature was developed and validated for risk classification and prognosis prediction based on The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) cohorts. Besides, the predictive ability of the IRRS in immunotherapy response was also determined. A total of 63 IRGs were identified, and 371 GC patients were stratified into two molecular subgroups with significantly different prognosis and immune characteristics. Then, an IRRS signature comprised of three IRGs (CENP8, NRP1, and SERPINE1) was constructed to predict the prognosis of GC patients in TCGA cohort. Importantly, external validation in multiple GEO cohorts further confirmed the universal applicability of the IRRS in distinct populations. Furthermore, we found that the IRRS was significantly correlated with patient's responsiveness to immunotherapy, GC patients with low IRRS are more likely to benefit from existing immunotherapy. The risk score could serve as a robust prognostic biomarker, provide therapeutic benefits for immunotherapy and may be helpful for clinical decision making in GC patients.

Cuproptosis-related DNA methylation signature predict prognosis and immune microenvironment in cutaneous melanoma

The prognosis for Cutaneous Melanoma (CM), a skin malignant tumor that is extremely aggressive, is not good. A recently identified type of controlled cell death that is intimately related to immunotherapy and the development of cancer is called cuproptosis. Using The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) database, we developed and validated a DNA-methylation located in cuproptosis death-related gene prognostic signature (CRG-located DNA-methylation prognostic signature) to predict CM’s prognosis. Kaplan–Meier analysis of our TCGA and GEO cohorts showed that high-risk patients had a shorter overall survival. The area under the curve (AUC) for the TCGA cohort was 0.742, while for the GEO cohort it was 0.733, according to the receiver operating characteristic (ROC) analysis. Furthermore, this signature was discovered as an independent prognostic indicator over CM patients based on Cox-regression analysis. Immunogenomic profiling indicated that majority immune-checkpoints got an opposite relationship with the signature, and patients in the group at low risk got higher immunophenoscore. Several immune pathways were enriched, according to functional enrichment analysis. In conclusion, a prognostic methylation signature for CM patients was established and confirmed. Because of its close relationship to the immune landscape, this signature may help clinicians make more accurate and individualized choices regarding therapy.

Prognositic value of anoikis and tumor immune microenvironment-related gene in the treatment of osteosarcoma.

Osteosarcoma is a highly aggressive primary malignant bone tumor commonly seen in children and adolescents, with a poor prognosis. Anchorage-dependent cell death (anoikis) has been proven to be indispensable in tumor metastasis, regulating the migration and adhesion of tumor cells at the primary site. However, as a type of programmed cell death, anoikis is rarely studied in osteosarcoma, especially in the tumor immune microenvironment. This study aims to clarify prognostic value of anoikis and tumor immune microenvironment-related gene in the treatment of osteosarcoma. Anoikis-related genes (ANRGs) were obtained from GeneCards. Clinical information and ANRGs expression profiles of osteosarcoma patients were sourced from the therapeutically applicable research to generate effective therapies and Gene Expression Omnibus (GEO) databases. ANRGs highly associated with tumor immune microenvironment were identified by the estimate package and the weighted gene coexpression network analysis (WGCNA) algorithm. Machine learning algorithms were performed to construct long-term survival predictive strategy, each sample was divided into high-risk and low-risk subgroups, which was further verified in the GEO cohort. Finally, based on single-cell RNA-seq from the GEO database, analysis was done on the function of signature genes in the osteosarcoma tumor microenvironment. A total of 51 hub ANRGs closely associated with the tumor microenvironment were identified, from which 3 genes (MERTK, BNIP3, S100A8) were selected to construct the prognostic model. Significant differences in immune cell activation and immune-related signaling pathways were observed between the high-risk and low-risk groups based on tumor microenvironment analysis (all P<0.05). Additionally, characteristic genes within the osteosarcoma microenvironment were identified in regulation of intercellular crosstalk through the GAS6-MERTK signaling pathway. The prognostic model based on ANRGs and tumor microenvironment demonstrate good predictive power and provide more personalized treatment options for patients with osteosarcoma.

Predictive model using four ferroptosis-related genes accurately predicts gastric cancer prognosis.

Gastric cancer (GC) is a common malignancy of the digestive system. According to global 2018 cancer data, GC has the fifth-highest incidence and the third-highest fatality rate among malignant tumors. More than 60% of GC are linked to infection with Helicobacter pylori (H. pylori), a gram-negative, active, microaerophilic, and helical bacterium. This parasite induces GC by producing toxic factors, such as cytotoxin-related gene A, vacuolar cytotoxin A, and outer membrane proteins. Ferroptosis, or iron-dependent programmed cell death, has been linked to GC, although there has been little research on the link between H. pylori infection-related GC and ferroptosis. To identify coregulated differentially expressed genes among ferroptosis-related genes (FRGs) in GC patients and develop a ferroptosis-related prognostic model with discrimination ability. Gene expression profiles of GC patients and those with H. pylori-associated GC were obtained from The Cancer Genome Atlas and Gene Expression Omnibus (GEO) databases. The FRGs were acquired from the FerrDb database. A ferroptosis-related gene prognostic index (FRGPI) was created using least absolute shrinkage and selection operator-Cox regression. The predictive ability of the FRGPI was validated in the GEO cohort. Finally, we verified the expression of the hub genes and the activity of the ferroptosis inducer FIN56 in GC cell lines and tissues. Four hub genes were identified (NOX4, MTCH1, GABARAPL2, and SLC2A3) and shown to accurately predict GC and H. pylori-associated GC. The FRGPI based on the hub genes could independently predict GC patient survival; GC patients in the high-risk group had considerably worse overall survival than did those in the low-risk group. The FRGPI was a significant predictor of GC prognosis and was strongly correlated with disease progression. Moreover, the gene expression levels of common immune checkpoint proteins dramatically increased in the high-risk subgroup of the FRGPI cohort. The hub genes were also confirmed to be highly overexpressed in GC cell lines and tissues and were found to be primarily localized at the cell membrane. The ferroptosis inducer FIN56 inhibited GC cell proliferation in a dose-dependent manner. In this study, we developed a predictive model based on four FRGs that can accurately predict the prognosis of GC patients and the efficacy of immunotherapy in this population.

Bioinformatics Analysis Reveals Prognostic Significance of the Macrophage Marker Gene Signature in Gastric Adenocarcinoma.

Gastric adenocarcinoma (GAC) is a malignant tumor with the highest incidence in the digestive system. Macrophages have been proven to play important roles in tumor microenvironment. Herein, single-cell RNA sequencing (scRNA-seq) profiles from the Gene Expression Omnibus (GEO) and bulk RNA-seq data from the Cancer Genome Atlas (TCGA) database were utilized to construct a macrophage marker gene signature (MMGS) to predict the prognosis of GAC patients. Subsequently, a risk score model based on the MMGS was built to predict the prognosis of GAC patients; further, this was validated in the GEO cohort. The risk score categorized patients into the high- and low-risk groups. A nomogram model based on the risk score and clinic-pathological characteristics was developed. Seven genes, ABCA1, CTHRC1, GADD45B, NPC2, PLTP, PRSS23, and RNASE1, were included in the risk score model. Patients with a low-risk score showed a better prognosis. The MMGS had good sensitivity and specificity for predicting the prognosis inGAC patients. The risk score was an independent prognostic factor. The constructed nomogram exhibited favorable predictability and reliability for predicting GAC prognosis. In conclusion, the risk score model based on the seven MMGSs performed well in the predicting prognosis of GAC patients. Our study may provide new insights into clinical decision-making for the personalized treatment of patients with gastric cancer (GC).

Identification of Novel Stemness-based Subtypes and Construction of a Prognostic Risk Model for Patients with Lung Squamous Cell Carcinoma.

Although cancer stem cells (CSCs) contribute to tumorigenesis, progression, and drug resistance, stemness-based classification and prognostic signatures of lung squamous cell carcinoma (LUSC) remain unclarified. This study attempted to identify stemness-based subtypes and develop a prognostic risk model for LUSC. Based on RNA-seq data from The Cancer Genome Atlas (TCGA), Gene-Expression Omnibus (GEO) and Progenitor Cell Biology Consortium (PCBC), mRNA expression-based stemness index (mRNAsi) was calculated by one-class logistic regression (OCLR) algorithm. A weighted gene coexpression network (WGCNA) was employed to identify stemness subtypes. Differences in mutation, clinical characteristics, immune cell infiltration, and antitumor therapy responses were determined. We constructed a prognostic risk model, followed by validations in GEO cohort, pan-cancer and immunotherapy datasets. LUSC patients with subtype C2 had a better prognosis, manifested by higher mRNAsi, higher tumor protein 53 (TP53) and Titin (TTN) mutation frequencies, lower immune scores and decreased immune checkpoints. Patients with subtype C2 were more sensitive to Imatinib, Pyrimethamine, and Paclitaxel therapy, whereas those with subtype C1 were more sensitive to Sunitinib, Saracatinib, and Dasatinib. Moreover, we constructed stemness-based signatures using seven genes (BMI1, CCDC51, CTNS, EIF1AX, FAM43A, THBD, and TRIM68) and found high-risk patients had a poorer prognosis in the TCGA cohort. Similar results were found in the GEO cohort. We verified the good performance of risk scores in prognosis prediction and therapy responses. The stemness-based subtypes shed novel insights into the potential roles of LUSC-stemness in tumor heterogeneity, and our prognostic signatures offer a promising tool for prognosis prediction and guide therapeutic decisions in LUSC.

DNMT3A promotes glioma growth and malignancy via TNF-α/NF-κB signaling pathway.

DNMT3A is the main molecule responsible for DNA methylation in cells. DNMT3A affects the progression of inflammation, degenerative diseases, and malignant tumors, and exhibits significant aberrantly expression in tumor tissues. Transcriptome data and relevant clinical information were downloaded from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), and Gene Expression Omnibus (GEO) datasets. Differential expression analysis and prognostic analysis were conducted based on above statistics. We constructed a clinical prognostic model and identified DNMT3A as an independent prognostic factor to accurately predict patient prognosis. Differential gene enrichment analysis revealed that DNMT3A affects the progression of glioma through multiple pathways, among which the tumor necrosis factor-α (TNF-α)/nuclear factor-kappa B (NF-κB) pathway shows a strong correlation. Immunological analysis also revealed a certain correlation between DNMT3A and tumor immunity. We demonstrated through gene editing that DNMT3A can affect the release of TNF-α in cells, thereby affecting the progression of glioma. Functional experiments have also demonstrated that DNMT3A plays a crucial role in tumors. RNA-sequencing and survival analyses of lower-grade glioma (LGG) patients in TCGA, CGGA, and GEO cohorts showed that high DNMT3A expression correlated with poor prognosis of LGG patients. Univariate and multivariate Cox regression analyses showed that DNMT3A expression was an independent prognostic indicator in LGG. The prognosis prediction nomogram with age, World Health Organization (WHO) grading, and DNMT3A expression showed reliable performance in predicting the 1-, 3-, and 5-year overall survival (OS) of LGG patients. Functional enrichment analysis, gene set enrichment analysis (GSEA), and ESTIMATE algorithm analyses showed that DNMT3A expression was associated with the tumor infiltration of immune cells and predicted response to immunotherapy in two immunotherapy cohorts of pan-cancer patients. Furthermore, short hairpin RNA (shRNA)-mediated knockdown of DNMT3A in the LGG cell lines suppressed proliferation, migration, and invasion of LGG cells by downregulating the TNF-α/NF-κB signaling pathway. Our data showed that DNMT3A was a potential prognostic biomarker in glioma. DNMT3A promoted proliferation and malignancy of LGG cells through the TNF-α/NF-κB signaling pathway. DNMT3A is a promising therapeutic target for treating patients with LGG.

A novel tumor-associated neutrophil gene signature for predicting prognosis, tumor immune microenvironment, and therapeutic response in breast cancer

Tumor-associated neutrophils (TANs) can promote tumor progression. This study aimed to investigate the molecular signature that predict the prognosis and immune response of breast cancer (BRCA) based on TAN-related gene (TANRG) expression data. The RNA-seq data of BRCA were gathered from The Cancer Genome Atlas (TCGA) and gene expression omnibus (GEO) datasets. Univariate Cox regression analysis and the least absolute shrinkage and selection operator for selecting prognostic genes. A neo-TAN-related risk signature was constructed by multivariate Cox regression analysis. Time-dependent receiver operating characteristic (ROC) curve analyses and Kaplan–Meier analyses were performed to validate the signature in GEO cohorts and the triple-negative breast cancer (TNBC) subtype. We constructed an independent prognostic factor model with 11 TANRGs. The areas under the ROC curve (AUCs) of the TCGA training cohorts for 3-, 5-, and 7-year overall survival were 0.72, 0.73, and 0.73, respectively. The AUCs of the GEO test cohorts for 3-, 5-, and 7-year overall survival were 0.83, 0.89, and 0.94 (GSE25066) and 0.67, 0.69, and 0.73 (GSE58812), respectively. The proportion of immune subtypes differed among the different risk groups. The IC50 values differed significantly between risk groups and can be used as a guide for systemic therapy. The prognostic model developed by TANRGs has excellent predictive performance in BRCA patients. In addition, this feature is closely related to the prediction of survival, immune activity and treatment response in BRCA patients.

Long noncoding RNA XIST: a novel independent prognostic biomarker for patients with ABC-DLBCL receiving R-CHOP treatment.

Approximately one-third of activated B-cell-like diffuse large B-cell lymphoma (ABC-DLBCL) cases were unresponsive to standard first-line therapy; thus, identifying biomarkers to evaluate therapeutic efficacy and assessing the emergence of drug resistance is crucial. Through early-stage screening, long noncoding RNA (lncRNA) X-inactive specific transcript (XIST) was found to be correlated with the R-CHOP treatment response. This study aimed to clarify the characteristics of XIST in ABC-DLBCL. The expression level of XIST in 161 patients with ABC-DLBCL receiving R-CHOP therapy was examined via RNA in situ hybridization, and the association between XIST expression and clinicopathological features, treatment response and prognosis was analyzed in the study cohort and validated in the Gene Expression Omnibus cohort. Cell biological experiments and bioinformatics analyses were conducted to reveal aberrant signaling. The proportion of complete response in patients with high XIST expression was lower than that in patients with low XIST expression (53.8% versus 77.1%) (P = 0.002). High XIST expression was remarkably associated with the characteristics of tumor progression and was an independent prognostic element for overall survival (P = 0.039) and progression-free survival (P = 0.027) in ABC-DLBCL. XIST was proven to be involved in m6A-related methylation and ATF6-associated autophagy. XIST knockdown repressed ABC-DLBCL cellular proliferation by regulating Raf/MEK/ERK signaling. High XIST expression was associated with ABC-DLBCL tumorigenesis and development and contributed to R-CHOP treatment resistance. XIST may be a promising signal to predict ABC-DLBCL prognosis.

Novel application of the ferroptosis-related genes risk model associated with disulfidptosis in hepatocellular carcinoma prognosis and immune infiltration.

Hepatocellular carcinoma (HCC) stands as the prevailing manifestation of primary liver cancer and continues to pose a formidable challenge to human well-being and longevity, owing to its elevated incidence and mortality rates. Nevertheless, the quest for reliable predictive biomarkers for HCC remains ongoing. Recent research has demonstrated a close correlation between ferroptosis and disulfidptosis, two cellular processes, and cancer prognosis, suggesting their potential as predictive factors for HCC. In this study, we employed a combination of bioinformatics algorithms and machine learning techniques, leveraging RNA sequencing data, mutation profiles, and clinical data from HCC samples in The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the International Cancer Genome Consortium (ICGC) databases, to develop a risk prognosis model based on genes associated with ferroptosis and disulfidptosis. We conducted an unsupervised clustering analysis, calculating a risk score (RS) to predict the prognosis of HCC using these genes. Clustering analysis revealed two distinct HCC clusters, each characterized by significantly different prognostic and immune features. The median RS stratified HCC samples in the TCGA, GEO, and ICGC cohorts into high-and low-risk groups. Importantly, RS emerged as an independent prognostic factor in all three cohorts, with the high-risk group demonstrating poorer prognosis and a more active immunosuppressive microenvironment. Additionally, the high-risk group exhibited higher expression levels of tumor mutation burden (TMB), immune checkpoints (ICs), and human leukocyte antigen (HLA), suggesting a heightened responsiveness to immunotherapy. A cancer stem cell infiltration analysis revealed a higher similarity between tumor cells and stem cells in the high-risk group. Furthermore, drug sensitivity analysis highlighted significant differences in response to antitumor drugs between the two risk groups. In summary, our risk prognostic model, constructed based on ferroptosis-related genes associated with disulfidptosis, effectively predicts HCC prognosis. These findings hold potential implications for patient stratification and clinical decision-making, offering valuable theoretical insights in this field.

Development and validation of glycosyltransferase related-gene for the diagnosis and prognosis of head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous cancer characterized by difficulties in early diagnosis and outcome prediction. Aberrant glycosylated structures produced by the aberrant expression of glycosyltransferases are prevalent in HNSCC. In this study, we aim to construct glycosyltransferase-related gene signatures with diagnostic and prognostic value to better stratify patients with HNSCC and improve their diagnosis and prognosis. Bioinformatic tools were used to process data of patients with HNSCC from The Cancer Genome Atlas (TCGA) database. The prognostic model was formatted using univariate and multivariate Cox regression methods, while the diagnostic signature was constructed using support vector machine (SVM) and LASSO analysis. The results were verified using the Gene Expression Omnibus (GEO) cohort. The tumor microenvironment and benefits of immune checkpoint inhibitor (ICI) therapy in subgroups defined by glycosyltransferase-related genes were analyzed. Molecular biology experiments, including western blotting, cell counting kit (CCK)-8, colony formation, wound healing, and Transwell assays, were conducted to confirm the oncogenic function of beta-1,4-galactosyltransferase 3 (B4GALT3) in HNSCC. We established a five-gene prognostic signature and a 15-gene diagnostic model. Based on the median risk score, patients with low risk had longer overall survival than those in the high-risk group, which was consistent with the results of the GEO cohort. The concrete results suggested that high-risk samples were related to a high tumor protein (TP)53 mutation rate, high infiltration of resting memory cluster of differentiation (CD)4 T cells, resting natural killer (NK) cells, and M0 macrophages, and benefited from ICI therapy. In contrast, the low-risk subgroup was associated with a low TP53 mutation rate; and high infiltration of naive B cells, plasma cells, CD8 T cells, and resting mast cells; and benefited less from ICI therapy. In addition, the diagnostic model had an area under curve (AUC) value of 0.997 and 0.978 in the training dataset and validation cohort, respectively, indicating the high diagnostic potential of the model. Ultimately, the depletion of B4GALT3 significantly hindered the proliferation, migration, and invasion of HNSCC cells. We established two new biomarkers that could provide clinicians with diagnostic, prognostic, and treatment guidance for patients with HNSCC.

Keratin gene signature expression drives epithelial-mesenchymal transition through enhanced TGF-β signaling pathway activation and correlates with adverse prognosis in lung adenocarcinoma

BackgroundLung adenocarcinoma (LUAD) stands as the foremost histological subtype of non-small-cell lung cancer, accounting for approximately 40% of all lung cancer diagnoses. However, there remains a critical unmet need to enhance the prediction of clinical outcomes and therapy responses in LUAD patients. Keratins (KRTs), serving as the structural components of the intermediate filament cytoskeleton in epithelial cells, play a crucial role in the advancement of tumor progression. This study investigated the prognostic significance of the KRT family gene and developed a KRT gene signature (KGS) for prognostic assessment and treatment guidance in LUAD. MethodsTranscriptome profiles and associated clinical details of LUAD patients were meticulously gathered from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. The KGS score was developed based on the expression of five prognostic KRT genes (KRT7, KRT8, KRT17, KRT18, and KRT80), and the upper quartile of the KGS score was chosen as the cutoff. The Kaplan-Meier method was evaluated to compare survival outcomes between KGS-high and KGS-low groups. The underlying mechanism was further investigated by GSEA, GSVA, and other bioinformatic algorithms. ResultsHigh expression of the KGS signature exhibited a robust association with poorer overall survival (OS) in the TCGA-LUAD dataset (HR: 1.81; 95% CI: 1.35–2.42, P = 0.00011). The association was further corroborated in three external GEO cohorts, including GSE31210 (HR: 3.31; 95% CI: 1.7–6.47, P = 0.00017), GSE72094 (HR: 1.95; 95% CI: 1.34–2.85, P = 0.00057) and GSE26939 (HR: 3.19; 95% CI: 1.74–5.84, P < 0.0001). Interestingly, KGS-high tumors revealed enrichments in TGF-β and WNT-β catenin signaling pathways, exhibited heightened activation of the epithelial-mesenchymal transition (EMT) pathway and proved intensified tumor stemness compared to their KGS-low counterparts. Additionally, KGS-high tumor cells exhibited increased sensitivity to several targeted agents, including gefitinib, erlotinib, lapatinib, and trametinib, in comparison to KGS-low cells. ConclusionThis study developed a KGS score that independently predicts the prognosis in LUAD. High expression of KGS score, accompanied by upregulation of TGF-β and WNT-β catenin signaling pathways, confers more aggressive EMT and tumor progression.

The pyroptosis-related signature predicts prognosis and influences the tumor immune microenvironment in dedifferentiated liposarcoma.

Dedifferentiated liposarcoma (DDL), a member of malignant mesenchymal tumors, has a high local recurrence rate and poor prognosis. Pyroptosis, a newly discovered programmed cell death, is tightly connected with the progression and outcome of tumor. The aim of this study was to explore the role of pyroptosis in DDL. We obtained the RNA sequencing data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression databases to identify different pyroptosis-related genes (PRGs) expression pattern. An unsupervised method for clustering based on PRGs was performed. Based on the result of cluster analysis, we researched clinical outcomes and immune microenvironment between clusters. The differentially expressed genes (DEGs) between the two clusters were used to develop a prognosis model by the LASSO Cox regression method, followed by the performance of functional enrichment analysis and single-sample gene set enrichment analysis. All of the above results were validated in the Gene Expression Omnibus (GEO) dataset. Forty-one differentially expressed PRGs were found between tumor and normal tissues. A consensus clustering analysis based on PRGs was conducted and classified DDL patients into two clusters. Cluster 2 showed a better outcome, higher immune scores, higher immune cells abundances, and higher expression levels in numerous immune checkpoints. DEGs between clusters were identified. A total of 5 gene signatures was built based on the DEGs and divided all DDL patients of the TCGA cohort into low-risk and high-risk groups. The low-risk group indicates greater inflammatory cell infiltration and better outcome. For external validation, the survival difference and immune landscape between the two risk groups of the GEO cohort were also significant. Receiver operating characteristic curves implied that the risk model could exert its function as an outstanding predictor in predicting DDL patients' prognoses. Our findings revealed the clinical implication and key role in tumor immunity of PRGs in DDL. The risk model is a promising predictive tool that could provide a fundamental basis for future studies and individualized immunotherapy.

Construction of the miRNA/Pyroptosis‐Related Molecular Regulatory Axis in Abdominal Aortic Aneurysm: Evidence From Transcriptome Data Combined With Multiple Machine Learning Approaches Followed by Experiment Validation

Background: Abdominal aortic aneurysm (AAA) represents a permanent and localized widening of the abdominal aorta, posing a potentially lethal risk of aortic rupture. Several recent studies have highlighted the role of pyroptosis, a pro‐inflammatory programed cell death, as critical molecular regulators in AAA occurrence, progression, and rupture. However, the potential effects of pyroptosis in AAA and its upstream microRNA (miRNA) have not been comprehensively clarified.Methods: Through a search of the gene expression omnibus (GEO) database, the expression profiles of mRNAs (GSE7084, GSE57691, and GSE98278) and miRNAs (GSE62179) and corresponding clinical features were downloaded, respectively. Expression profiles of 15 AAA and 10 normal vascular samples were consecutively collected for in vitro experimentation and subsequent analysis. Various machine learning techniques were employed to identify hub pyroptosis‐related genes (PRGs), leading to the development of a predictive model termed the PRG classifier. Quantitative real‐time‐polymerase chain reaction (qRT‐PCR), western blot (WB), and enzyme‐linked immunosorbent assay (ELISA) were used to confirm the expression of the hub PRGs. The diagnostic and predictive capabilities of the model were comprehensively evaluated in GEO and hospital cohorts. Then, the crucial immune cell infiltration and molecular pathways implicated in the initiation and rupture of AAA and their association with pyroptosis were explored. Lastly, a miRNA/hub pyroptosis‐related molecular regulatory axis was constructed using the TargetScan dataset, which was further explored through loss‐of‐function assays.Results: Differential analysis, enrichment score analysis, and principal component analysis (PCA) revealed that pyroptosis‐related molecules were significantly involved in the occurrence of AAA. Utilizing multiple machine learning algorithms, eight key PRGs (cysteinyl aspartate specific proteinase [CASP]1, infiltrating lymphocyte [IL]1B, IL18, IL6, NOD‐, LRR‐ and pyrin domain‐containing protein [NLRP]1, NLRP2, NLRP3, and tumor necrosis factor [TNF]) were integrated to establish a PRG classifier. Demonstrating robust diagnostic capabilities (area under curve [AUC] &gt; 0.90), the PRG classifier provided clinical insights across two GEO datasets and effectively differentiated small AAA from large AAA, elective stable AAA (eAAA), and ruptured AAA (rAAA), respectively. qRT‐PCR, WB, and ELISA verified the mRNA and protein expression of the hub PRGs. Notably, in hospital cohorts, a substantial positive link was unveiled between the PRG classifier and AAA risk factors (hypertension history, diastolic pressure, triglyceride levels, and aneurysm diameter). Furthermore, immune cell infiltration and functional enrichment analysis revealed significant associations of the PRG classifier/PRGs with M2 macrophage infiltration, activated dendritic cells, and enrichment scores of the cytosolic deoxyribonucleic acid (DNA) sensing pathway and tryptophan metabolism, potentially mediating AAA onset and rupture. Finally, based on 90 differentially expressed miRNAs (DEmiRNAs) and eight hub PRGs through TargetScan dataset, a hsa‐miR‐331‐3p/TNF regulatory axis was constructed, wherein upregulation of hsa‐miR‐331‐3p expression significantly reduced TNF and CASP1 protein levels.Conclusion: A predictive model (PRG classifier) incorporating eight PRGs through multiple machine learning algorithms was developed and validated. This model may stand as a potent tool for diagnosing AAA and assessing disease severity. The identification of the cytosolic DNA sensing pathway and the hsa‐miR‐331‐3p/TNF interaction axis may represent crucial targets for AAA treatment, offering deeper insights into its potential pathogenesis.

A Comprehensive Metabolism-Related Gene Signature Predicts the Survival of Patients with Acute Myeloid Leukemia.

(1) Background: Acute myeloid leukemia (AML) is a clonal malignancy with heterogeneity in genomics and clinical outcome. Metabolism reprogramming has been increasingly recognized to play an important role in the leukemogenesis and prognosis in AML. A comprehensive prognostic model based on metabolism signatures has not yet been developed. (2) Methods: We applied Cox regression analysis and the least absolute shrinkage and selection operator (LASSO) normalization to establish a metabolism-related prognostic gene signature based on glycolysis, fatty acid metabolism, and the tricarboxylic acid cycle gene signatures. The Cancer Genome Atlas-Acute Myeloid Leukemia-like (TCGA-LAML) cohort was set as the training dataset for model construction. Three independent AML cohorts (GSE37642, GSE10358, and GSE12417) combined from Gene Expression Omnibus (GEO) datasets and the Beat-AML dataset were retrieved as two validation sets to test the robustness of the model. The transcriptome data and clinic information of the cohorts were enrolled for the analysis. (3) Results: Divided by the median value of the metabolism risk score, the five-year overall survival (OS) of the high-risk and low-risk groups in the training set were 8.2% and 41.3% (p < 0.001), respectively. The five-year OS of the high-risk and low-risk groups in the combined GEO cohort were 25.5% and 37.3% (p = 0.002), respectively. In the Beat-AML cohort, the three-year OS of the high-risk and low-risk groups were 16.2% and 40.2% (p = 0.0035), respectively. The metabolism risk score showed a significantly negative association with the long-term survival of AML. Furthermore, this metabolism risk score was an independent unfavorable factor for OS by univariate analysis and multivariate analysis. (4) Conclusions: Our study constructed a comprehensive metabolism-related signature with twelve metabolism-related genes for the risk stratification and outcome prediction of AML. This novel signature might contribute to a better use of metabolism reprogramming factors as prognostic markers and provide novel insights into potential metabolism targets for AML treatment.

Prognostic signature of colorectal cancer based on uric acid-related genes

Colorectal cancer (CRC) is one of the deadliest cancers worldwide. Numerous studies have reported a correlation between uric acid (UA) level and CRC risk. Here, we investigated the role and prognostic value of UA-related genes in CRC progression. CRC-associated gene expression and clinical data were retrieved from The Cancer Genome Atlas (TCGA), and UA-related genes were identified by overlapping the TCGA and GeneCards databases. The Gene Ontology annotation, Kyoto Encyclopedia of Genes and Genomes pathway, and Molecular Signatures Database dataset were subjected to gene set enrichment analysis. A prognostic model was constructed using the univariate and multivariate COX regression and least absolute shrinkage and selection operator (LASSO) analyses and validated using the Gene Expression Omnibus cohort. Competing endogenous RNA network, CellMiner, and Human Protein Atlas were used to detect the signature of 13 UA-related genes in the prediction model. The expression of five potential UA-related genes in CRC cell lines was confirmed via qPCR. CIBERSORT was used to evaluate immune cell infiltration in the TCGA-CRC dataset. Thirteen highly prognostic UA-related genes were used to construct a prognostic model of CRC with risk score accuracy and predictive efficacy. Abundance of activated M0 macrophages, monocytes, CD8+ T cells, and natural killer cells positively correlated with the risk score. Five promising UA-related genes showed higher expression levels in CRC than in colonic cell lines. Thus, our model posits a direct relationship between UA-related genes and CRC risk, offering novel insights into diagnosis, prognosis, and treatment.