Development and validation of stable ferroptosis- and pyroptosis-related signatures in predicting prognosis and immune status in breast cancer

To develop and validate the predictive effects of stable ferroptosis- and pyroptosis-related features on the prognosis and immune status of breast cancer (BC). We retrieved as well as downloaded ferroptosis- and pyroptosis-related genes from the FerrDb and GeneCards databases. The minimum absolute contraction and selection operator (LASSO) algorithm in The Cancer Genome Atlas (TCGA) was used to construct a prognostic classifier combining the above two types of prognostic genes with differential expression, and the Integrated Gene Expression (GEO) dataset was used for validation. Seventeen genes presented a close association with BC prognosis. Thirteen key prognostic genes with prognostic value were considered to construct a new expression signature for classifying patients with BC into high- and low-risk groups. Kaplan–Meier analysis revealed a worse prognosis in the high-risk group. The receiver operating characteristic (ROC) curve and multivariate Cox regression analysis identified its predictive and independent features. Immune profile analysis showed that immunosuppressive cells were upregulated in the high-risk group, and this risk model was related to immunosuppressive molecules. We successfully constructed combined features of ferroptosis and pyroptosis in BC that are closely related to prognosis, clinicopathological and immune features, chemotherapy efficacy and immunosuppressive molecules. However, further experimental studies are required to verify these findings.     

丝氨酸蛋白酶抑制剂(SERPIN)家族相关基因在结肠腺癌中预后模型的建立及应用

应用生物信息学方法,构建结肠腺癌(COAD)丝氨酸蛋白酶抑制剂(SERPIN)家族相关基因预后模型。从TCGA数据库和GEO数据库下载结肠腺癌(COAD)转录组和临床数据,根据数据中SERPINs家族基因的表达量对COAD患者进行一致性聚类分析;将数据随机均分为训练集(Train)组和验证集(Test)组,基于两个亚型的差异基因,利用Train组进行COX回归和Lasso回归构建预后模型,根据模型风险评分中位值将样本分为高、低风险两组,绘制高低风险组患者生存曲线;通过ROC曲线评价模型预测能力;利用Test组数据验证模型;构建列线图,评估患者生存率模型预测值与实际值的一致性;并利用利用ESTIMATE算法和CIBERSORT算法评估风险评分和肿瘤微环境(TME)以及免疫浸润的相关性。通过34个SERPIN基因确定了两个亚型,基于2个亚型筛选出了436个预后相关分型差异基因,通过Lasso回归确定出了11个预后相关基因参与风险模型的构建,根据模型评分区分的高低风险组具有明显的生存差异,列线图可以准确预测1、3和5年生存率。肿瘤微环境分析和免疫浸润分析显示高风险评分组患者免疫活性差。SERPIN家族相关基因构建的风险评分模型能够预测COAD的预后,有利于进一步指导临床对COAD的诊治,提高患者生存率。     

Developing a risk scoring system based on immune-related lncRNAs for patients with gastric cancer

The immune system and the tumor interact closely during tumor development. Aberrantly expressed long non-coding RNAs (lncRNAs) may be potentially applied as diagnostic and prognostic markers for gastric cancer (GC). At present, the diagnosis and treatment of GC patients remain a formidable clinical challenge. The present study aimed to build a risk scoring system to improve the prognosis of GC patients. In the present study, ssGSEA was used to evaluate the infiltration of immune cells in GC tumor tissue samples, and the samples were split into a high immune cell infiltration group and a low immune cell infiltration group. About 1262 differentially expressed lncRNAs between the high immune cell infiltration group and the low immune cell infiltration group. About 3204 differentially expressed lncRNAs between GC tumor tissues and paracancerous tissues were identified. Then, 621 immune-related lncRNAs were screened using a Venn analysis based on the above results, and 85 prognostic lncRNAs were identified using a univariate Cox analysis. We constructed a prognostic signature using LASSO analysis and evaluated the predictive performance of the signature using ROC analysis. GO and KEGG enrichment analyses were performed on the lncRNAs using the R package, ‘clusterProfiler’. The TIMER online database was used to analyze correlations between the risk score and the abundances of the six types of immune cells. In conclusion, our study found that specific immune-related lncRNAs were clinically significant. These lncRNAs were used to construct a reliable prognostic signature and analyzed immune infiltrates, which may assist clinicians in developing individualized treatment strategies for GC patients.     

Exploration of the prognostic prediction value of the PANoptosis-based risk score and its correlation with tumor immunity in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a common cancer with high mortality worldwide. PANoptosis is a novel inflammatory programmed cell death modality with the characteristics of pyroptosis, apoptosis and necroptosis. It is necessary to explore PANoptosis-related genes in LUAD patients and offer evidence for prognosis prediction and therapeutic strategies. Single-cell RNA sequencing data and RNA expression profiles of LUAD patients from The Cancer Genome Atlas and Gene Expression Omnibus databases are used to screen PANoptosis-related differential genes for the construction of a risk model. Fifteen PANoptosis-related markers with prognostic value were identified by Least Absolute Shrinkage and Selection Operator (LASSO)–Cox regression analysis. Kaplan–Meier analysis and receiver operating characteristic curve analysis further demonstrated the significant predictive capability. Immune infiltration, Single Nucleotide Variants (SNV) mutations, and clinical drug susceptibility were analyzed. In conclusion, a risk model of 15 PANoptosis-related genes has significant value in prognostic prediction for LUAD and has potential to direct clinical therapeutic strategies during the treatment.     

Identification of an m6A-related lncRNA prognostic signature in colorectal cancer

Data sets of colorectal cancer (CRC) were obtained from The Cancer Genome Atlas (TCGA), three N6-methyladenosine (m6A) subtypes were identified using 21 m6A-related long noncoding RNAs (lncRNAs) and differential m6A subtypes of different CRC tumors were determined in this study to evaluate the m6A expression and the prognosis of patients with CRC. Subsequently, eight key lncRNAs were identified based on co-expression with 21 m6A-related genes in CRC tumors using the single-factor Cox and least absolute shrinkage and selection operator. Finally, an m6A-related lncRNA risk score model of CRC tumor was established using multifactor Cox regression based on the eight important lncRNAs and found to have a better performance in evaluating the prognosis of patients in the TCGA-CRC data set. TCGA-CRC tumor samples were divided based on the risk scores: high and low. Then, the clinical characteristics, tumor mutation load, and tumor immune cell infiltration difference between the high- and low-risk-score groups were explored, and the predictive ability of the risk score was assessed for immunotherapeutic benefits. We found that the risk score model can determine the overall survival, be a relatively independent prognostic indicator, and better evaluate the immunotherapeutic benefits for patients with CRC. This study provides data support for accurate immunotherapy in CRC.     

An outcome model for human bladder cancer: A comprehensive study based on weighted gene co‐expression network analysis

The precision evaluation of prognosis is crucial for clinical treatment decision of bladder cancer (BCa). Therefore, establishing an effective prognostic model for BCa has significant clinical implications. We performed WGCNA and DEG screening to initially identify the candidate genes. The candidate genes were applied to construct a LASSO Cox regression analysis model. The effectiveness and accuracy of the prognostic model were tested by internal/external validation and pan‐cancer validation and time‐dependent ROC. Additionally, a nomogram based on the parameter selected from univariate and multivariate cox regression analysis was constructed. Eight genes were eventually screened out as progression‐related differentially expressed candidates in BCa. LASSO Cox regression analysis identified 3 genes to build up the outcome model in E‐MTAB‐4321 and the outcome model had good performance in predicting patient progress free survival of BCa patients in discovery and test set. Subsequently, another three datasets also have a good predictive value for BCa patients' OS and DFS. Time‐dependent ROC indicated an ideal predictive accuracy of the outcome model. Meanwhile, the nomogram showed a good performance and clinical utility. In addition, the prognostic model also exhibits good performance in pan‐cancer patients. Our outcome model was the first prognosis model for human bladder cancer progression prediction via integrative bioinformatics analysis, which may aid in clinical decision‐making.     

Hub gene identification and prognostic model construction for isocitrate dehydrogenase mutation in glioma

Our study attempted to identify hub genes related to isocitrate dehydrogenase (IDH) mutation in glioma and develop a prognostic model for IDH-mutant glioma patients. In a first step, ten hub genes significantly associated with the IDH status were identified by weighted gene coexpression analysis (WGCNA). The functional enrichment analysis demonstrated that the most enriched terms of these hub genes were cadherin binding and glutathione metabolism. Three of these hub genes were significantly linked with the survival of glioma patients. 328 samples of IDH-mutant glioma were separated into two datasets: a training set (N = 228) and a test set (N = 100). Based on the training set, we identified two IDH-mutant subtypes with significantly different pathological features by using consensus clustering. A 31 gene-signature was identified by the least absolute shrinkage and selection operator (LASSO) algorithm and used for establishing a differential prognostic model for IDH-mutant patients. In addition, the test set was employed for validating the prognostic model, and the model was proven to be of high value in classifying prognostic information of samples. The functional annotation revealed that the genes related to the model were mainly enriched in nuclear division, DNA replication, and cell cycle. Collectively, this study provided novel insights into the molecular mechanism of IDH mutation in glioma, and constructed a prognostic model which can be effective for predicting prognosis of glioma patients with IDH-mutation, which might promote the development of IDH target agents in glioma therapies and contribute to accurate prognostication and management in IDH-mutant glioma patients.     

Novel prognostic model established for patients with head and neck squamous cell carcinoma based on pyroptosis-related genes

We aimed at establishing a risk – score model using pyroptosis-related genes to predict the prognosis of patients with head and neck squamous cell carcinoma (HNSCC). A total of 33 pyroptosis-related genes were selected. We then evaluated the data of 502 HNSCC patients and 44 normal patients from TCGA database. Gene expression was then profiled to detect differentially expressed genes (DEGs). Using the univariate, the least absolute shrinkage and selection operator (LASSO) Cox regression analyses, we generated a risk – score model. Tissue samples from neoplastic and normal sites of 44 HNSCC patients were collected. qRT-PCR were employed to analyze the mRNA level of the samples. Kaplan-Meier method was used to evaluate the overall survival rate (OS). Enrichment analysis was performed to elucidate the underlying mechanism of HNSCC patient's differentially survival status from the perspective of tumor immunology. 17 genes were categorized as DEGs. GSDME, IL-6, CASP8, CASP6, NLRP1 and NLRP6 were used to establish the risk – score model. Each patient's risk score in the TCGA cohort was calculated using the risk – score formula. The risk score was able to independently predict the OS of the HNSCC patients (P = 0.02). The OS analysis showed that the risk score model (P < 0.0001) was more reliable than single gene, a phenomenon verified by practical patient cohort. Additionally, enrichment analysis indicated more active immune activities in low-risk group than high-risk group. In conclusion, our risk – score model has provided novel strategy for the prediction of HNSCC patients’ prognosis.     

The landscape and prognostic value of immune characteristics in uterine corpus endometrial cancer

In the present study, we explored the clinical and immunological characteristics of 575 uterine corpus endometrial carcinoma (UCEC) samples obtained from The Cancer Genome Atlas (TCGA) using the ESTIMATE and CIBERSORT algorithms. First, Kaplan–Meier and univariate Cox regression analyses indicated that the immune cell score was a prognostic factor for overall survival (OS) and recurrence-free survival (RFS). Multivariate Cox regression analysis further revealed that the immune cell score was an independent prognostic factor for UCEC patients. Second, we investigated the correlation between the infiltration levels of 22 types of immune cells and the immune score. Survival analysis based on the 22 immune cell types showed that higher levels of regulatory T cell, activated NK cell, and follicular helper T-cell infiltration were associated with longer OS, while higher levels of CD8+ T cell and naive B-cell infiltration were associated with longer RFS. Next, we performed differential expression and prognosis analyses on 1534 immune-related genes and selected five from 14 candidate genes to construct a prognostic prediction model. The area under the receiver-operating characteristic (ROC) curve (AUC) for 3- and 5-year survival were 0.711 and 0.728, respectively. Further validation using a stage I–II subgroup showed similar results, presenting AUC values for 3- and five-year survival of 0.677 and 0.692, respectively. Taken together, the present study provides not only a deeper understanding of the relationship between UCEC and the immune landscape but also guidance for the future development of UCEC immunotherapy.     

Novel prognostic model established for patients with head and neck squamous cell carcinoma based on pyroptosis-related genes

We aimed at establishing a risk – score model using pyroptosis-related genes to predict the prognosis of patients with head and neck squamous cell carcinoma (HNSCC). A total of 33 pyroptosis-related genes were selected. We then evaluated the data of 502 HNSCC patients and 44 normal patients from TCGA database. Gene expression was then profiled to detect differentially expressed genes (DEGs). Using the univariate, the least absolute shrinkage and selection operator (LASSO) Cox regression analyses, we generated a risk – score model. Tissue samples from neoplastic and normal sites of 44 HNSCC patients were collected. qRT-PCR were employed to analyze the mRNA level of the samples. Kaplan-Meier method was used to evaluate the overall survival rate (OS). Enrichment analysis was performed to elucidate the underlying mechanism of HNSCC patient's differentially survival status from the perspective of tumor immunology. 17 genes were categorized as DEGs. GSDME, IL-6, CASP8, CASP6, NLRP1 and NLRP6 were used to establish the risk – score model. Each patient's risk score in the TCGA cohort was calculated using the risk – score formula. The risk score was able to independently predict the OS of the HNSCC patients (P = 0.02). The OS analysis showed that the risk score model (P < 0.0001) was more reliable than single gene, a phenomenon verified by practical patient cohort. Additionally, enrichment analysis indicated more active immune activities in low-risk group than high-risk group. In conclusion, our risk – score model has provided novel strategy for the prediction of HNSCC patients’ prognosis.     

A novel DNA methylation-based model that effectively predicts prognosis in hepatocellular carcinoma

Purpose: To build a novel predictive model for hepatocellular carcinoma (HCC) patients based on DNA methylation data.Methods: Four independent DNA methylation datasets for HCC were used to screen for common differentially methylated genes (CDMGs). Gene Ontology (GO) enrichment, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were used to explore the biological roles of CDMGs in HCC. Univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) Cox analysis were performed to identify survival-related CDMGs (SR-CDMGs) and to build a predictive model. The importance of this model was assessed using Cox regression analysis, propensity score-matched (PSM) analysis and stratification analysis. A validation group from the Cancer Genome Atlas (TCGA) was constructed to further validate the model.Results: Four SR-CDMGs were identified and used to build the predictive model. The risk score of this model was calculated as follows: risk score = (0.01489826 × methylation level of WDR69) + (0.15868618 × methylation level of HOXB4) + (0.16674959 × methylation level of CDKL2) + (0.16689301 × methylation level of HOXA10). Kaplan–Meier analysis demonstrated that patients in the low-risk group had a significantly longer overall survival (OS; log-rank P-value =0.00071). The Cox model multivariate analysis and PSM analysis identified the risk score as an independent prognostic factor (P<0.05). Stratified analysis results further confirmed this model performed well. By analyzing the validation group, the results of receiver operating characteristic (ROC) curve analysis and survival analysis further validated this model.Conclusion: Our DNA methylation-based prognosis predictive model is effective and reliable in predicting prognosis for patients with HCC.     

Identification and validation of autophagy-related prognostic signature for head and neck squamous cell carcinoma

BackgroundMany studies have demonstrated that autophagy plays a significant role in regulating tumor growth and progression. However, the effect of autophagy-related genes (ARGs) on the prognosis have rarely been analyzed in head and neck squamous cell carcinoma (HNSCC).MethodsWe obtained differentially expressed ARGs from HNSCC mRNA data in The Cancer Genome Atlas (TCGA) database. And then we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to explore the autophagy-related biological functions. The overall survival (OS)-related and disease specific survival (DSS)-related ARGs were identified by univariate Cox regression analyses. With these genes, we established OS-related and DSS-related risk signature by LASSO regression method, respectively. We validated the reliability of the risk signature with receiver operating characteristic (ROC) analysis, Kaplan-Meier survival curves, clinical correlation analysis, and nomogram. Then we analyzed relationships between risk signature and immune cell infiltration.ResultsWe established the prognostic signatures based on 14 ARGs for OS and 12 ARGs for DSS. The ROC curves, survival analysis, and nomogram validated the predictive accuracy of the models. Clinic correlation analysis showed that the risk group was closely related to Stage, pathological T stage, pathological N stage and human papilloma virus (HPV) subtype. Cox regression demonstrated that the risk score was an independent predictor for the prognosis of HNSCC patients. Furthermore, patients in low-risk score group exhibited higher immunescore and distinct immune cell infiltration than high-risk score group. And we further analysis revealed that the copy number alterations (CNAs) of ARGs-based signature affected the abundance of tumor-infiltrating immune cells.ConclusionIn this study, we identified novel autophagy-related signature for the prediction of OS and DSS in patients with HNSCC. Meanwhile, our study provides a novel sight to understand the role of autophagy and elucidate the important role of autophagy in tumor immune microenvironment (TIME) of HNSCC.     

Pyroptosis‐related molecular classification and immune microenvironment infiltration in breast cancer: A novel therapeutic target

The underlying role of pyroptosis in breast cancer (BC) remains unknown. Herein, we investigated the correlations of 33 pyroptosis‐related genes (PRGs) with immune checkpoints and immune cell infiltrations in BC patients based on The Cancer Genome Atlas cohort (n = 996) and Gene Expression Omnibus cohort (n = 3,262). Enrichment analysis revealed that these PRGs mainly functioned in pyroptosis, inflammasomes and regulation of autophagy pathway. Four prognostic independent PRGs (CASP9, TIRAP, GSDMC and IL18) were identified. Then, cluster 1/2 was recognized using consensus clustering for these four PRGs. Patients from cluster 1 had a favourable prognosis and diverse immune cell infiltrations. A nomogram was developed based on age, TNM stage, tumour subtype and pyroptosis score. Patients with the high‐risk group exhibited worse 5‐year OS, and the result was consistent in the external cohort. Additionally, high‐risk group patients were associated with downregulated immune checkpoint expression. Further analysis suggested that the high‐risk group patients were associated with a higher IC50 of paclitaxel, doxorubicin, cisplatin, methotrexate and vinorelbine. In summarizing, the pyroptosis score‐based nomogram might serve as an independent prognostic predictor and could guide medication for chemotherapy. Additionally, it may bring novel insight into the regulation of tumour immune microenvironment in BC and help to achieve precision immunotherapy.     

Construction of a five-gene prognostic model based on immune-related genes for the prediction of survival in pancreatic cancer

Purpose: To identify differentially expressed immune-related genes (DEIRGs) and construct a model with survival-related DEIRGs for evaluating the prognosis of patients with pancreatic cancer (PC).Methods: Six microarray gene expression datasets of PC from the Gene Expression Omnibus (GEO) and Immunology Database and Analysis Portal (ImmPort) were used to identify DEIRGs. RNA sequencing and clinical data from The Cancer Genome Atlas Program-Pancreatic Adenocarcinoma (TCGA-PAAD) database were used to establish the prognostic model. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analyses were applied to determine the final variables of the prognostic model. The median risk score was used as the cut-off value to classify samples into low- and high-risk groups. The prognostic model was further validated using an internal validation set of TCGA and an external validation set of {"type":"entrez-geo","attrs":{"text":"GSE62452","term_id":"62452"}}GSE62452.Results: In total, 142 DEIRGs were identified from six GEO datasets, 47 were survival-related DEIRGs. A prognostic model comprising five genes (i.e., ERAP2, CXCL9, AREG, DKK1, and IL20RB) was established. High-risk patients had poor survival compared with low-risk patients. The 1-, 2-, 3-year area under the receiver operating characteristic (ROC) curve of the model reached 0.85, 0.87, and 0.93, respectively. Additionally, the prognostic model reflected the infiltration of neutrophils and dendritic cells. The expression of most characteristic immune checkpoints was significantly higher in the high-risk group versus the low-risk group.Conclusions: The five-gene prognostic model showed reliably predictive accuracy. This model may provide useful information for immunotherapy and facilitate personalized monitoring for patients with PC.     

A six-gene-based signature for breast cancer radiotherapy sensitivity estimation

Breast cancer (BRCA) represents the most common malignancy among women worldwide with high mortality. Radiotherapy is a prevalent therapeutic for BRCA that with heterogeneous effectiveness among patients. Here, we proposed to develop a gene expression-based signature for BRCA radiotherapy sensitivity estimation. Gene expression profiles of BRCA samples from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) were obtained and used as training and independent testing dataset, respectively. Differential expression genes (DEGs) in BRCA samples compared with their paracancerous samples in the training set were identified by using the edgeR Bioconductor package. Univariate Cox regression analysis and LASSO Cox regression method were applied to screen optimal genes for constructing a radiotherapy sensitivity estimation signature. Nomogram combining independent prognostic factors was used to predict 1-, 3-, and 5-year OS of radiation-treated BRCA patients. Relative proportions of tumor infiltrating immune cells (TIICs) calculated by CIBERSORT and mRNA levels of key immune checkpoint receptors was adopted to explore the relation between the signature and tumor immune response. As a result, 603 DEGs were obtained in BRCA tumor samples, six of which were retained and used to construct the radiotherapy sensitivity prediction model. The signature was proved to be robust in both training and testing sets. In addition, the signature was closely related to the immune microenvironment of BRCA in the context of TIICs and immune checkpoint receptors’ mRNA levels. In conclusion, the present study obtained a radiotherapy sensitivity estimation signature for BRCA, which should shed new light in clinical and experimental research.     

Molecular pathways differentiate hepatitis C virus (HCV) recurrence from acute cellular rejection in HCV liver recipients

Acute cellular rejection (ACR) and hepatitis C virus (HCV) recurrence (HCVrec) are common complications after liver transplantation (LT) in HCV patients, who share common clinical and histological features, making a differential diagnosis difficult. Fifty-three liver allograft samples from unique HCV LT recipients were studied using microarrays, including a training set (n = 32) and a validation set (n = 19). Two no-HCV-ACR samples from LT recipients were also included. Probe set intensity values were obtained using the robust multiarray average method (RMA) method. Analysis of variance identified statistically differentially expressed genes (P ≤ 0.005). The limma package was used to fit the mixed-effects models using a restricted maximum likelihood procedure. The last absolute shrinkage and selection operator (LASSO) model was fit with HCVrec versus ACR as the dependent variable predicted. N-fold cross-validation was performed to provide an unbiased estimate of generalization error. A total of 179 probe sets were differentially expressed among groups, with 71 exclusive genes between HCVrec and HCV-ACR. No differences were found within ACR group (HCV-ACR vs. no-HCV-ACR). Supervised clustering analysis displayed two clearly independent groups, and no-HCV-ACR clustered within HCV-ACR. HCVrec-related genes were associated with a cytotoxic T-cell profile, and HCV-ACR-related genes were associated with the inflammatory response. The best-fitting LASSO model classifier accuracy, including 15 genes, has an accuracy of 100% in the training set. N-fold cross-validation accuracy was 78.1%, and sensitivity, specificity and positive and negative predictive values were 50.0%, 90.9%, 71.4% and 80.0%, respectively. Arginase type II (ARG2), ethylmalonic encephalopathy 1 (ETHE1), transmembrane protein 176A (TMEM176A) and TMEM176B genes were significantly confirmed in the validation set. A molecular signature capable of distinguishing HCVrec and ACR in HCV LT recipients was identified and validated.