Development and Validation of a Prognostic N6-Methyladenosine-Related Immune Gene Signature for Lung Adenocarcinoma.
Ontology highlight
ABSTRACT: Purpose
The prognostic value of an N6-methyladenosine (m6A) methylation-related immune gene signature for lung adenocarcinoma (LUAD) was investigated.Patients and methods
Gene expression and clinical phenotype data of LUAD patients were downloaded from The Cancer Genome Atlas database. A list of immune-related genes was retrieved from the InnateDB database. Correlation analysis, survival analysis, and univariate and multivariate Cox regression analyses were performed. After allocating patients into a high-risk or a low-risk group, the corresponding survival rates, immune microenvironment, expression of immune checkpoint genes, and modulation of Kyoto Encyclopedia of Genes and Genomes pathways were examined. Finally, the expression levels of prognostic biomarkers were assessed in the GSE126044 dataset.Results
Seven m6A-related immune prognostic genes were identified. High expression of PSMD10P1, DIDO1, ABCA5, and CIITA was associated with high survival rates, while that of PRC1, ZWILCH, and ANLN was associated with low survival rates. The high- and low-risk groups showed significant differences in terms of the abundance of six tumor-infiltrating immune cell types and expression of 12 immune checkpoint genes. The risk group acted as an independent prognostic factor (hazard ratio = 0.398, 95% confidence interval = 0.217-0.729, P = 0.003). Finally, the developed nomogram could predict most efficiently the 1-, 2-, and 3-year survival probability of LUAD patients with a C-index of 0.833.Conclusion
A seven-gene risk signature, associated with the immune microenvironment in LUAD, showed independent prognostic value.
SUBMITTER: Zhang L
PROVIDER: S-EPMC8643173 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundLung adenocarcinomas (LUAD) is the most common histological subtype of lung cancers. Tumor immune microenvironment (TIME) is involved in tumorigeneses, progressions, and metastases. This study is aimed to develop a robust immune-related signature of LUAD.MethodsA total of 1774 LUAD cases sourced from public databases were included in this study. Immune scores were calculated through ESTIMATE algorithm and weighted gene co-expression network analysis (WGCNA) was applied to identify immune-related genes. Stability selections and Lasso COX regressions were implemented to construct prognostic signatures. Validations and comparisons with other immune-related signatures were conducted in independent Gene Expression Omnibus (GEO) cohorts. Abundant infiltrated immune cells and pathway enrichment analyses were carried out, respectively, through ImmuCellAI and gene set enrichment analysis (GSEA).ResultsIn Cancer Genome Atlas (TCGA) LUAD cohorts, immune scores of higher levels were significantly associated with better prognoses (P < .05). Yellow (n = 270) and Blue (n = 764) colored genes were selected as immune-related genes, and after univariate Cox regression analysis (P < .005), a total of 133 genes were screened out for subsequent model constructions. A four-gene signature (ARNTL2, ECT2, PPIA, and TUBA4A) named IPSLUAD was developed through stability selection and Lasso COX regression. It was suggested by multivariate and subgroup analyses that IPSLUAD was an independent prognostic factor. It was suggested by Kaplan-Meier survival analysis that eight out of nine patients in high-risk groups had significantly worse prognoses in validation data sets (P < .05). IPSLUAD outperformed other signatures in two independent cohorts.ConclusionsA robust immune-related prognostic signature with great performances in multiple LUAD cohorts was developed in this study.
| S-EPMC7433810 | biostudies-literature
Project description:Lung cancer is the most commonly diagnosed cancer and the leading cause of cancer-related death. Among these, lung adenocarcinoma (LUAD) accounts for most cases. Due to the improvement of precision medicine based on molecular characterization, the treatment of LUAD underwent significant changes. With these changes, the prognosis of LUAD becomes diverse. N6-methyladenosine (m6A) is the most predominant modification in mRNAs, which has been a research hotspot in the field of oncology. Nevertheless, little has been studied to reveal the correlations between the m6A-related genes and prognosis in LUAD. Thus, we conducted a comprehensive analysis of m6A-related gene expressions in LUAD patients based on The Cancer Genome Atlas (TCGA) database by revealing their relationship with prognosis. Different expressions of the m6A-related genes in tumor tissues and non-tumor tissues were confirmed. Furthermore, their relationship with prognosis was studied via Consensus Clustering Analysis, Principal Components Analysis (PCA), and Least Absolute Shrinkage and Selection Operator (LASSO) Regression. Based on the above analyses, a m6A-based signature to predict the overall survival (OS) in LUAD was successfully established. Among the 479 cases, we found that most of the m6A-related genes were differentially expressed between tumor and non-tumor tissues. Six genes, HNRNPC, METTL3, YTHDC2, KIAA1429, ALKBH5, and YTHDF1 were screened to build a risk scoring signature, which is strongly related to the clinical features pathological stages (p < 0.05), M stages (p < 0.05), T stages (p?<?0.05), gender (p = 0.04), and survival outcome (p = 0.02). Multivariate Cox analysis indicated that risk value could be used as an independent prognostic factor, revealing that the m6A-related genes signature has great predictive value. Its efficacy was also validated by data from the Gene Expression Omnibus (GEO) database.
| S-EPMC7066421 | biostudies-literature
Project description:BackgroundThe incidence of stage I and stage II lung adenocarcinoma (LUAD) is likely to increase with the introduction of annual screening programs for high-risk individuals. We aimed to identify a reliable prognostic signature with immune-related genes that can predict prognosis and help making individualized management for patients with early-stage LUAD.MethodsThe public LUAD cohorts were obtained from the large-scale databases including 4 microarray data sets from the Gene Expression Omnibus (GEO) and 1 RNA-seq data set from The Cancer Genome Atlas (TCGA) LUAD cohort. Only early-stage patients with clinical information were included. Cox proportional hazards regression model was performed to identify the candidate prognostic genes in GSE30219, GSE31210 and GSE50081 (training set). The prognostic signature was developed using the overlapped prognostic genes based on a risk score method. Kaplan-Meier curve with log-rank test and time-dependent receiver operating characteristic (ROC) curve were used to evaluate the prognostic value and performance of this signature, respectively. Furthermore, the robustness of this prognostic signature was further validated in TCGA-LUAD and GSE72094 cohorts.ResultsA prognostic immune signature consisting of 21 immune-related genes was constructed using the training set. The prognostic signature significantly stratified patients into high- and low-risk groups in terms of overall survival (OS) in training data set, including GSE30219 (HR = 4.31, 95% CI 2.29-8.11; P = 6.16E-06), GSE31210 (HR = 11.91, 95% CI 4.15-34.19; P = 4.10E-06), GSE50081 (HR = 3.63, 95% CI 1.90-6.95; P = 9.95E-05), the combined data set (HR = 3.15, 95% CI 1.98-5.02; P = 1.26E-06) and the validation data set, including TCGA-LUAD (HR = 2.16, 95% CI 1.49-3.13; P = 4.54E-05) and GSE72094 (HR = 2.95, 95% CI 1.86-4.70; P = 4.79E-06). Multivariate cox regression analysis demonstrated that the 21-gene signature could serve as an independent prognostic factor for OS after adjusting for other clinical factors. ROC curves revealed that the immune signature achieved good performance in predicting OS for early-stage LUAD. Several biological processes, including regulation of immune effector process, were enriched in the immune signature. Moreover, the combination of the signature with tumor stage showed more precise classification for prognosis prediction and treatment design.ConclusionsOur study proposed a robust immune-related prognostic signature for estimating overall survival in early-stage LUAD, which may be contributed to make more accurate survival risk stratification and individualized clinical management for patients with early-stage LUAD.
| S-EPMC7542703 | 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
Project description:N6-methyladenosine (m6A) is the most prevalent type of RNA modification, and we hypothesized that patterns of m6A-related genes may be useful for estimating risk of lung adenocarcinoma (LUAD). An m6A-related gene set variation score (m6A-GSVS) was generated using RNA-sequencing data from LUAD patients in The Cancer Genome Atlas (TCGA). We investigated the association of m6A-GSVS with stemness, tumor mutational burden (TMB), expression of three immune checkpoints, levels of tumor-infiltrating lymphocytes (TILs), and patient prognosis. We found that m6A-GSVS was higher in LUAD than in healthy lung tissue, and it strongly correlated with stemness and TMB. Activated CD4 + T cells were more numerous in LUAD samples that had higher m6A-GSVS than in those with lower scores. Biological processes and pathways, including "Cell cycle," "DNA replication," and "RNA degradation," were significantly enriched in samples with high scores. Furthermore, m6A-GSVS was an independent prognostic indicator in LUAD. In conclusion, we proposed an m6A-GSVS in LUAD. It is a putative indicator for evaluating the ability to RNA m6A, an independent prognostic indicator and associated with tumor stemness.
| S-EPMC8295598 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common subtype of non-small cell lung cancer. Ferroptosis is a newly recognized process of cell death, which is different from other forms of cell death in morphology, biochemistry, and genetics, and has played a vital role in cancer biology. This study aimed to identify a ferroptosis-related gene signature associated with LUAD prognosis.MethodsDataset TCGA-LUAD which came from the TCGA portal was taken as the training cohort. GSE72094 and GSE68465 from the GEO database were treated as validation cohorts. Two hundred fifty-nine ferroptosis-related genes were retrieved from the FerrDb database. In the training cohort, Kaplan-Meier and univariate Cox analyses were conducted for preliminary screening of ferroptosis-related genes with potential prognostic capacity. These genes then entered into the LASSO Cox regression model, constructing a gene signature. The latter was then evaluated in the training and validation cohorts via Kaplan-Meier, Cox, and ROC analyses. In addition, the correlations between risk score and autophagy were examined by Pearson correlation coefficient. The analyses of GSEA and immune infiltrating were performed for better studying the function annotation of the gene signature and the character of each kind of immune cells played in the tumor microenvironment.ResultsA 15-gene signature was found from the training cohort and validated by Kaplan-Meier and Cox regression analyses, revealing its independent prognosis value in LUAD. Moreover, the ROC analysis was conducted, confirming a strong predictive ability that this signature owned for LUAD prognosis. One hundred fifty-one of 222 (68.01%) autophagy-related genes were discovered significantly correlated with risk scores. Analyses of GSEA and immune infiltration exhibited in detail the specific pathways that associate with the 15-gene signature and identified the crucial roles of resting mast cells and resting dendritic cells owned in the prognosis of the 15-gene signature.ConclusionIn this present study, a novel ferroptosis-related 15-gene signature (RELA, ACSL3, YWHAE, EIF2S1, CISD1, DDIT4, RRM2, PANX1, TLR4, ARNTL, LPIN1, HERPUD1, NCOA4, PEBP1, and GLS2) was built. It could accurately predict the prognosis of LUAD and was related to resting mast cells and resting dendritic cells, which provide potential for the personalized outcome prediction and the development of new therapies in LUAD population.
| S-EPMC8264775 | biostudies-literature
Project description:Background Studies have shown that the regulation of ferroptosis could be a new approach to cancer treatment and abnormal ferroptosis is closely associated with a dysregulated immune response. However, a combined signature with ferroptosis-related genes (FRGs) and immune-related genes (IRGs) is necessary to be constructed for predicting prognoses and guiding individualized precision therapy of lung adenocarcinoma (LUAD) patients. Methods In this study, based on the Cancer Genome Atlas (TCGA) cohort, prognosis-related FRGs and IRGs were first identified and incorporated into the Least Absolute Shrinkage and Selection Operator (LASSO)-Cox regression model to generate a combined signature of ferroptosis- and immune-related genes (CSFI) values to predict the overall survivals (OSs) of LUAD patients. And patients with LUAD from the Gene Expression Omnibus (GEO) database were applied for the validation set. Nomogram was constructed based on multivariate Cox regression analysis. Subsequently, ferroptosis, immunity, and gene mutation status of patients between the CSFI-high and -low groups were compared. Additionally, the enrichment pathways in CSFI-high and -low groups were explored by Gene Ontology (GO), Kyoto Gene and Genome Encyclopedia (KEGG) and Gene Set Enrichment Analysis (GSEA) analyses. Results As a result, the CSFI-low group showed a good prognosis instead of the CSFI-high group. CSFI was identified to be an independent prognosis factor for LUAD. In general, there were ferroptosis- and immune-suppressive states in CSFI-high patients. Notably, the mutation frequencies of TP53 were higher in CSFI-high patients. Conclusions In LUAD, CSFI which served as a novel classifier was offered for predicting the prognoses of patients and contributing to guiding personalized targeted therapy of patients. Therefore, based on these findings, it was believed that a synergistic treatment of ferroptosis and immunity would be more effective on LUAD patients with low CSFI values.
| S-EPMC9641098 | biostudies-literature
Project description:Aberrant regulation of m6A mRNA modification can lead to changes in gene expression, thus contributing to tumorigenesis in several types of solid tumors. In this study, by integrating analyses of m6A methylation and mRNA expression, we identified 84 m6A-regulated mRNAs in lung adenocarcinoma (LUAD). Although the m6A methylation levels of total RNA in LUAD patient tumor tissue were reduced, the majority (75.2%) of m6A-regulated mRNAs were hypermethylated. The m6A-hypermethylated mRNAs were mainly enriched in terms related to transcription factor activity. We established a 10-m6A-regulated-mRNA signature score system through least absolute shrinkage and selection operator Cox regression analysis, with its predictive value validated by Kaplan-Meier curve and time-dependent receiver operating characteristic curves. RFXAP and KHDRBS2 from the signature also exhibited an independent prognostic value. The co-expression and interaction network analyses demonstrated the strong correlation between m6A regulators and the genes in the signature, further supporting the results of the m6A methylation modification patterns. These findings highlight the potential utility of integrating multi-omics data (m6A methylation level and mRNA expression) to accurately obtain potential prognostic biomarkers, which may provide important insights into developing novel and effective therapies for LUAD.
| S-EPMC8377381 | biostudies-literature
Project description:Hepatocellular carcinoma (HCC) is among the most common types of cancer that threat the public health worldwide. N6-methyladenosine (m6A) RNA methylation is associated with cancer initiation and progression, and is dynamically regulated by m6A RNA methylation-associated genes. However, little is known about the expression status and the prognostic value of m6A-associated genes in HCC. The present study aimed to identify the expression profiling pattern and clinical significance of m6A-associated genes in HCC. Consensus clustering analysis was performed to identify the clusters of HCC with different clinical outcomes. A prognostic signature built by the least absolute shrinkage and selection operator Cox regression model was utilized to discover subtypes associated with different clinical outcomes of patients with HCC in the discovery cohort from The Cancer Genome Atlas. The differences between subgroups were characterized in terms of epigenetic dysregulation and somatic mutation frequencies. The International Cancer Genome Consortium cohort and two independent cohorts from the meta-Gene Expression Omnibus database were used for external validation. Most of the m6A-associated genes were upregulated and involved in the prognosis and malignancy of HCC. A four-gene prognostic signature revealed two HCC subtypes (namely, high- and low-risk group) that was associated with different clinical outcomes. Patients in the high-risk group were accompanied with increased epigenetic silencing and significant mutations in TP53 and FLG, while ALB was frequently mutated in the low-risk group. In conclusion, an m6A-based signature was constructed to predict the prognosis of patients with HCC, which may provide a tool for reliable prognosis assessment for clinicians, and aid clinical treatment decision-making.
| S-EPMC7859476 | biostudies-literature
Project description:BackgroundLung cancer is the most common malignant tumor in the world, and its prognosis is still not optimistic. The aim of this study was to establish an immune-related gene (IRG) prognostic index (IRGPI) for lung adenocarcinoma (LUAD) based on IRGs, and to explore the prognosis, molecular and immune features, and response to immune checkpoint inhibitor (ICI) therapy in IRGPI-classified different subgroups of LUAD.MethodsBased on the LUAD transcriptome RNA-sequencing data in TCGA database, the differentially expressed genes (DEGs) were selected. Subsequently, DEGs were intersected with IRGs to obtain differentially expressed immune-related genes (DEIRGs). Weighted gene co-expression network analysis (WGCNA) identified hub genes in DEIRGs. Finally, univariate and multivariate Cox regression analyses were used to build an IRGPI model. Subsequently, TCGA patients were divided into high- and low-risk groups, and the survival of patients in different groups was further analyzed. Besides, we validated the molecular and immune characteristics, relationship with immune checkpoints, angiogenesis-related genes, and immune subtypes distribution in different subgroups. Meanwhile, we further validated the response to ICI therapy in different subgroups.ResultsThe IRGPI was constructed based on 13 DEIRGs. Compared with the low-risk group, overall survival (OS) was lower in the high-risk group, and the high-risk score was independently associated with poorer OS. Besides, the high-risk score was associated with cell cycle pathway, high mutation rate of TP53 and KRAS, high infiltration of M0 macrophages, and immunosuppressive state, and these patients had poorer prognosis but the TIDE score of the high-risk group was lower than that of the other group, which means that the high-risk group could benefit more from ICI treatment. In contrast, the low-risk score was related to low mutation rate of TP53 and KRAS, high infiltration of plasma cells, and immunoactive state, and these patients had better prognosis but the low-risk group less benefit from ICI treatment based on the results of TIDE score.ConclusionsIRGPI is a prospective biomarker based on IRGs that can distinguish high- and low-risk groups to predict patient prognosis, help characterize the tumor immune microenvironment, and evaluate the benefit of ICI therapy in LUAD.
| S-EPMC10713328 | biostudies-literature