Back to Journals » Cancer Management and Research » Volume 10
Genetic variants in the exon region of versican predict survival of patients with resected early-stage hepatitis B virus-associated hepatocellular carcinoma
Authors Liu X, Han C, Liao X , Yu L, Zhu G , Su H, Qin W, Lu S, Ye X, Peng T
Received 8 January 2018
Accepted for publication 8 March 2018
Published 4 May 2018 Volume 2018:10 Pages 1027—1036
DOI https://doi.org/10.2147/CMAR.S161906
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 3
Editor who approved publication: Professor Harikrishna Nakshatri
Xiaoguang Liu,* Chuangye Han,* Xiwen Liao, Long Yu, Guangzhi Zhu, Hao Su, Wei Qin, Sicong Lu, Xinping Ye, Tao Peng
Department of Hepatobiliary Surgery, The First Affiliated Hospital of Guangxi Medical University, Nanning, Guangxi Zhuang Autonomous Region, China
*These authors contributed equally to this work
Background: The upregulated expression of versican (VCAN) promotes the proliferation, invasion, and metastasis of various types of human cancer cells, including hepatocellular carcinoma (HCC) cells.
Patients and methods: In this study, genetic variants in the exon region of VCAN were genotyped by DNA sequencing. Prognostic values of VCAN exon single nucleotide polymorphisms (SNPs) were assessed by Kaplan–Meier with the log-rank test, and uni- and multivariate Cox proportional hazard regression model.
Results: A total of 111 patients with resected hepatitis B virus-associated early-stage HCC were collected for genotyping VCAN exon SNPs using Sanger DNA sequencing. Haplotype analysis was performed using Haploview 4.2. Survival data were analyzed using Kaplan–Meier curves and Cox proportional hazards regression analyses. The rs2652098, rs309559, rs188703, rs160278, and rs160277 SNPs were significantly associated with overall patient survival (p<0.001, p=0.012, p=0.010, p=0.007, and p=0.007, respectively). Patients carrying the TAGTG haplotype had a poorer prognosis than those with the most common CGAAT haplotype, after adjusting for tumor size, tumor capsule, and regional invasion (adjusted hazard ratio [HR] =2.06, 95% CI: 1.27–3.34, p=0.003). Meanwhile, patients with the TAGTG haplotype and a larger tumor size or an incomplete tumor capsule had an increased risk of death, compared with the others (adjusted HR =3.00, 95% CI: 1.67–5.36, p<0.001; and adjusted HR = 1.99, 95% CI = 1.12–3.55, p = 0.02, respectively). The online database mining analysis showed that upregulated VCAN expression in HCC tissues was associated with a poor overall survival of 148 HCC patients.
Conclusion: Genetic variants in the exon region of VCAN were associated with overall survival in patients with resected early-stage hepatitis B virus-associated HCC, and may be a potential prognostic biomarker.
Keywords: versican, single nucleotide polymorphism, hepatitis B virus, hepatocellular carcinoma, survival
Introduction
Liver cancer is a lethal disease, and in 2012 alone, there were an estimated 782,500 new liver cancer cases and 745,500 liver cancer-related deaths occurred in the world.1 To date, 70%–90% of primary liver cancers are hepatocellular carcinoma (HCC). The highest incidence of HCC is in Asia, and half of HCC cases occur in China.2–4 The treatment options for HCC patients include hepatic resection, orthotopic liver transplantation, loco-regional chemotherapy, radiation therapy, and systemic or targeting therapy.5 Among them, hepatic resection is the first-line curative treatment, but it is only available for early-stage HCC patients; and the vast majority of HCC patients are not suitable for surgery due to hepatitis B virus (HBV)- or hepatitis C virus (HCV)-induced liver cirrhosis or poor liver function.6 Even after a hepatectomy, there is a 50%–60% chance of tumor recurrence, leading to a poor outcome for patients with HBV-associated HCC.6 Overall, a large tumor size, multiple tumor lesions, poor tumor cell differentiation, vascular invasion, and degree of inflammation have been associated with a poor HCC prognosis.7 Thus, identification and evaluation of novel biomarkers could help medical oncologists to select patients for more aggressive or targeted therapy in order to improve the outcome of HCC after surgery.
Versican (VCAN) is a large extracellular matrix proteoglycan that is expressed in majority human tissues. VCAN is localized on human chromosome 5q14.3, covering a genomic sequence of 109.4 kb.8 VCAN cDNA contains 15 exons and codes a protein with a molecular mass of >1000 kDa, which belongs to the aggrecan/versican proteoglycan family. VCAN protein functions to regulate cell adhesion, proliferation, migration, and angiogenesis; and it plays a central role in tissue morphogenesis and maintenance.9–14 Altered VCAN expression is closely related to proliferation, invasion, and metastasis in various types of human cancer cells, including HCC, and upregulated VCAN expression has been associated with a poor prognosis of cervical cancer patients.15 To date, VCAN has 4 splice variants, that is, V0, V1, V2, and V3, each of which appears to have distinct biological functions.16–18 Previous studies also have shown that VCAN genetic polymorphisms are associated with the development of gastric cancer and that altered VCAN expression is associated with outcomes in ovarian, breast, prostate, and gastric cancer.19–24 Versican V1 has been shown to promote HCC metastasis, and its expression has been associated with the prognosis of HCC patients after curative resection.9 VCAN contains 2 globular domains, that is, the G1 domain (N-terminal; G1-D), possessing the hyaluronan-binding property, and the G3 domain (C-terminal; G3-D), consisting of epidermal growth factor repeats.25 Another study has reported that the presence of the rs188703 variant allele A and the rs160277 variant allele T is a susceptibility genotype for gastric cancer, while the rs309559 variant allele G and the rs160278 variant allele A are also susceptibility factors in the development of gastric cancer.24 Therefore, in this study, we assessed 5 different polymorphisms in the exon region of VCAN and then associated them with the prognosis of resected early-stage HBV-associated HCC patients. These genetic polymorphisms in the VCAN exon region lead to changes in the VCAN protein amino acids; thus, we also performed an online database mining analysis to associate VCAN expression in HCC tissues with HCC patient survival. This study will provide insightful information regarding these VCAN single-nucleotide polymorphisms (SNPs) as a novel prognostic biomarker for HCC patients.
Patients and methods
Study population
The study protocol was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University (approval no. 2015KY-E-032). The written informed consents were obtained from all the participants enrolled in the study. From January 2005 to September 2013, we collected 111 patients with early-stage HBV-associated HCC among 1150 patients diagnosed with HBV-related HCC from The First Affiliated Hospital, Guangxi Medical University (Nanning, China). All of these 111 patients underwent a hepatectomy and had a confirmed HCC diagnosis based on histopathology. Tumor differentiation was assessed according to the Edmondson–Steiner grading system.26 Grade I was considered well-differentiated HCC, grades II and III were considered moderately differentiated HCC, and grade IV was considered poorly differentiated HCC. The inclusion criteria for patient enrollment into this study were as follows: 1) patients who tested positive for hepatitis B surface antigen, 2) patients with Barcelona Clinic Liver Cancer (BCLC) stage A,27,28 and 3) patients with Child–Pugh class A liver disease. The exclusion criteria were as follows: 1) patients with HCV infection and 2) patients with other cancers. The clinicopathological data were collected from the patients’ medical records and are summarized in Table 1.
In this study, radical hepatic resection was defined as the complete removal of the tumor lesion with an incisional margin ≥2 cm and without a visible tumor thrombus in the hepatic vein or in the extrahepatic lymph nodes. These patients were followed up in our outpatient clinic or through a telephone interview every 3 months, and the last follow-up was conducted on September 1, 2014. This study was approved by the Ethics Committee of The First Affiliated Hospital of Guangxi Medical University (Nanning, China).
SNP selection, DNA extraction, and genotyping
The SNPs in the VCAN exon region that are responsible for changes in amino acids were selected with a minor allele frequency >5% in the Han Chinese population from the 1000 Genomes Project (March 2012, build GRCh37/hg19). Five SNPs in the VCAN gene exon region (rs2652098, rs309559, rs188703, rs160278, and rs160277) were identified and assessed in this study.
Tissue samples from HCC lesions were obtained after surgical resection, immediately snap-frozen, and stored at −80°C until needed for DNA extraction. Genomic DNA from the HCC tissue samples was extracted using a TIANamp Genomic DNA Kit (Tiangen Biotech [Beijing] Co., Ltd., Beijing, China), according to the manufacturer’s protocol. VCAN genetic variations were genotyped by Sanger DNA sequencing on an ABI Prism 3100 (Applied Biosystems, Shanghai Sangon Biological Engineering Technology & Services Co., Ltd., Shanghai, China) after polymerase chain reaction amplification using different primers (Table 2). Ten samples were randomly selected for repeated polymerase chain reaction, the DNA sequences were used to validate the genotyping data, and the results showed 100% concordance.
Online mining of VCAN gene expression data
To associate and confirm VCAN expression with HCC patient survival, we searched the online database (GSE14520 [https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE14520] in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database) and obtained a cohort of 148 early-stage HBV-associated HCC patients (BCLC stage A) from Zhongshan Hospital, Fudan University (Shanghai, China).29,30 We used this cohort of patients and data as a validated set. These patients were similar to our cohort, with the same inclusion and exclusion criteria.
Statistical analysis
Hardy–Weinberg equilibrium for each SNP was estimated using a goodness-of-fit χ2-test with one degree of freedom. The linkage disequilibrium between SNPs of the VCAN gene and the haplotype block were conducted using Gabriel’s algorithm in Haploview version 4.2.31 For all of these 5 SNPs, 1 reasonable genetic model (dominant, additive, or recessive models) was used. Overall survival (OS) was defined from the date of the hepatectomy to death or the last follow-up (September 1, 2014). The univariate analysis between clinical features and OS was compared using the log-rank test, and those with p < 0.1 were entered into the Cox proportional hazard regression model. A p-value <0.05 was defined as statistically significant. All statistical analyses were performed using SPSS 16.0 (SPSS, Inc., Chicago, IL, USA).
Results
Patient characteristics and clinical predictors
In this study, we collected 111 cases of HBV-related HCC with BCLC stage A and Child–Pugh class A liver disease. Their clinicopathological characteristics are summarized in Table 1. In brief, the median patient age was 45 years old, there were 87.4% male patients, 67.6% of patients had a maximum tumor diameter ≥10 cm, and 84.7% of patients had liver cirrhosis. Moreover, 52.3% of patients had a low serum alpha-fetoprotein (AFP) level (≤400 ng/mL), 31.5% of patients had ever smoked tobacco, and 42.3% of patients drank alcohol. Seventy patients (63.1%) underwent a radical hepatic resection, 72 (64.9%) patients received a hepatic resection and adjuvant transarterial chemoembolization (TACE), and 95 (85.6%) patients did not receive adjuvant antiviral treatment.
The median follow-up period was 47 months (range: 10–117 months). During the follow-up, 69 (62.2%) patients had died, including 50 (45.0%) patients who died due to HCC recurrence. Tumor size, tumor capsule, and regional invasion were all associated with the OS of these patients (tumor size, p=0.020; tumor capsule, p=0.024; tumor regional invasion, p<0.001 with the log-rank test). Patients with a maximum tumor diameter ≥5 cm had a 94.0% increased risk of death (hazard ratio [HR] = 1.94, 95% CI: 1.10–3.45). Patients with an incomplete tumor capsule also had a higher risk of death (HR =1.72, 95% CI: 1.07–2.78). Moreover, patients who were treated with postoperative adjuvant antiviral therapy or adjuvant TACE did not show a decreased risk of death (adjuvant antiviral treatment, p=0.29; adjuvant TACE, p=0.12 with the log-rank test).
Associations of SNPs with overall patient survival
The Hardy–Weinberg equilibrium data for these 5 SNPs are listed in Table 2. All these 5 SNPs had 3 genotypes, and all of them were significantly associated with overall patient survival after adjusting for tumor size, tumor capsule, and regional invasion according to the Cox proportional hazard regression analyses (Table 3). The rs2652098 SNP was shown to be the most significantly associated with a shorter OS in patients who carried the T allele (CT/TT) than in those with the wild-type CC allele (adjusted HR =2.42, 95% CI: 1.49–3.94, under the dominant genetic model, Table 3).
Haplotype analysis of the association between SNPs and OS of HCC patients
Haplotype analysis of the SNPs rs2652098, rs309559, rs188703, rs160278, and rs160277 identified 1 haplotype block (Figure 1A). Haplo4 TAGTG was associated with a shorter median survival time, compared with those with the most common genotype CGAAT (34.0 vs. 83.0 months, p = 0.01 with the log-rank test; Figure 1B). After adjusting for tumor size, tumor capsule, and regional invasion, patients with haplotype TAGTG had a significantly increased risk of death, compared with those with the genotype CGAAT (adjusted HR =2.06, 95% CI: 1.27–3.34, p=0.003; Table 4).
Stratified analysis of the association of VCAN haplotypes with OS of HCC patients
Next, we carried out further stratified analyses by age, gender, AFP level, tumor size, tumor capsule, radical hepatic resection, adjuvant TACE, and adjuvant antiviral treatment to evaluate associations between the haplotype block of VCAN and OS of patients with HBV-associated HCC (Figure 2). Haplo1, Haplo2, or Haplo3 had a significant protective effect in male patients (adjusted HR =0.64, 95% CI: 0.42–0.96, Figure 2) and in patients with serum AFP levels >400 ng/mL (adjusted HR =0.35, 95% CI: 0.20–0.63, Figure 2) and adjuvant TACE (adjusted HR =0.58, 95% CI: 0.36–0.94). Haplo1, Haplo2, or Haplo3 also significantly decreased the death risk in patients without radical hepatic resection (adjusted HR =0.41, 95% CI: 0.20–0.83, Figure 2) and adjuvant antiviral therapy (adjusted HR =0.54, 95% CI: 0.36–0.82, Figure 2).
Combined effects of the VCAN haplotype as well as tumor size and capsule on OS of HCC patients
The combined effects of the VCAN haplotype and clinical features (tumor size and capsule) on the prediction of OS of patients with early-stage HBV-associated HCC are shown in Table 5. In brief, compared with patients with Haplo1, Haplo2, or Haplo3 and a small tumor size (≤5 cm), those with Haplo4 and a large tumor size (>5 cm) had a significantly increased risk of death (adjusted HR =3.00, 95% CI: 1.67–5.36, p<0.001, Table 5; Figure 3) after adjusting for tumor size, tumor capsule, and regional invasion (excluding the stratified factor in each stratum). Similarly, patients with Haplo4 and an incomplete tumor capsule had a significantly increased risk of death, compared with patients with Haplo1, Haplo2, or Haplo3 and a complete tumor capsule (adjusted HR =1.99, 95% CI: 1.12–3.55, p=0.020, Table 5; Figure 3).
Association of VCAN expression with OS of patients with early-stage HBV-associated HCC
Since these VCAN SNPs could alter the amino acids in VCAN protein, we retrieved data on the Fudan cohort from GSE14520 in the NCBI GEO database. These 148 cases of early-stage HBV-associated HCC (BCLC stage A) were divided into 2 groups according to the VCAN expression levels in tumor tissues, that is, the high VCAN expression group (n=76), in which the VCAN expression levels were above the median value, and a low VCAN expression group (n=72), in which the VCAN expression levels were below the median value. Our data showed that a low VCAN expression was associated with a better overall patient survival than a high VCAN expression (p=0.013 with the log-rank test; Figure 4).
Discussion
In the current study, we genotyped 5 different SNPs in the VCAN exon region and associated them with the clinicopathological and survival data of patients with resected early-stage HBV-associated HCC (classified as BCLC A and Child-Pugh class A). Our data showed that rs2652098, rs309559, rs188703, rs160278, and rs160277 were significantly associated with the clinicopathological and survival data of these patients. These 5 SNPs had strong linkage disequilibrium; however, compared with the other 4 VCAN SNPs, rs2652098 was more representative as a potential biomarker for OS of these patients after adjusting for the clinicopathological features. Moreover, haplotype analysis of the association between these 5 SNPs and overall patient survival showed that patients with the TAGTG haplotype had the worst prognosis. In addition, our online data mining analysis showed that a low VCAN expression in HCC tissues was associated with a better OS of HCC patients. Clinicopathological features, like tumor size, tumor capsule, and regional invasion, were also associated with the OS of these patients. Further studies will confirm the usefulness of these 5 VCAN SNPs in predicting the prognosis of resected HCC.
In the current study, we first retrieved SNP data on VCAN exonic polymorphisms with amino acid changes from the NCBI dbSNP database (GRCh37.p13) and found rs2652098 (C>T, Ser300Leu), rs309559 (A>G, Lys1516Arg), rs188703 (G>A, Arg1826His), rs160278 (T>A, Phe2301Tyr), and rs160277 (G>T, Asp2937Tyr) for genotyping. To date, there are no reports on these VCAN SNPs being associated with HCC outcome. Our current study demonstrated that these VCAN exonic polymorphisms were associated with the clinicopathological and survival data of patients with resected early-stage HBV-associated HCC. Haplo4 TAGTG (constituted by rs2652098T, rs309559A, rs188703G, rs160278T, and rs160277G) was associated with a poorer prognosis than those with the most common genotype CGAAT. A previous study has revealed that rs188703G and rs160277G can increase the susceptibility to intestinal-type gastric cancer.24 rs188703, rs160278, and rs160277 are localized on exon 8 of the VCAN gene that codes the GAG-beta domain (G3 domain) of VCAN protein. The GAG-beta domain is associated with enhanced cell proliferation and activation of the epidermal growth factor receptor (EGFR) signaling pathway;9 thus, blockade of the GAG-beta domain activity can reduce cell proliferation and EGFR pathway activation as well as reduce susceptibility to cancer development.9 Thus, several studies have indicated that targeting VCAN protein expression or activity could be used as a potential tumor therapeutic strategy.10,32–34 For example, VCAN G3 has been shown to enhance cancer cell resistance to chemotherapeutic agents, but a selective EGFR inhibitor has been demonstrated to inverse its resistance.10,11,32,35 Furthermore, VCAN G3 is able to enhance tumor angiogenesis by formation of a complex together with vascular endothelial growth factor via directly binding to fibronectin.11
Our current study also showed that HCC patients with Haplo4 and a larger tumor size or an incomplete tumor capsule had a poorer OS than those with Haplo1, Haplo2, or Haplo3 and a smaller tumor size or a complete tumor capsule. A previous study using PolyPhen prediction (Polymorphism Phenotyping at http://genetics.bwh.harvard.edu/pph/) has shown that rs160277 can affect the function or structure of VCAN protein through transcriptional regulation.36 Meanwhile, both rs188703 and rs160277 can affect the function or structure of the GAG-beta domain to interfere with cell proliferation and activation of the EGFR pathway. Overall, the use of these VCAN genetic polymorphisms combined with clinical characteristics might be better at predicting the survival of these early-stage HBV-associated HCC patients after a hepatectomy. However, VCAN, a large chondroitin sulfate proteoglycan, is a major component of the extracellular matrix.37,38 Expression of VCAN protein is associated with cell growth and differentiation;39 in addition, it is associated with invasion and metastasis of different tumors.11–14,40–45 VCAN also plays an important role in transforming growth factor (TGF)-β-mediated tumor cell invasion and metastasis.46,47 TGF-β-induced VCAN expression is able to activate the nuclear factor-κB signaling pathway and upregulate the expression of CD44, matrix metalloproteinase 9, and the hyaluronan-mediated motility receptor in ovarian cancer.47 VCAN G3 has been shown to promote growth, migration, and metastasis of mouse mammary tumor cells through the upregulation of EGFR signaling.48 Indeed, detection of VCAN expression can predict tumor metastasis and progression.49,50 Our current data also showed that increased VCAN expression was associated with a poorer OS of HCC patients. Recently, Xia et al found that FoxQ1 promotes HCC metastasis by transactivation of VCAN expression.51 In addition, they also demonstrated that VCAN expression contributes to a poorer overall patient survival and higher recurrence rates of HCC.51 Similarly, human cervical cancer patients with positive VCAN expression had lower OS rates.15,52 Thus, further studies will confirm whether detection of VCAN SNPs or protein expression is more useful in predicting the survival of HCC patients.
Of note, our study does have potential limitations. First, as an SNP study, the sample size was small and a fraction of the clinical data was missing. Second, there was a lack of functional characterizations between these SNPs and VCAN expression or activity.
Conclusion
Our current study provided a proof-of-principle to associate VCAN rs2652098, rs309559, rs188703, rs160278, and rs160277 SNPs with a poor OS of patients with resected early-stage HBV-associated HCC, even after adjusting for other clinicopathological features, such as tumor size, tumor capsule, or regional invasion. Our current data also confirmed that a low VCAN expression was associated with a better OS of HCC patients than a high VCAN expression. However, future studies are needed to validate our current findings and clarify the underlying mechanisms.
Acknowledgments
This work was supported in part by the National Nature Science Foundation of China (81560535, 81072321, 30760243, 30460143, and 30560133), 2009 Program for New Century Excellent Talents in University (NCET), Guangxi Nature Sciences Foundation (GuiKeGong 1104003A-7), and Guangxi Health Ministry Medicine Grant (Key-Scientific Research-Grant Z201018). The present study is also partly supported by Scientific Research Fund of the Health and Family Planning Commission of Guangxi Zhuang Autonomous Region (Z2016318). We would also acknowledge the support provided by the National Key Clinical Specialty Programs (General Surgery and Oncology) and the Key Laboratory of Early Prevention and Treatment for Regional High-Incidence-Tumor (Guangxi Medical University), Ministry of Education, China. The authors thank Minhao Peng, Lequn Li, Xiao Qin, Kaiyin Xiao, Xigang Chen, Bin Chen, Zhixiong Su, Ming Su, Zhang Wen, Jingning Lu, Ning Peng, Hai Zhu for providing part of HCC samples for this study, who are from the Department of Hepatobiliary Surgery, the First Affiliated Hospital of Guangxi Medical University. Thanks also go to Xue Qin (Department of Clinical Laboratory, The First Affiliated Hospital of Guangxi Medical University) for providing patient’s serum AFP assay results. We would like to acknowledge the researcher Jiaquan Li and Ying Gui from Guangxi Medical University for their contribution to specimen management. The authors thank the contributors of GSE14520 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE14520) for sharing the HCC data on open access. The authors also want to thank the Mingyi Li, Nianping Chen, Ming Chen, Xiaoyu Tan, Huilai Miao, Dong Dai, Ai Mo, Wei Dai, Hao Xu, Yongguang Yang, Manzhou Lin, Guohua Liu, who are from the Departments of Hepatobiliary Surgery, Affiliated Hospital of Guangdong Medical University, Zhanjiang, Guangdong 524001, China, for their contribution to manuscript revision. In addition, we would also like to acknowledge the helpful comments on this paper received from our reviewers.
Disclosure
The authors report no conflicts of interest in this work.
References
© 2018 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.