A new prognostic risk model based on autophagy-related genes in kidney renal clear cell carcinoma

Identification of molecular characteristics of FUT8 and alteration of core fucosylation in kidney renal clear cell cancer

BACKGROUND

Kidney renal clear cell cancer (KIRC) is a type of urological cancer that occurs worldwide. Core fucosylation (CF), as the most common post-translational modification, is involved in the tumorigenesis.

METHODS

The alterations of CF-related genes were summarized in pan-cancer. The "ConsensusClusterPlus" package was utilized to identify two CF-related KIRC subtypes. The "ssgsea" function was chosen to estimate the CF score, signaling pathways and cell deaths. Multiple algorithms were applied to assess immune responses. The "oncoPredict" was utilized to estimate the drug sensitivity. The IHC and subgroup analysis was performed to reveal the molecular features of FUT8. Single-cell RNA sequencing (scRNA-seq) data were scrutinized to evaluate the CF state.

RESULTS

In pan-cancer, there was a noticeable alteration in the expression of CF-related genes. In KIRC, two CF-related subtypes (i.e., C1, C2) were obtained. In comparison to C2, C1 exhibited a higher CF score and correlated with poorer overall survival. Additionally, the TME of C2 demonstrated increased activity in neutrophils, macrophages, myeloid dendritic cells, and B cells, alongside a higher presence of silent mast cells, NK cells, and endothelial cells. Compared to normal samples, higher expression of FUT8 is observed in KIRC. The mutation of SETD2 was more frequent in low-FUT8 samples while the mutation of DNAH9 was more frequent in high-FUT8 samples. scRNA-seq analyses revealed that the CF score was predominantly higher in endothelial cells and fibroblast cells.

CONCLUSIONS

Two CF-related subtypes with distinct prognosis and TME were identified in KIRC. FUT8 exhibited elevated expression in KIRC samples.

Construction and verification of a novel prognostic risk model for kidney renal clear cell carcinoma based on immunity-related genes

Background: Currently, there are no useful biomarkers or prognostic risk markers for the diagnosis of kidney renal clear cell carcinoma (KIRC), although recent research has shown that both, the onset and progression of KIRC, are substantially influenced by immune-associated genes (IAGs). Objective: This work aims to create and verify the prognostic value of an immune risk score signature (IRSS) based on IAGs for KIRC using bioinformatics and public databases. Methods: Differentially expressed genes (DEGs) related to the immune systems (IAGs) in KIRC tissues were identified from The Cancer Genome Atlas (TCGA) databases. The DEGs between the tumor and normal tissues were identified using gene ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses. Furthermore, a prognostic IRSS model was constructed and its prognostic and predictive performance was analyzed using survival analyses and nomograms. Kidney renal papillary cell carcinoma (KIRP) sets were utilized to further validate this model. Results: Six independent immunity-related genes (PAEP, PI3, SAA2, SAA1, IL20RB, and IFI30) correlated with prognosis were identified and used to construct an IRSS model. According to the Kaplan-Meier curve, patients in the high-risk group had significantly poorer prognoses than those of patients in the low-risk group in both, the verification set (p <0.049; HR = 1.84; 95% CI = 1.02-3.32) and the training set (p < 0.001; HR = 3.12, 95% CI = 2.23-4.37). The numbers of regulatory T cells (Tregs) were significantly positively correlated with the six immunity-related genes identified, with correlation coefficients were 0.385, 0.415, 0.399, 0.451, 0.485, and 0.333, respectively (p <0.001). Conclusion: This work investigated the association between immune infiltration, immunity-related gene expression, and severity of KIRC to construct and verify a prognostic risk model for KIRC and KIRP.

Comprehensive characterisation of immunogenic cell death in melanoma revealing the association with prognosis and tumor immune microenvironment

Increasing evidence has highlighted the critical functions of immunogenic cell death (ICD) within many tumors. However, the therapeutic possibilities and mechanism of utilizing ICD in melanoma are still not well investigated. Melanoma samples involved in our study were acquired from The Cancer Genome Atlas (TCGA) and the Gene Expression Omnibus (GEO) databases. First, pan-cancer analysis of ICD systematically revealed its expression characteristics, prognostic values, mutation information, methylation level, pathway regulation relationship in multiple human cancers. The non-negative matrix factorization clustering was utilized to separate the TCGA-melanoma samples into two subtypes (i.e. C1 and C2) with different prognosis and immune microenvironment based on the expression traits of ICD. Then, LASSO-Cox regression analysis was utilized to determine an ICD-dependent risk signature (ICDRS) based on the differentially expressed genes (DEGs) between the two subtypes. Principal component analysis and t-distributed stochastic neighbor embedding analysis of ICDRS showed that high- and low-risk subpopulations could be clearly distinguished. Survival analysis and ROC curves in the training, internal validation, and external validation cohorts highlighted the accurate prognosis evaluation of ICDRS. The obvious discrepancies of immune microenvironment between the different risk populations might be responsible for the different prognoses of patients with melanoma. These findings revealed the close association of ICD with prognosis and tumor immune microenvironment. More importantly, ICDRS-based immunotherapy response and targeted drug prediction might be beneficial to different risk subpopulations of patients with melanoma. The innotative ICDRS could function as a marker to determine the prognosis and tumor immune microenvironment in melanoma. This will aid in patient classification for individualized melanoma treatment.

Cuprotosis-Related Genes: Predicting Prognosis and Immunotherapy Sensitivity in Pancreatic Cancer Patients

Based on TCGA, GTEx, and TIMER databases and various bioinformatics analysis methods, the potential biological roles of cuprotosis-related genes in pancreatic cancer were deeply explored, and a predictive model for pancreatic cancer patients was constructed. We downloaded the RNA-Seq data and clinicopathological and predictive data of 179 pancreatic cancer tissues and 332 adjacent normal tissues from TCGA and GTEx databases. The differential expression of cuprotosis-related genes in pancreatic cancer tissue and adjacent normal tissue was analyzed, and the LASSO regression algorithm was used to construct a prediction model and verify the validity of the model prediction. Based on the LASSO regression algorithm, a predictive model composed of three genes LIPT1, LIAS, and DLAT was screened. The corresponding survival curves showed that the constructed prediction model could significantly distinguish the prognosis of pancreatic cancer patients, and the prognosis of patients in the high-risk group was worse (P = 0.00557). The ROC curve showed that the area under the curve of the predictive model for predicting the 4-, 5-, and 6-year survival rates in pancreatic cancer was 0.816, 0.836, and 0.956, respectively. The AUC value of this risk model was significantly higher than 0.7, which could more accurately predict the prognosis of pancreatic cancer patients. This study determined a risk-scoring model of cuprotosis-related genes, which can provide an essential basis for judging the prognosis of pancreatic cancer patients.

Prognostic Signature and Therapeutic Value Based on Membrane Lipid Biosynthesis-Related Genes in Breast Cancer

There is a need to improve diagnostic and therapeutic approaches to enhance the prognosis of breast cancer, the most common malignancy worldwide. Membrane lipid biosynthesis is a hot biological pathway in current cancer research. It is unclear whether membrane lipid biosynthesis is involved in the prognosis of BRCA. With LASSO regression, a 14-gene prediction model was constructed using data from the TCGA-BRCA cohort. The prediction model includes GPAA1, PIGF, ST3GAL1, ST6GALNAC4, PLPP2, ELOVL1, HACD1, SGPP1, PRKD2, VAPB, CERS2, SGMS2, ALDH3B2, and HACD3. BRCA patients from the TCGA-BRCA cohort were divided into two risk subgroups based on the model. Kaplan–Meier survival curves showed that patients with lower risk scores had significantly improved overall survival (P=2.49e − 09). In addition, risk score, age, stage, and TNM classification were used to predict mortality in BRCA patients. In addition, the 14 genes in the risk model were analyzed for gene variation, methylation level, drug sensitivity, and immune cell infiltration, and the miRNA-mRNA network was constructed. Afterward, the THPA website then analyzed the protein expression of 14 of these risk model genes in normal and pathological BRCA tissues. In conclusion, the membrane lipid biosynthesis-related risk model and nomogram can be used to predict BRCA clinical prognosis.

The Role of ERBB Signaling Pathway-Related Genes in Kidney Renal Clear Cell Carcinoma and Establishing a Prognostic Risk Assessment Model for Patients

Objective: We aimed to investigate the potential role of ERBB signaling pathway–related genes in kidney renal clear cell carcinoma (KIRC) and establish a new predictive risk model using various bioinformatics methods.

Methods: We downloaded the KIRC dataset and clinicopathological information from The Cancer Genome Atlas database. Univariate Cox analysis was used to identify essential genes significantly associated with KIRC progression. Next, we used the STRING website to construct a protein–protein interaction network of ERBB signaling pathway–related molecules. We then used the least the absolute shrinkage and selection operator (LASSO) regression analysis to build a predictive risk model for KIRC patients. Next, we used multiple bioinformatics methods to analyze the copy number variation, single-nucleotide variation, and overall survival of these risk model genes in pan-cancer. At last, we used the Genomics of Drug Sensitivity in Cancer to investigate the correlation between the mRNA expression of genes associated with this risk model gene and drug sensitivity.

Results: Through the LASSO regression analysis, we constructed a novel KIRC prognosis–related risk model using 12 genes: SHC1, GAB1, SOS2, SRC, AKT3, EREG, EIF4EBP1, ERBB3, MAPK3, transforming growth factor-alpha, CDKN1A, and PIK3CD. Based on this risk model, the overall survival rate of KIRC patients in the low-risk group was significantly higher than that in the high-risk group (p = 1.221 × 10⁻¹⁵). Furthermore, this risk model was associated with cancer metastasis, tumor size, node, stage, grade, sex, and fustat in KIRC patients. The receiver operating characteristic curve results showed that the model had better prediction accuracy. Multivariate Cox regression analysis showed that the model’s risk score was an independent risk factor for KIRC. The Human Protein Atlas database was used to validate the protein expression of risk model–associated molecules in tumors and adjacent normal tissues. The validation results were consistent with our previous findings.

Conclusions: We successfully established a prognostic-related risk model for KIRC, which will provide clinicians with a helpful reference for future disease diagnosis and treatment.

Development of a Novel Sphingolipid Signaling Pathway-Related Risk Assessment Model to Predict Prognosis in Kidney Renal Clear Cell Carcinoma

This study aimed to explore underlying mechanisms by which sphingolipid-related genes play a role in kidney renal clear cell carcinoma (KIRC) and construct a new prognosis-related risk model. We used a variety of bioinformatics methods and databases to complete our exploration. Based on the TCGA database, we used multiple R-based extension packages for data transformation, processing, and statistical analyses. First, on analyzing the CNV, SNV, and mRNA expression of 29 sphingolipid-related genes in various types of cancers, we found that the vast majority were protective in KIRC. Subsequently, we performed cluster analysis of patients with KIRC using sphingolipid-related genes and successfully classified them into the following three clusters with significant prognostic differences: Cluster 1, Cluster 2, and Cluster 3. We performed differential analyses of transcription factor activity, drug sensitivity, immune cell infiltration, and classical oncogenes to elucidate the unique roles of sphingolipid-related genes in cancer, especially KIRC, and provide a reference for clinical treatment. After analyzing the risk rates of sphingolipid-related genes in KIRC, we successfully established a risk model composed of seven genes using LASSO regression analysis, including SPHK1, CERS5, PLPP1, SGMS1, SGMS2, SERINC1, and KDSR. Previous studies have suggested that these genes play important biological roles in sphingolipid metabolism. ROC curve analysis results showed that the risk model provided good prediction accuracy. Based on this risk model, we successfully classified patients with KIRC into high- and low-risk groups with significant prognostic differences. In addition, we performed correlation analyses combined with clinicopathological data and found a significant correlation between the risk model and patient’s M, T, stage, grade, and fustat. Finally, we developed a nomogram that predicted the 5-, 7-, and 10-year survival in patients with KIRC. The model we constructed had strong predictive ability. In conclusion, we believe that this study provides valuable data and clues for future studies on sphingolipid-related genes in KIRC.