Dataset Information

Screening and verification of prognostic lncRNA markers related to immune infiltration in the metastasis of osteosarcoma.

ABSTRACT:

Background

We sought to screen and verify the long non-coding ribonucleic acids (lncRNAs) related to immune infiltration in metastatic osteosarcoma (OS).

Methods

We downloaded the RNA-sequencing expression data from The Cancer Genome Atlas (TCGA) database as the training data set. We downloaded the GSE39055 data set from the National Center for Biotechnology Information, Gene Expression Omnibus as the validation data set. The least absolute shrinkage and selection operator (LASSO) regression algorithm was used to screen the optimized lncRNA combinations. Kaplan-Meier curves were used to evaluate the associations between the lncRNAs and actual prognosis. The independent survival prognosis clinical factors were obtained by univariate and multivariate Cox analyses. A nomogram was established to explore the correlation between survival rate and risk information. The Tumor IMmune Estimation Resource was applied to estimate the composition of 6 subtypes of immune infiltration cells.

Results

In total, 1,398 lncRNAs and 14,631 messenger RNAs were screened from TCGA data set, and divided into the low and high immunity groups. The Estimation of STromal and Immune cells in MAlignant Tumour tissues using Expression data (ESTIMATE) scores differed significantly between the samples in the two groups. Additionally, 5 optimized lncRNA combinations were obtained using the LASSO algorithm. Risk factors including age, metastatic tumor, and risk-score (RS) were related to the prognosis of OS patients. The survival rates predicted by the nomogram model were consistent with the actual survival rates of OS patients. Finally, we found that RS was negatively correlated with the proportion of immune cells.

Conclusions

In total, 5 feature lncRNAs were identified as novel biomarkers for OS. Next, a RS nomogram model was constructed based on the 5 identified lncRNAs. This model predicted the survival rates and prognoses of OS patients well.

SUBMITTER: He C

PROVIDER: S-EPMC9552095 | biostudies-literature | 2022 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Screening and verification of prognostic lncRNA markers related to immune infiltration in the metastasis of osteosarcoma.

He Chuan C Wang Xiaolong X Ni Linying L Song Chunyu C Mai Wei W

Translational cancer research 20220901 9

<h4>Background</h4>We sought to screen and verify the long non-coding ribonucleic acids (lncRNAs) related to immune infiltration in metastatic osteosarcoma (OS).<h4>Methods</h4>We downloaded the RNA-sequencing expression data from The Cancer Genome Atlas (TCGA) database as the training data set. We downloaded the GSE39055 data set from the National Center for Biotechnology Information, Gene Expression Omnibus as the validation data set. The least absolute shrinkage and selection operator (LASSO) ...[more]

PMID: 36237238

Similar Datasets

Project description:BackgroundIncreasing evidence has shown that hypoxia microenvironment relates to tumor initiation and progression. However, no studies focus on the application of hypoxia-associated genes in predicting osteosarcoma patients' prognosis. This research aims to identify the hypoxia-associated genes related to osteosarcoma metastasis and construct a gene signature to predict osteosarcoma prognosis.MethodsThe differentially expressed messenger RNAs (DEmRNAs) related to osteosarcoma metastasis were identified from Therapeutically Applicable Research to Generate Effective Treatments (Target) database. Univariate and multivariate cox regression analyses were performed to develop the hypoxia-associated prognostic signature. The Kaplan-Meier (KM) survival analyses of patients with high and low hypoxia risk scores were conducted. The nomogram was constructed and the gene signature was validated in the external Gene Expression Omnibus (GEO) cohort. Single-sample gene set enrichment analysis (ssGSEA) was conducted to investigate the relationships between immune infiltration and gene signature.ResultsTwo genes, including decorin (DCN) and prolyl 4-hydroxylase subunit alpha 1 (P4HA1), were involved in the hypoxia-associated gene signature. In training and testing datasets, patients with high-risk scores showed lower survival rates and the gene signature was identified as the independent prognostic factor. Receiver operating characteristic (ROC) curves demonstrated the robustness of signature. Functional analyses of DEmRNAs among high- and low-risk groups revealed that immune-associated functions and pathways were significantly enriched. Furthermore, ssGSEA showed that five immune cells (DCs, macrophages, neutrophils, pDCs, and TIL) and three immune features (CCR, APC co inhibition, and Check-point) were down-regulated in the high-risk group.ConclusionThe current study established and validated a novel hypoxia-associated gene signature in osteosarcoma. It could act as a prognostic biomarker and serve as therapeutic guidance in clinical applications.

Project description:Osteosarcoma is the most common primary malignancy of the bones, and is associated with a high rate of metastasis and a poor prognosis. A tight association between the tumor microenvironment (TME) and osteosarcoma metastasis has been established. In the present study, the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) algorithm was applied to calculate the immune and stromal scores of patients with osteosarcoma based on data from The Cancer Genome Atlas database. A metagene approach and deconvolution method were used to reveal distinct TME landscapes in patients with osteosarcoma. Bioinformatics analysis was used to identify differentially expressed genes (DEGs) associated with metastasis and immune infiltration in osteosarcoma, and a risk model was constructed using the DEGs with potential prognostic significance. Subsequently, gene set enrichment and Spearman's correlation analyses were used to delineate the biological processes associated with these prognostic biomarkers. Finally, immunohistochemical (IHC) analysis was performed to evaluate the expression levels of immune infiltrates and prognostic biomarkers in clinical osteosarcoma tissues. The results of the ESTIMATE demonstrated that patients with non-metastatic osteosarcoma presented with higher immune/stromal scores and a more favorable prognosis compared with those with metastatic osteosarcoma. The TME landscapes in patients with osteosarcoma suggested that high levels of tumor-infiltrating immune cells (TIICs) may suppress metastasis. Increased numbers of CD56bright natural killer cells, immature B cells, M1 macrophages and neutrophils, and lower levels of M2 macrophages were observed in the non-metastatic tissues compared with those in the metastatic tissues. A total of 69 DEGs were identified to be associated with metastasis and immune infiltration in osteosarcoma. Of these, GATA3, LPAR5, EVI2B, RIAM and CFH exhibited prognostic potential and were highly expressed in non-metastatic osteosarcoma tissues based on the IHC analysis results. These biomarkers were involved in various immune-related biological processes and were positively associated with multiple TIICs and immune signatures. The risk model constructed using these prognostic biomarkers demonstrated high predictive accuracy for the prognosis of osteosarcoma. In conclusion, the present study proposed a five-biomarker prognostic signature for the prediction of metastasis and immune infiltration in patients with osteosarcoma.

Project description:BackgroundPyroptosis is closely related to cancer prognosis. In this study, we tried to construct an individualized prognostic risk model for hepatocellular carcinoma (HCC) based on within-sample relative expression orderings (REOs) of pyroptosis-related lncRNAs (PRlncRNAs).MethodsRNA-seq data of 343 HCC samples derived from The Cancer Genome Atlas (TCGA) database were analyzed. PRlncRNAs were detected based on differentially expressed lncRNAs between sample groups clustered by 40 reported pyroptosis-related genes (PRGs). Univariate Cox regression was used to screen out prognosis-related PRlncRNA pairs. Then, based on REOs of prognosis-related PRlncRNA pairs, a risk model for HCC was constructed by combining LASSO and stepwise multivariate Cox regression analysis. Finally, a prognosis-related competing endogenous RNA (ceRNA) network was built based on information about lncRNA-miRNA-mRNA interactions derived from the miRNet and TargetScan databases.ResultsHierarchical clustering of HCC patients according to the 40 PRGs identified two groups with a significant survival difference (Kaplan-Meier log-rank, p = 0.026). Between the two groups, 104 differentially expressed lncRNAs were identified (|log2(FC)|> 1 and FDR < 5%). Among them, 83 PRlncRNA pairs showed significant associations between their REOs within HCC samples and overall survival (Univariate Cox regression, p < 0.005). An optimal 11-PRlncRNA-pair prognostic risk model was constructed for HCC. The areas under the curves (AUCs) of time-dependent receiver operating characteristic (ROC) curves of the risk model for 1-, 3-, and 5-year survival were 0.737, 0.705, and 0.797 in the validation set, respectively. Gene Set Enrichment Analysis showed that inflammation-related interleukin signaling pathways were upregulated in the predicted high-risk group (p < 0.05). Tumor immune infiltration analysis revealed a higher abundance of regulatory T cells (Tregs) and M2 macrophages and a lower abundance of CD8 + T cells in the high-risk group, indicating that excessive pyroptosis might occur in high-risk patients. Finally, eleven lncRNA-miRNA-mRNA regulatory axes associated with pyroptosis were established.ConclusionOur risk model allowed us to determine the robustness of the REO-based PRlncRNA prognostic biomarkers in the stratification of HCC patients at high and low risk. The model is also helpful for understanding the molecular mechanisms between pyroptosis and HCC prognosis. High-risk patients may have excessive pyroptosis and thus be less sensitive to immune therapy.

Project description:BackgroundColonic adenocarcinoma (COAD) is a common gastrointestinal tract tumor, and its occurrence and progression are typically associated with genomic instability, tumor-suppressor gene and oncogene mutations, and tumor mutational load. N6-methyladenosine (m6A) modification of RNAs and long non-coding RNA (lncRNA) expression are important in tumorigenesis and progression. However, the regulatory roles of m6A-associated lncRNAs in the tumor microenvironment, stratification of prognosis, and immunotherapy are unclear.MethodsWe screened 43 prognostic lncRNAs linked to m6A and performed consistent molecular typing of COAD using consensus clustering. The single-sample Gene Set Enrichment Analysis and ESTIMATE algorithms were used to assess the immune characteristics of different subgroups. Covariation between methylation-related prognostic lncRNAs was eliminated by least absolute shrinkage and selection operator Cox regression. A nomogram was created and evaluated by combining the methylation-related prognostic lncRNA model with other clinical factors. The relationship between the prognostic model grouping and microsatellite instability, immunophenotype score, and tumor mutation burden was validated using R scripts. Finally, we used a linkage map to filter sensitive medicines to suppress the expression of high-risk genes. Three m6A-associated lncRNA modes were identified in 446 COAD specimens with different clinical endpoints and biological statuses. Risk scores were constructed based on the m6A-associated lncRNA signature genes. Patients with lower risk scores showed superior immunotherapy responses and clinical benefits compared to those with higher risk scores. Lower risk scores were also correlated with higher immunophenotype scores, tumor mutation burden, and mutation rates in significantly mutated genes (e.g., FAT4 and MUC16). Piperidolate, quinostatin, and mecamylamin were screened for their abilities to suppress the expression of high-risk genes in the model.ConclusionsQuantitative assessment of m6A-associated lncRNAs in single tumors can enhance the understanding of tumor microenvironment profiles. The prognostic model constructed using m6A-associated lncRNAs may facilitate prognosis and immunotherapy stratification of patients with COAD; finally, three drugs with potential therapeutic value were screened based on the model.

Dataset Information

Screening and verification of prognostic lncRNA markers related to immune infiltration in the metastasis of osteosarcoma.

Background

Methods

Results

Conclusions

Publications

Screening and verification of prognostic lncRNA markers related to immune infiltration in the metastasis of osteosarcoma.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets