scholarly journals Identification and Validation of a Prognostic Immune-Related Alternative Splicing Events Signature for Glioma

2021 ◽  
Vol 11 ◽  
Author(s):  
Minjie Wang ◽  
Zijie Zhou ◽  
Jianglin Zheng ◽  
Wenxuan Xiao ◽  
Jiameng Zhu ◽  
...  

BackgroundGlioma is the most common malignant brain tumor in adults, with its tumor-promoting immune microenvironment always being intricate to handle with. Amounts of evidence has accumulated to suggest that alternative splicing (AS) is related to tumor immune microenvironment. However, comprehensive analysis of immune-related AS events and their clinical significance are still lacking in glioma.MethodsAS events and transcriptome data of 653 glioma patients were downloaded online. ssGSEA was performed on transcriptome data of 653 patients to divided them into low, medium and high immune cell infiltration groups. Immune-related AS events were filtrated based on this grouping. Then lasso Cox regression analysis and multivariate Cox regression analysis were done to achieve an immune-related AS events prognostic signature for glioma. Kaplan-Meier analysis, ROC analyses, univariate Cox regression and multivariate Cox regression were performed to reveal the independent prognostic role of this signature. Meanwhile, a nomogram was constructed to achieved better prognostic value for glioma patients. Besides, functional enrichment analyses and correlation analyses with immune cells infiltration were used to validated the immune-related characteristic of this signature.Results36 immune-related AS events were achieved based on the grouping mentioned above. A nine-immune-related alternative splicing event signature was built for glioma patients. This signature showed an independent prognostic value and a nomogram containing gender, age, Karnofsky performance score, grade, IDH status, MGMT promoter status and risk score derived from the signature was constructed with a higher predictive ability for overall survival. Association with the infiltration of immune cell subtypes was validated and functional enrichment analysis found that the signature was mainly enriched in immune-related and pro-tumor functions.ConclusionOur research presented all immune-related AS events in glioma, identified an immune-related prognostic AS events risk model and a nomogram was constructed to predict the prognosis individually and more precisely. Tight connection was verified between this signature and clinical characteristics. Also, immune cells infiltration and immune checkpoints expression level were proved to link to risk scores, which enhanced the understanding of relationship between AS events and glioma immune microenvironment, firstly revealing the potential role of AS in immunotherapy of glioma.

2021 ◽  
Vol 12 ◽  
Author(s):  
Tao Lin ◽  
Hao Cheng ◽  
Da Liu ◽  
Lei Wen ◽  
Junlin Kang ◽  
...  

Since autophagy and the immune microenvironment are deeply involved in the tumor development and progression of Lower-grade gliomas (LGG), our study aimed to construct an autophagy-related risk model for prognosis prediction and investigate the relationship between the immune microenvironment and risk signature in LGG. Therefore, we identified six autophagy-related genes (BAG1, PTK6, EEF2, PEA15, ITGA6, and MAP1LC3C) to build in the training cohort (n = 305 patients) and verify the prognostic model in the validation cohort (n = 128) and the whole cohort (n = 433), based on the data from The Cancer Genome Atlas (TCGA). The six-gene risk signature could divide LGG patients into high- and low-risk groups with distinct overall survival in multiple cohorts (all p < 0.001). The prognostic effect was assessed by area under the time-dependent ROC (t-ROC) analysis in the training, validation, and whole cohorts, in which the AUC value at the survival time of 5 years was 0.837, 0.755, and 0.803, respectively. Cox regression analysis demonstrated that the risk model was an independent risk predictor of OS (HR > 1, p < 0.05). A nomogram including the traditional clinical parameters and risk signature was constructed, and t-ROC, C-index, and calibration curves confirmed its robust predictive capacity. KM analysis revealed a significant difference in the subgroup analyses’ survival. Functional enrichment analysis revealed that these autophagy-related signatures were mainly involved in the phagosome and immune-related pathways. Besides, we also found significant differences in immune cell infiltration and immunotherapy targets between risk groups. In conclusion, we built a powerful predictive signature and explored immune components (including immune cells and emerging immunotherapy targets) in LGG.


2021 ◽  
Vol 11 ◽  
Author(s):  
Liping Zhu ◽  
Zhiqiang Wang ◽  
Yilan Sun ◽  
Georgios Giamas ◽  
Justin Stebbing ◽  
...  

BackgroundAlternative splicing (AS) is a gene regulatory mechanism that drives protein diversity. Dysregulation of AS is thought to play an essential role in cancer initiation and development. This study aimed to construct a prognostic signature based on AS and explore the role in the tumor immune microenvironment (TIME) in lung adenocarcinoma.MethodsWe analyzed transcriptome profiling and clinical lung adenocarcinoma data from The Cancer Genome Atlas (TCGA) database and lists of AS-related and immune-related signatures from the SpliceSeq. Prognosis-related AS events were analyzed by univariate Cox regression analysis. Gene set enrichment analyses (GSEA) were performed for functional annotation. Prognostic signatures were identified and validated using univariate and multivariate Cox regression, LASSO regression, Kaplan–Meier survival analyses, and proportional hazards model. The context of TIME in lung adenocarcinoma was also analyzed. Gene and protein expression data of Cyclin-Dependent Kinase Inhibitor 2A (CDKN2A) were obtained from ONCOMINE and Human Protein Atlas. Splicing factor (SF) regulatory networks were visualized.ResultsA total of 19,054 survival-related AS events in lung adenocarcinoma were screened in 1,323 genes. Exon skip (ES) and mutually exclusive exons (ME) exhibited the most and fewest AS events, respectively. Based on AS subtypes, eight AS prognostic signatures were constructed. Patients with high-risk scores were associated with poor overall survival. A nomogram with good validity in prognostic prediction was generated. AUCs of risk scores at 1, 2, and 3 years were 0.775, 0.736, and 0.759, respectively. Furthermore, the prognostic signatures were significantly correlated with TIME diversity and immune checkpoint inhibitor (ICI)-related genes. Low-risk patients had a higher StromalScore, ImmuneScore, and ESTIMATEScore. AS-based risk score signature was positively associated with CD8+ T cells. CDKN2A was also found to be a prognostic factor in lung adenocarcinoma. Finally, potential functions of SFs were determined by regulatory networks.ConclusionTaken together, our findings show a clear association between AS and immune cell infiltration events and patient outcome, which could provide a basis for the identification of novel markers and therapeutic targets for lung adenocarcinoma. SF networks provide information of regulatory mechanisms.


2021 ◽  
Vol 14 (8) ◽  
pp. 1151-1159
Author(s):  
Chen-Lu Liao ◽  
◽  
Xing-Yu Sun ◽  
Qi Zhou ◽  
Min Tian ◽  
...  

AIM: To investigate the role of tumor microenvironment (TME)-related long non-coding RNA (lncRNA) in uveal melanoma (UM), probable prognostic signature and potential small molecule drugs using bioinformatics analysis. METHODS: UM expression profile data were downloaded from the Cancer Genome Atlas (TCGA) and bioinformatics methods were used to find prognostic lncRNAs related to UM immune cell infiltration. The gene expression profile data of 80 TCGA specimens were analyzed using the single sample Gene Set Enrichment Analysis (ssGSEA) method, and the immune cell infiltration of a single specimen was evaluated. Finally, the specimens were divided into high and low infiltration groups. The differential expression between the two groups was analyzed using the R package ‘edgeR’. Univariate, multivariate and Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analyses were performed to explore the prognostic value of TME-related lncRNAs. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses were also performed. The Connectivity Map (CMap) data set was used to screen molecular drugs that may treat UM. RESULTS: A total of 2393 differentially expressed genes were identified and met the criteria for the low and high immune cell infiltration groups. Univariate Cox analysis of lncRNA genes with differential expression identified 186 genes associated with prognosis. Eight prognostic markers of TME-included lncRNA genes were established as potentially independent prognostic elements. Among 269 differentially expressed lncRNAs, 69 were up-regulated and 200 were down-regulated. Univariate Cox regression analysis of the risk indicators and clinical characteristics of the 8 lncRNA gene constructs showed that age, TNM stage, tumor base diameter, and low and high risk indices had significant prognostic value. We screened the potential small-molecule drugs for UM, including W-13, AH-6809 and Imatinib. CONCLUSION: The prognostic markers identified in this study are reliable biomarkers of UM. This study expands our current understanding of the role of TME-related lncRNAs in UM genesis, which may lay the foundations for future treatment of this disease.


2020 ◽  
Vol 2020 ◽  
pp. 1-43
Author(s):  
Beilei Wu ◽  
Lijun Tao ◽  
Daqing Yang ◽  
Wei Li ◽  
Hongbo Xu ◽  
...  

Objective. Stromal cells and immune cells have important clinical significance in the microenvironment of colorectal cancer (CRC). This study is aimed at developing a CRC gene signature on the basis of stromal and immune scores. Methods. A cohort of CRC patients (n=433) were adopted from The Cancer Genome Atlas (TCGA) database. Stromal/immune scores were calculated by the ESTIMATE algorithm. Correlation between prognosis/clinical characteristics and stromal/immune scores was assessed. Differentially expressed stromal and immune genes were identified. Their potential functions were annotated by functional enrichment analysis. Cox regression analysis was used to develop an eight-gene risk score model. Its predictive efficacies for 3 years, 5 years, overall survival (OS), and progression-free survival interval (PFI) were evaluated using time-dependent receiver operating characteristic (ROC) curves. The correlation between the risk score and the infiltering levels of six immune cells was analyzed using TIMER. The risk score was validated using an independent dataset. Results. Immune score was in a significant association with prognosis and clinical characteristics of CRC. 736 upregulated and two downregulated stromal and immune genes were identified, which were mainly enriched into immune-related biological processes and pathways. An-eight gene prognostic risk score model was conducted, consisting of CCL22, CD36, CPA3, CPT1C, KCNE4, NFATC1, RASGRP2, and SLC2A3. High risk score indicated a poor prognosis of patients. The area under the ROC curves (AUC) s of the model for 3 years, 5 years, OS, and PFI were 0.71, 0.70, 0.73, and 0.66, respectively. Thus, the model possessed well performance for prediction of patients’ prognosis, which was confirmed by an external dataset. Moreover, the risk score was significantly correlated with immune cell infiltration. Conclusion. Our study conducted an immune-related prognostic risk score model, which could provide novel targets for immunotherapy of CRC.


2021 ◽  
Author(s):  
Feng Liu ◽  
Zewei Tu ◽  
Junzhe Liu ◽  
Xiaoyan Long ◽  
Bing Xiao ◽  
...  

Abstract Background: A role of DNAJC10 has been reported in several cancers, but its function in glioma is not clear. The purpose of this study was to investigate the prognostic role and the underlying functions of DNAJC10 in glioma.Methods: Reverse transcription and quantitative polymerase chain reaction and western blotting were performed to quantify the relative DNAJC10 mRNA and protein expressions of clinical samples. Wilcoxon rank sum tests were used to compare DNAJC10 expression between or among glioma subgroups with different clinicopathological features. The overall survival (OS) rates of glioma patients with different DNAJC10 expression were compared with the Kaplan-Meier method (two-sided log-rank test). The prognosis-predictive accuracy of the DNAJC10 was evaluated by time-dependent receiver operating characteristic (ROC) curves. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes annotations were conducted using the “clusterProfiler” package. Single-sample gene set enrichment analysis was used to estimate immune cell infiltrations and immune-related function levels. The independent prognostic role of DNAJC10 was determined by univariate and multivariate Cox regression analyses. A DNAJC10-based nomogram model was established using multivariate Cox regression in the R package “rms.” Results: Higher DNAJC10 expression was observed in gliomas. It was upregulated in tumors with higher World Health Organization grade, isocitrate dehydrogenase wild-type status, 1p/19q non-co-deletion, and methylguanine-DNA methyltransferase unmethylated gliomas. Patients with gliomas with higher DNAJC10 expression had poorer prognoses than those with low-DNAJC10 gliomas. The predictive accuracy of 1/3/5-year OS of DNAJC10 was stable and robust using a time-dependent ROC model. Functional enrichment analysis recognized that T cell activation and T cell receptor signaling were enriched in higher DNAJC10 gliomas. Immune cell and stromal cell infiltrations, tumor mutation burden, copy number alteration burden, and immune checkpoint genes were also positively correlated with glioma DNAJC10 expression. A DNAJ10-based nomogram model was established and showed strong prognosis-predictive ability.Conclusion: Higher DNAJC10 expression correlates with poor prognosis of patients with glioma and is a potential and useful prognostic biomarker.


2021 ◽  
Vol 8 ◽  
Author(s):  
Mingqin Ge ◽  
Jie Niu ◽  
Ping Hu ◽  
Aihua Tong ◽  
Yan Dai ◽  
...  

Objective: This study aimed to construct a prognostic ferroptosis-related signature for thyroid cancer and probe into the association with tumor immune microenvironment.Methods: Based on the expression profiles of ferroptosis-related genes, a LASSO cox regression model was established for thyroid cancer. Kaplan-Meier survival analysis was presented between high and low risk groups. The predictive performance was assessed by ROC. The predictive independency was validated via multivariate cox regression analysis and stratified analysis. A nomogram was established and verified by calibration curves. The enriched signaling pathways were predicted via GSEA. The association between the signature and immune cell infiltration was analyzed by CIBERSORT. The ferroptosis-related genes were validated in thyroid cancer tissues by immunohistochemistry and RT-qPCR.Results: A ferroptosis-related eight gene model was established for predicting the prognosis of thyroid cancer. Patients with high risk score indicated a poorer prognosis than those with low risk score (p = 1.186e-03). The AUCs for 1-, 2-, and 3-year survival were 0.887, 0.890, and 0.840, respectively. Following adjusting other prognostic factors, the model could independently predict the prognosis (p = 0.015, HR: 1.870, 95%CI: 1.132–3.090). A nomogram combining the signature and age was constructed. The nomogram-predicted probability of 1-, 3-, and 5-year survival approached the actual survival time. Several ferroptosis-related pathways were enriched in the high-risk group. The signature was distinctly associated with the immune cell infiltration. After validation, the eight genes were abnormally expressed between thyroid cancer and control tissues.Conclusion: Our findings established a prognostic ferroptosis-related signature that was associated with the immune microenvironment for thyroid cancer.


2020 ◽  
Vol 40 (5) ◽  
Author(s):  
Xiaodong Chen ◽  
Fen Tian ◽  
Peng Lun ◽  
Yugong Feng

Abstract Tumor-infiltrating immune cells play a decisive part in prognosis and survival. Until now, previous researches have not made clear about the diversity of cell types involved in the immune response. The objective of this work was to confirm the composition of tumor-infiltrating immune cells and their correlation with prognosis in meningiomas based on a metagene approach (known as CIBERSORT) and online databases. A total of 22 tumor-infiltrating immune cells were detected to determine the relationship between the immune infiltration pattern and survival. The proportion of M2 macrophages was more abundant in 68 samples, reaching more than 36%. Univariate Cox regression analysis displayed that the proportion of dendritic cells was obviously related to prognosis. Hierarchical clustering analysis identified two clusters by the method of within sum of squares errors, which exhibited different infiltrating immune cell composition and survival. To summarize, our results indicated that proportions of tumor-infiltrating immune cells as well as cluster patterns were associated with the prognosis, which offered clinical significance for research of meningiomas.


2021 ◽  
Vol 8 ◽  
Author(s):  
Zengyu Feng ◽  
Kexian Li ◽  
Yulian Wu ◽  
Chenghong Peng

Background: Accumulating evidence shows that the elevated expression of DCBLD2 (discoidin, CUB and LCCL domain-containing protein 2) is associated with unfavorable prognosis of various cancers. However, the correlation of DCBLD2 expression value with the diagnosis and prognosis of pancreatic ductal adenocarcinoma (PDAC) has not yet been elucidated. Methods: Univariate Cox regression analysis was used to screen robust survival-related genes. Expression pattern of selected genes was investigated in PDAC tissues and normal tissues from multiple cohorts. Kaplan–Meier (K–M) survival curves, ROC curves and calibration curves were employed to assess prognostic performance. The relationship between DCBLD2 expression and immune cell infiltrates was conducted by CIBERSORT software. Biological processes and KEGG pathway enrichment analyses were adopted to clarify the potential function of DCBLD2 in PDAC. Results: Univariate analysis, K–M survival curves and calibration curves indicated that DCBLD2 was a robust prognostic factor for PDAC with cross-cohort compatibility. Upregulation of DCBLD2 was observed in dissected PDAC tissues as well as extracellular vesicles from both plasma and serum samples of PDAC patients. Both DCBLD2 expression in tissue and extracellular vesicles had significant diagnostic value. Besides, DCBLD2 expression was correlated with infiltrating level of CD8+ T cells and macrophage M2 cells. Functional enrichment revealed that DCBLD2 might be involved in cell motility, angiogenesis, and cancer-associated pathways. Conclusion: Our study systematically analyzed the potential diagnostic, prognostic and therapeutic value of DCBLD2 in PDAC. All the findings indicated that DCBLD2 might play a considerably oncogenic role in PDAC with diagnostic, prognostic and therapeutic potential. These preliminary results of bioinformatics analyses need to be further validated in more prospective studies.


2021 ◽  
Vol 11 ◽  
Author(s):  
Lu Lu ◽  
Le-Ping Liu ◽  
Qiang-Qiang Zhao ◽  
Rong Gui ◽  
Qin-Yu Zhao

Lung adenocarcinoma (LUAD) is a highly heterogeneous malignancy, which makes prognosis prediction of LUAD very challenging. Ferroptosis is an iron-dependent cell death mechanism that is important in the survival of tumor cells. Long non-coding RNAs (lncRNAs) are considered to be key regulators of LUAD development and are involved in ferroptosis of tumor cells, and ferroptosis-related lncRNAs have gradually emerged as new targets for LUAD treatment and prognosis. It is essential to determine the prognostic value of ferroptosis-related lncRNAs in LUAD. In this study, we obtained RNA sequencing (RNA-seq) data and corresponding clinical information of LUAD patients from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database and ferroptosis-related lncRNAs by co-expression analysis. The best predictors associated with LUAD prognosis, including C5orf64, LINC01800, LINC00968, LINC01352, PGM5-AS1, LINC02097, DEPDC1-AS1, WWC2-AS2, SATB2-AS1, LINC00628, LINC01537, LMO7DN, were identified by Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis, and the LUAD risk prediction model was successfully constructed. Kaplan-Meier analysis, receiver operating characteristic (ROC) time curve analysis and univariate and multivariate Cox regression analysis and further demonstrated that the model has excellent robustness and predictive ability. Further, based on the risk prediction model, functional enrichment analysis revealed that 12 prognostic indicators involved a variety of cellular functions and signaling pathways, and the immune status was different in the high-risk and low-risk groups. In conclusion, a risk model of 12 ferroptosis related lncRNAs has important prognostic value for LUAD and may be ferroptosis-related therapeutic targets in the clinic.


2021 ◽  
Vol 8 ◽  
Author(s):  
Yong Chen ◽  
Fada Xia ◽  
Bo Jiang ◽  
Wenlong Wang ◽  
Xinying Li

Background: Epigenetic regulation, including DNA methylation, plays a major role in shaping the identity and function of immune cells. Innate and adaptive immune cells recruited into tumor tissues contribute to the formation of the tumor immune microenvironment (TIME), which is closely involved in tumor progression in breast cancer (BC). However, the specific methylation signatures of immune cells have not been thoroughly investigated yet. Additionally, it remains unknown whether immune cells-specific methylation signatures can identify subgroups and stratify the prognosis of BC patients.Methods: DNA methylation profiles of six immune cell types from eight datasets downloaded from the Gene Expression Omnibus were collected to identify immune cell-specific hypermethylation signatures (IC-SHMSs). Univariate and multivariate cox regression analyses were performed using BC data obtained from The Cancer Genome Atlas to identify the prognostic value of these IC-SHMSs. An unsupervised clustering analysis of the IC-SHMSs with prognostic value was performed to categorize BC patients into subgroups. Multiple Cox proportional hazard models were constructed to explore the role of IC-SHMSs and their relationship to clinical characteristics in the risk stratification of BC patients. Integrated discrimination improvement (IDI) was performed to determine whether the improvement of IC-SHMSs on clinical characteristics in risk stratification was statistically significant.Results: A total of 655 IC-SHMSs of six immune cell types were identified. Thirty of them had prognostic value, and 10 showed independent prognostic value. Four subgroups of BC patients, which showed significant heterogeneity in terms of survival prognosis and immune landscape, were identified. The model incorporating nine IC-SHMSs showed similar survival prediction accuracy as the clinical model incorporating age and TNM stage [3-year area under the curve (AUC): 0.793 vs. 0.785; 5-year AUC: 0.735 vs. 0.761]. Adding the IC-SHMSs to the clinical model significantly improved its prediction accuracy in risk stratification (3-year AUC: 0.897; 5-year AUC: 0.856). The results of IDI validated the statistical significance of the improvement (p < 0.05).Conclusions: Our study suggests that IC-SHMSs may serve as signatures of classification and risk stratification in BC. Our findings provide new insights into epigenetic signatures, which may help improve subgroup identification, risk stratification, and treatment management.


Sign in / Sign up

Export Citation Format

Share Document