- Primary research
- Open Access
Xanthine dehydrogenase as a prognostic biomarker related to tumor immunology in hepatocellular carcinoma
Cancer Cell International volume 21, Article number: 475 (2021)
Xanthine dehydrogenase (XDH) is a critical enzyme involved in the oxidative metabolism of purines, pterin and aldehydes and a central component of the innate immune system. However, the prognostic value of XDH in predicting tumor-infiltrating lymphocyte abundance, the immune response, and survival in different cancers, including hepatocellular carcinoma (HCC), is still unclear.
XDH expression was analyzed in multiple databases, including Oncomine, the Tumor Immune Estimation Resource (TIMER), the Kaplan–Meier plotter database, the Gene Expression Profiling Interactive Analysis (GEPIA) database, and The Cancer Genome Atlas (TCGA). XDH-associated transcriptional profiles were detected with an mRNA array, and the levels of infiltrating immune cells were validated by immunohistochemistry (IHC) of HCC tissues. A predictive signature containing multiple XDH-associated immune genes was established using the Cox regression model.
Decreased XDH mRNA expression was detected in human cancers originating from the liver, bladder, breast, colon, bile duct, kidney, and hematolymphoid system. The prognostic potential of XDH mRNA expression was also significant in certain other cancers, including HCC, breast cancer, kidney or bladder carcinoma, gastric cancer, mesothelioma, lung cancer, and ovarian cancer. In HCC, a low XDH mRNA level predicted poorer overall survival, disease-specific survival, disease-free survival, and progression-free survival. The prognostic value of XDH was independent of the clinical features of HCC patients. Indeed, XDH expression in HCC activated several immune-related pathways, including the T cell receptor, PI3K-AKT, and MAPK signaling pathways, which induced a cytotoxic immune response. Importantly, the microenvironment of XDHhigh HCC tumors contained abundant infiltrating CD8 + T cells but not exhausted T cells. A risk prediction signature based on multiple XDH-associated immune genes was revealed as an independent predictor in the TCGA liver cancer cohort.
These findings suggest that XDH is a valuable prognostic biomarker in HCC and other cancers and indicate that it may function in tumor immunology. Loss of XDH expression may be an immune evasion mechanism for HCC.
Hepatocellular carcinoma (HCC) is the seventh most common tumor and the third most common cause of cancer-related death worldwide . Patients with HCC are often diagnosed at advanced stages or are not eligible for curative treatments [2, 3]. Despite optimal treatment, the high recurrence rate of HCC remains a significant concern . The overall prognosis of HCC is poor, with a 5-year overall survival estimated at 10–18% .
The major risk factors for developing HCC include chronic hepatitis B or C virus infections, nonalcoholic steatohepatitis, and alcohol-related liver diseases . Indeed, HCC is strongly influenced by the immune system . The background of chronic inflamed livers leads to tumor development and is associated with an immune-rich contexture in the HCC microenvironment [8,9,10]. Recently, developed immune-based therapies for patients with advanced HCC represent a promising treatment option [5, 9]. However, the complex interface between inflammation, fibrosis, and the immune response involved in HCC pathogenesis is poorly understood .
Xanthine dehydrogenase (XDH) functions as a key regulator of purine metabolism [12, 13], inflammatory cascades , and the innate immune system . In fact, XDH is widely expressed in human tissues, with high levels in the liver [16,17,18,19,20]. Decreased XDH activity is believed to contribute to the development and progression of HCC [21, 22]. Indeed, it has been reported that decreases in XDH activity levels are associated with poor prognoses for cancers, including breast cancer , gastric cancer , ovarian cancer , non-small-cell lung cancer  and colorectal cancer . Despite these intriguing findings, how decreased activity or expression of XDH regulates tumor immunobiology in the development and progression of cancers, including HCC, remains poorly understood.
In this study, we found that the expression of XDH has significant prognostic implications in different types of tumors, including HCC. Moreover, the XDH-associated signaling pathway may regulate the immune response and tumor-infiltrating immune cells in HCC. Importantly, the prognostic value of the XDH-associated immune signature in HCC shed light on XDH might modulate tumor immunity in HCC.
Methods and materials
The Oncomine database was used to determine the expression level of the XDH gene in various types of cancers (https://www.oncomine.org/resource/login.html) . The screening criteria were as follows: P-value of 0.001, fold change of 1.5, and gene ranking of all. Kaplan–Meier plotter was used to assess the effect of 54,675 genes on survival using 10,461 cancer samples. The correlation between XDH expression and survival in breast, ovarian, lung and gastric cancers was analyzed by Kaplan–Meier plotter (http://kmplot.com/analysis/) . Hazard ratios (HRs) with 95% confidence intervals and log-rank P-values were also computed. The Tumor IMmune Estimation Resource (TIMER) is a comprehensive resource for the systematic analysis of immune infiltrates across diverse cancer types (https://cistrome.shinyapps.io/timer/) . TIMER applies a previously published statistical deconvolution method  to infer the abundance of tumor-infiltrating immune cells from gene expression profiles. The TIMER database assesses 10,897 samples from The Cancer Genome Atlas (TCGA) to estimate the abundance of immune infiltrates, including B cells, CD4 + T cells, CD8 + T cells, neutrophils, macrophages, and dendritic cells (DCs) . In addition, correlations between the expression levels of XDH and those of gene markers of tumor-infiltrating immune cells were explored via correlation modules. Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/index.html) is an interactive website that includes data for 9736 tumors and 8587 normal samples from the TCGA and Genotype-Tissue Expression (GTEx) projects and analyzes RNA sequencing expression . GEPIA was used to generate survival curves, including overall survival (OS) and disease-free survival (DFS) curves, based on gene expression with the log-rank test and the Mantel-Cox test in 25 different types of cancer. Gene expression correlation analysis was performed to obtain sets of TCGA expression data. The Spearman method was used to determine the correlation coefficient.
As in our previous report , MHCC-97H cells were transfected with plasmids for overexpressing XDH or a pEZ-Lv201 control vector (Genecopoeia, Guangzhou, China), which served as a negative control. RNA was extracted from the tumors using the Qiagen RNeasy Midi Kit according to the manufacturers' instructions. Complementary RNA (cRNA) was prepared according to the GeneChip Expression Analysis Technical Manual (Affymetrix), hybridized onto Affymetrix Human U133 Plus 2.0, and scanned by a GeneChip Scanner 3000 (Affymetrix). Differential gene expression analysis between treatment groups was calculated using the R limma package . For functional annotation, we used pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) and gene ontology (GO) databases as provided by the clusterprofiler package . P-values were corrected for multiple testing using the Benjamini–Hochberg method. The GSEA algorithm , which tests whether a gene set is significantly perturbed relative to all genes, was applied for analysis.
The preparation of cell lysates and the western blotting procedure were carried out as previously indicated [21, 37]. Briefly, whole-cell lysates were separated by 8–10% sodium dodecyl sulfate–polyacrylamide gel electrophoresis. Equal amounts of resolved proteins were transferred to polyvinylidene difluoride (PVDF) membranes. After incubation with QuickBlock blocking buffer (Beyotime Biotechnology, China, cat. no. P0233), the membranes were then probed with specific primary antibodies (Additional file 1: Table S1) and secondary antibodies (Beyotime Biotechnology, China). Images of the bands were acquired with an Amersham Imager 600 (GE Healthcare, Russellville, AR, USA).
The immunohistochemical staining procedure was conducted as previously described [21, 37, 38]. All patients provided their written informed consent. HCC liver samples were obtained according to a protocol approved by the ethics committee of Fudan University. Immunoreaction images were viewed and captured by Motic DSAssistant software (Motic VM V1 Viewer 2.0). To quantify the density of tumor-infiltrating immune cells, the three most representative areas of stroma were evaluated at 200 × magnification, and the mean value was adopted . For immune cell markers (CD3, CD4, CD8, CD20, CD68 and PD-1) (Additional file 1: Table S2), the counts of all positive cells by immunostaining were recorded in terms of cells/mm2. Two pathologists who were blinded to the patient outcomes independently evaluated the immunohistochemistry (IHC) staining of each sample.
Construction of prognostic signature
To develop an XDH-associated prognostic signature including multiple immune genes, a comprehensive list of immune-related genes was downloaded from the Immunology Database and Analysis Portal (ImmPort) database (https://immport.niaid.nih.gov). The Pearson coefficient between XDH and these immune genes was calculated, and an absolute value of coefficient > 0.3 with P < 0.05 was set as the identification criterion. Then, stepwise variable selection was performed with the Akaike information criterion in the Cox model . After the immune genes were chosen, the prognostic index, referred to as the risk score, was calculated as follows: risk score = β1x1 + β2x2 + … + βixi. In this formula, xi is the expression level of each gene, while βi is the risk coefficient of each gene derived from the Cox model . Multivariate analysis was performed for the risk score with adjustment for age, sex, T stage, N stage, M stage, and tumor-node-metastasis (TNM) stage. Time-dependent receiver operating characteristic (ROC) curves were adopted to determine the prognostic accuracy of the risk score using the timeROC package . Then, a nomogram was built by incorporating the clinical characteristics and the risk scores of HCC patients, and the prognostic value of the nomogram was evaluated . The nomogram was created via the rms package for R software. With the application of the bootstrap method (1000 replicates), a calibration curve was used to visualize the deviation of predicted probabilities from the actual values. The concordance index (C-index) was used to measure the predictive accuracy of the nomogram.
Survival curves were generated by Kaplan–Meier plotter. The results generated in Oncomine are displayed with P-values, fold changes, and ranks. The results of Kaplan–Meier plots and GEPIA are displayed with HR and P or Cox P-values from a log-rank test. The gene expression correlations were evaluated by Spearman’s correlation test and statistical evaluation. P-values < 0.05 were considered to indicate statistical significance. Significant changes are represented as follows: *p < 0.05; **p < 0.01; ***p < 0.001. Nonsignificant changes are labeled as ns.
Decreased expression of XDH mRNA in certain types of human cancer.
To profile XDH mRNA expression in tumor tissues and adjacent normal tissues across several types of cancer, we analyzed the Oncomine database and RNA-seq data of multiple cancers in the TCGA database. In the Oncomine database, lower XDH mRNA expression was detected in bladder cancer, breast cancer, colorectal cancer, leukemia, liver cancer, and lymphoma (Fig. 1A), while higher mRNA expression of XDH was observed in head and neck cancer (HNSC) and lung cancer in some datasets compared to that in normal tissues. The detailed results of XDH expression in different cancer types are summarized in Additional file 1: Table S3. Consistently, in TCGA, low XDH mRNA expression was significantly associated with breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), liver hepatocellular carcinoma (LIHC), and rectum adenocarcinoma (READ), while XDH mRNA expression was higher in HNSC, lung adenocarcinoma (LUAD), and lung squamous cell carcinoma (LUSC) tissues than in adjacent normal tissues (Fig. 1B). Additionally, in TCGA, XDH mRNA expression was lower in kidney chromophobe (KICH) and cholangiocarcinoma (CHOL) tumor tissues and higher in esophageal carcinoma (ESCA) and uterine corpus endometrial carcinoma (UCEC) than in adjacent normal tissues. However, XDH mRNA expression was higher in bladder urothelial carcinoma (BLCA). In summary, a decrease in the level of XDH mRNA was detected in human cancers originating from the liver, bladder, breast, colon, bile duct, kidney, and hematolymphoid system.
XDH mRNA expression is associated with prognosis.
To investigate the impact of XDH mRNA expression on outcomes in patients with cancer, the correlations between XDH expression and OS prognosis and disease progression were evaluated using Kaplan–Meier plotter data based on Affymetrix microarrays. In HCC, patients with a high XDH mRNA expression level had significantly prolonged OS (hazard ratio (HR) 0.55, 95% confidence interval (CI) 0.38 to 0.78, P = 0.00072) and slower disease progression (recurrence-free survival (RFS) HR 0.68, 95% CI 0.49 to 0.95, P = 0.024; progression-free survival (PFS) HR 0.69, 95% CI 0.5 to 0.0.94, P = 0.017; disease-specific survival (DSS) HR 0.53, 95% CI 0.33 to 0.83, P = 0.0045) than patients with a low XDH mRNA expression level (Fig. 2A–D). In addition, high XDH mRNA expression was also correlated with better OS in bladder carcinoma (HR 0.69, 95% CI 0.52 to 0.93, P = 0.015) (Fig. 2E, F) and ovarian cancer (HR 0.87, 95% CI 0.76 to 0.99, P = 0.031) (Fig. 2M, N) and better RFS in breast cancers (HR 0.77, 95% CI 0.69 to 0.86, P = 1.7e−6, Fig. 2G, H). However, high XDH mRNA expression was associated with worse OS in patients with gastric cancer and lung cancer (Fig. 2I–L). Furthermore, the prognostic potential of XDH mRNA expression in 25 types of human cancers was assessed via the GEPIA website (Additional file 1: Figure S1). Indeed, high XDH mRNA expression was associated with a better prognosis in terms of OS in LIHC and BLCA. Moreover, high XDH mRNA expression was correlated with a worse prognosis in terms of OS and DFS in mesothelioma (MESO) and in terms of DFS in LUSC. Together, these findings suggest that the level of expression of XDH mRNA has a significant impact on the prognosis of specific human cancers.
XDH mRNA expression correlates with the clinical characteristics of HCC patients.
We previously reported that decreased XDH mRNA expression is associated with aggressive HCC phenotypes . In humans, XDH is mainly expressed in hepatic tissues . Therefore, it is rational to consider liver cancers as a representative human cancer with weak XDH expression. To better understand the role of XDH expression in cancers, we investigated the relationship between XDH expression and the clinical characteristics of HCC patients in the Kaplan–Meier plotter database (Fig. 3). Indeed, low XDH mRNA expression correlated with poorer OS and PFS in males (OS: HR 0.57, P = 0.0132; PFS HR 0.56, P = 0.027), Asians (OS: HR 0.27, P = 0.0041; PFS HR 0.56, P = 0.016), alcohol consumers (OS: HR 0.23, P = 0.0014; PFS HR 0.53, P = 0.045), nonconsumers of alcohol (OS: HR 0.59, P = 0.024; PFS HR 0.65, P = 0.042), patients without vascular invasion (OS: HR 0.27, P = 0.0041; PFS HR 0.56, P = 0.016) and patients without hepatitis viral infection (OS: HR 0.63, P = 0.05; PFS HR 0.55, P = 0.013). Additionally, low XDH mRNA expression correlated with worse OS and PFS in HCC patients with grade 3 disease (OS: HR 0.39, P = 0.002; PFS: HR 0.55, P = 0.027). However, XDH mRNA expression was not significantly correlated with OS and PFS in stage 1 or stage 2 patients or patients with vascular invasion. Taken together, these results show the prognostic significance of XDH expression in HCC patients based on their clinical characteristics, particularly in those with advanced-stage HCC.
XDH mRNA expression is related to tumor-infiltrating immune cells.
Recently, several clinical trial results showed that immunotherapy is a promising therapy in HCC [5, 45]. In addition, the predictive value of tumor-infiltrating lymphocytes in lymph node metastases and survival has been validated in cancers, including HCC [46, 47], gastric cancer , and cutaneous melanoma . Next, we described the relationship between the expression of XDH and infiltrating immune cells in 40 cancer types, including HCC, using the TIMER database. The results demonstrated that XDH expression was significantly correlated with tumor purity in 18 types of cancer. In addition, XDH expression was associated with infiltrating levels of CD4+ T cells in 9 cancer types, B cells in 12 cancer types, CD8+ T cells in 16 cancer types, DCs in 22 cancer types, macrophages in 9 cancer types, and neutrophils in 19 cancer types (Additional file 1: Table S4). For HCC, we observed that the level of XDH expression was positively correlated with the infiltration levels of CD8+ T cells (r = 0.157, p=3.46e−03) and macrophages (r = 0.221, p = 3.40e−05) but negatively correlated with the infiltration levels of B cells (r = − 0.178, p = 8.68e−04) and myeloid DCs (r = − 0.179, p=8.66e−04; Fig. 4A). As there is no estimate of exhausted T cells in the TIMER algorithm, we selected several marker genes representing tumor-infiltrating exhausted T cells . The expression of XDH negatively correlated with that of programmed cell death protein 1 (PD1) (r = − 0.174, p = 1.18e−03) and cytotoxic T-lymphocyte-associated protein 4 (CTLA4) (r = − 0.212, p = 6.96e−05). Furthermore, a negative trend was also observed in T cell immunoglobulin and mucin domain-containing protein 3 (TIM3) and lymphocyte-activation gene 3 (LAG-3) (Fig. 4B). Together, these findings suggest that XDH may modulate the infiltration of immune cells into tumor tissues.
XDH correlates with the immune response in HCC.
To understand the mechanistic role of XDH, we performed in-depth analysis of the relationship between XDH expression and tumor-infiltrating immune cells. To profile XDH-related gene expression, XDH was overexpressed in an HCC cell line (MHCC-97H), and the cells were subjected to RNA array analysis. Differential gene expression analysis (overexpression group versus control group) was performed with the limma package. GO enrichment analysis showed that the immune response was upregulated by XDH (Fig. 5A). GSEA confirmed that XDH overexpression upregulated the purine metabolism pathway and activated several immune-related pathways, including the T cell receptor, PI3K-AKT, and MAPK signaling pathways (Fig. 5B, E). Indeed, western blot analysis confirmed that XDH overexpression induced the activation of the PI3K-Akt pathway (Fig. 5F; Additional file 1: Figure S2; Additional file 3). To clarify the correlation between XDH expression and immune cell markers in HCC tumor tissue, IHC analysis of 6 types of immune cells was performed (Fig. 5; Additional file 1: Figure S3; Additional file 2). As shown in Fig. 5I, the expression of XDH was positively correlated with the infiltration of CD8 + immune cells (r = 0.6354, p = 0.0025). Furthermore, there was a trend of a negative correlation between XDH and PD1 + immune cell infiltration (Fig. 5L). These findings suggest that the expression of XDH may trigger a cytotoxic immune response in HCC.
The prognostic implication of the XDH-associated immune signature in HCC
To further study the prognostic value of XDH, 218 immune-related genes were identified to be significantly associated with XDH expression by Pearson correlation analysis. The protein–protein network generated with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) online server showed a tight correlation (Additional file 1: Figure S3). Then, a stepwise multivariable Cox regression analysis was used to identify the prognostic value of the XDH-associated immune signature in HCC. Consequently, an optimal prognostic signature of 20 genes related to XDH in HCC was revealed. The biological functions of these genes are presented in Additional file 1: Table S5. The distribution of risk scores, survival status, and signature gene expression profiles for HCC were visualized (Fig. 6A). Compared with patients with high risk scores, patients with low risk scores had a significantly longer survival time, as indicated by Kaplan–Meier survival curves (log-rank test, P < 0.001; Fig. 6B). In addition, the risk score showed a strong discriminative ability for 3-year and 5-year OS (Fig. 6C). After adjusting for sex, age, and stage, multivariable Cox regression analysis showed that the risk score was an independent predictor of prognosis in HCC (HR = 1.35, 95% CI = 1.26–1.4, p < 0.001; Fig. 6D). Finally, we built a prognostic nomogram for HCC patients by determining weighted coefficients for risk score, stage, age, and sex. The calibration curves showed that the nomogram-predicted probability (solid line) well matched the ideal reference line for 3- and 5-year survival (Fig. 6E, F). In addition, the prognostic nomogram displayed good discrimination with a C-index of 0.73. These data show that the XDH-associated immune signature has probable prognostic value in HCC.
The recent success of immune-based therapies has revolutionized the HCC treatment armamentarium [51,52,53,54]. However, there are very few useful biomarkers available for the identification of sensitivity and resistance to checkpoint inhibitors and their combinations in HCC . In this study, we found that XDH mRNA levels correlated with the prognosis of several human cancers. In HCC, the downregulation of XDH was an independent survival predictor associated with worse prognosis. The XDH-related cell signaling pathway was associated with a cytotoxic immune response in HCC. Furthermore, XDH mRNA levels correlated with the numbers of tumor-infiltrated immune cells based on the levels of markers for different immune cell types in HCC. The prognostic implication of the XDH-associated immune signature in HCC supports that XDH is a potential immune biomarker for HCC.
In the present study, we first comprehensively analyzed the mRNA levels of XDH and their prognostic value in cancer tissues using multiple databases. We also found that XDH expression was significantly downregulated in certain types of cancer, while increased XDH expression was detected in some other cancers, such as HNSC, LUAD, and LUSC. Indeed, the XDH expression data were consistent in HCC tissues across different databases. Thus, the variability of XDH expression in different types of cancers may reflect differences in the underlying molecular and genetic mechanisms for cancer development and progression. Although the prognostic potential of XDH in several types of cancer has been reported, our work expands this knowledge of XDH in other cancers. Indeed, analysis of the GEPIA database revealed that low XDH expression was correlated with a worse prognosis in cancer types such as CHOL, lower grade glioma (LGG), LUAD, and thyroid carcinoma (THCA).
We hypothesized that XDH could play a significant role in regulating tumor immunology and therefore influence the outcomes of cancer patients. Analysis of the TIMER database revealed that XDH expression correlated with the infiltration status of immune cells in several cancer types, including HCC. In HCC, a positive correlation was found between XDH expression and CD8 + T-cell infiltration in the analysis of both the TIMER database and IHC staining of immune cell markers in HCC tumor tissues. Although there was a significant positive correlation between XDH expression and macrophage infiltration in the TIMER database, only a positive trend was found between XDH expression and CD68 + macrophage infiltration in this small group of HCC tissues. Indeed, Saidak et al. indicated that XDH expression is linked to an immune infiltrate in tumors . XDH expression initiates several immune-related pathways, including the T cell receptor, PI3K-AKT, and MAPK signaling pathways. Similarly, both Wang et al. and Peng et al. found that immune-related pathways in HCC were mainly involved in the MAPK signaling pathway and PI3K-AKT signaling pathway [56, 57]. Moreover, the expression of the exhausted T cell markers PD-1 and CTLA-4, which are critical inhibitory immune checkpoint proteins, negatively correlated with XDH expression. Thus, decreased expression of XDH may facilitate tumor invasion in HCC and possibly in other similar cancers with low XDH expression. However, this hypothesis warrants further investigation.
The development of well-verified signatures for cancer prognosis evaluation represents a critical step for the implementation of stratification strategies and personalized immunotherapy for cancer [58,59,60,61]. In this study, we constructed a nomogram for personalized prognosis prediction in HCC with a C-index of 0.73. The risk scores derived from the XDH-associated immune gene signatures were significantly associated with survival. Notably, most of the immune genes integrated into the prognostic signatures participate in the regulation of the activity of immune cells, highlighting the significance of cytotoxic activity in HCC. Indeed, several prognosis-related risk models have been established according to the TCGA liver cancer cohort [62,63,64,65,66,67]. For instance, among 6 developed models, the 10-gene model proposed by Zhao et al. reached the highest C-index of 0.715 . Similarly, Xu et al. established an 8-immune-gene prognostic signature with a C-index of 0.725 for HCC . Recently, Chen et al. found nine immune-related gene pairs (IRGPs) that could be used to determine the outcomes of HCC patients with a C-index of 0.755 by integrating three public datasets of HCC 67. Moreover, some other predication models with high C-index values have been reported . Nevertheless, none of those studies had risk models that reached area under the curve (AUC) values as high as 0.81 and 0.802 for 3-year and 5-year OS, respectively, as indicated by our model.
Our study has some limitations to be improved. First, this study is mainly based on data retrieved from public databases. Some in vitro or in vivo experiments are necessary to validate the immune evasion mechanism by which XDH contributes to the progression of HCC. Second, the sample size for some specific tumor types including HCC was too small to form a solid conclusion. Hence, further studies are necessary to verify the role of XDH in the regulation of immune-related pathways in HCC.
In summary, our results suggest that XDH is a potential independent prognostic biomarker for HCC. The XDH-associated cell signaling pathway may affect immune cell infiltration into the tumor microenvironment. In HCC, decreased XDH expression correlates with a reduced cytotoxic immune response, and our prognostic XDH-associated immune signature provides a valuable tool for precision therapy. These data may indicate that XDH plays a role in HCC tumor immunology via an immune evasion mechanism.
Tumor immune estimation resource
The Cancer Genome Atlas
Gene expression profiling interactive analysis
Kyoto Encyclopedia of Genes and Genomes
Receiver operating characteristic
Programmed cell death protein 1
Cytotoxic T-lymphocyte-associated protein 4
T cell immunoglobulin and mucin domain-containing protein 3
Lymphocyte-activation gene 3
Search Tool for the Retrieval of Interacting Genes/Proteins
Area under the curve
Bladder urothelial carcinoma
Breast invasive carcinoma
Cervical squamous cell carcinoma
Diffuse large b-cell lymphoma
Head and neck squamous cell carcinoma
Kidney renal clear cell carcinoma
Kidney renal papillary cell carcinoma
Liver hepatocellular carcinoma
Lung squamous cell carcinoma
Lower grade glioma
Ovarian serous cystadenocarcinoma
Pheochromocytoma and paraganglioma
Skin cutaneous melanoma
Testicular germ cell tumor
Uterine corpus endometrial carcinoma
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021. https://doi.org/10.3322/caac.21660.
Huang J, Yan L, Cheng Z, Wu H, Du L, Wang J, et al. A randomized trial comparing radiofrequency ablation and surgical resection for HCC conforming to the Milan criteria. Ann Surg. 2010;252(6):903–12. https://doi.org/10.1097/SLA.0b013e3181efc656.
Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, et al. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med. 2008;359(4):378–90. https://doi.org/10.1056/NEJMoa0708857.
El-Serag HB, Rudolph KL. Hepatocellular carcinoma: epidemiology and molecular carcinogenesis. Gastroenterology. 2007;132(7):2557–76. https://doi.org/10.1053/j.gastro.2007.04.061.
Qin S, Ren Z, Meng Z, Chen Z, Chai X, Xiong J, et al. Camrelizumab in patients with previously treated advanced hepatocellular carcinoma: a multicentre, open-label, parallel-group, randomised, phase 2 trial. Lancet Oncol. 2020;21(4):571–80. https://doi.org/10.1016/S1470-2045(20)30011-5.
Makarova-Rusher OV, Altekruse SF, McNeel TS, Ulahannan S, Duffy AG, Graubard BI, et al. Population attributable fractions of risk factors for hepatocellular carcinoma in the United States. Cancer. 2016;122(11):1757–65. https://doi.org/10.1002/cncr.29971.
Hou J, Zhang H, Sun B, Karin M. The immunobiology of hepatocellular carcinoma in humans and mice: Basic concepts and therapeutic implications. J Hepatol. 2020;72(1):167–82. https://doi.org/10.1016/j.jhep.2019.08.014.
Buonaguro L, Mauriello A, Cavalluzzo B, Petrizzo A, Tagliamonte M. Immunotherapy in hepatocellular carcinoma. Ann Hepatol. 2019;18(2):291–7. https://doi.org/10.1016/j.aohep.2019.04.003.
Pinato DJ, Guerra N, Fessas P, Murphy R, Mineo T, Mauri FA, et al. Immune-based therapies for hepatocellular carcinoma. Oncogene. 2020;39(18):3620–37. https://doi.org/10.1038/s41388-020-1249-9.
Johnston MP, Khakoo SI. Immunotherapy for hepatocellular carcinoma: current and future. World J Gastroenterol. 2019;25(24):2977–89. https://doi.org/10.3748/wjg.v25.i24.2977.
Keenan BP, Fong L, Kelley RK. Immunotherapy in hepatocellular carcinoma: the complex interface between inflammation, fibrosis, and the immune response. J Immunother Cancer. 2019;7(1):267. https://doi.org/10.1186/s40425-019-0749-z.
Battelli MG, Bolognesi A, Polito L. Pathophysiology of circulating xanthine oxidoreductase: new emerging roles for a multi-tasking enzyme. Biochim Biophys Acta. 2014;1842(9):1502–17. https://doi.org/10.1016/j.bbadis.2014.05.022.
Battelli MG, Polito L, Bortolotti M, Bolognesi A. Xanthine oxidoreductase in cancer: more than a differentiation marker. Cancer Med. 2016;5(3):546–57. https://doi.org/10.1002/cam4.601.
Ohtsubo T, Rovira II, Starost MF, Liu C, Finkel T. Xanthine oxidoreductase is an endogenous regulator of cyclooxygenase-2. Circ Res. 2004;95(11):1118–24. https://doi.org/10.1161/01.RES.0000149571.96304.36.
Ives A, Nomura J, Martinon F, Roger T, LeRoy D, Miner JN, et al. Xanthine oxidoreductase regulates macrophage IL1beta secretion upon NLRP3 inflammasome activation. Nat Commun. 2015;6:6555. https://doi.org/10.1038/ncomms7555.
Xu P, Huecksteadt TP, Harrison R, Hoidal JR. Molecular cloning, tissue expression of human xanthine dehydrogenase. Biochem Biophys Res Commun. 1994;199(2):998–1004. https://doi.org/10.1006/bbrc.1994.1328.
Wright RM, Vaitaitis GM, Wilson CM, Repine TB, Terada LS, Repine JE. cDNA cloning, characterization, and tissue-specific expression of human xanthine dehydrogenase/xanthine oxidase. Proc Natl Acad Sci USA. 1993;90(22):10690–4. https://doi.org/10.1073/pnas.90.22.10690.
Sato A, Nishino T, Noda K, Amaya Y, Nishino T. The structure of chicken liver xanthine dehydrogenase. cDNA cloning and the domain structure. J Biol Chem. 1995;270(6):2818–26. https://doi.org/10.1074/jbc.270.6.2818.
Harmon DB, Mandler WK, Sipula IJ, Dedousis N, Lewis SE, Eckels JT, et al. Hepatocyte-specific ablation or whole-body inhibition of xanthine oxidoreductase in mice corrects obesity-induced systemic hyperuricemia without improving metabolic abnormalities. Diabetes. 2019;68(6):1221–9. https://doi.org/10.2337/db18-1198.
Srivastava M, Kale RK. Radiomodfication of xanthine oxidoreductase system in the liver of mice by phenylmethylsulfonyl fluoride and dithiothreitol. Radiat Res. 2000;154(1):94–103. https://doi.org/10.1667/0033-7587(2000)154[0094:roxosi]2.0.co;2.
Chen GL, Ye T, Chen HL, Zhao ZY, Tang WQ, Wang LS, et al. Xanthine dehydrogenase downregulation promotes TGFbeta signaling and cancer stem cell-related gene expression in hepatocellular carcinoma. Oncogenesis. 2017;6(9): e382. https://doi.org/10.1038/oncsis.2017.81.
Sun Q, Zhang Z, Lu Y, Liu Q, Xu X, Xu J, et al. Loss of xanthine oxidoreductase potentiates propagation of hepatocellular carcinoma stem cells. Hepatology. 2020;71(6):2033–49. https://doi.org/10.1002/hep.30978.
Linder N, Lundin J, Isola J, Lundin M, Raivio KO, Joensuu H. Down-regulated xanthine oxidoreductase is a feature of aggressive breast cancer. Clin Cancer Res. 2005;11(12):4372–81. https://doi.org/10.1158/1078-0432.CCR-04-2280.
Linder N, Haglund C, Lundin M, Nordling S, Ristimaki A, Kokkola A, et al. Decreased xanthine oxidoreductase is a predictor of poor prognosis in early-stage gastric cancer. J Clin Pathol. 2006;59(9):965–71. https://doi.org/10.1136/jcp.2005.032524.
Linder N, Butzow R, Lassus H, Lundin M, Lundin J. Decreased xanthine oxidoreductase (XOR) is associated with a worse prognosis in patients with serous ovarian carcinoma. Gynecol Oncol. 2012;124(2):311–8. https://doi.org/10.1016/j.ygyno.2011.10.026.
Kim AW, Batus M, Myint R, Fidler MJ, Basu S, Bonomi P, et al. Prognostic value of xanthine oxidoreductase expression in patients with non-small cell lung cancer. Lung Cancer. 2011;71(2):186–90. https://doi.org/10.1016/j.lungcan.2010.05.006.
Linder N, Martelin E, Lundin M, Louhimo J, Nordling S, Haglund C, et al. Xanthine oxidoreductase—clinical significance in colorectal cancer and in vitro expression of the protein in human colon cancer cells. Eur J Cancer. 2009;45(4):648–55. https://doi.org/10.1016/j.ejca.2008.10.036.
Rhodes DR, Kalyana-Sundaram S, Mahavisno V, Varambally R, Yu J, Briggs BB, et al. Oncomine 3.0: genes, pathways, and networks in a collection of 18,000 cancer gene expression profiles. Neoplasia. 2007;9(2):166–80. https://doi.org/10.1593/neo.07112.
Nagy A, Munkacsy G, Gyorffy B. Pancancer survival analysis of cancer hallmark genes. Sci Rep. 2021;11(1):6047. https://doi.org/10.1038/s41598-021-84787-5.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–10. https://doi.org/10.1158/0008-5472.CAN-17-0307.
Li B, Severson E, Pignon JC, Zhao H, Li T, Novak J, et al. Comprehensive analyses of tumor immunity: implications for cancer immunotherapy. Genome Biol. 2016;17(1):174. https://doi.org/10.1186/s13059-016-1028-7.
Aran D, Sirota M, Butte AJ. Systematic pan-cancer analysis of tumour purity. Nat Commun. 2015;6:8971. https://doi.org/10.1038/ncomms9971.
Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45(W1):W98-102. https://doi.org/10.1093/nar/gkx247.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47. https://doi.org/10.1093/nar/gkv007.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102(43):15545–50. https://doi.org/10.1073/pnas.0506580102.
Chen GL, Li R, Chen XX, Wang J, Cao S, Song R, et al. Fra-2/AP-1 regulates melanoma cell metastasis by downregulating Fam212b. Cell Death Differ. 2021;28(4):1364–78. https://doi.org/10.1038/s41418-020-00660-4.
Chen GL, Luo Y, Eriksson D, Meng X, Qian C, Bauerle T, et al. High fat diet increases melanoma cell growth in the bone marrow by inducing osteopontin and interleukin 6. Oncotarget. 2016;7(18):26653–69. https://doi.org/10.18632/oncotarget.8474.
Wei L, Delin Z, Kefei Y, Hong W, Jiwei H, Yange Z. A classification based on tumor budding and immune score for patients with hepatocellular carcinoma. Oncoimmunology. 2020;9(1):1672495. https://doi.org/10.1080/2162402X.2019.1672495.
Choi I, Wells BJ, Yu C, Kattan MW. An empirical approach to model selection through validation for censored survival data. J Biomed Inform. 2011;44(4):595–606. https://doi.org/10.1016/j.jbi.2011.02.005.
Aguirre-Gamboa R, Gomez-Rueda H, Martinez-Ledesma E, Martinez-Torteya A, Chacolla-Huaringa R, Rodriguez-Barrientos A, et al. SurvExpress: an online biomarker validation tool and database for cancer gene expression data using survival analysis. PLoS ONE. 2013;8(9): e74250. https://doi.org/10.1371/journal.pone.0074250.
Saha-Chaudhuri P, Heagerty PJ. Non-parametric estimation of a time-dependent predictive accuracy curve. Biostatistics. 2013;14(1):42–59. https://doi.org/10.1093/biostatistics/kxs021.
Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. 2008;26(8):1364–70. https://doi.org/10.1200/JCO.2007.12.9791.
Stirpe F, Ravaioli M, Battelli MG, Musiani S, Grazi GL. Xanthine oxidoreductase activity in human liver disease. Am J Gastroenterol. 2002;97(8):2079–85. https://doi.org/10.1111/j.1572-0241.2002.05925.x.
Granier C, De Guillebon E, Blanc C, Roussel H, Badoual C, Colin E, et al. Mechanisms of action and rationale for the use of checkpoint inhibitors in cancer. ESMO Open. 2017;2(2): e213. https://doi.org/10.1136/esmoopen-2017-000213.
Yao W, He JC, Yang Y, Wang JM, Qian YW, Yang T, et al. The prognostic value of tumor-infiltrating lymphocytes in hepatocellular carcinoma: a systematic review and meta-analysis. Sci Rep. 2017;7(1):7525. https://doi.org/10.1038/s41598-017-08128-1.
Tian MX, Liu WR, Wang H, Zhou YF, Jin L, Jiang XF, et al. Tissue-infiltrating lymphocytes signature predicts survival in patients with early/intermediate stage hepatocellular carcinoma. BMC Med. 2019;17(1):106. https://doi.org/10.1186/s12916-019-1341-6.
Japanese Gastric Cancer Association. Japanese gastric cancer treatment guidelines 2014 (ver. 4). Gastric Cancer. 2017;20(1):1–19. https://doi.org/10.1007/s10120-016-0622-4.
Azimi F, Scolyer RA, Rumcheva P, Moncrieff M, Murali R, McCarthy SW, et al. Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma. J Clin Oncol. 2012;30(21):2678–83. https://doi.org/10.1200/JCO.2011.37.8539.
Jiang Y, Li Y, Zhu B. T-cell exhaustion in the tumor microenvironment. Cell Death Dis. 2015;6: e1792. https://doi.org/10.1038/cddis.2015.162.
Sangro B, Sarobe P, Hervas-Stubbs S, Melero I. Advances in immunotherapy for hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2021. https://doi.org/10.1038/s41575-021-00438-0.
Llovet JM, Montal R, Sia D, Finn RS. Molecular therapies and precision medicine for hepatocellular carcinoma. Nat Rev Clin Oncol. 2018;15(10):599–616. https://doi.org/10.1038/s41571-018-0073-4.
Nakano S, Eso Y, Okada H, Takai A, Takahashi K, Seno H. Recent advances in immunotherapy for hepatocellular carcinoma. Cancers. 2020. https://doi.org/10.3390/cancers12040775.
Pinter M, Scheiner B, Peck-Radosavljevic M. Immunotherapy for advanced hepatocellular carcinoma: a focus on special subgroups. Gut. 2021;70(1):204–14. https://doi.org/10.1136/gutjnl-2020-321702.
Saidak Z, Louandre C, Dahmani S, Sauzay C, Guedda S, Chauffert B, et al. A pan-cancer study of the transcriptional regulation of uricogenesis in human tumours: pathological and pharmacological correlates. Biosci Rep. 2018. https://doi.org/10.1042/BSR20171716.
Peng Y, Liu C, Li M, Li W, Zhang M, Jiang X, et al. Identification of a prognostic and therapeutic immune signature associated with hepatocellular carcinoma. Cancer Cell Int. 2021;21(1):98. https://doi.org/10.1186/s12935-021-01792-4.
Wang Z, Zhu J, Liu Y, Liu C, Wang W, Chen F, et al. Development and validation of a novel immune-related prognostic model in hepatocellular carcinoma. J Transl Med. 2020;18(1):67. https://doi.org/10.1186/s12967-020-02255-6.
Xie P, Ma Y, Yu S, An R, He J, Zhang H. Development of an immune-related prognostic signature in breast cancer. Front Genet. 2019;10:1390. https://doi.org/10.3389/fgene.2019.01390.
Hsu DS, Kim MK, Balakumaran BS, Acharya CR, Anders CK, Clay T, et al. Immune signatures predict prognosis in localized cancer. Cancer Invest. 2010;28(7):765–73. https://doi.org/10.3109/07357900903095755.
Zhou C, Li C, Yan F, Zheng Y. Identification of an immune gene signature for predicting the prognosis of patients with uterine corpus endometrial carcinoma. Cancer Cell Int. 2020;20(1):541. https://doi.org/10.1186/s12935-020-01560-w.
Song Q, Shang J, Yang Z, Zhang L, Zhang C, Chen J, et al. Identification of an immune signature predicting prognosis risk of patients in lung adenocarcinoma. J Transl Med. 2019;17(1):70. https://doi.org/10.1186/s12967-019-1824-4.
Zhao K, Xu L, Li F, Ao J, Jiang G, Shi R, et al. Identification of hepatocellular carcinoma prognostic markers based on 10-immune gene signature. 2020. Biosci Rep.
Qiao GJ, Chen L, Wu JC, Li ZR. Identification of an eight-gene signature for survival prediction for patients with hepatocellular carcinoma based on integrated bioinformatics analysis. PeerJ. 2019;7: e6548. https://doi.org/10.7717/peerj.6548.
Liu GM, Zeng HD, Zhang CY, Xu JW. Identification of a six-gene signature predicting overall survival for hepatocellular carcinoma. Cancer Cell Int. 2019;19:138. https://doi.org/10.1186/s12935-019-0858-2.
Wang Z, Teng D, Li Y, Hu Z, Liu L, Zheng H. A six-gene-based prognostic signature for hepatocellular carcinoma overall survival prediction. Life Sci. 2018;203:83–91. https://doi.org/10.1016/j.lfs.2018.04.025.
Zheng Y, Liu Y, Zhao S, Zheng Z, Shen C, An L, et al. Large-scale analysis reveals a novel risk score to predict overall survival in hepatocellular carcinoma. Cancer Manag Res. 2018;10:6079–96. https://doi.org/10.2147/CMAR.S181396.
Chen W, Ou M, Tang D, Dai Y, Du W. Identification and validation of immune-related gene prognostic signature for hepatocellular carcinoma. J Immunol Res. 2020;2020:5494858. https://doi.org/10.1155/2020/5494858.
Xu D, Wang Y, Zhou K, Wu J, Zhang Z, Zhang J, et al. Development and validation of a novel 8 immune gene prognostic signature based on the immune expression profile for hepatocellular carcinoma. Onco Targets Ther. 2020;13:8125–40. https://doi.org/10.2147/OTT.S263047.
The authors would like to thank the TCGA curators for the availability of the data.
This study was supported by the National Natural Science Foundation of China (No. 81802362). The funding agencies had no role in the design and conduct of the study; in the collection, analysis, and interpretation of the data; or in the preparation, review or approval of the manuscript.
Ethics approval and consent to participate
All patients provided their written informed consent. HCC liver samples were obtained according to a protocol approved by the ethics committee of Fudan University, Shanghai Cancer Center (No. 050432–4-1212B).
Consent for publication
Consent to publish was obtained from the participants.
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Correlation of xanthine dehydrogenase (XDH) expression with prognosis in 25 diverse types of cancer. Overall survival (OS) and disease-free survival (DFS) curves comparing the high and low XDH expression groups in the BLCA (A-B), BRCA (C-D), CESC (E-F), CHOL (G-H), COAD (I-J), ESCA (K-L), GBM (M-N), HNSC (O-P), KIRC (Q-R), KIRP (S-T), LIHC (U-V), LUSC (W-X), MESO (Y-Z), OV (AA-AB), PAAD (AC-AD), PRAD (AE-AF), READ (AG-AH), SARC (AI-AJ), SKCM (AK-AL), STAD (AM-AN), TGCT (AO-AP), THCA (AQ-AR), THYM (AS-AT), UCEC (AU-AV), and UCS (AW-AX) cohorts. Figure S2. Quantification of western blotting results in Fig. 5F. Figure S3. Representative images of CD3 (A), CD4 (B), CD8 (C), CD20 (D), CD68 (E), PD1 (F) and XDH (G) staining. Figure S4. Protein-protein network of 128 XDH-associated immune genes in LIHC produced by the STRING online server. Table S1. List of the primary antibodies and dilutions used in the study. Table S2. Antibody sources and staining conditions. Table S3. XDH expression in cancer tissues versus normal tissues in the bladder cancer, breast cancer, colorectal cancer, head and neck cancer, leukemia, liver cancer, lung cancer, and lymphoma datasets in the Oncomine database. Table S4. Correlation of XDH expression with immune infiltration levels in diverse types of cancer in the Tumor IMmune Estimation Resource (TIMER) database. Table S5. Functions of the genes included in the prognostic signature.
About this article
Cite this article
Lin, Z., Xie, YZ., Zhao, MC. et al. Xanthine dehydrogenase as a prognostic biomarker related to tumor immunology in hepatocellular carcinoma. Cancer Cell Int 21, 475 (2021). https://doi.org/10.1186/s12935-021-02173-7
- Hepatocellular carcinoma
- Exhaustion marker