A Five-lncRNAs Signature-Derived Risk Score Based on TCGA and CGGA for Glioblastoma: Potential Prospects for Treatment Evaluation and Prognostic Prediction

Accumulating studies have confirmed the crucial role of long non-coding RNAs (ncRNAs) as favorable biomarkers for cancer diagnosis, therapy, and prognosis prediction. In our recent study, we established a robust model which is based on multi-gene signature to predict the therapeutic efficacy and prognosis in glioblastoma (GBM), based on Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) databases. lncRNA-seq data of GBM from TCGA and CGGA datasets were used to identify differentially expressed genes (DEGs) compared to normal brain tissues. The DEGs were then used for survival analysis by univariate and multivariate COX regression. Then we established a risk score model, depending on the gene signature of multiple survival-associated DEGs. Subsequently, Kaplan-Meier analysis was used for estimating the prognostic and predictive role of the model. Gene set enrichment analysis (GSEA) was applied to investigate the potential pathways associated to high-risk score by the R package “cluster profile” and Wiki-pathway. And five survival associated lncRNAs of GBM were identified: LNC01545, WDR11-AS1, NDUFA6-DT, FRY-AS1, TBX5-AS1. Then the risk score model was established and shows a desirable function for predicting overall survival (OS) in the GBM patients, which means the high-risk score significantly correlated with lower OS both in TCGA and CGGA cohort. GSEA showed that the high-risk score was enriched with PI3K-Akt, VEGFA-VEGFR2, TGF-beta, Notch, T-Cell pathways. Collectively, the five-lncRNAs signature-derived risk score presented satisfactory efficacies in predicting the therapeutic efficacy and prognosis in GBM and will be significant for guiding therapeutic strategies and research direction for GBM.

Download Full-text

Development and Verification of Glutamatergic Synapse-Associated Prognosis Signature for Lower-Grade Gliomas

Frontiers in Molecular Neuroscience ◽

10.3389/fnmol.2021.720899 ◽

2021 ◽

Vol 14 ◽

Author(s):

Liguo Ye ◽

Yang Xu ◽

Ping Hu ◽

Long Wang ◽

Ji’an Yang ◽

...

Keyword(s):

Risk Score ◽

Cox Regression ◽

Gene Set Enrichment Analysis ◽

Normal Brain ◽

Lower Grade ◽

Glutamatergic Synapse ◽

Hub Genes ◽

Gene Set ◽

Genome Atlas

Background: Lower-grade glioma (LGG) is the most common histology identified in gliomas, a heterogeneous tumor that may develop into high-grade malignant glioma that seriously shortens patient survival time. Recent studies reported that glutamatergic synapses might play an essential role in the progress of gliomas. However, the role of glutamatergic synapse-related biomarkers in LGG has not been systemically researched yet.Methods: The mRNA expression data of glioma and normal brain tissue were obtained from The Cancer Genome Atlas database and Genotype-Tissue Expression, respectively, and the Chinese Glioma Genome Atlas database was used as a validation set. Difference analysis was performed to evaluate the expression pattern of glutamatergic synapse-related genes (GSRGs) in LGG. The least absolute shrinkage and selection operator (LASSO) Cox regression was applied to construct the glutamatergic synapse-related risk signature (GSRS), and the risk score of each LGG sample was calculated based on the coefficients and expression value of selected GSRGs. Univariate and multivariate Cox regression analyses were used to investigate the prognostic value of risk score. Immunity profile and single-sample gene set enrichment analysis (ssGSEA) were performed to explore the association between risk score and the characters of tumor microenvironment in LGG. Gene set variation analysis (GSVA) was performed to investigate the potential pathways related to GSRS. The HPA database and real-time PCR were used to identify the expression of hub genes identified in GSRS.Results: A total of 22 genes of 39 GSRGs were found differentially expressed among normal and LGG samples. Through the LASSO algorithm, 14-genes GSRS constructed were associated with the prognosis and clinicopathological features of patients with LGG. Furthermore, the risk score level was significantly positively correlated with the infiltrating level of immunosuppressive cells, including M2 macrophages and regulatory T cells. GSVA identified a series of cancer-related pathways related to GSRS, such as P13K-AKT and P53 pathways. Moreover, ATAD1, NLGN2, OXTR, and TNR, hub genes identified in GSRS, were considered as potential prognostic biomarkers in LGG.Conclusion: A 14-genes GSRS was constructed and verified in this study. We provided a novel insight into the role of GSRS in LGG through a series of bioinformatics methods.

Download Full-text

Clinicopathological signature and prognostic value of RNA:m5C methyltransferases in gliomas

10.21203/rs.2.22601/v1 ◽

2020 ◽

Author(s):

Peng Wang ◽

Kai Huang ◽

Miaojing Wu ◽

Qing Hu ◽

Chuming Tao ◽

...

Keyword(s):

Overall Survival ◽

Risk Score ◽

Prognostic Value ◽

Cox Regression ◽

Roc Curves ◽

Clinicopathological Features ◽

Gene Set Enrichment Analysis ◽

Consensus Clustering ◽

Score Model ◽

Genome Atlas

Abstract Background: Glioma is the most common primary intracranial tumor, accounting for the vast majority of intracranial malignant tumors. Aberrant expression of RNA:5-methylcytosine(m5C) methyltransferases has recently been the focus of research relating to the occurrence and progression of tumors. However, the prognostic value of RNA:m5C methyltransferases in glioma remains unclear. This study investigated RNA: m5C methyltransferase expression and defined its clinicopathological signature and prognostic value in gliomas. Methods: We systematically studied the RNA-sequence data of RNA:m5C methyltransferases underlying gliomas in the Chinese Glioma Genome Atlas (CGGA) and The Cancer Genome Atlas (TCGA) datasets and identified different subtypes using Consensus clustering analysis. Gene Ontology (GO) and Gene Set Enrichment analysis (GSEA) was used to annotate the function of these genes. Univariate Cox regression and the least absolute shrinkage and selection operator (LASSO) Cox regression algorithm analyses were performed to construct the risk score model. Kaplan-Meier method and Receiver operating characteristic (ROC) curves were used to assess the overall survival of glioma patients. Additionally, Cox proportional regression model analysis was developed to address the connections between the risk scores and clinical factors. Results: Consensus clustering of RNA:m5C methyltransferases identified three clusters of gliomas with different prognostic and clinicopathological features. Meanwhile, Functional annotations demonstrated that RNA:m5C methyltransferases were significantly associated with the malignant progression of gliomas. Thereafter, five RNA:m5C methyltransferase genes were screened to construct a risk score model which can be used to predict not only overall survival but also clinicopathological features in gliomas. ROC curves revealed the significant prognostic ability of this signature. In addition, Multivariate Cox regression analyses indicated that the risk score was an independent prognostic factor for glioma outcome. Conclusion: We demonstrated the role of RNA:m5C methyltransferases in the initiation and progression of glioma. We have expanded on the understanding of the molecular mechanism involved, and provided a unique approach to predictive biomarkers and targeted therapy.

Download Full-text

A Novel 10-Gene Signature Predicts Poor Prognosis in Low Grade Glioma

Technology in Cancer Research & Treatment ◽

10.1177/1533033821992084 ◽

2021 ◽

Vol 20 ◽

pp. 153303382199208

Author(s):

Wentao Liu ◽

Jiaxuan Zou ◽

Rijun Ren ◽

Jingping Liu ◽

Gentang Zhang ◽

...

Keyword(s):

Risk Score ◽

Poor Prognosis ◽

Cox Regression ◽

Gene Signature ◽

Low Grade Glioma ◽

The Cancer Genome Atlas ◽

Low Grade ◽

Normal Brain ◽

Cox Regression Analysis ◽

Genome Atlas

Aim: Low grade glioma (LGG) is a lethal brain cancer with relatively poor prognosis in young adults. Thus, this study was performed to develop novel molecular biomarkers to effectively predict the prognosis of LGG patients and finally guide treatment decisions. Methods: survival-related genes were determined by Kaplan-Meier survival analysis and multivariate Cox regression analysis using the expression and clinical data of 506 LGG patients from The Cancer Genome Atlas (TCGA) database and independently validated in a Chinese Glioma Genome Atlas (CGGA) dataset. A prognostic risk score was established based on a linear combination of 10 gene expression levels using the regression coefficients of the multivariate Cox regression models. GSEA was performed to analyze the altered signaling pathways between the high and low risk groups stratified by median risk score. Results: We identified a total of 1489 genes significantly correlated with patients’ prognosis in LGG. The top 5 protective genes were DISP2, CKMT1B, AQP7, GPR162 and CHGB, the top 5 risk genes were SP1, EYA3, ZSCAN20, ITPRIPL1 and ZNF217 in LGG. The risk score was predictive of poor overall survival and relapse-free survival in LGG patients. Pathways of small cell lung cancer, pathways in cancer, chronic myeloid leukemia, colorectal cancer were the top 4 most enriched pathways in the high risk group. SP1, EYA3, ZSCAN20, ITPRIPL1, ZNF217 and GPR162 were significantly up-regulated, while DISP2, CKMT1B, AQP7 were down-regulated in 523 LGG tissues as compared to 1141 normal brain controls. Conclusions: The 10-gene signature may become novel prognostic and diagnostic biomarkers to considerably improve the prognostic prediction in LGG.

Download Full-text

A glycolysis-based three-gene signature predicts survival in patients with lung squamous cell carcinoma

BMC Cancer ◽

10.1186/s12885-021-08360-z ◽

2021 ◽

Vol 21 (1) ◽

Author(s):

Guichuan Huang ◽

Jing Zhang ◽

Ling Gong ◽

Yi Huang ◽

Daishun Liu

Keyword(s):

Squamous Cell Carcinoma ◽

Cell Carcinoma ◽

Risk Score ◽

Squamous Cell ◽

Cox Regression ◽

Lung Squamous Cell Carcinoma ◽

Gene Signature ◽

Gene Set Enrichment Analysis ◽

Related Gene ◽

Cox Regression Analysis

Abstract Background Lung cancer is one of the most lethal and most prevalent malignant tumors worldwide, and lung squamous cell carcinoma (LUSC) is one of the major histological subtypes. Although numerous biomarkers have been found to be associated with prognosis in LUSC, the prediction effect of a single gene biomarker is insufficient, especially for glycolysis-related genes. Therefore, we aimed to develop a novel glycolysis-related gene signature to predict survival in patients with LUSC. Methods The mRNA expression files and LUSC clinical information were obtained from The Cancer Genome Atlas (TCGA) dataset. Results Based on Gene Set Enrichment Analysis (GSEA), we found 5 glycolysis-related gene sets that were significantly enriched in LUSC tissues. Univariate and multivariate Cox proportional regression models were performed to choose prognostic-related gene signatures. Based on a Cox proportional regression model, a risk score for a three-gene signature (HKDC1, ALDH7A1, and MDH1) was established to divide patients into high-risk and low-risk subgroups. Multivariate Cox regression analysis indicated that the risk score for this three-gene signature can be used as an independent prognostic indicator in LUSC. Additionally, based on the cBioPortal database, the rate of genomic alterations in the HKDC1, ALDH7A1, and MDH1 genes were 1.9, 1.1, and 5% in LUSC patients, respectively. Conclusion A glycolysis-based three-gene signature could serve as a novel biomarker in predicting the prognosis of patients with LUSC and it also provides additional gene targets that can be used to cure LUSC patients.

Download Full-text

A Hypoxia-related LncRNA Signature Predicts Survival of Primary Lower-grade glioma

10.21203/rs.3.rs-335140/v1 ◽

2021 ◽

Author(s):

Yanjia Hu ◽

Jing Zhang ◽

Jing Chen

Keyword(s):

High Risk ◽

Risk Score ◽

Cox Regression ◽

Risk Group ◽

Enrichment Analysis ◽

High Risk Group ◽

Gene Set Enrichment Analysis ◽

Lower Grade ◽

Prognostic Signature ◽

Gene Set Enrichment

Abstract Background Hypoxia-related long non-coding RNAs (lncRNAs) have been proven to play a role in multiple cancers and can serve as prognostic markers. Lower-grade gliomas (LGGs) are characterized by large heterogeneity. Methods This study aimed to construct a hypoxia-related lncRNA signature for predicting the prognosis of LGG patients. Transcriptome and clinical data of LGG patients were obtained from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA). LGG cohort in TCGA was chosen as training set and LGG cohorts in CGGA served as validation sets. A prognostic signature consisting of fourteen hypoxia-related lncRNAs was constructed using univariate and LASSO Cox regression. A risk score formula involving the fourteen lncRNAs was developed to calculate the risk score and patients were classified into high- and low-risk groups based on cutoff. Kaplan-Meier survival analysis was used to compare the survival between two groups. Cox regression analysis was used to determine whether risk score was an independent prognostic factor. A nomogram was then constructed based on independent prognostic factors and assessed by C-index and calibration plot. Gene set enrichment analysis and immune cell infiltration analysis were performed to uncover further mechanisms of this lncRNA signature. Results LGG patients with high risk had poorer prognosis than those with low risk in both training and validation sets. Recipient operating characteristic curves showed good performance of the prognostic signature. Univariate and multivariate Cox regression confirmed that the established lncRNA signature was an independent prognostic factor. C-index and calibration plots showed good predictive performance of nomogram. Gene set enrichment analysis showed that genes in the high-risk group were enriched in apoptosis, cell adhesion, pathways in cancer, hypoxia etc. Immune cells were higher in high-risk group. Conclusion The present study showed the value of the 14-lncRNA signature in predicting survival of LGGs and these 14 lncRNAs could be further investigated to reveal more mechanisms involved in gliomas.

Download Full-text

Development and Verification of the Hypoxia- and Immune-Associated Prognostic Signature for Pancreatic Ductal Adenocarcinoma

Frontiers in Immunology ◽

10.3389/fimmu.2021.728062 ◽

2021 ◽

Vol 12 ◽

Author(s):

Dongjie Chen ◽

Hui Huang ◽

Longjun Zang ◽

Wenzhe Gao ◽

Hongwei Zhu ◽

...

Keyword(s):

High Risk ◽

Pancreatic Ductal Adenocarcinoma ◽

Risk Score ◽

Cox Regression ◽

Pearson Correlation ◽

Risk Groups ◽

Low Risk ◽

Ductal Adenocarcinoma ◽

Prognostic Signature ◽

Score Model

We aim to construct a hypoxia- and immune-associated risk score model to predict the prognosis of patients with pancreatic ductal adenocarcinoma (PDAC). By unsupervised consensus clustering algorithms, we generate two different hypoxia clusters. Then, we screened out 682 hypoxia-associated and 528 immune-associated PDAC differentially expressed genes (DEGs) of PDAC using Pearson correlation analysis based on the Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression project (GTEx) dataset. Seven hypoxia and immune-associated signature genes (S100A16, PPP3CA, SEMA3C, PLAU, IL18, GDF11, and NR0B1) were identified to construct a risk score model using the Univariate Cox regression and the Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression, which stratified patients into high- and low-risk groups and were further validated in the GEO and ICGC cohort. Patients in the low-risk group showed superior overall survival (OS) to their high-risk counterparts (p < 0.05). Moreover, it was suggested by multivariate Cox regression that our constructed hypoxia-associated and immune-associated prognosis signature might be used as the independent factor for prognosis prediction (p < 0.001). By CIBERSORT and ESTIMATE algorithms, we discovered that patients in high-risk groups had lower immune score, stromal score, and immune checkpoint expression such as PD-L1, and different immunocyte infiltration states compared with those low-risk patients. The mutation spectrum also differs between high- and low-risk groups. To sum up, our hypoxia- and immune-associated prognostic signature can be used as an approach to stratify the risk of PDAC.

Download Full-text

A Novel Four-Gene Prognostic Signature for Prediction of Survival in Patients with Soft Tissue Sarcoma

Cancers ◽

10.3390/cancers13225837 ◽

2021 ◽

Vol 13 (22) ◽

pp. 5837

Author(s):

Changwu Wu ◽

Siming Gong ◽

Georg Osterhoff ◽

Nikolas Schopow

Keyword(s):

Soft Tissue ◽

Risk Score ◽

Cox Regression ◽

Soft Tissue Sarcomas ◽

Gene Signature ◽

Gene Set Enrichment Analysis ◽

The Cancer Genome Atlas ◽

Prognostic Models ◽

Malignant Tumours ◽

Free Survival

Soft tissue sarcomas (STS), a group of rare malignant tumours with high tissue heterogeneity, still lack effective clinical stratification and prognostic models. Therefore, we conducted this study to establish a reliable prognostic gene signature. Using 189 STS patients’ data from The Cancer Genome Atlas database, a four-gene signature including DHRS3, JRK, TARDBP and TTC3 was established. A risk score based on this gene signature was able to divide STS patients into a low-risk and a high-risk group. The latter had significantly worse overall survival (OS) and relapse free survival (RFS), and Cox regression analyses showed that the risk score is an independent prognostic factor. Nomograms containing the four-gene signature have also been established and have been verified through calibration curves. In addition, the predictive ability of this four-gene signature for STS metastasis free survival was verified in an independent cohort (309 STS patients from the Gene Expression Omnibus database). Finally, Gene Set Enrichment Analysis indicated that the four-gene signature may be related to some pathways associated with tumorigenesis, growth, and metastasis. In conclusion, our study establishes a novel four-gene signature and clinically feasible nomograms to predict the OS and RFS. This can help personalized treatment decisions, long-term patient management, and possible future development of targeted therapy.

Download Full-text

Immune‐related Long Noncoding RNA Signature for Patients with Cervical Cancer

10.21203/rs.3.rs-49217/v1 ◽

2020 ◽

Author(s):

Jianfeng Zheng ◽

Jinyi Tong ◽

Benben Cao ◽

Xia Zhang ◽

Zheng Niu

Keyword(s):

Cervical Cancer ◽

High Risk ◽

Risk Score ◽

Immune Cell ◽

Risk Group ◽

Risk Groups ◽

Low Risk ◽

Gene Set Enrichment Analysis ◽

Training Set ◽

Score Model

Abstract Background: Cervical cancer (CC) is a common gynecological malignancy for which prognostic and therapeutic biomarkers are urgently needed. The signature based on immune‐related lncRNAs(IRLs) of CC has never been reported. This study aimed to establish an IRL signature for patients with CC.Methods: The RNA-seq dataset was obtained from the TCGA, GEO, and GTEx database. The immune scores（IS）based on single-sample gene set enrichment analysis (ssGSEA) were calculated to identify the IRLs, which were then analyzed using univariate Cox regression analysis to identify significant prognostic IRLs. A risk score model was established to divide patients into low-risk and high-risk groups based on the median risk score of these IRLs. This was then validated by splitting TCGA dataset(n=304) into a training-set(n=152) and a valid-set(n=152). The fraction of 22 immune cell subpopulations was evaluated in each sample to identify the differences between low-risk and high-risk groups. Additionally, a ceRNA network associated with the IRLs was constructed.Results: A cohort of 326 CC and 21 normal tissue samples with corresponding clinical information was included in this study. Twenty-eight IRLs were collected according to the Pearson’s correlation analysis between immune score and lncRNA expression (P < 0.01). Four IRLs (BZRAP1-AS1, EMX2OS, ZNF667-AS1, and CTC-429P9.1) with the most significant prognostic values (P < 0.05) were identified which demonstrated an ability to stratify patients into low-risk and high-risk groups by developing a risk score model. It was observed that patients in the low‐risk group showed longer overall survival (OS) than those in the high‐risk group in the training-set, valid-set, and total-set. The area under the curve (AUC) of the receiver operating characteristic curve (ROC curve) for the four IRLs signature in predicting the one-, two-, and three-year survival rates were larger than 0.65. In addition, the low-risk and high-risk groups displayed different immune statuses in GSEA. These IRLs were also significantly correlated with immune cell infiltration. Conclusions: Our results showed that the IRL signature had a prognostic value for CC. Meanwhile, the specific mechanisms of the four-IRLs in the development of CC were ascertained preliminarily.

Download Full-text

Immune-related gene based prognostic signature for the risk stratification analysis of breast cancer

Current Bioinformatics ◽

10.2174/1574893616666211005110732 ◽

2021 ◽

Vol 16 ◽

Author(s):

Dongqing Su ◽

Qianzi Lu ◽

Yi Pan ◽

Yao Yu ◽

Shiyuan Wang ◽

...

Keyword(s):

Breast Cancer ◽

Regression Analysis ◽

Cancer Patients ◽

Risk Score ◽

Cox Regression ◽

Enrichment Analysis ◽

Gene Set Enrichment Analysis ◽

Breast Cancer Patients ◽

Cox Regression Analysis ◽

Score Model

Background: Breast cancer has plagued women for many years and caused many deaths around the world. Method: In this study, based on the weighted correlation network analysis, univariate Cox regression analysis and least absolute shrinkage and selection operator, 12 immune-related genes were selected to construct the risk score for breast cancer patients. The multivariable Cox regression analysis, gene set enrichment analysis and nomogram were also conducted in this study. Results: Good results were obtained in the survival analysis, enrichment analysis, multivariable Cox regression analysis and immune-related feature analysis. When the risk score model was applied in 22 breast cancer cohorts, the univariate Cox regression analysis demonstrated that the risk score model was significantly associated with overall survival in most of the breast cancer cohorts. Conclusion: Based on these results, we could conclude that the proposed risk score model may be a promising method, and may improve the treatment stratification of breast cancer patients in the future work.

Download Full-text

Identification and Validation of an Energy Metabolism-Related lncRNA-mRNA Signature for Lower-Grade Glioma

BioMed Research International ◽

10.1155/2020/3708231 ◽

2020 ◽

Vol 2020 ◽

pp. 1-14

Author(s):

Jingwei Zhao ◽

Le Wang ◽

Bo Wei

Keyword(s):

Energy Metabolism ◽

Risk Score ◽

Cox Regression ◽

The Cancer Genome Atlas ◽

Coexpression Network ◽

Normal Brain ◽

Lower Grade ◽

Prognostic Signature ◽

Lower Grade Glioma ◽

Genome Atlas

Energy metabolic processes play important roles for tumor malignancy, indicating that related protein-coding genes and regulatory upstream genes (such as long noncoding RNAs (lncRNAs)) may represent potential biomarkers for prognostic prediction. This study will develop a new energy metabolism-related lncRNA-mRNA prognostic signature for lower-grade glioma (LGG) patients. A GSE4290 dataset obtained from Gene Expression Omnibus was used for screening the differentially expressed genes (DEGs) and lncRNAs (DELs). The Cancer Genome Atlas (TCGA) dataset was used as the prognosis training set, while the Chinese Glioma Genome Atlas (CGGA) was for the validation set. Energy metabolism-related genes were collected from the Molecular Signatures Database (MsigDB), and a coexpression network was established between energy metabolism-related DEGs and DELs to identify energy metabolism-related DELs. Least absolute shrinkage and selection operator (LASSO) analysis was performed to filter the prognostic signature which underwent survival analysis and nomogram construction. A total of 1613 DEGs and 37 DELs were identified between LGG and normal brain tissues. One hundred and ten DEGs were overlapped with energy metabolism-related genes. Twenty-seven DELs could coexpress with 67 metabolism-related DEGs. LASSO regression analysis showed that 9 genes in the coexpression network were the optimal signature and used to construct the risk score. Kaplan-Meier curve analysis showed that patients with a high risk score had significantly worse OS than those with a low risk score (TCGA: HR=3.192, 95%CI=2.182‐4.670; CGGA: HR=1.922, 95%CI=1.431‐2.583). The predictive accuracy of the risk score was also high according to the AUC of the ROC curve (TCGA: 0.827; CGGA: 0.806). Multivariate Cox regression analyses revealed age, IDH1 mutation, and risk score as independent prognostic factors, and thus, a prognostic nomogram was established based on these three variables. The excellent prognostic performance of the nomogram was confirmed by calibration and discrimination analyses. In conclusion, our findings provided a new biomarker for the stratification of LGG patients with poor prognosis.

Download Full-text