Dataset Information

A Prognostic Signature of Glycolysis-Related Long Noncoding RNAs for Molecular Subtypes in the Tumor Immune Microenvironment of Lung Adenocarcinoma.

ABSTRACT:

Purpose

Long noncoding RNAs (lncRNAs) and glycolysis regulate multiple types of cancer. However, the prognostic roles and biological functions of glycolysis-related lncRNAs in lung adenocarcinoma (LUAD) remain unclear. In this study, we investigated the role of glycolysis-related lncRNAs in LUAD.

Patients and methods

We retrieved glycolysis-related genes from the Molecular Signatures Database and screened for prognostic glycolysis-related lncRNAs from The Cancer Genome Atlas.

Results

We identified three LUAD subtypes (clusters 1-3) by univariate Cox regression analysis and consensus clustering. Patients in cluster 1 had the best overall survival rates. Immune, stromal, and cytolytic-activity scores were the highest in cluster 1. The expression of immune checkpoint molecules (programmed cell death protein 1 and cytotoxic T-lymphocyte-associated protein 4) and other immune-related indicators was the highest in cluster 1, whereas that of epithelial cell biomarkers (Cadherin 1, Cadherin 2, and MET) was the lowest. Therefore, patients in cluster 1 may benefit from immunotherapy. Lasso-Cox regression and multivariate Cox regression analyses were used to select nine lncRNAs to build a robust prognostic model of LUAD. The area under the curve classifier values and a nomogram performed well in predicting survival times for patients with LUAD. The expression levels of nine lncRNAs were validated by quantitative reverse transcriptase-polymerase chain reaction analysis, and most of these lncRNAs were significantly related to immune-related mRNAs. Gene set enrichment analysis revealed that the high-risk group was enriched for cell cycle-related pathways and the low-risk group was enriched for pathways associated with immunity or immune-related diseases.

Conclusion

The LUAD subtypes and prognostic model developed here may help in clinical risk stratification, prognosis management, and treatment decisions for patients with LUAD.

SUBMITTER: Li N

PROVIDER: S-EPMC8637177 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Purpose: Our research developed immune-related long noncoding RNAs (lncRNAs) for risk stratification in cervical cancer (CC) and explored factors of prognosis, inflammatory microenvironment infiltrates, and chemotherapeutic therapies. Methods: The RNA-seq data and clinical information of CC were collected from the TCGA TARGET GTEx database and the TCGA database. lncRNAs and immune-related signatures were obtained from the GENCODE database and the ImPort database, respectively. We screened out immune-related lncRNA signatures through univariate Cox, LASSO, and multivariate Cox regression methods. We established an immune-related risk model of hub immune-related lncRNAs to evaluate whether the risk score was an independent prognostic predictor. The xCell and CIBERSORTx algorithms were employed to appraise the value of risk scores which are in competition with tumor-infiltrating immune cell abundances. The estimation of tumor immunotherapy response through the TIDE algorithm and prediction of innovative recommended medications on the target to immune-related risk model were also performed on the basis of the IC50 predictor. Results: We successfully established six immune-related lncRNAs (AC006126.4, EGFR-AS1, RP4-647J21.1, LINC00925, EMX2OS, and BZRAP1-AS1) to carry out prognostic prediction of CC. The immune-related risk model was constructed in which we observed that high-risk groups were strongly linked with poor survival outcomes. Risk scores varied with clinicopathological parameters and the tumor stage and were an independent hazard factor that affect prognosis of CC. The xCell algorithm revealed that hub immune-related signatures were relevant to immune cells, especially mast cells, DCs, megakaryocytes, memory B cells, NK cells, and Th1 cells. The CIBERSORTx algorithm revealed an inflammatory microenvironment where naive B cells (p < 0.01), activated dendritic cells (p < 0.05), activated mast cells (p < 0.0001), CD8+ T cells (p < 0.001), and regulatory T cells (p < 0.01) were significantly lower in the high-risk group, while macrophages M0 (p < 0.001), macrophages M2 (p < 0.05), resting mast cells (p < 0.0001), and neutrophils (p < 0.01) were highly conferred. The result of TIDE indicated that the number of immunotherapy responders in the low-risk group (124/137) increased significantly (p = 0.00000022) compared to the high-risk group (94/137), suggesting that the immunotherapy response of CC patients was completely negatively correlated with the risk scores. Last, we compared differential IC50 predictive values in high- and low-risk groups, and 12 compounds were identified as future treatments for CC patients. Conclusion: In this study, six immune-related lncRNAs were suggested to predict the outcome of CC, which is beneficial to the formulation of immunotherapy.

Project description:Purpose: The purpose of this study was to construct a novel risk scoring model with prognostic value that could elucidate tumor immune microenvironment of hepatocellular carcinoma (HCC). Samples and methods: Data were obtained through The Cancer Genome Atlas (TCGA) database. Univariate Cox analysis, least absolute shrinkage and selection operator (LASSO) analysis, and multivariate Cox analysis were carried out to screen for glycolysis-related long noncoding RNAs (lncRNAs) that could provide prognostic value. Finally, we established a risk score model to describe the characteristics of the model and verify its prediction accuracy. The receiver operating characteristic (ROC) curves of 1, 3, and 5 years of overall survival (OS) were depicted with risk score and some clinical features. ESTIMATE algorithm, single-sample gene set enrichment analysis (ssGSEA), and CIBERSORT analysis were employed to reveal the characteristics of tumor immune microenvironment in HCC. The nomogram was drawn by screening indicators with high prognostic accuracy. The correlation of risk signature with immune infiltration and immune checkpoint blockade (ICB) therapy was analyzed. After enrichment of related genes, active behaviors and pathways in high-risk groups were identified and lncRNAs related to poor prognosis were validated in vitro. Finally, the impact of MIR4435-2HG upon ICB treatment was uncovered. Results: After screening through multiple steps, four glycolysis-related lncRNAs were obtained. The risk score constructed with the four lncRNAs was found to significantly correlate with prognosis of samples. From the ROC curve of samples with 1, 3, and 5 years of OS, two indicators were identified with high prognostic accuracy and were used to draw a nomogram. Besides, the risk score significantly correlated with immune score, immune-related signature, infiltrating immune cells (i.e. B cells, etc.), and ICB key molecules (i.e. CTLA4,etc.). Gene enrichment analysis indicated that multiple biological behaviors and pathways were active in the high-risk group. In vitro validation results showed that MIR4435-2HG was highly expressed in the two cell lines, which had a significant impact on the OS of samples. Finally, we corroborated that MIR4435-2HG had intimate relationship with ICB therapy in hepatocellular carcinoma. Conclusion: We elucidated the crucial role of risk signature in immune cell infiltration and immunotherapy, which might contribute to clinical strategies and clinical outcome prediction of HCC.

Project description:AbstractReliable biomarkers are of great significance for the treatment and diagnosis of hepatocellular carcinoma (HCC). This study identified potential prognostic epithelial-mesenchymal transition related lncRNAs (ERLs) by the cancer genome atlas (TCGA) database and bioinformatics.The differential expression of long noncoding RNA (lncRNA) was obtained by analyzing the lncRNA data of 370 HCC samples in TCGA. Then, Pearson correlation analysis was carried out with EMT related genes (ERGs) from molecular signatures database. Combined with the univariate Cox expression analysis of the total survival rate of hepatocellular carcinoma (HCC) patients, the prognostic ERLs were obtained. Then use "step" function to select the optimal combination of constructing multivariate Cox expression model. The expression levels of ERLs in HCC samples were verified by real-time quantitative polymerase chain reaction.Finally, we identified 5 prognostic ERLs (AC023157.3, AC099850.3, AL031985.3, AL365203.2, CYTOR). The model showed that these prognostic markers were reliable independent predictors of risk factors (P value <.0001, hazard ratio [HR] = 2.400, 95% confidence interval [CI] = 1.667-3.454 for OS). In the time-dependent receiver operating characteristic analysis, this prognostic marker is a good predictor of HCC survival (area under the curve of 1 year, 2 years, 3 years, and 5 years are 0.754, 0.720, 0.704, and 0.662 respectively). We analyzed the correlation of clinical characteristics of these prognostic markers, and the results show that this prognostic marker is an independent factor that can predict the prognosis of HCC more accurately. In addition, by matching with the Molecular Signatures Database, we obtained 18 ERLs, and then constructed the HCC prognosis model and clinical feature correlation analysis using 5 prognostic ERLs. The results show that these prognostic markers have reliable independent predictive value. Bioinformatics analysis showed that these prognostic markers were involved in the regulation of EMT and related functions of tumor occurrence and migration.Five prognostic types of ERLs identified in this study can be used as potential biomarkers to predict the prognosis of HCC.

Project description:BackgroundOsteosarcoma is the most prevalent bone cancer that affects young adults and adolescents. It is the most frequent malignancy of the bone. In spite of the fact that complete surgical resection and chemotherapy have increased the overall survival of osteosarcoma patients considerably, the prognosis remains dismal in patients with recurring and/or metastasized osteosarcoma. Thus, finding predictive biomarkers representing osteosarcoma's biological variability may result in more effective treatment for osteosarcoma patients.MethodsIn this research, RNA data and clinical information were obtained from TARGET database. The risk score was calculated using a technique that incorporated both univariate and multivariate Cox regression. A variety of statistical methods were employed to assess the risk score's accuracy. These included ROC curves, nomograms, and Kaplan-Meier curves. Following that, bioinformatics studies were carried out in order to investigate the possible biological processes that influence the prognosis of osteosarcoma patients. GSEA was used to investigate the variations in pathway enrichment among the different groups of genes. To examine the disparities in the immune microenvironment, the analytical methods CIBERSORT and ssGSEA were employed.ResultsWe discovered three differentially expressed lncRNAs (RPARP-AS1, AC009159.3, and AC124312.3) that are linked to osteosarcoma prognosis. Kaplan-Meier analysis showed the presence of a signature of high-risk lncRNAs linked with a poor prognosis for osteosarcoma. Furthermore, the AUC of the lncRNAs signature was 0.773, indicating that they are useful in predicting osteosarcoma prognosis in certain cases. In predicting osteosarcoma prognosis, our risk assessment approach outperformed conventional clinicopathological characteristics. In the high-risk group of people, GSEA showed the presence of tumor-related pathways as well as immune-related pathways. Furthermore, TARGET revealed that immune-related functions such as checkpoint, T-cell coinhibition, and costimulation were significantly different between the high-risk and low-risk groups. LAIR1, LAG3, CD44, and CD22, as well as other immune checkpoints, were shown to be expressed differentially across the two risk groups.ConclusionThis study established that pyroptosis-derived lncRNAs had a significant predictive value for osteosarcoma patients' survival, indicating that they may be a viable target for future therapy.

Project description:Cholangiocarcinoma (CCA) is featured by common occurrence and poor prognosis. Autophagy is a biological process that has been extensively involved in the progression of tumors. Long noncoding RNAs (lncRNAs) have been discovered to be critical in diagnosing and predicting various tumors. It may be valuable to elaborate autophagy-related lncRNAs (ARlncRNAs) in CCA, and indeed, there are still few studies concerning the role of ARlncRNAs in CCA. Here, a prognostic ARlncRNA signature was constructed to predict the survival outcome of CCA patients. Through identification, three differentially expressed ARlncRNAs (DEARlncRNAs), including CHRM3.AS2, MIR205HG, and LINC00661, were screened and were considered predictive signatures. Furthermore, the overall survival (OS) of patients with high-risk scores was significantly lower than that of patients with low scores. Interestingly, the risk score was an independent factor for the OS of patients with CCA. Moreover, receiver operating characteristic (ROC) curve analysis showed that the screened and constructed prognosis signature for 1 year (AUC = 0.884), 3 years (AUC =0.759), and 5 years (AUC = 0.788) presented a high score of accuracy in predicting OS of CCA patients. Gene set enrichment analysis (GSEA) revealed that the three DEARlncRNAs were significantly enriched in CCA-related signaling pathways, including "pathways of basal cell carcinoma", "glycerolipid metabolism", etc. Quantitative real-time PCR (qRT-PCR) showed that expressions of CHRM3.AS2, MIR205HG, and LINC00661 were higher in CCA tissues than those in normal tissues, similar to the trends detected in the CCA dataset. Furthermore, Pearson's analysis reported an intimate correlation of the risk score with immune cell infiltration, indicating a predictive value of the signature for the efficacy of immunotherapy. In addition, the screened lncRNAs were found to have the ability to modulate the expression of mRNAs by interacting with miRNAs based on the established lncRNA-miRNA-mRNA network. In conclusion, our study develops a novel nomogram with good reliability and accuracy to predict the OS of CCA patients, providing a significant guiding value for developing tailored therapy for CCA patients.

Project description:BackgroundPatients with advanced clear cell renal cell carcinoma (ccRCC) have a poor prognosis and lack effective prognostic biomarkers. N6-methyladenosine-related lncRNAs (m6A-related long noncoding RNAs [lncRNAs]) have been confirmed to be associated with the development of multiple tumors, but its role in ccRCC is not clear.MethodsGene expression data and clinical information of ccRCC patients were extracted from The Cancer Genome Atlas Database. The prognostic m6A-related lncRNAs were obtained by Pearson's correlation analysis and univariate Cox regression analysis. Afterward, the cluster classification and its correlation with prognosis, clinical characteristics, and immunity were analyzed. LASSO regression was used to establish the prognostic risk model. The predictive performance of the prognostic model was evaluated and validated by survival analysis and receiver operating characteristic curve analysis, et al. The expression of immune checkpoints and immune cell infiltration in patients with different risks were systematically analyzed.ResultsA total of 27 prognostic m6A-related lncRNAs were identified. These m6A-related lncRNAs were differentially expressed between tumor and normal tissues. Among them, 24 high-risk m6A-related lncRNAs were overexpressed in Cluster 2 and correlated with poor prognosis, low stromal score, high expression of immune checkpoints, and immunosuppressive cells infiltration. Based upon, a prognostic risk model composed of seven m6A-related lncRNAs was constructed. After a series of analyses, it was proved that this model had good sensitivity and specificity, and could predict the prognosis of patients with different clinical stratification. The expression of PD-1, PD-L1, CTLA-4, LAG-3, TIM-3, and TIGIT were significantly increased in the high-risk patients, and there was a correlation between the risk score and immune cell infiltration.ConclusionsThe seven m6A-related lncRNAs prognostic risk signature showed reliable prognostic predictive power for ccRCC and was associated with the expression of immune checkpoints and immune cell infiltration. This seven m6A-related lncRNAs signature will be helpful in managing ccRCC and guiding individualized immunotherapy.

Project description:BackgroundRecently, a newly programmed cell death has been discovered, namely cuproptosis. It is considered a novel copper-dependent cell death model. Long non-coding RNA (lncRNA) influence the prognosis of bladder cancer. In this study, we established a scoring system based on 7 cuproptosis-related lncRNA to predict the prognosis and immune landscape of bladder cancer (BCa).MethodGene expression and clinical data of 431 tissues were downloaded from The Cancer Genome Atlas (TCGA), including 19 normal samples and 419 cancer samples. All samples were randomly categorized into train and test cohorts. Cuproptosis-related lncRNA were distinguished. Then we conduct univariate COX and multivariate COX regression, paralleled with LASSO regression to cultivate a cuproptosis-related lncRNA risk model. Kaplan-Meier curves, scatter diagram, C-index, ROC curves, nomogram, PCA analysis and univariate and multivariate Cox regression were used to test the accuracy of risk model and to predict patient survival. Additional, gene mutation status between high- and low-risk groups was calculated.GO and KEGG were used to access the DEGs (different expression genes)-related pathway.The ssGSEA and ESTIMATE algorithms were used to assess the immune function in different tumor samples. Besides, patient's response to immunotherapy and drug susceptibility were also been estimated.Results7 cuproptosis-related lncRNA (LINC01184, LINC00513, LINC02443, SMARCA5-AS1, BDNF-AS, SOD2-OT1, HYI-AS1) were selected to construct the risk model in the train cohort. This model can well predict the overall survival (OS) in test group and entire cohort with different stage. Despite no significant different is observed in gene mutation between high- and low-risk group, different immune infiltration, different survival and sensitivity to drugs are discovered.ConclusionWe established a novel cuproptosis-related lncRNA risk model which can predict the outcome and immunotherapy response with satisfactory predictive effects. This risk model can provide a new insight into prognostic evaluation and may have potential to guide comprehensive treatment in bladder cancer.

Project description:Background: Immune and stromal cells in the tumor microenvironment (TME) significantly contribute to the prognosis of lung adenocarcinoma; however, the TME-related immune prognostic signature is unknown. The aim of this study was to develop a novel immune prognostic model of the TME in lung adenocarcinoma. Methods: First, the immune and stromal scores among lung adenocarcinoma patients were determined using the ESTIMATE algorithm in accordance with The Cancer Genome Atlas (TCGA) database. Differentially expressed immune-related genes (IRGs) between high and low immune/stromal score groups were analyzed, and a univariate Cox regression analysis was performed to identify IRGs significantly correlated with overall survival (OS) among patients with lung adenocarcinoma. Furthermore, a least absolute shrinkage and selection operator (LASSO) regression analysis was performed to generate TME-related immune prognostic signatures. Gene set enrichment analysis was performed to analyze the mechanisms underlying these immune prognostic signatures. Finally, the functions of hub IRGs were further analyzed to delineate the potential prognostic mechanisms in comprehensive TCGA datasets. Results: In total, 702 intersecting differentially expressed IRGs (589 upregulated and 113 downregulated) were screened. Univariate Cox regression analysis revealed that 58 significant differentially expressed IRGs were correlated with patient prognosis in the training cohort, of which three IRGs (CLEC17A, INHA, and XIRP1) were identified through LASSO regression analysis. A robust prognostic model was generated on the basis of this three-IRG signature. Furthermore, functional enrichment analysis of the high-risk-score group was performed primarily on the basis of metabolic pathways, whereas analysis of the low-risk-score group was performed primarily on the basis of immunoregulation and immune cell activation. Finally, hub IRGs CLEC17A, INHA, and XIRP1 were considered novel prognostic biomarkers for lung adenocarcinoma. These hub genes had different mutation frequencies and forms in lung adenocarcinoma and participated in different signaling pathways. More importantly, these hub genes were significantly correlated with the infiltration of CD4+ T cells, CD8+ T cells, macrophages, B cells, and neutrophils. Conclusions: The robust novel TME-related immune prognostic signature effectively predicted the prognosis of patients with lung adenocarcinoma. Further studies are required to further elucidate the regulatory mechanisms of these hub IRGs in the TME and to develop new treatment strategies.