Establishing and validating a pathway prognostic signature in pancreatic cancer based on miRNA and mRNA sets using GSVA

Tumor hypoxia in immune infiltration and prognosis of bladder cancer

Background

Bladder cancer (BC) is the sixth most common cancer and the ninth leading cause of cancer death among men in the world. Previous studies have shown that tumor hypoxia plays an important role in the occurrence and development of BC, but the role of tumor hypoxia in the prognosis and immune infiltration of BC remains unclear. Our aim was to perform a bioinformatics analysis combined with a clinical analysis to explore the roles of hypoxia in BC.

Methods

We acquired datasets (GSE13507, GSE5287, and GSE1827) containing mRNA expression information from BC cohorts from the Gene Expression Omnibus (GEO) and measured the Hypoxia score using the Gene Set Variation Analysis (GSVA). Then we used X-tile method and log-rank test and Pearson's correlation test to analyze the relation among the Hypoxia score and the clinicopathological and immunological characteristics of BC and used stepwise Cox regression analysis to establish a Prognostic model.

Results

Hypoxia was found to be closely associated with tumor grade, pathological type, invasion, and prognosis of BC in our study. Moreover, we determined that hypoxia was closely related to the infiltration abundance of multiple immune cells through a correlation analysis, and the tumor immune cell infiltration was further found to be significantly associated with the tumor grade and tumor type of BC. Furthermore, we constructed several models based on the Hypoxia score and tumor immune infiltration with C-indexes ranging from 0.703 and 0.888, which showed good performance in predicting the prognosis of BC.

Conclusions

Our study showed that hypoxia plays an important role in the progression, prognosis, and tumor immune infiltration of BC. Our models based on hypoxia and tumor immune infiltration play a guiding role in the prognosis and treatment of BC patients.

AMIGO2 is a pivotal therapeutic target related to M2 polarization of macrophages in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a common kind of lethal cancer, with low early diagnostic rate and poor prognosis. In this study, we identified and verified the AMIGO2 with significant diagnostic and prognostic value in PDAC through LASSO regression combined with multiple machines learning methods, including RVM-RFE and Random Forest in TCGA and GEO datasets. The relevance between the expression of AMIGO2 and M2 polarization of macrophages was identified through pancancer, normal tissue, and cell lines data in TCGA, GTEx and CCLE datasets. The relevance between AMIGO2 and M2 polarization was then further identified in our local PDAC cohort. Finally, the role of AMIGO2 as cancer promoter and pivotal factor enrolled in M2 polarization was verified through siRNA transfection and M2 macrophages induction. These findings could facilitate diagnosis and treatment of PDAC. In addition, further research was deemed necessary on the deep mechanism between AMIGO2 and M2 polarization of macrophages in PDAC.

Signatures of tumor-associated macrophages correlate with treatment response in ovarian cancer patients

Ovarian cancer (OC) ranks as the second leading cause of death among gynecological cancers. Numerous studies have indicated a correlation between the tumor microenvironment (TME) and the clinical response to treatment in OC patients. Tumor-associated macrophages (TAMs), a crucial component of the TME, exert influence on invasion, metastasis, and recurrence in OC patients. To delve deeper into the role of TAMs in OC, this study conducted an extensive analysis of single-cell data from OC patients. The aim is to develop a new risk score (RS) to characterize the response to treatment in OC patients to inform clinical treatment. We first identified TAM-associated genes (TAMGs) in OC patients and examined the protein and mRNA expression levels of TAMGs by Western blot and PCR experiments. Additionally, a scoring system for TAMGs was constructed, successfully categorizing patients into high and low RS subgroups. Remarkably, significant disparities were observed in immune cell infiltration and immunotherapy response between the high and low RS subgroups. The findings revealed that patients in the high RS group had a poorer prognosis but displayed greater sensitivity to immunotherapy. Another important finding was that patients in the high RS subgroup had a higher IC50 for chemotherapeutic agents. Furthermore, further experimental investigations led to the discovery that THEMIS2 could serve as a potential target in OC patients and is associated with EMT (epithelial-mesenchymal transition). Overall, the TAMGs-based scoring system holds promise for screening patients who would benefit from therapy and provides valuable information for the clinical treatment of OC.

Analysis of m6A modulator-mediated methylation modification patterns and the tumor microenvironment in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer. In the development and progression of LUAD, epigenetic aberration plays a crucial role. However, the function of RNA N6-methyladenosine (m⁶A) modifications in the LUAD progression is unknown. The m6A regulator modification patterns in 955 LUAD samples were analyzed comprehensively. Patterns were systematically correlated with the tumor microenvironment (TME) cell-infiltration characteristics. Using principal component analysis algorithms, the m6Ascore was generated to quantify m⁶A modification patterns in individual tumors. Then, their values for predicting prognoses and therapeutic response in LUAD patients were assessed. Three distinct m⁶A modification patterns in LUAD were identified. Among them, the prognosis of m6Acluster C was the best, while the prognosis of m⁶Acluster A was the worst. Interestingly, the characterization of TME cell infiltration and biological behavior differed among the three patterns. To evaluate m⁶A modification patterns within individual tumors, an m⁶Ascore signature was constructed. The results showed that the high m⁶Ascore group was associated with a better prognosis; tumor somatic mutations and tumor microenvironment differed significantly between the high- and low- m⁶Ascore groups. Furthermore, in the cohort with anti-CTLA-4 treatment alone, patients with a high m⁶Ascore had higher ICI scores, which indicated significant therapeutic advantage and clinical benefits.

NXPH4 Used as a New Prognostic and Immunotherapeutic Marker for Muscle-Invasive Bladder Cancer

Background

One of the most common malignant tumors of the urinary system is muscle-invasive bladder cancer (MIBC). With the increased use of immunotherapy, its importance in the field of cancer is becoming abundantly evident. This study classifies MIBC according to GSVA score from the perspective of the GSEA immune gene set.

Methods

This study integrated the sequencing and clinical data of MIBC patients in TCGA and GEO databases, then scored the data using the GSVA algorithm, the CNMF algorithm was implemented to divide the subtypes of GEO and TCGA datasets, respectively, and finally screened and determined the key pathways in combination with clinical data. Simultaneously, LASSO Cox regression model was constructed based on key pathway genes to assess the model's predictive ability (ROC) and describe the immune landscape differences between high- and low-risk groups; key genes were further analyzed and verified in patient tissues.

Results

404 TCGA and 297 GEO datasets were divided into C1-3 groups (TCGA-C1:120/C2:152/C3:132; GEO- C1:112/C2:101/C3:84), of which TCGA-C2 (n = 152) subtype and GEO-C1 (n = 112) subtype had the worst prognosis. LASSO Cox regression model with ROC (train set = 0.718, test set = 0.667) could be constructed. When combined with the Cancer Immunome Atlas database, it was found that patients with high-risk scores were more sensitive to PD-1 inhibitor and PD-1 inhibitor combined with CTLA-4. NXPH4, as a key gene, plays a role in MIBC with tissue validation results show that nxph4 is highly expressed in tumor.

Conclusion

The immune gene score of MIBC data in TCGA and GEO databases was successfully evaluated using GSVA in this research. The lasso Cox expression model was successfully constructed by screening immune genes, the high-risk group had a worse prognosis and higher sensitivity to immunotherapy, PD-1 inhibitors or PD-1 combined with CTLA-4 inhibitors can be preferentially used in high-risk patients who are sensitive to immunotherapy, and NXPH4 may be a molecular target to adjust the effect of immunotherapy.

Comprehensive pan-cancer analysis identifies cellular senescence as a new therapeutic target for cancer: multi-omics analysis and single-cell sequencing validation

Although cellular senescence has long been recognized as an anti-tumor mechanism, mounting evidence suggests that in some circumstances, senescent cells promote tumor growth and malignancy spread. Therefore, research into the exact relationship between cellular senescence and tumor immunity is ongoing. We analyzed changes in the expression, copy number variation, single-nucleotide variation, methylation, and drug sensitivity of cellular senescence-related genes in 33 tumor types. The cellular senescence score was calculated using the single-sample gene-set enrichment analysis. The correlations between cellular senescence score and prognosis, tumor immune microenvironment (TIME), and expression of tumor immune-related genes were comprehensively analyzed. Single-cell transcriptome sequencing data were used to assess the activation state of cellular senescence in the tumor microenvironment (TME). The expression of cellular senescence-associated hub genes varied significantly across cancer types. In these genes, missense mutation was the major type of single nucleotide polymorphism, and heterozygous deletion and heterozygous amplification were the major types of copy number variation. Moreover, the cellular senescence pathway in tumors was sensitive to drugs such as XMD13-2, TPCA-1, methotrexate, and KIN001-102. Furthermore, the cellular senescence score was significantly higher in most cancer types, related to poor prognosis. The expression of immune checkpoint molecules such as NRP1, CD276, and CD44 was significantly correlated with the cellular senescence score. Monocyte cellular senescence was significantly higher in the TME of kidney renal clear cell carcinoma cells than in normal tissues. The findings of this study provide insights into the important role of cellular senescence in the TIME of human cancers and the effect of immunotherapy.

Identification and Validation of Three PDAC Subtypes and Individualized GSVA Immune Pathway-Related Prognostic Risk Score Formula in Pancreatic Ductal Adenocarcinoma Patients

Background

With the progress of precision medicine treatment in pancreatic ductal adenocarcinoma (PDAC), individualized cancer-related medical examination and prediction are of great importance in this high malignant tumor and tumor-immune microenvironment with changed pathways highly enrolled in the carcinogenesis of PDAC.

Methods

High-throughput data of pancreatic ductal adenocarcinoma were downloaded from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database. After batch normalization, the enrichment pathway and relevant scores were identified by the enrichment of immune-related pathway signature using gene set variation analysis (GSVA). Then, cancerous subtype in TCGA and GEO samples was defined through the NMF methods by cancertypes packages in R software, respectively. Subsequently, the significance between the characteristics of each TCGA sample and cancer type and the significant prognosis-related pathway with risk score formula is calculated through t-test and univariate Cox analysis. Next, the prognostic value of gained risk score formula and each significant prognosis-related pathway were validated in TCGA and GEO samples by survival analysis. The pivotal hub genes in the enriched significant prognosis-related pathway are identified and validated, and the TIMER database was used to identify the potential role of hub genes in the PDAC immune environment. The potential role of hub genes is promoting the transdifferentiation of cancer-associated fibroblasts.

Results

The enrichment pathway and relevant scores were identified by GSVA, and 3 subtypes of pancreatic ductal adenocarcinoma were defined in TCGA and GEO samples. The clinical stage, tumor node metastasis classification, and tumor grade are strongly relative to the subtype above in TCGA samples. A risk formula about GSVA significant pathway “GSE45365_WT_VS_IFNAR_KO_CD11B_DC_MCMV_INFECTION_DN ∗ 0.80 + HALLMARK_GLYCOLYSIS ∗ 16.8 + GSE19888_CTRL_VS_T_CELL_MEMBRANES_ACT_MAST_CELL_DN ∗ 14.4” was identified and validated in TCGA and GEO samples through survival analysis with significance. DCN, VCAN, B4GALT7, SDC1, SDC2, B3GALT6, B3GAT3, SDC3, GPC1, and XYLT2 were identified as hub genes in these GSVA significant pathways and validated in silico.

Conclusions

Three pancreatic ductal adenocarcinoma subtypes are identified, and an individualized GSVA immune pathway score-related prognostic risk score formula with 10 hub genes is identified and validated. The predicted function of the 10 upregulated hub genes in tumor-immune microenvironment was promoting the infiltration of cancer-associated fibroblasts. These findings will contribute to the precision medicine of pancreatic ductal adenocarcinoma treatment and tumor immune-related basic research.

Role of Calcium Signaling Pathway-Related Gene Regulatory Networks in Ischemic Stroke Based on Multiple WGCNA and Single-Cell Analysis

Background

This study is aimed at investigating the changes in relevant pathways and the differential expression of related gene expression after ischemic stroke (IS) at the single-cell level using multiple weighted gene coexpression network analysis (WGCNA) and single-cell analysis.

Methods

The transcriptome expression datasets of IS samples and single-cell RNA sequencing (scRNA-seq) profiles of cerebrovascular tissues were obtained by searching the Gene Expression Omnibus (GEO) database. First, gene pathway scoring was calculated via gene set variation analysis (GSVA) and was imported into multiple WGCNA to acquire key pathways and pathway-related hub genes. Furthermore, SCENIC was used to identify transcription factors (TFs) regulating these core genes using scRNA-seq data. Finally, the pseudotemporal trajectory analysis was used to analyse the role of these TFs on various cell types under hypoxic and normoxic conditions.

Results

The scores of 186 KEGG pathways were obtained via GSVA using microarray expression profiles of 40 specimens. WGCNA of the KEGG pathways revealed the two following pathways: calcium signaling pathway and neuroactive ligand-receptor interaction pathways. Subsequently, WGCNA of the gene expression matrix of the samples revealed the calcium signaling pathway-related genes (AC079305.10, BCL10, BCL2A1, BRE-AS1, DYNLL2, EREG, and PTGS2) that were identified as core genes via correlation analysis. Furthermore, SCENIC and pseudotemporal analysis revealed JUN, IRF9, ETV5, and PPARA score gene-related TFs. Jun was found to be associated with hypoxia in endothelial cells, whereas Irf9 and Etv5 were identified as astrocyte-specific TFs associated with oxygen concentration in the mouse cerebral cortex.

Conclusions

Calcium signaling pathway-related genes (AC079305.10, BCL10, BCL2A1, BRE-AS1, DYNLL2, EREG, and PTGS2) and TFs (JUN, IRF9, ETV5, and PPARA) were identified to play a key role in IS. This study provides a new perspective and basis for investigating the pathogenesis of IS and developing new therapeutic approaches.

The Role of microRNA in Pancreatic Cancer

MicroRNAs (miRNAs) are small ribonucleic acid molecules that play a key role in regulating gene expression. The increasing number of studies undertaken on the functioning of microRNAs in the tumor formation clearly indicates their important potential in oncological therapy. Pancreatic cancer is one of the deadliest cancers. The expression of miRNAs released into the bloodstream appears to be a good indicator of progression and evaluation of the aggressiveness of pancreatic cancer, as indicated by studies. The work reviewed the latest literature on the importance of miRNAs for pancreatic cancer development.

Comprehensive analysis of a chemokine- and chemokine receptor family-based signature for patients with lung adenocarcinoma

The clinical significance and comprehensive features of chemokines and their receptors in lung adenocarcinoma (LUAD) have not been clarified. We aimed to characterize the expression profiles of chemokine and chemokine receptor family members and construct a chemokine- and chemokine receptor-based prognosis signature. A total of 1511 patients with LUAD from seven independent cohorts were included in the study. The training set collected from The Cancer Genome Atlas (TCGA) database containing 468 cases. The validation was performed on the basis of six different cohorts downloaded from Gene Expression Omnibus (GEO) database. A five-chemokine- and chemokine receptor-(CXCL2, CXCL13, CCL26, CCL20, CX3CR1) based prognosis signature was constructed with TCGA dataset using LASSO Cox regression and Cox proportional hazards regression analysis. A multivariate analysis verified that this signature was an independent prognostic factor. The predictive value of this signature was further verified by other six independent cohorts and multiple clinical subtypes. We performed immune cell infiltration analysis and biological pathway analysis which provided more insight into this signature-related immune and inflammatory landscape and clarified the intrinsic molecular mechanism by which this signature could be used to predict clinical prognosis. Furthermore, we explored the close relationship between this signature and tumor mutation burden (TMB), neoantigen burden, PD-1, PD-L1, CTLA4, TIDE score, T cell-inflamed score. This signature provided a robust prognostic biomarker for LUAD and could serve as a predictor for immunotherapy response, which may be used as an important supplement to immunotherapy to achieve individualized tumor treatment by optimizing the prognostic management and immunotherapy for patients with LUAD.