Hexokinase 3 dysfunction promotes tumorigenesis and immune escape by upregulating monocyte/macrophage infiltration into the clear cell renal cell carcinoma microenvironment

Unlocking the potential: Targeting metabolic pathways in the tumor microenvironment for Cancer therapy

Cancer incidence and mortality are increasing and impacting global life expectancy. Metabolic reprogramming in the tumor microenvironment (TME) is intimately related to tumorigenesis, progression, metastasis and drug resistance. Tumor cells drive metabolic reprogramming of other cells in the TME through metabolic induction of cytokines and metabolites, and metabolic substrate competition. Consequently, this boosts tumor cell growth by providing metabolic support and facilitating immunosuppression and angiogenesis. The metabolic interplay in the TME presents potential therapeutic targets. Here, we focus on the metabolic reprogramming of four principal cell subsets in the TME: CAFs, TAMs, TILs and TECs, and their interaction with tumor cells. We also summarize medications and therapies targeting these cells' metabolic pathways, particularly in the context of immune checkpoint blockade therapy.

TCF3 as a multidimensional biomarker: oncogenicity, genomic alterations, and immune landscape in pan-cancer analysis

Transcription factor 3 (TCF3), a pivotal member of the TCF/LEF family, plays a critical role in tumorigenesis. Nonetheless, its impact on the tumor microenvironment (TME) and cancer phenotypes remains elusive. We perform an exhaustive analysis of TCF3 expression, DNA variation profiles, prognostic implications, and associations with the TME and immunological aspects. This study is based on a large-scale pan-cancer cohort, encompassing over 17,000 cancer patients from multiple independent datasets, validated by in vitro assays. Our results show that TCF3/4/7 exhibits differential expression patterns between normal and tumor tissues across pan-cancer analyses. Mutational analysis of TCF3 across diverse cancer types reveals the highest alteration rates in biliary tract cancer. Additionally, mutations and single nucleotide variants in TCF3/4/7 are found to exert varied effects on patient prognosis. Importantly, TCF3 emerges as a robust predictor of survival across all cancer cohorts and among patients receiving immune checkpoint inhibitors. Elevated TCF3 expression is correlated with more aggressive cancer subtypes, as validated by immunohistochemistry and diverse cohort data. Furthermore, TCF3 expression is positively correlated with intratumoral heterogeneity and angiogenesis. In vitro investigations demonstrate that TCF3 is involved in epithelial-mesenchymal transition, migration, invasion, and angiogenesis. These effects are likely mediated through the interaction of TCF3 with the NF-κB/MMP2 pathway, which is modulated by IL-17A in human uveal melanoma MUM2B cells. This study elucidates, for the first time, the significant associations of TCF3 with DNA variation profiles, prognostic outcomes, and the TME in multiple cancer contexts. TCF3 holds promise as a molecular marker for diagnosis and as a potential target for novel therapeutic strategies, particularly in uveal melanoma.

Hexokinase HK3-mediated O-GlcNAcylation of EP300: a key regulator of PD-L1 expression and immune evasion in ccRCC

Clear cell renal cell carcinoma (ccRCC) demonstrates enhanced glycolysis, critically contributing to tumor development. Programmed death-ligand 1 (PD-L1) aids tumor cells in evading T-cell-mediated immune surveillance. Yet, the specific mechanism by which glycolysis influences PD-L1 expression in ccRCC is not fully understood. Our research identified that the glycolysis-related gene (GRG) HK3 has a unique correlation with PD-L1 expression. HK3 has been identified as a key regulator of O-GlcNAcylation in ccRCC. O-GlcNAcylation exists on the serine 900 (Ser900) site of EP300 and can enhance its stability and oncogenic activity by preventing ubiquitination. Stably expressed EP300 works together with TFAP2A as a co-transcription factor to promote PD-L1 transcription and as an acetyltransferase to stabilize PD-L1 protein. Furthermore, ccRCC exhibits interactive dynamics with tumor-associated macrophages (TAMs). The uridine 5'-diphospho-N-acetylglucosamine (UDP-GlcNAc), which serves as a critical substrate for the O-GlcNAcylation process, facilitates TAMs polarization. In ccRCC cells, HK3 expression is influenced by IL-10 secreted by M2 TAMs. Our study elucidates that HK3-mediated O-GlcNAcylation of EP300 is involved in tumor immune evasion. This finding suggests potential strategies to enhance the efficacy of immune checkpoint blockade therapy.

PD1/PD-L1 blockade in clear cell renal cell carcinoma: mechanistic insights, clinical efficacy, and future perspectives

The advent of PD1/PD-L1 inhibitors has significantly transformed the therapeutic landscape for clear cell renal cell carcinoma (ccRCC). This review provides an in-depth analysis of the biological functions and regulatory mechanisms of PD1 and PD-L1 in ccRCC, emphasizing their role in tumor immune evasion. We comprehensively evaluate the clinical efficacy and safety profiles of PD1/PD-L1 inhibitors, such as Nivolumab and Pembrolizumab, through a critical examination of recent clinical trial data. Furthermore, we discuss the challenges posed by resistance mechanisms to these therapies and potential strategies to overcome them. We also explores the synergistic potential of combination therapies, integrating PD1/PD-L1 inhibitors with other immunotherapies, targeted therapies, and conventional modalities such as chemotherapy and radiotherapy. In addition, we examine emerging predictive biomarkers for response to PD1/PD-L1 blockade and biomarkers indicative of resistance, providing a foundation for personalized therapeutic approaches. Finally, we outline future research directions, highlighting the need for novel therapeutic strategies, deeper mechanistic insights, and the development of individualized treatment regimens. Our work summarizes the latest knowledge and progress in this field, aiming to provide a valuable reference for improving clinical efficacy and guiding future research on the application of PD1/PD-L1 inhibitors in ccRCC.

Exploring variations in glycolytic and gluconeogenic enzymes and isoforms across breast cancer cell lines and tissues

It is well established that tumour cells undergo metabolic changes to acquire biological advantage over normal cells with activation of the glycolytic pathway, a process termed "Warburg effect". Enzyme isoforms are alternative enzymatic forms with the same function but with different biochemical or epigenetic features. Moreover, isoforms may have varying impacts on different metabolic pathways. We challenge ourselves to analyse the glycolytic and gluconeogenic enzymes and isoforms in breast cancer, a complex and heterogeneous pathology, associated with high incidence and mortality rates especially among women. We analysed epithelial and tumour cell lines by RT-PCR and compared values to a publicly available database for the expression profile of normal and tumour tissues (Gepia) of enzymes and enzymatic isoforms from glycolytic and gluconeogenic pathways. Additionally, GeneMANIA was used to evaluate interactions, pathways, and attributes of each glycolytic/gluconeogenic steps. The findings reveal that the enzymes and enzymatic isoforms expressed in cell culture were somewhat different from those in breast tissue. We propose that the tumor microenvironment plays a crucial role in the expression of glycolytic and gluconeogenic enzymes and isoforms in tumour cells. Nonetheless, they not only participate in glycolytic and gluconeogenic enzymatic activities but may also influence other pathways, such as the Pentose-Phosphate-Pathway, TCA cycle, as well as other carbohydrate, lipid, and amino acid metabolism.

Identification of HK3 as a promising immunomodulatory and prognostic target in sepsis-induced acute lung injury

BACKGROUND

Sepsis is a life-threatening global disease with a significant impact on human health. Acute lung injury (ALI) has been identified as one of the primary causes of mortality in septic patients. This study aimed to identify candidate genes involved in sepsis-induced ALI through a comprehensive approach combining bioinformatics analysis and experimental validation.

METHODS

The datasets GSE65682 and GSE32707 obtained from the Gene Expression Omnibus database were merged to screen for sepsis-induced ALI related differentially expressed genes (DEGs). Functional enrichment and immune infiltration analyses were conducted on DGEs, with the construction of protein-protein interaction (PPI) networks to identify hub genes. In vitro and in vivo models of sepsis-induced ALI were used to study the expression and function of hexokinase 3 (HK3) using various techniques including Western blot, real-time PCR, immunohistochemistry, immunofluorescence, Cell Counting Kit-8, Enzyme-linked immunosorbent assay, and flow cytometry.

RESULTS

The results of bioinformatics analysis have identified HK3, MMP9, and S100A8 as hub genes with diagnostic and prognostic significance for sepsis-induced ALI. The HK3 has profound effects on sepsis-induced ALI and exhibits a correlation with immune regulation. Experimental results showed increased HK3 expression in lung tissue of septic mice, particularly in bronchial and alveolar epithelial cells. In vitro studies demonstrated upregulation of HK3 in lipopolysaccharide (LPS)-stimulated lung epithelial cells, with cytoplasmic localization around the nucleus. Interestingly, following the knockdown of HK3 expression, lung epithelial cells exhibited a significant decrease in proliferation activity and glycolytic flux, accompanied by an increase in cellular inflammatory response, oxidative stress, and cell apoptosis.

CONCLUSIONS

It was observed for the first time that HK3 plays a crucial role in the progression of sepsis-induced ALI and may be a valuable target for immunomodulation and therapy.Bioinformatics analysis identified HK3, MMP9, and S100A8 as hub genes with diagnostic and prognostic relevance in sepsis-induced ALI. Experimental findings showed increased HK3 expression in the lung tissue of septic mice, particularly in bronchial and alveolar epithelial cells. In vitro experiments demonstrated increased HK3 levels in lung epithelial cells stimulated with LPS, with cytoplasmic localization near the nucleus. Knockdown of HK3 expression resulted in decreased proliferation activity and glycolytic flux, increased inflammatory response, oxidative stress, and cell apoptosis in lung epithelial cells.

Metabolic alterations in hereditary and sporadic renal cell carcinoma

Kidney cancer is the seventh leading cause of cancer in the world, and its incidence is on the rise. Renal cell carcinoma (RCC) is the most common form and is a heterogeneous disease comprising three major subtypes that vary in their histology, clinical course and driver mutations. These subtypes include clear cell RCC, papillary RCC and chromophobe RCC. Molecular analyses of hereditary and sporadic forms of RCC have revealed that this complex and deadly disease is characterized by metabolic pathway alterations in cancer cells that lead to deregulated oxygen and nutrient sensing, as well as impaired tricarboxylic acid cycle activity. These metabolic changes facilitate tumour growth and survival. Specifically, studies of the metabolic features of RCC have led to the discovery of oncometabolites - fumarate and succinate - that can promote tumorigenesis, moonlighting functions of enzymes, and substrate auxotrophy owing to the disruption of pathways that enable the production of arginine and cholesterol. These metabolic alterations within RCC can be exploited to identify new therapeutic targets and interventions, in combination with novel approaches that minimize the systemic toxicity of metabolic inhibitors and reduce the risk of drug resistance owing to metabolic plasticity.

Unveiling the impact of tertiary lymphoid structures on immunotherapeutic responses of clear cell renal cell carcinoma

Tertiary lymphoid structures (TLS) are organized aggregates of immune cells that form under pathological conditions. However, the predictive value of TLS in clear cell renal cell carcinoma (ccRCC) for immunotherapies remains unclear. We comprehensively assessed the implications for prognosis and immunological responses of the TLS spatial and maturation heterogeneity in 655 ccRCC patients. A higher proportion of early-TLS was found in peritumoral TLS, while intratumoral TLS mainly comprised secondary follicle-like TLS (SFL-TLS), indicating markedly better survival. Notably, presence of TLS, especially intratumoral TLS and SFL-TLS, significantly correlated with better survival and objective reflection rate for ccRCC patients receiving anti-Programmed Cell Death Protein-1 (PD-1)/Programmed Cell Death-Ligand-1 (PD-L1) immunotherapies. In peritumoral TLS cluster, primary follicle-like TLS, the proportion of tumor-associated macrophages, and Treg infiltration in the peritumoral regions increased prominently, suggesting an immunosuppressive tumor microenvironment. Interestingly, spatial transcriptome annotation and multispectral fluorescence showed that an abundance of mature plasma cells within mature TLS has the capacity to produce IgA and IgG, which demonstrate significantly higher objective response rates and a superior prognosis for ccRCC patients subjected to immunotherapy. In conclusion, this study revealed the implications of TLS spatial and maturation heterogeneity on the immunological status and clinical responses, allowing the improvement of precise immunotherapies of ccRCC.

Heterogeneity in tertiary lymphoid structures predicts distinct prognosis and immune microenvironment of clear cell renal cell carcinoma

BACKGROUND

Tertiary lymphoid structures (TLS) are organized aggregates of immune cells that develop postnatally in non-lymphoid tissues and are associated with pathological conditions. TLS typically comprise B-cell follicles containing and are encompassed by T- cell zones and dendritic cells. The prognostic and predictive value of TLS in the tumor microenvironment (TME) as potential mediators of antitumor immunity have gained interest. However, the precise relationship between localization and maturation of TLS and the clinical outcome of their presence in clear cell renal cell carcinoma (ccRCC) is yet to be elucidated.

METHODS

Immunohistochemistry and multispectral fluorescence were used to evaluate the TLS heterogeneity along with TME cell-infiltrating characterizations. A thorough investigation of the prognostic implications of the TLS heterogeneity in 395 patients with ccRCC from two independent cohorts was conducted. Associations between TLS heterogeneity and immunologic activity were assessed by quantifying the immune cell infiltration.

RESULTS

Infiltrated TLS were identified in 34.2% of the ccRCC samples (N=395). These TLS were found to be tumor-proximal, tumor-distal, or both in 37.8%, 74.1%, and 11.9% of the TLS-positive cases, respectively. A higher proportion of early TLS was found in tumor-distal TLS (p=0.016), while tumor-proximal TLS primarily comprised secondary follicle-like structures (p=0.004). In the main study cohort (Fudan University Shanghai Cancer Center, N=290), Kaplan-Meier analyses revealed a significant correlation between the presence of tumor-proximal TLS and improved progression-free survival (PFS, p<0.001) and overall survival (OS, p=0.002). Conversely, the presence of tumor-distal TLS was associated with poor PFS (p=0.02) and OS (p=0.021). These findings were further validated in an external validation set of 105 patients with ccRCC. Notably, the presence of mature TLS (namely secondary follicle-like TLS, with CD23+ germinal center) was significantly associated with better clinical outcomes in patients with ccRCC. Furthermore, novel nomograms incorporating the presence of tumor-proximal TLS demonstrated remarkable predictability for the 8-year outcomes of resected ccRCC (area under the curve >0.80). Additionally, ccRCC samples with tumor-distal TLS enriched with primary follicle-like TLS exhibited higher programmed death-ligand 1 tumor-associated macrophages levels and regulatory T cells infiltration in the tumor-distal region, indicative of a suppressive TME.

CONCLUSION

This study for the first time elucidates the impact of TLS localization and maturation heterogeneities on the divergent clinical outcomes of ccRCC. The findings reveal that most TLS in ccRCC are located in the tumor-distal area and are associated with immature, immunosuppressive characterizations. Furthermore, our findings corroborate previous research demonstrating that tumor-proximal TLS were associated with favorable clinical outcomes.

Deciphering the implications of mitophagy-related signatures in clinical outcomes and microenvironment heterogeneity of clear cell renal cell carcinoma

BACKGROUND

The role of mitophagy in various cancer-associated biological processes is well recognized. Nonetheless, the comprehensive implications of mitophagy in clear cell renal cell carcinoma (ccRCC) necessitate further exploration.

METHODS

Based on the transcriptomic data encompassing 25 mitophagy-related genes (MRGs), we identified the distinct mitophage patterns in 763 ccRCC samples. Subsequently, a mitophage-related predictive signature with machine learning algorithms was constructed, designated as RiskScore, to quantify the individual mitophagy status in ccRCC patients. Employing multispectral immunofluorescence (mIF) and immunohistochemistry (IHC) staining, we detected the effect of PTEN-induced putative kinase 1 (PINK1) in the prognosis and immune microenvironment of ccRCC.

RESULTS

Our analysis initially encompassed a comprehensive assessment of the expression profiling, genomic variations, and interactions among the 25 MRGs in ccRCC. Subsequently, the consensus clustering algorithm was applied to stratify ccRCC patients into three clusters with distinct prognostic outcomes, tumor microenvironment (TME) characteristics, and underlying biological pathways. We screened eight pivotal genes (CLIC4, PTPRB, SLC16A12, ENPP5, FLRT3, HRH2, PDK4, and SCD5) to construct a mitophagy-related predictive signature, which showed excellent prognostic value for ccRCC patients. Moreover, patient subgroups divided by the RiskScore showed contrasting expression levels of immune checkpoints (ICPs), abundance of immune cells, and immunotherapy response. Additionally, a nomogram was established with robust predictive power integrating the RiskScore and clinical features. Notably, we observed that PINK1 expression markedly correlated with favorable treatment response and advanced maturation stages of tertiary lymphoid structures, which potentially shed light on enhancing anti-tumor immunity of ccRCC.

CONCLUSION

Collectively, this study initially developed a signature associated with mitophagy, which demonstrated an excellent ability to predict the clinical prognosis, TME characterization, and responsiveness to targeted therapy and immunotherapy for ccRCC patients. Of particular note is the pivotal role of PINK1 in mediating the treatment response and immune microenvironment for ccRCC patients.

PLAC8 is an innovative biomarker for immunotherapy participating in remodeling the immune microenvironment of renal clear cell carcinoma

Background

PLAC8 has been identified in the progression of various cancers by inducing tumorigenesis, immune response, chemotherapy resistance and metastasis. Nevertheless, the precise biological function of PLAC8 in renal cancer remains unknown.

Methods

We obtained the expression profile and associated clinical characteristics of patients diagnosed with clear cell renal cell carcinoma (ccRCC) from The Cancer Genome Atlas database. The biological behavior of specific cell lines was detected using Cell Counting Kit-8 (CCK-8), colony formation, and 5-ethynyl-2'-deoxyuridine (EdU) assay. A prognostic model was constructed based on PLAC8-related molecules through a machine-learning algorithm.

Results

We observed overexpression of PLAC8 in ccRCC patients. In addition, PLAC8 has been identified as being linked to unfavorable clinical characteristics and adverse prognosis outcomes. Biological enrichment analysis revealed the potential involvement of PLAC8 in cell cycle checkpoints, mitotic phase transformation, immunotherapy-predicted and reactive oxygen species (ROS) related pathways. In addition, immune analyses showed that PLAC8 was involved in remodeling the tumor microenvironment (TME) and affecting the effect of immunotherapy in ccRCC patients. In vitro experiments demonstrated a significant reduction in the proliferation, invasion and migration of renal cancer cells following the knockdown of PLAC8. Finally, LASSO logistics regression was applied to construct a prognosis model, which presented a favorable prediction ability on the prognosis of ccRCC.

Conclusion

Our results implied that PLAC8 may be a novel immunotherapy biomarker of ccRCC, which is a crucial molecule in remodeling the cancer microenvironment. PLAC8 can predict immunotherapy response and is expected to guide precise treatment.

HK3 stimulates immune cell infiltration to promote glioma deterioration

BACKGROUND

Glioma is the most common and lethal type of brain tumor, and it is characterized by unfavorable prognosis and high recurrence rates. The reprogramming of energy metabolism and an immunosuppressive tumor microenvironment (TME) are two hallmarks of tumors. Complex and dynamic interactions between neoplastic cells and the surrounding microenvironment can generate an immunosuppressive TME, which can accelerate the malignant progression of glioma. Therefore, it is crucial to explore associations between energy metabolism and the immunosuppressive TME and to identify new biomarkers for glioma prognosis.

METHODS

In our work, we analyzed the co-expression relationship between glycolytic genes and immune checkpoints based on the transcriptomic data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) and found the correlation between HK3 expression and glioma tumor immune status. To investigate the biological role of HK3 in glioma, we performed bioinformatics analysis and established a mouse glioblastoma (GBM) xenograft model.

RESULTS

Our study showed that HK3 significantly stimulated immune cell infiltration into the glioma TME. Tissue samples with higher HK3 expressive level showed increasing levels of immune cells infiltration, including M2 macrophages, neutrophils, and various subtypes of activated memory CD4+ T cells. Furthermore, HK3 expression was significantly increasing along with the elevated tumor grade, had a higher level in the mesenchymal subtype compared with those in other subtypes of GBM and could independently predict poor outcomes of GBM patients.

CONCLUSION

The present work mainly concentrated on the biological role of HK3 in glioma and offered a novel insight of HK3 regulating the activation of immune cells in the glioma microenvironment. These findings could provide a new theoretical evidence for understanding the metabolic molecular within the glioma microenvironment and identifying new therapeutic targets.

Tumor purity-associated genes influence hepatocellular carcinoma prognosis and tumor microenvironment

Introduction

Tumor purity takes on critical significance to the progression of solid tumors. The aim of this study was at exploring potential prognostic genes correlated with tumor purity in hepatocellular carcinoma (HCC) by bioinformatics analysis.

Methods

The ESTIMATE algorithm was applied for determining the tumor purity of HCC samples from The Cancer Genome Atlas (TCGA). The tumor purity-associated genes with differential expression (DEGs) were identified based on overlap analysis, weighted gene co-expression network analysis (WGCNA), and differential expression analysis. The prognostic genes were identified in terms of the prognostic model construction based on the Kaplan-Meier (K-M) survival analysis and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses. The expression of the above-described genes was further validated by the GSE105130 dataset from the Gene Expression Omnibus (GEO) database. We also characterized the clinical and immunophenotypes of prognostic genes. Gene set enrichment analysis (GSEA) was carried out for exploring the biological signaling pathway.

Results

A total of 26 tumor purity-associated DEGs were identified, which were involved in biological processes such as immune/inflammatory responses and fatty acid elongation. Ultimately, we identified ADCK3, HK3, and PPT1 as the prognostic genes for HCC. Moreover, HCC patients exhibiting higher ADCK3 expression and lower HK3 and PPT1 expressions had a better prognosis. Furthermore, high HK3 and PPT1 expressions and low ADCK3 expression resulted in high tumor purity, high immune score, high stromal score, and high ESTIMATE score. GSEA showed that the abovementioned prognostic genes showed a significant correlation with immune-inflammatory response, tumor growth, and fatty acid production/degradation.

Discussion

In conclusion, this study identified novel predictive biomarkers (ADCK3, HK3, and PPT1) and studied the underlying molecular mechanisms of HCC pathology initially.

Functions of Key Enzymes of Glycolytic Metabolism in Tumor Microenvironment

The tumor microenvironment (TME) plays a crucial role in tumor initiation, growth and metastasis. Metabolic enzymes involved in tumor glycolytic reprogramming, including hexokinase, pyruvate kinase, and lactate dehydrogenase, not only play key roles in tumorigenesis and maintaining tumor cell survival, but also take part in the modulation of the TME. Many studies have been devoted to the role of key glycolytic enzymes in the TME over the past decades. We summarize the studies on the role of glycolytic enzymes in the TME of these years and found that glycolytic enzymes remodel the TME primarily through regulating immune escape, angiogenesis, and affecting stromal cells and exosomes. Notably, abnormal tumor vascular system, peritumoral stromal cells, and tumor immunosuppressive microenvironment are important contributors to the failure of antitumor therapy. Therefore, we discuss the mechanisms of regulation by key glycolytic enzymes that may contribute to a promising biomarker for therapeutic intervention. We argue that targeting key glycolytic enzymes in combination with antiprogrammed cell death ligand 1 or antivascular endothelial growth factor could emerge as the more integrated and comprehensive antitumor treatment strategy.

CircZBTB44 promotes renal carcinoma progression by stabilizing HK3 mRNA structure

CircZBTB44 (hsa_circ_0002484) has been identified to be upregulated in renal cell carcinoma (RCC) tissues, while its role and contribution in RCC remain elusive. We confirmed the overexpression of circZBTB44 in RCC cells compared to normal kidney cell HK-2. CircZBTB44 knockdown suppressed the viability, proliferation, and migration of RCC cells and inhibited tumorigenesis in xenograft mouse models. Heterogeneous Nuclear Ribonucleoprotein C (HNRNPC) and Insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) are two RNA binding proteins of circZBTB44. HNRNPC facilitated the translocation of circZBTB44 from nuclei to cytoplasm via m6A modification, facilitating the interaction of IGF2BP3 and circZBTB44 in the cytoplasm of RCC cells. Furthermore, circZBTB44 upregulated Hexokinase 3 (HK3) expression by binding to IGF2BP3 in RCC cells. HK3 exerted oncogenic effects on RCC cell malignant behaviors and tumor growth. In the co-culture of RCC cells with macrophages, circZBTB44 promoted M2 polarization of macrophages by up-regulating HK3. In summary, HNRNPC mediated circZBTB44 interaction with IGF2BP3 to up-regulate HK3, promoting the proliferation and migration of RCC cells in vitro and tumorigenesis in vivo. The results of the study shed new light on the targeted therapy of RCC.

Aiding Cancer's "Sweet Tooth": Role of Hexokinases in Metabolic Reprogramming

Hexokinases (HKs) convert hexose sugars to hexose-6-phosphate, thus trapping them inside cells to meet the synthetic and energetic demands. HKs participate in various standard and altered physiological processes, including cancer, primarily through the reprogramming of cellular metabolism. Four canonical HKs have been identified with different expression patterns across tissues. HKs 1-3 play a role in glucose utilization, whereas HK 4 (glucokinase, GCK) also acts as a glucose sensor. Recently, a novel fifth HK, hexokinase domain containing 1 (HKDC1), has been identified, which plays a role in whole-body glucose utilization and insulin sensitivity. Beyond the metabolic functions, HKDC1 is differentially expressed in many forms of human cancer. This review focuses on the role of HKs, particularly HKDC1, in metabolic reprogramming and cancer progression.

Interfering biosynthesis by nanoscale metal-organic frameworks for enhanced radiation therapy

Radiation therapy (RT) is one of the most widely used cancer treatments. However, the vigorous biosynthesis of cancer cells plays an important role for RT resistance. Herein, we develop a hafnium-based nanoscale metal-organic frameworks (Hf-nMOFs) loaded with 3-bromopyruvate (3-BrPA) to overcome RT resistance and achieve favorable RT efficacy. The deposition of X-rays is greatly enhanced by Hf-nMOFs to induce stronger damage to DNA in RT. Simultaneously, as an inhibitor of glycolysis, the loaded 3-BrPA can reduce the supply of energy and interfere with the biosynthesis of proteins to decrease the DNA damage repair. As a result, the 3-BrPA@Hf-nMOFs (BHT) will overcome the RT resistance and enhance the curative effect of RT. Up and down-regulated genes as well as the related pathways in cellular metabolism and biosynthesis are well investigated to reveal the radiosensitization mechanism of BHT. In addition, the Hf element endows BHT with CT imaging capability to real-timely monitor the therapeutic process. Hence, the designed strategy of biosynthesis-targeted radiosensitization could decrease the doses of ionizing radiations and provide fresh perspectives on cancer treatment.

The m6A-regulation and single cell effect pattern in sunitinib resistance on clear cell renal cell carcinoma: Identification and validation of targets

Background: Sunitinib is the main target drug for clear cell renal cell carcinoma. However, the effect of sunitinib is often limited by acquired drug resistance.

Methods: The open-accessed data used in this study were obtained from different online public databases, which were analyzed using the R software. The RNA level of specific genes was detected using quantitative Real-Time PCR. Sunitinib-resistant cell lines were constructed based on protocol get from the previous study. Colony formation and Cell Counting Kit-8 assays were applied to detect cell proliferation ability.

Results: In this study, through publicly available data and high-quality analysis, we deeply explored the potential biological mechanisms that affect the resistance of sunitinib. Detailed, data from GSE64052, GSE76068 and The Cancer Genome Atlas were extracted. We identified the IFITM1, IL6, MX2, PCOLCE2, RSAD2 and SLC2A3 were associated with sunitinib resistance. Single-cell analysis, prognosis analysis and m6A regulatory network were conducted to investigate their role. Moreover, the MX2 was selected for further analysis, including its biological role and effect on the ccRCC microenvironment. Interestingly, we noticed that MX2 might be an immune-related gene that could affect the response rate of immunotherapy. Then, in vitro experiments validated the overexpression of MX2 in sunitinib-resistance cells. Colony formation assay indicated that the knockdown of MX2 could remarkably inhibit the proliferation ability of 786-O-Res and Caki-1-Res when exposed to sunitinib.

Conclusion: In summary, through publicly available data and high-quality analysis, we deeply explored the potential biological mechanisms that affect the resistance of sunitinib. MX2 was selected for further analysis, including its biological role and effect on the ccRCC microenvironment. Finally, in vitro experiments were used to validate its role in ccRCC.

Glioblastoma glycolytic signature predicts unfavorable prognosis, immunological heterogeneity, and ENO1 promotes microglia M2 polarization and cancer cell malignancy

Glioblastomas are the most malignant brain tumors, whose progress was promoted by aberrate aerobic glycolysis. The immune environment was highly engaged in glioblastoma formation, while its interaction with aerobic glycolysis remained unclear. Herein, we build a 7-gene Glycolytic Score (GS) by Elastic Net in the training set and two independent validating sets. The GS predicted malignant features and poor survival with good performances. Immune functional analyses and Cibersort calculation identified depressed T cells, B cells, natural killer cells immunity, and high immunosuppressive cell infiltration in the high-GS group. Also, high expressions of the immune-escape genes were discovered. Subsequently, the single-cell analyses validated the glycolysis-related immunosuppression. The functional results manifested the high-GS neoplastic cells' association with T cells, NK cells, and macrophage function regulation. The intercellular cross-talk showed strong associations between high-GS neoplastic cells and M2 macrophages/microglia in several immunological pathways. We finally confirmed that ENO1, the key gene of the GS, promoted M2 microglia polarization and glioblastoma cell malignant behaviors via immunofluorescence, clone formation, CCK8, and transwell rescue experiments. These results indicated the interactions between cancerous glycolysis and immunosuppression and glycolysis' role in promoting glioblastoma progression. Conclusively, we built a robust model and discovered strong interaction between GS and immune, shedding light on prognosis management improvement and therapeutic strategies development for glioblastoma patients.

Comprehensive analysis reveals signal and molecular mechanism of mitochondrial energy metabolism pathway in pancreatic cancer

Pancreatic cancer (PAAD) is one of the most malignant tumors with the worst prognosis. The abnormalities in the mitochondrial energy metabolism pathway are intimately correlated with the occurrence and progression of cancer. For the diagnosis and treatment of pancreatic cancer, abnormal genes in the mitochondrial energy metabolism system may offer new targets and biomarkers. In this study, we compared the dysregulated mitochondrial energy metabolism-associated pathways in PAAD based on pancreatic cancer samples in the Cancer Genome Atlas (TCGA) database and normal pancreas samples from the Genotype Tissue Expression project (GTEx) database. Then identified 32 core genes of mitochondrial energy metabolism pathway-related genes (MMRG) were based on the gene set enrichment analysis (GSEA). We found most of these genes were altered among different clinical characteristic groups, and showed significant prognostic value and association with immune infiltration, suggesting critical roles of MMRG involve tumor genesis of PAAD. Therefore, we constructed a four-gene (LDHA, ALDH3B1, ALDH3A1, and ADH6) prognostic biomarker after eliminating redundant factors, and confirming its efficiency and independence. Further analysis indicated the potential therapeutic compounds based on the mitochondrial energy metabolism-associated prognostic biomarker. All of the above analyses dissected the critical role of mitochondrial energy metabolism signaling in pancreatic cancer and gave a better understanding of the clinical intervention of PAAD.

Hexokinases in cancer and other pathologies

Glucose metabolism is indispensable for cell growth and survival. Hexokinases play pivotal roles in glucose metabolism through canonical functions of hexokinases as well as in immune response, cell stemness, autophagy, and other cellular activities through noncanonical functions. The aberrant regulation of hexokinases contributes to the development and progression of pathologies, including cancer and immune diseases.

Exploring the prognostic value of HK3 and its association with immune infiltration in glioblastoma multiforme

Background: Hexokinase 3 (HK3) is one of the key enzymes involved in glucose phosphorylation (the first step in most glucose metabolic pathways). Many studies have demonstrated the vital role of dysregulation of HK3 in several tumors. However, there is a need for in-depth characterization of the role of HK3 in glioblastoma multiforme (GBM). Methods: All data were sourced from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA). Kaplan-Meier analysis and univariate regression were applied for survival analysis. Gene set enrichment analysis (GSEA) was used for enrichment analysis. Tumor Immune Single Cell Hub (TISCH) database was applied for single-cell analysis. Tumor Immune Dysfunction and Exclusion (TIDE) analysis was applied to evaluate the immune response. Results: HK3 expression was upregulated in GBM and correlated with poor prognosis. The high HK3 expression group was primarily enriched in adaptive immune response, chemokine signaling pathway, and cytokine-cytokine receptor interaction. The high HK3 expression group showed significantly greater enrichment of the majority of immune cells and immune-related pathways. HK3 showed significant correlation with most immune cells, especially macrophages (p < .001, R = .81). TISCH analysis showed that HK3 was predominantly expressed in macrophages in most cancers. HK3 showed significant correlation with most immune-related genes, such as PD-1 (p < .001, R = .41), PDL-1 (p < .001, R = .27), and CTLA-4 (p < .001, R = .29). TIDE analysis revealed that the low HK3 expression group has a lower TIDE score and may benefit from immunotherapy. Drug sensitivity analysis showed that patients with high HK3 expression frequently showed drug resistance. Conclusion: HK3 was associated with poor prognosis and may serve as a biomarker of macrophages in GBM. HK3 was also associated with immune response and drug resistance. Our findings may provide novel insights for GBM immunotherapy.

A risk stratification and prognostic prediction model for lung adenocarcinoma based on aging-related lncRNA

To create a risk model of aging-related long non-coding RNAs (arlncRNAs) and determine whether they might be useful as markers for risk stratification, prognosis prediction, and targeted therapy guidance for patients with lung adenocarcinoma (LUAD). Data on aging genes and lncRNAs from LUAD patients were obtained from Human Aging Genomic Resources 3 and The Cancer Genome Atlas, and differential co-expression analysis of established differentially expressed arlncRNAs (DEarlncRNAs) was performed. They were then paired with a matrix of 0 or 1 by cyclic single pairing. The risk coefficient for each sample of LUAD individuals was obtained, and a risk model was constructed by performing univariate regression, least absolute shrinkage and selection operator regression analysis, and univariate and multivariate Cox regression analysis. Areas under the curve were calculated for the 1-, 3-, and 5-year receiver operating characteristic curves to determine Akaike information criterion-based cutoffs to identify high- and low-risk groups. The survival rate, correlation of clinical characteristics, malignant-infiltrating immune-cell expression, ICI-related gene expression, and chemotherapeutic drug sensitivity were contrasted with the high- and low-risk groups. We found that 99 DEarlncRNAs were upregulated and 12 were downregulated. Twenty pairs of DEarlncRNA pairs were used to create a prognostic model. The 1-, 3-, and 5-year survival curve areas of LUAD individuals were 0.805, 0.793, and 0.855, respectively. The cutoff value to classify patients into two groups was 0.992. The mortality rate was higher in the high-risk group. We affirmed that the LUAD outcome-related independent predictor was the risk score (p < 0.001). Validation of tumor-infiltrating immune cells and ICI-related gene expression differed substantially between the groups. The high-risk group was highly sensitive to docetaxel, erlotinib, gefitinib, and paclitaxel. Risk models constructed from arlncRNAs can be used for risk stratification in patients with LUAD and serve as prognostic markers to identify patients who might benefit from targeted and chemotherapeutic agents.

The significance of glycolysis in tumor progression and its relationship with the tumor microenvironment

It is well known that tumor cells rely mainly on aerobic glycolysis for energy production even in the presence of oxygen, and glycolysis is a known modulator of tumorigenesis and tumor development. The tumor microenvironment (TME) is composed of tumor cells, various immune cells, cytokines, and extracellular matrix, among other factors, and is a complex niche supporting the survival and development of tumor cells and through which they interact and co-evolve with other tumor cells. In recent years, there has been a renewed interest in glycolysis and the TME. Many studies have found that glycolysis promotes tumor growth, metastasis, and chemoresistance, as well as inhibiting the apoptosis of tumor cells. In addition, lactic acid, a metabolite of glycolysis, can also accumulate in the TME, leading to reduced extracellular pH and immunosuppression, and affecting the TME. This review discusses the significance of glycolysis in tumor development, its association with the TME, and potential glycolysis-targeted therapies, to provide new ideas for the clinical treatment of tumors.

Proteomic characteristics of saliva in patients with different subgroups of IgG4-RD

Background

Immunoglobulin G4-related disease (IgG4-RD) is a newly defined disease entity, with great heterogeneity among IgG4-RD subgroups with different organ involvement patterns. Identification of the proteomic characteristics of IgG4-RD subgroups will be critical for the understanding of the pathogenic mechanisms of IgG4-RD.

Method

In this study, we performed proteomic analysis using Tandem Mass Tags (TMT) technology with "high field" mass analyzer with improved resolution and sequencing speed to investigate the proteomic profile of saliva and plasma samples from ten untreated IgG4-RD patients and five healthy controls (HCs). Differentially expressed proteins (DEPs) were identified by "t test" function in R package. Functional enrichment analysis was used to investigate pathways enriched in IgG4-RD samples.

Results

Most salivary DEPs identified in IgG4-RD patients compared with HCs were mainly enriched in neutrophil mediated GO bioprocess. Within the comparisons between four IgG4-RD subgroups, more DEPs were identified in the comparison of Mikulicz group and Head and neck group. Among four subgroups of IgG4-RD, Head and neck group showed the most distinctive proteomic expression pattern when compared with HCs. Moreover, "Neutrophil mediated process" related GO bioprocess was commonly identified between comparisons of Mikulicz group and Head and neck group, Head and neck group and Retroperitoneal aorta group, Head and neck group and HCs, IgG4-RD patients with saliva gland involvement and those without saliva gland involvement. Key DEPs that involved in this GO bioprocess were identified. Besides, we performed proteomic analysis for plasma samples between ten IgG4-RD and five HCs and there were several DEPs identified overlapped in saliva and plasma.

Conclusion

We identified multiple processes/factors and several signaling pathways in saliva that may be involved in the IgG4-RD pathogenesis.

Integrated Analysis of the Role of Enolase 2 in Clear Cell Renal Cell Carcinoma

Enolase 2 (ENO2) has increasingly been documented in multiple cancers in recent years. However, the role of ENO2 in clear cell renal carcinoma (ccRCC) has not been fully explored. In the present study, open-access data were downloaded from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and the Human Protein Atlas (HPA) databases. All statistical analyses were performed in R and GraphPad Prism 8 softwares. Results showed that ENO2 was overexpressed in ccRCC tissues and cell lines and correlated with worse clinical features and prognosis. In vitro experiments indicated that the inhibition of ENO2 could hamper the malignant behaviors of ccRCC cells. Gene Set Enrichment Analysis showed that epithelial-mesenchymal transition, KRAS signaling, inflammatory response, angiogenesis, hypoxia, and WNT/β-catenin pathways were upregulated in the ENO2 high-expression group; whereas adipogenesis, DNA repair, and androgen response pathways were downregulated. Immune infiltration analysis indicated that patients with high ENO2 levels might have higher M2 macrophages and lower γβ T cells in the tumor microenvironment, which may account to some extent for the worse prognosis of ENO2. Moreover, it was found that patients with low and high ENO2 expression might be more sensitive to PD-1 therapy and CTLA-4 therapy, respectively. In addition, patients with high ENO2 expression showed lower sensitivity to common chemotherapy drugs for ccRCC, including axitinib, cisplatin, gemcitabine, pazopanib, sunitinib, and temsirolimus. Overall, these results suggest that ENO2 is a potential prognosis biomarker of ccRCC and could affect the malignant biological behavior of cancer cells, highlighting its value as a potential therapeutic target.

Proteomes of Residual Tumors in Curcumin-Treated Rats Reveal Changes in Microenvironment/Malignant Cell Crosstalk in a Highly Invasive Model of Mesothelioma

Curcumin exhibits both immunomodulatory properties and anticarcinogenic effects which have been investigated in different experimental tumor models and cancer types. Its interactions with multiple signaling pathways have been documented through proteomic studies on malignant cells in culture; however, in vivo approaches are scarce. In this study, we used a rat model of highly invasive peritoneal mesothelioma to analyze the residual tumor proteomes of curcumin-treated rats in comparison with untreated tumor-bearing rats (G1) and provide insights into the modifications in the tumor microenvironment/malignant cell crosstalk. The cross-comparing analyses of the histological sections of residual tumors from two groups of rats given curcumin twice on days 21 and 26 after the tumor challenge (G2) or four times on days 7, 9, 11 and 14 (G3), in comparison with G1, identified a common increase in caveolin-1 which linked with significant abundance changes affecting 115 other proteins. The comparison of G3 vs. G2 revealed additional features for 65 main proteins, including an increase in histidine-rich glycoprotein and highly significant abundance changes for 22 other proteins regulating the tumor microenvironment, linked with the presence of numerous activated T cells. These results highlight new features in the multiple actions of curcumin on tumor microenvironment components and cancer cell invasiveness.

The expression and prognostic value of transporter 1, ATP binding cassette subfamily B member in clear cell renal cell cancer with experimental validation

Transporter associated with antigen processing 1(TAP1) serves as a protein to transport antigenic peptides from the surface of the endoplasmic reticulum to the lumen of the endoplasmic reticulum when the antigens are presented by major histocompatibility complex type I (MHC-I), which has been identified to play a critical role in antigen presentation in innate immunity. In tumors, the role of TAP1 seems to remain controversial. On the one hand, given the role of TAP1 in antigen presentation, it is indicated that high TAP1 expression corresponds to the emergence of more neoantigens epitopes that facilitate the recognition for phagocytes, T cells and other cells. On the other hand, the genetic ablation of transporter associated with antigen processing (TAP) results in the presentation of new class I-restricted epitopes encoded in house-keeping products. Opposite result has been revealed by studies in other tumors suggest, which implies a more complex function of TAP1. Therefore, it's significant to clarify the role of TAP1 in clear cell renal cell carcinoma (ccRCC). In this study, we found the elevated expression levels in mRNA and protein of TAP1 in ccRCC tissues, which indicated a relatively worse prognosis. Transwell assay and Scratch assay in vitro demonstrated the promotive role of TAP1 in ccRCC migration as well as a significant role in metastasis. And the increased expression of TAP1 resulted in more immune cells infiltrated in cancer tissues. TAP1 was also demonstrated to be related to immune regulator genes, as gene set enrichment analysis (GSEA) indicated its significant role in immune regulation. The results of CancerSEA indicated the positive association of the high-level TAP1 expression with epithelial-mesenchymal transition (EMT) and the inverse association with Cell Cycle. The effective drugs were also predicted based on TAP1 expression, of which the high level was indeed associated with resistance to multiple drugs, but some effective drugs still identified based on high TAP1 expression. According to the analysis of various databases, the role of TAP1 in ccRCC was explored, especially in relationship of TAP1 with tumor microenvironment. These results indicate that TAP1 can serve as a potential target for treatment of ccRCC.

Reprogramming lipid metabolism as potential strategy for hematological malignancy therapy

Hematological malignancies are one of the most lethal illnesses that seriously threaten human life and health. Lipids are important constituents of various biological membranes and substances for energy storage and cell signaling. Furthermore, lipids are critical in the normal physiological activities of cells. In the process of the lethal transformation of hematological malignancies, lipid metabolism reprogramming meets the material and energy requirements of rapidly proliferating and dividing tumor cells. A large number of studies have shown that dysregulated lipid metabolism, commonly occurs in hematological malignancies, mediating the proliferation, growth, migration, invasion, apoptosis, drug resistance and immune escape of tumor cells. Targeting the lipid metabolism pathway of hematological malignancies has become an effective therapeutic approach. This article reviews the oncogenic mechanisms of lipid metabolism reprogramming in hematological malignancies, including fatty acid, cholesterol and phospholipid metabolism, thereby offering an insight into targeting lipid metabolism in the treatment of hematological malignancies.

Genomic alteration of MTAP/CDKN2A predicts sarcomatoid differentiation and poor prognosis and modulates response to immune checkpoint blockade in renal cell carcinoma

Sarcomatoid differentiation is a highly aggressive pathological characteristic of renal cell carcinoma (RCC) and is characterized by susceptibility to progression and extremely poor prognosis. In this study, we included all genomic alteration events that led to a loss of protein function of MTAP and CDKN2A, and enrolled 5,307 RCC patients with genomic sequencing data from Western and Chinese cohorts. Notably, MTAP/CDKN2AMUT occurred in the Chinese population ~2 times more frequently than in the Western cohort and showed significant co-mutation trends. We found significantly higher proportions of sarcomatoid-positive patients with MTAPMUT or CDKN2AMUT compared with MTAP/CDKN2A wild-type (WT) patients (P < 0.001). Of the 574 RCC samples from the FUSCC cohort and 3,563 RCC samples from 17 independent cohorts, the MTAP/CDKN2AMUT significantly predicted extremely poor outcomes (P < 0.0001). The Western cohort suggested a concordant relationship between MTAP/CDKN2AMUT and sarcomatoid differentiation in RCC. Moreover, although MTAP/CDKN2AMUT RCC may be insensitive to targeted therapy, the high degree of tumor heterogeneity and higher PD-L1 and CXCL13 expression characterizations reflected that MTAP/CDKN2A-deficient features could benefit from immunotherapy for patients with RCC. This study utilized RCC samples from large-scale, global, multicenter sequencing cohorts and first proved that MTAP/CDKN2A deficiency significantly correlates with sarcomatoid differentiation in RCC and predicts aggressive progression, poor prognosis, and primary resistance to targeted therapy and potential favorable responses to immune checkpoint blockade. Unlike conventional targeted therapies, emerging drugs such as immunotherapies or synthetic lethal PRMT5 inhibitors may become novel therapeutic options for patients with MTAP/CDKN2AMUT RCC.

Prognostic value, DNA variation and immunologic features of a tertiary lymphoid structure-related chemokine signature in clear cell renal cell carcinoma

BACKGROUND

The tumor microenvironment (TME) and tertiary lymphoid structures (TLS) affect the occurrence and development of cancers. How the immune contexture interacts with the phenotype of clear cell renal cell carcinoma (ccRCC) remains unclear.

METHODS

We identified and evaluated TLS clusters in ccRCC using machine learning algorithms and the 12-chemokine gene signature for TLS. Analyses for functional enrichment, DNA variation, immune cell distribution, association with independent clinicopathological features and predictive value of CXCL13 in ccRCC were performed.

RESULTS

We found a prominently enrichment of the 12-chemokine gene signature for TLS in patients with ccRCC compared with other types of renal cell carcinoma. We identified a prognostic value of CCL4, CCL5, CCL8, CCL19 and CXCL13 expression in ccRCC. DNA deletion of the TLS gene signature significantly predicted poor outcome in ccRCC compared with amplification and wild-type gene signature. We established TLS clusters (C1-4) and observed distinct differences in survival, stem cell-like characteristics, immune cell distribution, response to immunotherapies and VEGF-targeted therapies among the clusters. We found that elevated CXCL13 expression significantly predicted aggressive progression and poor prognosis in 232 patients with ccRCC in a real-world validation cohort.

CONCLUSION

This study described a 12-chemokine gene signature for TLS in ccRCC and established TLS clusters that reflected different TME immune status and corresponded to prognosis of ccRCC. We confirmed the dense presence of TILs aggregation and TLS in ccRCC and demonstrated an oncogenic role of CXCL13 expression of ccRCC, which help develop immunotherapies and provide novel insights on the long-term management of ccRCC.

Tumor-associated macrophage-derived chemokine CCL5 facilitates the progression and immunosuppressive tumor microenvironment of clear cell renal cell carcinoma

Background: Tumor-associated macrophages (TAMs) dominate the malignancy of cancers by perturbing the tumor microenvironment (TME). However, the clinical implications of heterogeneous subpopulations of TAMs in clear cell renal cell carcinoma (ccRCC) remain to be elucidated. Methods: We comprehensively evaluated the prognostic implications, biological behaviors, and immunogenomics features of the C-C Motif Chemokine Ligand 5 (CCL5) expression and CCL5+ TME in vitro and in 932 real-world ccRCC patients from testing and public validation cohorts. Flow cytometry was used to examine the functional patterns of CCL5+ TAMs with TME cell-infiltrating characterizations. Results: Our results identified distinct prognostic clusters with gradual changes in clinicopathological indicators based on CCL5 expression. Knockdown of CCL5 significantly restrained cell viability, migration capabilities of ccRCC cells, and the inhibits the proliferation and chemotaxis of THP1-derived TAMs. Mechanically, down-regulation of CCL5 arrested epithelial-mesenchymal transition by modulating the PI3K/AKT pathway in ccRCC cells. In ccRCC samples with CCL5 upregulation, the proportion of CCL5+ TAMs and PD-L1+ CD68+ TAMs were prominently increased, showing a typical suppressive tumor immune microenvironment (TIME). Besides, intra-tumoral CCL5+ TAMs showed distinct pro-tumorigenic TME features characterized by exhausted CD8+ T cells and increased expression of immune checkpoints. Furthermore, elevated CCL5+ TAMs infiltration was prominently associated with a dismal prognosis for patients with ccRCC. Conclusion: In conclusion, this study first revealed the predictive value of the chemokine CCL5 on the progression and TME of ccRCC. The intra-tumoral CCL5+ TAMs could be applied to comprehensively evaluate the prognostic patterns as well as unique TME characteristics among individuals, allowing for the identification of immunophenotypes and promotion of treatment efficiency for ccRCC.

Deciphering the role of miR-187-3p/LRFN1 axis in modulating progression, aerobic glycolysis and immune microenvironment of clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is one of the most common malignant genitourinary cancers with high recurrence risk worldwide. Recently, multi-omics data facilitate obtaining a molecular landscape of tumor development, and were implemented to affect pathogenesis, phenotype, and prognosis of ccRCC. In this study, after screening for differential expressed microRNAs based on multiply datasets, we tested expression levels and prognostic value of miR-187-3p in ccRCC samples, and transfected miR-187-3p mimics or negative controls into ccRCC cells. Up-regulation of miR-187-3p restrains proliferation, migration and promotes apoptosis ability in human ccRCC A498 and 786O cells. In addition, Luciferase reporter assay revealed that miR-187-3p directly targets LRFN1-3'-UTR and negatively modulates LRFN1 expression. LRFN1 rescues proliferation and invasion capacities after miR-187-3p mimic transfection in vitro and in subcutaneous xenograft models. We further performed deep-sequencing technology and bioinformatics analyses to evaluate the biological functions and potential clinical implications of LRFN1 expression in ccRCC. Interestingly, LRFN1 could serve as an independent and potential biomarker for prognosis in over 1000 patients with ccRCC from multiply independent cohorts. Besides, the up-regulated LRFN1 expression prominently promoted intra-tumoral heterogeneity and immune-infiltrating microenvironment, represented by elevated M2 macrophage infiltration, CD8+ T cells activity and PD-L1 expression. In conclusion, this study revealed the tumor-specific and immunological role of miR-187-3p/LRFN1 axis in the progression and reshaping of tumor immune microenvironment of ccRCC.

Risk coefficient model of necroptosis-related lncRNA in predicting the prognosis of patients with lung adenocarcinoma

Model algorithms were used in constructing the risk coefficient model of necroptosis-related long non-coding RNA in identifying novel potential biomarkers in the prediction of the sensitivity to chemotherapeutic agents and prognosis of patients with lung adenocarcinoma (LUAD). Clinic and transcriptomic data of LUAD were obtained from The Cancer Genome Atlas. Differently expressed necroptosis-related long non-coding RNAs got identified by performing both the univariate and co-expression Cox regression analyses. Subsequently, the least absolute shrinkage and selection operator technique was adopted in constructing the nrlncRNA model. We made a comparison of the areas under the curve, did the count of the values of Akaike information criterion of 1-year, 2-year, as well as 3-year receiver operating characteristic curves, after which the cut-off value was determined for the construction of an optimal model to be used in identifying high risk and low risk patients. Genes, tumor-infiltrating immune cells, clinical correlation analysis, and chemotherapeutic agents data of both the high-risk and low-risk subgroups were also performed. We identified 26 DEnrlncRNA pairs, which were involved in the Cox regression model constructed. The curve areas under survival periods of 1 year, 2 years, and 3 years of patients with LUAD were 0.834, 0.790, and 0.821, respectively. The cut-off value set was 2.031, which was used in the identification of either the high-risk or low-risk patients. Poor outcomes were observed in patients belonging to the high-risk group. The risk score was the independent predictor of the LUAD outcome (p < 0.001). The expression levels of immune checkpoint and infiltration of specific immune cells were anticipated by the gene risk model. The high-risk group was found to be highly sensitive to docetaxel, erlotinib, cisplatin, and paclitaxel. The model established through nrlncRNA pairs irrespective of the levels of expression could give a prediction on the LUAD patients’ prognosis and assist in identifying the patients who might gain more benefit from chemotherapeutic agents.

Comprehensive Analysis of the Oncogenic, Genomic Alteration, and Immunological Landscape of Cation-Chloride Cotransporters in Pan-Cancer

Background

Assessing the phenotypic diversity underlying tumor progression requires the identification of variations in the respective molecular interaction in the tumor microenvironment (TME). Despite emerging studies focusing on the association between cation-chloride cotransporters (CCCs) and carcinogenesis, direct evidence that CCCs (KCC2 and NKCC1) mediate tumor progression in pan-cancer remains unclear.

Methods

We conducted a comprehensive assessment of the expression, DNA variation profiles, and prognostic and immunologic implications of CCCs based on a large-scale pan-cancer population, including 10,967 cancer patients from the Cancer Genome Atlas, 9,162 cancer patients from Genomics Expression Omnibus, 48,834 cancer patients from 188 independent studies, and 356 cancer patients from three real-world cohorts.

Results

In this study, we first found that CCCs were highly expressed in most tumors, and prominently associated with prognosis. Kaplan-Meier analysis and Cox regression analysis revealed that KCC2 and NKCC1 significantly predicted survival for patients with pan-cancer, suggesting that CCCs have inconsistent tumorigenesis regulatory mechanisms in cancers. Next, we examined the DNA variation landscape of KCC2 and NKCC1 and their prognostic implications in pan-cancer. The results demonstrated that UCEC patients with somatic copy number variation (CNV) of NKCC1 received significantly better outcomes (p < 0.05). Besides emphasizing the clinical implications of CNV of CCCs for cancer patients, we found that NKCC1MUT could prominently prolong progression-free survival (p = 2.59e-04), disease-specific survival (p = 0.019), and overall survival (p = 0.034) compared with NKCC1WT cancer patients possibly via regulation of cell proliferation and oncogenic stress pathways. Additionally, KCC2 positively correlated with the levels of tumor-infiltrating macrophages and CD4+ T cells, but NKCC1 showed a significantly widely negative association with tumor-infiltrated lymphocytes, suggesting an immune-excluded TME in cancers. Similarly, expression of KCC2, rather than NKCC1, was positively correlated with the immune checkpoint molecules, indicating its role as an immune regulator in a wide variety of cancers. Finally, to verify our hypothesis and altered expression of CCCs, we performed IHC analysis and revealed the staining distribution in tumor and adjacent normal tissues of glioma, clear cell renal cell carcinoma, papillary cell renal cell carcinoma, and hepatocellular and breast cancer from three real-world cohorts, and validated prominently prognostic implications of CCCs in patients with clear cell renal cell carcinoma.

Conclusion

This study first comprehensively investigated the molecular and clinical role of CCCs, and illustrated the significant association among KCC2/NKCC1 expression, DNA variation profiles prognosis, and TME of pan-cancer. The pan-cancer findings provided an in-depth understanding of potential oncogenic and immunologic of differential expression and DNA alteration of KCC2/NKCC1 cancers.

An immune-related lncRNA risk coefficient model to predict the outcomes in clear cell renal cell carcinoma

Objective: Using model algorithms, we constructed an immune-related long non-coding RNAs (lncRNAs) risk coefficient model to predict outcomes for patients with clear cell renal cell carcinoma (ccRCC) to understand the infiltration of tumor immune cells and the sensitivity to immune-targeted drugs.

Methods: Open genes data were downloaded from The Cancer Genome Atlas and The Immunology Database and Analysis Portal, and immune-related lncRNAs were obtained through Pearson correlation analysis. R language software was used to obtain differentially expressed immune-related lncRNAs and immune-related lncRNA pairs. The model was constructed using least absolute shrinkage and selector operation regression analysis, and receiver operator characteristic curves were drawn. The Akaike information criterion was used to distinguish the high-risk from the low-risk group. We also conducted correlation analysis for the high- and low-risk subgroups.

Results: We identified 27 immune-related lncRNAs pairs, 16 of which were included in the model construction. After merging clinical data, the areas under the curve of 1 -year, 3-year, and 5-year survival times of ccRCC patients were 0.867, 0.832, and 0.838, respectively. Subgroup analyses were conducted according to the cut-off value. We found that the high-risk group was associated with poor outcomes. The risk score and tumor stage were independent predictors of the outcome of ccRCC. The risk model predicted specific immune cell infiltration, immune checkpoint gene expression levels, and high-risk groups more sensitive to sunitinib targeted therapy.

Conclusion: We obtained prognostic-related novel ccRCC markers and risk model that predicts the outcome of patients with ccRCC and helps identify those who can benefit from sunitinib.

Prognostic Value of Long Noncoding RNA DLEU2 and Its Relationship with Immune Infiltration in Kidney Renal Clear Cell Carcinoma and Liver Hepatocellular Carcinoma

Background

DLEU2 is a long noncoding RNA considered important in the progression of many cancers. However, correlations between DLEU2 and kidney renal clear cell carcinoma (KIRC) and liver hepatocellular carcinoma (LIHC) have rarely been reported.

Methods

We first analysed the expression of DLEU2 across cancers and the correlation between DLEU2 and the clinical features of KIRC and LIHC by using the “ggplot2” package in R and searched the Oncomine database and Timer website platform. We verified the expression of DLEU2 in the GEO dataset (GSE105261 and GSE45267). Receiver operating characteristic (ROC) curves were drawn using the “pROC” and “ggplot2” packages in R, and we constructed a DLEU2-based prognostic nomogram for KIRC and LIHC by using the “survival” and “rms” packages in R. Then, we analysed the correlation between DLEU2 expression and prognosis in R as well as the correlation between DLEU2 and immune cell infiltration in the TIMER database. Finally, we explored the causes of DLEU2 upregulation in the UCSC Xena and UALCAN databases.

Results

We found that DLEU2 was upregulated in many cancers, including KIRC and LIHC. Expression of DLEU2 is associated with tumour stage, grade, lymphatic metastasis, and distant metastasis in KIRC as well as alpha-fetoprotein (AFP), tumour stage, grade, lymphatic metastasis, and distant metastasis in LIHC. DLEU2 is an adverse factor for the prognosis of KIRC and LIHC. In addition, DLEU2 has moderate accuracy in diagnosing KIRC and LIHC and predicting their prognosis. Moreover, we found that expression of DLEU2 correlated positively with immune cell infiltration in KIRC and LIHC, and upregulation of DLEU2 in KIRC and LIHC suggests a poor prognosis based on immune cells analysis. Genetic and epigenetic analyses of DLEU2 indicate that copy number variations (CNVs) and methylation contribute to the upregulation of DLEU2.

Conclusion

The long noncoding RNA DLEU2 has the potential to predict the prognosis and immune infiltration of KIRC and LIHC.

Comprehensive Multi-Omics Identification of Interferon-γ Response Characteristics Reveals That RBCK1 Regulates the Immunosuppressive Microenvironment of Renal Cell Carcinoma

Interferon-gamma (IFN-γ) has a complex role in modulating the tumor microenvironment (TME) during renal cell carcinoma (RCC) development. To define the role of IFN-γ response genes in RCC progression, we characterized the differential gene expression, prognostic implications, and DNA variation profiles of selected IFN-γ response signatures, which exhibited a significant hazard ratio for the overall survival (OS) and progression-free survival (PFS) of papillary, chromophobia, and clear cell RCC (ccRCC) patients (n = 944). Prognostic nomograms were constructed to predict the outcomes for ccRCC patients, highlighting the prognostic implications of RANBP2-type and C3HC4-type zinc finger containing 1 (RBCK1). Interestingly, large-scale pan-cancer samples (n = 12,521) and three single-cell RNA datasets revealed that RBCK1 showed markedly differential expression between cancer and normal tissues and significantly correlated with tumor-infiltrating immune cells, tumor purity, and immune checkpoint molecules, such as PD-L1, CTLA-4, LAG-3, and TIGIT in pan-cancer samples. Notably, the TIDE score was significantly higher in the RBCK1high group compared with the RBCK1low group in both ccRCC and RCC cohorts. Besides, immunohistochemistry staining showed significantly elevated RBCK1 expression in tumors (n = 50) compared with kidney samples (n = 40) from a real-world cohort, Fudan University Shanghai Cancer Center (FUSCC, Shanghai). After RBCK1 expression was confirmed in ccRCC, we found a significantly decreased number of infiltrating CD4+ T cells, CD4+ FOXP3+ Treg cells, M1 macrophages, and CD56bight/dim NK cells in the immune-cold RBCK1high group. In addition to the distinct heterogeneous immune microenvironment, the increased expression of RBCK1 predicted a prominently worse prognosis than the RBCK1low group for 232 ccRCC patients in the FUSCC proteomic cohort. Furthermore, after transfected with siRNA in human ccRCC cells, extraordinarily decreased cell proliferation, migration capacities, and prominently elevated apoptosis tumor cell proportion were found in the siRNA groups compared with the negative control group. In conclusion, this study identified IFN-γ response clusters, which might be used to improve the prognostic accuracy of immune contexture in the ccRCC microenvironment. Immune-cold RBCK1high patients have pro-tumorigenic immune infiltration and significantly worse outcomes than RBCK1low patients based on results from multi-omics to real-world data. Our discovery of novel independent prognostic indicators for RCC highlights the association between tumor alterations and immune phenotype.

Malignant Tumor Purity Reveals the Driven and Prognostic Role of CD3E in Low-Grade Glioma Microenvironment

The tumor microenvironment (TME) contributes to the initiation and progression of many neoplasms. However, the impact of low-grade glioma (LGG) purity on carcinogenesis remains to be elucidated. We selected 509 LGG patients with available genomic and clinical information from the TCGA database. The percentage of tumor infiltrating immune cells and the tumor purity of LGG were evaluated using the ESTIMATE and CIBERSORT algorithms. Stromal-related genes were screened through Cox regression, and protein-protein interaction analyses and survival-related genes were selected in 487 LGG patients from GEO database. Hub genes involved in LGG purity were then identified and functionally annotated using bioinformatics analyses. Prognostic implications were validated in 100 patients from an Asian real-world cohort. Elevated tumor purity burden, immune scores, and stromal scores were significantly associated with poor outcomes and increased grade in LGG patients from the TCGA cohort. In addition, CD3E was selected with the most significant prognostic value (Hazard Ratio=1.552, P<0.001). Differentially expressed genes screened according to CD3E expression were mainly involved in stromal related activities. Additionally, significantly increased CD3E expression was found in 100 LGG samples from the validation cohort compared with adjacent normal brain tissues. High CD3E expression could serve as an independent prognostic indicator for survival of LGG patients and promotes malignant cellular biological behaviors of LGG. In conclusion, tumor purity has a considerable impact on the clinical, genomic, and biological status of LGG. CD3E, the gene for novel membrane immune biomarker deeply affecting tumor purity, may help to evaluate the prognosis and develop individual immunotherapy strategies for LGG patients. Evaluating the ratio of differential tumor purity and CD3E expression levels may provide novel insights into the complex structure of the LGG microenvironment and targeted drug development.

m6A Regulator-Mediated Methylation Modification Model Predicts Prognosis, Tumor Microenvironment Characterizations and Response to Immunotherapies of Clear Cell Renal Cell Carcinoma

BACKGROUND

This study aims to establish an N6-methyladenosine (m6A) RNA methylation regulators-mediated methylation model and explore its role in predicting prognostic accuracy of immune contexture and characterizations of clear cell renal cell carcinoma (ccRCC).

METHODS

The m6A modification subclasses (m6AMS) were identified by unsupervised cluster analysis and three clusters were determined by consensus clustering algorithm in a discovering cohort. Testing and real-world validation cohorts were used to identify predictive responses for immune checkpoint therapies (ICTs) of m6AMS.

RESULTS

Prognostic implications landscape of m6A regulators in cancers and its differential expression levels in ccRCC patients were identified. Based on discovering cohort, ccRCC were automatically divided into three m6AMS, and cluster 3 showed significant worse survival than cluster 1/2. Importantly, it was found that the immune checkpoint molecules expression was significantly elevated in cluster 3. Besides, m6A scoreLow group (cluster 1&2) have significantly elevated TIDE score compared with m6A scoreHigh group (cluster 3). There was conspicuous tertiary lymphoid tissue, aggressive phenotype, elevated glycolysis, expression of PD-L1, abundance of CD8+ T cells, CD4+ FOXP3+ Treg cells and TCRn immune cells infiltration in the high m6A score group. Interestingly, there are significantly increased patients with clinical benefit in m6A scoreHigh group in 368 patients receiving ICTs from testing IMvigor210 (n = 292) and validation FUSCC (n = 55) cohorts.

CONCLUSION

Our discovery highlights the relationship between tumor epigenetic heterogeneity and immune contexture. Immune-rejection cluster 3 has pro-tumorigenic immune infiltration, and shows significant clinical benefits for ccRCC patients receiving ICTs, enabling patient selection for future clinical treatment.