Pan-cancer analysis reveals that G6PD is a prognostic biomarker and therapeutic target for a variety of cancers

Abstract

Background: Despite accumulating evidence revealing that Glucose-6-phosphate dehydrogenase (G6PD) is highly expressed in many tumor tissues and plays a remarkable role in cancer tumorigenesis and progression, there is still a lack of G6PD pan-cancer analysis. This study was designed to analyze the expression status and prognostic significance of G6PD in pan-cancer.

Methods: G6PD expression data were obtained from multiple data resources including the Genotype-Tissue Expression, the Cancer Genome Atlas, and the Tumor Immunity Estimation Resource. These data were used to assess the G6PD expression, prognostic value, and clinical characteristics. The ESTIMATE algorithms were used to analyze the association between G6PD expression and immune-infiltrating cells and the tumor microenvironment. The functional enrichment analysis was also performed across pan-cancer. In addition, the GDSC1 database containing 403 drugs was utilized to explore the relationship between drug sensitivity and G6PD expression levels. Furthermore, we also performed clinical validation and in vitro experiments to further validate the role of G6PD in hepatocellular carcinoma (HCC) cells and its correlation with prognosis. The R software was used for statistical analysis and data visualization.

Results: G6PD expression was upregulated in most cancers compared to their normal counterparts. The study also revealed that G6PD expression was a prognostic indicator and high levels of G6PD expression were correlated with worse clinical prognosis including overall survival, disease-specific survival, and progression-free interval in multiple cancers. Furthermore, the G6PD level was also related to cancer immunity infiltration in most of the cancers, especially in KIRC, LGG, and LIHC. In addition to this, G6PD expression was positively related to pathological stages of KIRP, BRCA, KIRC, and LIHC. Functional analysis and protein-protein interactions network results revealed that G6PD was involved in metabolism-related activities, immune responses, proliferation, and apoptosis. Drug sensitivity analysis showed that IC50 values of most identified anti-cancer drugs were positively correlated with the G6PD expression. Notably, in vitro functional validation showed that G6PD knockdown attenuated the phenotypes of proliferation in HCC.

Conclusion: G6PD may serve as a potential prognostic biomarker for cancers and may be a potential therapeutic target gene for tumor therapy.

Keywords: G6PD; drug resistance; pan-cancer; survival analysis; tumor microenvironment.