Skip to main content
  • Primary research
  • Open access
  • Published:

Expression of the Wnt ligands gene family and its relationship to prognosis in hepatocellular carcinoma

Abstract

Background

The Wnt gene family members are known to participate regulating various normal and pathological processes including tumorigenesis. However, the association between Wnt ligands gene family and prognosis in hepatocellular carcinoma has not been systematically studied. Therefore, we evaluated the role of Wnt ligands gene family in hepatocellular carcinoma using publicly available data from The Cancer Genome Atlas (TCGA).

Methods

Clinical information and RNA-Seq mRNA expression data were derived from TCGA hepatocellular carcinoma cohort. Differences in overall survival (OS) and disease-free survival (DFS) between increased and decreased expression groups (defined by X-tile analyses) of Wnt ligands gene family were compared using Kaplan–Meier method and Cox regression model, with p-values calculated via log-rank test. Gene Set Enrichment Analysis (GSEA) was performed.

Results

Multivariate analysis adjusted for patient age, sex, BMI, tumor grade, and TMN stage revealed that Wnt1, Wnt3 and Wnt5B expressions were independent prognostic factors for OS and DFS (OS: HR = 0.58, P = 0.006; HR = 0.65, P = 0.03; HR = 0.56, P = 0.023, respectively; DFS: HR = 0.52, P < 0.001; HR = 1.93, P = 0.003; HR = 0.59, P = 0.011, respectively). Furthermore, expression of Wnt1 and Wnt5B was significantly associated with TMN stage (P = 0.02 and P = 0.03 for OS; P = 0.02 and P = 0.02 for DFS). GSEA showed that nucleotide excision repair was differentially enriched in Wnt1 low expression phenotype and aminoacyl trna biosynthesis and basal transcription factors were differentially enriched in Wnt5B low expression phenotype.

Conclusions

Our results identified associations of several Wnt ligands with prognosis of HCC patients, indicating that these genes could serve as prognostic biomarkers of HCC.

Background

Liver cancer is one of the most commonly diagnosed cancer and the second most frequent cause of cancer-related deaths globally, with 854,000 incident cases and 810,000 deaths in 2015 [1]. Hepatocellular carcinoma (HCC), as the most common pathological type, accounts for approximately 90% of primary liver cancers [2]. Chronic viral hepatitis (B and C), alcohol intake and aflatoxin exposure were well known underlying aetiologies. A series of therapies including liver resection, percutaneous ethanol injection, transcatheter arterial chemoembolization, microwave ablation, liver transplantation and systemic therapy were developed for patients with HCC. However, the mortality rate of HCC was always very high due to their disease being diagnosed at a late stage [3]. Detailed underlying mechanisms of development and progression of HCC were considered complicated and ambiguous. Studying the genes that play a key role in HCC development is crucial to identify disease biomarkers which could be of great use for diagnosis, prognostic prediction or even development of targeted drugs.

The Wnt family of 19 secreted glycoproteins have crucial roles in the regulation of diverse processes, such as embryogenesis, differentiation, and tumorigenesis through canonical dependent and non-canonical pathways [4, 5]. Numerous studies have reported that aberrant activation of Wnt signaling may contribute to the pathology of various types of cancer, including colon cancer [6], gastric cancer [7], HCC [8]. A previous study, using glutamine synthetase (encoded by canonical Wnt signaling target GLUL gene) as a sensitive and specific marker, showed that 36% HCCs displayed canonical Wnt activation [9]. Recently, Wnt3a expression in HCC is reported to be associated with the poorly-differentiated grade, liver cirrhosis, HBV infection, higher TNM stage, and a relatively shorter survival time [10]. Moreover, several previous investigations demonstrated that an increased Wnt1 expression was detected in human HCC tissue and human hepatoma cell lines and correlated with increased tumor recurrence after curative tumor resection [11, 12]. In addition to Wnt3a and Wnt1, increased expression of Wnt3, Wnt4, Wnt5a, and Wnt10b has also been revealed in HCC tumors and in the peritumoral liver tissues [8, 13, 14]. However, a comprehensive analysis of the association between the expression of Wnt ligands and clinicopathologic features of hepatocellular carcinoma (HCC) is lacking. Thus, the objective of the current study was to extensively evaluate the prognostic value of Wnt ligands in HCC based on data obtained from TCGA. To gain further insight into the biological pathways involved in HCC pathogenesis related Wnt regulatory network, GSEA was also performed.

Materials and methods

The level-3 expression data (RNA-seqV2) and clinicopathological data of 360 HCC patients and 50adjacent normal liver samples were downloaded from The Cancer Genome Atlas (TCGA, https://tcga-data.nci.nih.gov/tcga/) data portal. The clinicopathological characteristics of HCC patients, including age, sex, race, body mass index (BMI), tumor grade, tumor, node, metastasis (TNM) stage, overall survival time, overall survival status, disease free survival time and disease-free survival status, were collected. The methods of biospecimen collection, RNA isolation, and RNA sequencing were previously described by the Cancer Genome Atlas Research Network [15]. The data was processed according to the TCGA publication guidelines and data access policies.

Gene set enrichment analysis

To determine whether a priori defined set of genes shows statistically significant, consistent differences between two biological states (Increased expression vs. Decreased expression), Gene Set Enrichment Analysis (GSEA) was performed by the JAVA program (http://software.broadinstitute.org/gsea/downloads.jsp) using the MSigDB C2 KEGG pathways gene sets, which contains 186 gene sets. Normalized enrichment score (NES), nominal p value and false discovery rate (FDR) were used to quantify enrichment magnitude and statistical significance, respectively [16].

Statistical analysis

Data analysis was performed with SPSS (Version 22.0; IBM, New York, NY, USA). Differences between groups were calculated by using Chi square test or Fisher exact test. Box-plots were used to visualize expression differences for Wnt ligands between tumor and adjacent non-tumor tissues. The optimal cut-off values for Wnt ligands expression were determined by X-tile software (Version3.6.1, Yale University, New Haven, CT, USA) [17]. Subsequently, each expression level of Wnt ligands was divided into increased expression group and decreased expression group according to the optimal cut-off value. Chi square test was also applied to evaluate the association between Wnt ligands expression and the clinicopathologic features in HCC. The Kaplan–Meier survival analysis and log-rank test were used to compare differences in survival times. Univariate and multivariate survival analysis were performed using the Cox hazards regression model to analyze the independent parameters associated to the overall survival and disease-free survival of HCC patients. A P value less than 0.05 was considered to be statistically significant.

Results

Basic patient characteristics

Clinical information and Wnt ligands mRNA expression levels of 360 patients were obtained from TCGA. Among 360 patients, 310 patients were recorded with diseases free survival status. Detailed clinical characteristics of the 360 patients in the TCGA database are shown in Table 1. TNM stage was significantly associated with the OS and DFS (P < 0.001), but not sex, age, BMI, or race (all P > 0.05).

Table 1 Basic characteristics of TCGA HCC patients

Wnt ligands expression changes in HCC

All 19 members of the Wnt ligands family (Wnt1-Wnt16) had their expression analyzed in 360 HCC tumor tissues and 50 adjacent non-tumor tissues. The expression levels of Wnt2B, Wnt3A, Wnt6, Wnt8B and Wnt10B were significantly higher in primary liver tumor tissues than which in adjacent non-tumor tissues. However, the expression levels of Wnt2, Wnt5B, Wnt7A, Wnt7B, Wnt9A and Wnt11 were dramatically decreased in liver tumor tissues. Expression levels the remaining Wnt ligands, including Wnt1, Wnt3, Wnt4, Wnt5A, Wnt8A, Wnt9B, Wnt10A and Wnt16, remained insignificantly different between liver tumors tissues and adjacent non-tumor tissues. The box diagrams indicating the distribution of gene expression in HCC patients and adjacent normal for all Wnt ligands members are displayed in Additional file 1: Figure S1.

Multivariate analysis and survival outcomes

The expression levels of Wnt ligands were divided into increased expression group and decreased expression group according to the cut-off values determined by X-tile program. The univariate and multivariate analysis demonstrated that TNM stages (HR = 2.02, 95% CI 1.37–2.96, P < 0.001), decreased Wnt1 expression (HR = 0.58, 95% CI 0.39–0.85, P = 0.006), decreased Wnt3 expression (HR = 0.65, 95% CI 0.45–0.96, P = 0.030), decreased Wnt5B expression (HR = 0.56, 95% CI 0.34–0.93, P = 0.023), increased Wnt6 expression (HR = 1.62, 95% CI 1.05–2.51, P = 0.030) and increased Wnt8A expression (HR = 4.14, 95% CI 1.47–11.68, P = 0.007) were independent prognostic factors for overall survival (Table 2) and TNM stages (HR = 1.41, 95% CI 1.24–1.76, P < 0.001), decreased Wnt1 expression (HR = 0.52, 95% CI 0.35–0.75, P < 0.001), increased Wnt3 expression (HR = 1.93, 95% CI 1.24–2.98, P = 0.003), decreased Wnt5A expression (HR = 0.61, 95% CI 0.40–0.91, P = 0.017), decreased Wnt5B expression (HR = 0.59, 95% CI 0.39–0.89, P = 0.011) and increased Wnt8B expression (HR = 0.53, 95% CI 0.36–0.77, P = 0.001) were independent prognostic factors for disease-free survival (Additional file 2: Table S1). Consistent with multivariate analysis, Kaplan–Meier survival analysis showed that HCC patients with advanced TNM stage, decreased Wnt1 expression, decreased Wnt3 expression, decreased Wnt5B expression, increased Wnt6 expression and increased Wnt8A expression had a worse OS than the counterpart (all P < 0.05) (Fig. 1). Kaplan–Meier survival analysis for DFS was shown in Additional file 3: Figure S2).

Table 2 Univariate and multivariate analysis of overall survival using the Cox proportional hazard regression model
Fig. 1
figure 1

TMN stage (a) and expression of Wnt1, Wnt3, Wnt5B, Wnt6 and Wnt8A are associated with overall survival (b–f). Kaplan–Meier survival analysis and log-rank test were used to compare differences in overall survival between the groups classified using cut-off values determined by X-tile

Association with Wnt ligands expression and clinicopathologic variables

A total of 360 HCC patients with OS data and 310 HCC patients with DFS data were analyzed from TCGA. As showed in Table 3, low expression of Wnt1 and Wnt5B significantly correlated with the advanced TMN stage (Wnt1:P = 0.02; Wnt5B:P = 0.03). Low expression of Wnt3 correlated with histological grade (P = 0.04) but not TMN stage. Similar results were observed in the analysis for DFS. Low expressions of Wnt1 and Wnt5B correlated significantly with the advanced TMN stage (Wnt1:P = 0.02; Wnt5B:P = 0.02). Details were displayed in Additional file 4: Table S2).

Table 3 Correlation of clinicopathologic variables and expression of several specific Wnts significantly associated with OS

GSEA identifies a Wnt1-related and Wnt5B-related KEGG signaling pathway

To identify KEGG signaling pathways that are differentially activated in HCC, we conducted GSEA between increased and decreased Wnt1 and Wnt5B expression data sets. GSEA reveal significant differences (NOM p-val < 0.01) in enrichment of MSigDB Collection (c2.cp.kegg.v6.1.entrez). We selected the most significantly enriched signaling pathways based on their normalized enrichment score (NES) (Fig. 2 and Table 4). The Fig. 2 shows that nucleotide excision repair is differentially enriched in Wnt1 low expression phenotype and aminoacyl tRNA biosynthesis and basal transcription factors are differentially enriched in Wnt5B low expression phenotype.

Fig. 2
figure 2

Enrichment plots from gene set enrichment analysis (GSEA). GSEA results showing nucleotide excision repair is differentially enriched in Wnt1 decreased expression phenotype (a) and aminoacyl tRNA biosynthesis (b) and basal transcription factors (c) are differentially enriched in Wnt5B decreased expression phenotype (b, c)

Table 4 Gene sets enriched in decreased expression phenotype

Discussion

In this study, we investigated the association between Wnt ligand family genes and HCC. We observed that the mRNA expression levels of several specific Wnt ligand family genes, such asWnt1, Wnt3 and Wnt5B, are associated with distinct OS and DFS. Moreover, we found that expression of Wnt1 and Wnt5B significantly correlated with TMN stage. Thus, Wnt ligand family genes-especially Wnt1, Wnt3 and Wnt5B-may serve as prognostic biomarkers of HCC and represent possible oncogenes that could serve as therapeutic targets of HCC.

Wnt signals are known for regulation of diverse processes, including cell proliferation, survival, migration and polarization, embryonic development, specification of cell fate, and self-renewal in stem cells [18]. It is not surprising that Wnt pathway mutations are frequently observed in carcinomas. In the past decade, a large number of studies have been conducted to explore the role of Wnts and their downstream effectors in regulating cancer progression, including tumor initiation, tumor growth, cell senescence, cell death, differentiation and metastasis [18]. As activation of Wnt signals starts with the secretion of Wnt ligands, accumulating researches have investigated the expression profile of all 19 Wnt ligand genes indifferent cancer cell types, such as mammary carcinoma cell lines, human ovarian cancer cell lines and HCC cell lines [19,20,21,22]. Expectedly, Wnts and Wnt pathway components are frequently over- or under-expressed in different human malignant tumors. Interestingly, the expression patterns of Wnt signaling components can also serve as prognostic indicators of patient outcomes. For example, Wnt3a expression is reported to be significantly associated with poor prognosis of numerous cancers, including esophageal squamous cell carcinoma and HCC [10, 23]. In gastric cancer and ovarian carcinoma, increased level of Wnt5A protein was associated with high grade tumors and with decreased patient survival [24, 25], yet in colon cancer and HCC high level of Wnt5A protein correlated with increased patient survival [26, 27]. A recent study indicated that high Wnt2 expression in fibroblasts is associated with poor prognosis in human colorectal cancer [28]. However, a comprehensive analysis of the association between Wnts and liver cancer prognosis has not been performed.

In our study, it was found that increased expressions of Wnt 2B, Wnt3A, Wnt6, Wnt8B and Wnt10B in primary liver tumors and decreased expressions of Wnt2, Wnt5B, Wnt7A, Wnt7B, Wnt9A and Wnt11 in liver tumors tissues. These Wnts might have biomarker potential and could be utilized clinically in a diagnostic capacity. Interestingly, our study revealed that increased expressions of Wnt1, Wnt3, and Wnt5B and decreased expressions of Wnt6 and Wnt8A in HCC were associated with good OS probability. Increased expression of Wnt1, Wnt5A, Wnt5B and Wnt8B and decreased expression of Wnt3 were found to be associated with good DFS probability. The results showed that Wnt1, Wnt3 and Wnt5B were both independent prognostic factors affecting OS and DFS. Noteworthily, expression level of Wnt3 differed greatly for good OS and DFS probability. This reverse result could relate to the different cut-off values calculated by X-tile. As showed in Table 2 and Additional file 2: Table S1, the cut-off values of Wnt3 expression for OS and DFS were 37.97 and 98.6, respectively. Patients with low Wnt3 expression (< 37.97) had a worse OS and patients with low Wnt3 expression (< 98.6) had a better DFS. It could be explained that patients with Wnt3 expression level between 37.97 and 98.6 had a significantly better survival time and then contribute to this discrepancy.

Wnt1, one of the key ligands in β-catenin regulation, has been described for its prognostic role in several types of malignant tumors including non-small cell lung cancer [29], renal cell carcinoma [30], and colorectal cancer [31]. More importantly, a report by Lee et al. [11] found that high tumor Wnt1 expression was associated with increased hepatitis B virus (HBV)—related and hepatitis C virus (HCV)—related HCC recurrence after curative resection. The mechanism might correlate with increased nuclearβ-catenin accumulation accompanied by decreased membranous-cadherin expression and this seems to be consistent with the character of cancer cell metastasis. Surprisingly, our data showed that increased expression of Wnt1 was associated with good OS and DFS probability. The cause of this discrepancy between the previous study and our study remains unknown. It might be explained by the biological heterogeneity of HCC, which has an important impact on carcinogenesis and development.

Among various Wnts associated with the canonical Wnt/-catenin signaling pathway, Wnt3 was previously reported being frequently up-regulated in human cancers [8, 32, 33]. However, recent researches demonstrated that Wnt3 was down-regulated in some types of cancers and related to worse prognosis [34, 35]. With regards to HCC, Wnt3 expression was reported to be upregulated in human HCC compared to the adjacent peritumoral tissue [13]. In our study, the difference of Wnt3 expression remains insignificant between tumors tissue and adjacent non-tumor tissue. Moreover, since the cut-off values measured by X-tile differ too much for OS and DFS, the prognostic value of Wnt3 expression in HCC for a good OS and DFS was discrepant. Hence, prognostic value of Wnt3 expression in HCC remains further elucidation.

Wnt5B, an intermediately transforming or non-transforming Wnt family member, was reported to activate the noncanonical Wnt signals [36]. Several studies have demonstrated that Wnt5B is involved in the proliferation and migration of tumor cells and have critical role in tumor lymph angiogenesis and lymph node metastasis through the regulation of epithelial-mesenchymal transition (EMT) [37, 38]. Recently, an investigation from China reported that Wnt5b mRNA expression was significantly higher in hepatitis B virus-related HCC tissues than that of adjacent noncancerous tissues. Patients with up-regulated Wnt5b mRNA and protein had a shorter relapse-free survival [39]. However, our study indicated that increased expression of Wnt5b was significantly associated with better OS and DFS. Hence, similar to Wnt1, this discrepancy was observed in Wnt5B as well. We have a hypothesis that Wnt5B, who shows the greatest similarity with Wnt5A, may share a common Fzd receptor2 [40] and mediates similar Wnt5A effects of antagonizing Wnt signaling which may further inhibite HCC proliferation and migration [21, 37]. In addition, in our study, the optimal cut-off values for Wnt ligands expression were measured by X-tile, which presents substantial tumor subpopulations and shows the biological relationships between a biomarker and outcome. We believe the different methods to determine the cut-off value could also explain the discrepancy.

In recent years, a number of novel biomarkers, such as AFP [41, 42], glypican 3 [43], Osteopontin [44], CXCL1 [45], UQCRH [46], TIP30 [47], neutrophil–lymphocyte ratio [48], have emerged for diagnosing HCC and predicting patient outcome. Among these biomarkers, AFP is a well-known serum diagnostic and prognostic biomarker for HCC [49]. However, its prognostic value remains controversial. Several previous literature refuted the prognostic value of AFP in single, small HCC, and even for the prediction of HCC recurrence [50, 51]. Prognostic biomarkers’ use in the daily practice has never been endorsed by international guidelines. These frustrating results should not discourage the study of novel biomarkers and their translation at tissue level into prognostic and predictive indicators. In the present population, our work reveals a correlation between HCC outcome and Wnt ligand family genes. Hence, we postulate that the Wnt ligand family members may serve as potential serum biomarkers for prognosis of HCC.

Since the prognostic value of Wnt3 expression for OS and DFS is discrepant, we merely performed GSEA analysis to further investigate the potential mechanism of Wnt1 and Wnt5B in HCC. GSEA analysis showed that nucleotide excision repair is differentially enriched in Wnt1 low expression phenotype and aminoacyl tRNA biosynthesis and basal transcription factors are differentially enriched in Wnt5B low expression phenotype. Nucleotide excision repair (NER) is a repair system for many types of DNA damage, and therefore many types of genotoxic carcinogenic exposures, including ultraviolet light, products of organic combustion, metals, oxidative stress, etc. There have been a few reports demonstrating association of nucleotide excision repair with cancer [52, 53]. Ishikawa et al. [54] previously reported that the DNA repair system, especially the NER pathway, played a vital role in protection against human cancer. Moreover, Wang et al. [55] demonstrated that six polymorphisms of five genes involved in three steps of nucleotide excision repair pathways was associated with hepatocellular cancer risk. But the associations between Wnt1 and Wnt5B expression and nucleotide excision repair, aminoacyl tRNA biosynthesis and basal transcription factors were the first to be reported, and the regulatory mechanism needs to be further elucidated.

Nevertheless, the prediction of protein expression using mRNA was far from perfect. Because of limitations in our study design, the correlation between Wnts mRNA expression and Wnts protein expression could not be clearly assessed in this study. Further study in HCC is required.

Conclusion

Our study identifies that several Wnt ligand family genes, such as Wnt1, Wnt3 and Wnt5B, may be potential prognostic biomarkers of HCC. Moreover, the nucleotide excision repair, aminoacyl tRNA biosynthesis and basal transcription factors may be the key pathway regulated by Wnts in HCC. Due to the small sample size and incomplete clinical information in this study, further experimental validation should be performed to prove the biologic impact of Wnts.

References

  1. Global Burden of Disease Liver Cancer C, Akinyemiju T, Abera S, Ahmed M, Alam N, Alemayohu MA, et al. The burden of primary liver cancer and underlying etiologies from 1990 to 2015 at the global, regional, and national level: results from the global burden of disease study 2015. JAMA Oncol. 2017;3(12):1683–91.

    Article  Google Scholar 

  2. European Association for the Study of the Liver. Electronic address eee, European Association for the Study of the L. EASL clinical practice guidelines: management of hepatocellular carcinoma. J Hepatol. 2018;69(1):182–236.

    Article  Google Scholar 

  3. Allaire M, Nault JC. Advances in management of hepatocellular carcinoma. Curr Opin Oncol. 2017;29(4):288–95.

    Article  CAS  Google Scholar 

  4. Dale TC. Signal transduction by the Wnt family of ligands. Biochem J. 1998;329(Pt 2):209–23.

    Article  CAS  Google Scholar 

  5. Neth P, Ciccarella M, Egea V, Hoelters J, Jochum M, Ries C. Wnt signaling regulates the invasion capacity of human mesenchymal stem cells. Stem Cells. 2006;24(8):1892–903.

    Article  CAS  Google Scholar 

  6. Zhang X, Gaspard JP, Chung DC. Regulation of vascular endothelial growth factor by the Wnt and K-ras pathways in colonic neoplasia. Cancer Res. 2001;61(16):6050–4.

    CAS  Google Scholar 

  7. Yong X, Tang B, Xiao YF, Xie R, Qin Y, Luo G, et al. Helicobacter pylori upregulates Nanog and Oct4 via Wnt/beta-catenin signaling pathway to promote cancer stem cell-like properties in human gastric cancer. Cancer Lett. 2016;374(2):292–303.

    Article  CAS  Google Scholar 

  8. Bengochea A, de Souza MM, Lefrancois L, Le Roux E, Galy O, Chemin I, et al. Common dysregulation of Wnt/Frizzled receptor elements in human hepatocellular carcinoma. Br J Cancer. 2008;99(1):143–50.

    Article  CAS  Google Scholar 

  9. Audard V, Grimber G, Elie C, Radenen B, Audebourg A, Letourneur F, et al. Cholestasis is a marker for hepatocellular carcinomas displaying beta-catenin mutations. J Pathol. 2007;212(3):345–52.

    Article  CAS  Google Scholar 

  10. Pan LH, Yao M, Cai Y, Gu JJ, Yang XL, Wang L, et al. Oncogenic Wnt3a expression as an estimable prognostic marker for hepatocellular carcinoma. World J Gastroenterol. 2016;22(14):3829–36.

    Article  CAS  Google Scholar 

  11. Lee HH, Uen YH, Tian YF, Sun CS, Sheu MJ, Kuo HT, et al. Wnt-1 protein as a prognostic biomarker for hepatitis B-related and hepatitis C-related hepatocellular carcinoma after surgery. Cancer Epidemiol Biomarkers Prev. 2009;18(5):1562–9.

    Article  CAS  Google Scholar 

  12. Wei W, Chua MS, Grepper S, So SK. Blockade of Wnt-1 signaling leads to anti-tumor effects in hepatocellular carcinoma cells. Mol Cancer. 2009;8:76.

    Article  Google Scholar 

  13. Kim M, Lee HC, Tsedensodnom O, Hartley R, Lim YS, Yu E, et al. Functional interaction between Wnt3 and Frizzled-7 leads to activation of the Wnt/beta-catenin signaling pathway in hepatocellular carcinoma cells. J Hepatol. 2008;48(5):780–91.

    Article  CAS  Google Scholar 

  14. Yoshikawa H, Matsubara K, Zhou X, Okamura S, Kubo T, Murase Y, et al. WNT10B functional dualism: beta-catenin/Tcf-dependent growth promotion or independent suppression with deregulated expression in cancer. Mol Biol Cell. 2007;18(11):4292–303.

    Article  CAS  Google Scholar 

  15. Cancer Genome Atlas Research Network. Electronic address w b e, cancer genome atlas research N. Comprehensive and integrative genomic characterization of hepatocellular carcinoma. Cell. 2017;169(7):1327-41 e23.

    Google Scholar 

  16. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50.

    Article  CAS  Google Scholar 

  17. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10(21):7252–9.

    Article  CAS  Google Scholar 

  18. Anastas JN, Moon RT. WNT signalling pathways as therapeutic targets in cancer. Nat Rev Cancer. 2013;13(1):11–26.

    Article  CAS  Google Scholar 

  19. Vouyovitch CM, Perry JK, Liu DX, Bezin L, Vilain E, Diaz JJ, et al. WNT4 mediates the autocrine effects of growth hormone in mammary carcinoma cells. Endocr Relat Cancer. 2016;23(7):571–85.

    Article  CAS  Google Scholar 

  20. Yoshioka S, King ML, Ran S, Okuda H, MacLean JA, McAsey ME, et al. WNT7A regulates tumor growth and progression in ovarian cancer through the WNT/beta-catenin pathway. Mol Cancer Res. 2012;10(3):469–82.

    Article  CAS  Google Scholar 

  21. Yuzugullu H, Benhaj K, Ozturk N, Senturk S, Celik E, Toylu A, et al. Canonical Wnt signaling is antagonized by noncanonical Wnt5a in hepatocellular carcinoma cells. Mol Cancer. 2009;8:90.

    Article  Google Scholar 

  22. Wang W, Xu L, Liu P, Jairam K, Yin Y, Chen K, et al. Blocking Wnt secretion reduces growth of hepatocellular carcinoma cell lines mostly independent of beta-catenin signaling. Neoplasia. 2016;18(12):711–23.

    Article  CAS  Google Scholar 

  23. Oguma J, Ozawa S, Kazuno A, Nitta M, Ninomiya Y, Kajiwara H. Wnt3a expression is associated with poor prognosis of esophageal squamous cell carcinoma. Oncol Lett. 2018;15(3):3100–8.

    PubMed  Google Scholar 

  24. Kurayoshi M, Oue N, Yamamoto H, Kishida M, Inoue A, Asahara T, et al. Expression of Wnt-5a is correlated with aggressiveness of gastric cancer by stimulating cell migration and invasion. Cancer Res. 2006;66(21):10439–48.

    Article  CAS  Google Scholar 

  25. Peng C, Zhang X, Yu H, Wu D, Zheng J. Wnt5a as a predictor in poor clinical outcome of patients and a mediator in chemoresistance of ovarian cancer. Int J Gynecol Cancer. 2011;21(2):280–8.

    Article  Google Scholar 

  26. Dejmek J, Dejmek A, Safholm A, Sjolander A, Andersson T. Wnt-5a protein expression in primary dukes B colon cancers identifies a subgroup of patients with good prognosis. Cancer Res. 2005;65(20):9142–6.

    Article  CAS  Google Scholar 

  27. Wang L, Yao M, Fang M, Zheng WJ, Dong ZZ, Pan LH, et al. Expression of hepatic Wnt5a and its clinicopathological features in patients with hepatocellular carcinoma. Hepatobiliary Pancreat Dis Int. 2018;17(3):227–32.

    Article  Google Scholar 

  28. Kramer N, Schmollerl J, Unger C, Nivarthi H, Rudisch A, Unterleuthner D, et al. Autocrine WNT2 signaling in fibroblasts promotes colorectal cancer progression. Oncogene. 2017;36(39):5460–72.

    Article  CAS  Google Scholar 

  29. Huang CL, Liu D, Ishikawa S, Nakashima T, Nakashima N, Yokomise H, et al. Wnt1 overexpression promotes tumour progression in non-small cell lung cancer. Eur J Cancer. 2008;44(17):2680–8.

    Article  CAS  Google Scholar 

  30. Kruck S, Eyrich C, Scharpf M, Sievert KD, Fend F, Stenzl A, et al. Impact of an altered Wnt1/beta-catenin expression on clinicopathology and prognosis in clear cell renal cell carcinoma. Int J Mol Sci. 2013;14(6):10944–57.

    Article  Google Scholar 

  31. Wang J, Lu R, Fu X, Dan Z, Zhang YG, Chang X, et al. Novel regulatory roles of Wnt1 in infection-associated colorectal cancer. Neoplasia. 2018;20(5):499–509.

    Article  CAS  Google Scholar 

  32. Nakashima N, Liu D, Huang CL, Ueno M, Zhang X, Yokomise H. Wnt3 gene expression promotes tumor progression in non-small cell lung cancer. Lung Cancer. 2012;76(2):228–34.

    Article  Google Scholar 

  33. Wang HS, Nie X, Wu RB, Yuan HW, Ma YH, Liu XL, et al. Downregulation of human Wnt3 in gastric cancer suppresses cell proliferation and induces apoptosis. Onco Targets Ther. 2016;9:3849–60.

    Article  CAS  Google Scholar 

  34. Do Carmo NG, Sakamoto LH, Pogue R, Mascarenhas CD, Passos SK, Felipe MS, et al. Altered expression of PRKX, WNT3 and WNT16 in human nodular basal cell carcinoma. Anticancer Res. 2016;36(9):4545–51.

    Article  CAS  Google Scholar 

  35. Poppova L, Janovska P, Plevova K, Radova L, Plesingerova H, Borsky M, et al. Decreased WNT3 expression in chronic lymphocytic leukaemia is a hallmark of disease progression and identifies patients with worse prognosis in the subgroup with mutated IGHV. Br J Haematol. 2016;175(5):851–9.

    Article  CAS  Google Scholar 

  36. Kikuchi A, Yamamoto H. Tumor formation due to abnormalities in the beta-catenin-independent pathway of Wnt signaling. Cancer Sci. 2008;99(2):202–8.

    Article  CAS  Google Scholar 

  37. Gujral TS, Chan M, Peshkin L, Sorger PK, Kirschner MW, MacBeath G. A noncanonical Frizzled2 pathway regulates epithelial-mesenchymal transition and metastasis. Cell. 2014;159(4):844–56.

    Article  CAS  Google Scholar 

  38. Wang SH, Chang JS, Hsiao JR, Yen YC, Jiang SS, Liu SH, et al. Tumour cell-derived WNT5B modulates in vitro lymphangiogenesis via induction of partial endothelial-mesenchymal transition of lymphatic endothelial cells. Oncogene. 2017;36(11):1503–15.

    Article  CAS  Google Scholar 

  39. Xiang LY, Ou HH, Chen ZJ, Fang YH, Huang Y, Li XH, et al. Expression of Wnt5b in patients with HBV-related hepatocellular carcinoma and its clinical significance. Nan Fang Yi Ke Da Xue Xue Bao. 2017;37(8):1071–7.

    CAS  PubMed  Google Scholar 

  40. Saitoh T, Katoh M. Molecular cloning and characterization of human WNT5B on chromosome 12p13.3 region. Int J Oncol. 2001;19(2):347–51.

    CAS  PubMed  Google Scholar 

  41. Komorowski AL, Hsu CC, Julka KD, Vasavada B, Lin CC, Wang CC, et al. AFP role in predicting recurrence of hepatocellular carcinoma after living donor liver transplantation in HCV patients. Neoplasma. 2018;65(3):455–60.

    Article  Google Scholar 

  42. Bai DS, Zhang C, Chen P, Jin SJ, Jiang GQ. The prognostic correlation of AFP level at diagnosis with pathological grade, progression, and survival of patients with hepatocellular carcinoma. Sci Rep. 2017;7(1):12870.

    Article  Google Scholar 

  43. Liu H, Yang C, Lu W, Zeng Y. Prognostic significance of glypican-3 expression in hepatocellular carcinoma: a meta-analysis. Medicine (Baltimore). 2018;97(4):e9702.

    Article  CAS  Google Scholar 

  44. Cabiati M, Gaggini M, Cesare MM, Caselli C, De Simone P, Filipponi F, et al. Osteopontin in hepatocellular carcinoma: A possible biomarker for diagnosis and follow-up. Cytokine. 2017;99:59–65.

    Article  CAS  Google Scholar 

  45. Han KQ, Han H, He XQ, Wang L, Guo XD, Zhang XM, et al. Chemokine CXCL1 may serve as a potential molecular target for hepatocellular carcinoma. Cancer Med. 2016;5(10):2861–71.

    Article  CAS  Google Scholar 

  46. Park ER, Kim SB, Lee JS, Kim YH, Lee DH, Cho EH, et al. The mitochondrial hinge protein, UQCRH, is a novel prognostic factor for hepatocellular carcinoma. Cancer Med. 2017;6(4):749–60.

    Article  CAS  Google Scholar 

  47. Zhang X, Lv L, Ouyang X, Zhang S, Fang J, Cai L, et al. Association of TIP30 expression and prognosis of hepatocellular carcinoma in patients with HBV infection. Cancer Med. 2016;5(9):2180–9.

    Article  CAS  Google Scholar 

  48. Wang Y, Peng C, Cheng Z, Wang X, Wu L, Li J, et al. The prognostic significance of preoperative neutrophil-lymphocyte ratio in patients with hepatocellular carcinoma receiving hepatectomy: a systematic review and meta-analysis. Int J Surg. 2018;55:73–80.

    Article  Google Scholar 

  49. Chaiteerakij R, Addissie BD, Roberts LR. Update on biomarkers of hepatocellular carcinoma. Clin Gastroenterol Hepatol. 2015;13(2):237–45.

    Article  CAS  Google Scholar 

  50. Giannini EG, Marenco S, Borgonovo G, Savarino V, Farinati F, Del Poggio P, et al. Alpha-fetoprotein has no prognostic role in small hepatocellular carcinoma identified during surveillance in compensated cirrhosis. Hepatology. 2012;56(4):1371–9.

    Article  CAS  Google Scholar 

  51. Shim JH, Yoon DL, Han S, Lee YJ, Lee SG, Kim KM, et al. Is serum alpha-fetoprotein useful for predicting recurrence and mortality specific to hepatocellular carcinoma after hepatectomy? A test based on propensity scores and competing risks analysis. Ann Surg Oncol. 2012;19(12):3687–96.

    Article  Google Scholar 

  52. Vineis P, Manuguerra M, Kavvoura FK, Guarrera S, Allione A, Rosa F, et al. A field synopsis on low-penetrance variants in DNA repair genes and cancer susceptibility. J Natl Cancer Inst. 2009;101(1):24–36.

    Article  CAS  Google Scholar 

  53. Marteijn JA, Lans H, Vermeulen W, Hoeijmakers JH. Understanding nucleotide excision repair and its roles in cancer and ageing. Nat Rev Mol Cell Biol. 2014;15(7):465–81.

    Article  CAS  Google Scholar 

  54. Ishikawa T, Zhang SS, Qin X, Takahashi Y, Oda H, Nakatsuru Y, et al. DNA repair and cancer: lessons from mutant mouse models. Cancer Sci. 2004;95(2):112–7.

    Article  CAS  Google Scholar 

  55. Wang B, Xu Q, Yang HW, Sun LP, Yuan Y. The association of six polymorphisms of five genes involved in three steps of nucleotide excision repair pathways with hepatocellular cancer risk. Oncotarget. 2016;7(15):20357–67.

    PubMed  PubMed Central  Google Scholar 

Download references

Authors’ contributions

JJD and LY draft the manuscript and analyzed the data. KQS designed the study and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

None.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and/or analysed during the current study are available in the TCGA Research Network: http://cancergenome.nih.gov/.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Because the dataset in the current study was downloaded from TCGA, and data acquiring and application complied with the TCGA publication guidelines and data access policies, additional approval by an ethics committee and consent to participate were not needed.

Funding

This study was supported by Natural Science Foundation of Zhejiang Province (LY16H160047 and LY17H160057) and Provinces and Ministries Co-Contribution of Zhejiang, China (No.wkj-zj-1706).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ke-Qing Shi.

Additional files

Additional file 1: Figure S1.

Comparison of gene expression levels of all 19 members of Wnt family genes in tumor tissue and adjacent non-tumor tissue.

Additional file 2: Table S1.

Univariate and multivariate analysis of disease-free survival using the Cox proportional hazard regression model.

Additional file 3: Figure S2.

TMN stage (A) and expression of Wnt1, Wnt3, Wnt5A, Wnt5B and Wnt8B (B–F) are associated with disease-free survival. Kaplan–Meier survival analysis and log-rank test were used to compare differences in overall survival between the groups classified using cut-off values determined by X-tile.

Additional file 4: Table S2.

Correlation of clinicopathologic variables and expression of several specific Wnts significantly associated with DFS.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Dong, JJ., Ying, L. & Shi, KQ. Expression of the Wnt ligands gene family and its relationship to prognosis in hepatocellular carcinoma. Cancer Cell Int 19, 34 (2019). https://doi.org/10.1186/s12935-019-0743-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12935-019-0743-z

Keywords