Novel molecular signatures and potential therapeutics in renal cell carcinomas: Insights from a comparative analysis of subtypes

Genome-scale metabolic models in translational medicine: the current status and potential of machine learning in improving the effectiveness of the models

The genome-scale metabolic model (GEM) has emerged as one of the leading modeling approaches for systems-level metabolic studies and has been widely explored for a broad range of organisms and applications. Owing to the development of genome sequencing technologies and available biochemical data, it is possible to reconstruct GEMs for model and non-model microorganisms as well as for multicellular organisms such as humans and animal models. GEMs will evolve in parallel with the availability of biological data, new mathematical modeling techniques and the development of automated GEM reconstruction tools. The use of high-quality, context-specific GEMs, a subset of the original GEM in which inactive reactions are removed while maintaining metabolic functions in the extracted model, for model organisms along with machine learning (ML) techniques could increase their applications and effectiveness in translational research in the near future. Here, we briefly review the current state of GEMs, discuss the potential contributions of ML approaches for more efficient and frequent application of these models in translational research, and explore the extension of GEMs to integrative cellular models.

Personalized Driver Gene Prediction Using Graph Convolutional Networks with Conditional Random Fields

Cancer is a complex and evolutionary disease mainly driven by the accumulation of genetic variations in genes. Identifying cancer driver genes is important. However, most related studies have focused on the population level. Cancer is a disease with high heterogeneity. Thus, the discovery of driver genes at the individual level is becoming more valuable but is a great challenge. Although there have been some computational methods proposed to tackle this challenge, few can cover all patient samples well, and there is still room for performance improvement. In this study, to identify individual-level driver genes more efficiently, we propose the PDGCN method. PDGCN integrates multiple types of data features, including mutation, expression, methylation, copy number data, and system-level gene features, along with network structural features extracted using Node2vec in order to construct a sample-gene interaction network. Prediction is performed using a graphical convolutional neural network model with a conditional random field layer, which is able to better combine the network structural features with biological attribute features. Experiments on the ACC (Adrenocortical Cancer) and KICH (Kidney Chromophobe) datasets from TCGA (The Cancer Genome Atlas) demonstrated that the method performs better compared to other similar methods. It can identify not only frequently mutated driver genes, but also rare candidate driver genes and novel biomarker genes. The results of the survival and enrichment analyses of these detected genes demonstrate that the method can identify important driver genes at the individual level.

A Differential Transcriptional Regulome Approach to Unpack Cancer Biology: Insights on Renal Cell Carcinoma Subtypes

Cancer research calls for new approaches that account for the regulatory complexities of biology. We present, in this study, the differential transcriptional regulome (DIFFREG) approach for the identification and prioritization of key transcriptional regulators and apply it to the case of renal cell carcinoma (RCC) biology. Of note, RCC has a poor prognosis and the biomarker and drug discovery studies to date have tended to focus on gene expression independent from mutations and/or post-translational modifications. DIFFREG focuses on the differential regulation between transcription factors (TFs) and their target genes rather than differential gene expression and integrates transcriptome profiling with the human transcriptional regulatory network to analyze differential gene regulation between healthy and RCC cases. In this study, RNA-seq tissue samples (n = 1020) from the Cancer Genome Atlas (TCGA), including healthy and tumor subjects, were integrated with a comprehensive human TF-gene interactome dataset (1122603 interactions between 1289 TFs and 25177 genes). Comparative analysis of DIFFREG profiles, consisting of perturbed TF-gene interactions, from three common subtypes (clear cell RCC, papillary RCC and chromophobe RCC) revealed subtype-specific alterations, supporting the hypothesis that these signatures in the transcriptional regulome profiles may be considered potential biomarkers that may play an important role in elucidating the molecular mechanisms of RCC development and translating knowledge about the genetic basis of RCC into the clinic. In addition, these indicators may help oncologists make the best decisions for diagnosis and prognosis management.

Cross-Talks between RKIP and YY1 through a Multilevel Bioinformatics Pan-Cancer Analysis

Recent studies suggest that PEBP1 (also known as RKIP) and YY1, despite having distinct molecular functions, may interact and mutually influence each other's activity. They exhibit reciprocal control over each other's expression through regulatory loops, prompting the hypothesis that their interplay could be pivotal in cancer advancement and resistance to drugs. To delve into this interplay's functional characteristics, we conducted a comprehensive analysis using bioinformatics tools across a range of cancers. Our results confirm the association between elevated YY1 mRNA levels and varying survival outcomes in diverse tumors. Furthermore, we observed differing degrees of inhibitory or activating effects of these two genes in apoptosis, cell cycle, DNA damage, and other cancer pathways, along with correlations between their mRNA expression and immune infiltration. Additionally, YY1/PEBP1 expression and methylation displayed connections with genomic alterations across different cancer types. Notably, we uncovered links between the two genes and different indicators of immunosuppression, such as immune checkpoint blockade response and T-cell dysfunction/exclusion levels, across different patient groups. Overall, our findings underscore the significant role of the interplay between YY1 and PEBP1 in cancer progression, influencing genomic changes, tumor immunity, or the tumor microenvironment. Additionally, these two gene products appear to impact the sensitivity of anticancer drugs, opening new avenues for cancer therapy.

The novel m6A writer METTL5 as prognostic biomarker probably associating with the regulation of immune microenvironment in kidney cancer

Nowadays, among all urinary system cancers, the mortality of kidney cancer (KC) has risen to the first, and the incidence has been keeping on the third. Many recent studies have demonstrated that m6A modification regulated by the methyltransferases (writers) is closely related to the tumorigenesis of multiple cancers. In our previous study, we found that the methyltransferase METTL5 had a stronger association with the hazard ratio of KC more than most tumors, indicating its special function in carcinogenesis of KC. Until now, the expression, functions and mechanism of METTL5 in KC are still unclear. In this study, we analyzed the mRNA expression of METTL5 using the data sets from public databases, and revealed that the METTL5 expression was significantly up-regulated in tumor tissues of kidney renal clear cell carcinoma (KIRC) and kidney renal papillary cell carcinoma (KIRP) compared to normal tissues. Also, the METTL5 expression was correlated with the tumor stage and grade, indicating the potential involvement of METTL5 in tumor progression. Additionally, the higher expression of METTL5 predicted poorer prognosis of KIRC and KIRP patients. Subsequently, we revealed that the functions of METTL5 in KIRC might be related to immune modulation, because its co-expressed gene were enriched in immune-relevant pathways including Th17 cell differentiation, Th1 and Th2 cell differentiation, and phosphatidylinositol 3-kinase activity. Next, we disclosed that the METTL5 expression was correlated to the microenvironment score and immune score of KIRC and KIRP, and associated with the infiltration ratios of 25 types of immune cells. Besides, we demonstrated a wide difference of the METTL5's effect on the survival of patients with high and low immune infiltration, further suggesting METTL5 might affect tumor development via modulating the immune microenvironment. The findings of our study provide a novel potential prognostic biomarker and immune drug target for KC.

Potential Application of Pyroptosis in Kidney Renal Clear Cell Carcinoma Immunotherapy and Targeted Therapy

Renal cell carcinoma (RCC) is a type of cancer with an increasing rate of morbidity and mortality and is a serious threat to human health. The treatment of RCC, especially kidney renal clear cell carcinoma (KIRC), has always been the focus of clinical treatment. Using The Cancer Genome Atlas (TCGA) database as a starting point, we explored the feasibility of applying the pyroptosis mechanism to KIRC treatment by searching for cancer markers associated with pyroptosis and cancer treatment signatures. The obtained samples were clustered using unsupervised clustering analysis to define the different KIRC subtypes with different pyroptosis expression levels. Based on this, a gene expression analysis was performed to explore the carcinogenic mechanism that is markedly related to pyroptosis. The Genomics of Drug Sensitivity in Cancer database and single sample gene set enrichment analysis (ssGSEA) algorithm were used to analyze the different treatment methods of the current prominent KIRC to determine whether pyroptosis plays a role. Finally, LASSO regression was used to screen for related genes and construct a model to predict patient prognosis. The expression levels of GSDME, CASP3, CASP4, CASP5, CHMP3, and CHMP4C were incorporated into the model construction. After verification, the prediction accuracy of the 3-, 5-, 7- and 10 years survival rates of our prognostic model were 0.66, 0.701, 0.719, and 0.728, respectively. Through the above analysis, we demonstrated the feasibility of pyroptosis in the clinical treatment of KIRC and provided novel ideas and suggestions for the clinical treatment of KIRC.

Multiomics Data Integration Identifies New Molecular Signatures for Abdominal Aortic Aneurysm and Aortic Occlusive Disease: Implications for Early Diagnosis, Prognosis, and Therapeutic Targets

Cardiovascular disease (CVD) is the leading cause of death among adults in developed countries. Among CVDs, abdominal aortic aneurysm (AAA) and aortic occlusive disease (AOD) are of great public health importance because of the high mortality rate in the elderly population. Despite significant molecular insights into AAA and AOD, the molecular mechanisms of these diseases remain unclear, and the current lack of robust diagnostic and prognostic biomarkers requires novel approaches to biomarker discovery and molecular targeting. In this study, we performed a comparative analysis of genome-wide expression data from patients with large AAA (n = 29), small AAA (n = 20), AOD (n = 9), and controls (n = 10). Specifically, we identified the differentially expressed genes and associated molecular pathways and biological processes (BPs) in each disease. Using a systems science approach, these data were linked to comprehensive human biological networks (i.e., protein-protein interaction, transcriptional regulatory, and metabolic networks) to identify molecular signatures of the salient mechanisms of AAA and AOD. Significant alterations in lipid metabolism and valine, leucine, and isoleucine metabolism, as well as neurodegenerative diseases and sex differences in the pathogenesis of AAA and AOD were identified. In the presence of aneurysm, size-dependent changes in lipid metabolism were observed. In addition, molecules and signaling pathways related to immunity, inflammation, infectious disease, and oxidative phosphorylation were identified in common. The results of the comparative and integrative analyzes revealed important clues to disease mechanisms and reporter molecules at various levels that warrant future development as potential prognostic biomarkers and putative therapeutic targets.

Systems biomarkers for papillary thyroid cancer prognosis and treatment through multi-omics networks

The identification of biomolecules associated with papillary thyroid cancer (PTC) has upmost importance for the elucidation of the disease mechanism and the development of effective diagnostic and treatment strategies. Despite particular findings in this regard, a holistic analysis encompassing molecular data from different biological levels has been lacking. In the present study, a meta-analysis of four transcriptome datasets was performed to identify gene expression signatures in PTC, and reporter molecules were determined by mapping gene expression data onto three major cellular networks, i.e., transcriptional regulatory, protein-protein interaction, and metabolic networks. We identified 282 common genes that were differentially expressed in all PTC datasets. In addition, six proteins (FYN, JUN, LYN, PML, SIN3A, and RARA), two Erb-B2 receptors (ERBB2 and ERBB4), two cyclin-dependent receptors (CDK1 and CDK2), and three histone deacetylase receptors (HDAC1, HDAC2, and HDAC3) came into prominence as proteomic signatures in addition to several metabolites including lactaldehyde and proline at the metabolome level. Significant associations with calcium and MAPK signaling pathways and transcriptional and post-transcriptional activities of 12 TFs and 110 miRNAs were also observed at the regulatory level. Among them, six miRNAs (miR-30b-3p, miR-15b-5p, let-7a-5p, miR-130b-3p, miR-424-5p, and miR-193b-3p) were associated with PTC for the first time in the literature, and the expression levels of miR-30b-3p, miR-15b-5p, and let-7a-5p were found to be predictive of disease prognosis. Drug repositioning and molecular docking simulations revealed that 5 drugs (prochlorperazine, meclizine, rottlerin, cephaeline, and tretinoin) may be useful in the treatment of PTC. Consequently, we report here biomolecule candidates that may be considered as prognostic biomarkers or potential therapeutic targets for further experimental and clinical trials for PTC.

Integrated Analysis to Identify a Redox-Related Prognostic Signature for Clear Cell Renal Cell Carcinoma

The imbalance of the redox system has been shown to be closely related to the occurrence and progression of many cancers. However, the biological function and clinical significance of redox-related genes (RRGs) in clear cell renal cell carcinoma (ccRCC) are unclear. In our current study, we downloaded transcriptome data from The Cancer Genome Atlas (TCGA) database of ccRCC patients and identified the differential expression of RRGs in tumor and normal kidney tissues. Then, we identified a total of 344 differentially expressed RRGs, including 234 upregulated and 110 downregulated RRGs. Fourteen prognosis-related RRGs (ADAM8, CGN, EIF4EBP1, FOXM1, G6PC, HAMP, HTR2C, ITIH4, LTB4R, MMP3, PLG, PRKCG, SAA1, and VWF) were selected out, and a prognosis-related signature was constructed based on these RRGs. Survival analysis showed that overall survival was lower in the high-risk group than in the low-risk group. The area under the receiver operating characteristic curve of the risk score signature was 0.728 at three years and 0.759 at five years in the TCGA cohort and 0.804 at three years and 0.829 at five years in the E-MTAB-1980 cohort, showing good predictive performance. In addition, we explored the regulatory relationships of these RRGs with upstream miRNA, their biological functions and molecular mechanisms, and their relationship with immune cell infiltration. We also established a nomogram based on these prognostic RRGs and performed internal and external validation in the TCGA and E-MTAB-1980 cohorts, respectively, showing an accurate prediction of ccRCC prognosis. Moreover, a stratified analysis showed a significant correlation between the prognostic signature and ccRCC progression.

Differential Interactome Proposes Subtype-Specific Biomarkers and Potential Therapeutics in Renal Cell Carcinomas

Although many studies have been conducted on single gene therapies in cancer patients, the reality is that tumor arises from different coordinating protein groups. Unveiling perturbations in protein interactome related to the tumor formation may contribute to the development of effective diagnosis, treatment strategies, and prognosis. In this study, considering the clinical and transcriptome data of three Renal Cell Carcinoma (RCC) subtypes (ccRCC, pRCC, and chRCC) retrieved from The Cancer Genome Atlas (TCGA) and the human protein interactome, the differential protein–protein interactions were identified in each RCC subtype. The approach enabled the identification of differentially interacting proteins (DIPs) indicating prominent changes in their interaction patterns during tumor formation. Further, diagnostic and prognostic performances were generated by taking into account DIP clusters which are specific to the relevant subtypes. Furthermore, considering the mesenchymal epithelial transition (MET) receptor tyrosine kinase (PDB ID: 3DKF) as a potential drug target specific to pRCC, twenty-one lead compounds were identified through virtual screening of ZINC molecules. In this study, we presented remarkable findings in terms of early diagnosis, prognosis, and effective treatment strategies, that deserve further experimental and clinical efforts.