Construction of immune-related and prognostic lncRNA clusters and identification of their immune and genomic alterations characteristics in lung adenocarcinoma samples.
Ontology highlight
ABSTRACT:
SUBMITTER: Li J
PROVIDER: S-EPMC7288974 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundLung adenocarcinoma (LUAD) originates mainly from the mucous epithelium and glandular epithelium of the bronchi. It is the most common pathologic subtype of non-small cell lung cancer (NSCLC). At present, there is still a lack of clear criteria to predict the efficacy of immunotherapy. The 5-year survival rate for LUAD patients remains low.MethodsAll data were downloaded from The Cancer Genome Atlas (TCGA) database. We used Gene Set Enrichment Analysis (GSEA) database to obtain immune-related mRNAs. Immune-related lncRNAs were acquired by using the correlation test of the immune-related genes with R version 3.6.3 (Pearson correlation coefficient cor = 0.5, P < 0.05). The TCGA-LUAD dataset was divided into the testing set and the training set randomly. Based on the training set to perform univariate and multivariate Cox regression analyses, we screened prognostic immune-related lncRNAs and given a risk score to each sample. Samples were divided into the high-risk group and the low-risk group according to the median risk score. By the combination of Kaplan-Meier (KM) survival curve, the receiver operating characteristic (ROC) (AUC) curve, the independent risk factor analysis, and the clinical data of the samples, we assessed the accuracy of the risk model. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the differentially expressed mRNAs between the high-risk group and the low-risk group. The differentially expressed genes related to immune response between two risk groups were analyzed to evaluate the role of the model in predicting the efficacy and effects of immunotherapy. In order to explain the internal mechanism of the risk model in predicting the efficacy of immunotherapy, we analyzed the differentially expressed genes related to epithelial-mesenchymal transition (EMT) between two risk groups. We extracted RNA from normal bronchial epithelial cell and LUAD cells and verified the expression level of lncRNAs in the risk model by a quantitative real-time polymerase chain reaction (qRT-PCR) test. We compared our risk model with other published prognostic signatures with data from an independent cohort. We transfected LUAD cell with siRNA-LINC0253. Western blot analysis was performed to observed change of EMT-related marker in protein level.ResultsThrough univariate Cox regression analysis, 24 immune-related lncRNAs were found to be strongly associated with the survival of the TCGA-LUAD dataset. Utilizing multivariate Cox regression analysis, 10 lncRNAs were selected to establish the risk model. The K-M survival curves and the ROC (AUC) curves proved that the risk model has a fine predictive effect. The GO enrichment analysis indicated that the effect of the differentially expressed genes between high-risk and low-risk groups is mainly involved in immune response and intercellular interaction. The KEGG enrichment analysis indicated that the differentially expressed genes between high-risk and low-risk groups are mainly involved in endocytosis and the MAPK signaling pathway. The expression of genes related to the efficacy of immunotherapy was significantly different between the two groups. A qRT-PCR test verified the expression level of lncRNAs in LUAD cells in the risk model. The AUC of ROC of 5 years in the independent validation dataset showed that this model had superior accuracy. Western blot analysis verified the change of EMT-related marker in protein level.ConclusionThe immune lncRNA risk model established by us could better predict the prognosis of patients with LUAD.
| S-EPMC8044985 | biostudies-literature
Project description:BackgroundPancreatic ductal adenocarcinoma (PDAC) is a highly aggressive malignant tumor of the digestive system. Its grim prognosis is mainly attributed to the lack of means for early diagnosis and poor response to treatments. Genomic instability is shown to be an important cancer feature and prognostic factor, and its pattern and extent may be associated with poor treatment outcomes in PDAC. Recently, it has been reported that long non-coding RNAs (lncRNAs) play a key role in maintaining genomic instability. However, the identification and clinical significance of genomic instability-related lncRNAs in PDAC have not been fully elucidated.MethodsGenomic instability-derived lncRNA signature (GILncSig) was constructed based on the results of multiple regression analysis combined with genomic instability-associated lncRNAs and its predictive power was verified by the Kaplan-Meier method. And real-time quantitative polymerase chain reaction (qRT-PCR) was used for simple validation in human cancers and their adjacent non-cancerous tissues. In addition, the correlation between GILncSig and tumor microenvironment (TME) and epithelial-mesenchymal transition (EMT) was investigated by Pearson correlation analysis.ResultsThe computational framework identified 206 lncRNAs associated with genomic instability in PDAC and was subsequently used to construct a genome instability-derived five lncRNA-based gene signature. Afterwards, we successfully validated its prognostic capacity in The Cancer Genome Atlas (TCGA) cohort. In addition, via careful examination of the transcriptome expression profile of PDAC patients, we discovered that GILncSig is associated with EMT and an adaptive immunity deficient immune profile within TME.ConclusionsOur study established a genomic instability-associated lncRNAs-derived model (GILncSig) for prognosis prediction in patients with PDAC, and revealed the potential functional regulatory role of GILncSig.
| S-EPMC9486402 | biostudies-literature
Project description:Background: There is increasing evidence that long non-coding RNAs (lncRNAs) can be used as potential prognostic factors for cancer. This study aimed to develop a prognostic model for lung adenocarcinoma (LUAD) using angiogenesis-related long non-coding RNAs (lncRNAs) as potential prognostic factors. Methods: Transcriptome data from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) were analyzed to identify aberrantly expressed angiogenesis-related lncRNAs in LUAD. A prognostic signature was constructed using differential expression analysis, overlap analysis, Pearson correlation analysis, and Cox regression analysis. The model's validity was assessed using K-M and ROC curves, and independent external validation was performed in the GSE30219 dataset. Prognostic lncRNA-microRNA (miRNA)-messenger RNA (mRNA) competing endogenous RNA (ceRNA) networks were identified. Immune cell infiltration and mutational characteristics were also analyzed. The expression of four human angiogenesis-associated lncRNAs was quantified using quantitative real-time PCR (qRT-PCR) gene arrays. Results: A total of 26 aberrantly expressed angiogenesis-related lncRNAs in LUAD were identified, and a Cox risk model based on LINC00857, RBPMS-AS1, SYNPR-AS1, and LINC00460 was constructed, which may be an independent prognostic predictor for LUAD. The low-risk group had a significant better prognosis and was associated with a higher abundance of resting immune cells and a lower expression of immune checkpoint molecules. Moreover, 105 ceRNA mechanisms were predicted based on the four prognostic lncRNAs. qRT-PCR results showed that LINC00857, SYNPR-AS1, and LINC00460 were significantly highly expressed in tumor tissues, while RBPMS-AS1 was highly expressed in paracancerous tissues. Conclusion: The four angiogenesis-related lncRNAs identified in this study could serve as a promising prognostic biomarker for LUAD patients.
| S-EPMC10043447 | biostudies-literature
Project description:Necroptosis is a new form of cell death. Since the discovery that long non-coding RNAs can affect the proliferation of lung adenocarcinoma, much has been learned about it, yet those of necroptosis-related long non-coding RNAs (NRlncRNAs) in lung adenocarcinoma (LUAD) remain enigmatic. This study aims to explore novel biomarkers and therapeutic targets for LUAD. The LUAD data was downloaded from The Cancer Genome Atlas, and necroptosis-related genes were retrieved from published literature. Co-expression analysis, univariate Cox analysis, least absolute shrinkage and selection operator regression analysis were used to identify necroptosis-related prognostic long non-coding RNAs. A comprehensive evaluation of tumor immunity for necrosis-related features was performed, and we identified a 9-NRlncRNA signature. Kaplan-Meier and Cox regression analyses confirmed that the signature was an independent predictor of LUAD outcome in the test and train sets (all P < 0.05). The areas of 1-, 2-, and 3-year overall survival under the time-dependent receiver operating characteristics (ROC) curve (AUC) were 0.754, 0.746, and 0.720, respectively. The GSEA results showed that 9 NRlncRNAs were associated with multiple malignancy-associated and immunoregulatory pathways. Based on this model, we found that the immune status and level of response to chemotherapy and targeted therapy were significantly different in the low-risk group compared with the high-risk group. qRT-PCR assay revealed that 9 NRlncRNAs were involved in the regulation of tumor cell proliferation and may affect the expression of programmed cell death 1 (PD1) and CD28 at human immune checkpoints. Our results indicated that the novel signature involving 9 NRlncRNAs (AL031600.2, LINC01281, AP001178.1, AL157823.2, LINC01290, MED4-AS1, AC026355.2, AL606489.1, FAM83A-AS1) can predict the prognosis of LUAD and are associated with the immune response. This will provide new insights into the pathogenesis and development of therapies for LUAD.
| S-EPMC10560937 | biostudies-literature
Project description:Non-small-cell lung cancer (NSCLC) has been recognized as a highly heterogeneous disease with phenotypic and genotypic diversity in each subgroup. While never-smoker patients with NSCLC have been well studied through next generation sequencing, we have yet to recognize the potentially unique molecular features of young never-smoker patients with NSCLC. In this study, we conducted whole genome sequencing (WGS) to characterize the genomic alterations of 36 never-smoker Chinese patients, who were diagnosed with lung adenocarcinoma (LUAD) at 45 years or younger. Besides the well-known gene mutations (e.g., TP53 and EGFR), our study identified several potential lung cancer-associated gene mutations that were rarely reported (e.g., HOXA4 and MST1). The lung cancer-related copy number variations (e.g., EGFR and CDKN2A) were enriched in our cohort (41.7%, 15/36) and the lung cancer-related structural variations (e.g., EML4-ALK and KIF5B-RET) were commonly observed (22.2%, 8/36). Notably, new fusion partners of ALK (SMG6-ALK) and RET (JMJD1C-RET) were found. Furthermore, we observed a high prevalence (63.9%, 23/36) of potentially targetable genomic alterations in our cohort. Finally, we identified germline mutations in BPIFB1 (rs6141383, p.V284M), CHD4 (rs74790047, p.D140E), PARP1 (rs3219145, p.K940R), NUDT1 (rs4866, p.V83M), RAD52 (rs4987207, p.S346*), and MFI2 (rs17129219, p.A559T) were significantly enriched in the young never-smoker patients with LUAD when compared with the in-house noncancer database (p < 0.05). Our study provides a detailed mutational portrait of LUAD occurring in young never-smokers and gives insights into the molecular pathogenesis of this distinct subgroup of NSCLC.
| S-EPMC6175072 | biostudies-literature
Project description:BackgroundGenomic instability (GI) is a critical feature of cancer which plays a key role in the occurrence and development of pancreatic adenocarcinoma (PAAD). Long non-coding RNA (LncRNA) is an emerging prognostic biomarker because it is involved in regulating GI. Recently, researchers used such GI-related LncRNAs (GILncRNAs) to establish a prognostic signature for patients with cancer and helped in predicting the overall prognosis of the patients. However, it is evident that patients with PAAD still lack such prognostic signature constructed with GILncRNA.MethodsThe present study screened GILncRNAs from 83 patients with PAAD. Prognosis-related GILncRNAs were identified by univariate Cox regression analysis. The correlation coefficients of these GILncRNAs were obtained by multivariate Cox regression analysis and used to construct a signature. The signature in the present study was then assessed through survival analysis, mutation correlation analysis, independent prognostic analysis, and clinical stratification analysis in the training set and validated in the testing as well as all TCGA set. The current study performed external clinical relevance validation of the signature and validated the effect of AC108134.2 in GILncSig on PAAD using in vitro experiments. Finally, the function of GILncRNA signature (GILncSig) dependent on Gene Ontology enrichment analysis was explored and chemotherapeutic drug sensitivity analysis was also performed.ResultsResults of the present study found that a total of 409 GILncRNAs were identified, 5 of which constituted the prognostic risk signature in this study, namely, AC095057.3, AC108134.2, AC124798.1, AL606834.1, and AC104695.4. It was found that the signature of the present study was better than others in predicting the overall survival and applied to patients with PAAD of all ages, genders, and tumor grades. Further, it was noted that the signature of the current study in the GSE102238, was correlated with tumor length, and tumor stage of patients with PAAD. In vitro, functional experiments were used in the present study to validate that AC108134.2 is associated with PAAD genomic instability and progression. Notably, results of the pRRophetic analysis in the current study showed that the high-risk group possessed reverse characteristics and was sensitive to chemotherapy.ConclusionsIn conclusion, it was evident that the GILncSig used in the present study has good prognostic performance. Therefore, the signature may become a potential sensitive biological indicator of PAAD chemotherapy, which may help in clinical decision-making and management of patients with cancer.
| S-EPMC9012103 | biostudies-literature
Project description:BackgroundSuppressive tumor microenvironment is closely related to the progression and poor prognosis of lung adenocarcinoma (LUAD). Novel individual and universal immune-related biomarkers to predict the prognosis and immune landscape of LUAD patients are urgently needed. Two-gene pairing patterns could integrate and utilize various gene expression data.MethodsThe RNA-seq and relevant clinicopathological data of the LUAD project from the TCGA and well-known immune-related genes list from the ImmPort database were obtained. Co-expression analysis followed by an analysis of variance was performed to identify differentially expressed immune-related lncRNA (irlncRNA) (DEirlncRNA) between tumor and normal tissues. Two arbitrary DEirlncRNAs (DEirlncRNAs pair) in a tumor sample underwent pairwise comparison to generate a score (0 or 1). Next, Univariate analysis, Lasso regression and Multivariate analysis were used to screen survival-related DEirlncRNAs pairs and construct a prognostic model. The Acak information standard (AIC) values of the receiver operating characteristic (ROC) curve for 3 years are calculated to determine the cut-off point for high- or low-risk score. Finally, we evaluated the relationship between the risk score and overall survival, clinicopathological features, immune landscape, and chemotherapy efficacy.ResultsData of 54 normal and 497 tumor samples of LUAD were enrolled. After a strict screening process, 15 survival-independent-related DEirlncRNA pairs were integrated to construct a prognostic model. The AUC value of the 3-year ROC curve was 0.828. Kaplan-Meier analysis showed that patients with low risk lived longer than patients with high risk (p <0.001). Univariate and Multivariate Cox analysis suggested that the risk score was an independent factor of survival. The risk score was negatively associated with most tumor-infiltrating immune cells, immune score, and microenvironment scores. The low-risk group was correlated with increased expression of ICOS. The high-risk group had a connection with lower half inhibitory centration (IC50) of most chemotherapy drugs (e.g., etoposide, paclitaxel, vinorelbine, gemcitabine, and docetaxel) and targeted medicine-erlotinib, but with higher IC50 of methotrexate.ConclusionThe established irlncRNA pairs-based model is a promising prognostic signature for LUAD patients. Furthermore, the prognostic signature has great potential in the evaluation of tumor immune landscape and guiding individualized treatment regimens.
| S-EPMC8784752 | biostudies-literature
Project description:Bladder cancer (BLCA) is one of the most common urological cancer with increasing cases and deaths every year. In the present study, we aim to construct an immune-related prognostic lncRNA signature (IRPLS) in bladder cancer (BLCA) patients and explore its immunogenomic implications in pan-cancers. First, the immune-related differentially expressed lncRNAs (IRDELs) were identified by 'limma' R package and the score of IRPLS in every patient were evaluated by Cox regression. The dysregulation of IRDELs expression between cancer and para-cancer normal tissues was validated through RT-qPCR. Then, we further explore the biological functions of a novel lncRNA from IRPLS, RP11-89 in BLCA using CCK8 assay, Transwell assay and Apoptosis analysis, which indicated that RP11-89 was able to promote cell proliferation and invasive capacity while inhibits cell apoptosis in BLCA. In addition, we performed bioinformatic methods and RIP to investigate and validate the RP11-89/miR-27a-3p/PPARγ pathway in order to explore the mechanism. Next, CIBERSORT and ESTIMATE algorithm were used to evaluate abundance of tumour-infiltrating immune cells and scores of tumour environment elements in BLCA with different level of IRPLS risk scores. Finally, multiple bioinformatic methods were performed to show us the immune landscape of these four lncRNAs for pan-cancers. In conclusion, this study first constructed an immune-related prognostic lncRNA signature, which consists of RP11-89, PSORS1C3, LINC02672 and MIR100HG and might shed lights on novel targets for individualized immunotherapy for BLCA patients.
| S-EPMC8093983 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) is the most common subtype of lung cancer, but the prognosis of LUAD patients remains unsatisfactory. Here, we retrieved the RNA-seq data of LUAD cohort from The Cancer Genome Atlas (TCGA) database and then identified differentially expressed immune-related lncRNAs (DEirlncRNAs) between LUAD and normal controls. Based on a new method of cyclically single pairing along with a 0-or-1 matrix, we constructed a novel prognostic signature of 8 DEirlncRNA pairs in LUAD with no dependence upon specific expression levels of lncRNAs. This prognostic model exhibited significant power in distinguishing good or poor prognosis of LUAD patients and the values of the area under the curve (AUC) were all over 0.70 in 1, 3, 5 years receiver operating characteristic (ROC) curves. Moreover, the risk score of the model could serve as an independent prognostic factor for patients with LUAD. In addition, the risk model was significantly associated with clinicopathological characteristics, tumor-infiltrating immune cells, immune-related molecules and sensitivity of anti-tumor drugs. This novel signature of DEirlncRNA pairs in LUAD, which did not require specific expression levels of lncRNAs, might be used to guide the administration of patients with LUAD in clinical practice.
| S-EPMC8373953 | biostudies-literature
Project description:This study aimed to establish a prognostic risk model for lung adenocarcinoma (LUAD). We firstly divided 535 LUAD samples in TCGA-LUAD into high-, medium-, and low-immune infiltration groups by consensus clustering analysis according to immunological competence assessment by single-sample gene set enrichment analysis (ssGSEA). Profile of long non-coding RNAs (lncRNAs) in normal samples and LUAD samples in TCGA was used for a differential expression analysis in the high- and low-immune infiltration groups. A total of 1,570 immune-related differential lncRNAs in LUAD were obtained by intersecting the above results. Afterward, univariate COX regression analysis and multivariate stepwise COX regression analysis were conducted to screen prognosis-related lncRNAs, and an eight-immune-related-lncRNA prognostic signature was finally acquired (AL365181.2, AC012213.4, DRAIC, MRGPRG-AS1, AP002478.1, AC092168.2, FAM30A, and LINC02412). Kaplan-Meier analysis and ROC analysis indicated that the eight-lncRNA-based model was accurate to predict the prognosis of LUAD patients. Simultaneously, univariate COX regression analysis and multivariate COX regression analysis were undertaken on clinical features and risk scores. It was illustrated that the risk score was a prognostic factor independent from clinical features. Moreover, immune data of LUAD in the TIMER database were analyzed. The eight-immune-related-lncRNA prognostic signature was related to the infiltration of B cells, CD4+ T cells, and dendritic cells. GSEA enrichment analysis revealed significant differences in high- and low-risk groups in pathways like pentose phosphate pathway, ubiquitin mediated proteolysis, and P53 signaling pathway. This study helps to treat LUAD patients and explore molecules related to LUAD immune infiltration to deeply understand the specific mechanism.
| S-EPMC8301374 | biostudies-literature