Project description:This study was conducted to identify genes that are differentially expressed in paracancerous tissue and to determine the potential predictive value of selected gene panel. Gene transcriptome data of bladder tissue was downloaded from UCSC Xena browser and NCBI GEO repository, including GTEx (the Genotype-Tissue Expression project) data, TCGA (The Cancer Genome Atlas) data, and GEO (Gene Expression Omnibus) data. Differentially Expressed Genes (DEGs) analysis was performed to identify tumor-DEGs candidate genes, using the intersection of tumor-paracancerous DEGs genes and paracancerous-normal DEGs genes. The survival-related genes were screened by Kaplan-Meier (KM) survival analysis and univariable Cox regression with the cutoff criteria of KM < 0.05 and cox p-value < 0.05. The risk model was developed using Lasso regression. The clinical data were analyzed by univariate and multivariate Cox regression analysis. Gene Ontology (GO) and KEGG enrichment analysis were performed in the DEGs genes between the high-risk and low-risk subgroups. We identified six survival-related genes, EMP1, TPM1, NRP2, FGFR1, CAVIN1, and LATS2, found in the DEG analyses of both, tumor-paracancerous and paracancerous-normal differentially expressed data sets. Then, the patients were classified into two clusters, which can be distinguished by specific clinical characteristics. A three-gene risk prediction model (EMP1, FGFR1, and CAVIN1) was constructed in patients within cluster 1. The model was applied to categorize cluster 1 patients into high-risk and low-risk subgroups. The prognostic risk score was considered as an independent prognostic factor. The six identified survival-related genes can be used in molecular characterization of a specific subtype of bladder cancer. This subtype had distinct clinical features of T (topography), N (lymph node), stage, grade, and survival status, compared to the other subtype of bladder cancer. Among the six identified survival-related genes, three-genes, EMP1, FGFR1, and CAVIN1, were identified as potential independent prognostic markers for the specific bladder cancer subtype with clinical features described.
Project description:In recent years, genes associated with immunogenic cell death (ICD)-related genes have garnered significant interest as potential targets for immunotherapy. As a frontier in cancer treatment, immunotherapy has notably enhanced the therapeutic outcomes for cancer patients. However, since only a subset of patients benefits from this treatment approach, there is an imperative need for biomarker research to enhance patient sensitivity to immunotherapy. Expression of ICD-related genes and clinical patient data were sourced from The Cancer Genome Atlas (TCGA) database. Utilizing univariate Cox regression analysis, we constructed a signature for predicting the overall survival of colon adenocarcinoma (COAD) patients. A genomic feature analysis was performed, incorporating tumor mutation burden (TMB) and copy number variation (CNV). The immunological characteristics were analyzed via the ssGSEA and GSEA algorithms, with the resulting data visualized using R software (version 4.2.1). According to the univariate regression analysis for COAD, AIM2 emerged as the gene most significantly associated with overall survival among the 32 ICD-related genes in the TCGA dataset. Patients were divided into two groups based on high or low AIM2 expression, and genomic differences between the groups were explored. Patients expressing high levels of AIM2 had a higher TMB and a lower CNV. In addition, these patients had elevated immune checkpoint, immune cell, and immune function scores, thus indicating increased sensitivity to immunotherapy. TIDE analysis further confirmed that these patients were likely to respond more effectively to immunotherapy. Subclass mapping analysis corroborated our findings, demonstrating that patients with high AIM2 expression responded more positively to immunotherapy. Additionally, our study found that the suppression of AIM2 could significantly enhance the proliferation, invasion, and migration capabilities of colon cancer cells. In this research, we identified a novel prognostic signature suggesting that patients with higher AIM2 expression levels are more likely to respond favorably to immunotherapy.
Project description:Bladder cancer (BLCA) is the most common malignancy whose early diagnosis can ensure a better prognosis. However, the predictive accuracy of commonly used predictors, including patients' general condition, histological grade, and pathological stage, is insufficient to identify the patients who need invasive treatment. Autophagy is regarded as a vital factor in maintaining mitochondrial function and energy homeostasis in cancer cells. Whether autophagy-related genes (ARGs) can predict the prognosis of BLCA patients deserves to be investigated. Based on BLCA data retrieved from the Cancer Genome Atlas and ARGs list obtained from the Human Autophagy Database website, we identified prognosis-related differentially expressed ARGs (PDEARGs) through Wilcox text and constructed a PDEARGs-based prognostic model through multivariate Cox regression analysis. The predictive accuracy, independent forecasting capability, and the correlation between present model and clinical variables or tumor microenvironment were evaluated through R software. Enrichment analysis of PDEARGs was performed to explore the underlying mechanism, and a systematic prognostic signature with nomogram was constructed by integrating clinical variables and the aforementioned PDEARGs-based model. We found that the risk score generated by PDEARGs-based model could effectively reflect deteriorated clinical variables and tumor-promoting microenvironment. Additionally, several immune-related gene ontology terms were significantly enriched by PDEARGs, which might provide insights for present model and propose potential therapeutic targets for BLCA patients. Finally, a systematic prognostic signature with promoted clinical utility and predictive accuracy was constructed to assist clinician decision. PDEARGs are valuable prognostic predictors and potential therapeutic targets for BLCA patients.
Project description:Background: Every year, nearly 170,000 people die from bladder cancer worldwide. A major problem after transurethral resection of bladder tumor is that 40-80% of the tumors recur. Ferroptosis is a type of regulatory necrosis mediated by iron-catalyzed, excessive oxidation of polyunsaturated fatty acids. Increasing the sensitivity of tumor cells to ferroptosis is a potential treatment option for cancer. Establishing a diagnostic and prognostic model based on ferroptosis-related genes may provide guidance for the precise treatment of bladder cancer. Methods: We downloaded mRNA data in Bladder Cancer from The Cancer Genome Atlas and analyzed differentially expressed genes based on and extract ferroptosis-related genes. We identified relevant pathways and annotate the functions of ferroptosis-related DEGs using Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis and Gene Ontology functions. On the website of Search Tool for Retrieving Interacting Genes database (STRING), we downloaded the protein-protein interactions of DEGs, which were drawn by the Cytoscape software. Then the Cox regression analysis were performed so that the prognostic value of ferroptosis-related genes and survival time are combined to identify survival- and ferroptosis-related genes and establish a prognostic formula. Survival analysis and receiver operating characteristic curvevalidation were then performed. Risk curves and nomograms were generated for both groups to predict survival. Finally, RT-qPCR was applied to analyze gene expression. Results: Eight ferroptosis-related genes with prognostic value (ISCU, NFE2L2, MAFG, ZEB1, VDAC2, TXNIP, SCD, and JDP2) were identified. With clinical data, we established a prognostic model to provide promising diagnostic and prognostic information of bladder cancer based on the eight ferroptosis-related genes. RT-qPCR revealed the genes that were differentially expressed between normal and cancer tissues. Conclusion: This study found that the ferroptosis-related genes is associated with bladder cancer, which may serve as new target for the treatment of bladder cancer.
Project description:The tumor microenvironment (TME) is a complex system that plays an important role in tumor development and progression, but the current knowledge about its effect on bladder cancer (BC) is scarce. In this study, we performed a comprehensive analysis of the relationship between the TME and gene expression profiles to identify prognostic biomarkers for BC. The ESTIMATE algorithm was used to calculate immune and stromal scores of BC patients who were obtained from the Gene Expression Omnibus database. We found that the immune and stromal scores were associated with clinical characteristics and the prognosis of BC patients. Based on these scores, 104 immune-related differentially expressed genes were identified. Further, functional enrichment analysis revealed that these genes were mainly involved in the immune-related biological processes and signaling pathways. Three prognostic genes were then identified and used to establish a risk prediction model using Cox regression analyses. Kaplan-Meier survival analysis showed that the expression levels of COL1A1, COMP, and SERPINE2 significantly correlated with cancer-specific survival and overall survival of BC patients. Additionally, we validated the prognostic values of these genes using two independent cohorts from The Cancer Genome Atlas and Gene Expression Omnibus databases. Finally, the relationships between the three prognostic genes and several immune cells were evaluated using Tumor Immune Estimation Resource, indicating that the expression levels of COL1A1, COMP, and SERPINE2 correlated positively with the tumor infiltration levels of CD4+ T cells and macrophages. In conclusion, the current study comprehensively analyzed the TME and presented immune-related prognostic genes for BC, providing new insights into immunotherapeutic strategies for BC patients.
Project description:BackgroundThis study aimed to identify potential stemness-related targets in gastric cancer (GC) in order to support the development of new treatment strategies and improve patient survival.MethodsUsing the edgeR package, we identified stemness-related differentially expressed genes (DEGs) using GSE112631 and the stemness-related signaling pathways in the Gene Set Enrichment Analysis (GSEA) database. Lasso-penalized Cox regression analysis and multivariate Cox regression analysis tested by Akaike Information Criterion (AIC) were used to screen out survival genes in order to construct a prognostic model. We verified the accuracy of our prognostic model using a nomogram and receiver operating characteristic (ROC) curve analysis. Patients were divided into two groups based on the median risk score, and functional enrichment analysis was used to explore the differences between the two groups.ResultsEight genes were selected to establish a prognostic model of The Cancer Genome Atlas (TCGA) and a validation model of the GSE84437 dataset from the Genome Expression Omnibus (GEO). In both models, we found that the low risk score group had better overall survival (OS) than the high-risk score group. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways between the two risk groups were totally different.ConclusionsWe used eight stemness-related genes to build a prognostic model. The high-risk score group had a worse prognosis compared to the low-risk score group.
Project description:Bladder cancer is one of the most common genitourinary malignant cancers worldwide. Cell death processes, including apoptosis, ferroptosis, and necrosis, provide novel clinical and immunological insights promoting the management of precision medicine. Therefore, this study aimed to evaluate the transcriptomic profile of signatures in cell death pathways with significant prognostic implications in patients with bladder cancer from multiple independent cohorts (n = 1999). First, genes involved in apoptosis (n = 19), ferroptosis (n = 31), and necrosis (n = 6) were analyzed to evaluate the prognostic implications in bladder cancer. Significant genes were included to establish the cell-death index (CDI) of 36 genes that distinguished patients according to high and low risks. Survival analysis using the Kaplan-Meier curves clustered patients based on overall survival (18.8 vs. 96.7 months; hazard model [HR] = 3.12, P<00001). Cox proportional hazard model was significantly associated with a higher risk of mortality using 10 external independent cohorts in patients with CDIhigh (HR = 1.31, 95% CI: 1.04-1.62). To explore immune parameters associated with CDI, microenvironment cell-population-counter algorithms indicated increased intratumoral heterogeneity and macrophage/monocyte infiltration and CD8+ T cells in patients with CDIhigh group. Besides, the CDIhigh group showed an increased expression of the following immune checkpoints: CD276, PD-L1, CTLA-4, and T-cell exhaustion signatures. Cytokine expression analysis revealed the highest association of IL-9R, IL-17A, IL-17F, GDF7, and IFNW1 with the high-risk group. In addition, 42 patients with BCa receiving immunotherapies were enrolled from a real-world cohort, and expression patterns of three CDI hub genes (DRD5, SCL2A14, and IGF1) were detected using immunohistochemical staining. Patients with triple-negative staining of tumor tissues had significantly higher tumor-associated macrophage abundance, PD-L1 expression, predicted immunocompromised microenvironment, and prominently progressive progression (HR = 4.316, P = 0.0028). In conclusion, this study highlights the immunoevasive tumor microenvironment characterized by the higher tumor-associated macrophage infiltration with the presence of immune checkpoint and T-cell exhaustion genes in patients with BCa at CDIhigh risk who might suffer progression and be more suitable to benefit from immune checkpoint inhibitors or other immunotherapies.
Project description:Immunosenescence refers to the immune system undergoing a series of degenerative changes with advancing age and is tightly associated with the initiation and progression of cancers. However, the immunosenescence-related genes as critical biomarkers for bladder cancer (BLCA) have not been systematically analyzed. We retrieved the immunosenescence-related genes from the public database and verified their association with hallmarks of immunosenescence based on The Cancer Genome Atlas (TCGA) cohort. Through gene pairing, Lasso, and univariate Cox regression, an 8-gene pair model was constructed to evaluate the overall survival of BLCA, which was then validated in the training cohort (P < 0.001, n = 396), two external validation cohorts (P < 0.05, n = 165; P < 0.001, n = 224), and local samples (P < 0.05, n = 10). We also downloaded the clinical information and gene expression matrices of other 32 different cancers from TCGA. The established model showed significant predictive value for the prognosis in 15 cancers (P < 0.05). The risk model could also serve as a promising predictor for immunotherapeutic response, which has been verified by the TIDE algorithm (P < 0.05), IMvigor210 dataset (P < 0.01, n = 298), and other two datasets correlated with immunotherapy (P < 0.05, n = 56; P = 0.17, n = 27). The TCGA dataset, in vitro cell experiments, and pan-cancer analysis displayed that the gene signature was associated with cisplatin sensitivity (P < 0.05). Overall, we proposed a novel immunosenescence-related gene signature to predict prognosis, immunotherapeutic response, and cisplatin sensitivity of BLCA, which were validated in different independent cohorts, local samples, and pan-cancer analyses.
Project description:BackgroundBladder cancer (BC) is a commonly occurring malignant tumor of the urinary system, demonstrating high global morbidity and mortality rates. BC currently lacks widely accepted biomarkers and its predictive, preventive, and personalized medicine (PPPM) is still unsatisfactory. N6-methyladenosine (m6A) modification and non-coding RNAs (ncRNAs) have been shown to be effective prognostic and immunotherapeutic responsiveness biomarkers and contribute to PPPM for various tumors. However, their role in BC remains unclear.Methodsm6A-related ncRNAs (lncRNAs and miRNAs) were identified through a comprehensive analysis of TCGA, starBase, and m6A2Target databases. Using TCGA dataset (training set), univariate and least absolute shrinkage and selection operator (LASSO) regression analyses were performed to develop an m6A-related ncRNA-based prognostic risk model. Kaplan-Meier analysis of overall survival (OS) and receiver operating characteristic (ROC) curves were used to verify the prognostic evaluation power of the risk model in the GSE154261 dataset (testing set) from Gene Expression Omnibus (GEO). A nomogram containing independent prognostic factors was developed. Differences in BC clinical characteristics, m6A regulators, m6A-related ncRNAs, gene expression patterns, and differentially expressed genes (DEGs)-associated molecular networks between the high- and low-risk groups in TCGA dataset were also analyzed. Additionally, the potential applicability of the risk model in the prediction of immunotherapeutic responsiveness was evaluated based on the "IMvigor210CoreBiologies" data set.ResultsWe identified 183 m6A-related ncRNAs, of which 14 were related to OS. LASSO regression analysis was further used to develop a prognostic risk model that included 10 m6A-related ncRNAs (BAALC-AS1, MIR324, MIR191, MIR25, AC023509.1, AL021707.1, AC026362.1, GATA2-AS1, AC012065.2, and HCP5). The risk model showed an excellent prognostic evaluation performance in both TCGA and GSE154261 datasets, with ROC curve areas under the curve (AUC) of 0.62 and 0.83, respectively. A nomogram containing 3 independent prognostic factors (risk score, age, and clinical stage) was developed and was found to demonstrate high prognostic prediction accuracy (AUC = 0.83). Moreover, the risk model could also predict BC progression. A higher risk score indicated a higher pathological grade and clinical stage. We identified 1058 DEGs between the high- and low-risk groups in TCGA dataset; these DEGs were involved in 3 molecular network systems, i.e., cellular immune response, cell adhesion, and cellular biological metabolism. Furthermore, the expression levels of 8 m6A regulators and 12 m6A-related ncRNAs were significantly different between the two groups. Finally, this risk model could be used to predict immunotherapeutic responses.ConclusionOur study is the first to explore the potential application value of m6A-related ncRNAs in BC. The m6A-related ncRNA-based risk model demonstrated excellent performance in predicting prognosis and immunotherapeutic responsiveness. Based on this model, in addition to identifying high-risk patients early to provide them with focused attention and targeted prevention, we can also select beneficiaries of immunotherapy to deliver personalized medical services. Furthermore, the m6A-related ncRNAs could elucidate the molecular mechanisms of BC and lead to a new direction for the improvement of PPPM for BC.Supplementary informationThe online version contains supplementary material available at 10.1007/s13167-021-00259-w.