Genetic variants in the exon region of versican predict survival of patients with resected early-stage hepatitis B virus-associated hepatocellular carcinoma.

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.

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.

Table 1

Clinical features of patients with resected early-stage HBV-associated HCC

Variables	n=111	MST (months)	HR (95% CI)	p-value
Age (years)				0.846
≤45	63	47	1
>45	48	45	1.05 (0.65–1.68)
Gender				0.61
Male	97	45	1
Female	14	51	0.83 (0.40–1.73)
Smoking status				0.49
None	76	57	1
Ever	35	42	1.19 (0.72–1.97)
Drinking status				0.44
None	64	42	1
Ever	47	82	0.83 (0.51–1.34)
AFP (ng/mL)a				0.64
≤400	58	57	1
>400	47	41	1.12 (0.69–1.83)
Tumor size (cm)				0.020
≤5	36	80	1
>5	75	42	1.94 (1.10–3.45)
No. of tumors				0.92
Single (n=1)	106	45	1
Multiple (n>1)	5	65	0.94 (0.30–3.00)
Tumor capsule				0.024
Complete	59	88	1
Incomplete	52	38	1.72 (1.07–2.78)
Regional invasion				<0.001
Absent	101	74	1
Present	10	11	3.68 (1.90–7.13)
Tumor differentiationb				0.64
Well	9	79	1
Moderately	74	42	1.48 (0.59–3.72)
Poorly	2	40	0.96 (0.11–8.25)
Cirrhosis				0.63
Absent	17	82	1
Present	94	45	1.18 (0.60–2.31)
Radical liver resectionc				0.69
Yes	70	45	1
No	38	38	1.10 (0.67–1.81)
Adjuvant antiviral treatment				0.29
No	95	65	1
Yes	16	25	1.41 (0.74–2.69)
Adjuvant TACE				0.12
No	39	101	1
Yes	72	42	1.51 (0.89–2.57)

Notes:

AFP level information are unavailable for 6 patients.

Pathological diagnosis information are unavailable for 16 patients.

Radical hepatic resection information are unavailable for 3 patients. These P-values <0.05 were highlighted in bold font.

Abbreviations: AFP, alpha-fetoprotein; HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; HR, hazard ratio; MST, median survival time; TACE, transarterial chemoembolization.

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.

Table 2

Primer sequences and HWE of VCAN exon polymorphisms

SNPs	Primer sequences	Allelea	Function	HWE p-value
rs2652098	5′-AGAAGCTGCTACCTTGCCATGA-3′5′-TGTTTCCATTCCGTGCACACAC-3′	C>T	Ser300Leu	0.035
rs309559	5′-AGTTACATGGAAGCCTGAGACTTAC-3′5′-TTGGTCAATGGCAATCTCTCCTG-3′	A>G	Lys1516Arg	0.069
rs188703	5′-AGTCTGAAAGGGAAATGACAGATTCTAC-3′5′-AGCTACTTCCTTTTCGGCTTGAA-3′	G>A	Arg1826His	0.69
rs160278	5′-CCAACAATTCAAGAGTCAGATACTGAGC-3′5′-CCCACTACCTTCAAAGATGTCTGTT-3′	A>T	Phe2301Tyr	0.15
rs160277	5′-ACCTTCAGAAGATGATGGTAAACCT-3′5′-TCATTAGAAGTCTCCAGAAGGGAAA-3′	G>T	Asp2937Tyr	0.73

Note:

Common allele > the variant allele.

Abbreviations: HWE, Hardy–Weinberg equilibrium; SNP, single-nucleotide polymorphism; VCAN, versican.

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).

Table 3

Association of VCAN exon region SNPs with survival of patients with resected early-stage HBV-associated HCC

SNPs	Best fitting model	Death/total	Crude HR (95% CI)	Adjusted HR (95% CI)a	p-value (adjusted)a
rs2652098	CC	31/63	1	1
	CT	35/42	2.36 (1.45–3.85)	2.39 (1.46–3.90)	0.001
	TT	1/1	6.19 (0.819–46.81)	5.97 (0.77–46.54)	0.088
	Dominant		2.40 (1.48–3.90)	2.42 (1.49–3.94)	<0.001
rs309559	GG	12/29	1	1
	GA	29/44	1.96 (1.00–3.85)	2.09 (1.05–4.15)	0.036
	AA	25/34	2.59 (1.29–5.16)	2.50 (1.24–5.04)	0.011
	Recessive		2.21 (1.18–4.13)	2.26 (1.19–4.28)	0.012
rs188703	AA	6/17	1	1
	AG	31/50	2.12 (0.88–5.09)	2.86 (1.16–7.01)	0.022
	GG	31/43	3.07 (1.27–7.38)	3.37 (1.38–8.22)	0.008
	Additive		1.63 (1.14–2.33)	1.58 (1.12–2.24)	0.010
rs160278	AA	11/28	1	1
	AT	31/47	2.18 (1.09–4.33)	2.32 (1.15–4.68)	0.018
	TT	25/34	2.81 (1.38–5.73)	2.70 (1.30–5.49)	0.008
	Dominant		2.42 (1.26–4.62)	2.47 (1.28–4.77)	0.007
rs160277	TT	5/16	1	1
	TG	30/49	2.40 (0.93–6.20)	3.35 (1.26–8.89)	0.015
	GG	31/43	3.45 (1.34–8.90)	3.87 (1.47–10.13)	0.006
	Recessive		2.83 (1.14–7.06)	3.61 (1.42–9.19)	0.007

Note:

Adjusted for tumor size, tumor capsule, regional invasion. These P-values <0.05 were highlighted in bold font.

Abbreviations: HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; HR, hazard ratio; SNP, single-nucleotide polymorphism; VCAN, versican.

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).

Figure 1

Association of VCAN SNP haplotypes with OS of patients with resected early-stage HBV-associated HCC.

Notes: (A) Patterns of LD plot for 5 selected VCAN SNPs. The value within each diamond represents the pairwise LD (correlation, measured as r2 values) between these 5 SNPs defined by the upper left and the upper right sides of the diamond. One block was identified based on the LD patterns. (B) Kaplan–Meier curves of OS stratified by haplotype (TAGTG) and other haplotypes in patients with resected early-stage HBV-associated HCC.

Abbreviations: HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; LD, linkage disequilibrium; OS, overall survival; SNP, single nucleotide polymorphism; VCAN, versican.

Table 4

Associations between VCAN haplotypes and OS of patients with resected early-stage HBV-associated HCC

Variables	Haplotypes	Death/total (2n)	MST (months)	Adjusted HR (95% CI)a	p-value (adjusted)a
Haplo1	CGAAT	34/62	80	1
Haplo2	CGGAG	32/51	71	1.26 (0.78–2.04)	0.35
Haplo3	CAGTG	48/59	40	1.47 (0.94–2.28)	0.091
Haplo4	TAGTG	33/40	34	2.06 (1.27–3.34)	0.003

Notes:

Adjusted for tumor size, tumor capsule, and regional invasion. Haplotype analysis included rs2652098, rs309559, rs188703, rs160278, and rs160277 SNPs. These P-values <0.05 were highlighted in bold font.

Abbreviations: HBV- associated HCC, hepatitis B virus- associated hepatocellular carcinoma; HR, hazard ratio; MST, median survival time; OS, overall survival; VCAN, versican.

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).

Figure 2

Stratified analysis of the association between VCAN haplotypes (TAGTG haplotype vs. other haplotypes) and OS of patients with resected early-stage HBV-associated HCC.

Abbreviations: AFP, alpha-fetoprotein; HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; OS, overall survival; TACE, transarterial chemoembolization; VCAN, versican.

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).

Table 5

Joint-effects analysis of VCAN haplotypes and clinicopathological characteristics from patients with resected early-stage HBV-associated HCC

Group	Variables	Death/total	Crude HR (95% CI)	Adjusted HR (95% CI)a	p-value (adjusted)a
A	Haplo1/Haplo2/Haplo3 + Tumor size ≤5 cm	23/57	1	1
B	Haplo4 + Tumor size ≤5 cm	8/11	2.13 (0.95–4.77)	2.23 (1.00–5.00)	0.051
C	Haplo1/Haplo2/Haplo3 + Tumor size >5 cm	91/125	2.19 (1.39–3.47)	1.95 (1.22–3.11)	0.005
D	Haplo4 + Tumor size >5 cm	25/29	3.55 (2.01–6.28)	3.00 (1.67–5.36)	<0.001
a	Haplo1/Haplo2/Haplo3 + Complete tumor capsule	57/105	1	1
b	Haplo4 +Complete tumor capsule	17/22	1.95 (1.13–3.35)	1.91 (1.11–3.30)	0.02
c	Haplo1/Haplo2/Haplo3 + Incomplete tumor capsule	57/77	1.78 (1.22–2.59)	1.36 (0.92–2.03)	0.12
d	Haplo4 +Incomplete tumor capsule	16/18	2.61 (1.49–4.58)	1.99 (1.12–3.55)	0.02

Notes:

Adjusted for tumor size, tumor capsule, regional invasion (excluding the stratified factor in each stratum). Haplotype analysis included rs2652098, rs309559, rs188703, rs160278, and rs160277 SNPs. These P-values <0.05 were highlighted in bold font.

Abbreviations: HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; HR, hazard ratio; VCAN, versican.

Figure 3

Kaplan–Meier curve of joint-effects analysis in HBV-associated HCC.

Notes: (A) OS stratified by VCAN haplotypes and tumor size cut-off by 5 cm; (B) OS stratified by VCAN haplotypes and tumor capsule complete or not.

Abbreviations: HBV, hepatitis B virus; HCC, hepatocellular carcinoma; OS, overall survival; VCAN, versican.

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).

Figure 4

Kaplan–Meier curve prediction of poor OS of patients with resected early-stage HBV-associated HCC stratified by high VCAN expression.

Abbreviations: HBV-associated HCC, hepatitis B virus-associated hepatocellular carcinoma; OS, overall survival; VCAN, versican.

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.