Thioredoxin reductase 1 (TXNRD1) which is a selenocysteine-containing protein is overexpressed in many malignancies. Its role in the hepatocellular carcinoma (HCC) prognosis has not been investigated. In this study, we investigated whether TXNRD1 functions as an independent prognostic factor for HCC patients. We found TXNRD1 was overexpressed in HCC tissues and cells, immunohistochemical analysis suggested TXNRD1 was elevated in 57 of 120 (47.5%) clinical samples, and its level was increased with the increasing clinical stage. In addition, TXNRD1 expression was positively correlated with clinical stage (p=3.5e-5), N classification (p=4.4e-4), and M classification (p=0.037) of HCC patients. Kaplan-Meier analysis revealed that patients with high TXNRD1 expression had significantly shorter survival time than patients with low TXNRD1 expression. Multivariate analysis found TXNRD1 was an independent prognostic factor for HCC patients. In conclusion, our data suggested that TXNRD1 was a biomarker for the prognosis of patients with HCC.
State Key Projects on Infection Diseases of China2012ZX10002017-0052012ZX10002016-0232012ZX10002010-001-007National Natural Science Foundation of China813722438137057581301331Guangdong Natural Science Foundation2014A0303131312015A0303137562015A0303120132016A030313195Science and Technology Planning Project of Guangdong Province2014B0202280032014B0303010412015B020226004Science and Technology Planning Project of Guangzhou201400000001-31581000762013J4100081Guangzhou Municipal Health Bureau20141A01111720151A0111121. Introduction
HCC is one of the most common human tumors worldwide. According to the data published by the World Cancer Report 2014, the global incidence and mortality of HCC are 6% and 9%, respectively. About 746,000 patients died at 2012. Despite the advances in diagnosis and therapy, HCC still has poor prognosis and high recurrence rate, and looking for new diagnostic factors and therapeutic factors is critical for HCC therapy [1].
Selenoprotein TXNRD1 is a critical antioxidant enzyme catalyzing the NADPH-dependent reduction of thioredoxin to regulate cellular redox homeostasis [2]. It associated with cell proliferation, apoptosis, and transformation [3]. TXNRD1 is also upregulated in many human malignancies and functions as a prognostic factor for many tumors, such as oral squamous cell carcinomas [4], lung cancer [5, 6], breast cancer [7, 8], and astrocytomas [9]. Previous reports have shown that the hepatocarcinogenesis rate in TXNRD1 knockout mice was significantly increased, induced by carcinogen diethylnitrosamine (DEN), suggesting it could protect against chemically induced hepatocarcinogenesis through controlling the balance of cellular redox state, but it could not induce hepatocarcinogenesis [10]. Anticancer natural product gambogic acid interacts with TXNRD1 to inhibit its Trx-reduction activity and to increase the levels of reactive oxygen species to induce apoptosis in human HCC cell SMMC-7721; meanwhile, knockdown of TXNRD1 attenuates the gambogic acid cytotoxicity [11]. These findings suggested TXNRD1 might regulate the progression of HCC, but whether TXNRD1 could serve as a novel prognostic factor for HCC patients has not been studied until now. In this study, we determined the role of TXNRD1 in the prognosis for HCC patients, and we showed that TXNRD1 was upregulated in HCC cells and tissues and was an unfavorable prognostic factor for HCC patients.
2. Materials and Methods2.1. Cell Culture and Clinical Samples
Hepatocyte cell line LO2 cells were grown in RPMI-1640 medium (Gibco) supplemented with 10% fetal bovine serum (FBS, Gibco), 0.3% glutamine, 0.1 mM nonessential amino acids, 1 mM sodium pyruvate, and 10 μM β-mercaptoethanol (Gibco). HCC cell lines including Hep3B, BEL-7404, HepG2, Huh7, QGY-7703, MHCC97H, MHCC97L, BEL-7402, and HCCC-9801 were maintained in DMEM medium (Hyclone) supplemented with 10% FBS (Gibco). These cells were obtained from ATCC.
A cohort of 120 paraffin-embedded HCC samples was obtained from the Third Affiliated Hospital of Sun Yat-sen University (Guangzhou, China) with informed consent according to the Institute Research Ethics Committee; these patients were histologically and clinically diagnosed between July 2008 and November 2014. The detailed information was shown in Table 1, clinicopathological characteristics such as gender, age, clinical stage, and TNM classification were included.
Clinicopathological characteristics of HCC patients’ samples.
Number of cases (%)
Gender
Male
88 (73.3%)
Female
32 (26.7%)
Age (years)
≥45
76 (63.3%)
<45
44 (36.7%)
Clinical stage
I
30 (25.0%)
II
30 (25.0%)
III
34 (28.3%)
IV
26 (21.7%)
T classification
T1
39 (32.5%)
T2
33 (27.5%)
T3
18 (15.0%)
T4
30 (25.0%)
N classification
N0
98 (81.7%)
N1
22 (18.3%)
M classification
Yes
7 (5.8%)
No
113 (94.2%)
2.2. Western Blot
Cell lysates were extracted from cells using RIPA buffer (50 mM Tris (pH 7.4), 150 mM NaCl, 1% Triton X-100, 1% sodium deoxycholate, and 0.1% SDS supplemented with protease inhibitor cocktail (Roche)). Protein concentration was measured using BCA Protein Assay Kit (Pierce). Equivalent protein lysates were separated on 12% SDS-PAGE gels and transferred to PVDG membranes (Millipore) followed with blocking with 5% nonfat milk for 1 h at room temperature. After blocking, membranes were incubated with anti-TXNRD1 antibodies (1 : 1000, #15140, Cell Signaling) for 2 h at room temperature; then, the membranes were washed using TBST buffer and incubated with an anti-rabbit IgG, HRP-linked antibody (1 : 10000, #7074, Cell Signaling) for 1 h at room temperature. The band was detected using Amersham ECL Prime (GE Healthcare). The membranes were stripped and reprobed with an anti-β-actin antibody as the loading control.
2.3. Immunohistochemistry (IHC)
IHC was performed according to the standard method described previously [12, 13]. Two observers who were blinded to the clinical outcome evaluated and scored the degree of immunostaining of TXNRD1 independently. The staining index was calculated as the sum of the intensity of staining and the percentage of positively stained tumor cells. Scores for the intensity of staining were shown as follows: 0 (no staining), 1 (weak staining, light yellow), 2 (moderate staining, yellow brown), and 3 (strong staining, brown). Scores for the percentage of positively stained tumor cells were shown as follows: 0 (no positive tumor cells), 1 (<10% of the cells), 2 (10–50% of the cells), 3 (>50% of the cells). The staining index scores were 0, 1, 2, 3, 4, 6 and 9. The cut-off values for TXNRD1 expression were chosen based on the measure of heterogeneity using the log-rank test with respect to the overall survival. When SI score was ≥4, the tumor was considered to have high expression; otherwise, the tumor was considered to have low expression.
2.4. Statistical Analysis
All statistical analyses in this study were carried out using SPSS 13.0 (IBM). Chi-square test and Fisher’s test were used to determine the relationship between TXNRD1 expression and the clinicopathological characteristics, and Spearman correlation test was performed to calculate bivariate correlations between the study variables. The survival curves were plotted by the Kaplan-Meier method and compared by log-rank test. Univariate and multivariate analysis was carried out using Cox’s proportional hazards regression models. For all tests, a two-sided p value of less than 0.05 was considered to be statistically significant.
3. Results3.1. TXNRD1 Is Overexpressed in HCC Tissues and Cells
To investigate the prognostic value of TXNRD1, we first examined TXNRD1 expression in HCC tissues. We found TXNRD1 mRNA level was significantly upregulated in HCC tissues compared to liver nontumor tissues (p<0.0001), mRNA expression profiles of HCC and liver nontumor tissues came from GSE14520 (Figure 1(a)), we also used The Cancer Genome Atlas (TCGA) dataset to analyze TXNRD1 expression, and TXNRD1 mRNA level was also significantly upregulated in HCC tissues compared to normal live tissues (p<0.0001) (Figure 1(b)). Kaplan-Meier survival curves suggested patients with high TXNRD1 have poor outcome (p=0.016) (Figure 1(c)); this data also came from TCGA. Moreover, we determined TXNRD1 expression in HCC cells, and western blot assay suggested TXNRD1 protein level was upregulated in HCC cells compared to normal live cell LO2 (Figure 2(a)).
TXNRD1 was overexpressed in HCC tissues and patients with high TXNRD1 level had poor outcome. (a) The mRNA levels of TXNRD1 in HCC tissues and normal liver tissues, the data downloaded from GSE14520. (b) The mRNA levels of TXNRD1 in HCC tissues and normal liver tissues, the data downloaded from TCGA dataset. (c) Kaplan-Meier analysis for overall survival of HCC patients according to the level of TXNRD1. Data represent mean ± SD.
TXNRD1 expression was upregulated in HCC cells and was positively correlated with clinical stage. (a) Western blot analysis of TXNRD1 level in HCC cells and normal liver cell LO2. β-Actin was used as the loading control. (b) Representative photographs of TXNRD1 level in normal liver tissues, HCC samples in clinical stages II, IIIA, IIIB, and IV determined by IHC analysis. Original magnification: ×200.
We further examined TXNRD1 protein level in 120 paraffin-embedded archived HCC tissues using IHC. IHC analysis suggested 113 (94.2%) samples expressed TXNRD1 positively, and only 7 (5.8%) samples expressed TXNRD1 negatively. 63 samples (52.5%) had low TXNRD1 expression, 57 (47.5%) had high TXNRD1 expression (Table 2), and TXNRD1 was mainly located in the cytoplasm. These samples included 30 stage I samples, 30 stage II samples, 34 stage III samples, and 26 stage IV samples. IHC assay also showed TXNRD1 level was positively correlated with advancing clinical stage (Figure 2(b)); the detailed data was shown as follows: 23% (7/30) for stage I, 37% (11/30) for stage II, 53% (18/34) for stage III, and 81% (21/27) for stage IV (Table 3). This result suggested high TXNRD1 expression was associated with advancing clinical stage. These results suggested TXNRD1 was overexpressed in HCC tissues and might be a poor prognostic factor.
The expression of TXNRD1 in HCC.
Expression of TXNRD1
Negative
7 (5.8%)
Positive
113 (94.2%)
Low expression
63 (52.5%)
High expression
57 (47.5%)
Correlation between TXNRD1 expression and clinicopathological characteristics of HCC.
Characteristics
TXNRD1
Chi-square test p value
Fisher’s exact test p value
Low number of or no cases
High number of cases
Gender
Male
44
44
0.363
0.412
Female
19
13
Age (years)
≥45
40
36
0.970
1.000
<45
23
21
Clinical stage
I
23
7
3.5e-5
3.5e-5
II
19
11
III
16
18
IV
5
21
T classification
T1
24
15
0.265
0.263
T2
19
14
T3
8
10
T4
12
18
N classification
N0
59
39
4.4e-4
0.001
N1
4
18
M classification
No
62
51
0.037
0.052
Yes
1
6
3.2. Relationship between Clinicopathologic Features and TXNRD1 Expression in HCC
We investigated the relationship between TXNRD1 level and clinicopathological characteristics of HCC patients and found that TXNRD1 level was significantly correlated with clinical stage (p=3.5e-5) and N classification (p=4.4e-4 using Chi-square test, p=0.001 using Fisher’s exact test). TXNRD1 also was significantly correlated with M classification analyzed by Chi-square test (p=0.037) but was not significantly correlated with M classification analyzed by Fisher’s exact test (p=0.052). However, TXNRD1 level does not have significant correlation with gender, age, or T classification (Table 3). We used Spearman correlation analysis to confirm this result; Spearman correlation analysis showed TXNRD1 level was significantly correlated with clinical stage (r=0.405, p=5.3e-8), N classification (r=0.326, p=7.8e-7), and M classification (r=0.190, p=0.037), but TXNRD1 level was not significantly correlated with T classification (r=0.178, p=0.052) (Table 4). Taken together, high TXNRD1 level was positively correlated with clinical stage, N classification, and M classification.
Spearman correlation analysis between TXNRD1 and clinical pathologic factors.
Variables
TXNRD1 expression level
Spearman correlation
p value
Clinical stage
0.405
5.3e-8
T classification
0.178
0.052
N classification
0.326
7.8e-7
M classification
0.190
0.037
3.3. Patients with High TXNRD1 Level Have Poor Outcome
Using Kaplan-Meier analysis and the log-rank test, we observed that patients with low TXNRD1 levels had longer overall survival time (p<0.001, Figure 3(a)). Subgroup analyses found that patients with high TXNRD1 levels had poor prognosis in the clinical classification I-II (p<0.001) or in the clinical classification III-IV (p<0.05), suggesting TXNRD1 also was an unfavorable prognostic factor for clinical classifications I-II and III-IV (Figures 3(a) and 3(b)). Univariate Cox-regression analysis showed both clinical stage and TXNRD1 level were the significant poor prognostic factors (p=5.5e-9 and p=8.5e-8, resp.), and multivariate Cox-regression analysis found clinical stage and TXNRD1 level were independent prognostic factor for HCC patients (p=6.1e-7 and 5.4e-11, resp.) (Table 5). These findings revealed that TXNRD1 was an independent prognostic factor for HCC patients; patients with high TXNRD1 level had shorter survival time than those who had low TXNRD1 level.
Univariate and multivariate analyses of various prognostic parameters in patients with HCC Cox-regression analysis.
Univariate analysis
Multivariate analysis
Number of patients
p
Regression coefficient (SE)
p
Relative risk
95% confidence interval
Clinical stage
I
30
5.5e-9
1.807 (0.111)
6.1e-7
0.396
1.183–1.867
II
30
III
34
IV
26
Expression of TXNRD1
Low expression
63
8.5e-8
3.970 (0.230)
5.4e-11
1.071
1.792–4.755
High expression
57
Kaplan-Meier analysis to investigate the prognostic value of TXNRD1 level in overall survival and different clinical stage. (a) All patients. (b) Patients with clinical classification I-II. (c) Patients with clinical classification III-IV.
4. Discussion
In the present study, we investigated whether TXNRD1 could function as a prognostic factor for HCC patients. We first analyzed TXNRD1 expression using public database and found that TXNRD1 was significantly overexpressed in HCC tissues, and patients with high TXNRD1 level had poor outcome. Then, we further determined TXNRD1 expression in HCC cells and clinical specimens; TXNRD1 was overexpressed in HCC cells, and almost all collected specimens were TXNRD1-positive. Statistical analyses revealed that high TXNRD1 expression was positively correlated with clinical stage, N classification, and M classification. The patients who had low TXNRD1 level had longer survival time. Multivariate analysis showed TXNRD1 was an independent prognostic biomarker for HCC patients.
Thioredoxin system comprised by thioredoxin reductase (TrxR/TXNDR), thioredoxin (Trx), and NADPH plays an important role in maintaining intracellular redox homeostasis. There are three TXNDRs: TXNRD1 is located in cytosol, TXNRD2 is located in mitochondria, and TXNRD3 is located in testis [14]. IHC analysis also found TXNRD1 was mainly presented in cytosol, and the role and molecular mechanisms of TXNRD1 in HCC progression will be explored further.
TXNRD1 is involved in protecting from ROS by itself and via its function together with TXN to serve as an electron donor for peroxiredoxins or ribonucleotide reductase involved in DNA replication and repair [14]. TXNRD1 has been demonstrated to be overexpressed in many cancers, and ROS level also increased in many tumors, including HCC [15]. TXNRD1 promotes tumor growth, DNA replication, and tumorigenicity [16, 17], and knockdown of TXNRD1 also increases sensitivity of cancer cells to some chemotherapy drugs [18, 19]; these finding suggest it is a potential target for anticancer agents [20]. In HCC, ROS causes DNA damage and regulates p53, AP1, NF-AT, Nrf2/Maf, and multidrug resistance proteins to promote hepatocarcinogenesis and drug resistance [21]; we conferred that TXNRD1 might induce ROS generation to promote HCC development and drug resistance.
We found TXNRD1 could function as an unfavorable prognostic factor for HCC, while this finding is still to be replicated and to be verified in other patients’ population. In summary, we found TXNRD1 was upregulated in HCC tissues and cells, and TXNRD1 level was increased with increasing clinical stage. High TXNRD1 level was positively correlated with clinical stage, N classification, and M classification, and patients with high TXNRD1 level had poor outcome; TXNRD1 was a novel biomarker for the prognosis for HCC patients.
Ethical Approval
All procedures performed in the studies involving human participants were in accordance with the ethical standards of the institution or practice at which the studies were conducted.
Consent
Informed consent was obtained from all the individual participants included in the study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the State Key Projects on Infection Diseases of China (2012ZX10002017-005, 2012ZX10002016-023, and 2012ZX10002010-001-007), the National Natural Science Foundation of China (81372243, 81370575, and 81301331), Guangdong Natural Science Foundation (2014A030313131, 2015A030313756, 2015A030312013, and 2016A030313195), Science and Technology Planning Project of Guangdong Province (2014B020228003, 2014B030301041, and 2015B020226004), Science and Technology Planning Project of Guangzhou (201400000001-3, 158100076, and 2013J4100081), and Guangzhou Municipal Health Bureau (20141A011117, 20151A011112). These fundings were used to pay the salary of researchers and perform the experiments.
TorreL. A.BrayF.SiegelR. L.FerlayJ.Lortet-TieulentJ.JemalA.Global cancer statistics, 201220156528710810.3322/caac.212622-s2.0-84924271853NguyenP.AwwadR. T.SmartD. D. K.SpitzD. R.GiusD.Thioredoxin reductase as a novel molecular target for cancer therapy2006236216417410.1016/j.canlet.2005.04.0282-s2.0-33746430147ArnérE. S. J.HolmgrenA.Physiological functions of thioredoxin and thioredoxin reductase2000267206102610910.1046/j.1432-1327.2000.01701.x2-s2.0-0033775891IwasawaS.YamanoY.TakiguchiY.TanzawaH.TatsumiK.UzawaK.Upregulation of thioredoxin reductase 1 in human oral squamous cell carcinoma201125363764410.3892/or.2010.11312-s2.0-79551594518FernandesA. P.CapitanioA.SeleniusM.BrodinO.RundlöfA.-K.BjörnstedtM.Expression profiles of thioredoxin family proteins in human lung cancer tissue: correlation with proliferation and differentiation200955331332010.1111/j.1365-2559.2009.03381.x2-s2.0-69549116697SoiniY.KahlosK.NäpänkangasU.Kaarteenaho-WiikR.SäilyM.KoistinenP.PääkköP.HolmgrenA.KinnulaV. L.Widespread expression of thioredoxin and thioredoxin reductase in non-small cell lung carcinoma200176175017572-s2.0-0034901958CañasA.López-SánchezL. M.Valverde-EstepaA.HernándezV.FuentesE.Muñoz-CastañedaJ. R.López-PedreraC.De La Haba-RodríguezJ. R.ArandaE.Rodríguez-ArizaA.Maintenance of S-nitrosothiol homeostasis plays an important role in growth suppression of estrogen receptor-positive breast tumors2012146, article R15310.1186/bcr33662-s2.0-84870342895CadenasC.FranckensteinD.SchmidtM.GehrmannM.HermesM.GeppertB.SchormannW.MaccouxL. J.SchugM.SchumannA.WilhelmC.FreisE.IckstadtK.RahnenführerJ.BaumbachJ. I.SickmannA.HengstlerJ. G.Role of thioredoxin reductase 1 and thioredoxin interacting protein in prognosis of breast cancer2010123, article R4410.1186/bcr25992-s2.0-78049394036EsenH.ErdiF.KayaB.FeyziogluB.KeskinF.DemirL. S.Tissue thioredoxin reductase-1 expression in astrocytomas of different grades2015121345145810.1007/s11060-014-1661-52-s2.0-84925491384CarlsonB. A.YooM.-H.TobeR.MuellerC.Naranjo-SuarezS.HoffmannV. J.GladyshevV. N.HatfieldD. L.Thioredoxin reductase 1 protects against chemically induced hepatocarcinogenesis via control of cellular redox homeostasis20123391806181310.1093/carcin/bgs2302-s2.0-84865545696DuanD.ZhangB.YaoJ.LiuY.SunJ.GeC.PengS.FangJ.Gambogic acid induces apoptosis in hepatocellular carcinoma SMMC-7721 cells by targeting cytosolic thioredoxin reductase201469152510.1016/j.freeradbiomed.2013.12.0272-s2.0-84893502890ChenL.LiY.LinC. H.ChanT. H. M.ChowR. K. K.SongY.LiuM.YuanY.-F.FuL.KongK. L.QiL.LiY.ZhangN.TongA. H. Y.KwongD. L.-W.ManK.LoC. M.LokS.TenenD. G.GuanX.-Y.Recoding RNA editing of AZIN1 predisposes to hepatocellular carcinoma201319220921610.1038/nm.30432-s2.0-84873531277LuG.ShanT.HeS.RenM.ZhuM.HuY.LuX.ZhangD.Overexpression of CENP-H as a novel prognostic biomarker for human hepatocellular carcinoma progression and patient survival20133052238224410.3892/or.2013.26752-s2.0-84885066805HolmgrenA.LuJ.Thioredoxin and thioredoxin reductase: current research with special reference to human disease2010396112012410.1016/j.bbrc.2010.03.0832-s2.0-77952311188TrachoothamD.AlexandreJ.HuangP.Targeting cancer cells by ROS-mediated mechanisms: a radical therapeutic approach?20098757959110.1038/nrd28032-s2.0-67650071137YooM.-H.XuX.-M.CarlsonB. A.GladyshevV. N.HatfieldD. L.Thioredoxin reductase 1 deficiency reverses tumor phenotype and tumorigenicity of lung carcinoma cells200628119130051300810.1074/jbc.C6000122002-s2.0-33744958179YooM.-H.XuX.-M.CarlsonB. A.PattersonA. D.GladyshevV. N.HatfieldD. L.Targeting thioredoxin reductase 1 reduction in cancer cells inhibits self-sufficient growth and DNA replication2007210e111210.1371/journal.pone.00011122-s2.0-39749177838PoerschkeR. L.MoosP. J.Thioredoxin reductase 1 knockdown enhances selenazolidine cytotoxicity in human lung cancer cells via mitochondrial dysfunction201181221122110.1016/j.bcp.2010.09.0242-s2.0-78650173598TobeR.YooM.-H.FradejasN.CarlsonB. A.CalvoS.GladyshevV. N.HatfieldD. L.Thioredoxin reductase 1 deficiency enhances selenite toxicity in cancer cells via a thioredoxin-independent mechanism2012445342343010.1042/BJ201206182-s2.0-84864346767SmartD. K.OrtizK. L.MattsonD.BradburyC. M.BishtK. S.SieckL. K.BrechbielM. W.GiusD.Thioredoxin reductase as a potential molecular target for anticancer agents that induce oxidative stress200464186716672410.1158/0008-5472.CAN-03-39902-s2.0-4644313264KuoM. T.SavarajN.Roles of reactive oxygen species in hepatocarcinogenesis and drug resistance gene expression in liver cancers200645970170910.1002/mc.202402-s2.0-33748807185