AGR2 is clinically relevant in HCC
We retrieved two datasets, including the sorafenib-resistant dataset for Huh7 cells (GSE94550, [9]) and the Roessler Liver microarray dataset (GSE14520, [10]), and utilized > 2-fold change as the criterion for selection of sorafenib-modulated molecules in sorafenib-resistant Huh7 cells compared to parental cells (GSE94550) and > 1.2-fold as the criterion for choosing oncogenes related to survival in the bottom 25% vs. the top 25% of HCC patients from the Roessler Liver microarray dataset (GSE14520) to intersect the potential candidates for further study. According to the above analysis, the selected genes were validated and filtered more stringently and applied to evaluate highly significant molecules. After retrieving these datasets, 545 upregulated genes (sorafenib resistance vs. control > 2-fold) and 609 downregulated genes (sorafenib resistance vs. control < 2-fold) were identified in the sorafenib-resistant database in hepatoma Huh7 cells (GSE94550). In the Roessler Liver microarray (GSE14520), we used > 1.2-fold as the criterion for choosing oncogenes related to survival in the bottom 25% vs. the top 25% of HCC patients, and only 12 dysregulated genes were identified (Fig. 1A). Furthermore, we found that there were 4 potential candidates identified in the two intersecting microarray datasets, including neurotensin (NTS), AGR2, alpha-fetoprotein (AFP) and meprin A, alpha (MEP1A) (Fig. 1A, B). Previously, the AGR2 protein was demonstrated to be a part of the PDI family of ER proteins that mediate the formation of disulfide bonds and catalyze protein folding [32]. Moreover, AGR2 is highly expressed in numerous cancer types, including liver cancer [16]. Hrstka et al. showed that AGR2 expression can be used as a marker to predict poor prognosis in breast cancer via clinical studies [15]. Hence, we suggest that AGR2 may be a potential candidate target in sorafenib-treated HCC. Interestingly, AGR2 is the major gene highly correlated with the survival rate and sorafenib resistance in liver cancer; however, the relationship between AGR2 and sorafenib treatment in HCC has not been demonstrated. Therefore, AGR2 was selected for further study to investigate its molecular mechanism associated with sorafenib resistance and its physiological role in HCC.
The clinical significance of AGR2 expression in the Roessler Liver database (GSE14520) and our collected cohort was analyzed (Fig. 1C-F). Patients with lower expression (last 40%) of AGR2 had a better OS rate (log-rank P < 0.05; AGR2 high (top 40%): standard error, 2.738; 95% CI, 36.196–46.927; AGR2 low (last 40%): standard error, 2.432; 95% CI, 45.322–54.857) and RFS rate (log-rank P < 0.05; AGR2 high (top 40%): standard error, 2.766; 95% CI, 27.398–38.243; AGR2 low (last 40%): standard error, 2.645; 95% CI, 36.365–46.734) (Fig. 1C and D). The Roessler Liver microarray dataset (GSE14520) provides detailed clinical information; therefore, the correlation between AGR2 and various clinical parameters was analyzed statistically to define the role of sorafenib-regulated AGR2 in HCC progression (Tables 1–3). High AGR2 expression was significantly correlated with high AFP and ALT levels, a high predicted risk metastasis signature score (Fig. 1E, F and G), large main tumor size and more advanced pathological stages of HCC (Fig. 1H and I). A high AGR2 level in HCC was correlated with a high serum AFP level (P = 0.002) (Table 1). Univariate analysis showed that male sex (P = 0.009), tumor size (P = 0.045), CLIP score (P = 0.002), BCLC staging (P < 0.001), AJCC stage (P < 0.001), and high AGR2 level (P = 0.003) were significant predictors of worse RFS (Table 2). Multivariate analysis showed that male sex (P = 0.023, HR = 2.124, CI = 1.112–4.058), BCLC staging (P = 0.021, HR = 1.653, CI = 1.077–2.535), AJCC staging (P = 0.019, HR = 1.607, CI = 1.083–2.386) and a high AGR2 level (P = 0.010, HR = 1.572, CI = 1.112–2.224) were independently associated with RFS (Table 2). For OS, univariate analysis showed that cirrhosis (P = 0.023), AFP level (P = 0.011), tumor size (P = 0.001), CLIP score (P < 0.001), BCLC staging (P < 0.001), AJCC stage (P < 0.001), and high AGR2 level (P = 0.002) were significant predictors of worse OS (Table 3). Multivariate analysis showed that cirrhosis (P = 0.034, HR = 4.563, CI = 1.122–18.555), BCLC staging (P < 0.001, HR = 2.928, CI = 1.908–4.494), and a high AGR2 level (P = 0.008, HR = 1.735, CI = 1.154–2.609) were independently associated with OS (Table 3).
Table 1
Association of AGR2 level (Roessler liver array) with clinicopathologic indicators of hepatocellular carcinoma
Factors | Group | | AGR2 (mean ± SE) | P | |
Age | < 60 years | | 4.3208 ± 0.1474 | 0.900 | |
| ≥ 60 years | | 4.3935 ± 0.3107 | | |
Sex | Male | | 4.3334 ± 0.1422 | 0.926 | |
| Female | | 4.3525 ± 0.3850 | | |
Cirrhosis | Absent | | 3.6259 ± 0.3076 | 1.162 | |
| Present | | 4.3963 ± 0.1415 | | |
Serum AFP | < 300 | | 3.8669 ± 0.1446 | 0.002* | |
(ng/ml) | ≥ 300 | | 4.8537 ± 0.2261 | | |
Tumor size | < 5 cm | | 4.1165 ± 0.1462 | 0.055 | |
| ≥ 5 cm | | 4.7315 ± 0.2591 | | |
CLIP | 0–1 | | 4.2594 ± 0.1429 | 0.420 | |
| ≥ 2 | | 4.6445 ± 0.3413 | | |
BCLC | 0–1 | | 4.1908 ± 1.9342 | 0.094 | |
| ≥ 2 | | 4.8530 ± 2.4513 | | |
AJCC stage | I | | 4.1535 ± 1.8291 | 0.381 | |
| ≥II | | 4.4955 ± 2.2573 | | |
* P < 0.05. CLIP, Cancer of the Liver Italian Program score; BCLC, Barcelona Clinic Liver Cancer staging; AJCC, American Joint Committee on Cancer 2017; AFP, alpha-fetoprotein. |
Table 2
Prognostic significance of clinicopathologic indicators and AGR2 for recurrence-free survival in the Roessler liver array.
| RFS univariate | RFS multivariate |
Factor | Group | HR | 95% CI | P | | HR | 95% CI | P |
Age | < 60/≥60 years | 0.952 | 0.628–1.443 | 0.817 | | | | |
Sex | Female/Male | 2.359 | 1.238–4.493 | 0.009* | | 2.124 | 1.112–4.058 | 0.023* |
Cirrhosis | -/+ | 2.003 | 0.936–4.287 | 0.074 | | | | |
Serum AFP | < 300/≥300 ng/ml | 1.314 | 0.937–1.842 | 0.113 | | | | |
Tumor size | < 5/≥5 cm | 1.424 | 1.008–2.012 | 0.045* | | | | NS |
CLIP | 0–1/≥2 | 1.872 | 1.267–2.766 | 0.002* | | | | NS |
BCLC | 0–1/≥2 | 2.432 | 1.670–3.543 | < 0.001* | | 1.653 | 1.077–2.535 | 0.021* |
AJCC stage | I / ≥II | 2.063 | 1.405–3.029 | < 0.001* | | 1.607 | 1.083–2.386 | 0.019* |
AGR2 | Low/High | 1.691 | 1.202–2.378 | 0.003* | | 1.572 | 1.112–2.224 | 0.010* |
*P < 0.05. RFS, recurrence-free survival; CLIP, Cancer of the Liver Italian Program score; BCLC, Barcelona Clinic Liver Cancer staging; AJCC, American Joint Committee on Cancer 2017; AFP, alpha-fetoprotein. |
Table 3
Prognostic significance of clinicopathologic indicators and AGR2 for overall survival in the Roessler liver array.
| RFS univariate | RFS multivariate |
Factor | Group | HR | 95% CI | P | | HR | 95% CI | P |
Age | < 60/≥60 years | 0.990 | 0.972–1.008 | 0.990 | | | | |
Sex | Female/Male | 1.858 | 0.901–3.833 | 0.094 | | | | |
Cirrhosis | -/+ | 5.093 | 1.255–20.671 | 0.023* | | 4.563 | 1.122–18.555 | 0.034* |
Serum AFP | < 300/≥300 ng/ml | 1.686 | 1.126–2.527 | 0.011* | | | | |
Tumor size | < 5/≥5 cm | 1.960 | 1.309–2.933 | 0.001* | | | | |
CLIP | 0–1/≥2 | 2.811 | 1.832–4.313 | < 0.001* | | | | |
BCLC | 0–1/≥2 | 3.176 | 2.081–4.846 | < 0.001* | | 2.928 | 1.908–4.494 | < 0.001* |
AJCC stage | I / ≥II | 2.278 | 1.483-3.500 | < 0.001* | | | | |
AGR2 | Low/high | 1.919 | 1.282–2.873 | 0.002* | | 1.735 | 1.154–2.609 | 0.008* |
*P < 0.05. OS, overall survival; CLIP, Cancer of the Liver Italian Program score; BCLC, Barcelona Clinic Liver Cancer staging; AJCC, American Joint Committee on Cancer 2017; AFP, alpha-fetoprotein. |
Moreover, our HCC specimen cohort was also analyzed. We utilized qRT–PCR to examine the levels of AGR2, and a median level of 15.3 (39-△Ct) was defined as the cutoff to divide the HCC specimens into high and low expression of AGR2. Similar results were observed; high levels of AGR2 were related to significantly worse OS and RFS rates (Fig. 1J and K). Furthermore, the correlation between AGR2 expression and various clinical parameters from our collected cohort was also analyzed (Table 4–6). A higher AGR2 level in HCC was observed in female patients (P = 0.007) (Table 4). A high AGR2 level (39-△Ct ≥ 10.8) in HCC was significantly associated with worse OS (P = 0.016) (Fig. 1J) and RFS (P = 0.045) (Fig. 1K). Univariate analysis showed that cirrhosis (P = 0.022), AFP level (P = 0.033), vascular invasion (P = 0.003), AJCC stage (P = 0.003), and high AGR2 level (P = 0.047) were significant predictors of worse RFS (Table 5). Multivariate analysis showed that cirrhosis (P = 0.023, HR = 1.600, CI = 1.067–2.401), vascular invasion (P = 0.005, HR = 1.825, CI = 1.198–2.780), and a high AGR2 level (P = 0.043, HR = 1.662, CI = 1.017–2.716) were independently associated with RFS (Table 5). For OS, univariate analysis showed that cirrhosis (P = 0.020), AFP level (P = 0.012), vascular invasion (P < 0.001), AJCC stage (P = 0.007), and high AGR2 level (P = 0.018) were significant predictors of worse OS (Table 6). Multivariate analysis showed that cirrhosis (P = 0.017, HR = 1.655, CI = 1.093–2.505), vascular invasion (P = 0.001, HR = 2.200, CI = 1.392–3.476), and a high AGR2 level (P = 0.015, HR = 1.945, CI = 1.138–3.324) were independently associated with OS (Table 6).
Table 4
Association of AGR2 expression with clinicopathologic indicators of hepatocellular carcinoma
Factors | Group | | AGR2 (mean ± SE) | P | |
Age | < 60 years | | 15.8713 ± 0.7077 | 0.659 | |
| ≥ 60 years | | 16.2616 ± 0.4603 | | |
Sex | Male | | 15.4998 ± 0.4755 | 0.007* | |
| Female | | 17.6959 ± 0.5801 | | |
Hepatitis viral | Absent | | 15.9186 ± 0.6085 | 0.812 | |
infection | Present | | 16.2114 ± 0.4827 | | |
Cirrhosis | Absent | | 16.1498 ± 0.4687 | 0.989 | |
| Present | | 16.0949 ± 0.6798 | | |
Serum AFP | < 200 | | 15.8232 ± 0.5254 | 0.403 | |
(ng/ml) | ≥ 200 | | 16.5543 ± 0.5651 | | |
Tumor | W | | 14.4390 ± 1.5358 | 0.240 | |
differentiation | M-P | | 16.2353 ± 0.3985 | | |
Tumor size | < 5 cm | | 15.7190 ± 0.5735 | 0.482 | |
| ≥ 5 cm | | 16.3917 ± 0.5167 | | |
Vascular | Absent | | 15.5317 ± 0.6740 | 0.198 | |
invasion | Present | | 16.5114 ± 0.4624 | | |
AJCC stage | I | | 15.3584 ± 0.8662 | 0.241 | |
| ≥II | | 16.3968 ± 0.4242 | | |
* P < 0.05. Tumor differentiation by WHO; AJCC, American Joint Committee on Cancer 2017; AFP, alpha-fetoprotein. |
Table 5
Prognostic significance of clinicopathologic indicators and AGR2 for recurrence-free survival in the clinical cohort.
| RFS univariate | RFS multivariate |
Factor | Group | HR | 95% CI | P | | HR | 95% CI | P |
Age | < 60/≥60 years | 1.079 | (0.713–1.632) | 0.720 | | | | |
Sex | Male/female | 1.059 | (0.672–1.670) | 0.804 | | | | |
Viral infection | | | | 0.248 | | | | |
| No/B | 1.667 | (0.977–2.844) | 0.061 | | | | |
| No / C | 1.274 | (0.739–2.915) | 0.384 | | | | |
| No / B + C | 1.647 | (0.835–3.262) | 0.152 | | | | |
Cirrhosis | -/+ | 1.602 | (1.069–2.401) | 0.022* | | 1.600 | (1.067–2.401) | 0.023* |
Serum AFP | < 200/≥200 ng/ml | 1.548 | (1.037–2.312) | 0.033* | | | | NS |
Differentiation | W/M-P | 3.053 | (0.964–9.670) | 0.058 | | | | |
Tumor size | < 5/≥5 cm | 1.161 | (0.764–1.765) | 0.484 | | | | |
Vascular invasion | -/+ | 1.899 | (1.248–2.890) | 0.003* | | 1.825 | (1.198–2.780) | 0.005* |
AJCC stage | I / ≥II | 2.100 | (1.282–3.440) | 0.003* | | | | NS |
AGR2 | Low/high | 1.638 | (1.006–2.670) | 0.047* | | 1.662 | (1.017–2.716) | 0.043* |
*P < 0.05. DFS, disease-free survival; Tumor differentiation according to WHO system; AFP, alpha-fetoprotein; AJCC, American Joint Committee on Cancer 2017 |
Table 6
Prognostic significance of clinicopathologic indicators and AGR2 for overall survival in the clinical cohort.
| OS univariate | OS multivariate |
Factor | Group | HR | 95% CI | P | | HR | 95% CI | P |
Age | < 60/≥60 years | 0.975 | (0.638–1.490) | 0.907 | | | | |
Sex | Male/female | 1.407 | (0.906–2.186) | 0.936 | | | | |
Viral infection | | | | 0.281 | | | | |
| No/B | 1.546 | (0.890–2.685) | 0.122 | | | | |
| No / C | 1.014 | (0.573–1.794) | 0.962 | | | | |
| No / B + C | 1.433 | (0.704–2.916) | 0.322 | | | | |
Cirrhosis | -/+ | 1.632 | (1.080–2.466) | 0.020* | | 1.655 | (1.093–2.505) | 0.017* |
Serum AFP | < 200/≥200 ng/ml | 1.692 | (1.122–2.551) | 0.012* | | | | NS |
Differentiation | W/M-P | 1.330 | (0.851–2.079) | 0.210 | | | | |
Tumor size | < 5/≥5 cm | 1.101 | (0.720–1.683) | 0.658 | | | | |
Vascular invasion | -/+ | 2.302 | (1.458–3.635) | < 0.001* | | 2.200 | (1.392–3.476) | 0.001* |
AJCC stage | I/≥II | 2.076 | (1.224–3.521) | 0.007* | | | | NS |
AGR2 | Low/high | 1.894 | (1.114–3.221) | 0.018* | | 1.945 | (1.138–3.324) | 0.015* |
*P < 0.05. OS, overall survival; Tumor differentiation according to WHO system; AFP, alpha-fetoprotein; AJCC, American Joint Committee on Cancer 2017 |
Moreover, immunohistochemical staining was utilized to examine the expression of AGR2 in HCC tissues, and the results indicated that AGR2 was highly expressed in tumor tissues compared to normal tissues (Fig. 1L, M). Collectively, based on the evidence, we found that patients with lower expression of AGR2 have a better OS rate; hence, we suggest that AGR2 might play an oncogenic role in HCC progression and thus might be a useful prognostic marker for HCC progression.
Sorafenib decreases cell viability and increases cell apoptosis
First, to determine the effect of sorafenib on cell viability, the HCC cell lines were treated with various doses of sorafenib (5–10 µM, 24–48 h). Cell viability was significantly decreased with sorafenib treatment in a dose-dependent manner in J7, Hep3B, HepG2 and Huh7 cells according to the MTT assay (Fig. 2A-D). Moreover, flow cytometry was utilized to determine whether sorafenib influences HCC cell apoptosis. HepG2 and Huh7 cells were stimulated with 5 and 10 µM sorafenib for 24 h. Cell apoptosis was slightly induced with 5 µM sorafenib; however, the increase in the apoptosis rate compared with that in the control reached approximately 25% after 10 µM sorafenib stimulation in both HepG2 and Huh7 cells (Fig. 2E-H). Based on the results, we found that sorafenib can modulate HCC cell viability and apoptosis abilities; subsequently, we evaluated whether AGR2 is involved in the sorafenib-regulated phenotypes.
Sorafenib induces AGR2 secretion instead of transcriptional regulation
After we retrieved the dataset (GSE94550), AGR2 was induced in sorafenib-resistant Huh7 cells compared to parental cells (Fig. 1A). To further examine whether sorafenib could regulate AGR2 expression in HCC, RT–PCR and Western blotting were utilized in parental HCC cell lines treated with sorafenib. First, we found that the RNA level of AGR2 was not altered by sorafenib in J7, HepG2, Huh7 and Hep3B cells (Fig. 2I-L). However, the protein level of AGR2 in the cell lysate was unexpectedly decreased after sorafenib stimulation in these cells in a dose-dependent manner, as shown by Western blotting (Fig. 2M-P). Based on these contradictory results, sorafenib regulated the RNA and protein levels of AGR2 in parental HCC cells. Moreover, several reports have shown that AGR2 is an ER-resident protein that is also localized in the extracellular matrix, blood and urine[21, 22]. Therefore, we collected CM with 5 and 10 µM sorafenib treatments for 24 h. As expected, AGR2 was detected and induced by sorafenib in CM from Hep3B cells (Fig. 2Q, R). Based on the evidence, we suggest that sorafenib regulates AGR2 through posttranslational modification, not transcriptional regulation, in parental HCC cells.
Agr2 Plays A Role In Cell Viability And Apoptosis
To analyze the roles of AGR2 in HCC progression, we established AGR2-silenced Hep3B, HepG2 and Huh7 cells (Fig. 3A(a), B(a), C(a)), and cell viability was determined using the MTT assay. Cell viability was significantly decreased after AGR2 silencing, and the effect was more obvious when combined with sorafenib treatment (Fig. 3A(b), B(b), C(b)). The evidence suggests that AGR2 plays a role in promoting cancer progression. Additionally, we investigated whether AGR2 has effects on cell apoptosis. Flow cytometry analysis was utilized to demonstrate that sorafenib can induce significant cell apoptosis of approximately 10% compared with the control, which was more conspicuous after AGR2 silencing in the presence of sorafenib compared to the control (siNC) in Hep3B, HepG2 and Huh7 cells (Fig. 3D, E, F). Quantitative results shown in panels 3D(b), E(b), F(b). However, we found that AGR2 can be secreted into CM and detected by Western blotting (CM, Fig. 2Q, R). Previously, Fessart et al. reported that extracellular AGR2 is defined as an extracellular matrix pro-oncogenic regulator and makes cancer cells more aggressive [21]. Hence, we evaluated whether the addition of recombinant AGR2 protein regulates cell viability and apoptosis. Flow cytometry analysis indicated that 10 µM sorafenib can induce cell apoptosis up to 30–40% at 24 h in both HepG2 and Huh7 cells; however, the phenomenon can be reversed with 60 ng/ml recombinant AGR2 (rAGR2) (Fig. 3G, H). Based on the evidence, we suggest that AGR2 secreted into CM with sorafenib stimulation plays an oncogenic role in HCC cancer progression.
Sorafenib Induces Er Stress In Hcc
We have determined the functions and regulatory mechanisms of AGR2 by sorafenib. Based on a literature search, AGR2 has been demonstrated to be upregulated upon ER stress, and ER stress-related molecules such as protein kinase R (PKR)-like endoplasmic reticulum kinase (PERK), inositol-requiring enzyme 1 (IRE1) and activating transcription factor 6 (ATF6) are dysregulated in many types of cancers [16]. Therefore, the relationship between sorafenib and ER-related factors was determined. Among these molecules, we found that phospho-IRE1α (p-IRE1α) was upregulated with 10 µM sorafenib treatment in HepG2 and Huh7 cells using Western blotting (Fig. 4A). Moreover, X-box binding protein 1 (XBP1) has been reported as a unique transcription factor modulating gene expression for ERAD and promoting protein folding[33]. In the UPR, IRE1α is activated via oligomerization and autophosphorylation, followed by the activation of its endoribonuclease to cleave and splice XBP1. The activated IRE1α endoribonuclease can remove 26 nucleotides from the intron of XBP1, leading XBP1 from preform XBP1 (XBP1u: unspliced) into activated XBP1 (XBP1 s: spliced) [34]. Therefore, RT–PCR was used to detect the status of the IRE-1 downstream factor XBP1. We found that XBP1 was spliced from inactive XBP1u to active XBP1s after stimulation with sorafenib for 24 and 48 h in HepG2 and Huh7 cells by RT–PCR (Fig. 4B, C). Subsequently, we wanted to verify whether AGR2 could play a role in the process from inactive XBP1u to active XBP1 s. As expected, the sorafenib-induced XBP1s levels were more robust in HepG2 and Huh7 cells (Fig. 4D). In contrast, we further verified whether added recombinant AGR2 (rAGR2) could modulate the splicing of inactive XBP1u to activate XBP1s. Similarly, the levels of XBP1s induced by 10 µM sorafenib were reduced after stimulation with 60 ng/ml rAGR2 in HepG2 and Huh7 cells (Fig. 4E). Collectively, the evidence indicated that sorafenib induced HCC ER stress via the IRE1α-XBP1 cascade through AGR2 regulation.
AGR2 possesses diverse roles
To determine whether AGR2 plays a diverse role in resistant sublines compared to sorafenib-sensitive HCC cells, HepG2 sorafenib-resistant (HepG2-SR) and Huh7 sorafenib-resistant (Huh7-SR) cells were established (Fig. 5A, B). We applied 7 µM sorafenib to the culture medium to generate sorafenib-resistant cell lines for the following experiments. We found that cell viability was decreased by approximately 50% after sorafenib treatment (7 µM) in parental HepG2 and Huh7 cells (indicated as PC) using the MTT assay; however, sorafenib only reduced cell viability 10 ~ 20% in resistant cells (indicated as SR) (Fig. 5A, B). The evidence indicated that these SR cell lines are protected against sorafenib challenge. We found that AGR2 was related to resistance and upregulated by sorafenib in the GSE94550 dataset (Fig. 1A, B). Therefore, we wanted to verify whether the regulation could be observed in our sorafenib-resistant cells. Expectedly, AGR2 was highly expressed intracellularly in sorafenib-resistant HepG2 and Huh7 cells compared with parental cells (Fig. 5C). We also found that sorafenib reduced intracellular AGR2 levels in HepG2-SR and Huh7-SR cells (Fig. 5D), and the tendency was similar to that of sorafenib-regulated AGR2 in HepG2-PCs and Huh7-PCs (Fig. 2M-P). Moreover, we also found that sorafenib can induce AGR2 secretion in CM in HepG2-SR and Huh7-SR cells compared with controls. The effect was stronger in resistant cells than in parental cells (Fig. 5E). Collectively, the evidence indicated that AGR2 induction in both the cell lysate and CM of sorafenib-resistant cells was more robust than that in parental cells in the presence and absence of sorafenib. This evidence might indicate why SR cells are protected against sorafenib.
Sorafenib-resistant cells modulate ER stress and reduce cell apoptosis
To determine whether sorafenib-resistant cells can resist the effect of sorafenib-induced cell apoptosis, flow cytometry analysis was utilized to demonstrate that cell apoptosis was induced after stimulation with various doses of sorafenib for 24 h in HepG2 (Fig. 6A, B) and Huh7 (Fig. 6C, D) parental cells. The ratio of cells undergoing sorafenib-induced apoptosis reached approximately 30% among parental cells; however, the effect was not observed in sorafenib-resistant cells (Fig. 6A-D). Based on the evidence, we suggest that these resistant cells are protected against sorafenib toxicity, which prolongs cancer cell viability. Furthermore, we wanted to analyze whether AGR2 plays a vital role in inducing HCC resistance to sorafenib. Therefore, we silenced AGR2 in HepG2 SR and Huh7 SR cells and performed an apoptosis assay. The evidence indicated that sorafenib can induce cell apoptosis, but the effect was more robust after AGR2 knockdown (Fig. 6E-H). According to these findings, we speculated that AGR2 plays a critical role in inducing sorafenib resistance in HCC.
We found that the ER stress-related molecule p-IRE1α was induced with sorafenib treatment; subsequently, we determined whether p-IRE1α regulation could occur in sorafenib-resistant cells. The levels of p-IRE1α were decreased after sorafenib stimulation in sorafenib-resistant (SR) HepG2 and Huh7 cells, and the effect was inconsistent in parental cells (PC) by Western blot (Fig. 7A, B). Moreover, we found that the sorafenib-induced changes in the levels of spliced XBP1 (XBP1s), a downstream factor of IRE1α, were abolished or attenuated in HepG2-SR and Huh7-SR cells according to RT–PCR, and this was not observed in parental HepG2 and Huh7 cells (Fig. 7C).
Furthermore, we evaluated whether the effect could be mediated by AGR2 in sorafenib-resistant cells. The data indicated that the expression of XBP1s was induced after silencing AGR2 expression and was more elevated after stimulation with sorafenib (Fig. 7D). The evidence implied that AGR2 is more essential under cell stress conditions. Collectively, the contradictory evidence in sorafenib-sensitive and -resistant cells may imply that cells can resist sorafenib toxicity by modulating ER stress-related molecules to decrease cellular ER stress.
Based on the evidence, we proposed that sorafenib reduces cell viability and induces cell apoptosis via downregulation of AGR2 in the cell lysate and increased secretion in CM, which induces ER stress via upregulation of p-IRE1α and spliced XBP1 in HCC. However, the phenomenon of sorafenib-induced apoptosis was abolished in sorafenib-resistant cells, and this effect may occur through increased AGR2 expression in cell lysates and downregulation of p-IRE1α and spliced XBP1. Overall, AGR2 might modulate ER stress to protect cells from sorafenib toxicity and extend cell survival (Fig. 8).