Dataset Information

The role of m6A-related genes in the prognosis and immune microenvironment of pancreatic adenocarcinoma.

ABSTRACT:

Background

Pancreatic adenocarcinoma (PAAD) is among the most lethal diseases and has a dismal prognosis; however, efficient treatment is currently limited. Several studies have observed epigenetic variation during tumorigenesis, suggesting the potential role of RNA methylation, especially N6-methyladenosine (m6A) modification, as a novel epigenetic modification mediating PAAD prognosis.

Methods

The expression levels of m6A-related genes were downloaded from The Cancer Genome Atlas-Pancreatic Adenocarcinoma (TCGA) and Genotype-Tissue Expression (GTEx) projects, and the findings were validated in four Expression Omnibus (GEO) datasets. A predictive model was constructed using a lasso regression and evaluated by a survival analysis and receiver operating characteristic curve. Consensus clustering identified two distinct subgroups with different immune activity signatures based on the expression pattern of m6A-related genes. The relationship between the mutation state of m6A-related genes and infiltration of immune cells was established and visualized using Tumor Immune Estimation Resource (https://cistrome.shinyapps.io/timer/).

Results

Fourteen of twenty-one m6A-related genes were differentially expressed between PAAD and normal tissues in TCGA-GTEx cohort. Among these genes, HNRNPC, IGF2BP2 and YTHDF1 were further validated in four GEO datasets. Moreover, an m6A-based model exhibited moderate accuracy in predicting overall survival in PAAD samples. Additionally, potential m6A modification targets were screened by selecting genes from a set of 23,391 genes that not only harbored the most m6A-modified sites but also showed a robust correlation with PAAD survival. Moreover, we correlated the expression level of m6A-related genes with the immune microenvironment of pancreatic cancer for the first time. Specifically, both arm-level gain and deletion of ALKBH5 decreased the infiltration of CD8+T cells (P < 0.05 and P < 0.01, respectively).

Conclusion

Collectively, our findings suggest a novel anticancer strategy for restoring balanced RNA methylation in tumor cells and guide clinical physicians in developing a new practical approach for considering the impact of related genes on prognosis.

SUBMITTER: Tang R

PROVIDER: S-EPMC7528816 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive, fatal tumor. N6-methylandenosine (m6A) methylation is the major epigenetic modification of RNA including lncRNAs. The roles of m6A-related lncRNAs in PDAC have not been fully clarified. This study aims to assess gene signatures and prognostic value of m6A-related lncRNAs in PDAC. The Cancer Genome Atlas (TCGA) dataset and the International Cancer Genome Consortium (ICGC) dataset were explored to identify m6A-related lncRNAs. Univariate, least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression were performed to construct the m6A-related lncRNAs prognostic riskscore (m6A-LPR) model to predict the overall survival (OS) in the TCGA training cohort. Kaplan-Meier curve with log-rank test and receiver operating characteristic (ROC) curve were used to evaluate the prognostic value of the m6A-LPR. Furthermore, the robustness of the m6A-LPR was further validated in the ICGC cohort. Tumor immunity was evaluated using ESTIMATE and CIBERSORT algorithms. A total of 262 m6A-related lncRNAs were identified in two datasets. In the TCGA training cohort, 28 prognostic m6A-related lncRNAs were identified and the m6A-LPR including four m6A-related lncRNAs was constructed. The m6A-LPR was able to identify high-risk patients with significantly poorer OS and accurately predict OS in both the TCGA training cohort and the ICGC validation cohort. Analysis of tumor immunity revealed that high-risk groups had remarkably lower stromal, immune, and ESTIMATE scores. Moreover, high-risk groups were associated with significantly higher levels of plasma B cells and resting NK cells infiltration, and lower levels of infiltrating resting memory CD4 T cells, monocytes, and resting mast cells. Our study proposed a robust m6A-related prognostic signature of lncRNAs for predicting OS in PDAC, which provides some clues for further studies focusing on the mechanism process underlying m6A modification of lncRNAs.

Project description:BackgroundPancreatic ductal adenocarcinoma (PDAC) is a malignancy characterized by challenging early diagnosis and poor prognosis. It is believed that coagulation has an impact on the tumor microenvironment of PDAC. The aim of this study is to further distinguish coagulation-related genes and investigate immune infiltration in PDAC.MethodsWe gathered two subtypes of coagulation-related genes from the KEGG database, and acquired transcriptome sequencing data and clinical information on PDAC from The Cancer Genome Atlas (TCGA) database. Using an unsupervised clustering method, we categorized patients into distinct clusters. We investigated the mutation frequency to explore genomic features and performed enrichment analysis, utilizing Gene Ontology (GO) and Kyoto Encyclopedia of Genes (KEGG) to explore pathways. CIBERSORT was used to analyze the relationship between tumor immune infiltration and the two clusters. A prognostic model was created for risk stratification, and a nomogram was established to assist in determining the risk score. The response to immunotherapy was assessed using the IMvigor210 cohort. Finally, PDAC patients were recruited, and experimental samples were collected to validate the infiltration of neutrophils using immunohistochemistry. In addition, and identify the ITGA2 expression and function were identified by analyzing single cell sequencing data.ResultsTwo coagulation-related clusters were established based on the coagulation pathways present in PDAC patients. Functional enrichment analysis revealed different pathways in the two clusters. Approximately 49.4% of PDAC patients experienced DNA mutation in coagulation-related genes. Patients in the two clusters displayed significant differences in terms of immune cell infiltration, immune checkpoint, tumor microenvironment and TMB. We developed a 4-gene prognostic stratified model through LASSO analysis. Based on the risk score, the nomogram can accurately predict the prognosis in PDAC patients. We identified ITGA2 as a hub gene, which linked to poor overall survival (OS) and short disease-free survival (DFS). Single-cell sequencing analysis demonstrated that ITGA2 was expressed by ductal cells in PDAC.ConclusionsOur study demonstrated the correlation between coagulation-related genes and the tumor immune microenvironment. The stratified model can predict the prognosis and calculate the benefits of drug therapy, thus providing the recommendations for clinical personalized treatment.

Project description:N6-methylandrostenedione (m6A) methylation plays a very important role in the development of malignant tumors. The immune system is the key point in the progression of tumors, particularly in terms of tumor treatment and drug resistance. Tumor immunotherapy has now become a hot spot and a new approach for tumor treatment. However, as far as the stomach adenocarcinoma (STAD) is concerned, the in-depth research is still a gap in the m6A-associated immune markers. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases is extremely important for our research, where we obtained gene mutation, gene expression data and relevant clinical information of STAD patients. Firstly, the samples from GEO were used as external validation groups, while the TCGA samples were divided into a training group and an internal validation group randomly. Using the way of Single factor COX-LASSO- and multi-factor Cox to construct the prognostic model. Then, all samples were subjected to cluster analysis to generate high and low expression groups of immune gene. Meanwhile, we also collected the correlation between these types and tumor microenvironment. On this basis, a web version of the dynamic nomogram APP was developed. In addition, we performed microenvironmental correlation, copy number variation and mutation analyses for model genes. The prognostic model for STAD developed here demonstrated a very strong predictive ability. The results of cluster analysis manifested that the immune gene low expression group had lower survival rate and higher degree of immune infiltration. Therefore, the immune gene low expression group was associated with lower survival rates and a higher degree of immune infiltration. Gene set enrichment analysis suggested that the potential mechanism might be related to the activation of immunosuppressive functions and multiple signaling pathways. Correspondingly, the web version of the dynamic nomogram APP produced by the DynNom package has successfully achieved rapid and accurate calculation of patient survival rates. Finally, the multi-omics analysis of model genes further enriched the research content. Interference of RAB19 was confirmed to facilitate migration of STAD cells in vitro, while its overexpression inhibited these features. The prognostic model for STAD constructed in this study is accurate and efficient based on multi-omics analysis and experimental validation. Additionally, the results of the correlation analysis between the tumor microenvironment and m6Ascore are the basics of further exploration of the pathophysiological mechanism in STAD.

Project description:Pancreatic adenocarcinoma (PAAD) is one of the deadliest malignancies. Aging is described as the degeneration of physiological function, which is complexly correlated with cancer. It is significant to explore the influences of aging-related genes (ARGs) on PAAD. Based on The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) datasets, we used univariate Cox regression analysis and acquired eight differentially expressed ARGs with prognostic values. Two molecular subtypes were identified based on these ARGs to depict PAAD patients' overall survival (OS) and immune microenvironments preliminarily. Cluster 1 had a poor OS as well as a worse immune microenvironment. Through least absolute shrinkage and selection operator (LASSO) regression analysis, we constructed a seven-ARG risk signature based on the TCGA dataset and verified it in Gene Expression Omnibus (GEO) and International Cancer Genome Consortium (ICGC) to predict the prognoses, immune microenvironments, signal pathways, tumor mutations, and drug sensitivity of PAAD patients. The high-risk group possessed an unfavorable OS compared with that of the low-risk group. We also verified the independence and clinical availability of the risk signature by Cox regression analyses and the establishment of a nomogram, respectively. The higher risk score was associated with several clinical factors such as higher grade and advanced tumor stage as well as lower immunoscore and cluster 1. The negative associations of risk scores with immune, stroma, and estimate scores proved the terrible immune microenvironment in the high-risk group. Relationships between risk score and immune checkpoint gene expression as well as signal pathways provided several therapeutic targets. PAAD patients in the low-risk group possessed lower tumor mutations as well as a higher susceptibility to axitinib and vorinostat. The high-risk group bore a higher TMB and cisplatin and dasatinib may be better options. We used immunohistochemistry and qPCR to confirm the expression of key ARGs with their influences on OS. In conclusion, we identified two ARG-mediated molecular subtypes and a novel seven-ARG risk signature to predict prognoses, immune microenvironments, signal pathways, tumor mutations, and drug sensitivity of PAAD patients.

Project description:BackgroundLung adenocarcinoma (LUAD) is the most common subtype of lung cancer which typically exhibits a diverse progression trajectory. Our study sought to explore the cell differentiation trajectory of LUAD and its clinical relevance.MethodsUtilizing a single-cell RNA-sequencing dataset (GSE117570), we identified LUAD cells of distinct differential status along with differentiation-related genes (DRGs). DRGs were applied to the analysis of bulk-tissue RNA-sequencing dataset (GSE72094) to classify tumors into different subtypes, whose clinical relevance was further analyzed. DRGs were also applied to gene co-expression network analysis (WGCNA) using another bulk-tissue RNA-sequencing dataset (TCGA-LUAD). Genes from modules that demonstrated a significant correlation with clinical traits and were differentially expressed between normal tissue and tumors were identified. Among these, genes with significant prognostic relevance were used for the development of a prognostic nomogram, which was tested on TCGA-LUAD dataset and validated in GSE72094. Finally, CCK-8, EdU, cell apoptosis, cell colony formation, and Transwell assays were used to verify the functions of the identified genes.ResultsFour clusters of cells with distinct differentiation status were characterized, whose DRGs were predominantly correlated with pathways of immune regulation. Based on DRGs, tumors could be clustered into four subtypes associated with distinct immune microenvironment and clinical outcomes. DRGs were categorized into four modules. A total of nine DRGs (SFTPB, WFDC2, HLA-DPA1, TIMP1, MS4A7, HLA-DQA1, VCAN, KRT8, and FABP5) with most significant survival-predicting power were integrated to develop a prognostic model, which outperformed the traditional parameters in predicting clinical outcomes. Finally, we verified that knockdown of WFDC2 inhibited proliferation, migration, and invasion but promoted the apoptosis of A549 cells in vitro.ConclusionThe cellular composition and cellular differentiation status of tumor mass can predict the clinical outcomes of LUAD patients. It also plays an important role in shaping the tumor immune microenvironment.

Project description:Backgroundm6A is the most prevalent and abundant form of mRNA modifications and is closely related to tumor proliferation, differentiation, and tumorigenesis. In this study, we try to conduct an effective prediction model to investigated the function of m6A RNA methylation modulators in pancreatic adenocarcinoma and estimated the potential association between m6A RNA methylation modulators and tumor microenvironment infiltration for optimization of treatment.MethodsExpression of 28 m6A RNA methylation modulators and clinical data of patients with pancreatic adenocarcinoma and normal samples were obtained from TCGA and GTEx database. Differences in the expression of 28 m6A RNA methylation modulators between tumour (n = 40) and healthy (n = 167) samples were compared by Wilcoxon test. LASSO Cox regression was used to select m6A RNA methylation modulators to analyze the relationship between expression and clinical characteristics by univariate and multivariate regression. A risk score prognosis model was conducted based on the expression of select m6A RNA methylation modulators. Bioinformatics analysis was used to explore the association between the m6Ascore and the composition of infiltrating immune cells between high and low m6Ascore group by CIBERSORT algorithm. Evaluation of m6Ascore for immunotherapy was analyzed via the IPS and three immunotherapy cohort. Besides, the biological signaling pathways of the m6A RNA methylation modulators were examined by gene set enrichment analysis (GSEA).ResultsExpression of 28 m6A RNA methylation modulators were upregulated in patients with PAAD except for MTEEL3. An m6Ascore prognosis model was established, including KIAA1429, IGF2BP2, IGF2BP3, METTL3, EIF3H and LRPPRC was used to predict the prognosis of patients with PAAD, the high risk score was an independent prognostic indicator for pancreatic adenocarcinoma, and a high risk score presented a lower overall survival. In addition, m6Ascore was related with the immune cell infiltration of PAAD. Patients with a high m6Ascore had lower infiltration of Tregs and CD8+T cells but a higher resting CD4+ T infiltration. Patients with a low m6Ascore displayed a low abundance of PD-1, CTLA-4 and TIGIT, however, the IPS showed no difference between the two groups. The m6Ascore applied in three immunotherapy cohort (GSE78220, TCGA-SKCM, and IMvigor210) did not exhibit a good prediction for estimating the patients' response to immunotherapy, so it may need more researches to figure out whether the m6A modulator prognosis model would benefit the prediction of pancreatic patients' response to immunotherapy.ConclusionModulators involved in m6A RNA methylation were associated with the development of pancreatic cancer. An m6Ascore based on the expression of IGF2BP2, IGF2BP3, KIAA1429, METTL3, EIF3H and LRPPRC is proposed as an indicator of TME status and is instrumental in predicting the prognosis of pancreatic cancer patients.

Project description:BackgroundCuproptosis plays a crucial role in cancer, and different subtypes of cuproptosis have different immune profiles in prostate adenocarcinoma (PRAD). This study aimed to investigate immune genes associated with cuproptosis and develop a risk model to predict prognostic characteristics and chemotherapy/immunotherapy responses of patients with PRAD.MethodsThe CIBERSORT algorithm was used to evaluate the immune and stromal scores of patients with PRAD in The Cancer Genome Atlas (TCGA) cohort. Validation of differentially expressed genes DLAT and DLD in benign and malignant tissues by immunohistochemistry, and the immune-related genes of DLAT and DLD were further screened. Univariable Cox regression were performed to select key genes. Least absolute shrinkage and selection operator (LASSO)-Cox regression analyse was used to develop a risk model based on the selected genes. The model was validated in the TCGA, Memorial Sloan-Kettering Cancer Center (MSKCC) and Gene Expression Omnibus (GEO) datasets, as well as in this study unit cohort. The genes were examined via functional enrichment analysis, and the tumor immune features, tumor mutation features and copy number variations (CNVs) of patients with different risk scores were analysed. The response of patients to multiple chemotherapeutic/targeted drugs was assessed using the pRRophetic algorithm, and immunotherapy was inferred by the Tumor Immune Dysfunction and Exclusion (TIDE) and immunophenoscore (IPS).ResultsCuproptosis-related immune risk scores (CRIRSs) were developed based on PRLR, DES and LECT2. High CRIRSs indicated poor overall survival (OS), disease-free survival (DFS) in the TCGA-PRAD, MSKCC and GEO datasets and higher T stage and Gleason scores in TCGA-PRAD. Similarly, in the sample collected by the study unit, patients with high CRIRS had higher T-stage and Gleason scores. Additionally, higher CRIRSs were negatively correlated with the abundance of activated B cells, activated CD8+ T cells and other stromal or immune cells. The expression of some immune checkpoints was negatively correlated with CRIRSs. Tumor mutational burden (TMB), mutant-allele tumor heterogeneity (MATH) and copy number variation (CNV) scores were all higher in the high-CRIRS group. Multiple chemotherapeutic/targeted drugs and immunotherapy had better responsiveness in the low-CRIRS group.ConclusionOverall, lower CRIRS indicated better response to treatment strategies and better prognostic outcomes.