Gene Expression Omnibus Cohorts Research Articles

Immunogenic cell death-related gene landscape predicts the overall survival and immune infiltration status of ovarian cancer.

Background: Ovarian cancer (OC) is the most troubling malignant tumor of the female reproductive system. It has a low early diagnosis rate and a high tumor recurrence rate after treatment. Immunogenic cell death (ICD) is a unique form of regulated cell death that can activate the adaptive immune system through the release of DAMPs and cytokines in immunocompromised hosts and establish long-term immunologic memory. Therefore, this study aims to explore the prognostic value and underlying mechanisms of ICD-related genes in OC on the basis of characteristics. Methods: The gene expression profiles and related clinical information of OC were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. ICD-related genes were collected from the Genecards database. ICD-related prognostic genes were obtained by intersecting ICD-related genes with the OC prognostic-related genes that were analyzed in the TCGA database. Functional enrichment, genetic mutation, and immune infiltration correlation analyses were further performed to identify underlying mechanisms. Subsequently, we developed a TCGA cohort-based prognostic risk model that included a nine-gene signature through univariate and multivariate Cox regression and LASSO regression analyses. Meanwhile, external validation was performed on two sets of GEO cohorts and the TCGA training cohort for three other common tumors in women. In addition, a nomogram was established by integrating clinicopathological features and ICD-related gene signature to predict survival probability. Finally, functional enrichment and immune infiltration analyses were performed on the two risk subgroups. Results: By utilizing nine genes (ERBB2, RB1, CCR7, CD38, IFNB1, ANXA2, CXCL9, SLC9A1, and SLAMF7), we constructed an ICD-related prognostic signature. Subsequently, patients were subdivided into high- and low-risk subgroups in accordance with the median value of the risk score. In multivariate Cox regression analyses, risk score was an independent prognostic factor (hazard ratio = 2.783; p < 0.01). In the TCGA training cohort and the two GEO validation cohorts, patients with high-risk scores had worse prognosis than those with low-risk scores (p < 0.05). The time-dependent receiver operating characteristic curve further validated the prognostic power of the gene signature. Finally, gene set enrichment analysis indicated that multiple oncological pathways were significantly enriched in the high-risk subgroup. By contrast, the low-risk subgroup was strongly related to the immune-related signaling pathways. Immune infiltration analysis further illustrated that most immune cells showed higher levels of infiltration in the low-risk subgroup than in the high-risk subgroup. Conclusion: We constructed a novel ICD-related gene model for forecasting the prognosis and immune infiltration status of patients with OC. In the future, new ICD-related genes may provide novel potential targets for the therapeutic intervention of OC.

Identification and validation of a tyrosine metabolism-related prognostic prediction model and characterization of the tumor microenvironment infiltration in hepatocellular carcinoma.

Hepatocellular carcinoma (HCC) is an aggressive and heterogeneous disease characterized by high morbidity and mortality. The liver is the vital organ that participates in tyrosine catabolism, and abnormal tyrosine metabolism could cause various diseases, including HCC. Besides, the tumor immune microenvironment is involved in carcinogenesis and can influence the patients' clinical outcomes. However, the potential role of tyrosine metabolism pattern and immune molecular signature is poorly understood in HCC. Gene expression, somatic mutations, copy number variation data, and clinicopathological information of HCC were downloaded from The Cancer Genome Atlas (TCGA) database. GSE14520 from the Gene Expression Omnibus (GEO) databases was used as a validation dataset. We performed unsupervised consensus clustering of tyrosine metabolism-related genes (TRGs) and classified patients into distinct molecular subtypes. We used ESTIMATE algorithms to evaluate the immune infiltration. We then applied LASSO Cox regression to establish the TRGs risk model and validated its predictive performance. In this study, we first described the alterations of 42 TRGs in HCC cohorts and characterized the clinicopathological characteristics and tumor microenvironmental landscape of the two distinct subtypes. We then established a tyrosine metabolism-related scoring system and identified five TRGs, which were highly correlated with prognosis and representative of this gene set, namely METTL6, GSTZ1, ADH4, ADH1A, and LCMT1. Patients in the high-risk group had an inferior prognosis. Univariate and multivariate Cox proportional hazards regression analysis also showed that the tyrosine metabolism-related signature was an independent prognostic indicator. Besides, receiver operating characteristic curve (ROC) analysis demonstrated the predictive accuracy of the TRGs signature that could reliably predict 1-, 3-, and 5-year survival in both TCGA and GEO cohorts. We also got consistent results by performing clone formation and invasion analysis, and immunohistochemical (IHC) assays. Moreover, we also discovered that the TRGs signature was significantly associated with the different immune landscapes and therapeutic drug sensitivity. Our comprehensive analysis revealed the potential molecular signature and clinical utilities of TRGs in HCC. The model based on five TRGs can accurately predict the survival outcomes of HCC, improving our knowledge of TRGs in HCC and paving a new path for guiding risk stratification and treatment strategy development for HCC patients.

Identification of cuproptosis-related subtypes and development of a prognostic signature in colorectal cancer

Cuproptosis, a novel form of copper-mediated regulated cell death, participates in tumor progression. However, the role of cuproptosis-related genes (CRGs) in colorectal cancer (CRC) remains unclear. We aimed to investigate the cuproptosis subtypes and build a predictive model to improve the prognosis of patients with CRC. Gene expression data were downloaded from the TCGA database to identify distinct molecular subtypes using a non-negative matrix factorization algorithm. A robust and efficient prognostic signature was constructed by performing multivariate Cox regression analysis and further validated using the Gene Expression Omnibus cohort. Based on the gene expression matrix of CRC, the abundance of infiltrating immune cells and tumour microenvironment scores were calculated using the CIBERSORT and ESTIMATE algorithms, respectively. The pRRophetic algorithm was used to predict the sensitivity of the patients to different chemotherapy drugs. Two distinct molecular subtypes were identified based on 41 CRGs, with subtype C1 being characterized by an advanced clinical stage and worse overall survival. A prognostic signature was constructed based on the DEGs between the two cuproptosis subtypes, and its predictive ability was validated in an external database. Patients with CRC who belonged to the low-risk group had significantly higher survival rates than those who belonged to the high-risk group. Additionally, it remained a valid prognostic indicator in strata of age, sex, tumor location, and TNM stage, and its significance persisted after the multivariate Cox regression analysis. By further analyzing the prognostic signature, a higher immune score was observed in the low-risk group, which presented a better prognosis. AKT.inhibitor.VIII, doxorubicin, lenalidomide, and tipiparnib were more sensitive in the high-risk score group. A highly accurate nomogram was constructed to improve clinical application of the risk score. Compared with an ideal nomogram, our model, consisting of clinicopathological features, performed well in predicting patient survival. In conclusion, our study provides new ways and perspectives for the prediction of the prognosis of patients with CRC and guide more effective treatment regimens.

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on T-cell marker genes to predict prognosis and therapeutic response in lung squamous cell carcinoma.

Cancer immunotherapy is an increasingly successful strategy for treating patients with advanced or conventionally drug-resistant cancers. T cells have been proved to play important roles in anti-tumor and tumor microenvironment shaping, while these roles have not been explained in lung squamous cell carcinoma (LUSC). In this study, we first performed a comprehensive analysis of single-cell RNA sequencing (scRNA-seq) data from the gene expression omnibus (GEO) database to identify 72 T-cell marker genes. Subsequently, we constructed a 5-gene prognostic signature in the training cohort based on the T-cell marker genes from the cancer genome atlas (TCGA) database, which was further validated in the testing cohort and GEO cohort. The areas under the receiver operating characteristic curve at 1-, 3-, and 5-years were 0.614, 0.713 and 0.702 in the training cohort, 0.669, 0.603 and 0.645 in the testing cohort, 0.661, 0.628 and 0.590 in the GEO cohort, respectively. Furthermore, we created a highly reliable nomogram to facilitate clinical application. Gene set enrichment analysis showed that immune-related pathways were mainly enriched in the high-risk group. Tumor immune microenvironment indicated that high-risk group exhibited higher immune score, stromal score, and immune cell infiltration levels. Moreover, genes of the immune checkpoints and human leukocyte antigen family were all overexpressed in high-risk group. Drug sensitivity revealed that low-risk group was sensitive to 8 chemotherapeutic drugs and high-risk group to 4 chemotherapeutic drugs. In short, our study reveals a novel prognostic signature based on T-cell marker genes, which provides a new target and theoretical support for LUSC patients.

Comprehensive Molecular Analyses of an SLC Family-Based Model in Stomach Adenocarcinoma.

Background: Solute carrier (SLC) family members are crucial in transporting amino acids across membranes. Amino acids are indispensable for both cancer and immune cells. However, the clinical significance of amino acid transporting SLC members in stomach adenocarcinoma (STAD) remains unclear. This study aimed to develop an SLC family-based model to predict the prognosis and the response of STAD patients to immunotherapy. Methods: A total of 1239 tumor cases were obtained from online databases. The training set (n = 371) consisted of RNA sequencing profiles obtained from The Cancer Genome Atlas (TCGA), while those from Gene Expression Omnibus (GEO) were used as the test set. Subsequently, the clinical characteristics and immune profiles were investigated, and potential immunotherapy response prediction values of the model were assessed. Results: Based on the TCGA cohort, an SLC family-based model was developed using multivariate Cox analysis. All tumor cases were stratified into high- and low-risk groups considering the SLC model. High-risk patients had a worse overall survival (OS) than low-risk patients, consistent with the results of GEO cohorts. Comprehensive analyses revealed that the high-risk group was correlated with aggressiveness-related pathways, whereas the low-risk group had better T helper cell infiltration and stronger immunotherapy response. Compared to the high-risk group, the low-risk group presented increased PD-L1 and tumor mutation burden. Conclusion: This SLC family-based model has the potential to predict the prognosis and immunotherapy outcomes of STAD patients. The survival of patients in the low-risk group was greatly prolonged, and the patients may benefit more from immunotherapy.

Characterization of pyruvate metabolism and citric acid cycle patterns predicts response to immunotherapeutic and ferroptosis in gastric cancer

BackgroundGastric cancer is one of the most common malignancies of the digestive system with a high lethal rate. Studies have shown that inherited and acquired mutations in pyruvate metabolism and citric acid cycle (P-CA) enzymes are involved in tumorigenesis and tumor development. However, it is unclear how different P-CA patterns affect the tumor microenvironment (TME), which is critical for cancer progression.MethodsThis study mainly concentrated on investigating the role of the P-CA patterns in multicellular immune cell infiltration of GC TME. First, the expression levels of P-CA regulators were profiled in GC samples from The Cancer Genome Atlas and Gene Expression Omnibus cohorts to construct a consensus clustering analysis and identify three distinct P-CA clusters. GSVA was conducted to reveal the different biological processes in three P-CA clusters. Subsequently, 1127 cluster-related differentially expressed genes were identified, and prognostic-related genes were screened using univariate Cox regression analysis. A scoring system was then set up to quantify the P-CA gene signature and further evaluate the response of the patients to the immunotherapy.ResultsWe found that GC patients in the high P-CA score group had a higher tumor mutational burden, higher microsatellite instability, and better prognosis. The opposite was observed in the low P-CA score group. Interestingly, we demonstrated P-CA gene cluster could predict the sensitivity to immunotherapy and ferroptosis-induced therapy.ConclusionCollectively, the P-CA gene signature in this study exhibits potential roles in the tumor microenvironment and predicts the response to immunotherapeutic. The identification of these P-CA patterns may significantly accelerate the strategic development of immunotherapy for GC.

Correlation between Ferroptosis-Related Gene Signature and Immune Landscape, Prognosis in Breast Cancer.

Breast cancer (BC) is the most commonly diagnosed cancer and second leading cause of cancer-related death in women worldwide. Ferroptosis, an iron-dependent newly discovered mode of cell death, can be induced by lenaltinib and plays an important role in the biological behaviors of BC. Therefore, the prognostic value of ferroptosis-related genes (FRGs) in BC warrants further investigation. FRG expression profiles and clinical data were downloaded from The Cancer Genome Atlas (TCGA) database and Gene Expression Omnibus (GEO). Immune-related pathways were found in the functional analysis. Significant differences in enrichment scores for immune cells were observed. Some patients from TCGA-BRCA were included as the training cohort. A six-gene prediction signature was constructed with the least absolute shrinkage and selection operator Cox regression. This model was validated in the rest of the TCGA-BRCA and GEO cohort. The expressions of the six FRGs were verified with real-time quantitative polymerase chain reaction and immunohistochemistry in the Human Protein Atlas. Relapse or metastasis was more likely in the high-risk group. Risk score was an independent predictor of disease-free survival. Collectively, the ferroptosis-related risk model established in this study may serve as an effective tool to predict the prognosis in BC.

Comprehensive analysis of PTPN gene family revealing PTPN7 as a novel biomarker for immuno-hot tumors in breast cancer.

Background: The non-receptor protein tyrosine phosphatase (PTPN) gene family has been considered to be involved in the oncogenesis and development of multiple cancers. However, its prognostic utility and immunological relevance in breast cancer (BrCa) have not been clarified. Methods: A transcriptional level interpretation of the expressions and prognostic values was analyzed using the data from The Cancer Genome Atlas (TCGA) cohort. In addition, GO and DAVID pinpoint the functional enrichment of PTPNs. Moreover, the immune correlations of PTPN7 in BrCa and pan-cancer were further investigated based on the TCGA cohort and were testified using the in-house and the Gene Expression Omnibus (GEO) cohorts. Results: For systematic analysis of the PTPN family, we found that the expression levels of PTPN1, PTPN6, PTPN7, PTPN18, PTPN20, and PTPN22 was promoted in tumor tissues while comparing with paraneoplastic tissues during our study. We further investigated their functions and protein-protein interactions (PPI), and these results strongly suggested that PTPN family was associated with protein dephosphorylation. Next, we performed an immunological relevance analysis and found that PTPN7 was correlated with immune infiltration, suggesting a stronger association of PTPN7 with immuno-hot tumors in BrCa. In addition, results from the in-house cohort confirmed the positive correlation between PTPN7 and PD-L1. The pan-cancer analysis revealed that PTPN7 was related to PD-L1 and CTLA-4 expression in almost all cancer types. Finally, the predictive value of PTPN7 for immunotherapy was significant in two independent GEO cohorts. Conclusion: In conclusion, this is the first extensive research on the correlation between PTPN family expression and immune characterization in BrCa. As results, PTPN7 expression is associated with immuno-hot tumors and could be a promising predictive biomarker for immunotherapy in not only BrCa but multiple cancers.

Research into the biological differences and targets in lung cancer patients with diverse immunotherapy responses.

BackgroundImmunotherapy has gradually become an important therapy option for lung cancer patients.MethodsThe Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were responsible for all the public data.ResultsIn our study, we firstly identified 22 characteristic genes of NSCLC immunotherapy response using the machine learning algorithm. Molecule subtyping was then conducted and two patient subtypes were identified Cluster1 and Cluster2. Results showed that Cluster1 patients had a lower TIDE score and were more sensitive to immunotherapy in both TCGA and combined GEO cohorts. Biological enrichment analysis showed that pathways of epithelial-mesenchymal transition (EMT), apical junction, KRAS signaling, myogenesis, G2M checkpoint, E2F targets, WNT/β-catenin signaling, hedgehog signaling, hypoxia were activated in Cluster2 patients. Genomic instability between Cluster1 and Cluster2 patients was not significantly different. Interestingly, we found that female patients were more adaptable to immunotherapy. Biological enrichment revealed that compared with female patients, pathways of MYC target, G2M checkpoints, mTORC1 signaling, MYC target, E2F target, KRAS signaling, oxidative phosphorylation, mitotic spindle and P53 pathway were activated. Meanwhile, monocytes might have a potential role in affecting NSCLC immunotherapy and underlying mechanism has been explored. Finally, we found that SEC14L3 and APCDD1L were the underlying targets affecting immunotherapy, as well as patients survival.ConclusionsThese results can provide direction and guidance for future research focused on NSCLC immunotherapy.

A novel lactate metabolism-related signature predicts prognosis and tumor immune microenvironment of breast cancer.

Background: Lactate, an intermediate product of glycolysis, has become an essential regulator of tumor maintenance, development, and metastasis. Lactate can drive tumors by changing the microenvironment of tumor cells. Because of lactate’s important role in cancer, we aim to find a novel prognostic signature based on lactate metabolism-related genes (LMRGs) of breast cancer (BC). Methods: RNA-sequencing data and clinical information of BC were enrolled from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. We obtained LMRGs from the Molecular Signature Database v7.4 and articles, and then we compared candidate genes with TCGA data to get differential genes. Univariate analysis and most minor absolute shrinkage and selector operator (LASSO) Cox regression were employed to filter prognostic genes. A novel lactate metabolism-related risk signature was constructed using a multivariate Cox regression analysis. The signature was validated by time-dependent ROC curve analyses and Kaplan–Meier analyses in TCGA and GEO cohorts. Then, we further investigated in depth the function of the model’s immune microenvironment. Results: We constructed a 3-LMRG-based risk signature. Kaplan–Meier curves confirmed that high-risk score subgroups had a worse prognosis in TCGA and GEO cohorts. Then a nomogram to predict the probability of survival for BC was constructed. We also performed Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway function analysis. The function analysis showed that the lactate metabolism-related signature was significantly related to immune response. A significant correlation was observed between prognostic LMRGs and tumor mutation burden, checkpoints, and immune cell infiltration. An mRNA–miRNA network was built to identify an miR-203a-3p/LDHD/LYRM7 regulatory axis in BC. Conclusion: In conclusion, we constructed a novel 3-LMRG signature and nomogram that can be used to predict the prognosis of BC patients. In addition, the signature is closely related to the immune microenvironment, which may provide new insight into future anticancer therapies.

Peripheral immune signature resembles tumor microenvironment and predicts clinical outcomes in head and neck squamous cell carcinoma.

The contour of the tumor immune microenvironment (TIME) is very important for tumor prognostic prediction but hard to be characterized in clinical practice. It is unclear practice whether the peripheral immune signature (pIS) reflects the TIME as a feasible prognostic indicator for head and neck squamous cell carcinoma (HNSCC) patients. Here, we enrolled 599 HNSCC patients from three domestic institutes to explore the relationship between the pIS and survival. The peripheral neutrophil-to-lymphocyte ratio (pNLR) was screened out as a significant prognostic variable through multivariable COX regression analyses. An inverse correlation between pNLR and survival was found in the data of these 599 patients. Meanwhile, the bulk tumor RNA-seq data of 913 cases were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases to identify the prognosis-associated TIME features. The TIME feature was consistent to the finding of clinical data, in which high tissue NLR predicted a poor prognosis. Differentially expressed immune-related gene (DEIRG) enrichment analysis also showed a trend that the gene sets in patients with a good prognosis were enriched in lymphocyte-related functions, while those with a poor prognosis were enriched in neutrophil-related functions. At the same time, the well prediction performance of our model based on DEIRGs was verified in both TCGA and GEO cohorts. Finally, the correlation between pIS and the TIME was confirmed in a small independent cohort of 30 HNSCC patients. A positive correlation was confirmed prospectively between the pNLR and the TIME pattern in our independent cohort. Our findings provide evidence that the pNLR is a feasible prognostic signature that reflects the TIME patterns to some extent in HNSCC.

Development and validation of a prognosis prediction model based on 18 endoplasmic reticulum stress-related genes for patients with lung adenocarcinoma.

BackgroundEndoplasmic reticulum (ER) stress had a crucial impact on cell survival, proliferation, and metastasis in various cancers. However, the role of ER stress in lung adenocarcinoma remains unclear.MethodGene expression and clinical data of lung adenocarcinoma (LUAD) samples were extracted from The Cancer Genome Atlas (TCGA) and three Gene Expression Omnibus (GEO) datasets. ER stress score (ERSS) was constructed based on hub genes selected from 799 ER stress-related genes by least absolute shrinkage and selection operator (LASSO) regression. A Cox regression model, integrating ERSS and the TNM stage, was developed to predict overall survival (OS) in TCGA cohort and was validated in GEO cohorts. Gene set enrichment analysis (GSEA), single-sample GSEA (ssGSEA), and gene mutation analyses were performed to further understand the molecular features of ERSS. The tumor immune infiltration was evaluated by ESTIMATE, CIBERSORT, and xCell algorithms. The receiver operating characteristic (ROC) curves were used to evaluate the predictive value of the risk model. p< 0.05 was considered statistically significant.ResultsOne hundred fifty-seven differentially expressed genes (DEGs) were identified between tumor and para-carcinoma tissues, and 45 of them significantly correlated with OS. Next, we identified 18 hub genes and constructed ERSS by LASSO regression. Multivariate analysis demonstrated that higher ERSS (p< 0.0001, hazard ratio (HR) = 3.8, 95%CI: 2.8–5.2) and TNM stage (p< 0.0001, HR = 1.55, 95%CI: 1.34–1.8) were independent predictors for worse OS. The prediction model integrating ERSS and TNM stage performed well in TCGA cohort (area under the curve (AUC) at five years = 0.748) and three GEO cohorts (AUC at 5 years = 0.658, 0.717, and 0.739). Pathway enrichment analysis showed that ERSS significantly correlated with unfolded protein response. Meanwhile, pathways associated with the cell cycle, growth, and metabolism were significantly enriched in the high ERSS group. Patients with SMARCA4, TP53, and EGFR mutations showed significantly higher ERSS (p = 4e−04, 0.0027, and 0.035, respectively). Tissues with high ERSS exhibited significantly higher infiltration of M1 macrophages, activated dendritic cells, and lower infiltration of CD8+ T cells and B cells, which indicate an activated tumor antigen-presenting but suppressive immune response status.ConclusionWe developed and validated an ER stress-related risk model that exhibited great predictive value for OS in patients with LUAD. Our work also expanded the understanding of the role of ER stress in LUAD.

Development and Validation of a Novel Circadian Rhythm-Related Signature to Predict the Prognosis of the Patients with Hepatocellular Carcinoma.

Hepatocellular carcinoma (HCC) is one of the most important causes of cancer-related deaths and remains a major public health challenge worldwide. Considering the extensive heterogeneity of HCC, more accurate prognostic models are imperative. The circadian genes regulate the daily oscillations of key biological processes, such as nutrient metabolism in the liver. Circadian rhythm disruption has recently been recognized as an independent risk factor for cancer. In this study, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) were compared and 248 differentially expressed genes (DEGs) of the circadian rhythm were identified. HCC was classified into two subtypes based on these DEGs. The prognostic value of each circadian rhythm-associated gene (CRG) for survival was assessed by constructing a multigene signature from TCGA cohort. A 6-gene signature was created by applying the least absolute shrinkage and selection operator (LASSO) Cox regression method, and all patients in TCGA cohort were divided into high- and low-risk groups according to their risk scores. The survival rate of patients with HCC in the low-risk group was significantly higher than that in the high-risk group (p < 0.001). The patients with HCC in the Gene Expression Omnibus (GEO) cohort were also divided into two risk subgroups using the risk score of TCGA cohort, and the overall survival time (OS) was prolonged in the low-risk group (p = 0.012). Based on the clinical characteristics, the risk score was an independent predictor of OS in the patients with HCC. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that multiple metabolic pathways, cell cycle, etc., were enhanced in the high-risk group. Using the metabolic pathway single-sample gene set enrichment analysis (ssGSEA), it was found that the metabolic pathways in the high- and low-risk groups between TCGA and GEO cohorts were altered essentially in the same way. In conclusion, the circadian genes play an important role in HCC metabolic rearrangements and can be further used to predict the prognosis the patients with HCC.

Establish immune-related gene prognostic index for esophageal cancer

Background: Esophageal cancer is a tumor type with high invasiveness and low prognosis. As immunotherapy has been shown to improve the prognosis of esophageal cancer patients, we were interested in the establishment of an immune-associated gene prognostic index to effectively predict the prognosis of patients. Methods: To establish the immune-related gene prognostic index of esophageal cancer (EC), we screened 363 upregulated and 83 downregulated immune-related genes that were differentially expressed in EC compared to normal tissues. By multivariate Cox regression and weighted gene coexpression network analysis (WGCNA), we built a prognostic model based on eight immune-related genes (IRGs). We confirmed the prognostic model in both TCGA and GEO cohorts and found that the low-risk group had better overall survival than the high-risk group. Results: In this study, we identified 363 upregulated IRGs and 83 downregulated IRGs. Next, we found a prognostic model that was constructed with eight IRGs (OSM, CEACAM8, HSPA6, HSP90AB1, PCSK2, PLXNA1, TRIB2, and HMGB3) by multivariate Cox regression analysis and WGCNA. According to the Kaplan–Meier survival analysis results, the model we constructed can predict the prognosis of patients with esophageal cancer. This result can be verified by the Gene Expression Omnibus (GEO). Patients were divided into two groups with different outcomes. IRGPI-low patients had better overall survival than IRGPI-high patients. Conclusion: Our findings indicated the potential value of the IRGPI risk model for predicting the prognosis of EC patients.

An immune-related gene prognostic risk index for pancreatic adenocarcinoma.

ObjectiveOur goal is to construct an immune-related gene prognostic risk index (IRGPRI) for pancreatic adenocarcinoma (PAAD), and to clarify the immune and molecular features in IRGPRI-defined PAAD subgroups and the benefit of immune checkpoint inhibitors (ICIs) therapy.MethodThrough differential gene expression analysis, weighted gene co-expression network analysis (WGCNA), and univariate Cox regression analysis, 16 immune-related hub genes were identified using the Cancer Genome Atlas (TCGA) PAAD dataset (n = 182) and immune gene set. From these genes, we constructed an IRGPRI with the Cox regression method and the IRGPRI was verified based on the Gene Expression Omnibus (GEO) dataset (n = 45). Then, we analyzed the immune and molecular features and the benefit of ICI therapy in IRGPRI-defined subgroups.ResultsFive genes, including S100A16, CD40, VCAM1, TNFRSF4 and TRAF1 were used to construct IRGPRI. As with the results of the GEO cohort, the overall survival (OS) was more favorable in low IRGPRI patients versus high IRGPRI patients. The composite results pointed out that low IRGPRI was associated with immune response-related pathways, high level of CTLA4, low KRAS and TP53 mutation rate, more infiltration of activated memory CD4+ T cells, CD8+ T cells, and more benefits from ICIs therapy. In comparison, high IRGPRI was associated with cancer-related pathways, low expression of CTLA4, high KRAS and TP53 mutation rate, more infiltration of M2 macrophages, and less benefit from ICIs therapies.ConclusionThis IRGPRI is an encouraging biomarker to define the prognosis, immune and molecular features, and benefits from ICIs treatments in PAAD.

Establishment of a prognostic signature for lung adenocarcinoma by integration of 7 pyroptosis-related genes and cross-validation between the TCGA and GEO cohorts: A comprehensive bioinformatics analysis.

Pyroptosis-related genes (PRGs) have been reported to be associated with prognosis of lung adenocarcinoma (LUAD). Until now, the relationship of PRGs to the prognosis of LUAD patients and its underlying mechanisms have been poorly elucidated. Using The Cancer Genome Atlas (TCGA) LUAD cohort, a prior bioinformatics analysis constructed a prognostic signature incorporating 5 PRGs (NLRP7, NLRP1, NLRP2, NOD1, and CASP6) for predicting prognosis of LUAD patients. However, it has not been validated by the Gene Expression Omnibus (GEO) LUAD cohort yet. We implemented a modified bioinformatics analysis to, respectively, construct one prognostic signature with the TCGA cohort and with the GEO cohort and attempted to perform cross-validations by the GEO cohort and the TCGA cohort alternately in turn. Univariate and multivariate Cox regression analysis screened PRGs and constructed 2 prognostic signatures with the TCGA and GEO cohorts. All LUAD samples were classified into high- and low-risk groups according to the median risk score that was generated by regression formula. Kaplan-Meier survival analysis compared the overall survival rate between the 2 risk groups, and receiver operating characteristic curve analysis evaluated predictive performance of the 2 signatures. Additionally, risk score, combined with clinicopathological features, was subjected to multivariate Cox regression analysis, to evaluate independent prognostic value of the 2 signatures. Finally, the 2 signatures received cross-validations by the GEO and TCGA cohorts, alternately. The TCGA cohort yielded a 3-gene signature (PYCARD, NLRP1, and NLRC4), whereas the GEO cohort built a 7-gene signature (SCAF11, NOD1, NLRP2, NLRP1, GPX4, CASP8, and AIM2) for predicting the prognosis of LUAD patients. Multivariate analysis proved independent prognostic value of risk score in the TCGA cohort (hazard ratio, = 1.939,; P = 8.43 × 10−4) and the GEO cohort (hazard ratio, = 2.291,; P = 4.34 × 10−9). Cross-validations confirmed prognostic value for the 7-gene signature from the GEO cohort by the TCGA cohort but not for the 3-gene signature from the TCGA cohort by the GEO cohort. We develop and validate a 7-gene prognostic signature (SCAF11, NOD1, NLRP2, NLRP1, GPX4, CASP8, and AIM2) with independent prognostic value for patients with LUAD.

Immunometabolism characteristics and a potential prognostic risk model associated with TP53 mutations in breast cancer.

TP53, a gene with high-frequency mutations, plays an important role in breast cancer (BC) development through metabolic regulation, but the relationship between TP53 mutation and metabolism in BC remains to be explored. Our study included 1,066 BC samples from The Cancer Genome Atlas (TCGA) database, 415 BC cases from the Gene Expression Omnibus (GEO) database, and two immunotherapy cohorts. We identified 92 metabolic genes associated with TP53 mutations by differential expression analysis between TP53 mutant and wild-type groups. Univariate Cox analysis was performed to evaluate the prognostic effects of 24 TP53 mutation-related metabolic genes. By unsupervised clustering and other bioinformatics methods, the survival differences and immunometabolism characteristics of the distinct clusters were illustrated. In a training set from TCGA cohort, we employed the least absolute shrinkage and selection operator (LASSO) regression method to construct a metabolic gene prognostic model associated with TP53 mutations, and the GEO cohort served as an external validation set. Based on bioinformatics, the connections between risk score and survival prognosis, tumor microenvironment (TME), immunotherapy response, metabolic activity, clinical characteristics, and gene characteristics were further analyzed. It is imperative to note that our model is a powerful and robust prognosis factor in comparison to other traditional clinical features and also has high accuracy and clinical usefulness validated by receiver operating characteristic (ROC) and decision curve analysis (DCA). Our findings deepen our understanding of the immune and metabolic characteristics underlying the TP53 mutant metabolic gene profile in BC, laying a foundation for the exploration of potential therapies targeting metabolic pathways. In addition, our model has promising predictive value in the prognosis of BC.

Ferroptosis-Related lncRNA Signature Correlates with the Prognosis, Tumor Microenvironment, and Therapeutic Sensitivity of Esophageal Squamous Cell Carcinoma.

Esophageal squamous cell carcinoma (ESCC) is the most prevalent form of esophageal cancer in China and is closely associated with malignant biological characteristics and poor survival. Ferroptosis is a newly discovered iron-dependent mode of cell death that plays an important role in the biological behavior of ESCC cells. The clinical significance of ferroptosis-related long noncoding RNAs (FRLs) in ESCC remains unknown and warrants further research. The current study obtained RNA sequencing profiles and corresponding clinical data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, and FRLs were obtained through coexpression analysis. Consensus clustering was employed to divide the subjects into clusters, and immune-associated pathways were identified by functional analysis. The current study observed significant differences in the enrichment scores of immune cells among different clusters. Patients from TCGA-ESCC database were designated as the training cohort. A ten-FRL prediction signature was established using the least absolute shrinkage and selection operator Cox regression model and validated using the GEO cohort and our own independent validation database. Real-time quantitative polymerase chain reaction was used to verify the expression of the ten FRLs, and the ssGSEA analysis was employed to evaluate their function. In addition, the IMvigor database was used to assess the predictive value of the signature in terms of immunotherapeutic responses. Multivariate Cox and stratification analyses revealed that the ten-FRL signature was an independent predictor of the overall survival (OS). Patients with ESCC in the high-risk group displayed worse survival, a characteristic tumor immune microenvironment, and low immunotherapeutic benefits compared to those in the low-risk group. Collectively, the risk model established in this study could serve as a promising predictor of prognosis and immunotherapeutic response in patients with ESCC.

Identification of a Seven-Differentially Expressed Gene-Based Recurrence-Free Survival Model for Melanoma Patients.

Melanoma is a malignant tumor that originates in melanocytes of the skin or mucous membrane, which has a high mortality rate and worse prognosis. Therefore, perspective prognosis evaluation seems more important for patients' treatment. Gene expression profiles of melanoma were downloaded from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases, respectively. 130 consistent differentially expressed genes (DEGs) were identified between melanoma and nevus tissues from two GEO cohorts. Prognostic genes were identified by univariate analysis, and 20 of them were regarded to be associated with the recurrence-free survival (RFS) of melanoma patients. Then, the LASSO Cox regression analysis chose seven of them to establish a seven-DEG-based RFS predicting signature. We demonstrated that this model was more powerful to predict RFS risk than other individual clinical features and was able to independently predict the RFS outcomes in different subsets of patients. We attempted to search for the underlying mechanisms by analyzing the coexpression genes of the seven candidates, and the pathway enrichment analyses indicated that immune response-related pathways might play a critical role in melanoma progression. Finally, we establish a robust seven-DEG-based RFS predicting signature, which will facilitate the personalized treatment of melanoma patients.

Lactate Metabolism-Associated lncRNA Pairs: A Prognostic Signature to Reveal the Immunological Landscape and Mediate Therapeutic Response in Patients With Colon Adenocarcinoma.

BackgroundLactate metabolism is critically involved in the tumor microenvironment (TME), as well as cancer progression. It is important to note, however, that lactate metabolism-related long non-coding RNAs (laRlncRNAs) remain incredibly understudied in colon adenocarcinoma (COAD).MethodsA gene expression profile was obtained from the Cancer Genome Atlas (TCGA) database to identify laRlncRNA expression in COAD patients. A risk signature with prognostic value was identified from TCGA and Gene Expression Omnibus (GEO) cohort based on laRlncRNA pairs by the least absolute shrinkage and selection operator (LASSO) and Cox regression analyses. Quantitative real-time polymerase chain reaction (qRT-PCR) and functional experiments were carried out to verify the expression of laRlncRNAs in COAD. The relationship of laRlncRNA pairs with immune landscape as well as the sensitivity of different therapies was explored.ResultsIn total, 2378 laRlncRNAs were identified, 1,120 pairs of which were studied to determine their prognostic validity, followed by a risk signature established based on the screened 5 laRlncRNA pairs. The laRlncRNA pairs-based signature provided a better overall survival (OS) prediction than other published signatures and functioned as a prognostic marker for COAD patients. According to the calculated optimal cut-off point, patients were divided into high- and low-risk groups. The OS of COAD patients in the high-risk group were significantly shorter than that of those in the low-risk group (P=4.252e-14 in the TCGA cohort and P=2.865-02 in the GEO cohort). Furthermore, it remained an effective predictor of survival in strata of gender, age, TNM stage, and its significance persisted after univariate and multivariate Cox regressions. Additionally, the risk signature was significantly correlated with immune cells infiltration, tumor mutation burden (TMB), microsatellite instability (MSI) as well as immunotherapeutic efficacy and chemotherapy sensitivity. Finally, one of the laRlncRNA, LINC01315, promotes proliferation and migration capacities of colon cancer cells.ConclusionThe newly identified laRlncRNAs pairs-based signature exhibits potential effects in predicting prognosis, deciphering patients’ immune landscape, and mediating sensitivity to immunotherapy and chemotherapy. Findings in our study may provide evidence for the role of laRlncRNAs pairs as novel prognostic biomarkers and potentially individualized therapy targets for COAD patients.