Identification of a Hypoxia-Related Gene Model for Predicting the Prognosis and Formulating the Treatment Strategies in Kidney Renal Clear Cell Carcinoma

Deciphering the prognostic and therapeutic significance of BAG1 and BAG2 for predicting distinct survival outcome and effects on liposarcoma

Liposarcoma (LPS) is the second most common kind of soft tissue sarcoma, and a heterogeneous malignant tumor derived from adipose tissue. Up to now, the prognostic value of BAG1 or BAG2 in LPS has not been defined yet. Expression profiling data of LPS patients were collected from TCGA and GEO database. Survival curves were plotted to verify the outcome differences of patients based on BAG1 or BAG2 expression. Univariate and multivariate Cox regression models were used to analyze the prognostic ability of BAG1 or BAG2. Chaperone's regulators BAG1 and BAG2 were identified as prognostic biomarkers for LPS patients, which exhibited distinct expression patterns and survival outcome prediction performances. Patients with high BAG2 expression and/or low BAG1 expression had worse prognosis. Enrichment analysis showed that BAG1 was involved in negative regulation of TGF-β signaling. Low expression of BAG1 was associated with high abundance of regulatory T cells (Tregs). The 2-gene signature model further confirmed the improved risk assessment performance of BAG1 and BAG2: high risk patients displayed poor prognosis. BAG1 and BAG2 are supposed to be potential prognostic biomarkers for LPS and have impacts on liposarcomagenesis and immune infiltration in distinctive manners, which may function as potential therapy targets (BAG1 agonists/BAG2 inhibitors) for LPS.

AURKB promotes immunogenicity and immune infiltration in clear cell renal cell carcinoma

BACKGROUND

Chromatin regulators (CRs) are capable of causing epigenetic alterations, which are significant features of cancer. However, the function of CRs in controlling Clear Cell Renal Cell Carcinoma (ccRCC) is not well understood. This research aims to discover a CRs prognostic signature in ccRCC and to elucidate the roles of CRs-related genes in tumor microenvironment (TME).

METHODS

Expression profiles and relevant clinical annotations were retrieved from the Cancer Genome Atlas (TCGA) and UCSC Xena platform for progression-free survival (PFS) data. The R package "limma" was used to identify differentially expressed CRs. A predictive model based on five CRs was developed using LASSO-Cox analysis. The model's predictive power and applicability were validated using K-M curves, ROC curves, nomograms, comparisons with other models, stratified survival analyses, and validation with the ICGC cohort. GO and GSEA analyses were performed to investigate mechanisms differentiating low and high riskScore groups. Immunogenicity was assessed using Tumor Mutational Burden (TMB), immune cell infiltrations were inferred, and immunotherapy was evaluated using immunophenogram analysis and the expression patterns of human leukocyte antigen (HLA) and checkpoint genes. Differentially expressed CRs (DECRs) between low and high riskScore groups were identified using log2|FC|> 1 and FDR < 0.05. AURKB, one of the high-risk DECRs and a component of our prognostic model, was selected for further analysis.

RESULTS

We constructed a 5 CRs signature, which demonstrated a strong capacity to predict survival and greater applicability in ccRCC. Elevated immunogenicity and immune infiltration in the high riskScore group were associated with poor prognosis. Immunotherapy was more effective in the high riskScore group, and certain chemotherapy medications, including cisplatin, docetaxel, bleomycin, and axitinib, had lower IC50 values. Our research shows that AURKB is critical for the immunogenicity and immune infiltration of the high riskScore group.

CONCLUSION

Our study produced a reliable prognostic prediction model using only 5 CRs. We found that AURKB promotes immunogenicity and immune infiltration. This research provides crucial support for the development of prognostic biomarkers and treatment strategies for ccRCC.

Cuproptosis gene-related, neural network-based prognosis prediction and drug-target prediction for KIRC

BACKGROUND

Kidney renal clear cell carcinoma (KIRC), as a common case in renal cell carcinoma (RCC), has the risk of postoperative recurrence, thus its prognosis is poor and its prognostic markers are usually based on imaging methods, which have the problem of low specificity. In addition, cuproptosis, as a novel mode of cell death, has been used as a biomarker to predict disease in many cancers in recent years, which also provides an important basis for prognostic prediction in KIRC. For postoperative patients with KIRC, an important means of preventing disease recurrence is pharmacological treatment, and thus matching the appropriate drug to the specific patient's target is also particularly important. With the development of neural networks, their predictive performance in the field of medical big data has surpassed that of traditional methods, and this also applies to the field of prognosis prediction and drug-target prediction.

OBJECTIVE

The purpose of this study is to screen for cuproptosis genes related to the prognosis of KIRC and to establish a deep neural network (DNN) model for patient risk prediction, while also developing a personalized nomogram model for predicting patient survival. In addition, sensitivity drugs for KIRC were screened, and a graph neural network (GNN) model was established to predict the targets of the drugs, in order to discover potential drug action sites and provide new treatment ideas for KIRC.

METHODS

We used the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC) database, and DrugBank database for our study. Differentially expressed genes (DEGs) were screened using TCGA data, and then a DNN-based risk prediction model was built and validated using ICGC data. Subsequently, the differences between high- and low-risk groups were analyzed and KIRC-sensitive drugs were screened, and finally a GNN model was trained using DrugBank data to predict the relevant targets of these drugs.

RESULTS

A prognostic model was built by screening 10 significantly different cuproptosis-related genes, the model had an AUC of 0.739 on the training set (TCGA data) and an AUC of 0.707 on the validation set (ICGC data), which demonstrated a good predictive performance. Based on the prognostic model in this paper, patients were also classified into high- and low-risk groups, and functional analyses were performed. In addition, 251 drugs were screened for sensitivity, and four drugs were ultimately found to have high sensitivity, with 5-Fluorouracil having the best inhibitory effect, and subsequently their corresponding targets were also predicted by GraphSAGE, with the most prominent targets including Cytochrome P450 2D6, UDP-glucuronosyltransferase 1A, and Proto-oncogene tyrosine-protein kinase receptor Ret. Notably, the average accuracy of GraphSAGE was 0.817 ± 0.013, which was higher than that of GAT and GTN.

CONCLUSION

Our KIRC risk prediction model, constructed using 10 cuproptosis-related genes, had good independent prognostic ability. In addition, we screened four highly sensitive drugs and predicted relevant targets for these four drugs that might treat KIRC. Finally, literature research revealed that four drug-target interactions have been demonstrated in previous studies and the remaining targets are potential sites of drug action for future research.

A novel hypoxia-associated gene signature for prognosis prediction in head and neck squamous cell carcinoma

BACKGROUND

Head and neck squamous cell carcinoma (HNSCC) is the most common malignant tumor of head and neck, which seriously threatens human life and health. However, the mechanism of hypoxia-associated genes (HAGs) in HNSCC remains unelucidated. This study aims to establish a hypoxia-associated gene signature and the nomogram for predicting the prognosis of patients with HNSCC.

METHODS

Previous literature reports provided a list of HAGs. The TCGA database provided genetic and clinical information on HNSCC patients. First, a hypoxia-associated gene risk model was constructed for predicting overall survival (OS) in HNSCC patients and externally validated in four GEO datasets (GSE27020, GSE41613, GSE42743, and GSE117973). Then, immune status and metabolic pathways were analyzed. A nomogram was constructed and assessed the predictive value. Finally, experimental validation of the core genes was performed by qRT-PCR.

RESULTS

A HNSCC prognostic model was constructed based on 8 HAGs. This risk model was validated in four external datasets and exhibited high predictive value in various clinical subgroups. Significant differences in immune cell infiltration levels and metabolic pathways were found between high and low risk subgroups. The nomogram was highly accurate for predicting OS in HNSCC patients.

CONCLUSIONS

The 8 hypoxia-associated gene signature can serve as novel independent prognostic indicators in HNSCC patients. The nomogram combining the risk score and clinical stage enhanced predictive performance in predicting OS compared to the risk model and clinical characteristics alone.

A novel risk model of three SUMOylation genes based on RNA expression for potential prognosis and treatment sensitivity prediction in kidney cancer

Introduction: Kidney cancer is one of the most common and lethal urological malignancies. Discovering a biomarker that can predict prognosis and potential drug treatment sensitivity is necessary for managing patients with kidney cancer. SUMOylation is a type of posttranslational modification that could impact many tumor-related pathways through the mediation of SUMOylation substrates. In addition, enzymes that participate in the process of SUMOylation can also influence tumorigenesis and development. Methods: We analyzed the clinical and molecular data which were obtanied from three databases, The Cancer Genome Atlas (TCGA), the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC), and ArrayExpress. Results: Through analysis of differentially expressed RNA based on the total TCGA-KIRC cohort, it was found that 29 SUMOylation genes were abnormally expressed, of which 17 genes were upregulated and 12 genes were downregulated in kidney cancer tissues. A SUMOylation risk model was built based on the discovery TCGA cohort and then validated successfully in the validation TCGA cohort, total TCGA cohort, CPTAC cohort, and E-TMAB-1980 cohort. Furthermore, the SUMOylation risk score was analyzed as an independent risk factor in all five cohorts, and a nomogram was constructed. Tumor tissues in different SUMOylation risk groups showed different immune statuses and varying sensitivity to the targeted drug treatment. Discussion: In conclusion, we examined the RNA expression status of SUMOylation genes in kidney cancer tissues and developed and validated a prognostic model for predicting kidney cancer outcomes using three databases and five cohorts. Furthermore, the SUMOylation model can serve as a biomarker for selecting appropriate therapeutic drugs for kidney cancer patients based on their RNA expression.

A comprehensive analysis of the FOX family for predicting kidney renal clear cell carcinoma prognosis and the oncogenic role of FOXG1

Previous studies have confirmed that the forkhead box (FOX) superfamily of transcription factors regulates tumor progression and metastasis in multiple cancer. The purpose of this study was to develop a model based on FOX family genes for predicting kidney renal clear cell carcinom (KIRC) prognosis. We downloaded the transcriptional profiles and clinical data of KIRC patients from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) datasets. To build a new prognosis model, we screened prognosis-related FOX family genes using Lasso regression and Multivariate Cox regression analyses. Receiver operating characteristic (ROC) curves were used to evaluate model performance. Additionally, a prognostic nomogram was developed using clinical information and selected genes to improve the accuracy of prognostic prediction. We also investigated whether prognosis-related FOX family genes are related to the immune response in KIRC. Finally, we validated the oncogenic role of FOXG1 in KIRC using an in vitro tumor function assay. Six prognosis-related FOX family genes were screened: FOXO1, FOXM1, FOXK2, FOXG1, FOXA1, and FOXD1. The ROC curves indicated that our model was capable of making accurate predictions for 1-, 3-, and 5-year overall survival (OS). The nomogram further improved the accuracy of prognostic predictions. In addition, compared to those in patients with low-risk scores, high-risk scores predicted a decreased level of immune cell infiltration and a lower immune response rate. Moreover, the results of in vitro studies confirmed that FOXG1 supports the proliferation and invasion of KIRC.

Clinical significance and immune landscape of cuproptosis-related lncRNAs in kidney renal clear cell carcinoma: a bioinformatical analysis

Background

Kidney renal clear cell carcinoma (KIRC) is considered an immunogenic tumor. Cuproptosis is a newly identified copper-induced regulated cell death that relies on mitochondria respiration. Long noncoding RNAs (lncRNAs) have emerged as significant players in tumorigenesis and metastasis. However, there is a huge knowledge gap on the prognostic role of cuproptosis-related lncRNAs in KIRC. And, the clinical value of them is still unknown. Here, we aimed to develop a cuproptosis-related lncRNA prognostic signature in KIRC.

Methods

The messenger RNA (mRNA)/lncRNA expression profiles and the clinical information including age, gender, tumor stage, grade, and overall survival (OS) were acquired from The Cancer Genome Atlas (TCGA) database. The included KIRC samples were further randomly assigned into training (n=258) or testing (n=257) data sets. We performed Pearson correlation analysis to identify the cuproptosis-related lncRNAs and then constructed the prognostic signature using Cox regression analysis and LASSO algorithm. Subsequently, Kaplan-Meier survival analysis, a nomogram, and receiver operating characteristic (ROC) curve were performed to assess the predictive performance of the signature. Moreover, the immune characteristics and drug sensitivity related to the signature were also explored.

Results

The signature comprised 7 cuproptosis-related lncRNAs. The patients with a low-risk score had superior OS compared with those with a high-risk score. The survival rates of the high- and low-risk groups were 44.96% and 83.72% (P<0.001). The area under the curve (AUC) value for 1-, 3-, 5-year survival rate reached 0.814, 0.762 and 0.825, respectively. In addition, a nomogram was also generated; the AUC was 0.785 for risk score, higher than that for age (0.593), gender (0.489), grade (0.679), and stage (0.721). The high-risk group had more enriched immune- and tumor-related genes. Patients with low-risk scores were more sensitive to immunotherapy and the small molecular drugs GSK1904529A, tipifarnib, BX-912, FR-180204, and GSK1070916. Meanwhile, the high-risk group tended to be more sensitive to pyrimethamine, MS-275, and CGP-60474.

Conclusions

Collectively, we constructed a cuproptosis-related lncRNA prognostic signature with a higher predictive accuracy compared to multiple clinicopathological parameters, which may provide vital guidance for therapeutic strategies in KIRC. Combination of more prognostic biomarkers may further improve the accuracy.

Comprehensive analysis of a homeobox family gene signature in clear cell renal cell carcinoma with regard to prognosis and immune significance

The homeobox (HOX) family genes have been linked to multiple types of tumors, while their effect on malignant behaviors of clear cell renal cell carcinoma (ccRCC) and clinical significance remains largely unknown. Here, we comprehensively analyzed the expression profiles and prognostic value of HOX genes in ccRCC using datasets from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) databases. We developed a prognostic signature comprising eight HOX genes (HOXB1, HOXA7, HOXB5, HOXD8, HOXD9, HOXB9, HOXA9, and HOXA11) for overall survival prediction in ccRCC and it allowed patients to be subdivided into high- and low-risk groups. Kaplan-Meier survival analysis in all the internal and external cohorts revealed significant difference in clinical outcome of patients in different risk groups, indicating the satisfactory predictive power of the signature. Additionally, we constructed a prognostic nomogram by integrating signature-derived risk score and clinical factors such as gender, age, T and M status, which might be helpful for clinical decision-making and designing tailored management schedules. Immunological analysis revealed that the regulatory T cells (Tregs) infiltrated differently between the two subgroups in both TCGA and ICGC cohorts. ssGSEA method showed that the enrichment scores for mast cells were significantly lower in high-risk group compared with the low-risk group, which was consistent in both TCGA and ICGC cohorts. As for the related immune function, the enrichment scores of APC co-inhibition, para-inflammation, and type II IFN response were consistently lower in high-risk group in both cohorts. Of the eight HOX genes, the mRNA and protein levels of HOXD8 were downregulated in ccRCC than that in normal tissues, and decreased expression of HOXD8 was associated with increased tumor grade and stage, and lymph node metastasis. Survival analysis revealed that lower expression of HOXD8 predicted worse overall survival in ccRCC. In conclusion, our HOX gene-based signature was a favorable indicator to predict the prognosis of ccRCC cases and associated with immune cell infiltration. HOXD8 might be a tumor suppressor gene in ccRCC and a potential predictor of tumor progression.

Telomere-related gene risk model for prognosis and drug treatment efficiency prediction in kidney cancer

Kidney cancer is one of the most common urological cancers worldwide, and kidney renal clear cell cancer (KIRC) is the major histologic subtype. Our previous study found that von-Hippel Lindau (VHL) gene mutation, the dominant reason for sporadic KIRC and hereditary kidney cancer-VHL syndrome, could affect VHL disease-related cancers development by inducing telomere shortening. However, the prognosis role of telomere-related genes in kidney cancer has not been well discussed. In this study, we obtained the telomere-related genes (TRGs) from TelNet. We obtained the clinical information and TRGs expression status of kidney cancer patients in The Cancer Genome Atlas (TCGA) database, The International Cancer Genome Consortium (ICGC) database, and the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database. Totally 353 TRGs were differential between tumor and normal tissues in the TCGA-KIRC dataset. The total TCGA cohort was divided into discovery and validation TCGA cohorts and then using univariate cox regression, lasso regression, and multivariate cox regression method to conduct data analysis sequentially, ten TRGs (ISG15, RFC2, TRIM15, NEK6, PRKCQ, ATP1A1, ELOVL3, TUBB2B, PLCL1, NR1H3) risk model had been constructed finally. The kidney patients in the high TRGs risk group represented a worse outcome in the discovery TCGA cohort (p<0.001), and the result was validated by these four cohorts (validation TCGA cohort, total TCGA cohort, ICGC cohort, and CPTAC cohort). In addition, the TRGs risk score is an independent risk factor for kidney cancer in all these five cohorts. And the high TRGs risk group correlated with worse immune subtypes and higher tumor mutation burden in cancer tissues. In addition, the high TRGs risk group might benefit from receiving immune checkpoint inhibitors and targeted therapy agents. Moreover, the proteins NEK6, RF2, and ISG15 were upregulated in tumors both at the RNA and protein levels, while PLCL1 and PRKCQ were downregulated. The other five genes may display the contrary expression status at the RNA and protein levels. In conclusion, we have constructed a telomere-related genes risk model for predicting the outcomes of kidney cancer patients, and the model may be helpful in selecting treatment agents for kidney cancer patients.