Haptoglobin expression correlates with tumor differentiation and five-year overall survival rate in hepatocellular carcinoma

Abstract

Elevated serum haptoglobin (Hp) is identified as a prognostic marker in multiple types of solid tumors, which is correlated with poor prognosis. HCC is one of the major causes of cancer deaths in worldwide, which remains poor prognosis and is clinically urgent for discovering early diagnostic markers. However, except for serum Hp, the correlation of tumor Hp expression with hepatocellular carcinoma (HCC) progression is still unclear. In this study, we evaluated and identified the tissue Hp expression as a prognostic marker to predict the survival rate of HCC patients. To evaluate the prognostic value of Hp expression for HCC, two cohorts were enrolled in our study, including total 130 matched pair tissue sections (both adjacent non-tumorous and tumor tissue derived from same patient) of HCC patients from Changhua Christian Hospital (CCH) and total 316 RNA-seq data with clinical information of HCC patients from The Cancer Genome Atlas (TCGA) database. In contrast to other types of cancers, HCC tumor tissues have lower Hp protein expression in CCH cohort and have lower Hp mRNA expression in TCGA cohort as compared with adjacent non-tumorous tissues (p < 0.001). Moreover, lower Hp expression is significantly correlated with different stages of HCC cancer differentiation in CCH cohort (one-way ANOVA, p < 0.001). Most importantly, lower Hp expression is highly correlated with poor five-year overall survival rate in TCGA cohort (p < 0.01). Based on our data, we conclude that tissue Hp expression positively correlates with better HCC tumor differentiation and increased five-year overall survival rate of HCC patients. The results indicated that tissue Hp is potentially a prognostic marker for HCC patients. Our findings may further provide a new insight of effective treatments along with biopsy diagnosis of HCC patients.

About the Authors:

Chun-San Tai

Contributed equally to this work with: Chun-San Tai, Yan-Ren Lin

Affiliations Department of Biological Science and Technology, National Chiao Tung University, Hsinchu, Taiwan, Institute of Molecular Medicine and Bioengineering, National Chiao Tung University, Hsinchu, Taiwan, Program of Interdisciplinary Neuroscience, National Chiao Tung University, Hsinchu, Taiwan

Yan-Ren Lin

Contributed equally to this work with: Chun-San Tai, Yan-Ren Lin

Affiliations Department of Emergency Medicine, Changhua Christian Hospital, Changhua, Taiwan, School of Medicine, Kaohsiung Medical University, Kaohsiung, Taiwan, School of Medicine, Chung Shan Medical University, Taichung, Taiwan

Tsung-Han Teng

Affiliation: Institute of Molecular Medicine and Bioengineering, National Chiao Tung University, Hsinchu, Taiwan

Ping-Yi Lin

Affiliations Transplant Medicine & Surgery Research Centre, Changhua Christian Hospital, Changhua, Taiwan, Department of Medical Research, China Medical University Hospital, Taichung, Taiwan

Siang-Jyun Tu

Affiliation: Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu, Taiwan

Chih-Hung Chou

Affiliation: Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu, Taiwan

Ya-Rong Huang

Affiliation: Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu, Taiwan

Wei-Chih Huang

Affiliation: Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu, Taiwan

Shun-Long Weng

Affiliations Department of Medicine, Mackay Medical College, New Taipei City, Taiwan, Department of Obstetrics and Gynecology, Hsinchu MacKay Memorial Hospital, Hsinchu, Taiwan, Mackay Junior College of Medicine, Nursing and Management, Taipei, Taiwan

Hsien-Da Huang

* E-mail: [email protected] (WLC); [email protected] (YLC); [email protected] (HDH)

Affiliations Department of Biological Science and Technology, National Chiao Tung University, Hsinchu, Taiwan, Institute of Bioinformatics and Systems Biology, National Chiao Tung University, Hsinchu, Taiwan

Yao-Li Chen

* E-mail: [email protected] (WLC); [email protected] (YLC); [email protected] (HDH)

Affiliations School of Medicine, Kaohsiung Medical University, Kaohsiung, Taiwan, Department of General Surgery, Changhua Christian Hospital, Changhua, Taiwan

Wen Liang Chen

* E-mail: [email protected] (WLC); [email protected] (YLC); [email protected] (HDH)

Affiliation: Department of Biological Science and Technology, National Chiao Tung University, Hsinchu, Taiwan

ORCID http://orcid.org/0000-0003-4308-3128

Introduction

Haptoglobin (Hp) is a secreted acute phase protein that is mainly produced in liver[1]. Elevated levels of Hp in plasma and locally in tumor tissue is closely associated with different clinical conditions, especially in various types of cancers such as ovarian cancer, colorectal cancer, pancreatic cancer, breast, and lung cancer[2–6]. Previous study reported the serum Hp could be regarded as an independent prognostic factor for epithelial ovarian cancer patients since patients with higher serum levels of Hp had significantly poor survival rates[7]. Additionally, Sun et al. demonstrated that serum Hp is highly expressed in colorectal cancer and is correlated with poor survival rate[8]. Moreover, it has been shown that the increased serum levels of Hp in non-small cell lung cancer (NSCLC) is associated with advanced TNM stages and distant metastasis[9]. These findings suggest the potential role of Hp in cancer prognosis and biomarker, which is worth and clinically urgent for further investigations.

Liver cancer is the second leading cause of cancer mortality and the fifth incidence rate among men and the ninth incidence rate in women[10]. Among all types of liver cancer, hepatocellular carcinoma (HCC) is the most commonly diagnosed type, representing 83% of all cases[11]. Currently, the critical issue for HCC is the poor five-year overall survival rate of only 30%-40%, which HCC patients often suffered postoperative recurrence and poor response to systemic chemotherapeutic treatments[12]. Despite the improvement in therapeutic treatments such as, radiotherapy, drugs, surgical resections or liver transplantation, the high incidence of tumor recurrence and metastasis still lead to poor prognosis of HCC patients[13]. Thus, identification of promising diagnostic markers as a direction for effective clinical treatments is worth for development. We here attempt to discover a prognostic marker in tissue level, which can be utilized with the biopsy for HCC diagnosis and, at the same time, provide more information for further therapeutic strategies to improve the overall survival rate of HCC patients. Unfortunately, besides the current pathological interpretation of tissue sections, there is still lack of promising tissue markers in predicting the outcome of HCC patients. Since the clinical significance of Hp is reported by previous studies and it is also highly expressed in liver, we hypothesized that Hp might possessed with a critically prognostic role in HCC development. Moreover, the differences of tissue Hp levels between tumor and adjacent non-tumorous parts (in the same patient) that may reflect the tendency of Hp secretion (under different cancer-related stress), which had never been well analyzed for prognostic value. Therefore, in this study, we aim to firstly identify tissue Hp level difference (between adjacent non-tumorous tissues and tumor tissues) as a prognostic marker for HCC outcomes using two cohorts, including prospective patients sample collections (n = 130) from Changhua Christian Hospital (CCH) and Hp mRNA analysis (n = 316) from The Cancer Genome Atlas (TCGA) database.

In this study, we enrolled a total number of 130 HCC patients from Changhua Christian Hospital (CCH) and 316 HCC patients from The Cancer Genome Atlas (TCGA) database. By using the immunohistochemistry staining for measuring the Hp protein expression of CCH tissue sections and computational methods analyzing for the Hp mRNA expression from TCGA datasets, we verified the tissue Hp expression between adjacent non-tumorous tissues and tumor tissues. We further analyzed the correlation of Hp expression with HCC cancer cell differentiation and five-year overall survival rate. Based on these data, we concluded that Hp serves as a prognostic marker for HCC patients. Higher tissue Hp expression is correlated with well-differentiated HCC cancer cells, which finally reflects on improved five-year overall survival rate.

Materials and methods

Changhua Christian Hospital (CCH) HCC patients

A total of 130 HCC patients, confirmed by histopathologic examination, were prospectively collected from Changhua Christian Hospital (Changhua, Taiwan) between 2011 to 2013. The study was approved by the institutional review board of the Changhua Christian hospital and written informed consent was obtained from all patients (IRB ID: 130704). The clinicopathological characteristics of enrolled HCC patients (gender, age at primary diagnosis, pathological stage, tumor size, pathological TNM classification, HCC cancer cell differentiation, hepatitis B history, liver cirrhosis history and alpha-fetoprotein) were shown in Table 1.

[Figure omitted. See PDF.]

Table 1. Clinicopathological characteristics of CCH HCC patients.

https://doi.org/10.1371/journal.pone.0171269.t001

Preparation of tissue core array

Surgically resected matched pair tissue sections (both adjacent non-tumorous and tumor tissues were originated from the same patient) were collected from 130 CCH HCC patients and all of the pathological sections were reviewed and reported by senior pathologists. Briefly, embedded tissue cores were manually obtained from the tissue blocks by a 10-gauge syringe needle. The cores in demand were placed in a warm cast, which contains melted paraffin wax, and were arranged into an 8 × 12 matrix with an orientated marker from removing the core at the upper-left corner. Each slide contained 96 cores and each patient had two to five cores of their matched pair specimens, depending on the availability of tissue blocks. The positions of each core were recorded on reference sheets to facilitate data acquisition.

The Cancer Genome Atlas (TCGA) HCC patients and RNA-seq data analysis

RNA sequencing (RNA-seq) data and patients’ clinical and demographic information of HCC patients from TCGA Liver Hepatocellular Carcinoma (LIHC) dataset were downloaded from the TCGA data portal (https://tcga-data.nci.nih.gov/). The results presented here are based on the data generated by the TCGA Research Network: (http://cancergenome.nih.gov/). We selected Illumina HiSeq 2000 RNA Sequencing Version 2 platform for RNA-seq data and we selected 316 HCC tumor tissues and 40 adjacent non-tumorous tissues for further analysis, which are excluded the inconsistent features of HCC patients such as NA pathological stage and zero or NA information in days to death and days to last follow up. The normalized results for mRNA expression, the level 3 expression data from RNA-seq data, of a gene in TCGA Liver Hepatocellular Carcinoma (LIHC) dataset were analyzed by TCGA Research Network through using MapSplice to perform the alignment and RSEM to perform the quantitation. To identify the correlation of selected genes with Hp in either low tumor Hp or high tumor Hp HCC tumor tissues, we calculated the Pearson correlation coefficient and its confidence intervals and p-values for each gene by R, version 3.3.1[14].

Immunohistochemistry staining and assessment

Immunohistochemistry (IHC) staining for Hp was performed by biotin-streptavidin- peroxidase procedures using a Vectastain ABC Kit (Vector Laboratories, Burlingame, CA). Fresh tumor tissues and adjacent non-tumorous liver tissue from HCC patients in tissue core arrays were fixed in formalin and embedded in paraffin. Each sliced serial section is 4 μm and were mounted on glass slides, de-paraffinized with xylene and dehydrated with alcohol according to the manufacturer’s instructions. Antigen retrieval was conducted by heating the slice in a microwave for 2 minutes at 900W then incubated the slides with 0.3% H₂O₂ solution in methanol for 30 minutes to block endogenous peroxidase. After rinsing the slides with PBS for three times, the tissues were treated with 20 μg/mL of protease K in the buffer containing 50 mM Tris, 1mM EDTA, 0.5% Triton X-100, pH 8.0 for 15 minutes at room temperature to expose the immunoreactive parts of Hp. Subsequently, the rest of the immunostaining procedures were performed according to the above-mentioned description. Each stained core was scored by Dr. Tsung-Han Teng and Dr. Yan-Ren Lin on the percentage of positively stained (0–100%) and stain intensity (negative, 0; weak, +1; moderate, +2; and strong +3). Total 130 matched pair tissues had two cores available, for which the mean percentage and maximum stain intensity was averaged from two scores. Besides, according to WHO rules and previous references[15, 16], senior pathologists reviewed the stages of cancer cell differentiation of all the tissue sections.

Statistical analysis

All data analyses, statistical significances and figures, except for Fig 2B and 2C, were performed by R, version 3.3.1[14]. Fig 2B and 2C were plotted by GraphPad Prism version 6.0 (GraphPad Software, La Jolla, CA, USA). Statistical significances between two groups of each sample were analyzed with two-sided Wilcoxon rank-sum test. For multi-groups of samples, the statistical significance was analyzed with one-way ANOVA. The correlation analysis of selected genes with Hp is calculated by the Pearson correlation coefficient. Survival rate analysis was presented in Kaplan-Meier curves and statistical significance is assessed using the log-rank test. The optimal cut point for survival rate analysis is determined by analyzing p-values form continuous Hp expression of all the HCC patients, which is based on PrognoScan database[17]. For all statistical analysis, p-values of less than 0.05 were considered as statistically significant.

Results

Tissue Hp protein expression level is correlated with HCC cancer differentiation in CCH HCC patients

To compare Hp protein expression levels between HCC tumor tissues and adjacent non-tumorous tissues, we analyzed 130 matched pairs of HCC tissue sections, including the pathological stage (stage I, II, III, and IV) with the Hp protein expression level by immunohistochemistry (IHC) staining. The clinicopathological characteristics of 130 HCC patients from CCH cohort are shown in Table 1. Under 100X magnification, the positive Hp staining is presented in adjacent non-tumorous tissue, yet largely absent in tumor tissue (Fig 1A, 100X magnification). Intense Hp staining can be clearly observed in adjacent non-tumorous cells under 400X magnification. In contrast, the tumor portion was presented in negative Hp staining (Fig 1A, 400X magnification). Comparing between matched pair tumor tissues and adjacent non-tumorous tissues of 130 HCC patients, the Hp expression is significantly decreased in tumor tissues than adjacent non-tumorous tissues (p < 0.001, Fig 1B). These results indicated that in HCC patients, tissue Hp protein expression levels are reduced in tumor. To further investigate the prognostic value of Hp in HCC cancer progression, we analyzed the correlation between Hp expression levels and HCC cancer cell differentiation. As shown in Fig 1C, among different levels of cancer differentiation subgroups and adjacent non-tumorous tissues, Hp expression was significantly correlated with cell differentiation (one-way ANOVA, p < 0.001). Hp protein expression is significantly higher in adjacent non-tumorous cells than other three subgroups of HCC cancer cell differentiation (adjacent non-tumorous vs. well differentiation, p < 0.05; adjacent non-tumorous vs. moderate differentiation, p < 0.001; adjacent non-tumorous vs. poor differentiation, p < 0.001) (Fig 1C). It is worth noting that well-differentiated cancer cells show significantly higher Hp expression than poor-differentiated cancer cells (p < 0.05) (Fig 1C). In the schematic diagram of Fig 1D, we demonstrated that HCC cancer cells with better differentiation phenotype show higher Hp expression (one-way ANOVA, p < 0.001). The gradient increase of Hp expression is from poor-differentiated cancer cells to well-differentiated cancer cells to the outer adjacent non-tumorous tissue (Fig 1D). By zooming into the junction of tumor part and adjacent non-tumorous part, clear boundary was observed and IHC staining of Hp expression confirmed the results (Fig 1E).

[Figure omitted. See PDF.]

Fig 1. Hp protein expression and cancer differentiation analysis of CCH HCC patients.

(A) An example of Hp IHC staining in adjacent non-tumorous tissues and tumor tissues under 100X and 400X magnification. (B) Hp protein expression analysis between matched pair adjacent non-tumorous tissues and tumor tissues of HCC patients from CCH (p < 0.001). (C) Boxplot of Hp protein expression in adjacent non-tumorous tissues and three stages HCC cancer differentiation (well, moderate, and poor) (p < 0.001). (D) The schematic diagram interpreted the correlation between Hp expression and HCC cancer differentiation (p < 0.001). (E) An example of IHC staining for measuring the Hp protein expression at the junction area between adjacent non-tumorous tissue and tumor tissue.

https://doi.org/10.1371/journal.pone.0171269.g001

[Figure omitted. See PDF.]

Fig 2. Hp gene expression and five-year overall survival rate analysis in TCGA HCC patients.

(A) The boxplot demonstrated Hp gene expression level between adjacent non-tumorous tissues and tumor tissues in TCGA HCC patients (p < 0.001). (B) The Hp expression plot showed a total of 316 HCC patients that are ordered by Hp mRNA expression values (RPKM) from the lowest to the highest. The X-axis represents the accumulative number of HCC patients and the Y-axis represents Hp expression value in RPKM. The dashed line shows optimal cut point which dichotomized HCC patients into high (red) and low (green) Hp expression groups. (C) The p-value plot displayed the survival significance of each potential cut point of Hp expression measurement. According to each cut point, HCC patients are dichotomized into high and low Hp expression groups, and are calculated by log-rank test for survival significances. The X-axis represents the accumulative number of HCC patients on the order same as the Hp expression plot in panel B and the Y-axis represents raw p-values on a log10 scale. The optimal cut point is colored in red and is indicated by dashed line. (D) The Kaplan-Meier plot of HCC survival analysis. HCC patients were dichotomized into low Hp expression group (n = 34) and high Hp expression group (n = 282). The Kaplan-Meier survival curve showed high tumor Hp expression group has increased five-year overall survival rate, comparing with low tumor Hp expression group in TCGA HCC patients (p < 0.001).

https://doi.org/10.1371/journal.pone.0171269.g002

Hp expression is correlated with five-year overall survival rate of TCGA HCC patients

To further confirm our clinical results, we next analyzed RNA-seq data of HCC patients in TCGA database. Total 316 HCC patients from TCGA database and the detailed clinical pathological characteristics are displayed in Table 2. As expected, the results showed the similar trend with CCH cohort, which Hp expression is decreased in tumor tissues and also highly expressed in adjacent non-tumorous tissues in TCGA HCC patients (p < 0.001, Fig 2A). To evaluate the prognostic value of Hp in HCC patients, we next performed overall survival analysis to all (n = 316) HCC patients from TCGA library. TCGA HCC patients were dichotomized into two groups, low tumor Hp expression (n = 34) and high tumor Hp expression (n = 282), through the optimal cut point, which is determined by calculating p-values of continuous Hp expression of TCGA HCC patients (Fig 2B and 2C). As shown in Fig 2B and 2C, HCC patients were ordered by Hp expression from the lowest to the highest (Fig 2B) and were calculated by log-rank test for survival significances (Fig 2C). Surprisingly, the result showed that HCC patients with higher tumor Hp expression (n = 282) were correlated with higher overall survival rate, when comparing to patients with lower Hp expression (n = 34) (log-rank test, p < 0.001) (Fig 2D).

[Figure omitted. See PDF.]

Table 2. Clinicopathological characteristics of TCGA HCC patients.

https://doi.org/10.1371/journal.pone.0171269.t002

Hp expression is negatively correlated with poor cancer differentiation markers

To verify our results that lower Hp expression is correlated with poor HCC cancer cell differentiation, we selected several markers (TWIST1, LAMB1, THY1, EZH2, SALL4, and TCF3) that are closely associated with the activation of poor HCC cancer differentiation and analyzed the Pearson correlation coefficient of each selected marker to Hp. Based on the RNA-seq data and previous results of survival analysis, we again divided all stages of HCC patients into two groups, low tumor Hp expression (n = 34) and high tumor Hp expression (n = 282). We then calculated the Pearson correlation coefficient of each marker to Hp in low tumor Hp expression and high tumor Hp expression respectively (Table 3). In low tumor Hp expression group, TWIST1 (r = -0.54, p < 0.01), LAMB1 (r = -0.58, p < 0.001), and THY1 (r = -0.48, p < 0.01) showed negative correlation with Hp. In high tumor Hp expression group, EZH2 (r = -0.53, p < 0.001), SALL4 (r = -0.4, p < 0.001), TCF3 (r = -0.45, p < 0.001) were also negatively correlated with Hp.

[Figure omitted. See PDF.]

Table 3. The Pearson correlation coefficient of selected cancer differentiation markers with Hp.

https://doi.org/10.1371/journal.pone.0171269.t003

Discussion

The clinical significance of Hp is associated with the progression of diverse malignancies such as autoimmune disorders, diabetes, cardiovascular diseases, and, especially, cancers[18]. Elevated serum Hp expression is often observed in various types of cancer, which is identified as a poor prognostic marker[8, 19, 20]. Previous study showed that serum Hp is highly expressed in HCC patients as compared with liver cirrhosis (LC) patients, suggesting that Hp could be another potential biomarker complementary to alpha-fetoprotein[21]. However, the clinical significance of tumor Hp expression in HCC patients is not yet been revealed. In this study, we first reported the alteration of Hp expression in tumor tissues as compared with adjacent non-tumorous tissues in HCC patients, which can be regarded as a novel prognostic marker for HCC patients.

To the beginning, we examined tissue Hp expression in the two enrolled cohorts, CCH and TCGA HCC patients. As shown in Fig 1A and 1B, we found significantly higher tissue Hp expression in adjacent non-tumorous tissues than in tumor tissues by IHC staining of tissues sections from CCH HCC patients. In addition, as shown in Fig 2A, we validated our clinical observation by analyzing RNA-seq data in TCGA database, which showed the same trend as the tissue sections from CCH HCC patients. Hp expression is often raised in various malignancies. However, our study showed opposite results. We proposed that the reasonable explanation of such result is because of the detected location of Hp expression. Previous studies mostly focused on the serum level of Hp expression, which is in contrast to our observation on tissue levels of Hp expression. Since Hp is a secreted acute phase protein that can also be slightly secreted by other organs[1], except for liver, under inflammation responses, serum Hp expression may gathered from other places instead of originating from the cancer tissues. Furthermore, chronic inflammation in tumor microenvironment reduced the expression level of pro-inflammatory cytokines thereby affecting Hp expression in HCC tumor tissue, since the expressional alteration of Hp expression in hepatocyte in regulated by pro-inflammatory cytokines including IL-6, IL-1β, and TNF-α[22]. In our study, we firstly proved the lower Hp expression in tumor tissues, comparing with adjacent non-tumorous tissues, by directly look into tissue levels of Hp expression through protein and mRNA level.

Given this intriguing results of much higher Hp expression in adjacent non-tumorous tissues than tumor tissues, we next endeavored to assess the potential role of Hp in HCC prognosis. Surprisingly, statistically significant correlation of Hp expression with three stages of cancer differentiation and adjacent non-tumorous cells was discovered (Fig 1C). The gradient increasing Hp expression from poor-differentiated cancer cells to well-differentiated cancer cells to peripheral adjacent non-tumorous cells is illustrated in schematic diagram (Fig 1D). We found that higher tissue Hp expression is correlated with better cancer cell differentiation phenotype. Additionally, the observation of elevated Hp expression in adjacent non-tumorous tissues could be elaborated by the biological stimulation of Hp under inflammation, infections, tissue damages and even stress[23, 24]. Therefore, we proposed that intense Hp expression in adjacent non-tumorous tissues is stimulated by extremely extruding stress from growing tumor (Fig 1D). Furthermore, it might be the reason of elevated serum levels of Hp expression is often detected in diverse cancers in previous reports. IHC staining further showed the intense Hp expression in adjacent non-tumorous tissues and a clear boundary between adjacent non-tumorous tissues and tumor tissues representing well-differentiated cancer cells (Fig 1E). Cancer cell differentiation is closely linked with patients’ prognosis and overall survival rate, since poor differentiated cancer cells are more aggressive than well-differentiated cancer cells, resulting in cancer metastasis and tumor recurrence. As we discovered the statistically significance of higher Hp expression in well-differentiated cancer cells, we next proceeded to used tumor Hp expression as the factor in analyzing five-year overall survival rate to assess whether Hp could serve as a prognostic marker in HCC patients. Through analyzing the p-value of each HCC patients, which is ordered by Hp expression, the result provided the optimal cut point for survival analysis and displayed the potential of Hp as a prognostic value in HCC patients (Fig 2B and 2C). In agreement with previous results, Kaplan-Meier overall survival analysis of HCC patients proved that HCC patients with high tumor Hp expression result in improved five-year overall survival rate. By contrast, HCC patients with low tumor Hp expression showed a shorter life expectancy (Fig 2D). Based on these evidences, we revealed that higher Hp expression in tumor tissues is not only correlated with better cancer cell maturation, but also correlated with increasing five-year overall survival rate. We thus demonstrated that the higher tissue Hp expression is a potentially prognostic marker with the better outcome for HCC patients.

Lastly, not merely conclude by the clinical observation and survival analysis, we further attempted to explore the groups of genes, which are related to promoting poor cancer cell differentiation, is correlated with Hp. The six selected genes (TWIST1, LAMB1, THY1, EZH2, SALL4, and TCF3) are involved in cell differentiation, HCC cancer stem cells (CSCs) markers, and CSCs involved pathways. The negative Pearson correlation coefficient demonstrated negative correlations of these selected genes with Hp in either low tumor Hp expression group (n = 34) or high tumor Hp expression group (n = 282) (Table 3). TWIST1 is one of the key transcription factor in activating epithelial-mesenchymal transition (EMT), which is a crucial process for HCC cell de-differentiation[25]. According to previous study, Hp expression is associated with EMT in buccal cancer patients[26]. Another selected marker for cell differentiation is LAMB1, which is in the family of extracellular matrix glycoprotein and is implicated with cell differentiation. Up-regulation of LAMB1 showed poor differentiation and aggressive cancer phenotypes such as migration and invasion[27, 28]. THY1(CD90) is one of the representative HCC Cancer Stem Cells (CSCs) markers. CSCs represents undifferentiated cells that bear stem-like properties and differentiation ability, which contributes to poor prognosis and recurrence[29]. Thus, the negative correlation of TWIST1 (r = -0.54, p < 0.01), LAMB1 (r = -0.58, p < 0.001), and THY1 (r = -0.48, p < 0.01) in low tumor Hp expression groups could be a signature of poor cancer differentiation also the prognostic role of Hp in HCC cancer cell differentiation. On the other hand, negative correlation of EZH2 (r = -0.53, p < 0.001), SALL4 (r = -0.4, p < 0.001), TCF3 (r = -0.45, p < 0.001) with Hp in high tumor Hp expression group could also be an indication of relative well cancer differentiation. EZH2 and SALL4 are two critical genes in regulating CSC properties. Previous study indicated that EZH2 closely regulates differentiation of hepatic stem/progenitor cells[30, 31]. SALL4 is also reported in increasing stem cell properties and causing poor prognosis[31, 32]. Last but not least, TCF3 is one of the components in Wnt/β-catenin signaling pathway, which plays important role in activating liver cancer stem cells[33]. Taken together, the correlation of these selected genes (TWIST1, LAMB1, THY1, EZH2, SALL4, and TCF3) not only as a signature of poor cancer differentiation in low tumor Hp expression group and relative well cancer differentiation in high tumor Hp expression group, but also demonstrates the correlation of Hp expression and cancer differentiation.

To conclude, our data suggested that Hp expression could serve as a novel prognostic marker for HCC patients, which is correlated with cancer cell differentiation and five-year overall survival rate. Our findings contribute in several clinical significances, suggesting Hp can potentially predict HCC prognosis by cancer cell differentiation and then provide information for adequate clinical managements.

Acknowledgments

The authors thank to Ms. Li Shu Chao for technical assistance and contribution for our study.

Author Contributions

1. Conceptualization: CST YRL HDH YLC WLC.

2. Data curation: CST SJT CHC YRH WCH SLW HDH WLC.

3. Formal analysis: CST YRL THT PYL SJT CHC YRH WCH SLW HDH YLC WLC.

4. Funding acquisition: YRL SLW.

5. Investigation: CST YRL THT PYL SJT CHC YRH WCH SLW HDH YLC WLC.

6. Methodology: CST CHC WLC HDH.

7. Project administration: THT CHC SLW YLC HDH WLC.

8. Resources: THT CHC SLW YLC HDH WLC.

9. Software: CST SJT CHC YRH WCH SLW HDH WLC.

10. Supervision: CHC YLC HDH WLC.

11. Validation: CST YRL THT PYL SJT CHC YRH WCH SLW HDH YLC WLC.

12. Visualization: CST YRL THT PYL SJT CHC YRH WCH SLW HDH YLC WLC.

13. Writing – original draft: CST YRL PYL CHC HDH YLC WLC.

14. Writing – review & editing: CST YRL THT PYL SJT CHC YRH WCH SLW HDH YLC WLC.

Citation: Tai C-S, Lin Y-R, Teng T-H, Lin P-Y, Tu S-J, Chou C-H, et al. (2017) Haptoglobin expression correlates with tumor differentiation and five-year overall survival rate in hepatocellular carcinoma. PLoS ONE 12(2): e0171269. https://doi.org/10.1371/journal.pone.0171269

References

1. Galicia G, Ceuppens JL. Haptoglobin Function and Regulation in Autoimmune Diseases: INTECH Open Access Publisher; 2011.

2. Okuyama N, Ide Y, Nakano M, Nakagawa T, Yamanaka K, Moriwaki K, et al. Fucosylated haptoglobin is a novel marker for pancreatic cancer: a detailed analysis of the oligosaccharide structure and a possible mechanism for fucosylation. Int J Cancer. 2006;118(11):2803–8. pmid:16385567

3. Garibay-Cerdenares OL, Hernandez-Ramirez VI, Osorio-Trujillo JC, Gallardo-Rincon D, Talamas-Rohana P. Haptoglobin and CCR2 receptor expression in ovarian cancer cells that were exposed to ascitic fluid: exploring a new role of haptoglobin in the tumoral microenvironment. Cell Adh Migr. 2015;9(5):394–405. PubMed Central PMCID: PMCPMC4955374. pmid:26211665

4. Hoagland LFt, Campa MJ, Gottlin EB, Herndon JE 2nd, Patz EF Jr. Haptoglobin and posttranslational glycan-modified derivatives as serum biomarkers for the diagnosis of nonsmall cell lung cancer. Cancer. 2007;110(10):2260–8. pmid:17918261

5. Takeda Y, Shinzaki S, Okudo K, Moriwaki K, Murata K, Miyoshi E. Fucosylated haptoglobin is a novel type of cancer biomarker linked to the prognosis after an operation in colorectal cancer. Cancer. 2012;118(12):3036–43. pmid:22006099

6. Carlsson MC, Cederfur C, Schaar V, Balog CI, Lepur A, Touret F, et al. Galectin-1-binding glycoforms of haptoglobin with altered intracellular trafficking, and increase in metastatic breast cancer patients. PLoS One. 2011;6(10):e26560. PubMed Central PMCID: PMCPMC3196588. pmid:22028908

7. Zhao C, Annamalai L, Guo C, Kothandaraman N, Koh SC, Zhang H, et al. Circulating haptoglobin is an independent prognostic factor in the sera of patients with epithelial ovarian cancer. Neoplasia. 2007;9(1):1–7. PubMed PMID: 17325738; PubMed Central PMCID: PMCPMC1803031. pmid:17325738

8. Sun L, Hu S, Yu L, Guo C, Sun L, Yang Z, et al. Serum haptoglobin as a novel molecular biomarker predicting colorectal cancer hepatic metastasis. Int J Cancer. 2016;138(11):2724–31. pmid:26756179

9. Lu J, Wang Y, Yan M, Feng P, Yuan L, Cai Y, et al. High serum haptoglobin level is associated with tumor progression and predicts poor prognosis in non-small cell lung cancer. Oncotarget. 2016.

10. Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86. pmid:25220842

11. Farazi PA, DePinho RA. Hepatocellular carcinoma pathogenesis: from genes to environment. Nat Rev Cancer. 2006;6(9):674–87. pmid:16929323

12. Dong P, Ma L, Liu L, Zhao G, Zhang S, Dong L, et al. CD86(+)/CD206(+), Diametrically Polarized Tumor-Associated Macrophages, Predict Hepatocellular Carcinoma Patient Prognosis. Int J Mol Sci. 2016;17(3):320. PubMed Central PMCID: PMCPMC4813183. pmid:26938527

13. Han K, Kim JH, Ko GY, Gwon DI, Sung KB. Treatment of hepatocellular carcinoma with portal venous tumor thrombosis: A comprehensive review. World J Gastroenterol. 2016;22(1):407–16. PubMed Central PMCID: PMCPMC4698503. pmid:26755886

14. Team RC. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2015. URL http://www.R-project.org. 2016.

15. Brunt EM. Histopathologic features of hepatocellular carcinoma. Clinical Liver Disease. 2012;1(6):194–9.

16. Schlageter M, Terracciano LM, D'Angelo S, Sorrentino P. Histopathology of hepatocellular carcinoma. World J Gastroenterol. 2014;20(43):15955–64. PubMed Central PMCID: PMCPMC4239483. pmid:25473149

17. Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genomics. 2009;2:18. PubMed Central PMCID: PMCPMC2689870. pmid:19393097

18. Quaye IK. Haptoglobin, inflammation and disease. Transactions of the royal society of tropical medicine and hygiene. 2008;102(8):735–42. pmid:18486167

19. Tabassum U, Reddy O, Mukherjee G. Elevated serum haptoglobin is associated with clinical outcome in triple-negative breast cancer patients. Asian Pac J Cancer Prev. 2012;13(9):4541–4. PubMed PMID: 23167376. pmid:23167376

20. Haas M, Kern C, Kruger S, Michl M, Modest DP, Giessen C, et al. Assessing novel prognostic serum biomarkers in advanced pancreatic cancer: the role of CYFRA 21–1, serum amyloid A, haptoglobin, and 25-OH vitamin D3. Tumour Biol. 2015;36(4):2631–40. pmid:25472579

21. Shu H, Kang X, Guo K, Li S, Li M, Sun L, et al. Diagnostic value of serum haptoglobin protein as hepatocellular carcinoma candidate marker complementary to alpha fetoprotein. Oncol Rep. 2010;24(5):1271–6. PubMed PMID: 20878120. pmid:20878120

22. Bicho MC, da Silva AP, Medeiros R, Bicho M. The role of haptoglobin and its genetic polymorphism in cancer: a review. Acute Phase Proteins. 2013.

23. Gruys E, Toussaint MJ, Niewold TA, Koopmans SJ. Acute phase reaction and acute phase proteins. J Zhejiang Univ Sci B. 2005;6(11):1045–56. PubMed Central PMCID: PMCPMC1390650. pmid:16252337

24. Dobryszycka W. Biological functions of haptoglobin—new pieces to an old puzzle. Eur J Clin Chem Clin Biochem. 1997;35(9):647–54. PubMed PMID: 9352226. pmid:9352226

25. Giannelli G, Koudelkova P, Dituri F, Mikulits W. Role of epithelial to mesenchymal transition in hepatocellular carcinoma. J Hepatol. 2016.

26. Lee CC, Ho HC, Chien SH, Hsiao SH, Hung SK, Huang TT, et al. Association of acute phase protein-haptoglobin, and epithelial-mesenchymal transition in buccal cancer: a preliminary report. Clin Chem Lab Med. 2013;51(2):429–37. pmid:23093274

27. Ho DW, Yang ZF, Yi K, Lam CT, Ng MN, Yu WC, et al. Gene expression profiling of liver cancer stem cells by RNA-sequencing. PLoS One. 2012;7(5):e37159. PubMed Central PMCID: PMCPMC3351419. pmid:22606345

28. Yu G-R, Kim S-H, Park S-H, Cui X-D, Xu D-Y, Yu H-C, et al. Identification of molecular markers for the oncogenic differentiation of hepatocellular carcinoma. Experimental and Molecular Medicine. 2007;39(5):641. pmid:18059140

29. Yamashita T, Wang XW. Cancer stem cells in the development of liver cancer. The Journal of clinical investigation. 2013;123(5):1911–8. pmid:23635789

30. Aoki R, Chiba T, Miyagi S, Negishi M, Konuma T, Taniguchi H, et al. The polycomb group gene product Ezh2 regulates proliferation and differentiation of murine hepatic stem/progenitor cells. Journal of hepatology. 2010;52(6):854–63. pmid:20395008

31. Chiba T, Iwama A, Yokosuka O. Cancer stem cells in hepatocellular carcinoma: Therapeutic implications based on stem cell biology. Hepatol Res. 2016;46(1):50–7. pmid:26123821

32. Oikawa T, Kamiya A, Zeniya M, Chikada H, Hyuck AD, Yamazaki Y, et al. Sal-like protein 4 (SALL4), a stem cell biomarker in liver cancers. Hepatology. 2013;57(4):1469–83. pmid:23175232

33. Wang R, Sun Q, Wang P, Liu M, Xiong S, Luo J, et al. Notch and Wnt/beta-catenin signaling pathway play important roles in activating liver cancer stem cells. Oncotarget. 2016;7(5):5754–68. PubMed Central PMCID: PMCPMC4868719. pmid:26735577