scholarly journals Systemic characterization of alternative splicing related to prognosis and immune infiltration in malignant mesothelioma

BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Jinzhi Lai ◽  
Hainan Yang ◽  
Tianwen Xu

Abstract Background Malignant mesothelioma (MM) is a relatively rare and highly lethal tumor with few treatment options. Thus, it is important to identify prognostic markers that can help clinicians diagnose mesothelioma earlier and assess disease activity more accurately. Alternative splicing (AS) events have been recognized as critical signatures for tumor diagnosis and treatment in multiple cancers, including MM. Methods We systematically examined the AS events and clinical information of 83 MM samples from TCGA database. Univariate Cox regression analysis was used to identify AS events associated with overall survival. LASSO analyses followed by multivariate Cox regression analyses were conducted to construct the prognostic signatures and assess the accuracy of these prognostic signatures by receiver operating characteristic (ROC) curve and Kaplan–Meier survival analyses. The ImmuCellAI and ssGSEA algorithms were used to assess the degrees of immune cell infiltration in MM samples. The survival-related splicing regulatory network was established based on the correlation between survival-related AS events and splicing factors (SFs). Results A total of 3976 AS events associated with overall survival were identified by univariate Cox regression analysis, and ES events accounted for the greatest proportion. We constructed prognostic signatures based on survival-related AS events. The prognostic signatures proved to be an efficient predictor with an area under the curve (AUC) greater than 0.9. Additionally, the risk score based on 6 key AS events proved to be an independent prognostic factor, and a nomogram composed of 6 key AS events was established. We found that the risk score was significantly decreased in patients with the epithelioid subtype. In addition, unsupervised clustering clearly showed that the risk score was associated with immune cell infiltration. The abundances of cytotoxic T (Tc) cells, natural killer (NK) cells and T-helper 17 (Th17) cells were higher in the high-risk group, whereas the abundances of induced regulatory T (iTreg) cells were lower in the high-risk group. Finally, we identified 3 SFs (HSPB1, INTS1 and LUC7L2) that were significantly associated with MM patient survival and then constructed a regulatory network between the 3 SFs and survival-related AS to reveal potential regulatory mechanisms in MM. Conclusion Our study provided a prognostic signature based on 6 key events, representing a better effective tumor-specific diagnostic and prognostic marker than the TNM staging system. AS events that are correlated with the immune system may be potential therapeutic targets for MM.

2021 ◽  
Author(s):  
BO SONG ◽  
Lijun Tian ◽  
Fan Zhang ◽  
Zheyu Lin ◽  
Boshen Gong ◽  
...  

Abstract Background: Thyroid cancer (TC) is the most common endocrine malignancy worldwide. The incidence of TC is high and increasing worldwide due to continuous improvements in diagnostic technology. TC is still often overtreated due to a lack of reliable diagnostic biomarkers. Therefore, determining accurate prognostic predictions to stratify TC patients is important.Methods: Raw data were downloaded from the TCGA database, and pairwise comparisons were applied to identify differentially expressed immune-related lncRNA (DEirlncRNA) pairs. Then, we used univariate Cox regression analysis and a modified Lasso algorithm on these pairs to construct a risk assessment model for TC. Next, TC patients were assigned to high- and low-risk groups based on the optimal cutoff score of the model for the 1-year ROC curve. We evaluated the signature in terms of prognostic independence, predictive value, immune cell infiltration, ICI-related molecules and small-molecule inhibitor efficacy. Results: We identified 30 DEirlncRNA pairs through Lasso regression, and 14 pairs served as the novel predictive signature. The high-risk group had a significantly poorer prognosis than the low-risk group. Cox regression analysis revealed that this immune-related signature can predict prognosis independently and reliably for TC. With the CIBERSORT algorithm, we found an association between the signature and immune cell infiltration. Additionally, several immune checkpoint inhibitor (ICI)-related molecules, such as PD-1 and PD-L1, showed a negative correlation with the high-risk group. We further found that some commonly used small-molecule inhibitors, such as sunitinib, were related to this new signature. Conclusions: We constructed a prognostic immune-related lncRNA signature that can predict TC patient survival without considering the technical bias of different platforms, and this signature also sheds light on TC overall prognosis and novel clinical treatments, such as ICB therapy and small molecular inhibitors.


2021 ◽  
Vol 2021 ◽  
pp. 1-12
Author(s):  
Rui Wang ◽  
Wenxuan Bu ◽  
Yang Yang

Multiple myeloma (MM) is the second most commonly diagnosed hematological malignancy. Understanding the basic mechanisms of the metabolism in MM may lead to new therapies that benefit patients. We collected the gene expression profile data of GSE39754 and performed differential analysis. Furthermore, identify the candidate genes that affect the prognosis of the differentially expressed genes (DEGs) related to the metabolism. Enrichment analysis is used to identify the biological effects of candidate genes. Perform coexpression analysis on the verified DEGs. In addition, the candidate genes are used to cluster MM into different subtypes through consistent clustering. Use LASSO regression analysis to identify key genes, and use Cox regression analysis to evaluate the prognostic effects of key genes. Evaluation of immune cell infiltration in MM is by CIBERSORT. We identified 2821 DEGs, of which 348 genes were metabolic-related prognostic genes and were considered candidate genes. Enrichment analysis revealed that the candidate genes are mainly related to the proteasome, purine metabolism, and cysteine and methionine metabolism signaling pathways. According to the consensus clustering method, we identified the two subtypes of group 1 and group 2 that affect the prognosis of MM patients. Using the LASSO model, we have identified 10 key genes. The prognosis of the high-risk group identified by Cox regression analysis is worse than that of the low-risk group. Among them, PKLR has a greater impact on the prognosis of MM, and the prognosis of MM patients is poor when the expression is high. In addition, the level of immune cell infiltration in the high-risk group is higher than that in the low-risk group. In the summary, metabolism-related genes significantly affect the prognosis of MM patients through the metabolic process of MM patients. PKLR may be a prognostic risk factor for MM patients.


2021 ◽  
Author(s):  
Qijun Xie ◽  
Yuwei Zhang ◽  
Wu Huang ◽  
Haoran Wang ◽  
Fang Liu

Abstract Glioblastoma (GBM) is the most common and aggressive primary tumor of the central nervous system with high recurrence and extremely poor prognosis. Multiple recent studies have indicated a pivotal correlation between GBM prognosis and immune-related risk signature. Nevertheless, the potential value of endothelial cells (ECs) immune-related genes (EIRGs) in prognosis, immune infiltration, and their correlation with therapeutic response to immunotherapy and TMZ chemotherapy remain obscure, especially in GBM. Here, we screened out 11 EIRGs after intersecting the identified 59 GBM ECs related prognostic genes and the identified 438 immune-related prognostic genes. A prognostic-related 6-EIRGs signature was established through univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) Cox regression analysis, and patients in the high-risk group were significantly worse overall survival (OS) compared to those in the low-risk group. Additionally, univariate and multivariate Cox regression analysis confirmed that risk score was an independent predictor of OS in patients with GBM. The nomogram which comprised age, gender, IDH mutation status, radiation therapy, and risk score yielded a strong predictive ability of 0.5-, 1-, and 2-year OS for GBM patients. Our results demonstrate that the EIRGs signature, which is associated with immune cell infiltration, may play a regulatory role in the immunobiological process of TIME. Prognostic-related 6-EIRGs signature is a promising classification index for predicting the drug sensitivity to immunotherapy and TMZ chemotherapy, suggesting that EIRGs signature may serve as a biomarker to stratify patients who will benefit from immunotherapy and chemotherapy.


2021 ◽  
Author(s):  
Shaopei Ye ◽  
Wenbin Tang ◽  
Ke Huang

Abstract Background: Autophagy is a biological process to eliminate dysfunctional organelles, aggregates or even long-lived proteins. . Nevertheless, the potential function and prognostic values of autophagy in Wilms Tumor (WT) are complex and remain to be clarifed. Therefore, we proposed to systematically examine the roles of autophagy-associated genes (ARGs) in WT.Methods: Here, we obtained differentially expressed autophagy-related genes (ARGs) between healthy and Wilms tumor from Therapeutically Applicable Research To Generate Effective Treatments(TARGET) and The Cancer Genome Atlas (TCGA) database. The functionalities of the differentially expressed ARGs were analyzed using Gene Ontology. Then univariate COX regression analysis and multivariate COX regression analysis were performed to acquire nine autophagy genes related to WT patients’ survival. According to the risk score, the patients were divided into high-risk and low-risk groups. The Kaplan-Meier curve demonstrated that patients with a high-risk score tend to have a poor prognosis.Results: Eighteen DEARGs were identifed, and nine ARGs were fnally utilized to establish the FAGs based signature in the TCGA cohort. we found that patients in the high-risk group were associated with mutations in TP53. We further conducted CIBERSORT analysis, and found that the infiltration of Macrophage M1 was increased in the high-risk group. Finally, the expression levels of crucial ARGs were verifed by the experiment, which were consistent with our bioinformatics analysis.Conclusions: we emphasized the clinical significance of autophagy in WT, established a prediction system based on autophagy, and identified a promising therapeutic target of autophagy for WT.


2021 ◽  
Author(s):  
Sijia Li ◽  
Hongyang Zhang ◽  
Wei Li

Abstract Background: The purpose of our study is establishing a model based on ferroptosis-related genes predicting the prognosis of patients with head and neck squamous cell carcinoma (HNSCC).Methods: In our study, transcriptome and clinical data of HNSCC patients were from The Cancer Genome Atlas, ferroptosis-related genes and pathways were from Ferroptosis Signatures Database. Differentially expressed genes (DEGs) were screened by comparing tumor and adjacent normal tissues. Functional enrichment analysis of DEGs, protein-protein interaction network and gene mutation examination were applied. Univariate Cox regression analysis and least absolute shrinkage and selection operator (LASSO) regression were used to identified DEGs. The model was constructed by multivariate Cox regression analysis and verified by Kaplan-Meier analysis. The relationship between risk scores and other clinical features was also analyzed. Univariate and multivariate Cox analysis was used to verified the independence of our model. The model was evaluated by receiver operating characteristic analysis and calculation of the area under the curve (AUC). A nomogram model based on risk score, age, gender and TNM stages was constructed.Results: We analyzed data including 500 tumor tissues and 44 adjacent normal tissues and 259 ferroptosis-related genes, then obtained 73 DEGs. Univariate Cox regression analysis screened out 16 genes related to overall survival, and LASSO analysis fingered out 12 of them with prognostic value. A risk score model based on these 12 genes was constructed by multivariate Cox regression analysis. According to the median risk score, patients were divided into high-risk group and low-risk group. The survival rate of high-risk group was significantly lower than that of low-risk group in Kaplan-Meier curve. Risk scores were related to T and grade. Univariate and multivariate Cox analysis showed our model was an independent prognostic factor. The AUC was 0.669. The nomogram showed high accuracy predicting the prognosis of HNSCC patients.Conclusion: Our model based on 12 ferroptosis-related genes performed excellently in predicting the prognosis of HNSCC patients. Ferroptosis-related genes may be promising biomarkers for HNSCC treatment and prognosis.


2022 ◽  
Vol 20 (1) ◽  
Author(s):  
Junyu Huo ◽  
Ge Guan ◽  
Jinzhen Cai ◽  
Liqun Wu

Abstract Background Stromal cells in tumor microenvironment could promote immune escape through a variety of mechanisms, but there are lacking research in the field of gastric cancer (GC). Methods We identified differential expressed immune-related genes (DEIRGs) between the high- and low-stromal cell abundance GC samples in The Cancer Genome Atlas and GSE84437 datasets. A risk score was constructed basing on univariate cox regression analysis, LASSO regression analysis, and multivariate cox regression analysis in the training cohort (n=772). The median value of the risk score was used to classify patients into groups with high and low risk. We conducted external validation of the prognostic signature in four independent cohorts (GSE26253, n=432; GSE62254, n=300; GSE15459, n=191; GSE26901, n=109) from the Gene Expression Omnibus (GEO) database. The immune cell infiltration was quantified by the CIBERSORT method. Results The risk score contained 6 genes (AKT3, APOD, FAM19A5, LTBP3, NOV, and NOX4) showed good performance in predicting 5-year overall survival (OS) rate and 5-year recurrence-free survival (RFS) rate of GC patients. The risk death and recurrence of GC patients growing with the increasing risk score. The patients were clustered into three subtypes according to the infiltration of 22 kinds of immune cells quantified by the CIBERSORT method. The proportion of cluster A with the worst prognosis in the high-risk group was significantly higher than that in the low-risk group; the risk score of cluster C subtype with the best prognosis was significantly lower than that of the other two subtypes. Conclusion This study established and validated a robust prognostic model for gastric cancer by integrated analysis 1804 samples of six centers, and its mechanism was explored in combination with immune cell infiltration characterization.


BMC Cancer ◽  
2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Xuelong Wang ◽  
Bin Zhou ◽  
Yuxin Xia ◽  
Jianxin Zuo ◽  
Yanchao Liu ◽  
...  

Abstract Background DNA methylation alteration is frequently observed in Lung adenocarcinoma (LUAD) and may play important roles in carcinogenesis, diagnosis, and prognosis. Thus, this study aimed to construct a reliable methylation-based nomogram, guiding prognostic classification screening and personalized medicine for LUAD patients. Method The DNA methylation data, gene expression data and corresponding clinical information of lung adenocarcinoma samples were extracted from The Cancer Genome Atlas (TCGA) database. Differentially methylated sites (DMSs) and differentially expressed genes (DEGs) were obtained and then calculated correlation by pearson correlation coefficient. Functional enrichment analysis and Protein-protein interaction network were used to explore the biological roles of aberrant methylation genes. A prognostic risk score model was constructed using univariate Cox and LASSO analysis and was assessed in an independent cohort. A methylation-based nomogram that included the risk score and the clinical risk factors was developed, which was evaluated by concordance index and calibration curves. Result We identified a total of 1362 DMSs corresponding to 471 DEGs with significant negative correlation, including 752 hypermethylation sites and 610 hypomethylation sites. Univariate cox regression analysis showed that 59 DMSs were significantly associated with overall survival. Using LASSO method, we constructed a three-DMSs signature that was independent predictive of prognosis in the training cohort. Patients in high-risk group had a significant shorter overall survival than patients in low-risk group classified by three-DMSs signature (log-rank p = 1.9E-04). Multivariate cox regression analysis proved that the three-DMSs signature was an independent prognostic factor for LUAD in TCGA-LUAD cohort (HR = 2.29, 95%CI: 1.47–3.57, P = 2.36E-04) and GSE56044 cohort (HR = 2.16, 95%CI: 1.19–3.91, P = 0.011). Furthermore, a nomogram, combining the risk score with clinical risk factors, was developed with C-indexes of 0.71 and 0.70 in TCGA-LUAD and GSE56044 respectively. Conclusions The present study established a robust three-DMSs signature for the prediction of overall survival and further developed a nomogram that could be a clinically available guide for personalized treatment of LUAD patients.


2021 ◽  
Vol 12 ◽  
Author(s):  
Tao Han ◽  
Zhifan Zuo ◽  
Meilin Qu ◽  
Yinghui Zhou ◽  
Qing Li ◽  
...  

Background: Although low-grade glioma (LGG) has a good prognosis, it is prone to malignant transformation into high-grade glioma. It has been confirmed that the characteristics of inflammatory factors and immune microenvironment are closely related to the occurrence and development of tumors. It is necessary to clarify the role of inflammatory genes and immune infiltration in LGG.Methods: We downloaded the transcriptome gene expression data and corresponding clinical data of LGG patients from the TCGA and GTEX databases to screen prognosis-related differentially expressed inflammatory genes with the difference analysis and single-factor Cox regression analysis. The prognostic risk model was constructed by LASSO Cox regression analysis, which enables us to compare the overall survival rate of high- and low-risk groups in the model by Kaplan–Meier analysis and subsequently draw the risk curve and survival status diagram. We analyzed the accuracy of the prediction model via ROC curves and performed GSEA enrichment analysis. The ssGSEA algorithm was used to calculate the score of immune cell infiltration and the activity of immune-related pathways. The CellMiner database was used to study drug sensitivity.Results: In this study, 3 genes (CALCRL, MMP14, and SELL) were selected from 9 prognosis-related differential inflammation genes through LASSO Cox regression analysis to construct a prognostic risk model. Further analysis showed that the risk score was negatively correlated with the prognosis, and the ROC curve showed that the accuracy of the model was better. The age, grade, and risk score can be used as independent prognostic factors (p < 0.001). GSEA analysis confirmed that 6 immune-related pathways were enriched in the high-risk group. We found that the degree of infiltration of 12 immune cell subpopulations and the scores of 13 immune functions and pathways in the high-risk group were significantly increased by applying the ssGSEA method (p < 0.05). Finally, we explored the relationship between the genes in the model and the susceptibility of drugs.Conclusion: This study analyzed the correlation between the inflammation-related risk model and the immune microenvironment. It is expected to provide a reference for the screening of LGG prognostic markers and the evaluation of immune response.


2021 ◽  
Vol 2021 ◽  
pp. 1-14
Author(s):  
Qing Ma ◽  
Kai Geng ◽  
Ping Xiao ◽  
Lili Zeng

Background. Non-small-cell lung cancer (NSCLC) is a prevalent malignancy with high mortality and poor prognosis. The radiotherapy is one of the most common treatments of NSCLC, and the radiotherapy sensitivity of patients could affect the individual prognosis of NSCLC. However, the prognostic signatures related to radiotherapy response still remain limited. Here, we explored the radiosensitivity-associated genes and constructed the prognostically predictive model of NSCLC cases. Methods. The NSCLC samples with radiotherapy records were obtained from The Cancer Genome Atlas database, and the mRNA expression profiles of NSCLC patients from the GSE30219 and GSE31210 datasets were obtained from the Gene Expression Omnibus database. The Weighted Gene Coexpression Network Analysis (WGCNA), univariate, least absolute shrinkage and selection operator (LASSO), multivariate Cox regression analysis, and nomogram were conducted to identify and validate the radiotherapy sensitivity-related signature. Results. WGCNA revealed that 365 genes were significantly correlated with radiotherapy response. LASSO Cox regression analysis identified 8 genes, including FOLR3, SLC6A11, ALPP, IGFN1, KCNJ12, RPS4XP22, HIST1H2BH, and BLACAT1. The overall survival (OS) of the low-risk group was better than that of the high-risk group separated by the Risk Score based on these 8 genes for the NSCLC patients. Furthermore, the immune infiltration analysis showed that monocytes and activated memory CD4 T cells had different relative proportions in the low-risk group compared with the high-risk group. The Risk Score was correlated with immune checkpoints, including CTLA4, PDL1, LAG3, and TIGIT. Conclusion. We identified 365 genes potentially correlated with the radiotherapy response of NSCLC patients. The Risk Score model based on the identified 8 genes can predict the prognosis of NSCLC patients.


2020 ◽  
Author(s):  
Xuelong Wang ◽  
Bin Zhou ◽  
Yuxin Xia ◽  
Jianxin Zuo ◽  
Yanchao Liu ◽  
...  

Abstract Background DNA methylation alteration is frequently observed in Lung adenocarcinoma (LUAD) and may play important roles in carcinogenesis, diagnosis, and prognosis. Thus, this study aimed to construct a reliable methylation-based nomogram, guiding prognostic classification screening and personalized medicine for LUAD patients. Method: The DNA methylation data, gene expression data and corresponding clinical information of lung adenocarcinoma samples were extracted from The Cancer Genome Atlas (TCGA) database. Differentially methylated sites (DMSs) and differentially expressed genes (DEGs) were obtained and then calculated expression correlation by pearson correlation coefficient. Functional enrichment analysis and Protein-protein interaction network were used to explore the biological roles of aberrant methylation genes. A prognostic risk score model was constructed using univariate Cox and LASSO analysis and was assessed in an independent cohort. A methylation-based nomogram that included the risk score and the clinical risk factors was developed, which was evaluated by concordance index and calibration curves. Result We identified a total of 1362 DMSs corresponding to 471 DEGs with significant negative correlation, including 752 hypermethylation sites and 610 hypomethylation sites. Univariate cox regression analysis showed that 59 DMSs were significantly associated with overall survival. Using LASSO method, we constructed a three-DMSs signature that was independent predictive of prognosis in the training cohort. Patients in high-risk group had a significant shorter overall survival than patients in low-risk group classified by three-DMSs signature (log-rank p = 1.9E-04). Multivariate cox regression analysis proved that the three-DMSs signature was an independent prognostic factor for LUAD in TCGA-LUAD cohort (HR = 2.29, 95%CI: 1.47–3.57, P = 2.36E-04) and GSE56044 cohort (HR = 2.16, 95%CI: 1.19–3.91, P = 0.011). Furthermore, a nomogram, combining the risk score with clinical risk factors, was developed with C-indexes of 0.71 and 0.70 in TCGA-LUAD and GSE56044 respectively. Conclusions The present study established a robust three-DMSs signature for the prediction of overall survival and further developed a nomogram that could be a clinically available guide for personalized treatment of LUAD patients.


Sign in / Sign up

Export Citation Format

Share Document