Dataset Information

Construction of an Immune-Related Six-lncRNA Signature to Predict the Outcomes, Immune Cell Infiltration, and Immunotherapy Response in Patients With Hepatocellular Carcinoma.

ABSTRACT:

Background

Hepatocellular carcinoma (HCC) is one of the world's most lethal malignant tumors with a poor prognosis. Growing evidence has been demonstrating that immune-related long non-coding RNAs (lncRNAs) are relevant to the tumor microenvironment (TME) and can help assess the effects of immunotherapy and evaluate one's prognosis. This study aims to identify an immune-related lncRNA signature for the prospective assessment of the immunotherapy and prognosis in HCC.

Method

We downloaded HCC RNA-seq data and clinical information from The Cancer Genome Atlas (TCGA) project database. We first used ESTIMATE to evaluate the TME. Then, we conducted a cox regression analysis to construct a prognostic signature and the riskScore. We then applied the univariate Cox regression, multivariate Cox regression, principal components analysis (PCA), receiver operating characteristic (ROC) curve, and stratification analyses to confirm our previous assessments. Afterward, we employed a gene set enrichment analysis (GSEA) to explore the biological processes and pathways. Besides, we used CIBERSORT to estimate the abundance of tumor-infiltrating immune cells (TIICs). Furthermore, we investigated the relationship between the immune-related lncRNA signature and immune checkpoint genes. Finally, we used the quantitative real-time polymerase chain reaction (qRT-PCR) assays to demonstrate the expression of the six lncRNAs.

Results

We identified six immune-related lncRNAs - MSC-AS1, AC145207.5, SNHG3, AL365203.2, AL031985.3, NRAV - which show the ability to stratify patients into high-risk and low-risk groups with significantly different survival rates. The univariate Cox regression, multivariate Cox regression, ROC, and stratification analyses confirmed that the immune-related six-lncRNA signature was a novel independent prognostic factor in HCC patients. The high-risk group and low-risk group illustrated contrasting distributions in PCA. The GSEA suggested that the immune-related six-lncRNA signature was involved in the immune-related biological processes and pathways. Besides, the immune-related six-lncRNA signature was associated with the infiltration of immune cells. Furthermore, it was linked with the expression of critical immune genes and could predict immunotherapy's clinical response. Finally, the qRT-PCR demonstrated that the six lncRNAs were significantly differentially expressed in HCC cell lines and normal hepatic cell lines.

Conclusion

In summary, we identified an immune-related six-lncRNA signature that can predict the outcomes, immune cell infiltration, and immunotherapy response in patients with hepatocellular carcinoma.

SUBMITTER: Zhou P

PROVIDER: S-EPMC8283691 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundColon cancer remains one of the most common malignancies across the world. Thus far, a biomarker, which can comprehensively predict the survival outcomes, clinical characteristics, and therapeutic sensitivity, is still lacking.MethodsWe leveraged transcriptomic data of colon cancer from the existing datasets and constructed immune-related lncRNA (irlncRNA) pairs. After integrating with clinical survival data, we performed differential analysis and identified 11 irlncRNAs signature using Lasso regression analysis. We next plotted the 1-, 5-, and 10-year curve lines of receiver operating characteristics, calculated the areas under the curve, and recognized the optimal cutoff point. Then, we validated the pair-risk model in terms of the survival outcomes of the patients involved. Moreover, we tested the reliability of the model for predicting tumor aggressiveness and therapeutic susceptibility of colon cancer. Additionally, we reemployed the 11 of irlncRNAs involved in the pair-risk model to construct an expression-risk model to predict the prognostic outcomes of the patients involved.ResultsWe recognized a total of 377 differentially expressed irlncRNAs (DEirlcRNAs), including 28 low-expressed and 349 high-expressed irlncRNAs in colon cancer patients. After performing a univariant Cox analysis, we identified 115 risk irlncRNAs that were significantly correlated with survival outcomes of patients involved. By taking the overlap of the DEirlcRNAs and the risk irlncRNAs, we ultimately recognized 55 irlncRNAs as core irlncRNAs. Then, we established a Cox HR model (pair-risk model) as well as an expression HR model (exp-risk model) based on 11 of the 55 core irlncRNAs. We found that both of the two models significantly outperformed the commonly used clinical characteristics, including age, T, N, and M stages when predicting survival outcomes. Moreover, we validated the pair-risk model as a potential tool for studying the tumor microenvironment of colon cancer and drug susceptibility. Additionally, we noticed that combinational use of the pair-risk model and the exp-risk model yielded a more robust approach for predicting the survival outcomes of patients with colon cancer.ConclusionsWe recognized 11 irlncRNAs and created a pair-risk model and an exp-risk model, which have the potential to predict clinical characteristics of colon cancer, either solely or conjointly.

Project description:BackgroundGenomic instability (GI) plays a crucial role in the development of various cancers including hepatocellular carcinoma. Hence, it is meaningful for us to use long non-coding RNAs related to genomic instability to construct a prognostic signature for patients with HCC.MethodsCombining the lncRNA expression profiles and somatic mutation profiles in The Cancer Genome Atlas database, we identified GI-related lncRNAs (GILncRNAs) and obtained the prognosis-related GILncRNAs through univariate regression analysis. These lncRNAs obtained risk coefficients through multivariate regression analysis for constructing GI-associated lncRNA signature (GILncSig). ROC curves were used to evaluate signature performance. The International Cancer Genomics Consortium (ICGC) cohort, and in vitro experiments were used for signature external validation. Immunotherapy efficacy, tumor microenvironments, the half-maximal inhibitory concentration (IC50), and immune infiltration were compared between the high- and low-risk groups with TIDE, ESTIMATE, pRRophetic, and ssGSEA program.ResultsFive GILncRNAs were used to construct a GILncSig. It was confirmed that the GILncSig has good prognostic evaluation performance for patients with HCC by drawing a time-dependent ROC curve. Patients were divided into high- and low-risk groups according to the GILncSig risk score. The prognosis of the low-risk group was significantly better than that of the high-risk group. Independent prognostic analysis showed that the GILncSig could independently predict the prognosis of patients with HCC. In addition, the GILncSig was correlated with the mutation rate of the HCC genome, indicating that it has the potential to measure the degree of genome instability. In GILncSig, LUCAT1 with the highest risk factor was further validated as a risk factor for HCC in vitro. The ESTIMATE analysis showed a significant difference in stromal scores and ESTIMATE scores between the two groups. Multiple immune checkpoints had higher expression levels in the high-risk group. The ssGSEA results showed higher levels of tumor-antagonizing immune cells in the low-risk group compared with the high-risk group. Finally, the GILncSig score was associated with chemotherapeutic drug sensitivity and immunotherapy efficacy of patients with HCC.ConclusionOur research indicates that GILncSig can be used for prognostic evaluation of patients with HCC and provide new insights for clinical decision-making and potential therapeutic strategies.

Project description:Ferroptosis is an iron-dependent cell death process that plays important regulatory roles in the occurrence and development of cancers, including hepatocellular carcinoma (HCC). Moreover, the molecular events surrounding aberrantly expressed long non-coding RNAs (lncRNAs) that drive HCC initiation and progression have attracted increasing attention. However, research on ferroptosis-related lncRNA prognostic signature in patients with HCC is still lacking. In this study, the association between differentially expressed lncRNAs and ferroptosis-related genes, in 374 HCC and 50 normal hepatic samples obtained from The Cancer Genome Atlas (TCGA), was evaluated using Pearson's test, thereby identifying 24 ferroptosis-related differentially expressed lncRNAs. The least absolute shrinkage and selection operator (LASSO) algorithm and Cox regression model were used to construct and validate a prognostic risk score model from both TCGA training dataset and GEO testing dataset (GSE40144). A nine-lncRNA-based signature (CTD-2033A16.3, CTD-2116N20.1, CTD-2510F5.4, DDX11-AS1, LINC00942, LINC01224, LINC01231, LINC01508, and ZFPM2-AS1) was identified as the ferroptosis-related prognostic model for HCC, independent of multiple clinicopathological parameters. In addition, the HCC patients were divided into high-risk and low-risk groups according to the nine-lncRNA prognostic signature. The gene set enrichment analysis enrichment analysis revealed that the lncRNA-based signature might regulate the HCC immune microenvironment by interfering with tumor necrosis factor α/nuclear factor kappa-B, interleukin 2/signal transducers and activators of transcription 5, and cytokine/cytokine receptor signaling pathways. The infiltrating immune cell subtypes, such as resting memory CD4(+) T cells, follicular helper T cells, regulatory T cells, and M0 macrophages, were all significantly different between the high-risk group and the low-risk group as indicated in Spearman's correlation analysis. Moreover, a substantial increase in the expression of B7H3 immune checkpoint molecule was found in the high-risk group. Our findings provided a promising insight into ferroptosis-related lncRNAs in HCC and a personalized prediction tool for prognosis and immune responses in patients.

Project description:Increasing evidence has demonstrated that pyroptosis, a type of inflammatory programmed cell death, plays an important role in the pathogenesis and progression of gastric cancer. However, it remains unclear whether pyroptosis-related long non-coding RNAs (lncRNAs) can be used to predict the diagnosis and prognosis of gastric adenocarcinoma. This study aimed to evaluate and test the role of the lncRNA signature associated with pyroptosis as a prognostic tool for stomach adenocarcinoma (STAD) and to ascertain their immune value. Relative RNA-sequencing data were extracted from The Cancer Genome Atlas database (TCGA), and data preprocessing was performed for STAD. Pearson correlation analysis was used to determine whether lncRNAs were significantly correlated with pyroptosis based on 23 genes related to pyroptosis. Univariate Cox regression and least absolute shrinkage and selection operator(LASSO) analyses were both adopted to select features and establish the pyroptosis-related lncRNA (PRL) prognostic signature. Kaplan-Meier(KM) survival analysis of the different risk groups was conducted according to the risk scores. We further examined the functional enrichment, tumor microenvironment, and landscape of mutation status among the different risk groups, and these analyses further explained the reasons for the differences in the prediction as well as survival value of the different risk groups. Four lncRNAs, including HAND2-AS1, LINC01354, RP11-276H19.1, and PGM5-AS1, were involved in the PRL signature and used to split STAD patients into two risk groups. Overall survival time(OS) was significantly higher in the low-risk group than in the high-risk group in both the training and validation groups. Functional enrichment analysis was further employed to analyze differentially expressed genes in high- and low-risk groups to identify potential molecular functions and pathways associated with pyroptosis in the gastric cancer microenvironment. Protein-protein interaction (PPI) and Friends analysis identified hub genes that may play a key role in differentially expressed genes in high- and low-risk groups. In addition, there were remarkable discrepancies between the different risk groups in the tumor stage (P < 0.01) and histologic grade (P < 0.05). Furthermore, drug-susceptibility testing indicated potential sensitive chemotherapeutic drugs for each risk group. This study is the first to establish and validate STAD-associated PRLs that can effectively guide the prognosis and the immune microenvironment in STAD patients and provide evidence for the development of molecularly targeted therapies related to pyroptosis.

Project description:Breast cancer (BC) is among the most universal malignant tumors in women worldwide. Aging is a complex phenomenon, caused by a variety of factors, that plays a significant role in tumor development. Consequently, it is crucial to screen for prognostic aging-related long non-coding RNAs (lncRNAs) in BC. The BC samples from the breast-invasive carcinoma cohort were downloaded from The Cancer Genome Atlas (TCGA) database. The differential expression of aging-related lncRNAs (DEarlncRNAs) was screened by Pearson correlation analysis. Univariate Cox regression, LASSO-Cox analysis, and multivariate Cox analysis were performed to construct an aging-related lncRNA signature. The signature was validated in the GSE20685 dataset from the Gene Expression Omnibus (GEO) database. Subsequently, a nomogram was constructed to predict survival in BC patients. The accuracy of prediction performance was assessed through the time-dependent receiver operating characteristic (ROC) curves, Kaplan-Meier analysis, principal component analyses, decision curve analysis, calibration curve, and concordance index. Finally, differences in tumor mutational burden, tumor-infiltrating immune cells, and patients' response to chemotherapy and immunotherapy between the high- and low-risk score groups were explored. Analysis of the TCGA cohort revealed a six aging-related lncRNA signature consisting of MCF2L-AS1, USP30-AS1, OTUD6B-AS1, MAPT-AS1, PRR34-AS1, and DLGAP1-AS1. The time-dependent ROC curve proved the optimal predictability for prognosis in BC patients with areas under curves (AUCs) of 0.753, 0.772, and 0.722 in 1, 3, and 5 years, respectively. Patients in the low-risk group had better overall survival and significantly lower total tumor mutational burden. Meanwhile, the high-risk group had a lower proportion of tumor-killing immune cells. The low-risk group could benefit more from immunotherapy and some chemotherapeutics than the high-risk group. The aging-related lncRNA signature can provide new perspectives and methods for early BC diagnosis and therapeutic targets, especially tumor immunotherapy.