Next Article in Journal
Single-Cell Heterogeneity of Cutaneous T-Cell Lymphomas Revealed Using RNA-Seq Technologies
Previous Article in Journal
Genomic Database Analysis of Uterine Leiomyosarcoma Mutational Profile
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Clinicopathological and Prognostic Significance of Nrf2 and Keap1 Expression in Hepatocellular Carcinoma

1
Department of Pathology, Seoul National University Hospital, Seoul National University College of Medicine, Seoul 03080, Korea
2
Department of Pathology, Seoul National University Bundang Hospital, Seoul National University College of Medicine, Seongnam 13620, Korea
3
Department of Pathology, Gangnam Severance Hospital, Yonsei University College of Medicine, Seoul 06273, Korea
4
Department of Pathology and Translational Genomics, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul 06351, Korea
5
Department of Internal Medicine and Liver Research Institute, Seoul National University College of Medicine; Biomedical Research Institute, Center for Medical Innovation, Seoul National University Hospital, Seoul 03080, Korea
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(8), 2128; https://doi.org/10.3390/cancers12082128
Submission received: 16 June 2020 / Revised: 16 July 2020 / Accepted: 28 July 2020 / Published: 31 July 2020

Abstract

:
Nuclear factor E2-related factor2 (Nrf2) activation is associated with both cytoprotective effects and malignant behavior of cancer cells. This study aimed to evaluate the clinicopathological implications of the expression of Nrf2, pNrf2, and its regulator Keap1 in human hepatocellular carcinomas (HCCs). Tissue microarrays consisting of 285 surgically resected HCCs were immunohistochemically stained with pNrf2, Nrf2, Keap1, stemness-related markers (keratin 19 (K19), epithelial cell adhesion molecule (EpCAM)), carbonic anhydrase IX (CAIX), epithelial–mesenchymal transition (EMT)-related markers (ezrin, uPAR, E-cadherin), and p53, and the results were correlated with the clinicopathological features. pNrf2 expression was significantly associated with increased proliferative activity, as well as EpCAM, ezrin, p53, and CAIX expression and E-cadherin loss (p < 0.05, all). Strong cytoplasmic Nrf2 expression was associated with CAIX and ezrin expression (p < 0.05, both). Keap1 was associated with increased proliferative activity, portal vein invasion, EMT-related markers, and p53 expression in CAIX-negative HCCs (p < 0.05, all). Both pNrf2 and cytoplasmic Nrf2 expression were associated with decreased overall survival (p < 0.05, both), and cytoplasmic Nrf2 expression was an independent predictor of decreased overall survival on multivariate analysis (hazard ratio 4.15, p < 0.001). Both pNrf2 and cytoplasmic Nrf2 expression were associated with poor survival and aggressive behavior of HCC. In addition, Keap1 expression was also associated with aggressive HCC behavior in CAIX-negative HCCs, suggesting that Keap1 expression should be interpreted in the context of hypoxia status.

1. Introduction

Hepatocellular carcinoma (HCC) is the fourth leading cause of cancer-related mortality worldwide [1]. The genomic landscape of HCC was unraveled during the past decade, with the most common genetic alterations being mutations in TP53, CTNNB1, and TERT promoter genes [2,3,4]; following the identification of distinct morphological subtypes of HCCs that partly correlate with the molecular features, there was increasing interest in the heterogeneity of HCC in recent years [5]. This heterogeneity has important clinical implications, as understanding the morphomolecular features of the different subtypes of HCC will lead to increased treatment options and provide background for drug development.
Of the less frequent molecular alterations of HCC, somatic mutations in NFE2L2 and KEAP1 are observed in 3–6% and 2–8% of HCCs, respectively [4,6]. Activating NFE2L2 mutations and inactivating KEAP1 mutations were demonstrated to play important roles in tumor progression by preventing the degradation of nuclear factor E2-related factor2 (Nrf2), and the Nrf2/Keap1 pathway recently attracted attention as a potential therapeutic target for HCC [3,7,8]. Under physiological conditions, Nrf2 has cytoprotective functions; oxidative stress induces the phosphorylation of Nrf2, and phospho-Nrf2 (pNrf2) is subsequently released from Keap1 to enter the nucleus, bind to antioxidant response elements, and ultimately trigger the expression of cytoprotective genes [6,9,10]. However, although this cytoprotective mechanism initially prevents carcinogenesis, the overexpression of Nrf2 in cancer cells may result in protection of cancer cells from oxidative stress, enhancement of cancer cell proliferation and survival, and resistance to therapy [8,11,12,13,14].
Although the prognostic and clinical implication of Nrf2 and Keap1 protein expression was recently demonstrated in various malignant neoplasms, previous studies reported conflicting results and there are limited data on HCC [15,16,17,18,19,20,21,22]. In this study, we evaluated the prognostic significance and the clinicopathological features of HCCs expressing Nrf2, pNrf2, and Keap1, and we correlated the findings with the expression status of markers associated with hypoxia (carbonic anhydrase IX; CAIX), stemness (keratin 19; K19, epithelial cell adhesion molecule; EpCAM) and epithelial–mesenchymal transition (EMT) (ezrin, urokinase-type plasminogen activator receptor (uPAR) and E-cadherin).

2. Results

The clinicopathological features of the 285 HCC cases are summarized in Table 1. The most common etiology was hepatitis B virus (HBV) infection (203/285, 71.2%). Multiplicity was present in 49/285 (17.2%) cases, and 205/285 (72.0%) cases were poorly differentiated (Edmondson–Steiner grade 3 or 4). Microvascular and portal vein invasion was seen in 113/285 (39.6%) and 19/285 (6.7%) cases, respectively. After a follow up period of 66.5 months (range, 0–163 months), 161/285 (56.5%) cases recurred and 54/285 (18.9%) patients died.

2.1. pNrf2 Expression in HCCs and Clinicopathological Correlation

The immunohistochemical stain results and clinicopathological data are demonstrated in Figure 1 and Table 2. pNrf2 expression was observed in 73/285 (25.6%) cases, and all positive cases demonstrated strong nuclear labeling. High pNrf2 expression was associated with significantly higher proliferative activity (mitotic index: 17.4 ± 20.0% vs. 10.1 ± 13.4%, respectively, p < 0.001; Ki-67 labeling index (12.9 ± 12.2 vs. 5.1 ± 7.5/10 high-power field (HPF), respectively, p < 0.001) compared with pNrf2-low HCCs. pNrf2-high HCCs demonstrated significantly more frequent expression of EpCAM (41/73 (56.2%) vs. 84/212 (39.6%), respectively, p = 0.020), CAIX (27/72 (37.5%) vs. 48/205 (23.4%), respectively, p = 0.030), ezrin (43/72 (59.7%) vs. 78/212 (36.8%), respectively, p = 0.001), E-cadherin loss (56/73 (76.7%) vs. 128/212 (60.4%), respectively, p = 0.015), and p53 overexpression (26/70 (37.1%) vs. 40/195 (20.5%), respectively, p = 0.009), compared with pNrf2-low HCCs (Figure 2). K19 expression was also slightly more frequent in pNrf2-high HCCs; however, the difference was not statistically significant. pNrf2-high status was significantly associated with decreased overall survival (OS) on univariate analysis (hazard ratio (HR), 1.89: 95% confidence interval (CI), 1.04–3.46; p = 0.038) (Table 3, Figure 3). There was no significant difference in disease-free survival (DFS) according to pNrf2 status. However, analysis of an independent cohort of 293 surgically resected HCCs also revealed a significant association between pNrf2 expression and decreased OS (p = 0.001) and DFS (p = 0.029) (Figure S1).

2.2. Cytoplasmic Nrf2 Expression in HCCs and Clinicopathological Correlation

Cytoplasmic Nrf2 expression was observed in 48/285 (16.8%) of HCCs. There was no correlation between pNrf2 and Nrf2 expression status. Nrf2-high HCCs were smaller in size compared with Nrf2-low HCCs (3.1 ± 1.8 cm vs. 4.4 ± 2.9 cm, respectively, p = 0.005). Cytoplasmic Nrf2 expression was associated with increased expression of CAIX (20/48 (41.7%) vs. 55/229 (24.0%), respectively, p = 0.019) and Ezrin (31/48 (64.6%) vs. 90/236 (38.1%), respectively, p = 0.001) compared to Nrf2-low HCCs. Notably, Nrf2-high status was significantly associated with decreased OS in HCC patients on univariate analysis (HR, 3.33; 95% CI, 1.74–6.37; p < 0.001) (Table 3, Figure 3).

2.3. Keap1 Expression in HCCs and Clinicopathological Correlation

Keap1 was expressed in the cytoplasm of 161/285 (56.5%) HCCs. Interestingly, a significant positive correlation was observed between pNrf2 and Keap1 expression (p = 0.020). Keap1 expression was associated with increased expression of uPAR (29/124 (23.4%) vs. 55/159 (34.6%), respectively, p = 0.049). In addition, more frequent p53 expression (47/147 (32.0%) vs. 19/118 (16.1%), respectively, p = 0.004) and higher proliferative activity (Ki-67 labeling index: 8.1 ± 10.2% vs. 5.7 ± 8.5%, respectively, p = 0.035; mitotic index: 13.1 ± 16.9/10 HPF vs. 10.5 ± 13.6/10 HPF, respectively, p = 0.056) were noted in Keap1-high HCCs, compared to Keap1-low HCCs.
The positive correlation between pNrf2 and Keap1 and the association between Keap1 and increased proliferative activity in this study were unexpected, as Keap1 is known to suppress Nrf2 activation and, thus, expected to have tumor-suppressive effects. We, thus, analyzed the clinicopathological features of Keap1-positive HCCs separately in CAIX-positive (n = 75) and CAIX-negative HCCs (n = 202). Of note, we found that Keap1 expression was associated with increased Ki-67 labeling index (7.4 ± 10.3% vs. 4.6 ± 7.1%, respectively, p = 0.008), more frequent portal vein invasion (1/202 (0.9%) vs. 6/75 (13.6%), respectively, p = 0.024), ezrin expression (3/91 (3.3%) vs. 13/110 (11.8%), respectively, p = 0.035), uPAR expression (15/91 (16.5%) vs. 33/110 (30.0%), respectively, p = 0.031), and p53 expression (31/202 (29.8%) vs. 14/75 (34.1%), respectively, p = 0.009) only in CAIX-negative HCCs, while there was no association between Keap1 status and the clinicopathological factors in CAIX-positive HCCs (Table S1). In addition, although there were no differences in patient survival according to Keap1 expression status, there was interestingly a tendency for decreased DFS for Keap1-high HCCs in the CAIX-negative HCC group (p = 0.062).

2.4. Multivariate Survival Analysis

Multivariate survival analysis was performed with clinicopathological and immunohistochemical parameters that demonstrated statistical significance on univariate analysis, including pNrf2 and cytoplasmic Nrf2. Other parameters associated with decreased OS on univariate analysis included tumor multiplicity (HR, 2.81; 95% CI, 1.49–5.31; p = 0.001), tumor size (≥5 cm vs. <5 cm) (HR, 2.40; 95% CI, 1.32–4.38; p = 0.004), portal vein invasion (HR, 3.29; 95% CI, 1.39–7.80; p = 0.007), the presence of extrahepatic metastasis (HR, 4.46; 95% CI, 2.45–8.13; p < 0.001), and expression of CAIX (HR, 1.85; 95% CI, 1.02–3.39; p = 0.044) and ezrin (HR, 3.16; 95% CI, 1.67–5.97; p < 0.001). Multivariate analysis revealed cytoplasmic Nrf2 expression to be a significant independent prognostic factor for OS (HR, 4.35; 95% CI, 2.11–8.97; p < 0.001), along with tumor multiplicity (HR, 1.98; 95% CI, 1.02–3.86; p = 0.045), presence of extrahepatic metastasis (HR, 3.75; 95% CI, 1.91–7.36; p < 0.001), and ezrin expression (HR, 2.41; 95% CI, 1.25–4.64; p = 0.009) (Table 3). As for DFS, tumor multiplicity (HR, 1.49; 95% CI, 1.02–2.18; p = 0.042) and the presence of extrahepatic metastasis (HR, 4.48; 95% CI, 3.15–6.35; p < 0.001) were the only significant independent predictors on multivariate analysis, and none of the immunohistochemical markers were associated with poor DFS.

3. Discussion

In this study, we examined the expression of Nrf2, pNrf2, and Keap1 in HCCs, and we correlated the findings with the clinicopathological parameters and expression status of stemness, hypoxia, and EMT-related markers. The association between Nrf2 expression and poor prognosis was previously demonstrated in various solid malignancies, including in two meta-analysis studies [18,21]; however, the prognostic value of Nrf2/Keap1 in HCCs was only reported in two studies so far [16,22]. Zhang et al. demonstrated an association between Nrf2 expression in HCCs and decreased survival, poor histological differentiation, increased proliferation, and more frequent metastasis [22]. Although we observed a significant association between strong cytoplasmic Nrf2 expression in HCC with decreased OS and EMT-related marker expression, other clinicopathological factors such as proliferative activity were not significantly influenced by Nrf2 expression status. The different antibodies used may partly explain the differences; Zhang et al. described nucleocytoplasmic Nrf2 expression with predominantly nuclear localization in their study, while Nrf2 expression in the current study was cytoplasmic. As it is the nuclear form of Nrf2 that is recognized to be the active pro-tumorigenic form, we also examined nuclear pNrf2 expression in HCC, and we found that pNrf2-high HCCs were associated with poor prognosis, increased proliferative activity, and more frequent expression of markers associated with stemness, EMT, and hypoxia. The association between pNrf2 expression and decreased survival in HCCs is in line with a previous study by Chen et al., although we did not find the same expected positive correlation between pNrf2 and Nrf2 expression or an inverse correlation between pNrf2 and Keap1 expression that they observed [16]. It should be noted that Chen et al. reported a larger proportion of pNrf2-positive HCCs (50%) which could be explained by the different cut-offs; we evaluated pNrf2 using a quantitative approach and defined “pNrf2-high” HCCs (~25%) as those with strong nuclear labeling exceeding 10%, while the former study adopted a semiquantitative approach for both pNrf2 and Keap1, counting any labeling for pNrf2 and Keap1 (in the nucleus and cytoplasm, respectively) that was more than focal weak positive as “positive” [16]. Nevertheless, despite the differences in interpretation methods, it is interesting that pNrf2 expression was associated with CAIX, EpCAM, K19, and ezrin expression in this study.
A subset of HCCs that are morphologically compatible with HCC but express immunohistochemical markers associated with stemness (e.g., K19, EpCAM, and SAL-like protein 4 (SALL4)) recently received increasing attention, as these HCCs demonstrate more aggressive behavior compared to conventional HCCs, expression of EMT-related markers, poor prognosis, and resistance to locoregional and systemic treatment [23,24,25]. HCCs expressing stemness-related markers were shown to more frequently result in residual tumors after transarterial chemoembolization (TACE), compared to conventional HCCs, and the residual tumors after TACE often expressed stemness-related markers that were co-localized with CAIX expression, thus suggesting a relationship between TACE-induced hypoxia and stemness [26]. In this study, we found positive correlations between pNrf2 expression and hypoxia, stemness, and EMT-related markers, suggesting a possible role for Nrf2 in the link between tumor hypoxia, stemness, and aggressive behavior of HCC. pNrf2 may be upregulated in HCCs in response to oxidative stress, and pNrf2 may further enhance the aggressive behavior of hypoxic HCCs; however, functional studies would be necessary for validation.
Although Keap1 expression status was not correlated with clinicopathological factors or patient survival, Keap1-high HCCs more frequently demonstrated p53 expression and higher proliferative activity. In addition, a positive correlation was found between Keap1 and pNrf2 expression status. These results are somewhat counterintuitive, as Keap1 was described to be a protein that suppresses Nrf2 activation. Interestingly, when we examined the transcriptomics data provided by The Human Protein Atlas (https://www.proteinatlas.org/) and UALCAN (http://ualcan.path.uab.edu/), which are based on The Cancer Genome Atlas (TCGA) data, we found that KEAP1 messenger RNA (mRNA) expression was in fact higher in HCCs than in non-neoplastic liver tissues, and high KEAP1 mRNA was associated with poor prognosis in HCCs (Figure S2) [27,28]. An association between Keap1 expression and poor prognosis was also demonstrated in endometrial cancers, and the authors postulated that the overexpression of Keap1 may reflect Keap1 induction in the setting of oxidative stress [15]. On the other hand, when we analyzed the Keap1 expression status separately in CAIX-positive and CAIX-negative HCCs, we interestingly found that Keap1 expression was associated with more frequent portal vein invasion, higher proliferative activity, frequent p53 expression, and a tendency for decreased DFS in CAIX-negative HCCs, while Keap1 was not associated with any of the clinicopathological or immunohistochemical factors in the CAIX-positive HCCs. Instead, Nrf2 and/or pNrf2 were the determinants of aggressive behavior and poor survival in the setting of CAIX-positive HCCs. Thus, it is possible that, under hypoxic stress, Nrf2 is released from Keap1 to contribute to the aggressive behavior of HCC, while, under normoxic conditions, Keap1 is also a determinant of aggressive behavior. In addition, the association between strong cytoplasmic Nrf2 expression and decreased OS in this study suggests that, although the nuclear pNrf2 is the active form of Nrf2, strong overexpression of Nrf2 in the cytoplasm is also abnormal and has prognostic implications.
There are a few limitations to this study. Firstly, although we analyzed a large number of HCCs, this is a retrospective cohort study using formalin-fixed paraffin-embedded tissue samples and the results are based on immunohistochemical results. As all the HCCs enrolled in this study were surgically resected HCCs, it is possible that early well-differentiated HCCs—which are often treated by loco-regional modalities—were less likely to be included in this study, resulting in a biased cohort predominantly consisting of moderately to poorly differentiated HCCs. Functional studies would need to be performed to further evaluate the functional relationships among Nrf2/Keap1, hypoxia, and EMT. Secondly, as the majority of the HCCs in this study have an HBV-related etiology, it would be necessary to perform a similar analysis on an independent cohort of HCCs from regions where HBV is not as predominant.

4. Materials and Methods

4.1. Patients

A total of 285 cases of primary HCCs that were surgically resected between 2003 and 2012 at Seoul National University Bundang Hospital (Seongnam, Korea) were evaluated in this study. Clinicopathological information was obtained from the patients’ electronic medical records, pathology database, and a review of the archived glass slides, and it included sex, age at operation, underlying etiology (including hepatitis B virus (HBV), hepatitis C virus (HCV), alcohol and non-alcoholic fatty liver disease (NAFLD)), tumor size, multiplicity, tumor differentiation (Edmonson–Steiner grade), presence of microvascular and portal vein invasion, mitotic index, extrahepatic metastasis, and pathological T stages according to the American Joint Committee on Cancer TNM staging system (eighth edition). The mitotic index was expressed as the total number of mitotic figures in 10 high-power fields (HPFs). The clinical follow-up information, including the dates of surgery, disease recurrence, death, and most recent follow-up, was retrieved from the electronic medical records. Overall survival (OS) was defined as the number of months between surgery and patient death, and disease-free survival (DFS) was defined as the number of months between surgery and disease recurrence (local recurrence and/or distant metastasis). In addition, a separate independent tissue microarray cohort of 293 resected HCCs from Seoul National University Hospital was evaluated to further examine the prognostic significance of pNrf2 expression. This study was approved by the Institutional Review Board of Seoul National University Hospital and Seoul National University Bundang Hospital (joint IRB# H-1808-073-965), and patient consent was waived due to the retrospective nature of this study.

4.2. Tissue Samples and Immunohistochemistry

Tissue microarrays with tissue cores measuring 2 mm in diameter were constructed from formalin-fixed paraffin-embedded HCC tissue, and they were arranged in recipient tissue microarray blocks using a trephine apparatus (Superbiochips Laboratories, Seoul, Korea). Immunohistochemistry was performed on 4-μm-thick tissue sections using antibodies for pNrf2, Nrf2, Keap1, K19, EpCAM, E-cadherin, CAIX, p53, and Ki-67, using an automated platform (Ventana BenchMark GX, Ventana Medical Systems, Oro Valley, AZ, USA). Immunohistochemical staining for ezrin and uPAR was performed manually, with antigen retrieval using the citrate buffer (pH 6.0) and the Envision Kit (Dako, Glostrup, Denmark). Detailed information on the antibodies is summarized in Table 4.
For Nrf2 and Keap1, the intensity of cytoplasmic expression in the tumor cells was recorded as 0 (negative), 1+ (weak), 2+ (moderate), and 3+ (strong). Nrf2 expression was cytoplasmic, and strong positive (3+) staining was designated as “Nrf2-high” (overexpression). “Keap1-high” was defined as moderate to strong cytoplasmic staining in the tumor cells, and the staining intensity was similar to or more intense than that seen in non-neoplastic hepatocytes, cholangiocytes, and stromal cells. pNrf2 staining was evaluated in a quantitative manner, using the ImageJ software (National Institutes of Health, Bethesda, MD, USA); the percentage of nuclear pNrf2 labeling was calculated by dividing the number of positively stained tumor cell nuclei by the total number of tumor cell nuclei, in an HPF (400× magnification). “pNrf2-high” was defined as pNrf2 labeling in ≥10% of tumor cells. K19, EpCAM, CAIX, ezrin, and uPAR expression was designated positive when there was cytoplasmic (K19, CAIX, ezrin) or membranous (EpCAM and uPAR) expression in ≥5% of the tumor cells. Loss of E-cadherin expression was defined as complete absence of membranous E-cadherin expression in the tumor cells. p53 expression was defined as positive when strong nuclear labeling was present in more than 10% HCC cells. Two pathologists (K.L. and H.K.) independently evaluated the immunohistochemical stains and reached an agreement using a multihead microscope when there were discrepancies in the interpretation results.

4.3. Statistical Analysis

All data were analyzed using IBM SPSS statistics software version 25.0 (SPSS Inc., Chicago, IL, USA) and GraphPad Prism software version 7.0 (San Diego, CA, USA), and a p-value <0.05 was considered statistically significant. The associations between the immunohistochemical stain results and clinicopathological parameters were assessed using the χ2 test, Fisher exact test, and Student’s t-test, as deemed appropriate. OS and DFS were analyzed using the Kaplan–Meier method and the log-rank test. Multivariate analysis using the Cox regression model was performed for clinicopathological and immunohistochemical parameters with p < 0.05 on univariate analysis.

5. Conclusions

In summary, we demonstrate that pNrf2 expression in HCC was associated with increased proliferative activity, decreased survival, and more frequent expression of hypoxia, stemness, and EMT-related markers, suggesting a possible role for pNrf2 as a marker of aggressive behavior and poor prognosis in HCC. In addition, contrary to previous knowledge, strong cytoplasmic Nrf2 expression was also an independent predictor of decreased survival in HCC, and Keap1 expression was associated with aggressive HCC behavior in CAIX-negative tumors, suggesting that Keap1 expression in HCCs should be interpreted in the context of the metabolic status. Further prognostic studies in independent cohorts and functional studies would be required to validate the current findings.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/8/2128/s1. Figure S1: Kaplan-Meier survival curves demonstrating decreased overall survival (A) and disease-free survival (B) in pNrf2-high HCCs in an independent cohort of HCC (n = 293), Figure S2: Increased KEAP1 mRNA expression in HCCs compared to non-neoplastic liver tissue (A) and decreased overall survival in KEAP1-high HCCs (B) in the LIHC cohort from The Cancer Genome Atlas (graphs generated from UALCAN, http://ualcan.path.uab.edu/); Table S1. Clinicopathologic features according to Keap1 expression status in CAIX-positive and CAIX-negative hepatocellular carcinomas.

Author Contributions

Writing—original draft, K.L. (Kiryang Lee) and H.K.; writing—review and editing, K.L. (Kiryang Lee) and H.K.; conceptualization, K.L. (Kiryang Lee) and H.K.; resources, K.L. (Kiryang Lee), S.K., H.K., S.J.Y., Y.L. (Yangkyu Lee), H.L., Y.L. (Youngeun Lee), H.P., J.H.N., S.A., and K.L. (Kyoungbun Lee); formal analysis, K.L. (Kiryang Lee), S.K., and H.K.; software, K.L. (Kiryang Lee) and S.K.; visualization, K.L. (Kiryang Lee); supervision, H.K.; funding acquisition, H.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (NRF-2019R1A2C2010056, NRF-2016R1D1A1A09919042).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ferlay, J.; Colombet, M.; Soerjomataram, I.; Mathers, C.; Parkin, D.M.; Piñeros, M.; Znaor, A.; Bray, F. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int. J. Cancer 2019, 144, 1941–1953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fujimoto, A.; Furuta, M.; Totoki, Y.; Tsunoda, T.; Kato, M.; Shiraishi, Y.; Tanaka, H.; Taniguchi, H.; Kawakami, Y.; Ueno, M.; et al. Whole-genome mutational landscape and characterization of noncoding and structural mutations in liver cancer. Nat. Genet. 2016, 48, 500–509. [Google Scholar] [CrossRef] [PubMed]
  3. Guichard, C.; Amaddeo, G.; Imbeaud, S.; Ladeiro, Y.; Pelletier, L.; Ben Maad, I.; Calderaro, J.; Bioulac-Sage, P.; Letexier, M.; Degos, F.; et al. Integrated analysis of somatic mutations and focal copy-number changes identifies key genes and pathways in hepatocellular carcinoma. Nat. Genet. 2012, 44, 694. [Google Scholar] [CrossRef] [PubMed]
  4. Cancer Genome Atlas Research Network. Comprehensive and Integrative Genomic Characterization of Hepatocellular Carcinoma. Cell 2017, 169, 1327–1341. [Google Scholar] [CrossRef] [PubMed]
  5. Torbenson, M.S.; Ng, I.O.L.; Park, Y.N.; Roncalli, M.; Sakamoto, M. Hepatocellular carcinoma. In WHO Classification of Tumours Editorial Board, 5th ed.; Digestive System Tumours, Ed.; International Agency for Research on Cancer: Lyon, France, 2019; pp. 229–239. [Google Scholar]
  6. Zucman-Rossi, J.; Villanueva, A.; Nault, J.C.; Llovet, J.M. Genetic Landscape and Biomarkers of Hepatocellular Carcinoma. Gastroenterology 2015, 149, 1226–1239. [Google Scholar] [CrossRef] [Green Version]
  7. Sporn, M.B.; Liby, K.T. NRF2 and cancer: The good, the bad and the importance of context. Nat. Rev. Cancer 2012, 12, 564–571. [Google Scholar] [CrossRef]
  8. Suzuki, T.; Motohashi, H.; Yamamoto, M. Toward clinical application of the Keap1-Nrf2 pathway. Trends Pharm. Sci. 2013, 34, 340–346. [Google Scholar] [CrossRef]
  9. Chun, Y.; Ye-xiong, T.; Guang-zhen, Y.; Jian, Z.; Yu-fei, P.; Chen, L.; Jing, F.; Yao, C.; Zhi-wen, D.; Li-wei, D.; et al. Gankyrin has an antioxidative role through the feedback regulation of Nrf2 in hepatocellular carcinoma. J. Exp. Med. 2016, 213, 859–875. [Google Scholar]
  10. Kahroba, H.; Shirmohamadi, M.; Hejazi, M.S.; Samadi, N. The Role of Nrf2 signaling in cancer stem cells: From stemness and self-renewal to tumorigenesis and chemoresistance. Life Sci. 2019, 239, 116986. [Google Scholar] [CrossRef]
  11. Jaramillo, M.C.; Zhang, D.D. The emerging role of the Nrf2-Keap1 signaling pathway in cancer. Genes Dev. 2013, 27, 2179–2191. [Google Scholar] [CrossRef] [Green Version]
  12. Ngo, H.K.C.; Kim, D.H.; Cha, Y.N.; Na, H.K.; Surh, Y.J. Nrf2 Mutagenic Activation Drives Hepatocarcinogenesis. Cancer Res. 2017, 77, 4797–4808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Raghunath, A.; Sundarraj, K.; Arfuso, F.; Sethi, G.; Perumal, E. Dysregulation of Nrf2 in Hepatocellular Carcinoma: Role in Cancer Progression and Chemoresistance. Cancers 2018, 10, 481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, C.; Cheng, L.; Wu, H.; He, P.; Zhang, Y.; Yang, Y.; Chen, J.; Chen, M. Activation of the KEAP1-NRF2-ARE signaling pathway reduces oxidative stress in Hep2 cells. Mol. Med. Rep. 2018, 18, 2541–2550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ahtikoski, A.M.; Kangas, J.; Salonen, R.; Puistola, U.; Karihtala, P. Cytoplasmic Keap1 Expression Is Associated With Poor Prognosis in Endometrial Cancer. Anticancer. Res. 2019, 39, 585–590. [Google Scholar] [CrossRef]
  16. Chen, J.; Yu, Y.; Ji, T.; Ma, R.; Chen, M.; Li, G.; Li, F.; Ding, Q.; Kang, Q.; Huang, D.; et al. Clinical implication of Keap1 and phosphorylated Nrf2 expression in hepatocellular carcinoma. Cancer Med. 2016, 5, 2678–2687. [Google Scholar] [CrossRef]
  17. Cho, H.Y.; Kim, K.; Kim, Y.B.; Kim, H.; No, J.H. Expression Patterns of Nrf2 and Keap1 in Ovarian Cancer Cells and their Prognostic Role in Disease Recurrence and Patient Survival. Int. J. Gynecol. Cancer 2017, 27, 412–419. [Google Scholar] [CrossRef]
  18. Guo, Y.; Shen, L. Overexpression of NRF 2 is correlated with prognoses of patients with malignancies: A meta-analysis. Thorac. Cancer 2017, 8, 558–564. [Google Scholar] [CrossRef]
  19. Isohookana, J.; Haapasaari, K.M.; Soini, Y.; Karihtala, P. Keap1 expression has independent prognostic value in pancreatic adenocarcinomas. Diagn. Pathol. 2015, 10, 28. [Google Scholar] [CrossRef] [Green Version]
  20. Lingzhi, P.; Yaoxing, D.; Yongbin, Z.; Diande, Z.; Yanghui, W.; Jibin, Y.; Mingxu, D. Expression of Nrf2 and NQO1 in human gastric cancer and their clinical significance. Int. J. Clin. Exp. Pathol. 2016, 9, 1635–1643. [Google Scholar]
  21. Wang, L.; Zhang, C.; Qin, L.; Xu, J.; Li, X.; Wang, W.; Kong, L.; Zhou, T.; Li, X. The prognostic value of NRF2 in solid tumor patients: A meta-analysis. Oncotarget 2018, 9, 1257–1265. [Google Scholar] [CrossRef] [Green Version]
  22. Zhang, M.; Zhang, C.; Zhang, L.; Yang, Q.; Zhou, S.; Wen, Q.; Wang, J. Nrf2 is a potential prognostic marker and promotes proliferation and invasion in human hepatocellular carcinoma. BMC Cancer 2015, 15, 531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Olivier, G.; Mina, K.; Johannes, B.; Bart, S.; Carl, J.; Francesca, L.; Aezam, K.; Jasper, W.; Léon, K.; Anne, D.; et al. Keratin 19: A key role player in the invasion of human hepatocellular carcinomas. Gut 2014, 63, 674–685. [Google Scholar]
  24. Kawai, T.; Yasuchika, K.; Ishii, T.; Katayama, H.; Yoshitoshi, E.Y.; Ogiso, S.; Kita, S.; Yasuda, K.; Fukumitsu, K.; Mizumoto, M.; et al. Keratin 19, a Cancer Stem Cell Marker in Human Hepatocellular Carcinoma. Clin. Cancer Res. 2015, 21, 3081–3091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kim, H.; Choi, G.H.; Na, D.C.; Ahn, E.Y.; Kim, G.I.; Lee, J.E.; Cho, J.Y.; Yoo, J.E.; Choi, J.S.; Park, Y.N. Human hepatocellular carcinomas with “Stemness”-related marker expression: Keratin 19 expression and a poor prognosis. Hepatology 2011, 54, 1707–1717. [Google Scholar] [CrossRef] [PubMed]
  26. Rhee, H.; Nahm, J.H.; Kim, H.; Choi, G.H.; Yoo, J.E.; Lee, H.S.; Koh, M.J.; Park, Y.N. Poor outcome of hepatocellular carcinoma with stemness marker under hypoxia: Resistance to transarterial chemoembolization. Mod. Pathol. 2016, 29, 1038–1049. [Google Scholar] [CrossRef]
  27. Uhlen, M.; Fagerberg, L.; Hallstrom, B.M.; Lindskog, C.; Oksvold, P.; Mardinoglu, A.; Sivertsson, A.; Kampf, C.; Sjöstedt, E.; Asplund, A.; et al. Proteomics. Tissue-based map of the human proteome. Science 2015, 347, 1260419. [Google Scholar] [CrossRef]
  28. Chandrashekar, D.S.; Bashel, B.; Balasubramanya, S.A.H.; Creighton, C.J.; Ponce-Rodriguez, I.; Chakravarthi, B.V.; Varambally, S. UALCAN: A Portal for Facilitating Tumor Subgroup Gene Expression and Survival Analyses. Neoplasia 2017, 19, 649–658. [Google Scholar] [CrossRef]
Figure 1. Expression of phosphorylated nuclear factor E2-related factor2 (pNrf2), Nrf2, and Keap1 in hepatocellular carcinomas. (a) Nuclear pNrf2 expression in 63% of tumor cells; (b) high cytoplasmic Nrf2 expression in tumor cells; (c) high cytoplasmic Keap1 expression in tumor cells. Original magnification ×200 (a), ×400 (b,c).
Figure 1. Expression of phosphorylated nuclear factor E2-related factor2 (pNrf2), Nrf2, and Keap1 in hepatocellular carcinomas. (a) Nuclear pNrf2 expression in 63% of tumor cells; (b) high cytoplasmic Nrf2 expression in tumor cells; (c) high cytoplasmic Keap1 expression in tumor cells. Original magnification ×200 (a), ×400 (b,c).
Cancers 12 02128 g001
Figure 2. Comparison of the Ki-67 labeling indices (A,D,G), mitotic indices (B,E,H), and immunohistochemical characteristics (C,F,I) of hepatocellular carcinomas according to pNrf2 (AC), Nrf2 (DF), and Keap1 (GI) expression status; * p < 0.05, ** p < 0.01; HPF: high-power field.
Figure 2. Comparison of the Ki-67 labeling indices (A,D,G), mitotic indices (B,E,H), and immunohistochemical characteristics (C,F,I) of hepatocellular carcinomas according to pNrf2 (AC), Nrf2 (DF), and Keap1 (GI) expression status; * p < 0.05, ** p < 0.01; HPF: high-power field.
Cancers 12 02128 g002
Figure 3. Kaplan–Meier survival curves comparing the overall survival (AC) and disease-free survival (DF) of hepatocellular carcinomas according to pNrf2 (A,D), Nrf2 (B,E), and Keap1 (C,F) expression status.
Figure 3. Kaplan–Meier survival curves comparing the overall survival (AC) and disease-free survival (DF) of hepatocellular carcinomas according to pNrf2 (A,D), Nrf2 (B,E), and Keap1 (C,F) expression status.
Cancers 12 02128 g003
Table 1. Summary of the clinicopathological features of hepatocellular carcinomas (n = 285).
Table 1. Summary of the clinicopathological features of hepatocellular carcinomas (n = 285).
ParametersFrequency (%) or Median (Range)
Sex, n (%)
Male216 (75.8)
Female69 (24.2)
Age at operation, median (range, years)63 (29–87)
Etiology, n (%)
Hepatitis B virus (HBV)202 (70.9)
Hepatitis C virus (HCV)25 (8.8)
HBV + HCV1 (0.3)
Alcohol14 (4.9)
Non-alcoholic fatty liver disease12 (4.2)
Uncertain etiology31 (10.9)
Tumor size, median (range, cm)3.2 (0.9–17)
Multiplicity present, n (%)49 (17.2)
Edmondson–Steiner grade, n (%)
Grade I3 (1.0)
Grade II77 (27.0)
Grade III171 (60.0)
Grade IV34 (12.0)
Microvascular invasion present, n (%)113 (39.6)
Portal vein invasion present, n (%)19 (6.7)
Pathologic T category (AJCC 8th edition), n (%)
pT1a37 (13.0)
pT1b101 (35.4)
pT2112 (39.3)
pT311 (3.9)
pT424 (8.4)
Extrahepatic metastasis present, n (%)55 (19.3)
Cirrhosis in background liver present, n (%)152 (533.3)
Recurrence on follow-up present, n (%)161 (56.5)
Status at last follow-up, n (%)
Alive128 (44.9)
Deceased of disease44 (15.4)
Deceased of other cause10 (3.5)
Follow up loss103 (36.2)
Overall survival, median (range, months) 66.5 (0–163)
Disease-free survival, median (range, months) 32.5 (0–163)
AJCC, American Joint Committee on Cancer.
Table 2. Summary of the clinicopathological and immunohistochemical features of hepatocellular carcinomas according to pNrf2, Nrf2 and Keap1 status (n = 285), n (%).
Table 2. Summary of the clinicopathological and immunohistochemical features of hepatocellular carcinomas according to pNrf2, Nrf2 and Keap1 status (n = 285), n (%).
ParameterpNrf2p-ValueNrf2p-ValueKeap1p-Value
Low pNrf2High pNrf2Low Nrf2High Nrf2Low Keap1High Keap1
Frequency212 (74.3)73 (25.6)-237 (83.2)48 (16.8)-124 (43.5)161 (56.5)-
Sex (male/female)159 (75.0)/53 (25.0)57 (78.1)/16 (21.9)0.638179 (75.5)/58 (24.5)37 (77.1)/11 (22.9)1.00094 (75.8)/30 (24.2)122 (75.8)/39 (24.2)1.000
Age ≥60 years95 (44.8)28 (38.4)0.411100 (42.2)23 (47.9)0.54255 (44.7)68 (42.2)0.809
HBV-related etiology150 (70.8)53 (72.6)0.881167 (70.5)36 (75.0)0.60286 (69.4)117 (72.7)0.598
Size (cm, mean ± SD)4.2 ± 3.04.0 ± 2.40.2744.4 ± 2.93.1 ± 1.80.0054.4 ± 3.13.9 ± 2.50.105
Multiplicity37 (17.5)12 (16.4)1.00041 (17.3)8 (16.7)1.00023 (18.5)26 (16.1)0.636
E–S grade III and IV155 (73.1)50 (68.5)0.454173 (73.0)32 (66.7)0.38292 (74.2)113 (70.2)0.507
Microvascular invasion78 (36.8)35 (47.9)0.09896 (40.5)17 (35.4)0.62847 (37.9)66 (41.0)0.627
Portal vein invasion12 (5.7)7 (9.6)0.27817 (7.2)2 (4.2)0.75012 (9.7)7 (4.3)0.094
pT category 0.146 0.323 0.166
pT1a31 (14.6)6 (8.2) 33 (13.9)4 (8.3) 18 (14.5)19 (11.8)
pT1b75 (35.4)26 (35.6) 83 (35.0)18 (37.5) 46 (37.1)55 (34.2)
pT276 (35.8)36 (49.3) 89 (37.6)23 (47.9) 41 (33.1)71 (44.1)
pT310 (4.7)1 (1.4) 11 (4.6)0 (0) 4 (3.2)7 (4.3)
pT420 (9.4)4 (5.5) 21 (8.9)3 (6.3) 15 (12.1)9 (5.6)
Recurrence118 (55.7)43 (58.9)0.784133 (56.1)28 (58.3)0.74864 (51.6)97 (60.2)0.183
Extrahepatic metastasis37 (17.6)18 (24.7)0.22950 (21.2)5 (9.1)0.10924 (19.5)31 (19.4)1.000
Underlying cirrhosis108 (52.9)44 (60.3)0.338125 (54.1)27 (58.7)0.62865 (55.1)87 (54.7)1.000
K19 expression38 (17.9)16 (21.9)0.49042 (17.7)12 (25.0)0.23318 (14.5)36 (22.4)0.127
EpCAM expression84 (39.6)41 (56.2)0.020103 (43.5)22 (45.8)0.87349 (39.5)76 (47.2)0.229
CAIX expression48 (23.4)27 (37.5)0.03055 (24.0)20 (41.7)0.01931 (25.4)44 (28.4)0.589
Ezrin expression78 (36.8)43 (59.7)0.00190 (38.1)31 (64.6)0.00148 (38.7)73 (45.6)0.277
uPAR expression58 (27.5)26 (36.1)0.18070 (29.7)14 (29.8)1.00029 (23.4)55 (34.6)0.049
E-cadherin loss128 (60.4)56 (76.7)0.015153 (64.6)31 (64.6)1.00081 (65.3)103 (64.0)0.901
p53 overexpression40 (20.5)26 (37.1)0.00953 (24.1)13 (28.9)0.57019 (16.1)47 (32.0)0.004
Mitotic index (/10 HPF, mean ± SD)10.1 ± 13.417.4 ± 20.0<0.00111.9 ± 15.612.2 ± 16.10.62310.5 ± 13.613.1 ± 16.90.056
Ki-67 labeling index (%, mean ± SD)5.1 ± 7.512.9 ± 12.2<0.0017.2 ± 9.66.7 ± 9.60.9355.7 ± 8.58.1 ± 10.20.035
HR: hazard ratio, CI: confidence interval, HBV: hepatitis B virus, E–S grade: Edmondson–Steiner grade, K19: keratin 19, EpCAM: epithelial cell adhesion molecule, CAIX: carbonic anhydrase IX; uPAR: urokinase-type plasminogen activator receptor, HPF: high-power field.
Table 3. Summary of the survival analysis results (n = 285).
Table 3. Summary of the survival analysis results (n = 285).
ParameterOverall SurvivalDisease-Free Survival
Univariate AnalysisMultivariate AnalysisUnivariate AnalysisMultivariate Analysis
ParametersHR (95% CI)p-ValueHR (95% CI)p-ValueHR (95% CI)p-ValueHR (95% CI)p-Value
Sex (male/female)0.742 (0.345–1.597)0.446--0.702 (0.480–1.029)0.069--
Age (≥60/<60 years)0.960 (0.522–1.765)0.894--0.848 (0.621–1.157)0.299--
Underlying HBV1.160 (0.586–2.296)0.671--1.129 (0.801–1.590)0.488--
Size (≥5cm/<5cm)2.400 (1.315–4.382)0.0041.298 (0.650–2.594)0.4601.418 (1.012–1.988)0.0431.085 (0.765–1.540)0.647
Multiplicity2.812 (1.488–5.312)0.0012.156 (1.093–4.253)0.0271.991 (1.375–2.884)<0.0011.487 (1.015–2.180)0.042
E–S grade III and IV0.951 (0.497–1.818)0.879--1.095 (0.778–1.542)0.601--
Microvascular invasion1.477 (0.813–2.682)0.200--1.308 (0.962–1.779) 0.087--
Portal vein invasion3.289 (1.386–7.804)0.0072.435 (0.980–6.052)0.0551.573 (0.853–2.903)0.147--
Extrahepatic metastasis4.458 (2.445–8.130)<0.0013.256 (1.742–6.083)<0.0014.799 (3.411–6.753)<0.0014.475 (3.152–6.354)<0.001
Underlying cirrhosis1.282 (0.702–2.340)0.419--1.031 (0.758–1.404)0.844--
K191.860 (0.957–3.614)0.067--1.209 (0.826–1.771)0.329--
EpCAM1.727 (0.951-3.137)0.073--0.889 (0.613–1.292)0.538--
CAIX1.854 (1.015–3.386)0.0441.402 (0.738–2.664)0.3021.043 (0.743–1.464)0.809--
Ezrin3.160 (1.674–5.966)<0.0012.893 (1.527–5.483)0.0011.090 (0.766–1.552)0.632--
uPAR1.492 (0.812–2.743)0.198 1.058 (0.756–1.480)0.742
E-cadherin loss1.037 (0.556–1.936)0.909 0.853 (0.622–1.169)0.323
p53 overexpression1.777 (0.961–3.288)0.067--1.219 (0.854–1.740)0.276--
Nrf23.329 (1.738–6.374)<0.0014.151 (2.025–8.508)<0.0011.315 (0.879–1.968)0.183--
pNrf21.894 (1.037–3.458)0.0381.181 (0.950–3.450)0.0710.933 (0.661–1.317)0.693--
Keap11.061 (0.582–1.937)0.846--1.290 (0.944–1.762)0.109--
HBV: hepatitis B virus, HR: hazard ratio, CI: confidence interval, SD: standard deviation, E–S grade: Edmondson–Steiner grade, K19: keratin 19, EpCAM: epithelial cell adhesion molecule, CAIX: carbonic anhydrase IX, uPAR: urokinase-type plasminogen activator receptor, HPF: high-power field.
Table 4. Summary of the antibodies used in this study.
Table 4. Summary of the antibodies used in this study.
AntibodyDescriptionSourceDilutionMethod
Nrf2Rabbit monoclonal (EP1808Y)Abcam 1:50Autostainer
(Cambridge, UK)
pNrf2Rabbit monoclonal (EP1809Y)Abcam 1:100Autostainer
(Cambridge, UK)
Keap1Rabbit polyclonalProteintech 1:300Autostainer
(Manchester, UK)
Keratin 19 (K19)Mouse monoclonal (RCK108)Dako1:200Autostainer
(Glostrup, Denmark)
EpCAMMouse monoclonal (VU-1D9)Millipore1:1500Autostainer
(Temecula, CA, USA)
EzrinMouse monoclonal (3C12)Abcam 1:100Manual (pH 6.0)
(Cambridge, UK)
uPARMouse monoclonal (R-4)Abcam1:40Manual (pH 6.0)
(Cambridge, UK)
E-cadherinMouse monoclonal (36B5)Novocastra1:100Autostainer
(Newcastle, UK)
CAIX Rabbit polyclonalAbcam 1:500Autostainer
(Cambridge, UK)
p53Mouse monoclonal (DO-7)Dako1:1000Autostainer
(Glostrup, Denmark)
Ki-67Mouse monoclonal (MIB-1)Dako1:100Autostainer
(Glostrup, Denmark)

Share and Cite

MDPI and ACS Style

Lee, K.; Kim, S.; Lee, Y.; Lee, H.; Lee, Y.; Park, H.; Nahm, J.H.; Ahn, S.; Yu, S.J.; Lee, K.; et al. The Clinicopathological and Prognostic Significance of Nrf2 and Keap1 Expression in Hepatocellular Carcinoma. Cancers 2020, 12, 2128. https://doi.org/10.3390/cancers12082128

AMA Style

Lee K, Kim S, Lee Y, Lee H, Lee Y, Park H, Nahm JH, Ahn S, Yu SJ, Lee K, et al. The Clinicopathological and Prognostic Significance of Nrf2 and Keap1 Expression in Hepatocellular Carcinoma. Cancers. 2020; 12(8):2128. https://doi.org/10.3390/cancers12082128

Chicago/Turabian Style

Lee, Kiryang, Seunghye Kim, Yangkyu Lee, Hyejung Lee, Youngeun Lee, Hyunjin Park, Ji Hae Nahm, Soomin Ahn, Su Jong Yu, Kyoungbun Lee, and et al. 2020. "The Clinicopathological and Prognostic Significance of Nrf2 and Keap1 Expression in Hepatocellular Carcinoma" Cancers 12, no. 8: 2128. https://doi.org/10.3390/cancers12082128

APA Style

Lee, K., Kim, S., Lee, Y., Lee, H., Lee, Y., Park, H., Nahm, J. H., Ahn, S., Yu, S. J., Lee, K., & Kim, H. (2020). The Clinicopathological and Prognostic Significance of Nrf2 and Keap1 Expression in Hepatocellular Carcinoma. Cancers, 12(8), 2128. https://doi.org/10.3390/cancers12082128

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop