Identification of Metabolism-Associated Prostate Cancer Subtypes and Construction of a Prognostic Risk Model

A circadian rhythm-related lncRNA signature correlates with prognosis and tumor immune microenvironment in head and neck squamous cell carcinoma

OBJECTIVE

To investigate circadian rhythm-associated long non-coding RNA (lncRNA) signatures in predicting prognosis, metabolism, and immune infiltration in Head and Neck Squamous Cell Carcinoma (HNSC).

METHODS

HNSC samples were collected from the TCGA database. A signature was constructed using Cox regression and Least Absolute Shrinkage and Selection Operator (LASSO) methods. The immune cell infiltration was analyzed using CIBERSORT, ssGSEA, and MCPcounter. The RT-qPCR was used to detect the expression of signature lncRNAs.

RESULTS

A signature comprising 8 lncRNAs was constructed. The constructed signature demonstrated good prognostic prediction capability for HNSC. A nomogram encompassing risk score accurately predicted the long-term OS probability of HNSC. The infiltration levels of T cell, B cell and Macrophages were significantly higher in the high-risk group than in the low-risk group. Cluster analysis showed that the signature lncRNAs could classify the HNSC samples into two clusters. The RT-qPCR suggested that the expression of lncRNAs in signature was consistent with the data in TCGA.

CONCLUSION

The circadian rhythm-associated lncRNA signature has potential as a prognostic indicator for HNSC. It exhibits associations with metabolism, immune microenvironment, and drug sensitivity, thereby providing valuable insights for informing the treatment of HNSC.

Integrating single-cell RNA-sequencing and bulk RNA-sequencing data to explore the role of mitophagy-related genes in prostate cancer

Prostate cancer (PCa) is the most common malignancy of the male urinary system. Mitophagy, as a type of autophagy, can remove damaged mitochondria in cells. Mitophagy-related genes (MRGs) have been shown to play critical roles in the development of PCa. To this end, based on the comprehensive analysis of RNA-seq and scRNA-seq data of PCa samples and their controls, this paper identified PCa subtypes and constructed a prognostic model. In this paper, we downloaded scRNA-seq and RNA-seq data from Gene Expression Omnibus (GEO) and TCGA database. Based on the R package "Seurat" to process the scRNA-seq data, a total of five cell types were identified. Each cell population was scored based on the R package "AUCell" and using the intersection genes between MRGs and each cell population. The B cell population was then identified as a high-scoring cell population. Differentially expressed genes in RNA-seq data were identified based on the R package "limma" and intersected with previously intersected genes. Then, based on univariate Cox regression analysis and Lasso-Cox regression analysis, the prognostic genes were screened, and the risk model was constructed (composed of ADH5, CAT, BCAT2, DCXR, OGT, and FUS). The model is validated on internal and external test sets. Independent prognostic analysis identified age, N stage, and risk score as independent prognostic factors. This paper's risk models and prognostic genes can provide a reference for developing novel therapeutic targets for PCa.

Elevated ADH5 expression suggested better prognosis in kidney renal clear cell carcinoma (KIRC) and related to immunity through single-cell and bulk RNA-sequencing

BACKGROUND

Despite the rapid advances in modern medical technology, kidney renal clear cell carcinoma (KIRC) remains a challenging clinical problem in urology. Researchers urgently search for useful markers to break through the therapeutic conundrum due to its high lethality. Therefore, the study explores the value of ADH5 on overall survival (OS) and the immunology of KIRC.

METHODS

The gene expression matrix and clinical information on ADH5 in the TCGA database were validated using external databases and qRT-PCR. To confirm the correlation between ADH5 and KIRC prognosis, univariate/multivariate Cox regression analysis was used. We also explored the signaling pathways associated with ADH5 in KIRC and investigated its association with immunity.

RESULTS

The mRNA and protein levels showed an apparent downregulation of ADH5 in KIRC. Correlation analysis revealed that ADH5 was directly related to histological grade, clinical stage, and TMN stage (p < 0.05). Univariate and multivariate Cox regression analysis identified ADH5 as an independent factor affecting the prognosis of KIRC. Enrichment analysis looked into five ADH5-related signaling pathways. The results showed no correlation between ADH5 and TMB, TNB, and MSI. From an immunological perspective, ADH5 was found to be associated with the tumor microenvironment, immune cell infiltration, and immune checkpoints. Lower ADH5 expression was associated with greater responsiveness to immunotherapy. Single-cell sequencing revealed that ADH5 is highly expressed in immune cells.

CONCLUSION

ADH5 could be a promising prognostic biomarker and a potential therapeutic target for KIRC. Besides, it was found that KIRC patients with low ADH5 expression were more sensitive to immunotherapy.

Systematic review and integrated analysis of prognostic gene signatures for prostate cancer patients

Prostate cancer (PC) is one of the most common cancers in men and becoming the second leading cause of cancer fatalities. At present, the lack of effective strategies for prognosis of PC patients is still a problem to be solved. Therefore, it is significant to identify potential gene signatures for PC patients' prognosis. Here, we summarized 71 different prognostic gene signatures for PC and concluded 3 strategies for signature construction after extensive investigation. In addition, 14 genes frequently appeared in 71 different gene signatures, which enriched in mitotic and cell cycle. This review provides extensive understanding and integrated analysis of current prognostic signatures of PC, which may help researchers to construct gene signatures of PC and guide future clinical treatment.

Expressions of glucose transporter genes are diversely attenuated and significantly associated with prostate cancer progression

Prostate cancer is a health-threaten disease in men worldwide, however, lacking is the reliable biomarkers for patient management. Aberrant metabolic events including glucose metabolism are involved in prostate cancer progression. To examine the involvement of glucose metabolic pathways in prostate cancer, we analyzed the expression profiles of glucose transporter family genes using multiple RNA-seq datasets. Our results showed that three SLC2A family genes (SLC2A4/5/9) were significantly downregulated in primary prostate cancers compared to their benign compartments. These down-regulated expressions were inversely correlated with their gene promoter methylation and genome abnormalities. Among these three SLC2A genes, only SLC2A4 showed a significantly reverse correlation with all clinicopathological parameters, including TNM stage, disease relapse, Gleason score, disease-specific survival, and progression-free interval. In addition, the expression levels of these three genes were strongly correlated with anti-cancer immune cell filtration in primary prostate cancers. In a group of patients with early-onset prostate cancers, SLC2A4 also showed a strong negative correlation with multiple clinicopathological parameters, such as tumor mutation burden, biochemical relapse, pre-surgical PSA levels, and Gleason score but a positive correlation with progression-free interval after surgery. In metastatic castration-resistant prostate cancers (CRPC), SLC2A9 gene expression but not SLC2A4 or SLC2A5 genes showed a significant correlation with androgen receptor (AR) activity score and neuroendocrinal (NE) activity score. Meanwhile, SLC2A2/9/13 expression was significantly elevated in CRPC tumors with neuroendocrinal features compared to those without NE features. On the other hand, SLC2A10 and SlC2A12 gene expression were significantly reduced in NEPC tumors compared to CRPC tumors. Consistently, SLC2A10/12 expression levels were significantly reduced in castrated animals carrying the LuCaP35 xenograft models. Survival outcome analysis revealed that SLC2A4 expression in primary tumors is a favorable prognostic factor and SLC2A6 is a worse prognostic factor for disease-specific survival and progression-free survival in prostate cancer patients. In conclusion, our results suggest that SLC2A4/6 expressions are strong prognostic factors for prostate cancer progression and survival. The significance of SLC2A2/9/13 over-expression during NEPC progression needs more investigation.

Prognostic Significance of Amino Acid Metabolism-Related Genes in Prostate Cancer Retrieved by Machine Learning

Prostate cancer is among the leading cancers according to both incidence and mortality. Due to the high molecular, morphological and clinical heterogeneity, the course of prostate cancer ranges from slow growth that usually does not require immediate therapeutic intervention to aggressive and fatal disease that spreads quickly. However, currently available biomarkers cannot precisely predict the course of a disease, and novel strategies are needed to guide prostate cancer management. Amino acids serve numerous roles in cancers, among which are energy production, building block reservoirs, maintenance of redox homeostasis, epigenetic regulation, immune system modulation and resistance to therapy. In this article, by using The Cancer Genome Atlas (TCGA) data, we found that the expression of amino acid metabolism-related genes is highly aberrant in prostate cancer, which holds potential to be exploited in biomarker design or in treatment strategies. This change in expression is especially evident for catabolism genes and transporters from the solute carrier family. Furthermore, by using recursive partitioning, we confirmed that the Gleason score is strongly prognostic for progression-free survival. However, the expression of the genes SERINC3 (phosphatidylserine and sphingolipids generation) and CSAD (hypotaurine generation) can refine prognosis for high and low Gleason scores, respectively. Therefore, our results hold potential for novel prostate cancer progression biomarkers.

An Investigation of the Prognostic Role of Genes Related to Lipid Metabolism in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) has become a prevalent malignancy, and its incidence and mortality rate are increasing worldwide. Accumulating evidence has indicated that lipid metabolism-related genes (LMRGs) are involved in the occurrence and development of HNSCC. This study investigated the latent association of lipid metabolism with HNSCC and established a prognostic signature based on LMRGs. A prognostic risk model composed of eight differentially expressed LMRGs (PHYH, CYP4F8, INMT, ELOVL6, PLPP3, BCHE, TPTE, and STAR) was constructed through The Cancer Genome Atlas database. Then, ELOVL6 expression was validated in oral squamous cell carcinoma (OSCC), which is a common type of HNSCC, by immunohistochemical analysis. ELOVL6 expression in the OSCC II/III group was significantly higher than that in the other three groups (normal, dysplasia, and OSCC I), and OSCC patients with high ELOVL6 expression had poorer survival than those with low ELOVL6 expression. In summary, the LMRG-based prognostic feature had prognostic predictive capacity. ELOVL6 may be a potential prognostic factor for HNSCC patients.

Development and Validation of a Novel Circadian Rhythm-Related Signature to Predict the Prognosis of the Patients with Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is one of the most important causes of cancer-related deaths and remains a major public health challenge worldwide. Considering the extensive heterogeneity of HCC, more accurate prognostic models are imperative. The circadian genes regulate the daily oscillations of key biological processes, such as nutrient metabolism in the liver. Circadian rhythm disruption has recently been recognized as an independent risk factor for cancer. In this study, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) were compared and 248 differentially expressed genes (DEGs) of the circadian rhythm were identified. HCC was classified into two subtypes based on these DEGs. The prognostic value of each circadian rhythm-associated gene (CRG) for survival was assessed by constructing a multigene signature from TCGA cohort. A 6-gene signature was created by applying the least absolute shrinkage and selection operator (LASSO) Cox regression method, and all patients in TCGA cohort were divided into high- and low-risk groups according to their risk scores. The survival rate of patients with HCC in the low-risk group was significantly higher than that in the high-risk group (p < 0.001). The patients with HCC in the Gene Expression Omnibus (GEO) cohort were also divided into two risk subgroups using the risk score of TCGA cohort, and the overall survival time (OS) was prolonged in the low-risk group (p = 0.012). Based on the clinical characteristics, the risk score was an independent predictor of OS in the patients with HCC. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that multiple metabolic pathways, cell cycle, etc., were enhanced in the high-risk group. Using the metabolic pathway single-sample gene set enrichment analysis (ssGSEA), it was found that the metabolic pathways in the high- and low-risk groups between TCGA and GEO cohorts were altered essentially in the same way. In conclusion, the circadian genes play an important role in HCC metabolic rearrangements and can be further used to predict the prognosis the patients with HCC.

Metabolic-related gene signature model forecasts biochemical relapse in primary prostate cancer

Metabolism plays an important role in the pathogenesis of prostate cancer (PCa). Hence, we explored candidate metabolic-related genes attributed to biochemical relapse (BCR) of PCa. Gene expression profile and clinical parameters were downloaded from GSE70769 as a "training set". Using univariate Cox and LASSO-COX regression models, risk scores (RSs) were constructed. Kaplan-Meier (K-M) survival and time-dependent receiver operating characteristic (t-ROC) curves were employed. Univariate and multivariate Cox models were utilized to validate prognostic factors for biochemical relapse-free survival (BCRFS). Nomogram was plotted to facilitate clinical application. The dataset obtained from GSE70768 served as "validation set". RSs were constructed by using 7 metabolic-related genes. RSs could significantly predict 1, 3, 5-year BCRFS (AUCs for training set: 0.810-0.836; AUC for validation set: 0.673-0.827). Nomograms could effectively predicted BCRFS (training set: C-index=0.831; validation set: C-index=0.737). RSs model is an independent prognostic factor for BCR, holding greater predictive value than traditional clinicopathological parameters. Clinical Relevance- We built the prognostic nomogram based on metabolic-related gene signatures and clinicopathological features. The nomogram might further optimize biochemical relapse risk stratification for prostate cancer patients with crucial accuracy.

Predicting biochemical-recurrence-free survival using a three-metabolic-gene risk score model in prostate cancer patients

BACKGROUND

Biochemical recurrence (BCR) after initial treatment, such as radical prostatectomy, is the most frequently adopted prognostic factor for patients who suffer from prostate cancer (PCa). In this study, we aimed to construct a prognostic model consisting of gene expression profiles to predict BCR-free survival.

METHODS

We analyzed 70 metabolic pathways in 152 normal prostate samples and 494 PCa samples from the UCSC Xena dataset (training set) via gene set enrichment analysis (GSEA) to select BCR-related genes and constructed a BCR-related gene risk score (RS) model. We tested the power of our model using Kaplan-Meier (K-M) plots and receiver operator characteristic (ROC) curves. We performed univariate and multivariate analyses of RS using other clinicopathological features and established a nomogram model, which has stronger prediction ability. We used GSE70770 and DFKZ 2018 datasets to validate the results. Finally, we performed differential expression and quantitative real-time polymerase chain reaction analyses of the UCSC data for further verification of the findings.

RESULTS

A total of 194 core enriched genes were obtained through GSEA, among which 16 BCR-related genes were selected and a three-gene RS model based on the expression levels of CA14, LRAT, and MGAT5B was constructed. The outcomes of the K-M plots and ROC curves verified the accuracy of the RS model. We identified the Gleason score, pathologic T stage, and RS model as independent predictors through univariate and multivariate Cox analyses and constructed a nomogram model that presented better predictability than the RS model. The outcomes of the validation set were consistent with those of the training set. Finally, the results of differential expression analyses support the effectiveness of our model.

CONCLUSION

We constructed an RS model based on metabolic genes that could predict the prognosis of PCa patients. The model can be easily used in clinical applications and provide important insights into future research on the underlying mechanism of PCa.

Metabolic Phenotyping in Prostate Cancer Using Multi-Omics Approaches

Prostate cancer (PCa), one of the most frequently diagnosed cancers among men worldwide, is characterized by a diverse biological heterogeneity. It is well known that PCa cells rewire their cellular metabolism to meet the higher demands required for survival, proliferation, and invasion. In this context, a deeper understanding of metabolic reprogramming, an emerging hallmark of cancer, could provide novel opportunities for cancer diagnosis, prognosis, and treatment. In this setting, multi-omics data integration approaches, including genomics, epigenomics, transcriptomics, proteomics, lipidomics, and metabolomics, could offer unprecedented opportunities for uncovering the molecular changes underlying metabolic rewiring in complex diseases, such as PCa. Recent studies, focused on the integrated analysis of multi-omics data derived from PCa patients, have in fact revealed new insights into specific metabolic reprogramming events and vulnerabilities that have the potential to better guide therapy and improve outcomes for patients. This review aims to provide an up-to-date summary of multi-omics studies focused on the characterization of the metabolomic phenotype of PCa, as well as an in-depth analysis of the correlation between changes identified in the multi-omics studies and the metabolic profile of PCa tumors.

DNA Methylation Modification Regulator-Mediated Molecular Clusters and Tumor Metabolic Characterization in Prostate Cancer

Background

An increasing number of studies have indicated a close link between DNA methylation and tumor metabolism. However, the overall influence of DNA methylation on tumor metabolic characteristics in prostate cancer (PCa) remains unclear.

Methods

We first explored the subtypes of DNA methylation modification regulators and tumor metabolic features of 1,205 PCa samples using clustering analysis and gene set variation analysis based on the mRNA levels of DNA methylation modification regulators. A DNA methylation-related score (DMS) was calculated using principal component analysis and the DNA methylation modification-related gene signatures to quantify DNA methylation characteristics. We then performed a meta-analysis to identify the hazard ratio of DMS in the six cohorts. In addition, a nomogram was drawn using univariate and multivariate Cox analyses based on the DMS and clinical variables. Finally, a drug sensitivity analysis of the DMS was performed based on the genomics of drug sensitivity in cancer datasets.

Results

Three PCa clusters showing different DNA methylation modification patterns and tumor metabolic features were identified. A DMS system was established to quantify the characteristics of DNA methylation modification. PCa samples showed a differential metabolic landscape between the high and low DMS groups. The prognostic value of the DMS and nomogram was independently validated in multiple cohorts. A high DMS was associated with increases in the tumor mutation burden, copy number variation, and microsatellite instability; high tumor heterogeneity; and poor prognosis. Finally, DMS was closely related to different types of antitumor treatment.

Conclusion

Improving the understanding of tumor metabolism by characterizing DNA methylation modification patterns and using the DMS may help clinicians predict prognosis and aid in more personalized antitumor therapy strategies for PCa.

Metabolic characterization and metabolism-score of tumor to predict the prognosis in prostate cancer

Tumor metabolism patterns have been reported to be associated with the prognosis of many cancers. However, the metabolic mechanisms underlying prostate cancer (PCa) remain unknown. This study aimed to explore the metabolic characteristics of PCa. First, we downloaded mRNA expression data and clinical information of PCa samples from multiple databases and quantified the metabolic pathway activity level using single-sample gene set enrichment analysis (ssGSEA). Through unsupervised clustering and principal component analyses, we explored metabolic characteristics and constructed a metabolic score for PCa. Then, we independently validated the prognostic value of our metabolic score and the nomogram based on the metabolic score in multiple databases. Next, we found the metabolic score to be closely related to the tumor microenvironment and DNA mutation using multi-omics data and ssGSEA. Finally, we found different features of drug sensitivity in PCa patients in the high/low metabolic score groups. In total, 1232 samples were analyzed in the present study. Overall, an improved understanding of tumor metabolism through the characterization of metabolic clusters and metabolic score may help clinicians predict prognosis and aid the development of more personalized anti-tumor therapeutic strategies for PCa.

A Robust Prognostic Gene Signature Based on eRNAs-Driven Genes in Prostate Cancer

Prostate cancer (PCa) is the second most common malignancy in men, but its exact pathogenetic mechanisms remain unclear. This study explores the effect of enhancer RNAs (eRNAs) in PCa. Firstly, we screened eRNAs and eRNA -driven genes from The Cancer Genome Atlas (TCGA) database, which are related to the disease-free survival (DFS) of PCa patients;. screening methods included bootstrapping, Kaplan-Meier (KM) survival analysis, and Pearson correlation analysis. Then, a risk score model was established using multivariate Cox analysis, and the results were validated in three independent cohorts. Finally, we explored the function of eRNA-driven genes through enrichment analysis and analyzed drug sensitivity on datasets from the Genomics of Drug Sensitivity in Cancer database. We constructed and validated a robust prognostic gene signature involving three eRNA-driven genes namely MAPK15, ZNF467, and MC1R. Moreover, we evaluated the function of eRNA-driven genes associated with tumor microenvironment (TME) and tumor mutational burden (TMB), and identified remarkable differences in drug sensitivity between high- and low-risk groups. This study identified a prognostic gene signature, which provides new insights into the role of eRNAs and eRNA-driven genes while assisting clinicians to determine the prognosis and appropriate treatment options for patients with PCa.

Artificial Intelligence in Bulk and Single-Cell RNA-Sequencing Data to Foster Precision Oncology

Artificial intelligence, or the discipline of developing computational algorithms able to perform tasks that requires human intelligence, offers the opportunity to improve our idea and delivery of precision medicine. Here, we provide an overview of artificial intelligence approaches for the analysis of large-scale RNA-sequencing datasets in cancer. We present the major solutions to disentangle inter- and intra-tumor heterogeneity of transcriptome profiles for an effective improvement of patient management. We outline the contributions of learning algorithms to the needs of cancer genomics, from identifying rare cancer subtypes to personalizing therapeutic treatments.

Metabolic Classification and Intervention Opportunities for Tumor Energy Dysfunction

A comprehensive view of cell metabolism provides a new vision of cancer, conceptualized as tissue with cellular-altered metabolism and energetic dysfunction, which can shed light on pathophysiological mechanisms. Cancer is now considered a heterogeneous ecosystem, formed by tumor cells and the microenvironment, which is molecularly, phenotypically, and metabolically reprogrammable. A wealth of evidence confirms metabolic reprogramming activity as the minimum common denominator of cancer, grouping together a wide variety of aberrations that can affect any of the different metabolic pathways involved in cell physiology. This forms the basis for a new proposed classification of cancer according to the altered metabolic pathway(s) and degree of energy dysfunction. Enhanced understanding of the metabolic reprogramming pathways of fatty acids, amino acids, carbohydrates, hypoxia, and acidosis can bring about new therapeutic intervention possibilities from a metabolic perspective of cancer.