scholarly journals A Novel Signature to Predict Thyroid Cancer Prognosis and Immune Landscape Using Immune-Related LncRNA Pairs

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.

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 ◽  
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-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):  
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 ◽  
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.


2021 ◽  
Vol 8 ◽  
Author(s):  
Ji Yin ◽  
Xiaohui Li ◽  
Caifeng Lv ◽  
Xian He ◽  
Xiaoqin Luo ◽  
...  

Background: Long non-coding RNA (lncRNA) plays a significant role in the development, establishment, and progression of head and neck squamous cell carcinoma (HNSCC). This article aims to develop an immune-related lncRNA (irlncRNA) model, regardless of expression levels, for risk assessment and prognosis prediction in HNSCC patients.Methods: We obtained clinical data and corresponding full transcriptome expression of HNSCC patients from TCGA, downloaded GTF files to distinguish lncRNAs from Ensembl, discerned irlncRNAs based on co-expression analysis, distinguished differentially expressed irlncRNAs (DEirlncRNAs), and paired these DEirlncRNAs. Univariate Cox regression analysis, LASSO regression analysis, and stepwise multivariate Cox regression analysis were then performed to screen lncRNA pairs, calculate the risk coefficient, and establish a prognosis model. Finally, the predictive power of this model was validated through the AUC and the ROC curves, and the AIC values of each point on the five-year ROC curve were calculated to select the maximum inflection point, which was applied as a cut-off point to divide patients into low- or high-risk groups. Based on this methodology, we were able to more effectively differentiate between these groups in terms of survival, clinico-pathological characteristics, tumor immune infiltrating status, chemotherapeutics sensitivity, and immunosuppressive molecules.Results: A 13-irlncRNA-pair signature was built, and the ROC analysis demonstrated high sensitivity and specificity of this signature for survival prediction. The Kaplan–Meier analysis indicated that the high-risk group had a significantly shorter survival rate than the low-risk group, and the chi-squared test certified that the signature was highly related to survival status, clinical stage, T stage, and N stage. Additionally, the signature was further proven to be an independent prognostic risk factor via the Cox regression analyses, and immune infiltrating analyses showed that the high-risk group had significant negative relationships with various immune infiltrations. Finally, the chemotherapeutics sensitivity and the expression level of molecular markers were also significantly different between high- and low-risk groups.Conclusion: The signature established by paring irlncRNAs, with regard to specific expression levels, can be utilized for survival prediction and to guide clinical therapy in HNSCC.


2021 ◽  
Vol 8 ◽  
Author(s):  
Haige Zheng ◽  
Huixian Liu ◽  
Yumin Lu ◽  
Hengguo Li

Background: Head and neck squamous cell carcinoma (HNSCC) is a highly heterogeneous tumor with a high incidence and poor prognosis. Therefore, effective predictive models are needed to evaluate patient outcomes and optimize treatment.Methods: Robust Rank Aggregation (RRA) method was used to identify highly robust differentially-expressed genes (DEGs) between HNSCC and normal tissue in 9 GEO and TCGA datasets. Univariate Cox regression analysis and Lasso Cox regression analysis were performed to identify DEGs related to the Overall survival (OS) and to construct a prognostic gene signature (HNSCCSig). External validation was performed using GSE65858 dataset. Moreover, comprehensive bioinformatics analyses were used to identify the association between HNSCCSig and tumor immune environment.Results: A total of 257 reliable DEGs were identified by differentially analysis result of TCGA and GSE65858 datasets. The HNSCCSig including 7 mRNAs (SLURP1, SCARA5, CLDN10, MYH11, CXCL13, HLF, and ITGA3) were developed and validated to identify high-risk group who had a worse OS than low-risk group in TCGA and GSE65858 datasets. Cox regression analysis showed that the HNSCCSig could independently predict OS in both the TCGA and the GSE65858 datasets. Further research demonstrated that the infiltration bundance of CD8 + T cells, B cells, neutrophils, and NK cells were significantly lower in the high-risk group. A nomogram was also constructed by combining the HNSCCSig and clinical characters.Conclusion: We established and validated the HNSCCSig consisting of SLURP1, SCARA5, CLDN10, MYH11, CXCL13, HLF, and ITGA3. A nomogram combining HNSCCSig and some clinical parameters was constructed to identify high-risk HNSCC-patients with poor prognosis.


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.


2020 ◽  
Author(s):  
Kankan Zhao ◽  
Mengchuan Wang ◽  
Houlong Kang ◽  
Aiguo Wu

Abstract Background: Our study aimed to identify immune related long non-coding RNAs (LncRNAs) to serve as potential prognostic indicators and immune therapeutic targets in patients with colon cancer.Methods: The Cancer Genome Atlas (TCGA) and Molecular Signatures Databases (MSigDB) database were used to identify immune related lncRNAs in patients with colon cancer. The least absolute shrinkage and selection operator (LASSO) analysis, and multivariate Cox proportional hazards regression analysis were employed to screen prognostic lncRNAs and construct immune-related multi-lncRNA signature. We used time-dependent receiver operating characteristic curve to assess the performance of this signature in colon cancer by calculating the area under the curve (AUC). Univariate and multivariate Cox regression analysis were performed to verify the independence of the prognostic value of this signature in colon cancer.Results: Five immune related lncRNAs (AC025575.2, AL161729.4, ELFN1-AS1, LINC00513, MIR210HG) were found to be significantly associated with overall survival (OS) of patients with colon cancer. Then, we developed a five immune-related lncRNA signature. According to this signature, patients were ranked into a high risk group (n = 208) and a low risk group (n = 209). Kaplan-Meier curve and log-rank method showed that patients in high risk group had worse OS than patients in low risk group (P = 5.5644e-05). AUC for predicting 3 year survival and 5 year survival was 0.776 and 0.762 respectively, which indicated good performance of this signature. Finally, this five immune-related lncRNA signature was demonstrated to be independently associated with prognosis of patients with colon cancer.Conclusion: We developed a five immune-related lncRNA signature as a prognostic biomarker for patients with colon cancer.


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.


Sign in / Sign up

Export Citation Format

Share Document