Integrated genomic analysis of proteasome alterations across 11,057 patients with 33 cancer types: clinically relevant outcomes in framework of 3P medicine

Development of an m6A subtype classifier to guide precision therapy for patients with bladder cancer

Purpose: Bladder cancer (BLCA) is a highly heterogeneous tumor. We aim to construct a classifier from the perspective of N6-methyladenosine methylation (m6A) to identify patients with different prognostic risks and treatment responsiveness for precision therapy. Methods: Data on gene expression profile, mutation, and clinical characteristics were mainly obtained from the TCGA-BLCA cohort. Unsupervised clustering was performed to construct m6A subtypes. The tumor microenvironment (TME) landscapes were explored by using ssGSEA, ESTIMATE, and MCPcounter algorithms. K-M survival curves and Cox regression analysis were used to demonstrate the significance of m6A subtypes in predicting prognosis. pRRophetic, oncoPredict, and TIDE algorithms were used to evaluate responsiveness to antitumor therapy. A classifier of m6a subtypes was finally developed based on random forest and artificial neural network (ANN). Results: The two m6A subtypes have significantly different m6A-related gene expression profiles and mutational landscapes. TME analysis showed a higher level of stromal and Inhibitory immune components in subtype B compared with subtype A. The m6A subtype is a clinically independent prognostic predictor of BLCA, subtype B has a poorer prognosis. Drug sensitivity analysis showed that subtype B has lower IC50 values and AUC values for cisplatin and docetaxel. Efficacy assessment showed significantly poorer radiotherapy efficacy and lower immunotherapy responsiveness in subtype B. We finally constructed an ANN classifier to accurately classify BLCA patients into two m6A subtypes. Conclusion: Our study developed a classifier for identifying subtypes with different m6A characteristics, and BLCA patients with different m6A subtypes have significantly different prognosis and responsiveness to antitumor therapy.

A new gene signature associated with disulfidptosis that forecasts myasthenia gravis and suggests infiltration of immune microenvironment in thymoma patients

Introduction

The mechanism of thymoma-associated myasthenia gravis (TAMG) is currently unknown, although patients with TAMG experience more severe myasthenic symptoms and have worse prognoses compared to regular thymoma patients. The objective of this research is to create a transcriptome map of TAMG using genes linked to disulfidptosis, detect possible biomarkers, and examine the disparities in the tumor immune microenvironment (TIME) among different thymoma patients. The findings will offer valuable knowledge for personalized treatment alternatives.

Methods

Thymoma samples' RNA-seq data, along with their corresponding clinical data, were acquired from the TCGA database using methods. Next, genes and disulfidptosis-related lncRNAs(DRLs) were chosen through correlation analysis. Then, a prediction model of TAMG was established by LASSO regression. Subsequent to that, an analysis of the mutation data, the tumor mutational burden (TMB), and the assessment of immune and stromal elements within the tumor microenvironment were conducted.

Results

A total of 87 patients diagnosed with thymoma were included in the study, with 29 of them having TAMG. We discovered a group of 325 lncRNAs in this sample that showed significant associations with genes related to disulfidptosis, with 25 of them displaying significantly altered expression. Moreover, utilizing LASSO regression, we constructed a predictive model incorporating 11 DRLs. The analysis revealed an area under the curve (AUC) of 0.934 (CI 0.879-0.989), a cut-off value of 0.797, along with a sensitivity of 82.8 % and specificity of 93.1 %. Furthermore, we examined the TIME in both the high-risk and low-risk groups, and observed noteworthy disparities in B cells, T cells, and APC among the two groups (p < 0.05).

Conclusion

This research offers the initial examination of genes associated with disulfidptosis and TAMG through genomic and transcriptomic analysis. Furthermore, a strong risk forecasting model was created and the significance of TIME in TAMG was also clarified. The discoveries offer significant understanding into the molecular processes of TAMG and present possible indicators for categorizing risk.

Regulator of G protein signaling-1 regulates immune infiltration and macrophage polarization in clear cell renal cell carcinoma

OBJECTIVE

To better understand how to clear cell renal cell cancer (ccRCC) is affected by the regulator of G protein signaling-1 (RGS1), its effect on immune infiltration, macrophage polarization, tumor proliferation migration, and to explore whether RGS1 may serve as a marker and therapeutic target for ccRCC.

PATIENTS AND METHODS

In this study, a total of 20 surgical specimens of patients with pathological diagnosis of ccRCC admitted to the Department of Urology of the Second Affiliated Hospital of Anhui Medical University from November 2021 to June 2022 were selected for pathological and protein testing, while the expression of RGS1 in tumors, immune infiltration, and macrophage polarization, particularly M2 macrophage linked to the development of tumor microenvironment (TME), were combined with TGCA database and GO analysis. We also further explored and studied the expression and function of RGS1 in TME, investigated how RGS1 affected tumor growth, migration, apoptosis, and other traits, and initially explored the signaling pathways and mechanisms that RGS1 may affect.

RESULTS

RGS1 was found to be expressed at higher quantities in ccRCC than in normal cells or tissues, according to bioinformatics analysis and preliminary experimental data from this work. Using the TCGA database and GO analysis to describe the expression of RGS1 in a range of tumors, it was found that ccRCC had a much higher level of RGS1 expression than other tumor types. The results of gene enrichment analysis indicated that overexpression of RGS1 may be associated with immune infiltration. The outcomes of in vitro tests revealed that RGS1 overexpression in ccRCC did not significantly alter the proliferation and migration ability of ccRCC, but RGS1 overexpression promoted apoptosis in ccRCC. By in vitro co-culture experiments, RGS1 overexpression inhibited M2 macrophage polarization and also suppressed the Jagged-1/Notch signaling pathway.

CONCLUSIONS

RGS1 is highly expressed in ccRCC, while overexpression of RGS1 may increase immune infiltration in the TME and reduce the polarization of M2 macrophages while promoting apoptosis in ccRCC.

Ubiquitinomics revealed disease- and stage-specific patterns relevant for the 3PM approach in human sigmoid colon cancers

Objective

The patients with sigmoid colorectal cancer commonly show high mortality and poor prognosis. Increasing evidence has demonstrated that the ubiquitinated proteins and ubiquitination-mediated molecular pathways influence the growth and aggressiveness of colorectal cancer. It emphasizes the scientific merits of quantitative ubiquitinomics in human sigmoid colon cancer. We hypothesize that the ubiquitinome and ubiquitination-mediated pathway networks significantly differ in sigmoid colon cancers compared to controls, which offers the promise for in-depth insight into molecular mechanisms, discovery of effective therapeutic targets, and construction of reliable biomarkers in the framework of predictive, preventive, and personalized medicine (PPPM; 3P medicine).

Methods

The first ubiquitinome analysis was performed with anti-K-ε-GG antibody beads (PTMScan ubiquitin remnant motif [K-ε-GG])-based label-free quantitative proteomics and bioinformatics to identify and quantify ubiquitination profiling between sigmoid colon cancer tissues and para-carcinoma tissues. A total of 100 human sigmoid colon cancer samples that included complete clinical information and the corresponding gene expression data were obtained from The Cancer Genome Atlas (TCGA). Ubiquitination was the main way of protein degradation; the relationships between differentially ubiquitinated proteins (DUPs) and their differently expressed genes (DEGs) and between DUPs and their differentially expressed proteins (DEPs) were analyzed between cancer tissues and control tissues. The overall survival of those DUPs was obtained with Kaplan-Meier method.

Results

A total of 1249 ubiquitinated sites within 608 DUPs were identified in human sigmoid colon cancer tissues. KEGG pathway network analysis of these DUPs revealed 35 statistically significant signaling pathways, such as salmonella infection, glycolysis/gluconeogenesis, and ferroptosis. Gene Ontology (GO) analysis of 608 DUPs revealed that protein ubiquitination was involved in 98 biological processes, 64 cellular components, 51 molecule functions, and 26 immune system processes. Protein-protein interaction (PPI) network of 608 DUPs revealed multiple high-combined scores and co-expressed DUPs. The relationship analysis between DUPs and their DEGs found 4 types of relationship models, including DUP-up (increased ubiquitination level) and DEG-up (increased gene expression), DUP-up and DEG-down (decreased gene expression), DUP-down (decreased ubiquitination level) and DEG-up, and DUP-down and DEG-down. The relationship analysis between DUPs and their DEPs found 4 types of relationship models, including DUP-up and DEP-up (increased protein expression), DUP-up and DEP-down (decreased protein expression), DUP-down and DEP-up, and DUP-down and DEP-down. Survival analysis found 46 overall survival-related DUPs in sigmoid colon cancer, and the drug sensitivity of overall survival-related DUPs were identified.

Conclusion

The study provided the first differentially ubiquitinated proteomic profiling, ubiquitination-involved signaling pathway network changes, and the relationship models between protein ubiquitination and its gene expression and between protein ubiquitination and its protein expression, in human sigmoid colon cancer. It offers the promise for deep insights into molecular mechanisms of sigmoid colon cancer, and discovery of effective therapeutic targets and biomarkers for patient stratification, predictive diagnosis, prognostic assessment, and personalized treatment in the context of 3P medicine.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-023-00328-2.

Quantitative phosphoproteomics reveals molecular pathway network alterations in human early-stage primary hepatic carcinomas: potential for 3P medical approach

Objective

Hepatic carcinoma is one of the most common types of malignant tumors in the digestive system, and its biological characteristics determine its high rate of metastasis and recurrence after radical resection, leading to a poor prognosis for patients. Increasing evidence demonstrates that phosphoproteins and phosphorylation-mediated molecular pathways influence the occurrence and development of hepatic carcinoma. It is urgent need to develop early-stage biomarkers for improving diagnosis, therapy, medical service, and prognostic assessment. We hypothesize that phosphoproteome and phosphorylation-mediated signaling pathway networks significantly differ in human early-stage primary hepatic carcinomas relative to control liver tissues, which will identify the key differentially phosphorylated proteins and phosphorylation-mediated signaling pathway network alterations in human early-stage primary hepatic carcinoma to innovate predictive diagnosis, prognostic assessment, and personalized medical services and progress beyond the state of the art in the framework of predictive, preventive, and personalized medicine (PPPM).

Methods

Tandem mass tag (TMT)-based quantitative proteomics coupled with TiO2 enrichment of phosphopeptides was used to identify phosphorylation profiling, and bioinformatics was used to analyze the pathways and biological functions of phosphorylation profiling between early-stage hepatic carcinoma tissues and tumor-adjacent normal control tissues. Furthermore, the integrative analysis with transcriptomic data from TCGA database obtained differently expressed genes (DEGs) corresponding to differentially phosphorylated proteins (DPPs) and overall survival (OS)-related DPPs.

Results

A total of 1326 phosphopeptides derived from 858 DPPs in human early-stage primary hepatic carcinoma were identified. KEGG pathway network analysis of 858 DPPs revealed 33 statistically significant signaling pathways, including spliceosome, glycolysis/gluconeogenesis, B-cell receptor signaling pathway, HIF-1 signaling pathway, and fatty acid degradation. Gene Ontology (GO) analysis of 858 DPPs revealed that protein phosphorylation was involved in 57 biological processes, 40 cellular components, and 37 molecular functions. Protein-protein interaction (PPI) network constructed multiple high-combined scores and co-expressed DPPs. Integrative analysis of transcriptomic data and DPP data identified 105 overlapped molecules (DPPs; DEGs) between hepatic carcinoma tissues and control tissues and 125 OS-related DPPs. Overlapping Venn plots showed 14 common molecules among datasets of DPPs, DEGs, and OS-related DDPs, including FTCD, NDRG2, CCT2, PECR, SLC23A2, PNPLA7, ANLN, HNRNPM, HJURP, MCM2, STMN1, TCOF1, TOP2A, and SSRP1. The drug sensitivities of OS-related DPPs were identified, including LMOD1, CAV2, UBE2E2, RAPH1, ANXA5, HDLBP, CUEDC1, APBB1IP, VCL, SRSF10, SLC23A2, EPB41L2, ESR1, PLEKHA4, SAFB2, SMARCAD1, VCAN, PSD4, RDH16, NOP56, MEF2C, BAIAP2L2, NAGS, SRSF2, FHOD3, and STMN1.

Conclusions

Identification and annotation of phosphoproteomes and phosphorylation-mediated signaling pathways in human early-stage primary hepatic carcinoma tissues provided new directions for tumor prevention and treatment, which (i) helps to enrich phosphorylation functional research and develop new biomarkers; (ii) enriches phosphorylation-mediated signaling pathways to gain a deeper understanding of the underlying mechanisms of early-stage primary hepatic carcinoma; and (iii) develops anti-tumor drugs that facilitate targeted phosphorylated sites. We recommend quantitative phosphoproteomics in early-stage primary hepatic carcinoma, which offers great promise for in-depth insight into the molecular mechanism of early-stage primary hepatic carcinoma, the discovery of effective therapeutic targets/drugs, and the construction of reliable phosphorylation-related biomarkers for patient stratification, predictive diagnosis, prognostic assessment, and personalized medical services in the framework of PPPM.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13167-023-00335-3.

Construction and comprehensive analysis of a novel prognostic signature associated with pyroptosis molecular subtypes in patients with pancreatic adenocarcinoma

Background

Treatment of cancer with pyroptosis is an emerging strategy. Molecular subtypes based on pyroptosis-related genes(PRGs) seem to be considered more conducive to individualized therapy. It is meaningful to construct a pyroptosis molecular subtypes-related prognostic signature (PMSRPS) to predict the overall survival (OS) of patients with pancreatic adenocarcinoma(PAAD) and guide treatment.

Methods

Based on the transcriptome data of 23 PRGs, consensus clustering was applied to divide the TCGA and GSE102238 combined cohort into three PRGclusters. Prognosis-related differentially expressed genes(DEGs) among PRGclusters were subjected to LASSO Cox regression analysis to determine a PMSRPS. External cohort and in vitro experiments were conducted to verify this PMSRPS. The CIBERSORT algorithm, the ESTIMATE algorithm and the Immunophenoscore (IPS) were used to analyze the infiltrating abundance of immune cells, the tumor microenvironment (TME), and the response to immunotherapy, respectively. Wilcoxon analysis was used to compare tumor mutational burden (TMB) and RNA stemness scores (RNAss) between groups. RT-qPCR and in vitro functional experiments were used for evaluating the expression and function of SFTA2.

Results

Based on three PRGclusters, 828 DEGs were obtained and a PMSRPS was subsequently constructed. In internal and external validation, patients in the high-risk group had significantly lower OS than those in the low-risk group and PMSRPS was confirmed to be an independent prognostic risk factor for patients with PAAD with good predictive performance. Immune cell infiltration abundance and TME scores indicate patients in the high-risk group have typical immunosuppressive microenvironment characteristics. Analysis of IPS suggests patients in the high-risk group responded better to novel immune checkpoint inhibitors (ICIs) than PD1/CTLA4. The high-risk group had higher TMB and RNAss. In addition, 10 potential small-molecule compounds were screened out. Finally, we found that the mRNA expression of SFTA2 gene with the highest risk coefficient in PMSRPS was significantly higher in PAAD than in paracancerous tissues, and knockdown of it significantly delayed the progression of PAAD.

Conclusions

PMSRPS can well predict the prognosis, TME and immunotherapy response of patients with PAAD, identify potential drugs, and provide treatment guidance based on individual needs.

Machine Learning Identifies Pan-Cancer Landscape of Nrf2 Oxidative Stress Response Pathway-Related Genes

BACKGROUND

Oxidative stress produced a large amount of reactive oxygen species (ROS), which played a pivotal role in balanced ability and determining cell fate. The activated Nrf2 signaling pathway that responds to the excessive ROS regulated the expressions of antiapoptotic proteins, antioxidative enzymes, drug transporters, and detoxifying factors.

METHODS

The Nrf2 signaling pathway-related genes that had a direct relationship with Nrf2, including ATF4, BACH1, CREBBP, CUL3, EIF2AK3, EP300, FOS, FOSL1, GSK3B, JUN, KEAP1, MAF, MAFF, MAFG, MAFK, MAPK1, MAPK3, MAPK7, MAPK8, MAPK9, PIK3CA, PRRT2, and RIT1, were selected to do a systematic pan-cancer analysis. The relationship of Nrf2 signaling pathway-related gene expressions with tumor mutation burden, microsatellite status, clinical characteristics, immune system, cancer stemness index, and drug sensitivity was calculated by the Spearson correlation analysis across 11,057 subjects representing 33 cancer types. The prognosis models in lung squamous carcinoma, breast cancer, and stomach cancer were constructed with the Cox multivariate regression analysis and least absolute shrinkage and selection operator (Lasso) regression.

RESULTS

Many Nrf2 signaling pathway-related genes were differently expressed between tumor and normal tissues. PIK3CA showed high mutation rate in pan-cancer. The expressions of Nrf2 signaling pathway-related genes were significantly related to tumor mutation burden, copy number variant, microsatellite instability score, survival rate, pathological stage, immune phenotype, immune score, immune cell, cancer stemness index, and drug sensitivity. The prognosis models were significantly associated with survival rate in lung squamous carcinoma, breast cancer, and stomach cancer; and the prognosis model-based riskscore was significantly associated with clinicopathological characteristics of each cancer.

CONCLUSIONS

The study provided a comprehensive pan-cancer landscape of Nrf2 pathway-related genes. Based on the same Nrf2 pathway-related genes, the different prognosis models were constructed for different types of cancers.