Abstract

Background and Objectives. Lymph node metastasis (LNM) is common in hepatocellular carcinoma (HCC). In order to intervene HCC LNM in advance, we developed a prediction nomogram based on serum long noncoding RNA (lncRNA). Methods. Serum samples from 242 HCC patients were gathered and randomly enrolled into the training and validation cohorts. LncRNAs screened out from microarray were quantified with qRT-PCR. Univariate and multivariate analyses were applied for screening independent risk factors. A prediction nomogram was ultimately developed for HCC LNM. The nomogram was estimated by discrimination and calibration tests in the validation cohort. The effects of the candidate lncRNA on the malignant phenotypes of HCC cells were further explored by wound healing assay and colony formation assay. Results. ENST00000418803, lnc-ZNF35-4:1, lnc-EPS15L1-2:1, BCLC stage, and vascular invasion were selected as components of the nomogram according to the adjusted multivariate analysis. The nomogram effectively predicted the HCC LNM risk among the cohorts with suitable calibration fittings and displayed high discrimination with C-index of 0.89 and 0.85. Moreover, the abnormally high expression of lnc-EPS15L1-2:1 in HCC cell lines showed significant carcinogenic effects. Conclusions. The noninvasive nomogram may provide more diagnostic basis for treatments of HCC. The biomarkers identified can bring new clues to basic researches.

1. Introduction

Liver cancer is the second leading cause of cancer-related death worldwide, with approximately 850, 000 new cases occurring each year [1]. Among all the cases of liver cancer, hepatocellular carcinoma (HCC) accounts for about 90% [2]. Due to the poor prognosis of HCC, its mortality ranks third among all cancer types and shows an upward trend on a global scale [3, 4]. Compared with other prognostic factors, metastasis is one of the main reasons for cancer-related death [5]. Thereinto, lymph node metastasis (LNM) is common within extrahepatic metastases, second only to lung metastasis [6]. It was reported that the intraoperative detection rate of HCC LNM was 0.75%–7.5%, while the detection rate during autopsy was as high as 30.3% [6]. In the previous follow-up, we have found that around 10.3% of HCC patients after hepatectomy will develop LNM [7]. Although external beam radiotherapy is an effective means of treating HCC LNM, the prognosis of patients with LNM is worse than that of patients without LNM [8, 9].

Until now, the mechanisms underlying the development of LNM in HCC are poorly understood. Screening for biomarkers associated with HCC LNM can effectively identify high-risk patients, thereby taking steps to prevent disease progression after curative resection. In the previous studies, we have identified multiple biomarkers that are closely related to HCC LNM, such as hypoxia inducible factor-1α (HIF-1α), vascular endothelial growth factor (VEGF), and matrix metalloproteinase-2 (MMP-2)[7]. At this stage, we urgently hope to find out noninvasive and specific diagnostic markers to prevent the occurrence of HCC LNM.

Researches have found that lncRNAs are involved in multiple steps of cancer development, and they can be used as sensitive biomarkers to predict cancer metastasis [10, 11]. Moreover, lncRNAs have been identified in body fluids [1214] and are emerging as novel biomarkers for disease and as targets for disease intervention [15, 16]. Based on the expression profile of serum lncRNAs in HCC LNM [17], we quantitatively analyzed the candidate lncRNAs in the serum samples of 242 patients and evaluated the correlation between the lncRNAs and HCC LNM. The biomarkers with potential predictive efficacy were finally screened out, and a prediction nomogram for HCC LNM based on serum lncRNAs was constructed for the first time in combination with clinicopathological parameters.

By applying prediction nomogram for noninvasive risk assessment, HCC patients at high risk for LNM could receive prophylactic radiotherapy towards regional lymph nodes, which could effectively reduce the incidence of LNM and prolong the overall survival of HCC patients [1820]. Moreover, the abnormal expression and function of the HCC LNM-related serum lncRNA at cellular level may provide a basis for the mechanism study.

2. Materials and Methods

2.1. Patients and Serum Samples

The enrolled patients received treatments at Zhongshan Hospital, Fudan University, from 2012 to 2016 and were pathologically diagnosed after surgery. The blood samples from the enrolled patients were collected at the time of surgery and the serums were subsequently extracted with centrifugation (1000×g, 10 min, 4°C). The serum samples were than stored at -80°C immediately. All the enrolled patients had complete medical history and clinicopathological data with follow-up and regular clinical examinations and were randomly assigned to the training and validation cohorts in a ratio of 2 : 1. The grouping and exclusion of the enrolled patients are shown in Figure 1. Among the included clinical indicators, the grade of tumor differentiation was determined using the Edmondson grading system and tumor size was measured depending on the maximum diameter of the tumor specimen. The vascular invasion was evaluated with microscopic examination of resected specimens.

All the procedures performed in this study involving human participants were in accordance with the ethical standards of the institutional and national research committee.

2.2. Follow-Up and Postoperative Treatments

All the enrolled patients received routine follow-up every 3 months after surgery, which contains clinical examinations and medical history collection. Clinical examinations are composed of abdominal ultrasound examination, liver function test, and blood routine measurement, and all the examinations were performed by doctors who were blind to this study. When HCC LNM was suspected, CT scanning or MRI was carried out immediately for diagnosis. Once diagnosed with HCC LNM, patients were treated with radiotherapy for regional lymph nodes.

2.3. RNA Isolation and qRT-PCR

In order to relatively quantify the target lncRNAs, cel-mir-39-3p and GAPDH were used as the external and inner reference, respectively. Total RNAs were isolated from the serum samples using TRIzol LS reagent according to the manufacturer’s instructions. Specifically, 0.75 ml TRIzol LS reagent (Life technologies, Carlsbad, CA, US) was added into 0.25 ml serum sample, and 50 fmol of cel-miR-39 was then added for normalization. Cell-derived RNAs were isolated using TRIzol reagent according to the manufacturer’s instructions. After assessing the purity of the extracted total RNA, the process of reverse transcription was performed using Prime Script RT reagent Kit (Takara Bio, Shiga, Japan). The primer sequences were listed in Table 1. QRT-PCR analyses were conducted using SYBR® Premix Ex Taq™ (Takara Bio, Shiga, Japan) in the 7500 Real-Time PCR System (Applied Biosystems). The relative expression levels of the candidate lncRNAs were normalized against reference standards by utilizing method.

2.4. Cell Culture

Human hepatocyte cell line QSG-7701 and human HCC cell lines SMMC-7721, HuH-7, HepG2, MHCC-97H, MHCC-97L, and HCC-LM3 were cultured in DMEM containing 10% fetal bovine serum in a constant temperature incubator at 5% CO2, 37°C.

2.5. Fluorescence In Situ Hybridization

The fixed cells were infiltrated with citrate buffer at room temperature and placed in the cell punching solution for 10 minutes. The cells were immersed in prehybridized working solution at 42°C in the dark and incubated for 2 hours in a wet box. The sample was washed with 0.2×sodium citrate buffer and the hybridization working fluid containing lncRNA probe was then added for overnight rest. After removing the cover glass, the sample was washed at 37°C with 2×sodium citrate buffer, 0.2×sodium citrate buffer, and PBS-T solution. The nuclei were stained with DAPI for 10 minutes at room temperature and then encapsulated. Finally, the cells were observed and photographed by confocal laser scanning microscopy in a field in which the cell division phase was observed.

2.6. Lentivirus Transfection

An appropriate amount of lentivirus suspension (moi=30) was added at a cell density of 2×105 /well and cultured at 5% CO2, 37°C. After 48 hours, the expression of fluorescent protein TurboGFP was detected under the microscope and the transfection efficiency was confirmed. The stable transfectant was screened by using 2 μg/ml concentration of puromycin.

2.7. Wound Healing Assay

The back side of the 6-well plate was marked with a horizontal line and uniformly plated at a cell density of 5×105 per well. After cellular fusion, the vertical lines perpendicular to the marking lines were made with pipette tip in the cell distribution area. The 6-well plate was then washed three times with PBS buffer to ensure the complete removal of exfoliated cells, and the medium with low concentration of serum (2%) was also added. The treated 6-well plate was cultured in a constant temperature incubator at 5% CO2 and 37°C and photographed under a microscope at 0, 24, 48, and 72 hours.

2.8. Colony Formation Assay

Four hundred and eight hundred cells were laid in a 6 cm culture dish and cultured in a constant temperature incubator at 5% CO2 and 37°C. Two weeks later, visible clones were formed in the culture dish and more than 50 cells were found in most single clonal colonies. After removing the culture medium, PBS was gently washed along the medial wall of the dish three times, and 4% polyformaldehyde was used to fix the cells for 30 minutes. After one hour of continuous dyeing with crystal violet dye solution, rinse was performed along the side wall of the dish until clear dark blue colony spots appeared. The colonies containing more than 50 cells were counted three times in each dish, and the average value was taken as the final clone number. Quantitative analysis was carried out by GraghPad Prism software, and the clone formation rate was equal to the number of colonies/inoculated cells.

2.9. Statistical Analysis

By using SPSS 22.0 software (Chicago, IL, USA), GraphPad Prism 6.01 (La Jolla, CA, USA), Image J (National Institutes of Health, NIH), MedCalc Statistical Software 18.2.1 (Ostend, Belgium), and R software (version 3.2.3), all the statistical data were shown from at least three separate experiments. Specifically, the correlations between clinicopathological features and HCC LNM were analyzed by Pearson’s correlation test and Fisher’s exact test. Univariate analysis and multivariate Cox regression analysis (confounding variables that affect the regression coefficient by more than 10% were adjusted) were performed to determine independent risk factors of HCC LNM. The model was finally determined for nomogram with a backward step-down selection process. Estimates of the cutoff value and the model discriminatory ability were measured using time-depended ROC curves. Moreover, calibration curves originating from Hosmer-Lemeshow (H-L) test were applied for evaluation of the nomogram. The P < 0.05 were normally considered statistically significant. Vascular invasion was also included as an independent risk factor with slightly higher P value (P = 0.055) and appropriate hazard ratio.

3. Results

3.1. Patients’ Background Data

Characteristics of all the enrolled HCC patients (157 in the training cohort and 85 in the validation cohort) used to establish a prediction model are summarized in Table 2. Patients of the training cohort were observed until October 2016, and the median follow-up time was 55 months (range, 9-76 months). For the validation cohort, observation lasted until September 2016 and the median follow-up time was 45 months (range, 10-76 months). During the follow-up time, 20 patients (12.7%) in the training cohort and 9 patients (10.6%) in the validation cohort developed LNM.

3.2. Screening of Serum Biomarkers with Microarray

Through previous microarray analyses on serum samples, we have identified 235 lncRNAs that were differentially expressed among LNM group and non-LNM group as reported [17]. Among all the significant differential lncRNAs, five had fold change > 2 and FDR values < 0.05 with high expression abundance in serum samples. Lnc-GALR2-1:1 was downregulated with a FDR value of 0.034; ENST00000418803 was downregulated with a FDR value of 0.037; lnc-ZNF35-4:1 was downregulated with a FDR of 0.042; lnc-CAMKK2-3:2 was upregulated with a FDR of 0.037; and lnc-EPS15L1-2:1 was upregulated with a FDR of 0.037. We therefore conducted qRT-PCR to examine the expression of the candidate lncRNAs in serum samples from the training cohort and further evaluated the ability of these potential biomarkers for predicting HCC LNM.

3.3. Analysis of Serum Candidate lncRNAs Expression

With cel-mir-39-3p as an external reference, qRT-PCR was carried out to quantify the relative expression of the candidate lncRNAs in the serum samples. We first assessed the relevance between the serum lncRNA expression and diagnosis of HCC LNM based on the ROC curve and determined the cutoff value for judging lncRNA expression level simultaneously. Relative expressions of lnc-GALR2-1:1 and lnc-CAMKK2-3:2 in serum samples from the enrolled patients were verified to have no statistical significance for predicting HCC LNM. Besides, the remaining three lncRNAs were defined as high- or low-expression level based on the maximum value of Youden’s index in the ROC analysis [21]. As showed in Figure 2(a), the AUC of ENST00000418803, lnc-EPS15L1-2:1, and lnc-ZNF35-4:1 were 0.753 (95% CI: 0.678-0.819), 0.721 (95% CI: 0.644-0.790), and 0.766 (95% CI: 0.692-0.830), respectively, which indicated a considerable distinguishing power to HCC LNM. To further determine the diagnostic relevance between the above three lncRNAs and HCC LNM, Kaplan-Meier and Log-rank tests were performed among all patients from the cohorts as Figures 2(b), 2(c), and 2(d) show. The distinctions between non-LNM and LNM groups in the training cohort divided by the optimum cutoff values were displayed in Figure 3. The low expression of ENST00000418803 and lnc-ZNF35-4:1 and the high expression of lnc-EPS15L1-2:1 showed significant correlations with the occurrence of LNM.

In the training cohort of 157 patients, low ENST00000418803 expression was found in 51 of 157 patients (32.5%), high lnc-EPS15L1-2:1 expression in 68 (43.3%), and low lnc-ZNF35-4:1 expression in 55 (35.0%). In the validation cohort, low ENST00000418803 expression was found in 30 of 85 patients (35.3%), high lnc-EPS15L1-2:1 expression in 33 (38.8%), and low lnc-ZNF35-4:1 expression in 27 (31.8%).

3.4. Significant Predictors of Lymph Node Metastasis in Hepatocellular Carcinoma

For the training cohort, 19 clinicopathological features consisted of age, gender, HCV-Ab, HBsAg, a-fetoprotein, tumor differentiation, Child-Pugh score, intrahepatic metastasis, tumor size, vascular invasion, BCLC staging, ALT, γ-GT, liver cirrhosis, tumor number, and distant metastasis and the expression of the serum lncRNAs mentioned above was considered for the univariate analysis with Cox proportional hazards regression. The associations of the included variables with HCC LNM in the training cohort are summarized in Table 3. Through univariate analysis, vascular invasion (P = 0.010), BCLC stage (P = 0.023), ENST00000418803 (P < 0.001), lnc-EPS15L1-2:1 (P < 0.001), and lnc-ZNF35-4:1 (P = 0.001) were selected to be significantly associated with LNM in HCC patients, whereas age (P = 0.286), gender (P = 0.187), HBsAg (P = 0.248), HCV-Ab (P = 0.997), a-fetoprotein (P = 0.539), tumor differentiation (P = 0.611), Child-Pugh score (P = 0.998), intrahepatic metastasis (P = 0.365), tumor size (P = 0.509), alanine aminotransferase (ALT) (P = 0.456), γ‐glutamyltransferase (γ-GT) (P = 0.878), liver cirrhosis (P = 0.756), tumor number (P = 0.585), and distant metastasis (P = 0.752) displayed no significant association with LNM. Statistically significant variables were further adopted for multivariate analysis. By adjusted multivariate analysis, the following five variables were found to be independent risk factors for LNM in HCC: vascular invasion (P = 0.055, HR: 2.5, 95% CI: 1.0 ~ 6.5), BCLC stage (P = 0.007, HR: 4.2, 95% CI: 1.5 ~ 12.0), ENST00000418803 (P < 0.001, HR: 0.2, 95% CI: 0.1 ~ 0.5), lnc-EPS15L1-2:1 (P = 0.001, HR: 8.7, 95% CI: 2.5 ~ 30.6), and lnc-ZNF35-4:1 (P = 0.025, HR: 0.3, 95% CI: 0.1 ~ 0.9).

3.5. Construction of Prediction Nomogram for LNM in HCC

As shown in Figure 4, the following five independent risk variables from multivariate cox regression analyses were selected into the visible nomogram to predict the risk of LNM: high lnc-EPS15L1-2:1 expression has the highest score of 100; low ENST00000418803 has the score of 73; low lnc-ZNF35-4:1 has the score of 51; BCLC stage; and vascular invasion was scored as 54 and 34, respectively. Considering the time distribution of LNM occurrence and the time-dependent AUC in the both cohorts, we determined the 29 months’ time (75% quartile) as an observation point, in which the nomogram has the best prediction performance and application value. The sum of score from all the included risk factors can further correspond to the risk assessment of LNM occurring within 29 months. The prediction nomogram demonstrated a good accuracy with stable and favourable time-dependent AUC among the training cohorts (Figure 5(a)). Calibration curves revealed a suitable calibration between the predictive LNM risk and the observed LNM risk as well (Figure 5(c)). Harrell’s C-index of the stepwise selected model for LNM prediction was 0.89, which indicated a sound discrimination ability.

3.6. Validation for the Predictive Value of the lncRNA-Based Nomogram

As Figure 5(b) shows, the trend of AUC was stable and satisfied in the interval of 15 to 31 months’ observation time point, which revealed a good prediction performance for LNM within the optimum time (29 months’ time). The calibration curve noted that the nomogram was well calibrated with a favorable fitting between the observed and predicted risk (Figure 5(d)). Furthermore, Harrell’s C-index 0.85 of the model demonstrated good discrimination in the validation step.

3.7. Role of Lnc-EPS15L1-2:1 in Lymph Node Metastasis of HCC Cells

The aforementioned analyses showed that the high expression of serum lnc-EPS15L1-2:1 had a strong correlation with the occurrence of lymph node metastasis in HCC patients. However, the source of free lnc-EPS15L1-2:1 and the biological functions associated with HCC LNM remain unknown.

We performed qRT-PCR on the expression of lnc-EPS15L1-2:1 in human hepatocyte line QSG-7701 and human HCC cell lines SMMC-7721, Huh-7, HepG2, MHCC-97H, MHCC-97L, and HCC-LM3. Among these, SMMC-7721, HuH-7, MHCC-97H, and HCC-LM3 cell lines have high invasion and metastasis characteristics, which are highly malignant. On the contrary, HepG2 and MHCC-97L cell lines show low malignancy in invasion and metastasis. As shown in Figure 6, the expression of lnc-EPS15L1-2:1 in HCC cells was higher than that in normal hepatocytes and was significantly elevated in highly malignant cell lines. Fluorescence in situ hybridization (FISH) further confirmed the subcellular distribution of lnc-EPS15L1-2:1 which was mainly cytoplasmic (Figure 7).

The formation of tumor metastases depends on both high invasion and migration potential and strong clone-forming ability. Therefore, we performed wound healing assay and colony formation assay after the overexpression of SMMC-7721 cell line by lentiviral transfection. As shown in Figures 8 and 9, lnc-EPS15L1-2:1 significantly enhanced the migration and clonality of HCC cells.

4. Discussion

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors worldwide. The incidence of LNM in extrahepatic metastases of HCC is approximately 33.8%, and the overall survival of untreated HCC patients with LNM is only about three months [6, 22]. Lymph node metastasis is a clear prognostic factor in treatments of cancer and has a significant impact on long-term survival of patients [23]. However, imaging technique is still not sensitive to the initial diagnosis of lymph node micrometastasis. It is necessary to build a specific model that predicts LNM risk in HCC patients for preventive intervention.

In general, the free nucleic acid in the circulation of tumor patients is mainly derived from the frequent apoptosis and necrosis of cells [24, 25]. Therefore, the abnormal expression of some cell-free RNA can partly reflect the expression profile of cancer cells, which is significantly related to the malignant process of cancer [26]. Aberrant expressions of lncRNAs have been reported to participate in diverse biological processes in cancer, including LNM [2729]. Furthermore, lncRNAs are becoming a type of potential biomarkers for disease prediction [12, 30]. In the current study, we therefore sought to identify sensitive serum lncRNAs as noninvasive biomarkers to predict LNM in HCC.

We constructed and then validated a novel prediction nomogram based on three lncRNAs to predict LNM for HCC patients after hepatectomy. Our lncRNA-based nomogram incorporates the following five independent prognostic factors: BCLC stage, vascular invasion, low expression of ENST00000418803 and lnc-ZNF35-4:1, and high expression of lnc-EPS15L1-2:1. Using a linear predictor of 1.7 as the optimum cutoff value, 18.5% (n = 29) patients in the training cohort are identified as high-risk group, and 48.3% (n = 14) patients in the high-risk group developed LNM within the follow-up period. When applying the nomogram to the verification step, 14.1% (n = 12) patients in the validation group were allocated to the high-risk group, and 41.7% (n = 5) of whom were diagnosed with LNM. To our knowledge, this is the first report of serum lncRNA-based prediction nomogram in HCC, especially with respect to LNM. Utilizing this nomogram, posthepatectomy HCC patients can be accurately classified into groups with low and high risk of LNM at the early stage. The nomogram could serve as a valuable criterion for determining optimal treatment strategies for HCC patients.

We have two conjectures about the source of the lncRNAs that have potential predictive effects on HCC LNM. One possibility is that the abnormal cell-free lncRNAs may originate from the transformation process of premetastatic niches mediated by circulating tumor cells or tumor exosomes. NONCODE database displayed the notion that lnc-EPS15L1-2:1 are significantly highly expressed in lymph node compared with other tissues [31]. Another possibility is that changes in the expression profile of lncRNAs in tumor cells promote metastatic propensity. Cis analyses of protein-coding genes adjacent to the lncRNAs loci uncovered their possible roles in tumor malignant phenotype. Specifically, Kruppel-like factor 2 (KLF2), associated gene of lnc-EPS15L1-2:1 from Cis and Trans analyses, was reported as a terminal component of tumor proliferation and metastasis related pathway axis [3234]. And ENST00000418803 associated gene Sad1 and UNC84 domain containing 2 (SUN2) act as novel suppressors in cancer [35, 36].

There is a strong correlation between the high expression of serum lnc-EPS15L1-2:1 and the occurrence of HCC LNM. Our further experiments confirmed that lnc-EPS15L1-2:1 is highly expressed in HCC cells compared with normal hepatocytes and is associated with the malignancy of HCC cells. Moreover, lnc-EPS15L1-2:1, mainly distributed in the cytoplasm, significantly promotes the migration and clonality of HCC cells. We are about to investigate the specific role of lnc-EPS15L1-2:1 in HCC LNM in the next stage of research.

Some limitations should also be acknowledged as follows. As it was a retrospective cohort study, and because of the limited quantity of patients involved, the results need to be further validated in a large scale of prospective study.

In conclusion, the nomogram based on serum ENST00000418803, lnc-ZNF35-4:1, lnc-EPS15L1-2:1, BCLC stage, and vascular invasion has good predictive performance for HCC LNM, and HCC patients at high risk of LNM may benefit from this. Moreover, as an independent risk factor, overexpression of lnc-EPS15L1-2:1 may mediate the occurrence of HCC LNM at the cellular level, which needs further verification. This study will provide new ideas for the clinical diagnosis and mechanism research of HCC LNM.

Abbreviations

HCC:Hepatocellular carcinoma
LNM:Lymph node metastasis
lncRNA:Long noncoding RNA
qRT-PCR:Quantitative reverse transcription polymerase chain reaction
HIF-1α:Hypoxia inducible factor-1α
VEGF:Vascular endothelial growth factor
MMP2:Matrix metalloproteinase-2
BCLC:Barcelona Clinic Liver Cancer
AFP:Alpha-fetoprotein
ALT:Alanine aminotransferase
γ-GT:γ-glutamyltransferase
ROC curve:Receiver operating characteristic curve
AUC:Area under the curve.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

Zhao-Chong Zeng and Zuo-Lin Xiang contributed equally to this work.

Acknowledgments

This study is supported by Natural Science Foundation of Shanghai (Grant no. 17ZR1405300), Science and Technology supporting project of Shanghai (Grant no. 17411962600), Shanghai Municipal Human Resources and Social Security Bureau (Grant no. Q2016-019), and Pudong New Area Science and Technology Development Fund (Grant no. PKJ2018-Y02).