The pan-cancer landscape of aldo-keto reductase1B10 reveals that its expression is diminished in gastric cancer

Article type: Research Article

Keywords: AKR1B10, pan-cancer, progression, tumor microenvironment, immunity

Authors: Anqi Wu, Hao Li, Mengnan Gao, Juan Liang, Jiaqi Huang, Jaume Farrés, Deliang Cao, Guoqing Li

Affiliations: Department of Clinical Research Center, The Second Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, China; Hunan Province Key Laboratory of Basic and Clinical Pharmacological Research on Gastrointestinal Tumors, The Second Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, China; Department of Pathology, The Second Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, China; Department of Gastroenterology, The Second Affiliated Hospital, Hengyang Medical School, University of South China, Hengyang, China; Department of Biochemistry and Molecular Biology, Universitat Autònoma de Barcelona, Barcelona, Spain; Hunan Province Key Laboratory of Cancer Cellular and Molecular Pathology, Hengyang Medical School, Cancer Research Institute, University of South China, Hengyang, China

License: Copyright © 2024 Wu, Li, Gao, Liang, Huang, Farrés, Cao and Li CC BY 4.0 This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Article links: DOI: 10.3389/fimmu.2024.1488042  |  PubMed: 39712017  |  PMC: PMC11659136

Relevance: Moderate: mentioned 3+ times in text

Full text: PDF (12.6 MB)

Introduction

Aldo-keto reductase 1B10 (AKR1B10), also known as aldose reductase-like 1 (ARL-1), is a member of aldo/keto reductase (AKR) superfamily. AKR1B10 along with AKR1B1 and AKR1B15 together belong to the human AKR1B subfamily. These genes cluster at chromosome 7q33, and share more than 60% amino acid sequence identity (ref. 1). AKR1Bs catalyze the reduction of aromatic and aliphatic aldehydes to alcoholic forms with NADPH as a co-enzyme (ref. 1–ref. 5). These enzymes show certain substrates overlaps, but with clear variation of the catalytic efficiencies. For instance, AKR1B10 specifically catalyzes retinal (all-trans-retinaldehyde) with high efficiency to retinol (ref. 3), but its glucose reductase activity (ref. 1–ref. 3), prostaglandin F synthase activity (ref. 6), 17β-hydroxysteroid dehydrogenase activity for estrone and 4-androstene-3,17-dione (ref. 4, ref. 5) are relatively limited. AKR1B10 also efficiently reduces endogenous carbonyl compounds, including lipid peroxidation-derived cytotoxic aldehydes (ref. 2), isoprenyl aldehydes (ref. 7), and xenobiotics, such as drugs (ref. 8) and polycyclic aromatic hydrocarbon derivatives (ref. 9). The active site of AKR1B10 is the α/β-barrel tertiary structure which is consisted of 8 α-helices and β-sheets (ref. 10).

Under physiological conditions, AKR1B10 is constitutively abundant in small intestine and colon (ref. 1, ref. 11), which can protect mucosa epithelial cells from lesions of the reactive aldehydes in digested food through reductive detoxification. Pathologically, abnormal AKR1B10 expression is associated with various skin diseases (ref. 12–ref. 14), precancerous diseases (ref. 15–ref. 22) and cancerous diseases. However, the role of AKR1B10 in cancer development and progression is complex. AKR1B10 expression is elevated and correlated with poor prognosis in several types of cancers like hepatocellular carcinoma (ref. 23), breast carcinoma (ref. 24–ref. 26), lung cancer (ref. 27–ref. 29) and oral cancer (ref. 30–ref. 32), but is quite opposite in colorectal cancer (ref. 33–ref. 37). AKR1B10 may promote cancer progression through fatty acid synthesis, but prevent cancer via regulating autophagy and FGF1. Therefore, a special function of AKR1B10 in an indicated type of cancer may be context-dependent. A comprehensive analysis of AKR1B10 in pan-cancers is certainly warranted.

In this study, we compared the expression of AKR1B10 in multiple types of normal tissues and cancers, evaluated the relationship of AKR1B10 expression with DNA mutation, DNA methylation, cancer clinicopathology, cancer microenvironment and immune cell infiltration, and survivals in pan-cancers. Our study for the first time presents a landscape of AKR1B10 in pan-cancers and discusses the context-dependent function of AKR1B10 in cancers.

Material and methods

Data collection

All of the data analyzed in this article were public available. The RNA sequencing data of health tissues was downloaded from Genotype-Tissue Expression (GTEx) portal (https://www.gtexportal.org/) which collects organs from autopsies. The RNA sequencing data, simple nucleotide variations, copy number variations, as well as clinical data of TCGA projects were downloaded from GDC portal (https://portal.gdc.cancer.gov/), which contains datasets of adrenocortical carcinoma (ACC), bladder urothelial carcinoma (BLCA), breast invasive carcinoma (BRCA), cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), cholangiocarcinoma (CHOL), colon adenocarcinoma (COAD), lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), esophageal carcinoma (ESCA), glioblastoma multiforme (GBM), head and neck squamous cell carcinoma (HNSC), kidney chromophobe (KICH), kidney renal papillary cell carcinoma (KIRP), kidney renal clear cell carcinoma (KIRC), acute myeloid leukemia (LAML), brain lower grade glioma (LGG), liver hepatocellular carcinoma (LIHC), lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), mesothelioma (MESO), ovarian serous cystadenocarcinoma (OV), pancreatic adenocarcinoma (PAAD), pheochromocytoma and paraganglioma (PCPG), prostate adenocarcinoma (PRAD), rectum adenocarcinoma (READ), sarcoma (SARC), stomach adenocarcinoma (STAD), skin cutaneous melanoma (SKCM), testicular germ cell tumors (TGCT), thyroid carcinoma (THCA), thymoma (THYM), uterine corpus endometrial carcinoma (UCEC), uterine carcinosarcoma (UCS), uveal melanoma (UVM).

The immunohistochemistry (IHC) results of AKR1B10 in cancers and normal tissues were downloaded from the Human Protein Atlas (HPA) database (https://www.proteinatlas.org/). STAD enrolled samples from stomach (T-6300), stomach upper (T-62350) and stomach lower (T-63700) Adenocarcinoma; COAD enrolled samples from colon (T-67000) adenocarcinoma; LIHC enrolled samples from liver (T-56000) hepatocellular carcinoma; PAAD enrolled pancreas (T-59000) adenocarcinoma; UCEC enrolled samples from endometrium (T-84000) adenocarcinoma; CESC enrolled samples from cervix (T-83000) squamous cell carcinoma.

The validation datasets (GSE107850, GSE7696, GSE41272, GSE8894 and GSE57303) were available in gene expression omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/).

Analysis of AKR1B10 expression and clinical features

The association of AKR1B10 expression with stages or grades in each cancer from TCGA database was analyzed by TISIDB (http://cis.hku.hk/TISIDB/), which is an integrated repository portal for tumor-immune system interactions.

The prognosis information including overall survival (OS), progression-free interval (PFI) and disease-free interval (DFI) from clinical follow-up data were integrated with AKR1B10 expression from RNA sequencing data. The cutoff value of AKR1B10 expression was determined by the optimal calculation result. The prognostic value of AKR1B10 was determined by Kaplan-Meier analysis.

Mutation analysis

The frequency and site of mutations (missense mutations, nonsense mutations, frame shift inserts, frame shift deletions and splice sites) in each cancer based on the simple nucleotide variation dataset processed by MuTect2 software (ref. 38) were integrated with AKR1B10 protein domain information acquired from maftools (version 2.2.10) in R software. The tumor mutation burden (TMB)scores were calculated by tmb function of maftool.

The copy number variation (CNV) datasets processed by GISTIC software (ref. 39) were integrated with log2(x+0.001) transformed AKR1B10 expression from RNA-seq datasets. A total of 26 types of cancers were analyzed with R software (version 3.6.4) when cancers with sample number less than 3 were excluded.

Visualization of methylation data and single cell RNA sequencing data

The correlation between methylation of AKR1B10 promoter and FPKM in each cancer was analyzed by DNMIVD (ref. 40) which is an interactive visualization database for DNA methylation profiles of diverse human cancers including TCGA databases.

The expression of AKR1B10 among different organs and cell types from single cell RNA sequencing were interrogated and visualized by Human Cell Landscape (http://bjs.zju.edu.cn/HCL/) (ref. 41), which contains 60 human tissue samples from 8 systems.

Microenvironment analysis

The stromal, immune and ESTIMATE scores of each cancer sample were calculated by ESTIMATE package (version 1.0.13) in R software based on the gene expression matrix after data normalization and transformation mentioned above. While the purity and homologous recombination deficiency (HRD) data of each cancer were downloaded from previous published literature (ref. 42). The expression of AKR1B10 was integrated with purity, stromal scores, immune scores, ESTIMATE scores, and HRD data, and then the correlation was analyzed with corr.test function of psych package (version 2.1.6) in R software.

Immune cell infiltration analysis

Normalized pan-cancer data were annotated with gene symbols, and the expression was transformed with log2(x+0.001). The infiltration scores of B cells, T cells CD4, T cells CD8, neutrophils, macrophages and DCs in cancers were calculated with TIMER (ref. 43) algorithm from IOBR package (version 0.99.9) (ref. 44) in R software. Moreover, the infiltration scores of B_cell_naive, B_cell_memory, plasma_cells, T_cells_CD8, T_cells_CD4_naive, T_cells_CD4_memory_resting, T_cells_CD4_memory_activated, T_cells_follicular_helper, T_cells_regulatory_(Tregs), T_cells_gamma_delta, NK_cells_activated, Monocytes, Macrophages_M0, Macrophages_M1, Macrophages_M2, Dendritic_cells_resting, Dendritic_cells_activated, Mast_cells_resting, Mast_cells_activated, Eosinophils and Neutrophils in cancers were determined by deconvo_CIBERSORT (ref. 45) algorithm from IOBR package in R software. The association between AKR1B10 expression and immune cell infiltration scores in each cancer was analyzed with corr.test function of psych package (version 2.1.6).

Correlation of AKR1B10 with immune related gene expression

The expression of 41 chemokine genes, 18 chemokine receptors, 21 MHC genes, 24 immunoinhibitors, 46 immunostimulators and 24 immune checkpoint inhibitors (ref. 42) was extracted from normalized pan-cancer RNA-seq data, the correlation of which with AKR1B10 in each cancer was analyzed with corr.test function of psych package (version 2.1.6) after log2(x+0.001) transformed.

Functional analysis

The function of AKR1B10 in pan-cancers was analyzed with gene set enrichment analysis (GSEA) software (version 4.2.3). RNA-seq data of 33 cancers from TCGA project were loaded into GSEA software. C2.cp.kegg.v7.5.1.symbols was chosen as a gene set database, and the number of permutations was set to 1000. The no_collapse was used to gene symbols, and phenotype was set for permutation type. Enriched result lists were first filtered with threshold of |NES|≥1.5, NOM p-value<0.05, and FDR q-value<0.25. Filtered results were further presented with a heatmap of NES.

Multiplex immunofluorescence staining

Gastric cancer tissue microarray (A962) was purchased from WeiAoBio (Shanghai, China). Multiplex immunofluorescence staining procedure on paraffin-embedded tissue section was followed by the instruction of supplier (Absin, Cat No. abs50013). In briefly, the tissue section was deparaffinize, rehydration, antigen retrieval by microwave, serum blocking, primary antibody incubation, secondary antibody incubation, fluorescence staining signal amplification (multiple rounds), staining nuclear with DAPI and mounting. The primary antibodies were AKR1B10 (Abnova, Cat No. H00057016-M01, 1:200), CD68 (Absin, Cat No. 171440, 1:200), CD19 (Abcam, Cat No. 134114, 1:100), PanCK (Absin, Cat No. 123684, 1:400). The staining with the horseradish peroxidase conjugated secondary antibody and the TSA tyramide signal amplification (520-TSA, 620-TSA, 570-TAS, 690-TSA) were supplied within the staining kit. Multispectral images were analyzed. The percentage of positive staining cell was calculated as number of positive cells divide total cell number (P). For the intensity (I), extensive staining was recorded as 3, moderate staining was recorded as 2, weak staining was recorded as 1. $\text{IHC} \text{score}={\displaystyle \sum (P\times I)=(weak percentage\times 1)+(moderate percentage\times 2)+(strong percetage\times 3)}$.

Statistical analysis

The difference of AKR1B10 expression between cancers and corresponding normal tissues was analyzed by unpaired t-test. The significance between two groups was analyzed by unpaired Wilcoxon Rank Sum and Signed Rank tests. The significant difference among multiple groups was analyzed by the Kruskal test. All of the correlations were presented as spearman’s correlation coefficient.

Result

Tissue-specific expression of AKR1B10

To demonstrate the expression of AKR1B10 in human normal tissues, we draw an anatomy map based on AKR1B10 expression from GTEx databases. The result showed that the AKR1B10 highly expressed in stomach and bladder, followed by esophagus, colon, liver, skin and adrenal gland, and AKR1B10 is barely expressed in brain and lymph node (f1). Overall, the gastrointestinal system showed a high-medium AKR1B10 expression, while development gastrointestinal system, genitourinary system, and gynecological system showed a medium-low AKR1B10 expression. The single-cell RNA sequencing (ScRNA Seq) data from Human Cell Landscape (HCL) also showed that the AKR1B10 is highly expressed in the jejunum, stomach, esophagus and transverse-colon (f1).

Aberrant AKR1B10 expression in pan-cancers

AKR1B10 expression is significantly altered in various types of cancers. In general, AKR1B10 mRNA down-regulated in most of gastrointestinal cancer, genitourinary cancer and PCPG, but up-regulated in most of developmental gastrointestinal cancer, some of gynecological cancer, glioma and lung cancer when compared with normal tissues from TCGA databases (f2). Consistently, AKR1B10 protein decreased in STAD and COAD, but increased in LIHC, PAAD, UCEC and CESC when compared with their individual normal tissues (f2). Moreover, we analyzed the association of AKR1B10 mRNA expression with clinical features. The results showed AKR1B10 mRNA expression was significantly correlated with grades of HNSC and KIRC (f2). AKR1B10 was also significantly correlated with stages of KIRC, KIRP, PAAD and ACC (f2).

Both genetic and epigenetic alteration can result in gene expression change. Thus, we analyzed the relationship of DNA mutation, copy number variation (CNV) and DNA methylation with AKR1B10 expression. Mutation analysis indicated that scatter missense mutations (41), nonsense mutations (3), frame shift inserts (1), frame shift deletions (3) and splice sites (1) were observed within the coding sequence of AKR1B10 in 17/33 cancers, of which UCEC had the highest incidence rate of 3.4% (f3). CNV analysis showed that there are no distinct change in AKR1B10 expression among different types of CNV in most cancers (f3). While the methylation analysis showed that AKR1B10 expression was significantly negatively associated with DNA methylation in 12 types of cancer. Specifically, LUSC, CHOL and BLCA showed a moderate correlation (0.4< r <0.6), while LUAD, HNSC, ESCA, LIHC, PAAD, BRCA and CESC showed a weak correlation (0.2< r <0.4) (f4). Taken together, epigenetic might contribute more to the alteration of AKR1B10 expression in pan-cancer.

Potential role of AKR1B10 in tumor microenvironment

Tumor microenvironment (TME) is the non-cancerous cells surround a tumor including fibroblasts and immune cells, which are reported to play a role in tumor growth and drug resistance Purity and ESTIMATE (including stromal score, immune score and ESTIMATE score) algorithm are novel computational methods to infer tumor purity from bulk cell genomic analysis. Purity patterns differ by cancer type (Supplementary Figure 1A). High purity is correlated with shorter OS of PCPG, ACC, ESCA, LUAD, READ,CHOL, CESC and SARC. It is also correlated with longer OS of UVM, MESO, STAD, PAAD, KIRC, BLCA and TGCT (Supplementary Figure 1B). High stromal score is correlated with shorter OS in LGG,GBM, LAML, UVM, MESO, LUSC, STAD, COAD, KIRC, KIRP and BLCA. It is also correlated with longer OS in THCA, ACC, SKCM, LUAD, CHOL, LIHC and SARC (Supplementary Figure 2). High immune score is correlated with shorter OS in LGG, LAML, UVM, ESCA and KIRC. It is also correlated with longer OS in PCPG, THCA, ACC, SKCM, MESO, LUAD, HNSC, COAD, CHOL, LIHC, KICH, BRCA and CESC (Supplementary Figure 3). High ESTIMATE score is correlated with shorter OS in LGG, GBM, LAML, UVM, LUSC, ESCA, STAD, PAAD, KIRC, KIRP and TGCT. It is also correlated with longer OS in PCPG, THCA, ACC, SKCM, LUAD, CHOL, LIHC, CESC and SARC (Supplementary Figure 4).

From the HCL, we knew the AKR1B10 expressed in the epithelial cells, goblet cells, basal cells, enterocytes and enterocyte progenitors, gastric chief cells and hepatocytes, as well as some T cells, B cells, fibroblasts and endothelial cells in normal tissues (f1). The function of AKR1B10 in cancer have attracted the attention of scientists, the role of AKR1B10 in TME remains elusive. Herein, we analyzed the association of AKR1B10 expression with tumor purity (f5). There were 7 cancers stood out, of which the MESO showed positive correlation (r>0.2), but the GBM, THYM, THCA, READ, SKCM, and KIRP showed negative correlation (r<-0.2). Further analysis with ESTIMATE algorithm showed the LGG, GBM, THCA, ACC and UVM were top 5 cancers positively associated with ESTIMATE scores, immune scores and stromal scores, but the CHOL and LUSC were top 2 cancers negatively associated with ESTIMATE scores, immune scores and stromal scores (f5). The correlations in LGG, GBM, CHOL and non-small cell lung cancer (NSCLC) were validated with GEO datasets (Supplementary Figure 5). Taken LGG, GBM, LUSC and CHOL as examples, we find the AKR1B10 expression is elevated in purity low group, stromal score high group, immune score high group and ESTIMATE score high group in LGG and GBM (f5); AKR1B10 expression is also elevated in the purity high group, stromal score low group, immune score low group and ESTIMATE score low group in LUSC and CHOL (f5). Together these data may imply the potential role of AKR1B10 in TME.

Potential role of ARK1B10 in immune responses

We analyzed AKR1B10 expression and immune cells infiltrations with both TIMER and CIBERSORT algorithms. TIMER analysis showed the AKR1B10 expression was associated with multiple immune cells, especially CD4+ T cell, neutrophils and dendritic cells, and these associations were positive in GBM, THCA, THYM and PRAD, but negative in UVM, MESO, LUSC (f6). CIBERSORT analysis showed that the association was mainly focused on T cells CD4+ memory resting, T cells regulatory (Tregs), and dendritic cell resting (f6). AKR1B10 exhibited a broad correlation with immunomodulators involved in 5 immune pathways, among which the correlation with MHC pathway was highest (f7). Our results also supported an association of AKR1B10 expression with dendritic cells and T cells, indicating that AKR1B10 may be involved in antigen presentation.

Potential role of AKR1B10 in immune therapy responses

The immune therapy based on immune checkpoint inhibitors (ICI) is a novel and promising anti-tumor therapy mode. Our analysis showed the correlation between AKR1B10 expression and various ICIs in LGG, GBM, PCPG, THCA, ACC, UVM, KICH, and TGCT, but no correlations in THYM, MESO, LUSC, ESCA, and CHOL (f7).

Tumor mutation burden (TMB) was defined as the number of gene mutations (e.g., base substitutions, insertions and deletions in coding area) per million bases in somatic cells. It is an indicator of mutation frequencies. The higher TMB was usually associated with the favored ICI response regardless of PD-L1 status in various types of tumors. Our analysis showed AKR1B10 expression was significantly associated with TMB in 12 types of cancers, with HNSC correlation coefficient larger than 0.2, THYM, SKCM, ESCA and READ correlation coefficient smaller than -0.2 (f7). Together, these results indicate that AKR1B10 may be a potential predictor and target of immune therapies.

Homologous recombination deficiency (HRD) refers to HRR dysfunction caused by HRR related gene mutations or epigenetic inactivation. HRD often presents in various tumor and produces specific, quantifiable and stable genomic alterations. Clinical research showed HRD status was related to the sensitivity of PARP inhibitors and platinum chemotherapy. We found that AKR1B10 expression significantly correlated with HRD in 9 types of cancer, with ESCA, KIRP and PCPG correlation coefficient lager than 0.2 (f7).

Functional analysis of AKR1B10 in pan-cancer cells

We further explored AKR1B10 function in pan-cancer cells through GSEA. We noticed that AKR1B10 was enriched with cancer cell features and pathways, as well as apoptosis, cell cycle, autophagy, and ECM interaction processes in some cancers. Further studies revealed that AKR1B10 was most remarkably and positively enriched with multiple metabolism processes, including drugs, hormones, saccharides, and amino acids in pan-cancer (f8). These results indicate a potential role of AKR1B10 in cancer development and progression.

The prognostic role of AKR1B10 in pan-cancers

Subsequently, we evaluated the prognostic value of AKR1B10 in multiple types of cancers through Kaplan-Meier analysis. As shown in f9, the increased AKR1B10 expression was significantly related to poor overall survival (OS) in GBM, THCA, UVM, SKCM, MESO, LIHC, PAAD, KICH, KIRC, KIRP, and PRAD, but better OS in ACC, COAD, and READ. In addition, elevated AKR1B10 expression was associated with shorter PFI in GBM, SKCM, STAD, PAAD, KIRC, KIRP, TGCT, and BRCA, as well as decreased DFI in ACC, STAD, LIHC, BLCA, and TGCT. On the contrary, increased AKR1B10 expression was associated with longer PFI in ACC, COAD, and CESC, and also the prolonged DFI in THCA, BRCA, and CESC.

Validation of AKR1B10 expression and its association with immune in gastric cancer

Our previous analysis showed stomach was the organ with highest AKR1B10 expression. Some studies indicated that decreased AKR1B10 expression promote M2 macrophage marker expression in THP1 macrophage cells (ref. 46), promote NK cell kill activity on liver cancer cell Huh7 (ref. 47), and knock down of AKR1B8 transgenic mice (homologous of AKR1B10 in mice) show a immunity disorder in colon (ref. 48), however, the correlation between AKR1B10 expression and immune infiltration in STAD remains elusive. We evaluated AKR1B10 expression in paired normal and gastric cancer tissues, and the correlation of AKR1B10 expression with immune cells in gastric cancer (f10). The results showed that AKR1B10 expression was dramatically decreased in gastric cancer compared to normal tissues (f10). Meantime, the CD19 expression (B cell marker) was also decreased in gastric cancer (f10), and CD68 (Macrophage marker) expression elevated in gastric cancer (f10). The PanCK (Epithelial marker) expression difference is not significant between gastric cancer and normal tissue (f10).

As AKR1B10 expression almost diminished in gastric cancer, we evaluated the relationship between AKR1B10 and immunity by dividing samples with AKR1B10 expression into AKR1B10 high group with H score > 10, AKR1B10 medium group with 10 > H score > 0.5 and AKR1B10 low group with H score < 0.5. We found CD68 show an elevated trend as AKR1B10 expression decreased (f10), which is consistent with the TIMER analysis (f6). However, either expression difference of CD19 or PanCK is clear (f10).

Discussion

In physiological conditions, AKR1B10 eliminates various endogenous and exogenous cytotoxic carbonyl compounds and mediates long-chain fatty acid/lipid synthesis, playing an important role in intestinal hemostasis. However, the role of AKR1B10 in different cancers is complex and remains elusive. In this study, we systematically analyzed the expression pattern, genetic alterations, and epigenetic alterations in 33 cancer types. We evaluated the association of AKR1B10 expression with components of tumor microenvironment and clinical outcomes. Our study provides a landscape of AKR1B10 in pan-cancers and suggests a potential role of AKR1B10 in tumor microenvironment remodeling. This study also takes gastric cancer as a role model to examine the AKR1B10 expression and its association with TME.

In literature, AKR1B10 mRNA was detected at a high level in adrenal gland and intestinal tract (ref. 11), and the AKR1B10 protein was expressed in small intestine, colon, and liver (ref. 1). AKR1B10 was reported to be upregulated in liver cancer (ref. 23, ref. 49–ref. 51), lung cancer (ref. 27, ref. 28), breast cancer (ref. 25), oral cancer (ref. 30–ref. 32), pancreatic cancer (ref. 52), cervical cancer (ref. 53) and papillary renal cell carcinoma (ref. 54). AKR1B10 was reported to suppress tumor development and progression in colorectal cancer (ref. 33–ref. 37) and gastric cancer (ref. 55), but remains controversial in nasopharyngeal cancer (ref. 56, ref. 57) and esophageal cancer (ref. 58, ref. 59). In this pan-cancer analysis we found that AKR1B10 is increased in GBM, LUAD, LUSC, CHOL, LIHC, PAAD, KIRP, UCEC and CESC, but decreased in STAD, COAD, READ, KICH, BLCA and PRAD. However, the BRCA and HNSC (including oral cancer and nasopharyngeal cancer) didn’t show significant alterations in AKR1B10 expression, which suggests the analysis with a refined cancer subtype may be helpful to clarify the AKR1B10 expression in special pathological conditions. As for prognostic values, AKR1B10 expression was reported to be correlated with poor prognosis in liver cancer (ref. 50), oral cancer (ref. 30–ref. 32), breast cancer (ref. 25), lung cancer (ref. 29), but with favored prognosis in CRC (ref. 34, ref. 35). This pan-cancer analysis suggested that AKR1B10 also associated with poor prognosis in GBM, SKCM, STAD, PAAD, KIRC, KIRP, and TGCT, but with favored prognosis in ACC and CESC.

Tumor microenvironment underlies the basis of tumor growth, which contains kinds of cell and molecular components, forming a complex interaction net. Although the single cell RNA-Sequencing revealed AKR1B10 expression was mainly restrained in epithelial cells, we found that AKR1B10 expression was significantly associated with tumor purity, stromal cells, various immune cells infiltration (e.g., neutrophils, dendritic cells and CD4+ T cells) and immune gene expression. Cancer-associated-fibroblasts (CAF) were the main component of stroma, and participated in bidirectional regulation of tumorigenesis (ref. 60). Immune cells like CD4+ T cells and neutrophils were also reported to be dual functional in tumor immunity. CD4+ T cells contains a few subpopulations including Th1, Th2, Th9, Th17, follicular helpher T (TFH) cells and T cells regulatory (Treg). Among them, stable Tregs, which express foxp3, suppress immune responses by preventing Tconv cell activation, inducing B cell apoptosis, and restraining the antigen-presentation of DCs, while fragile Tregs exert an antitumor effect through secreting IFN-γ, and unstable Tregs was not suppressive (ref. 61). Tumor-associated neutrophils (TANs) are an important component of TME as well. Two types of TANs are identified in TME. The antitumor N1 kills tumor cells via oxidative stress and inflammatory factors. Pro-tumor N2 promotes tumor growth and dissemination via remodeling TME or forming neutrophils extracellular traps (NETs) (ref. 62). Neutrophils are novel potential targets of anti-cancer therapy (ref. 63). Dendritic cells (ref. 64) and immune checkpoints (ref. 65) both play crucial roles in immune responses, and are dazzling stars in tumor immunotherapy. One of the possible mechanisms linking AKR1B10 expression in epithelial cells with TME components may be the metabolism re-programming. For example, AKR1B10 may create an immunosuppressive TME by promoting long-chain fatty acid/lipid synthesis. These results suggest that AKR1B10 may be a promising regulator of TME immune status. Interestingly, our data showed for the first time that AKR1B10 expression was significantly associated with DNA methylation in pan-cancers. Further study is warranted.

This study verified AKR1B10 expression and immune infiltration in gastric cancer, and our results consistently demonstrated that AKR1B10 expression was diminished in gastric cancer accompanied with CD19+ B cell proportion decline and CD68+ macrophage increase. However, the function and prognostic value of AKR1B10 is still obscure, and the relationship with immune infiltration is largely unknown. For example, Yao et al. report that high AKR1B10 expression is associated with smaller tumor size, less metastasis and high survival rate of gastric cancer (ref. 55), but Ahmed et al. exhibit that high AKR1B10 is correlated with high metastasis, tumor progression, poor responses to chemotherapy and low survival rate in gastric cancer (ref. 66). Two independent studies show AKR1B10 may induce drug resistance of gastric cancer (ref. 67, ref. 68), but the other report that AKR1B10 inhibit proliferation and migration of gastric cancer cells (ref. 69). The B cells are major components of tertiary lymphoid structures in gastric cancer, and associated with favored prognosis (ref. 70). Tumor associated macrophage contains M1 (anti-tumor) and M2 (pro-tumor) macrophage, deletion of AK1B10 results in a decrease to M1 polarization markers and an increase to M2 polarization markers (ref. 46). Considering the essential biological function of AKR1B10 and its downregulation in gastric cancer, we would like to consider AKR1B10 as a protective factor in normal gastric tissues. It is possible that AKR1B10 is important for the gastric cancer development, and relate to macrophage infiltration.

Conclusion

This study provided a landscape of AKR1B10 in pan-cancers, including its expression and potential roles in TME, tumor immunity, and survival. The data suggested that AKR1B10 may be a potential biomarker and therapeutic target for multiple cancers. This study particularly revealed the diminished expression of AKR1B10 in gastric cancer, and discussed the potential role of AKR1B10 in the gastric cancer TME.

References

1. D Cao, ST Fan, SS Chung. Identification and characterization of a novel human aldose reductase-like gene.. J Biol Chem. (, 1998. [DOI]
2. S Endo, M Toshiyuki, M Hiroaki, O Chisato, S Midori, K Yukio. Kinetic studies of AKR1B10, human aldose reductase-like protein: endogenous substrates and inhibition by steroids.. Arch Biochem Biophys. (, 2009. [DOI | PubMed]
3. J Giménez-Dejoz, MH Kolář, FX Ruiz, I Crespo, A Cousido-Siah, A Podjarny. Substrate specificity, inhibitor selectivity and structure-function relationships of aldo-keto reductase 1B15: A novel human retinaldehyde reductase.. PloS One. (, 2015. [DOI | PubMed]
4. J Giménez-Dejoz, S Weber, Á Fernández-Pardo, G Möller, J Adamski, S Porté. Engineering aldo-keto reductase 1B10 to mimic the distinct 1B15 topology and specificity towards inhibitors and substrates, including retinoids and steroids.. Chem Biol Interact. (, 2019. [DOI]
5. S Weber, JK Salabei, G Möller, E Kremmer, A Bhatnagar, J Adamski. Aldo-keto Reductase 1B15 (AKR1B15): a mitochondrial human aldo-keto reductase with activity toward steroids and 3-keto-acyl-CoA conjugates.. J Biol Chem. (, 2015. [DOI]
6. Z Kabututu, M Manin, JC Pointud, T Maruyama, N Nagata, S Lambert. Prostaglandin F2alpha synthase activities of aldo-keto reductase 1B1, 1B3 and 1B7.. J Biochem. (, 2009. [DOI]
7. T Matsunaga, Y Wada, S Endo, M Soda, O El-Kabbani, A Hara. Aldo-keto reductase 1B10 and its role in proliferation capacity of drug-resistant cancers.. Front Pharmacol. (, 2012. [DOI | PubMed]
8. OS Bains, TA Grigliatti, RE Reid, KW Riggs. Naturally occurring variants of human aldo-keto reductases with reduced in vitro metabolism of daunorubicin and doxorubicin.. J Pharmacol Exp Ther. (, 2010. [DOI]
9. AM Quinn, RG Harvey, TM Penning. Oxidation of PAH trans-dihydrodiols by human aldo-keto reductase AKR1B10.. Chem Res Toxicol. (, 2008. [DOI]
10. TM Penning. The aldo-keto reductases (AKRs): Overview.. Chem Biol Interact. (, 2015. [DOI]
11. DJ Hyndman, TG Flynn. Sequence and expression levels in human tissues of a new member of the aldo-keto reductase family.. Biochim Biophys Acta. (, 1998. [DOI | PubMed]
12. VN Sumantran, P Mishra, R Bera, N Sudhakar. Microarray analysis of differentially-expressed genes encoding CYP450 and phase II drug metabolizing enzymes in psoriasis and melanoma.. Pharmaceutics. (, 2016. [DOI]
13. N Jumper, T Hodgkinson, G Arscott, Y Har-Shai, R Paus, A Bayat. The aldo-keto reductase AKR1B10 is up-regulated in keloid epidermis, implicating retinoic acid pathway dysregulation in the pathogenesis of keloid disease.. J Invest Dermatol. (, 2016. [DOI]
14. CT Soares, LRV Fachin, APF Trombone, PS Rosa, CC Ghidella, AFF Belone. Potential of AKR1B10 as a biomarker and therapeutic target in type 2 leprosy reaction.. Front Med (Lausanne). (, 2018. [DOI | PubMed]
15. Y Shen, J Ma, R Yan, H Ling, X Li, W Yang. Impaired self-renewal and increased colitis and dysplastic lesions in colonic mucosa of AKR1B8-deficient mice.. Clin Cancer Res. (, 2015. [DOI]
16. F Zeng, Y Zhang, X Han, M Zeng, Y Gao, J Weng. Predicting non-alcoholic fatty liver disease progression and immune deregulations by specific gene expression patterns.. Front Immunol. (, 2020. [DOI | PubMed]
17. J Starmann, M Fälth, W Spindelböck, KL Lanz, C Lackner, K Zatloukal. Gene expression profiling unravels cancer-related hepatic molecular signatures in steatohepatitis but not in steatosis.. PloS One. (, 2012. [DOI | PubMed]
18. H Tsuzura, T Genda, S Sato, A Murata, Y Kanemitsu, Y Narita. Expression of aldo-keto reductase family 1 member b10 in the early stages of human hepatocarcinogenesis.. Int J Mol Sci. (, 2014. [DOI]
19. A Bitter, P Rümmele, K Klein, BA Kandel, JK Rieger, AK Nüssler. Pregnane X receptor activation and silencing promote steatosis of human hepatic cells by distinct lipogenic mechanisms.. Arch Toxicol. (, 2015. [DOI]
20. P Pettinelli, BM Arendt, A Teterina, I McGilvray, EM Comelli, SK Fung. Altered hepatic genes related to retinol metabolism and plasma retinol in patients with non-alcoholic fatty liver disease.. PloS One. (, 2018. [DOI | PubMed]
21. O Govaere, S Cockell. Transcriptomic profiling across the nonalcoholic fatty liver disease spectrum reveals gene signatures for steatohepatitis and fibrosis.. Sci Transl Med. (, 2020. [DOI]
22. M Kanno, K Kawaguchi, M Honda, R Horii, H Takatori, T Shimakami. Serum aldo-keto reductase family 1 member B10 predicts advanced liver fibrosis and fatal complications of nonalcoholic steatohepatitis.. J Gastroenterol. (, 2019. [DOI]
23. Z Scuric, SC Stain, WF Anderson, JJ Hwang. New member of aldose reductase family proteins overexpressed in human hepatocellular carcinoma.. Hepatology. (, 1998. [DOI]
24. J Ma, R Yan, X Zu, JM Cheng, K Rao, DF Liao. Aldo-keto reductase family 1 B10 affects fatty acid synthesis by regulating the stability of acetyl-CoA carboxylase-alpha in breast cancer cells.. J Biol Chem. (, 2008. [DOI]
25. J Ma, DX Luo, C Huang, Y Shen, Y Bu, S Markwell. AKR1B10 overexpression in breast cancer: association with tumor size, lymph node metastasis and patient survival and its potential as a novel serum marker.. Int J Cancer. (, 2012. [DOI]
26. KA Reddy, PU Kumar, M Srinivasulu, B Triveni, K Sharada, A Ismail. Overexpression and enhanced specific activity of aldoketo reductases (AKR1B1 & AKR1B10) in human breast cancers.. Breast. (, 2017. [DOI]
27. AK MacLeod, L Acosta-Jimenez, PJ Coates, M McMahon, FA Carey, T Honda. Aldo-keto reductases are biomarkers of NRF2 activity and are co-ordinately overexpressed in non-small cell lung cancer.. Br J Cancer. (, 2017. [DOI]
28. S Fukumoto, N Yamauchi, H Moriguchi, Y Hippo, A Watanabe, J Shibahara. Overexpression of the aldo-keto reductase family protein AKR1B10 is highly correlated with smokers’ non-small cell lung carcinomas.. Clin Cancer Res. (, 2005. [DOI]
29. JJ Hung, YC Yeh, WH Hsu. Prognostic significance of AKR1B10 in patients with resected lung adenocarcinoma.. Thorac Cancer. (, 2018. [DOI]
30. CY Fang, YH Lin, CL Chen. Overexpression of AKR1B10 predicts tumor recurrence and short survival in oral squamous cell carcinoma patients.. J Oral Pathol Med. (, 2019. [DOI]
31. HH Ko, SL Cheng. Expression of AKR1B10 as an independent marker for poor prognosis in human oral squamous cell carcinoma.. Head Neck. (, 2017. [DOI]
32. HH Ko, HH Peng, SJ Cheng. Increased salivary AKR1B10 level: Association with progression and poor prognosis of oral squamous cell carcinoma.. Head Neck. (, 2018. [DOI]
33. Y Yao, X Wang, D Zhou, H Li, H Qian, J Zhang. Loss of AKR1B10 promotes colorectal cancer cells proliferation and migration via regulating FGF1-dependent pathway.. Aging (Albany NY). (, 2020. [DOI]
34. B Taskoparan, EG Seza, S Demirkol, S Tuncer, M Stefek, AO Gure. Opposing roles of the aldo-keto reductases AKR1B1 and AKR1B10 in colorectal cancer.. Cell Oncol. (, 2017. [DOI]
35. T Ohashi, M Idogawa, Y Sasaki, H Suzuki, T Tokino. AKR1B10, a transcriptional target of p53, is downregulated in colorectal cancers associated with poor prognosis.. Mol Cancer Res. (, 2013. [DOI]
36. ES Kropotova, RA Tychko, OL Zinov’eva, AF Zyrianova, SL Khankin, VL Cherkes. Downregulation of AKR1B10 gene expression in colorectal cancer.. Mol Biol (Mosk). (, 2010. [DOI]
37. OL Zinovieva, EN Grineva, GS Krasnov, DS Karpov, AO Zheltukhin, AV Snezhkina. Treatment of cancer cells with chemotherapeutic drugs results in profound changes in expression of genes encoding aldehyde-metabolizing enzymes.. J Cancer. (, 2019. [DOI]
38. R Beroukhim, CH Mermel, D Porter, G Wei, S Raychaudhuri, J Donovan. The landscape of somatic copy-number alteration across human cancers.. Nature. (, 2010. [DOI | PubMed]
39. CH Mermel, SE Schumacher, B Hill, ML Meyerson, R Beroukhim, G Getz. GISTIC2.0 facilitates sensitive and confident localization of the targets of focal somatic copy-number alteration in human cancers.. Genome Biol. (, 2011. [DOI | PubMed]
40. W Ding, J Chen, G Feng, G Chen, J Wu, Y Guo. DNMIVD: DNA methylation interactive visualization database.. Nucleic Acids Res. (, 2020. [DOI | PubMed]
41. X Han, Z Zhou, L Fei, H Sun, R Wang, Y Chen. Construction of a human cell landscape at single-cell level.. Nature. (, 2020. [DOI]
42. V Thorsson, LG David, DB Scott, W Denise, SB Dante, OY Tai-Hsien. The immune landscape of cancer.. Immunity. (, 2018. [DOI | PubMed]
43. T Li, J Fan, B Wang, N Traugh, Q Chen, JS Liu. TIMER: A web server for comprehensive analysis of tumor-infiltrating immune cells.. Cancer Res. (, 2017. [DOI]
44. D Zeng, Y Zilan, S Rongfang, Y Guangchuang, W Jiani, X Yi. IOBR: multi-omics immuno-oncology biological research to decode tumor microenvironment and signatures.. Front Immunol. (, 2021. [DOI | PubMed]
45. AM Newman, CL Liu, MR Green. Robust enumeration of cell subsets from tissue expression profiles.. Nature Methods. (, 2015. [DOI]
46. Y Wu, Y Hao, Q Zhuang, X Ma, C Shi. AKR1B10 regulates M2 macrophage polarization to promote the Malignant phenotype of gastric cancer.. Biosci Rep. (, 2023. [DOI]
47. T Wu, Y Ke, H Tang, C Liao, J Li, L Wang. Fidarestat induces glycolysis of NK cells through decreasing AKR1B10 expression to inhibit hepatocellular carcinoma.. Mol Ther Oncolytics. (, 2021. [DOI]
48. X Wang, K Ramina, S Yi, C Yu, L Minglin, Z Yun. Impaired barrier function and immunity in the colon of aldo-keto reductase 1B8 deficient mice.. Front Cell Dev Biol. (, 2021. [DOI | PubMed]
49. X Ye, C Li, X Zu, M Lin, Q Liu, J Liu. A large-scale multicenter study validates aldo-keto reductase family 1 member B10 as a prevalent serum marker for detection of hepatocellular carcinoma.. Hepatology. (, 2019. [DOI]
50. J Shi, L Chen, Y Chen, Y Lu, X Chen, Z Yang. Aldo-Keto Reductase Family 1 Member B10 (AKR1B10) overexpression in tumors predicts worse overall survival in hepatocellular carcinoma.. J Cancer. (, 2019. [DOI]
51. JK DiStefano, B Davis. Diagnostic and prognostic potential of AKR1B10 in human hepatocellular carcinoma.. Cancers (Basel). (, 2019. [DOI]
52. YT Chung, KA Matkowskyj, H Li, H Bai, W Zhang, MS Tsao. Overexpression and oncogenic function of aldo-keto reductase family 1B10 (AKR1B10) in pancreatic carcinoma.. Mod Pathol. (, 2012. [DOI]
53. H Yoshitake, M Takahashi, H Ishikawa, M Nojima, H Iwanari, A Watanabe. Aldo-keto reductase family 1, member B10 in uterine carcinomas: a potential risk factor of recurrence after surgical therapy in cervical cancer.. Int J Gynecol Cancer. (, 2007. [DOI]
54. A Ooi, JC Wong, D Petillo, D Roossien, V Perrier-Trudova, D Whitten. An antioxidant response phenotype shared between hereditary and sporadic type 2 papillary renal cell carcinoma.. Cancer Cell. (, 2011. [DOI]
55. HB Yao, Y Xu, LG Chen, TP Guan, YY Ma, XJ He. AKR1B10, a good prognostic indicator in gastric cancer.. Eur J Surg Oncol. (, 2014. [DOI]
56. YC He, Y Shen, Y Cao, FQ Tang, DF Tian, CF Huang. Overexpression of AKR1B10 in nasopharyngeal carcinoma as a potential biomarker.. Cancer biomark. (, 2016. [DOI]
57. Y Guo, W Luo, Z Hu, J Li, X Li, H Cao. Low expression of Aldo-keto reductase 1B10 is a novel independent prognostic indicator for nasopharyngeal carcinoma.. Cell Biosci. (, 2016. [DOI | PubMed]
58. J Breton, MC Gage, AW Hay, JN Keen, CP Wild, C Donnellan. Proteomic screening of a cell line model of esophageal carcinogenesis identifies cathepsin D and aldo-keto reductase 1C2 and 1B10 dysregulation in Barrett’s esophagus and esophageal adenocarcinoma.. J Proteome Res. (, 2008. [DOI]
59. DJ Nancarrow, AD Clouston, BM Smithers, DC Gotley, PA Drew, DI Watson. Whole genome expression array profiling highlights differences in mucosal defense genes in Barrett’s esophagus and esophageal adenocarcinoma.. PloS One. (, 2011. [DOI | PubMed]
60. T Liu, C Han, S Wang, P Fang, Z Ma, L Xu. Cancer-associated fibroblasts: an emerging target of anti-cancer immunotherapy.. J Hematol Oncol. (, 2019. [DOI | PubMed]
61. AE Overacre-Delgoffe, DAA Vignali. Treg fragility: A prerequisite for effective antitumor immunity?. Cancer Immunol Res. (, 2018. [DOI]
62. S Jaillon, A Ponzetta. Neutrophil diversity and plasticity in tumour progression and therapy.. Nat Rev Cancer. (, 2020. [DOI | PubMed]
63. IL Linde, TR Prestwood, J Qiu, G Pilarowski, MH Linde, X Zhang. Neutrophil-activating therapy for the treatment of cancer.. Cancer Cell. (, 2023. [DOI | PubMed]
64. TA Patente, MP Pinho, AA Oliveira, GCM Evangelista, PC Bergami-Santos, JAM Barbuto. Human dendritic cells: their heterogeneity and clinical application potential in cancer immunotherapy.. Front Immunol. (, 2018. [DOI | PubMed]
65. G Morad, BA Helmink, P Sharma, JA Wargo. Hallmarks of response, resistance, and toxicity to immune checkpoint blockade.. Cell. (, 2022. [DOI | PubMed]
66. SMU Ahmed, ZN Jiang, ZH Zheng, Y Li, XJ Wang, X Tang. AKR1B10 expression predicts response of gastric cancer to neoadjuvant chemotherapy.. Oncol Lett. (, 2019. [DOI]
67. Y Morikawa, C Kezuka, S Endo, A Ikari, M Soda, K Yamamura. Acquisition of doxorubicin resistance facilitates migrating and invasive potentials of gastric cancer MKN45 cells through up-regulating aldo-keto reductase 1B10.. Chem Biol Interact. (, 2015. [DOI]
68. T Matsunaga, A Suzuki, C Kezuka, N Okumura, K Iguchi, I Inoue. Aldo-keto reductase 1B10 promotes development of cisplatin resistance in gastrointestinal cancer cells through down-regulating peroxisome proliferator-activated receptor-γ-dependent mechanism.. Chem Biol Interact. (, 2016. [DOI]
69. X Shao, J Wu, S Yu, Y Zhou, C Zhou. AKR1B10 inhibits the proliferation and migration of gastric cancer via regulating epithelial-mesenchymal transition.. Aging (Albany NY). (, 2021. [DOI]
70. D Bruni, HK Angell. The immune contexture and Immunoscore in cancer prognosis and therapeutic efficacy.. Nat Rev Cancer. (, 2020. [DOI]