Comprehensive Analysis of APA Events and Their Association With Tumor Microenvironment in Lung Adenocarcinoma.
Ontology highlight
ABSTRACT: Background
Alternative polyadenylation (APA) is a pervasive posttranscriptional mechanism regulating gene expression. However, the specific dysregulation of APA events and its potential biological or clinical significance in lung adenocarcinoma (LUAD) remain unclear.Methods
Here, we collected RNA-Seq data from two independent datasets: GSE40419 (n = 146) and The Cancer Genome Atlas (TCGA) LUAD (n = 542). The DaPars algorithm was employed to characterize the APA profiles in tumor and normal samples. Spearman correlation was used to assess the effects of APA regulators on 3' UTR changes in tumors. The Cox proportional hazard model was used to identify clinically relevant APA events and regulators. We stratified 512 patients with LUAD in the TCGA cohort through consensus clustering based on the expression of APA factors.Findings
We identified remarkably consistent alternative 3' UTR isoforms between the two cohorts, most of which were shortened in LUAD. Our analyses further suggested that aberrant usage of proximal polyA sites resulted in escape from miRNA binding, thus increasing gene expression. Notably, we found that the 3' UTR lengths of the mRNA transcriptome were correlated with the expression levels of APA factors. We further identified that CPSF2 and CPEB3 may serve as key regulators in both datasets. Finally, four LUAD subtypes according to different APA factor expression patterns displayed distinct clinical results and oncogenic features related to tumor microenvironment including immune, metabolic, and hypoxic status.Interpretation
Our analyses characterize the APA profiles among patients with LUAD and identify two key regulators for APA events in LUAD, CPSF2 and CPEB3, which could serve as the potential prognostic genes in LUAD.
SUBMITTER: Zhang Y
PROVIDER: S-EPMC8012674 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Lung adenocarcinoma (LUAD) is the most common malignancy, leading to more than 1 million related deaths each year. Due to low long-term survival rates, the exploration of molecular mechanisms underlying LUAD progression and novel prognostic predictors is urgently needed to improve LUAD treatment. In our study, cancer-specific differentially expressed genes (DEGs) were identified using the robust rank aggregation (RRA) method between tumor and normal tissues from six Gene Expression Omnibus databases (GSE43458, GSE62949, GSE68465, GSE115002, GSE116959, and GSE118370), followed by a selection of prognostic modules using weighted gene co-expression network analysis. Univariate Cox regression, least absolute shrinkage and selection operator (LASSO), and multivariate Cox regression analyses were applied to identify nine hub genes (CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6) that constructed a prognostic risk model. The RNA expressions of nine hub genes were validated in tumor and normal tissues by RNA-sequencing and single-cell RNA-sequencing, while immunohistochemistry staining from the Human Protein Atlas database showed consistent results in the protein levels. The risk model revealed that high-risk patients were associated with poor prognoses, including advanced stages and low survival rates. Furthermore, a multivariate Cox regression analysis suggested that the prognostic risk model could be an independent prognostic factor for LUAD patients. A nomogram that incorporated the signature and clinical features was additionally built for prognostic prediction. Moreover, the levels of hub genes were related to immune cell infiltration in LUAD microenvironments. A CMap analysis identified 13 small molecule drugs as potential agents based on the risk model for LUAD treatment. Thus, we identified a prognostic risk model including CBFA2T3, CR2, SEL1L3, TM6SF1, TSPAN32, ITGA6, MAPK11, RASA3, and TLR6 as novel biomarkers and validated their prognostic and predicted values for LUAD.
| S-EPMC8440811 | biostudies-literature
Project description:Background: Trace elements such as copper are essential for human health. Recently the journal Nat Rev Cancer has put forward the concept of Cuproplasia, a way of promoting tumor growth through reliance on copper. We attempted to conduct a comprehensive analysis of Cuproplasia-related genes in lung adenocarcinoma (LUAD) to explore the mechanism of action of Cuproplasia-related genes in LUAD. Method: Transcriptome data and clinical information of LUAD were obtained from TCGA-LUAD and GSE31210, and prognostic models of Cuproplasia-related genes were constructed and verified by regression analysis of GSVA, WGCNA, univariate COX and lasso. The signal pathways affected by Cuproplasia-related genes were analyzed by GO, KEGG and hallmarK pathway enrichment methods. Five immunocell infiltration algorithms and IMVIGOR210 data were used to analyze immune cell content and immunotherapy outcomes in the high-low risk group. Results: In the results of WGCNA, BROWN and TURQUOISE were identified as modules closely related to Cuproplasia score. In the end, lasso regression analysis established a Cuproplasia-related signature (CRS) based on 24 genes, and the prognosis of high-risk populations was worse in TCGA-LUAD and GSE31210 datasets. The enrichment analysis showed that copper proliferation was mainly through chromosome, cell cycle, dna replication, g2m checkpoint and other pathways. Immunoinfiltration analysis showed that there were differences in the content of macrophages among the four algorithms. And IMVIGOR210 found that the lower the score, the more effective the immunotherapy was. Conclusion: The Cuproplasia related gene can be used to predict the prognosis and immunotherapy outcome of LUAD patients, and may exert its effect by affecting chromosome-related pathways and macrophages.
| S-EPMC10540310 | biostudies-literature
Project description:Lung adenocarcinoma (LUAD) was the first one all over the world. RAB11FIP1 was found to be expressed differently in a critical way among different cancers. However, the prognostic value and immune infiltration of RAB11FIP1 expression in LUAD are unclear. In this study, the expression of RAB11FIP1 in LUAD was investigated in the Oncomine, TCGA, GEO, and UALCAN databases. Kaplan-Meier analysis was chosen to compare the association between RAB11FIP1 expression and overall survival (OS) in LUAD patients. The dataset of TCGA was used to analyze the pertinence between RAB11FIP1 and clinicpathological factors. GO, KEGG, and network analysis of protein-protein interactions (PPI) were conducted to investigate the potential mechanism of RAB11FIP1. In the end, the relevance of RAB11FIP1 to cancer-immune infiltrates was investigated. RAB11FIP1 was found to be down-regulated by tumors compared with adjacent normal tissue in multiple LUAD cohorts. RAB11FIP1 is an independent prognostic factor in lung adenocarcinoma. There was a high correlation between low RAB11FIP1 in tumors and worse OS in LUAD. Functional network analysis suggested that RAB11FIP1 was associated with multiple pathways. Besides, the expression of RAB11FIP1 was closely related to the infiltration levels of B cell, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells. RAB11FIP1 expression in LUAD occurred with a variety of immune markers. Our findings suggest that RAB11FIP1 is related to prognosis and immune infiltrates in LUAD.
| S-EPMC8576257 | biostudies-literature
Project description:Background:Prognostic genes in the tumor microenvironment play an important role in immune biological processes and the response of cancer to immunotherapy. Thus, we aimed to assess new biomarkers that are associated with immune/stromal cells in lung adenocarcinomas (LUAD) using the ESTIMATE algorithm, which also significantly affects the prognosis of cancer. Methods:The RNA sequencing (RNA-Seq) and clinical data of LUAD were downloaded from the the Cancer Genome Atlas (TCGA ). The immune and stromal scores were calculated for each sample using the ESTIMATE algorithm. The LUAD gene chip expression profile data and the clinical data (GSE37745, GSE11969, and GSE50081) were downloaded from the Gene Expression Omnibus (GEO) for subsequent validation analysis. Differentially expressed genes were calculated between high and low score groups. Univariate Cox regression analysis was performed on differentially expressed genes (DEGs) between the two groups to obtain initial prognosis genes. These were verified by three independent LUAD cohorts from the GEO database. Multivariate Cox regression was used to identify overall survival-related DEGs. UALCAN and the Human Protein Atlas were used to analyze the mRNA /protein expression levels of the target genes. Immune cell infiltration was evaluated using the Tumor Immune Estimation Resource (TIMER) and CIBERSORT methods, and stromal cell infiltration was assessed using xCell. Results:In this study, immune scores and stromal scores are significantly associated with the clinical characteristics of LUAD, including T stage, M stage, pathological stage, and overall survival time. 530 DEGs (18 upregulated and 512 downregulated) were found to coexist in the difference analysis with the immune scores and stromal scores subgroup. Univariate Cox regression analysis showed that 286 of the 530 DEGs were survival-related genes (p < 0.05). Of the 286 genes initially identified, nine prognosis-related genes (CSF2RB, ITK, FLT3, CD79A, CCR4, CCR6, DOK2, AMPD1, and IGJ) were validated from three separate LUAD cohorts. In addition, functional analysis of DEGs also showed that various immunoregulatory molecular pathways, including regulation of immune response and the chemokine signaling pathways, were involved. Five genes (CCR6, ITK, CCR4, DOK2, and AMPD1) were identified as independent prognostic indicators of LUAD in specific data sets. The relationship between the expression levels of these genes and immune genes was assessed. We found that CCR6 mRNA and protein expression levels of LUAD were greater than in normal tissues. We evaluated the infiltration of immune cells and stromal cells in groups with high and low levels of expression of CCR6 in the TCGA LUAD cohort. In summary, we found a series of prognosis-related genes that were associated with the LUAD tumor microenvironment.
| S-EPMC7382940 | biostudies-literature
Project description:Because of immunotherapy failure in lung adenocarcinoma, we have tried to find new potential biomarkers for differentiating different tumor subtypes and predicting prognosis. We identified two subtypes based on tumor microenvironment-related genes in this study. We used seven methods to analyze the immune cell infiltration between subgroups. Further analysis of tumor mutation load and immune checkpoint expression among different subgroups was performed. The least absolute shrinkage and selection operator Cox regression was applied for further selection. The selected genes were used to construct a prognostic 14-gene signature for LUAD. Next, a survival analysis and time-dependent receiver operating characteristics were performed to verify and evaluate the model. Gene set enrichment analyses and immune analysis in risk groups was also performed. According to the expression of genes related to the tumor microenvironment, lung adenocarcinoma can be divided into cold tumors and hot tumors. The signature we built was able to predict prognosis more accurately than previously known models. The signature-based tumor microenvironment provides further insight into the prediction of lung adenocarcinoma prognosis and may guide individualized treatment.
| S-EPMC9283441 | biostudies-literature
Project description:BackgroundSuccinate dehydrogenase (SDH), one of the key enzymes in the tricarboxylic acid cycle, is mainly found in the mitochondria. SDH consists of four subunits encoding SDHA, SDHB, SDHC, and SDHD. The biological function of SDH is significantly related to cancer progression. Colorectal cancer (CRC) is one of the most common malignant tumors globally, whose most common histological subtype is colon adenocarcinoma (COAD). However, the correlation between SDH factors and COAD remains unclear.MethodsThe data on pan-cancer was obtained from The Cancer Genome Atlas (TCGA) database. Kaplan-Meier survival analysis showed the prognostic ability of SDHs. The cBioPortal database reflected genetic variations of SDHs. The correlation analysis was conducted between SDHs and mitochondrial energy metabolism genes (MMGs) and the protein-protein interaction (PPI) network was built. Consequently, Univariate and Multivariate Cox Regression Analysis on SDHs and other clinical characteristics were conducted. A nomogram was established. The ssGSEA analysis visualized the association between SDHs and immune infiltration. Immunophenoscore (IPS) explored the correlation between SDHs and immunotherapy, and the correlation between SDHs and targeted therapy was investigated through Genomics of Drug Sensitivity in Cancer. Finally, qPCR and immunohistochemistry detected SDHs' expression.ResultsAfter assessing SDHs differential expression in pan-cancer, we found that SDHB, SDHC, and SDHD benefit COAD patients. The cBioPortal database demonstrated that SDHA was the top gene in mutation frequency rank. Correlation analysis mirrored a strong link between SDHs and MMGs. We formulated a nomogram and found that SDHB, SDHC, SDHD, and clinical characteristics correlated with COAD patients' survival. For T helper cells, Th2 cells, and Tem, SDHA, SDHB, SDHC, and SDHD were significantly enriched in the high expression group. Moreover, COAD patients with high SDHA expression were more suitable for immunotherapy. And COAD patients with different SDHs' expression have different sensitivity to targeted drugs. Further verifying the gene and protein expression levels of SDHs, we found that the tissues were consistent with the bioinformatics analysis.ConclusionsOur study analyzed the expression and prognostic value of SDHs in COAD, explored the pathway mechanisms involved, and the immune cell correlations, indicating that SDHs might be biomarkers for COAD patients.
| S-EPMC10025334 | biostudies-literature
Project description:BackgroundThe N6-methyladenosine (m6A) RNA modification is the most prevalent and abundant type found in eukaryotic cells. It plays a crucial role in the initiation and progression of cancers. In this study, we aimed to comprehensively investigate the landscape of m6A regulators and their association with tumor microenvironment (TME), immunotherapeutic strategies in colon adenocarcinoma (COAD).ResultsThe differential expression, mutation, CNV frequency and prognostic value of 27 m6A regulators were systematically analyzed in COAD. Patients were classified into two clusters based on m6A regulators through consistent clustering analysis, with cluster A showing significant survival benefits. Most of the m6A regulators were negatively correlated with immune cells, except for WTAP, IGF2BP3, FTO, ALKBH5, which showed a positive correlation. We developed an m6A scoring system to calculate the m6Ascore for each patient. Patients with a high-m6Ascore had a better outcome, with the AUC of 0.775. An independent cohort of 416 COAD patients acquired from GSE38832 database was used to validate the prognosis prediction ability of m6Ascore. Moreover, the m6Ascore was negatively correlated with infiltration of anti-tumor immune cells. Additionally, patients with a high-m6Ascore responded better to anti-PD1 and anti-CTLA4 therapies, and those with MSI-H had a higher m6Ascore. Finally, we investigated the value of m6Ascore in predicting the response of patients to 15 commonly used drugs.ConclusionsWe comprehensively analyzed m6A regulators in COAD, including RNA expression, CNV changes, mutations and their correlation with TME. Our results showed that the m6A scoring system had significant predictive power for the prognosis of COAD patients, potentially leading to new personalized immunotherapy strategies.
| S-EPMC10422724 | biostudies-literature
Project description:Immune checkpoint blockade (ICB) therapies that target programmed cell death 1 (PD1) and PD1 ligand 1 (PDL1) have demonstrated promising benefits in lung adenocarcinoma (LUAD), and tumor mutational burden (TMB) is the most robust biomarker associated with the efficacy of PD-1-PD-L1 axis blockade in LUAD, but the assessment of TMB by whole-exome sequencing (WES) is rather expensive and time-consuming. Although targeted panel sequencing has been developed and approved by the US Food and Drug Administration (FDA) to estimate TMB, we found that its predictive accuracy for ICB response was significantly lower than WES in LUAD. Given that previous studies were mainly focusing on genomic variations to explore surrogate biomarkers of TMB, we turned to examine the transcriptome-based correlation with TMB in this study. Combining three immunotherapeutic cohorts with two independent The Cancer Genome Atlas (TCGA) datasets, we revealed that the expression of mutS homolog 2 (MSH2), one of the most crucial genes involved in DNA mismatch repair (MMR) pathway, was the strongest feature associated with increased TMB in multivariate analysis. Furthermore, MSH2 expression also displayed a significantly positive correlation with smoking signature while an inverse association with MMR deficiency (MMRd) signature in LUAD. More importantly, high expression of MSH2 markedly correlated with increased PD-L1 expression and CD8+ T cell infiltration, both suggesting a prominent immunotherapy-responsive microenvironment in LUAD. Notably, detecting MSH2 expression is much easier, faster, and cheaper than TMB in clinical practice. Taken together, this study demonstrates the association of MSH2 expression with TMB and the immune microenvironment in LUAD. MSH2 expression may be developed as a potential surrogate biomarker of TMB to identify ICB responders in LUAD.
| S-EPMC7046689 | biostudies-literature
Project description:Lung adenocarcinoma has entered into an era of immunotherapy with the development of immune checkpoint inhibitors (ICIs). The identification of immune subtype is crucial to prolonging survival in patients. The tumor microenvironment (TME) and metabolism have a profound impact on prognosis and therapy. The majority of previous studies focused on only one aspect, while both of them are essential to the understanding of tumorigenesis and development. We hypothesized that lung adenocarcinoma can be stratified into immune subgroups with alterations in the TME infiltration. We aimed to explore the "TME-Metabolism-Risk" patterns in each subtypes and the mechanism behind. Glycolysis and cholesterol were selected for the analysis of metabolic states based on the first half of the study. Bioinformatic analysis was performed to investigate the transcriptomic and clinical data integrated by three lung adenocarcinoma cohorts (GSE30219, GSE31210, GSE37745, N = 415). The results were validated in an independent cohort (GSE50081, N = 127). In total, 415 lung adenocarcinoma samples were integrated and analyzed. Four major immune subtypes were indentified using bioinformatic analysis. Subtype NC1, characterized by a high level of glycolysis, with extremely low microenvironment cell infiltration. Subtype NC2, characterized by the "Silence" and "Cholesterol biosynthesis Predominant" metabolic states, with a middle degree infiltration of microenvironment cell. Subtype NC3, characterized by the lack of "Cholesterol biosynthesis Predominant" metabolic state, with abundant microenvironment cell infiltration. Subtype NC4, characterized by "Mixed" metabolic state, with a relatively low microenvironment cell infiltration. Least absolute shrinkage and selection operator (LASSO) regression and multivariate analyses were performed to calculate the risk of each sample, and we attempted to find out the potential immune escape mechanism in different subtypes. The result revealed that the lack of immune cells infiltration might contribute to the immune escape in subtypes NC1 and NC4. NC3 was characterized by the high expression of immune checkpoint molecules and fibroblasts. NC2 had defects in activation of innate immune cells. There existed an obviously survival advantage in subtype NC2. Gene set enrichment analysis (GSEA) and Gene Ontology analysis indicated that the PI3K-AKT-mTOR, TGF-β, MYC-related pathways might be correlated with this phenomenon. In addition, some differentially expressed genes (DEGs) were indentified in subtype NC3, which might be potential targets for survival phenotype transformation.
| S-EPMC8193848 | biostudies-literature
Project description:BackgroundLung adenocarcinoma (LUAD) is the most common histological type of lung cancers, which is the primary cause of cancer-related mortality worldwide. Growing evidence has suggested that tumor microenvironment (TME) plays a pivotal role in tumorigenesis and progression. Hence, we investigate the correlation of TME related genes with LUAD prognosis.MethodThe information of LUAD gene expression data was obtained from The Cancer Genome Atlas (TCGA). According to their immune/stromal scores calculated by the ESTIMATE algorithm, differentially expressed genes (DEGs) were identified. Then, we performed univariate Cox regression analysis on DEGs to obtain genes that are apparently bound up with LUAD survival (SurGenes). Functional annotation and protein-protein interaction (PPI) was also conducted on SurGenes. By validating the SurGenes with data sets of lung cancer from the Gene Expression Omnibus (GEO), 106 TME related SurGenes were generated. Further, intersection analysis was executed between the 106 TME related SurGenes and hub genes from PPI network, PTPRC and CD19 were obtained. Gene Set Enrichment Analysis and CIBERSORT analysis were performed on PTPRC and CD19. Based on the TCGA LUAD dataset, we conducted factor analysis and Step-wise multivariate Cox regression analysis for 106 TME related SurGenes to construct the prognostic model for LUAD survival prediction. The LUAD dataset in GEO (GSE68465) was used as the testing dataset to confirm the prognostic model. Multivariate Cox regression analysis was used between risk score from the prognostic model and clinical parameters.ResultA total of 106 TME related genes were collected in our research totally, which were markedly correlated with the overall survival (OS) of LUAD patient. Bioinformatics analysis suggest them mainly concentrated on immune response, cell adhesion, and extracellular matrix. More importantly, among 106 TME related SurGenes, PTPRC and CD19 were highly interconnected nodes among PPI network and correlated with immune activity, exhibiting significant prognostic potential. The prognostic model was a weighted linear combination of the 106 genes, by which the low-OS LUAD samples could be separated from the high-OS samples with success. This model was also able to rebustly predict the situation of survival (training set: p-value < 0.0001, area under the curve (AUC) = 0.649; testing set: p-value = 0.0009, AUC = 0.617). By combining with clinical parameters, the prognostic model was optimized. The AUC achieved 0.716 for 3 year and 0.699 for 5 year.ConclusionA series of TME-related prognostic genes were acquired in this research, which could reflect immune disorders within TME, and PTPRC and CD19 show the potential to be an indicator for LUAD prognosis and tumor microenvironment modulation. The prognostic model constructed base on those prognostic genes presented a high predictive ability, and may have clinical implications in the overall survival prediction of LUAD.
| S-EPMC7811787 | biostudies-literature