Research Article | Open Access
Chengkun Yang, Chuangye Han, Xiangkun Wang, Xiwen Liao, Xiaoguang Liu, Wei Qin, Long Yu, Guangzhi Zhu, Hao Su, Sicong Lu, Zhiwei Chen, Tingdong Yu, Zhen Liu, Ketuan Huang, Zhengtao Liu, Yu Liang, Jianlu Huang, Tao Peng, "Clinical Implication and the Hereditary Factors of NM23 in Hepatocellular Carcinoma Based on Bioinformatics Analysis and Genome–Wide Association Study", Journal of Oncology, vol. 2018, Article ID 6594169, 18 pages, 2018. https://doi.org/10.1155/2018/6594169
Clinical Implication and the Hereditary Factors of NM23 in Hepatocellular Carcinoma Based on Bioinformatics Analysis and Genome–Wide Association Study
NM23 expression is closely associated with hepatocellular carcinoma (HCC) recurrence, but the hereditary factors influencing NM23 levels are unknown. Using public database, the diagnostic value of NM23 in HCC was investigated. A total of 424 hepatitis B virus- (HBV-) related HCC patients were enrolled to perform a genome–wide association study for identifying candidate variants associated with NM23 expression level. Additionally, a logistic regression model, haplotypes, and survival analysis were performed in the subsequent analysis. We identified high NM23 expression levels that have a diagnostic accuracy in HCC tissues and had a poor recurrence-free survival in HBV-related HCC patients. Variants near Psoriasis susceptibility 1 candidate 1 (PSORS1C1) and StAR related lipid transdomain containing 3 (STARD3) are associated with NM23 expression. The PSORS1C1 haplotype TGCACA and the STARD3 haplotype GG have favorable cumulative effects on NM23 expression. Further, variants in PSORS1C1 were associated with either overall survival (rs556285588, rs3095301, and rs3131003) only or overall survival and recurrence-free survival (rs560052000 and rs541820233) both in HCC patients. Our findings suggested that variants at the PSORS1C1 and STARD3 loci play an important role in NM23 regulation. Moreover, variants in PSORS1C1 are potential biomarkers for the prediction of postoperative clinical outcomes in HBV-related HCC patients. Thus, variants in PSORS1C1 and STARD3 are associated with NM23 expression and clinical outcomes of HBV-related HCC patients, which may be regarded as potential biomarkers for this disease.
Primary liver cancer (PLC) is a common malignant neoplasm, with an estimated 854 000 incident cases and 810 000 deaths globally in 2015, contributing to 20 578 000 disability–adjusted life-years . Hepatocellular carcinoma (HCC) comprises 85%–90% of PLC, and the 5-year survival after resection for early–stage HCC ranges from 17 to 53% with recurrence rate as high as 70% [2, 3]. HCC is a multifactorial disease involving a complex interplay between genetic and environmental factors. Epidemiological studies indicated that the major etiological factors affecting HCC include hepatitis B virus (HBV), hepatitis C virus (HCV), aflatoxin exposure, excessive alcohol intake, liver flukes, and cirrhosis [4, 5]. The onset and development of HCC are generally considered to be the consequence of a multistepped process involving the activation of oncogenes and inactivation of tumor suppressor genes.
NME/NM23 nucleoside diphosphate kinase 1 (NM23), also known as NME1, is first reported as an antitumor–metastasis gene that was correlated with metastasis in murine melanoma . NM23 protein expression was associated with cell–cell adhesion, cell migration, proliferation, and invasion depth [7, 8]. Several studies suggested that NM23 expression was inversely proportional to the aggressive metastatic behavior of melanoma as well as gastric, colon, and breast carcinomas . Metastasis is the major cause of morbidity and mortality in individuals with HCC. It has also been reported that the expression of NM23 in tumor tissues is correlated with the occurrence of metastasis and length of survival of HCC patients . Moreover, some studies reported that NM23 expression was upregulated in HCC neoplastic tissue as compared to nontumor tissue [11, 12]. Wei–lu et al.  showed that transcatheter arterial chemoembolization (TACE) enhanced the expression of NM23 in HCC patients. However, little is known about the genetic determinants of NM23 expression in HCC.
Recently, genome–wide association studies (GWAS) have become an efficient method to study the molecular genetics of HCC development and progression . In the present study, we investigated the diagnostic value of NM23 in HCC and performed a GWAS to explore the association between genetic variants and NM23 expression in HBV-related HCC, aiming to identify a novel therapeutic target for NM23 regulation. The NM23 isoforms (NM23-H1 or NM23-H2, also called NME1 or NME2) are heterogeneous in the process of metastasis of HCC; however, NME1 is recommended as one of immunohistochemical markers associated with biological properties of HCC in China. Thus, our study focuses on the expression of this specific isoform.
2. Materials and Methods
2.1. Evaluation of Diagnostic Value for NM23 Expression in HCC
The NM23 specific isoform expression in HCC was obtained from GEO (http://www.ncbi.nlm.nih.gov/geo/) and Oncomine (https://www.oncomine.com/). The criteria used to determine study eligibility were as follows: (1) human species; (2) histopathology confirmed with HCC; (3) availability NME1 expression in HCC and paracancerous; (4) use of prospective or retrospective cohort design with a clearly defined source population and justify all excluded eligible cases; and (5) selection of the latest and most complete study to avoid duplication. NM23 expression in tumor and nontumor tissues was presented as mean and standard deviation (SD) and compared by Student’s test. The receiver operating characteristic (ROC) curve was performed to identify the diagnostic value of NM23 in patients with HCC. The area under the curve (AUC) value was calculated for evaluating the predictive accuracy and discriminative ability of ROC.
2.2. Study Population
A total of 424 patients were enrolled at the First Affiliated Hospital of Guangxi Medical University (Guangxi, China) from 2005 to 2013. All HBV-related HCC subjects were histopathologically confirmed after hepatectomy. This study was approved by the Ethics Committee of the First Affiliated Hospital of Guangxi Medical University.
All tumor tissues were immunohistochemically stained for NME1 by full–time pathologists according to routine processes and the guidelines for standardized pathological diagnosis of primary liver cancer. Supersensitive reagents and mouse anti-human NME1 monoclonal antibody (clone OTI4G3) and streptavidin-peroxidase anti-human kit were purchased from OriGene (Beijing OriGene Technologies, Inc., China). The immunohistochemical staining was carried out by following the manufacturer’s instructions of the kit. For the negative control, the primary antibody was replaced by phosphatebuffered saline/Tween. Positive staining for the NM23 protein appeared as yellow–brown particles.
The criteria used to analyze NM23 expression were based on the number and staining intensity of the stained cells [16, 17]. Briefly, a mean percentage of positive tumor cells was determined in at least five areas at ×400 magnification (50–250 cancer cells per area) and assigned scores as follows: 0: ≤5%; 1: 6%–25%; 2: 26%–50%; 3: 51%–75%; and 4: ≥76%. For convenience of assessment, the intensity of immunostaining was scored following a quantitative principle of proportion: 0, negative, equal to the negative control; 1, weak, cytoplasmic stain slightly darker than the negative control; 2, moderate, defined as an intensity between 1 and 3; 3, strong staining, darker than the positive control. Each sample was processed along with a negative and positive control tissue as references. The sum of the staining intensity and staining extent scores (0–7) was used as the final staining score; that is, a final staining score of 0–1, 2–3, 4–5, or 6–7 was considered to be negative (–), weak (+), moderate (++), or strong (+++) expression, respectively (the representative staining of NM23 in HCC tissues is shown in Figure 1(a)). The results were examined independently and checked collectively by two pathologists who were blinded to the clinicopathological variables. Consequently, the total number of negative, weak, moderate, and strong NM23 expression samples was 43, 285, 48, and 48, respectively.
2.4. DNA Extraction and Genotyping
DNA was extracted from HCC specimens (including paracancer tissues) using the TIANamp Genomic DNA Kit (Tiangen Biotech, Beijing, China) according to the manufacturer’s instructions. DNA yield and purity were measured by the NanoDrop2000 system (Thermo Fisher Scientific, Waktham, MA, USA).
We genotyped all samples on the Illumina Human Exome BeadChip-12-1_A, which includes 242,901 markers focused on protein-altering variants. Genotype calling was carried out using Genotyping Module v1.0 in GenomeStudio. A total of 50 samples (over 10%) were randomly selected using random number table and sequenced for candidate loci by ABI Prism 3100 (Applied Biosystems, Shanghai Sangon Biological Engineering Technology & Services, Shanghai, China), yielding a 100% concordance rate with the genotyped variants (the SNP Sequenced primers were showed in Supplementary Table S1).
2.5. Patients Follow-Up
All patients were followed up after discharge until death or the last time of follow-up (September 2014) by personal or family contacts. The median follow-up time for the 424 patients was 42 months (ranging from 4 to 125 months), and the median survival time (MST) was 30 months. Overall survival (OS) and recurrence-free survival (RFS) were calculated.
2.6. Statistical Analysis
2.6.1. Study Design and Analysis
We performed a GWAS following the process displayed in the flowchart (Figure 1(b)). A large set of single nucleotide polymorphisms (SNPs) was identified using the GWAS array. Then, candidate SNPs were selected using association and pathway analyses. Finally, for these candidate SNPs, the association with the clinical outcomes of HCC patients was evaluated by survival analysis in SPSS.
2.6.2. Quality Control (QC)
Population stratification was estimated by a principal components analysis (PCA), implemented by the EPACTS package in MATLAB 7.0 (MathWorks, Natick, MA, USA). A quantile–quantile (Q–Q) plot was used to evaluate the potential impact of population stratification. QC filtering was performed with PLINK version 1.07, R v. 3.0.1 (https://www.r-project.org/) and the EIGENSOFT package (http://genetics.med.harvard.edu/reich/Reich_Lab/Software.html) as previously described [18, 19]. Finally, candidate SNPs with acceptable quality were analyzed in a discovery GWAS for NM23 expression.
2.7.1. Association Analysis
A linear regression model was used to test the association of quantitative trait (intensity of NM23 expression) with the SNPs that passed QC using the EPACTS package version 3.2.6 , adjusting for age, gender, race, body mass index (BMI), smoking status, drinking status, preoperative TACE, Barcelona Clinic Liver Cancer (BCLC) stage, liver cirrhosis, and Edmondson classification as covariates, and an additive model for allelic effect was assumed. The p–value threshold for the GWAS was 0.01. Baseline variables in intensity of NM23 expression were evaluated using the Chi–square test or Fisher’s exact test. Afterward, we performed ordinal logistic regression analysis to assess the cumulative effect of baseline variables and candidate SNP genotypes on NM23 expression by computing odds ratios (ORs) and 95% confidence intervals (CIs). The reliability of all results was assessed by the test of parallel lines. The correlation of mRNA expression of candidate genes and NM23 was measured by correlation analysis. Differences of mRNA expression of genes between HCC and adjacent normal tissues were verified by Student test.
2.7.2. Bioinformatics Analysis
We selected all genes containing candidate SNPs to search for gene–gene interactions. The signaling pathway network diagram was performed with GeneMANIA Software . Moreover, local linkage disequilibrium (LD) and recombination patterns nearby candidate SNPs were analyzed using LocusZoom  to create regional association plots. The LD values between SNPs were analyzed by the Haploview 4.2 program. The SNPs expression Quantitative Trait Locus (eQTL) mapping was analyzed in GTEx portal (https://www.gtexportal.org).
2.8. Survival Analysis
The Kaplan–Meier lifetable method was used to calculate OS and RFS rates, and differences in survival rates were estimated using a generalized log-rank test. The association between candidate gene and clinical outcomes of HCC patients was analyzed using Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/)  based on the Cancer Genome Atlas (TCGA) database. A Cox proportional hazard regression model was used to calculate hazard ratios (HR) and 95% CIs. All statistical analyses were two–sided and performed using SPSS version 24.0 (SPSS, Chicago, IL, USA). A P–value of less than 0.05 was considered statistically significant. The survival curves were depicted by GraphPad 7.01 (GraphPadSoftware, Inc., Sandiego, CA, USA).
3.1. Public Database Analysis
A total of 28 datasets were enrolled in this study from GEO and Oncomine (Supplementary Table S2). NM23 expression was elevated in HCC tissues as compared to non–tumor tissues in 26 of 28 datasets (Figure 2). The ROC analysis of NM23 expression in HCC datasets indicates high NM23 mRNA level had high accuracy in distinguishing tumor from non–tumor tissues (the AUC of the ROC curves in most of datasets were >0.70, Figure 3).
3.2.1. Baseline Characteristics
The groups were similar with respect to most of the characteristics (P>0.05, Table 1). Age, BMI, regional invasion, and antiviral therapies were significantly different (P<0.05, Table 1) based on the results of the Chi–square test or Fisher’s exact test. The univariate ordinal regression results show BMI≤25 (OR=1.75, 95% CI=1.05–2.93, and Table 1) might be a favorable factor for positive expression of NM23, but ethnic Han ancestry (OR=0.62, 95% CI=0.41–0.95), serum alpha–fetoprotein (AFP)≤400 ng/ml (OR=0.61, 95% CI=0.40–0.93), absence of regional invasion (OR=0.40, 95% CI=0.23–0.67), and antiviral therapies (OR=0.59, 95% CI=0.39–0.90) may associate with lower NM23 expression.
Note. OR for univariate ordinal logistic regression analysis; when the P–value of the OR in the test of parallel lines was less than 0.05, multivariate logistic regression analysis was used but is not shown in the table. †P–value for Chi-square test or Fisher’s exact test. ‡HR and P–values are for univariate survival analysis. PVTT was classified using the radiographic results of enhanced computed tomography scanning as follows: vp1=PVTT in distal to second-order portal branches; vp2=PVTT in second–order portal branches; vp3=PVTT in first–order branches; and vp4=PVTT in the main trunk .
Abbreviations. GWAS: genome–wide association study; OR: odds ratio; 95% CI: 95% confidence intervals; MST: median survival time; MRT: median recurrence time; HR, hazard ratio; Ref.: reference; NA: not applicable; TACE: transcatheter arterial chemoembolization; BMI: body mass index; AFP: alpha–fetoprotein; BCLC: Barcelona Clinic Liver Cancer; PVTT: portal vein tumor thrombus.
A total of 408 HBV-related HCC patients and 21,529 SNPs were included in further analysis after QC filtering. A PCA plot demonstrated that there were no outliers in this study population (Figure 1(c)). The genomic inflation factor (λ) in this study was 1.004 (Figure 1(d)).
3.2.3. Association Analysis
Based on the results of the GWAS displayed in a Manhattan plot (Figure 1(e)), as well as mRNA expression and LD analysis, we identified that SNPs in/near the candidate genes “Psoriasis susceptibility 1 candidate 1” (PSORS1C1) and “StAR related lipid trans domain containing 3” (STARD3) were associated with the expression of NM23 in HBV-related HCC (Supplementary Table S3). The association analysis between these eight SNPs and NM23 expression showed that rs560052000, rs541820233, rs556285588, rs3131003, rs3095301, and rs3095302 were strongly associated with a cumulative effect on NM23 expression (Supplementary Table S4). Moreover, rs11869286–CC was associated with high NM23 expression, whereas rs1877031–GG was associated with lower NM23 expression and rs1877031–AG was associated with it when adjusted than rs1877031–AA.
3.2.4. Pathway Analysis and Correlation Analysis in mRNA
The signaling pathway network showed that NM23 may interact with STARD3 (Figure 4(a)). Additionally, PSORS1C1 may interact with NM23 through T–cell lymphoma invasion and metastasis 1 (TIAM1) and TP53. Further, we used data from the Gene Expression Omnibus (GEO accession: GSE14520) to analyze the mRNA expression of PSORS1C1, STARD3, and NM23 between HCC and adjacent normal tissues. Downregulation of PSORS1C1 was observed in the tumor tissues while STARD3 and NM23 gene expression increased, when compared with the adjacent normal tissues (Figure 4(b)). Correlation analysis was performed to account for the relationship among PSORS1C1, STARD3, and NM23 expression. There was a statistically significant negative correlation between PSORS1C1 and NM23 expression (r=-0.163, P=0.001, Figure 4(c)), but STARD3 and NM23 were positively correlated (r=0.259, P=3.01 × 10−8, Figure 4(d)). PSORS1C1 expression was negatively related to STARD3 expression (r=0.230, P=9.42 × 10−7, Figure 4(e)).
3.2.5. LD and Haplotype Analysis
Haplotype analysis revealed that six SNPs of PSORS1C1 and the two SNPs in STARD3 showed a strong LD block in Haploview (Figures 4(f) and 4(g)). The regional association plots (Figures 4(f) and 4(g)) show the negative log P–values and LD patterns of these SNPs in the combined analysis. The PSORS1C1 haplotypes CAGGTG and CAGACA and the STARD3 haplotypes (CA+GA+CG) had a lower cumulative effect on the expression of NM23 (OR=0.44, 95% CI=0.30–0.65; OR=0.56, 95%=0.32–0.99; OR=0.65, 95% CI=0.47–0.91, respectively) (Table 2), compared to the PSORS1C1 haplotype TGCACA and the STARD3 haplotype GG, respectively.
Note. ★OR, †HR and P–value are adjusted for age, gender, race, smoking status, drinking status, BMI, BCLC stage, Child–Pugh stage, preoperative serum AFP level, TACE status before hepatectomy, pathological grade, cirrhosis, intrahepatic metastasis, PVTT, regional invasion, radical resection, and use of antiviral therapies. Others include CGCACA and CGGGCA. The group of all others is merged CAGGTG, CAGACA, CGCACA, and CGGGCA. Others include CA, GA, CG, and the total number of patients harbour GA and CG is 6.
Abbreviations. OR: odds ratio; HR: hazard ratio; 95% CI: 95% confidence intervals; Ad: adjusted.
3.2.6. eQTL Analysis
The eQTL mapping was performed in a total number of 153 samples in GTEx portal. Rs1877031 and rs11869286 presented a negative eQTL relationship in liver samples (both normalized effect size=-0.25, P=2.1 × 10−7) while rs3095302 and rs3131003 had a positive eQTL association in liver tissues (both normalized effect size=0.54, P=1.3 × 10−6, Supplementary Figure S1).
3.3. Survival Analysis
3.3.1. Distribution of Patient Characteristics and Clinical Outcomes Analysis
MST was calculated using the Kaplan–Meier method with a log-rank test in the different subgroups of baseline variables (Table 1). There were statistical differences among the subgroups in Child–Pugh classification A, BCLC stage A, BCLC stage B, single tumor nodes, tumor size, absence of intrahepatic metastasis, PVTT, radical resection, and insufficient of antiviral therapies in OS of HCC patients. In the RFS, Child–Pugh classification A, BCLC stage A, absence of intrahepatic metastasis, and PVTT had a longer time until recurrence (Table 1).
3.3.2. Association between Candidate Genes and NM23 and Their Complication in HCC Patients
The prognostic value of PSORS1C1, STARD3, and NM23 was evaluated in GEPIA and GSE14520. The cutoff value was set at median in STARD3 and NM23 and quartile in PSORS1C1. High STARD3 and NM23 expression level in patients with HCC had an unfavorable OS and RFS (Figure 5). Besides, high PSORS1C1 expression in HCC patients was a risk factor in OS (HR=1.8, P=0.00023; Figure 5). In HBV-related HCC patients of Guangdong cohort, high NM23 expression level was associated with RFS (Figure 5; HR=1.47, 95% CI=1.01-2.13, Supplementary Table S5).
3.3.3. Association of PSORS1C1 and STARD3 SNPs with Clinical Outcomes in HBV-Related HCC Patients
After hepatic resection, there was a significant difference in the MST of PSORS1C1 SNPs rs541820233, rs556285588, rs560052000, and rs3095301 and when merged two of the genotypes in rs541820233–TT+CC, rs556285588–AA+GG, rs560052000–CC+GG, rs3095301–TC+CC, and rs3131003–AA+AG (Figures 6(a)–6(i), Supplementary Table S5). rs541820233 genotype TC (=0.64, 95% CI=0.47–0.87), rs556285588 genotype AG (=0.62, 95% CI=0.45–0.86), and rs560052000 genotype CG (=0.61, 95% CI=0.44–0.84) were associated with better survival. Additionally, when comparing the groups that were merged the homozygote genotypes, we found that people who harbored heterozygote genotype rs541820233–TC (=0.69, 95% CI=0.50–0.96) and rs560052000 CG (=0.67, 95% CI=0.48–0.93) may have a favorable prognosis in terms of RFS after hepatic resection (Figures 6(j)–6(k), Supplementary Table S6).
3.3.4. Stratification Analysis of PSORS1C1 Associated with Clinical Outcomes
A Cox proportional hazard regression model was applied to perform stratified analysis and further assessed the relationship between the SNPs with OS. We found that PSORS1C1 SNPs rs541820233–TC, rs556285588–AG, and rs560052000–CG were protective factors for most of the clinicopathological and oncological features analyzed (Supplementary Figure S2).
3.3.5. Association of PSORS1C1 and STARD3 Haplotypes with Clinical Outcomes
The association between PSORS1C1 and STARD3 haplotypes and prognosis after hepatic resection was analyzed through a Cox proportional hazard regression model. However, there was no statistically difference between PSORS1C1 haplotypes CAGGTG, CAGACA, and CGCACA+CGGGCA, when compared to TGCACA (Table 2). Similarly, STARD3 haplotypes were not significantly different (CA+GA+CG; =0.97, 95% CI=0.76–1.24; =1.18, 95% CI=0.74–1.88) versus the GG haplotype (Table 2).
In this study, we investigated the diagnostic value of NM23 and performed a GWAS to explore the association between hereditary factors with NM23 expression in HBV-related HCC patients in Guangxi. High NM23 expression showed a precise discrimination in HCC patients, with the AUC of the ROC curves in most of datasets more than 0.70. We detected variants in candidate genes PSORS1C1 and STARD3 fell into two strong LD blocks, with a cumulative effect on the expression intensity of NM23. Interestingly, we found that PSORS1C1 variants associated with clinical outcomes in HBV-related HCC patients. In addition, multivariate Cox proportional hazard regression model analysis demonstrated that PSORS1C1 SNPs rs541820233 and rs560052000 were associated with the OS and RFS of HCC patients, and heterozygous genotypes at these variants were associated with lower risk, implying that these SNPs may be independent prognostic indicators.
HCC is the third leading cause of cancer-related deaths worldwide . Chronic HBV infection affects over 350 million people worldwide and remains one of the leading causes of cirrhosis, liver failure and HCC . Different virus infections may cause variability and complexity of HCC, like HBV and HCV, resulting in different oncology characteristics. Some independent studies suggest that different viral proteins have critical roles in regulating NM23 functions and changing biological activities in cancer progression [26–28]. Run et al.  showed that the disease-free survival in HCC patients with negative NM23 expression was significantly poorer than that in patients with positive NM23 expression. By contrast, report by YB Liu et al. indicated the high NM23 expression in the group with high tendency to metastasis and recurrence and in patients with metastasis or recurrence during the follow-up . In our study, high NM23 expression in HBV-related patients harbored an unfavorable RFS. However, some critical issues as well as host factors, including immunosuppression, somatic mutations, genetic predisposition, and exposure to carcinogens, have important contributory roles . Our result suggested that variants of two immune disease-related genes were associated with NM23 expression in patients with HBV-related HCC.
PSORS1C1 is located at 6p21.3, near the major histocompatibility complex (MHC) class I region. A regional association plot showed that PSORS1C1 lies nearby human leukocyte antigen (HLA)–B and HLA–C. HLA is an important component of the MHC region. One study  reported that PSORS1C1 was associated with HLA–independent systemic sclerosis. Additionally, PSORS1C1 was found to be in strong linkage disequilibrium with the HLA–DQB1 haplotype . In our previous report, we identified that HLA–DQB1 variants associated with OS in HCC patients . The HLA gene family has also been associated with HBV-related HCC in recent studies [33, 34]. Our study showed that PSORS1C1 variants are associated with NM23 expression. Further, we found that the rs560052000–GC, rs541820233–AG, rs556285588–TC, rs3131003–TT+TC, and rs3095301–AA+AG genotypes were associated with better OS in our study subjects. Stratified analysis demonstrated that rs560052000–GC, rs541820233–AG, and rs556285588–TC are protective genotypes in HBV-related HCC patients, having lower HRs for most clinicopathological factors. Although no associations were found in rs3095302 with survival outcomes, all six PSORS1C1 SNPs were in strong linkage disequilibrium. The other haplotypes (CAGGTG, CAGACA, CGCACA, and CGGGCA) had a lower cumulative effect on the expression of NM23 when compared with the TGCACA haplotype. The LD of PSORS1C1 with the HLA region may play a critical role. Tschiedel et al.  reported that NM23 was identified as a novel HLA–A32 restricted tumor–associated antigen in chronic myeloid leukemia. Alternatively, the result of our pathway analysis showed that PSORS1C1 may interact with NM23 through TIAM1 and TP53, which are associated with metastasis. We infer that PSORS1C1 variants in HBV-related HCC patients result in immunosuppression and further cause the HBV to interfere with NM23 expression. However, how PSORS1C1 affects NM23 expression remains unknown, and more experiments are needed to investigate the specific mechanism.
A correlation between STARD3 and NM23 expression was found in the GEO database. STARD3 maps to chromosome 17q11-q12 and encodes a member of a subfamily of lipid trafficking proteins that are characterized by a C–terminal steroidogenic acute regulatory domain and an N–terminal metastatic lymph node 64 domain. Studies have reported that STARD3 expression plays a role in focal adhesion kinase, and correlates with adhesive capacity and prognosis in breast cancer patients [36, 37]. The results of our pathway analysis showed that NM23 and STARD3 may be coexpressed. She et al.  showed that NM23 overexpression reduces the phosphorylation of focal adhesion kinase, mediating the invasive process of hepatocarcinoma cells. These may have a coordinated impact on modulating the cell adhesive capability via upregulated expression, due to mutual gene interactions. We also found that STARD3 variants rs1877031 and rs11869286 were partly associated with NM23 expression. rs1877031 genotype TC has been reported to promote histogenesis in gastric cancer, and STARD3 haplotype CCCT (rs9972882, rs881844, rs11869286, and rs1877031) conferred a protective effect on susceptibility to gastric cancer . In our haplotype analysis, STARD3 (rs11869286, rs1877031) haplotype CA, GA, and CG had a lower cumulative effect on the expression of NM23 (OR=0.65, 95% CI=0.47–0.91), compared to haplotype GG. Moreover, STARD3 colocatizes with EGFR, which may influence clinical phenotypes associated with EGFR by affecting its expression and amplification . Mandai et al.  reported that EGFR is correlated with NM23 expression. According to these studies, STARD3 variants may play a role in regulating NM23 expression through candidate SNPs’ LD effects and effects on EGFR expression, but the specific mechanisms require further investigation.
Several limitations of this study warrant discussion. First, our sample size is modest, as is common to many pharmacogenomic GWAS, and additional studies with larger sample sizes and multiple centers are needed to clarify our results. In addition, because the subjects evaluated in this study included minority subjects, racial heterogeneity may also represent a major limitation of the study. However, we accounted for this by including race, age, and gender as covariates in our GWAS model, and based on the low genomic inflation factor and the Q–Q plot, there is no evidence of population stratification. Finally, our research is preliminary, and further mechanistic and functional studies should be undertaken to discern the potential role of variants near PSORS1C1 and STARD3.
In summary, we identified high NM23 mRNA level offering high diagnostic ability for the discrimination of HCC and demonstrated that genetic variants near PSORS1C1 and STARD3 are associated with NM23 expression in HBV-related HCC. Moreover, variants near PSORS1C1 (rs560052000, rs541820233, rs556285588, rs3095301, and rs3131003) are associated with the clinical outcomes acting as potential biomarkers for the prediction of postoperative patients with HBV-related HCC. The associations and molecular mechanisms of NM23 regulation merit further research.
All the data supporting our findings can be found in the Results and Supplementary materials section of the paper. Please contact corresponding author for more data on reasonable request.
Conflicts of Interest
The authors declare no conflicts of interest.
This work was supported in part by the National Nature Science Foundation of China (Grant nos. 81560535, 81072321, 30760243, 30460143, and 30560133), 2009 Program for New Century Excellent Talents in University (NCET), Guangxi Nature Sciences Foundation (no. GuiKeGong 1104003A-7), and Guangxi Health Ministry Medicine Grant (Key-Scientific Research–Grant Z201018). The authors would also like to acknowledge the support by the National Key Clinical Specialty Programs (General Surgery & Oncology) and the Key Laboratory of Early Prevention & Treatment for Regional High–Incidence–Tumor (Guangxi Medical University), Ministry of Education, China. Thanks are also due to TCGA group, GTEx, GEO, and Oncomine that database provided relevant data.
Supplementary Table S1. Primers of polymerase chain reaction for candidate SNPs. Supplementary Table S2. NM23 expression and ROC analysis of enrolled HCC datasets in this study. Supplementary Table S3. Candidate SNPs in HBV-related HCC patients. Supplementary Table S4. Association between genotypes of candidate SNPs with NM23 expression. Supplementary Table S5. The association between Gene expression levels and clinical outcomes of HBV-related HCC patients. Supplementary Table S6. Association between PSORS1C1 and STARD3 SNPs with clinical outcomes of HBV-related HCC patients after hepatic resection. Supplementary Figure S1. Candidate SNPs eQTL analysis in GTEx portal. Supplementary Figure S2. Stratified analysis on association of rs541820233–TC , rs556285588–AG , and rs560052000–GC with outcomes in HBV-related HCC patients. The HRs are indicated for overall survival. The figure is stratified by favorable and adverse strata. (Supplementary Materials )
- C. Fitzmaurice, T. Akinyemiju, S. Abera 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 Oncology, vol. 3, no. 12, pp. 1683–1691, 2017.
- K. A. McGlynn, L. Tsao, A. W. Hsing, S. S. Devesa, and J. F. Fraumeni Jr., “International trends and patterns of primary liver cancer,” International Journal of Cancer, vol. 94, no. 2, pp. 290–296, 2001.
- G. L. Grazi, G. Ercolani, F. Pierangeli et al., “Improved results of liver resection for hepatocellular carcinoma on cirrhosis give the procedure added value,” Annals of Surgery, vol. 234, no. 1, pp. 71–78, 2001.
- R. N. Aravalli, C. J. Steer, and E. N. K. Cressman, “Molecular mechanisms of hepatocellular carcinoma,” Hepatology, vol. 48, no. 6, pp. 2047–2063, 2008.
- H. B. El-Serag, “Hepatocellular carcinoma,” The New England Journal of Medicine, vol. 365, no. 12, pp. 1118–1127, 2011.
- P. S. Steeg, G. Bevilacqua, R. Pozzatti, L. A. Liotta, and M. E. Sobel, “Altered expression of NM23, a gene associated with low tumor metastatic potential, during Adenovirus 2 Ela inhibition of experimental metastasis,” Cancer Research, vol. 48, no. 22, pp. 6550–6554, 1988.
- M. Boissan, O. De Wever, F. Lizarraga et al., “Implication of metastasis suppressor NM23-H1 in maintaining adherens junctions and limiting the invasive potential of human cancer cells,” Cancer Research, vol. 70, no. 19, pp. 7710–7722, 2010.
- H.-Z. Ren, J.-S. Wang, G.-Q. Pan et al., “Comparative proteomic analysis of β-catenin-mediated malignant progression of esophageal squamous cell carcinoma,” Diseases of the Esophagus, vol. 23, no. 2, pp. 175–184, 2010.
- S. Kapitanović, T. Čačev, M. Berković et al., “nm23-H1 expression and loss of heterozygosity in colon adenocarcinoma,” Journal of Clinical Pathology, vol. 57, no. 12, pp. 1312–1318, 2004.
- R. An, J. Meng, Q. Shi et al., “Expressions of nucleoside diphosphate kinase (nm23) in tumor tissues are related with metastasis and length of survival of patients with hepatocellular carcinoma,” Biomedical and Environmental Sciences, vol. 23, no. 4, pp. 267–272, 2010.
- N. Marino, J. Nakayama, J. W. Collins, and P. S. Steeg, “Insights into the biology and prevention of tumor metastasis provided by the Nm23 metastasis suppressor gene,” Cancer and Metastasis Reviews, vol. 31, no. 3-4, pp. 593–603, 2012.
- J. Cui, B.-W. Dong, P. Liang, X.-L. Yu, and D.-J. Yu, “Construction and clinical significance of a predictive system for prognosis of hepatocellular carcinoma,” World Journal of Gastroenterology, vol. 11, no. 20, pp. 3027–3033, 2005.
- W. Lu, Y.-H. Li, X.-F. He, Y.-Y. Mo, and Z.-Y. Zhu, “Transcatheter arterial chemoembolization enhances expression of Nm23-H1 and TIMP-2 in the tumor tissue of patients with hepatocellular carcinoma,” Hepato-Gastroenterology, vol. 58, no. 106, pp. 558–564, 2011.
- J. U. Marquardt, P. R. Galle, and A. Teufel, “Molecular diagnosis and therapy of hepatocellular carcinoma (HCC): an emerging field for advanced technologies,” Journal of Hepatology, vol. 56, no. 1, pp. 267–275, 2012.
- I. Ikai, Y. Yamamoto, N. Yamamoto et al., “Results of hepatic resection for hepatocellular carcinoma invading major portal and/or hepatic veins,” Surgical Oncology Clinics of North America, vol. 12, no. 1, pp. 65–75, 2003.
- S. Radović, M. Dorić, A. Hukić, M. Babić, S. Kuskunović, and N. Spahović, “Immunohistochemical expression and significance of NM23 suppressor protein in primary gastric adenocarcinoma,” Bosnian Journal of Basic Medical Sciences, vol. 13, no. 2, pp. 72–77, 2013.
- Y. Guan-Zhen, C. Ying, N. Can-Rong, W. Guo-Dong, Q. Jian-Xin, and W. Jie-Jun, “Reduced protein expression of metastasis-related genes (nm23, KISS1, KAI1 and p53) in lymph node and liver metastases of gastric cancer,” International Journal of Clinical and Experimental Pathology, vol. 88, no. 3, pp. 175–183, 2007.
- C. Han, X. Liao, W. Qin et al., “EGFR and SYNE2 are associated with p21 expression and SYNE2 variants predict post-operative clinical outcomes in HBV-related hepatocellular carcinoma,” Scientific Reports, vol. 6, 2016.
- C. K. Yang, T. D. Yu, C. Y. Han et al., “Genome-wide association study of mki67 expression and its clinical implications in hbv-related hepatocellular carcinoma in southern china,” in Proceedings of the Cellular Physiology Biochemistry International Journal of Experimental Cellular Physiology Biochemistry Pharmacology, vol. 42, p. 1342, 2017.
- S. Lee, G. R. Abecasis, M. Boehnke, and X. Lin, “Rare-variant association analysis: study designs and statistical tests,” American Journal of Human Genetics, vol. 95, no. 1, pp. 5–23, 2014.
- D. W. Farley, S. L. Donaldson, O. Comes et al., “The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function,” Nucleic Acids Research, vol. 38, no. 2, pp. W214–W220, 2010.
- R. J. Pruim, R. P. Welch, S. Sanna et al., “LocusZoom: regional visualization of genome-wide association scan results,” Bioinformatics, vol. 26, no. 18, pp. 2336-2337, 2010.
- Z. Tang, C. Li, B. Kang, G. Gao, C. Li, and Z. Zhang, “GEPIA: A web server for cancer and normal gene expression profiling and interactive analyses,” Nucleic Acids Research, vol. 45, no. 1, pp. W98–W102, 2017.
- S. F. Altekruse, K. A. McGlynn, and M. E. Reichman, “Hepatocellular carcinoma incidence, mortality, and survival trends in the United States from 1975 to 2005,” Journal of Clinical Oncology, vol. 27, no. 9, pp. 1485–1491, 2009.
- A. Schweitzer, J. Horn, R. T. Mikolajczyk, G. Krause, and J. J. Ott, “Estimations of worldwide prevalence of chronic hepatitis B virus infection: a systematic review of data published between 1965 and 2013,” The Lancet, vol. 386, no. 10003, pp. 1546–1555, 2015.
- R. Kaul, M. Murakami, K. Lan, T. Choudhuri, and E. S. Robertson, “EBNA3C can modulate the activities of the transcription factor necdin in association with metastasis suppressor protein Nm23-H1,” Journal of Virology, vol. 83, no. 10, pp. 4871–4883, 2009.
- Z. Qin, L. Dai, B. Toole, E. Robertson, and C. Parsons, “Regulation of Nm23-H1 and cell invasiveness by Kaposi's sarcoma-associated herpesvirus,” Journal of Virology, vol. 85, no. 7, pp. 3596–3606, 2011.
- L. Khera, C. Paul, and R. Kaul, “Hepatitis C Virus E1 protein promotes cell migration and invasion by modulating cellular metastasis suppressor Nm23-H1,” Virology, vol. 506, pp. 110–120, 2017.
- Y.-B. Liu, S.-L. Gao, X.-P. Chen et al., “Expression and significance of heparanase and nm23-H1 in hepatocellular carcinoma,” World Journal of Gastroenterology, vol. 11, no. 9, pp. 1378–1381, 2005.
- J. B. Liao, “Cancer issue: Viruses and human cancer,” Yale Journal of Biology Medicine, vol. 79, p. 115, 2006.
- Y. Allanore, M. Saad, P. Dieudé et al., “Genome-Wide scan identifies TNIP1, PSORS1C1, and RHOB as novel risk loci for systemic sclerosis,” PLoS Genetics, vol. 7, no. 7, Article ID e1002091, 2011.
- X. Liu, L. Yu, C. Han et al., “Polymorphisms of HLA-DQB1 predict survival of hepatitis B virus-related hepatocellular carcinoma patients receiving hepatic resection,” Clinics and Research in Hepatology and Gastroenterology, vol. 40, no. 6, pp. 739–747, 2016.
- L. Hu, X. Zhai, J. Liu et al., “Genetic variants in human leukocyte antigen/DP-DQ influence both hepatitis B virus clearance and hepatocellular carcinoma development,” Hepatology, vol. 55, no. 5, pp. 1426–1431, 2012.
- G. Xia, W. Liu, X. Zhang et al., “Genetic polymorphism of hla-dq confers susceptibility to hepatitis b virus-related hepatocellular carcinoma: A case-control study in han population in china,” in Proceedings of the Tumour Biology the Journal of the International Society for Oncodevelopmental Biology Medicine, p. 1, 2016.
- S. Tschiedel, C. Gentilini, T. Lange et al., “Identification of NM23-H2 as a tumour-associated antigen in chronic myeloid leukaemia,” Leukemia, vol. 22, no. 8, pp. 1542–1550, 2008.
- W. Cai, L. Ye, J. Sun, R. E. Mansel, and W. G. Jiang, “Expression of MLN64 influences cellular matrix adhesion of breast cancer cells, the role for focal adhesion kinase,” International Journal of Molecular Medicine, vol. 25, no. 4, pp. 573–580, 2010.
- B. Vassilev, H. Sihto, S. Li et al., “Elevated levels of StAR-related lipid transfer protein 3 alter cholesterol balance and adhesiveness of breast cancer cells: Potential mechanisms contributing to progression of HER2-positive breast cancers,” The American Journal of Pathology, vol. 185, no. 4, pp. 987–1000, 2015.
- S. She, B. Xu, M. He, X. Lan, and Q. Wang, “Nm23-H1 suppresses hepatocarcinoma cell adhesion and migration on fibronectin by modulating glycosylation of integrin beta1,” Journal of Experimental & Clinical Cancer Research, vol. 29, no. 1, article no. 93, 2010.
- Y. Qiu, Z.-Y. Zhang, W.-D. Du et al., “Association analysis of ERBB2 amplicon genetic polymorphisms and STARD3 expression with risk of gastric cancer in the chinese population,” Gene, vol. 535, no. 2, pp. 225–232, 2014.
- R. W. Glynn, N. Miller, and M. J. Kerin, “17q12-21 - The pursuit of targeted therapy in breast cancer,” Cancer Treatment Reviews, vol. 36, no. 3, pp. 224–229, 2010.
- M. Mandai, I. Konishi, M. Koshiyama et al., “Expression of Metastasis-related nm23-Hl and nm23-H2 Genes in Ovarian Carcinomas: Correlation with Clinicopathology, EGFR, c-erbB-2, and c-erbB-3 Genes, and Sex Steroid Receptor Expression1,” Cancer Research, vol. 54, no. 7, pp. 1825–1830, 1994.
Copyright © 2018 Chengkun Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.