Comprehensive analysis of immune related lncRNAs in the tumor microenvironment of stage II–III colorectal cancer
Ontology highlight
ABSTRACT: Background
Long non-coding RNAs (lncRNAs) associated with immunological function have increasingly been found to act as effective prognostic biomarkers of the overall survival (OS) of colorectal cancer (CRC) patients. We sought to identify a signature of immune-related lncRNAs that offered value as a tool for the prospective prognostic evaluation of patients with stage II–III CRC. Methods
The clinical and gene expression data of CRC patients in The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases was obtained and separated into a training cohort composed of 202 samples, a test cohort of 124 samples from the GSE72970 dataset, and a validation cohort of 91 samples from the GSE143985 dataset. Results
We firstly evaluated intratumoral immune cell infiltration by conducting a Single-sample gene set enrichment analyses (ssGSEA) analysis to separate patient tumors into those with low immune cell infiltration and those with high immune cell infiltration. We then compared lncRNA and mRNA expression profiles between these two tumor types, leading us to focus on eight lncRNAs identified within the resultant mRNA-lncRNA co-expression network. Multivariate Cox regression models were then utilized to detect an immune-associated lncRNA signature that offered value for prognostic model construction. Functional analyses revealed this lncRNA signature to be associated with key immunological pathways including the JAK-STAT signaling, T cell receptor signaling, and Rap1 signaling pathways. Conclusions
Together, our results suggest that our immune-related 4 lncRNA signature can reliably predict stage II–III CRC patient prognosis, thereby guiding efforts to better understand this disease and to effectively treat it.
SUBMITTER: Zhang Z
PROVIDER: S-EPMC8576234 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundN6-methyladenosine (m6A) modification and long non-coding RNAs (lncRNAs) play pivotal roles in gastric cancer (GC) progression. The emergence of immunotherapy in GC has created a paradigm shift in the approaches of treatment, whereas there is significant heterogeneity with regard to degree of treatment responses, which results from the variability of tumor immune microenvironment (TIME). How the interplay between m6A and lncRNAs enrolling in the shaping of TIME remains unclear.MethodsThe RNA sequencing and clinical data of GC patients were collected from TCGA database. Pearson correlation test and univariate Cox analysis were used to screen out m6A-related lncRNAs. Consensus clustering method was implemented to classify GC patients into two clusters. Survival analysis, the infiltration level of immune cells, Gene set enrichment analysis (GSEA) and the mutation profiles were analyzed and compared between two clusters. A competing endogenous RNA (ceRNA) network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied for the identification of pathways in which m6A-related lncRNAs enriched. Then least absolute shrinkage and selection operator (LASSO) COX regression was implemented to select pivotal lncRNAs, and risk model was constructed accordingly. The prognosis value of the risk model was explored. In addition, the response to immune checkpoint inhibitors (ICIs) therapy were compared between different risk groups. Finally, we performed qRT-PCR to detect expression patterns of the selected lncRNAs in the 35 tumor tissues and their paired adjacent normal tissues, and validated the prognostic value of risk model in our cohort (N = 35).ResultsThe expression profiles of 15 lncRNAs were included to cluster patients into 2 subtypes. Cluster1 with worse prognosis harbored higher immune score, stromal score, ESTIMATE score and lower mutation rates of the genes. Different immune cell infiltration patterns were also displayed between the two clusters. GSEA showed that cluster1 preferentially enriched in tumor hallmarks and tumor-related biological pathways. KEGG pathway analysis found that the target mRNAs which m6A-related lncRNAs regulated by sponging miRNAs mainly enriched in vascular smooth muscle contraction, cAMP signaling pathway and cGMP-PKG signaling pathway. Next, eight lncRNAs were selected by LASSO regression algorithm to construct risk model. Patients in the high-risk group had poor prognoses, which were consistent in our cohort. As for predicting responses to ICIs therapy, patients from high-risk group were found to have lower tumor mutation burden (TMB) scores and account for large proportion in the Microsatellite Instability-Low (MSI-L) subtype. Moreover, patients had distinct immunophenoscores in different risk groups.ConclusionOur study revealed that the interplay between m6A modification and lncRNAs might have critical role in predicting GC prognosis, sculpting TIME landscape and predicting the responses to ICIs therapy.
| S-EPMC8943990 | biostudies-literature
Project description:Cellular subpopulations within the colorectal tumor microenvironment (TME) include CD3+ and CD8+ lymphocytes, CD68+ and CD163+ macrophages, and tumor buds (TBs), all of which have known prognostic significance in stage II colorectal cancer. However, the prognostic relevance of their spatial interactions remains unknown. Here, by applying automated image analysis and machine learning approaches, we evaluate the prognostic significance of these cellular subpopulations and their spatial interactions. Resultant data, from a training cohort retrospectively collated from Edinburgh, UK hospitals (n?=?113), were used to create a combinatorial prognostic model, which identified a subpopulation of patients who exhibit 100% survival over a 5-year follow-up period. The combinatorial model integrated lymphocytic infiltration, the number of lymphocytes within 50-?m proximity to TBs, and the CD68+/CD163+ macrophage ratio. This finding was confirmed on an independent validation cohort, which included patients treated in Japan and Scotland (n?=?117). This work shows that by analyzing multiple cellular subpopulations from the complex TME, it is possible to identify patients for whom surgical resection alone may be curative.
| S-EPMC7229187 | biostudies-literature
Project description:Background: Globally, pancreatic adenocarcinoma (PAAD) is a common and highly devastating gastrointestinal malignancy that seriously threatens human health. Pyroptosis refers to an emerging form of programmed cell death that has been discovered in recent years, and studies have demonstrated that long non-coding RNA (lncRNA) may act as a moderator in the pyroptosis process of cancer cells. However, relevant explorations about lncRNAs and pyroptosis are still insufficient in PAAD. Therefore, our research is designed to make a comprehensive analysis of the potential values of pyroptosis-related lncRNAs in PAAD. Methods: By integrating the RNA-sequencing, somatic mutation, and copy number variation (CNV) datasets, as well as the clinicopathological features, we established and validated a risk signature based on pyroptosis-related lncRNAs, and comprehensively analyzed its clinical significance and the potential connection with the tumor immune microenvironment (TIME). Consequences: The genetic variation landscape displayed that the somatic mutations were rare while CNV changes were general and mainly concentrated on copy number amplification of these 52 pyroptosis-related genes. Subsequently, a risk signature consisting of 10 lncRNAs (TRAF3IP2-AS1, LINC00519, LINC01133, LINC02251, AC005332.6, AL590787.1, AC090114.2, TRPC7-AS1, MIR223HG, and MIR3142HG) was constructed and patients were divided into different subgroups according to the median risk score; patients with high-risk scores presented worse outcomes compared to those with low-risk scores in the training, testing, and entire cohorts. Furthermore, patients at low-risk scores possessed a higher infiltration abundance of immune cells compared with high-risk patients, which was consistent with the expression levels of lncRNAs between the high/low-risk groups. Drug sensitivity analysis showed that low-risk scores were related to anti-cancer agents like AICAR and Axitinib, whereas high-risk scores were connected with certain drugs such as AUY922. These results demonstrated that our risk signature could be used for prognosis prediction; additionally, it was also related to the TIME that might act as a potential indicator to instruct immunotherapeutic strategies. Conclusion: This work explored the significance of the risk model constructed by pyroptosis-related lncRNAs in prognosis prediction and its internal link with the immune microenvironment of PAAD. The results are expected to assist in the diagnosis, prognostic assessment, and management of patients with PAAD.
| S-EPMC9296806 | biostudies-literature
Project description:Background: Accumulating evidence shows that pyroptosis plays a crucial role in hepatocellular carcinoma (HCC). However, the relationship between pyroptosis-related long non-coding RNAs (lncRNAs) and HCC tumor characteristics remains enigmatic. We aimed to explore the predictive effect of pyroptosis-related lncRNAs (PRLs) in the prognosis of HCC. Methods: We comprehensively analyzed the role of the PRLs in the tumor microenvironment and HCC prognosis by integrating genomic data from patients of HCC. Consensus clustering analysis of PRLs was applied to identify HCC subtypes. A prognostic model was then established with a training cohort from The Cancer Genome Atlas (TCGA) using univariate and least absolute shrinkage and selection operator (LASSO) Cox regression analysis. Further, we evaluated the accuracy of this predictive model using a validation set. We predicted IC50s of commonly used chemotherapeutic and targeted drugs through the R package pRRophetic. Results: Based on pyroptosis-related lncRNAs, a prognostic risk signature composed of seven PRLs (MKLN1AS, AL031985.3, SNHG4, GHRLOS, AC005479.2, AC099850.4, and AC026412.3) was established. For long-term prognosis of HCC patients, our model shows excellent accuracy to forecast overall survival of HCC individuals both in training set and testing set. We found a significant correlation between clinical features and the risk score. Patients in the high-risk group had tumor characteristics associated with progression such as aggressive pathological grade and stage. Besides that, gene set enrichment analysis (GSEA) showed that cell cycle and focal adhesion were significantly enriched in the high-risk group. Conclusion: The association of the risk model constituted by these seven pyroptosis-related lncRNAs with clinical prognosis, tumor microenvironment, chemotherapy and small molecule drugs was evaluated. Our study provides strong evidence for individualized prediction of prognosis, shedding light on immunotherapy in HCC patients.
| S-EPMC9087742 | biostudies-literature
Project description:Background: Abnormal expression of lncRNA is closely related to the occurrence and metastasis of osteosarcoma. The tumor immune microenvironment (TIM) is considered to be an important factor affecting the prognosis and treatment of osteosarcoma. This study aims to explore the effect of immune-related lncRNAs (IRLs) on the prognosis of osteosarcoma and its relationship with the TIM. Methods: Ninety-five osteosarcoma samples from the TARGET database were included. Iterative LASSO regression and multivariate Cox regression analysis were used to screen the IRLs signature with the optimal AUC. The predict function was used to calculate the risk score and divide osteosarcoma into a high-risk group and low-risk group based on the optimal cut-off value of the risk score. The lncRNAs in IRLs signature that affect metastasis were screened for in vitro validation. Single sample gene set enrichment analysis (ssGSEA) and ESTIMATE algorithms were used to evaluate the role of TIM in the influence of IRLs on osteosarcoma prognosis. Results: Ten IRLs constituted the IRLs signature, with an AUC of 0.96. The recurrence and metastasis rates of osteosarcoma in the high-risk group were higher than those in the low-risk group. In vitro experiments showed that knockdown of lncRNA (AC006033.2) could increase the proliferation, migration, and invasion of osteosarcoma. ssGSEA and ESTIMATE results showed that the immune cell content and immune score in the low-risk group were generally higher than those in the high-risk group. In addition, the expression levels of immune escape-related genes were higher in the high-risk group. Conclusion: The IRLs signature is a reliable biomarker for the prognosis of osteosarcoma, and they alter the prognosis of osteosarcoma. In addition, IRLs signature and patient prognosis may be related to TIM in osteosarcoma. The higher the content of immune cells in the TIM of osteosarcoma, the lower the risk score of patients and the better the prognosis. The higher the expression of immune escape-related genes, the lower the risk score of patients and the better the prognosis.
| S-EPMC8529014 | biostudies-literature
Project description:Background: Hepatocellular carcinoma (HCC) remains the most prevalent gastrointestinal malignancy worldwide, with robust drug resistance to therapy. N7-methylguanosine (m7G) mRNA modification has been significantly related to massive human diseases. Considering the effect of m7G-modified long non-coding RNAs (lncRNAs) in HCC progression is unknown, the study aims at investigating a prognostic signature to improve clinical outcomes for patients with HCC. Methods: Two independent databases (TCGA and ICGC) were used to analyze RNAseq data of HCC patients. First, co-expression analysis was applied to obtain the m7G-related lncRNAs. Moreover, consensus clustering analysis was employed to divide HCC patients into clusters. Then, using least absolute shrinkage and selection operator-Cox regression analysis, the m7G-related lncRNA prognostic signature (m7G-LPS) was first tested in the training set and then confirmed in both the testing and ICGC sets. The expression levels of the nine lncRNAs were further confirmed via real-time PCR in cell lines, principal component analysis, and receiver operating characteristic curve. The m7G-LPS could divide HCC patients into two different risk groups with the optimal risk score. Then, Kaplan-Meier curves, tumor mutation burden (TMB), therapeutic effects of chemotherapy agents, and expressions of immune checkpoints were performed to further enhance the availability of immunotherapeutic treatments for HCC patients. Results: A total of 1465 lncRNAs associated with the m7G genes were finally selected from the TCGA database, and through the univariate Cox regression, the expression levels of 22 m7G-related lncRNAs were concerning HCC patients' overall survival (OS). Then, the whole patients were grouped into two subgroups, and the OS in Cluster 1 was longer than that of patients in Cluster 2. Furthermore, nine prognostic m7G-related lncRNAs were identified to conduct the m7G-LPS, which were further verified. A prognostic nomogram combined age, gender, HCC grade, stage, and m7G-LPS showed strong reliability and accuracy in predicting OS in HCC patients. Finally, immune checkpoint expression, TMB, and several chemotherapy agents were remarkably associated with risk scores. More importantly, the OS of the TMB-high patients was the worst among the four groups. Conclusion: The prognostic model we established was validated by abundant algorithms, which provided a new perspective on HCC tumorigenesis and thus improved individualized treatments for patients.
| S-EPMC9445240 | biostudies-literature
Project description:Distant metastasis is the major causes of death in colorectal cancer (CRC) patients. In order to identify genes influencing the prognosis of patients with CRC, we compared gene expression in primary tumors with and without distant metastasis using an oligonucleotide microarray. We also examined the expression of the candidate gene in 100 CRC patients by quantitative real-time reverse transcription PCR and studied the relationship between its expression and the prognosis of patients with CRC. As a result, we identified MUC12 as a candidate gene involved in metastasis processes by microarray analysis. Quantitative real-time reverse transcription PCR showed that MUC12 expression was significantly lower in cancer tissues than in adjacent normal tissues (P < 0.001). In stage II and stage III CRC, patients with low expression showed worse disease-free survival (P = 0.038). Multivariate analysis disclosed that MUC12 expression status was an independent prognostic factor in stage II and stage III CRC (relative risk, 9.532; 95% confidence interval, 2.303-41.905; P = 0.002). This study revealed the prognostic value of MUC12 expression in CRC patients. Moreover, our result suggests MUC12 expression is a possible candidate gene for assessing postoperative adjuvant therapy for CRC patients. Total of 111 microarray datasets (77 for LCM samples, and 17 pairs for homogenized samples from tumor and adjacent tissues) were normalized using robust multi-array average (RMA) method under R 2.6.2 statistical software together with BioConductor package, as described previously. Then, the gene expression levels were log2-transformed, and 62 control probe sets were removed for further analysis. In order to identify a set of genes associated with development of metastatic recurrence, we performed Wilcoxon rank-sum test for gene expression differences of 54,613 probe sets between recurrence and non-recurrence groups. Similarly, Wilcoxon singed-rank test was conducted to select genes which showed significant expression difference between tumor and adjacent tissue. Then, we selected a set of genes that satisfied both of above two criteria.
2010-03-05 | E-GEOD-18105 | biostudies-arrayexpress
Project description:Distant metastasis is the major causes of death in colorectal cancer (CRC) patients. In order to identify genes influencing the prognosis of patients with CRC, we compared gene expression in primary tumors with and without distant metastasis using an oligonucleotide microarray. We also examined the expression of the candidate gene in 100 CRC patients by quantitative real-time reverse transcription PCR and studied the relationship between its expression and the prognosis of patients with CRC. As a result, we identified MUC12 as a candidate gene involved in metastasis processes by microarray analysis. Quantitative real-time reverse transcription PCR showed that MUC12 expression was significantly lower in cancer tissues than in adjacent normal tissues (P < 0.001). In stage II and stage III CRC, patients with low expression showed worse disease-free survival (P = 0.038). Multivariate analysis disclosed that MUC12 expression status was an independent prognostic factor in stage II and stage III CRC (relative risk, 9.532; 95% confidence interval, 2.303-41.905; P = 0.002). This study revealed the prognostic value of MUC12 expression in CRC patients. Moreover, our result suggests MUC12 expression is a possible candidate gene for assessing postoperative adjuvant therapy for CRC patients.
2010-02-04 | GSE18105 | GEO
Project description:With the increasing incidence of colorectal cancer (CRC) and continued difficulty in treating it using immunotherapy, there is an urgent need to identify an effective immune-related biomarker associated with the survival and prognosis of patients with this disease. DNA methylation plays an essential role in maintaining cellular function, and changes in methylation patterns may contribute to the development of autoimmunity, aging, and cancer. In this study, we aimed to identify a novel immune-related methylated signature to aid in predicting the prognosis of patients with CRC. We investigated DNA methylation patterns in patients with stage II/III CRC using datasets from The cancer genome atlas (TCGA). Overall, 182 patients were randomly divided into training (n = 127) and test groups (n = 55). In the training group, five immune-related methylated CG sites (cg11621464, cg13565656, cg18976437, cg20505223, and cg20528583) were identified, and CG site-based risk scores were calculated using univariate Cox proportional hazards regression in patients with stage II/III CRC. Multivariate Cox regression analysis indicated that methylated signature was independent of other clinical parameters. The Kaplan-Meier analysis results showed that CG site-based risk scores could significantly help distinguish between high- and low-risk patients in both the training (P = 0.000296) and test groups (P = 0.022). The area under the receiver operating characteristic curve in the training and test groups were estimated to be 0.771 and 0.724, respectively, for prognosis prediction. Finally, stratified analysis results suggested the remarkable prognostic value of CG site-based risk scores in CRC subtypes. We identified five methylated CG sites that could be used as an efficient overall survival (OS)-related biomarker for stage II/III CRC patients.
| S-EPMC8273301 | biostudies-literature
Project description:Cuproptosis was characterized as a novel type of programmed cell death. Recently, however, the role of cuproptosis-related long noncoding RNAs (CRLs) in tumors has not yet been studied. Identifying a predictive CRL signature in hepatocellular carcinoma (HCC) and investigating its putative molecular function were the goals of this work. Initially, Pearson's test was used to assess the relationship between lncRNAs and cuproptosis-associated genes obtained from HCC data of The Cancer Genome Atlas (TCGA). By implementing differential expression and univariate Cox analysis, 61 prognostic CRLs were subsequent to the least absolute shrinkage and selection operator (LASSO) Cox regression analysis. A prognostic risk score model was then constructed to evaluate its ability to predict patients' survival when combined with clinicopathological parameters in HCC. The five-lncRNA prognostic signature categorized the HCC patients into high- and low-risk groups. The low-risk group exhibited more sensitivity to elesclomol than the high-risk one. Surprisingly, distinct mitochondrial metabolism pathways connected to cuproptosis and pivotal immune-related pathways were observed between the two groups via gene set enrichment analysis (GSEA). Meanwhile, there were substantial differences between the high-risk group and the low-risk group in terms of tumor-infiltrating immune cells (TIICs). Furthermore, a positive relationship was shown between the risk score and the expression of immune checkpoints. Additionally, differential expression of the five lncRNAs was confirmed in our own HCC samples and cell lines via RT-qPCR. Finally, in vitro assays confirmed that WARS2-AS1 and MKLN1-AS knockdown could sensitize HCC cells to elesclomol-induced cuproptosis. Overall, our predictive signature may predict the prognosis of HCC patients in an independent manner, give a better understanding of how CRLs work in HCC, and offer therapeutic reference for patients with HCC.
| S-EPMC9846072 | biostudies-literature