Abstract

Background. Circulating long noncoding RNAs (lncRNAs) have been demonstrated to serve as diagnostic biomarkers for various cancers. We aimed to elucidate the diagnostic efficacy of eight serum lncRNAs HULC, MALAT1, Linc00152, PTENP1, PTTG3P, SPRY4-IT1, UBE2CP3, and UCA1 and their combinations for the diagnosis of hepatocellular carcinoma (HCC). Methods. A total of 129 patients with HCC, 49 patients with liver cirrhosis, 27 patients with chronic hepatitis B, and 93 healthy controls were enrolled in this study. The levels of serum lncRNAs were assessed by quantitative real-time polymerase chain reaction. The correlations between serum lncRNAs and clinicopathological characteristics were further analyzed. The receiver operating characteristic (ROC) curve and area under curve (AUC) were utilized to estimate the diagnostic capacity of serum lncRNAs and their combination with AFP for HCC. A logistic regression model was performed to establish a multiple-lncRNA panel. Results. The levels of serum HULC, MALAT1, Linc00152, PTTG3P, SPRY4-IT1, UBE2CP3, and UCA1 were significantly higher in HCC patients than in patients with benign liver diseases and healthy controls, whereas serum PTENP1 was significantly decreased in HCC patients compared with healthy participants. Positive correlations between serum Linc00152 and GGT, serum PTTG3P and GGT, and serum SPRY4-IT1 and ALT were noted in HCC patients. ROC analysis revealed that all these lncRNAs had a significantly predictive value for HCC except for PTENP1. The best performance of single lncRNA was obtained by Linc00152 with an AUC of 0.877. When combined with AFP, the combination of Linc00152 and AFP gained the highest accuracy, yielding an AUC of 0.906. Through logistic regression analysis, the panel consisting of serum linc00152, UCA1, and AFP provided the greatest predictive ability, obtaining an AUC of 0.912 with 82.9% sensitivity and 88.2% specificity. Conclusion. The panel of serum Linc00152, UCA1, and AFP demonstrates a novel and noninvasive biomarker with relatively high sensitivity and specificity for HCC diagnosis.

1. Introduction

Hepatocellular carcinoma (HCC) is one of the most common cancers globally with relatively high morbidity and mortality [1]. Despite great progresses made in HCC therapy in recent years, the prognosis of HCC patients still remains poor due to being diagnosed at an advanced stage and high rate of recurrence [2]. Currently, histological examination of the tumor tissue serves as the standard procedure for definitive diagnosis of HCC, while imaging examinations including ultrasound and computer tomography are the common supplement procedures [3]. However, histopathology may cause hemorrhage and tissue damage, and imaging examinations are expensive and lack adequate sensitivity and specificity. For serum biomarkers, alpha fetoprotein (AFP) is one of the most commonly used biomarkers for HCC screening and diagnosis clinically [4]. Nevertheless, the sensitivity and specificity of AFP are unsatisfactory, especially for patients with early-stage HCC [5]. Therefore, it is of great value to develop novel and noninvasive biomarkers with improved sensitivity and specificity for HCC diagnosis and prognosis.

Accumulating evidence highlights the critical roles of long noncoding RNA (lncRNA) in the initiation and progression of HCC. The lncRNAs, which are defined as noncoding RNA nucleotides in length, have been confirmed to play crucial roles in various biological processes involved in hepatocarcinogenesis, including cell growth, angiogenesis, and metastasis [6]. In addition, previous reports have confirmed the stable presence of lncRNAs in various body fluids such as plasma, serum, urine, gastric juice, and saliva [7]. Several circulating lncRNAs have been implied to be promising markers with high accuracies and efficiencies for the diagnosis and prognosis of HCC [810]. Examples include HULC, MALAT1, Linc00152, UCA1, SPRY4-IT1, and UBE2CP3 [9, 11]. However, the diagnostic value of many lncRNAs for HCC reported by different researchers remains controversial and needs further investigation. Moreover, the diagnostic efficacy of only a small number of circulating lncRNAs has been confirmed. The majority of circulating lncRNAs involved in HCC still remains largely unclear.

Previous studies have reported that serum PTENP1 can distinguish patients with gastric cancer or clear cell renal cell carcinoma from healthy controls [12, 13] and PTENP1 functions as a vital tumor suppressor for HCC [14]. Moreover, we previously found PTTG3P plays an important role in HCC progression and can be detectable in the serum of HCC patients [15]. Therefore, in this study, we initially enrolled PTENP1, PTTG3P, and other six lncRNAs (HULC, MALAT1, Linc00152, SPRY4-IT1, UBE2CP3, and UCA1) which have been already reported to be promising serum markers for HCC [9, 11, 15]. This study is aimed at evaluating the serum levels of these eight lncRNAs among patients with HCC, chronic hepatitis B (CHB), and liver cirrhosis (LC) and healthy participants and at assessing the diagnostic values of these lncRNAs and their combination with AFP, thus to identify a panel of circulating lncRNAs as a tool for the diagnosis of HCC.

2. Materials and Methods

2.1. Patients

A total of 298 individuals were enrolled in the study (129 patients with HCC, 27 patients with CHB, 49 patients with LC, and 93 healthy controls). All patients were recruited from the First Affiliated Hospital of Fujian Medical University between September 2017 and May 2019. The patients with HCC were histologically confirmed. Patients with CHB were diagnosed according to the Clinical Practice Guidelines on the management of hepatitis B virus infection of EASL in 2017. Patients with LC were verified according to the American Association for the Study of Liver Diseases Practice Guidelines. Clinical characteristics of patients were obtained from the medical records retrospectively, including age, gender, HBV surface antigen (HBsAg), tumor number, and tumor size. As a healthy control, 93 individuals who performed their annual health check at the hospital and did not have any liver diseases or other cancerous diseases were recruited. Written informed consent was obtained from each participant, and the study was approved by the Ethics Committee of the First Affiliated Hospital of Fujian Medical University (2018[048]).

2.2. Sample Collection

Peripheral blood samples were collected in a separate vacuum tube from all participants before surgery, pharmacological intervention, or chemotherapy. The serum was isolated by centrifugation at 3,000 rpm for 10 min. After that, the isolated serum was immediately transferred and then stored at −80°C for further analysis.

2.3. RNA Isolation and Quantitative Real-Time PCR (qRT-PCR)

Total RNA was extracted from serum samples using a Hipure Liquid RNA Kit (Cat# R4163-03, Magen, Guangzhou, China). The RNA quantity and purity were evaluated via the NanoDrop One spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The purified RNA was reversely transcribed into cDNA using the M-MLV Reverse Transcriptase (Cat# M1701, Promega, Madison WI, USA). Then, the levels of lncRNAs were assessed by qRT-PCR using TB Green™ Premix Ex Taq (Cat# RR420A, Takara, Dalian, China) which was performed on the QuantStudio Real Time PCR system (Applied Biosystems, Foster City, CA, USA). All reactions were performed with the following conditions: 95°C for 30 seconds, 45 cycles of 95°C for 5 seconds, and 60°C for 30 seconds. The specificity of the PCR products was ensured by melting curve analysis following each reaction. The relative expression of each lncRNA was determined using the 2−ΔΔCt method with GAPDH as the endogenous control for data normalization. The primers used are listed in Table 1.

2.4. Clinical Chemistry Analysis and Serum AFP Detection

Serum alanine aminotransferase (ALT), aspartic transaminase (AST), gamma-glutamyl transpeptidase (GGT), total protein (TP), and albumin (ALB) were measured on the ADVIA 2400 clinical chemistry analyzer (Siemens, Germany). Serum AFP was detected on the Roche Cobas analyzer (Roche Diagnostics, Mannheim, Germany).

2.5. Statistical Analysis

Statistical analysis was carried out by SPSS 23.0 software (SPSS Inc., Chicago, IL, USA) and GraphPad prism 7.0 software (GraphPad Software, Inc., La Jolla, CA, USA). Specifically, the Kolmogorov-Smirnov test was performed to analyze the normal distribution of the variables. If the variables were normally distributed, they were presented as the (SD) and analyzed by two-tailed Student’s -test or one-way ANOVA. Besides, for nonparametric data, the Mann-Whitney test was used for comparisons between groups. Pearson’s correlation coefficient test was applied to judge the correlations between two variables with normal distribution while Spearman’s correlation coefficient test was used for variables with abnormal distribution. Values of the area under curve (AUC) obtained from receiver operating characteristic (ROC) curve were performed to evaluate the diagnostic performance of HULC, PTENP1, MALAT1, Linc00152, PTTG3P, UBE2CP3, SPRY4-IT1, UCA1, AFP, and their combinations. The multi-lncRNA panel was established by logistic regression analysis. A two-tailed value of <0.05 was considered statistically significant.

3. Results

3.1. Patient Characteristics

A total of 298 individuals including 129 patients with HCC, 27 patients with CHB, 49 patients with LC, and 93 healthy controls were enrolled into this study. Table 2 lists the detailed demographic information of all participants. The median age was 59 years (range, 23-88 years) for patients with HCC, 58 years (range, 35-87 years) for patients with LC, 52 years (range, 33-71 years) for patients with CHB, and 55 years (21-79 years) for healthy controls. The HCC cohort consisted of 111 males and 18 females. For patients with LC, 32 patients were male and 17 were female. The CHB group contained 21 males and 6 females. As for healthy controls, 69 participants were male and 24 were female.

3.2. Serum Levels of lncRNAs in Patients with HCC, LC, and CHB and Healthy Controls

As shown in Figure 1, significantly high levels of lncRNA HULC, MALAT1, Linc00152, PTTG3P, SPRY4-IT1, UBE2CP3, and UCA1 were observed in patients with HCC when compared with healthy controls, CHB patients, or LC patients. Serum PTENP1 was markedly decreased in HCC patients than in healthy controls, whereas there was no significant difference with respect to serum PTENP1 between HCC patients and CHB patients or LC patients (Figure 1(d)). Elevated levels of HULC, Linc00152, SPRY4-IT1, UBE2CP3, and UCA1 and the declined level of PTENP1 were observed in patients with LC relatively to healthy controls (Figure 1). Surprisingly, the levels of all eight lncRNAs showed no significance between the LC and CHB groups (Figure 1). Patients with CHB had higher levels of MALAT1, SPRY4-IT1, UBE2CP3, and UCA1 than healthy controls while the levels of other four lncRNAs exhibited no statistical difference (Figure 1).

3.3. Correlations between Serum lncRNAs and Clinicopathological Characteristics

We further investigated the associations between the serum levels of eight lncRNAs and clinicopathological features in patients with HCC. Our results showed that serum Linc00152 were positively correlated with GGT (, ; Table 3, Figures 2(a) and 2(b)). Positive correlations were also observed between PTTG3P and GGT (, ; Table 4, Figures 2(c) and 2(d)). Furthermore, HCC patients with a lower level of serum PTENP1 tended to have a higher level of AFP (, Table 4). Serum levels of serum SPRY4-IT1 were positively associated with ALT (, ; Table 4, Figures 2(e) and 2(f)). HCC patients older than 55 years tended to have a lower level of SPRY4-IT1 and UCA1 (Tables 4 and 5, Figures 2(g) and 2(h)). However, we did not observe any associations between serum lncRNA and gender, HBsAg, tumor size, tumor number, AST, TP, and ALB.

3.4. Diagnosis Efficacy of Each Serum lncRNA in Patients with HCC

Figures 3(a) and 3(b) show the diagnostic value (ROC curves) of 8 lncRNAs in patients with HCC. When HCC patients were tested against healthy controls, the AUC was 0.796 (0.734-0.858) for HULC, 0.768 (0.706-0.830) for MALAT1, 0.895 (0.854-0.936) for Linc00152, 0.602 (0.526-0.678) for PTENP1, 0.785 (0.723-0.847) for PTTG3P, 0.808 (0.750-0.866) for SPRY4-IT1, 0.812 (0.754-0.870) for UBE2CP3, 0.858 (0.810-0.907) for UCA1, and 0.862 (0.815-0.909) for AFP, respectively (Figure 3(a) and Table 6). Among these markers, Linc00152 demonstrated the highest accuracy in predicting HCC from healthy participants with 78.3% sensitivity and 89.2% specificity (Figure 3(a) and Table 6).

We further explored whether the relative levels of the 8 lncRNAs could discriminate HCC patients from patients with CHB and LC. As shown in Figure 3(b) and Table 6, ROC analysis showed that AFP (AUC: 0.811; 95% CI: 0.761-0.861), HULC (AUC: 0.756; 95% CI: 0.702-0.810), MALAT1 (AUC: 0.733; 95% CI: 0.676-0.790), Linc00152 (AUC: 0.877; 95% CI: 0.835-0.918), PTTG3P (AUC: 0.768; 95% CI: 0.715-0.820), SPRY4-IT1 (AUC: 0.768; 95% CI: 0.715-0.821), UBE2CP3 (AUC: 0.756; 95% CI: 0.702-0.810), and UCA1 (AUC: 0.809; 95% CI: 0.761-0.857) had a significantly predictive value for distinguishing HCC patients from CHB and LC patients and HC. However, PTENP1 (AUC: 0.530; 95% CI: 0.465-0.596) did not show any diagnosis efficacy to predicting HCC from CHB patients, LC patients, and healthy controls (, Figure 3(b) and Table 6). The best diagnostic value was obtained with Linc00152 with 81.4% sensitivity and 82.8% specificity (Figure 3(b) and Table 6).

3.5. Diagnosis Efficacy of the Combination of Serum lncRNAs with AFP and Establishment of Diagnostic Panel

Given the fact that AFP is one of the most common markers used for HCC screening and diagnosis, we also evaluated the combinations of these seven lncRNAs with AFP for the diagnosis of HCC. The AUC value, 95% CI, sensitivity, and specificity for the combinations of serum lncRNAs and AFP to discriminate HCC from CHB patients, LC patients, and healthy individuals are listed in Table 7. The AUC values of the combination of seven lncRNAs with AFP to predicting HCC were higher than that of AFP alone (Figure 3(c) and Table 7). The combination of Linc00152 and AFP gained the highest accuracy, yielding an AUC of 0.906 (95% CI: 0.870-0.942) with 85.3% sensitivity and 83.4% specificity (Figure 3(c) and Table 7).

Based on these data, a stepwise regression model was established to get the strongest panel for predicting HCC. Finally, two lncRNAs (Linc00152 and UCA1) and AFP were enrolled. The final formula was listed as follows: . The panel of Linc00152, UCA1, and AFP provided the greatest predictive ability, obtaining an AUC of 0.912 (95% CI: 0.878-0.945) with 82.9% sensitivity and 88.2% specificity (Figure 3(d) and Table 7).

4. Discussion

HCC is an aggressive tumor with high rates of recurrence and metastasis, which leads to the poor prognosis. Therefore, it is of great value to find an effective biomarker for the early diagnosis and therapy of HCC. Previous studies have confirmed that lncRNAs are detectable and also stable in serum of cancer patients, making it possible to be utilized as diagnostic biomarkers for various cancers [16, 17]. Currently, several lncRNAs have been reported to be dysregulated in the serum or plasma of HCC patients. For instance, HULC is highly expressed in tumor tissues as well as plasma of HCC patients, and plasma HULC achieves a fine diagnostic accuracy in predicting HCC [18]. However, clinical use of circulating lncRNAs for HCC still remains controversial due to inconsistent results reported by different researchers and some technical issues, such as using single or multiple lncRNAs.

In this study, we selected eight lncRNAs (HULC, MALAT1, Linc00152, PTENP1, PTTG3P, SPRY4-IT1, UBE2CP3, and UCA1) from previously published HCC-related studies [11, 14, 15, 1820] and assessed their expressions in serum of 129 HCC patients, 49 patients with LC, 27 patients with CHB, and 93 healthy controls. Our data revealed that seven lncRNAs (HULC, MALAT1, Linc00152, PTTG3P, SPRY4-IT1, UBE2CP3, and UCA1) were markedly high expressed in serum of HCC patients compared with healthy controls, which are similar to some previous studies [11, 18, 2123]. These results lead to an indication that these lncRNAs might be useful candidates for the screening and diagnosis of HCC. The elevation of circulating lncRNAs may result from the increased release of nucleic acids from necrosis tumor cells or increased secretion of extracellular vesicles by tumor cells into the blood [17, 24]. PTENP1 is a reported tumor suppressor for HCC [14]. Here, we revealed that serum PTENP1 was significantly downregulated in HCC patients. Moreover, MALAT1, SPRY4-IT1, UBE2CP3, and UCA1 were upregulated in CHB patients, indicating they may be involved in the early period of liver injury. The increased level of all these lncRNAs was observed in LC patients except for MALAT1. Given the fact that LC increases the risk of HCC with an annual incidence between 2 and 4% [3], these seven lncRNAs might play important roles in the early stage of HCC initiation. Intriguingly, plasma MALAT1 is reported to be upregulated in CHB patients and also LC patients [25]. In our study, significant elevation of serum MALAT1 is noted in patients with CHB while only an increasing trend is observed in LC patients. The limited sample size may be possible to be blamed. As for PTTG3P, our previous study has demonstrated the vital roles of PTTG3P in promoting cell growth and metastasis in HCC [15]. Here, we provide the first evidence that the serum level of PTTG3P is markedly high in HCC patients relatively to patients with LC and CHB and healthy controls. This finding provides a basis for the application of PTTG3P as a novel diagnostic biomarker of HCC.

Next, we also investigated the correlations between these circulating lncRNAs and clinicopathological features. A significant association between serum Linc00152 and GGT as well as serum PTTG3P and GGT was demonstrated in patients with HCC. GGT is an enzyme which can be induced by intrahepatic obstruction or directly synthesized by liver cancer cells in HCC [26]. The diagnostic value of GGT has been deeply studied, and elevated GGT has been demonstrated to predict poor survival for HCC [27]. Accordingly, serum linc00152 and PTTG3P might be associated with the prognosis for HCC patients. Our data also revealed an inverse trend between serum PTENP1 and AFP. Furthermore, we found that serum SPRY4-IT1 was positively correlated with ALT. Clinically, ALT is a useful indicator of liver injury. When liver cells are damaged, ALT will leak out from the injured cells into the blood, resulting in an increased serum ALT level [26]. Thus, it is suggested that serum SPRY4-IT1 may also be induced by liver injury. Besides, serum SPRY4-IT1 and UCA1 tended to decrease along with age in patients with HCC. In another word, younger patients tended to have higher levels of serum SPRY4-IT1 and UCA1. The same tendency of UCA1 along with age is reported by Kamel and his colleagues [28]. Some researchers found that young HCC patients have a larger tumor size and more advanced tumor stage at the time of diagnosis and thus show a worse prognosis than elder ones [29]. Therefore, we speculate that high serum levels of SPRY4-IT1 and UCA1 may serve as unfavorable prognostic factors for HCC.

To further investigate the diagnostic values of these lncRNAs, ROC curve analysis was performed. The AUC of the single lncRNA in our study is similar to some previous studies [9, 18]. We firstly investigate the diagnostic values of serum PTENP1 and PTTG3P for HCC. Serum PTTG3P was demonstrated to obtain a fine diagnostic accuracy in distinguishing HCC from LC patients, CHB patients, and healthy controls while the diagnostic potential of serum PTENP1 seems unsatisfactory. The AUC of Linc00152 was greater than those of any other circulating lncRNAs and the currently used tumor marker AFP, indicating remarkable efficacy of Linc00152 in predicting HCC from either benign liver diseases or healthy controls. Nowadays, serum AFP is the most widely used biomarker for HCC diagnosis [4]. However, its sensitivity and specificity are limited. A cohort study of 665 HCC patients carried out by Agopian et al. reported that serum AFP levels do not elevate in 31.3% of patients with HCC [30]. Therefore, we further investigate whether the combination of these lncRNAs and AFP could improve the diagnostic efficacy of single AFP. Our data revealed that the AUC of all the combinations was greater than single AFP. The best diagnostic performance was obtained in the combined use of linc00152 and AFP with an AUC of 0.906.

Previously, most of the HCC-related lncRNAs focus on single molecule when exploring the potential of lncRNA as novel biomarkers. However, the diagnostic accuracy of any single index is limited. Several studies have already developed several multiple biomarker panels to increase diagnostic accuracy, such as the lncRNA panel established by lncRNA-LET, PVT1, PANDAR, PTENP1, and linc00963 for the diagnosis of clear cell renal cell carcinoma [13] and a panel consisting of two lncRNAs (SOX2OT and ANRIL) and three tumor markers (CEA, CYFRA21-1, and SCCA) for the detection of non-small-cell lung cancer [31]. To elevate diagnostic capacity, we identified and combined three markers (linc00152, UCA1, and AFP) by a stepwise selection model, resulting in an AUC of 0.912. The diagnostic sensitivity and specificity of this panel were higher than any single markers. Interestingly, linc00152 and UCA1 have previously been linked with cancers by others. Accumulating evidence support an important role of linc00152 in glioma [32], bladder cancer [33], and HCC [34]. Linc00152 has been demonstrated to promote cell cycle progression in HCC through the miR-193a/b-3p/CCND1 axis [35]. Additionally, elevated serum linc00152 is found in patients with gastric cancer [36], non-small-cell lung cancer [37], and also HCC. The AUC for linc00152 to predicting HCC from healthy control reported by Li et al. [18] is 0.85, which is similar to our results. Similarly, frequent upregulation of UCA1 is observed in various cancers, including lung cancer, bladder cancer, gastric cancer, liver cancer, and colorectal cancer [38]. Xiao et al. [39] reported that UCA1 regulates Snail2 expression by effectively sponging miR-203 to promote HCC progression. Moreover, several studies provide strong evidence that serum UCA1 is a potential marker for HCC diagnosis and prognosis [23, 28, 40]. Consistent with our findings, Zheng et al. [23] revealed that serum UCA1 could differentiate patients with HCC from healthy participants and benign liver diseases. These studies further strengthen our findings that circulating lncRNA linc00152 and UCA1 have high potential to act as HCC biomarkers.

There are a few limitations to our study. First, we did not validate the diagnostic efficacy of the panel in another set. Second, our study was performed at a single center with relatively limited sample size. Thus, further investigations and large-scale multicenter studies are strongly recommended to fully validate this novel panel. Third, apart from AFP, AFPL3 and PIVKA-II are also common markers for HCC diagnosis. The comparisons between the diagnostic values of these serum lncRNAs and AFPL3 and PIVKA-II for HCC are also recommended for a further study.

5. Conclusion

In summary, we established a diagnostic panel consisting of serum lncRNA linc00152, UCA1, and tumor marker AFP with superior sensitivity and specificity to diagnose HCC. This panel might serve as a novel and noninvasive biomarker for HCC diagnosis.

Abbreviations

AFP:Alpha fetoprotein
ALB:Albumin
ALT:Alanine aminotransferase
AST:Aspartic transaminase
AUC:Area under curve
CHB:Chronic hepatitis B
CI:Confidence interval
GGT:Gamma-glutamyl transpeptidase
HBsAg:HBV surface antigen
HC:Healthy control
HCC:Hepatocellular carcinoma
HULC:Highly upregulated in liver cancer
LC:Liver cirrhosis
lncRNA:Long noncoding RNA
MALAT1:Metastasis-associated lung adenocarcinoma transcript 1
PTENP1:Phosphatase and tensin homolog pseudogene 1
PTTG3P:Pituitary tumor-transforming 3, pseudogene
qRT-PCR:Quantitative real-time polymerase chain reaction
ROC:Receiver operating characteristic curve
SPRY4-IT1:SPRY4 intronic transcript 1
TP:Total protein
UBE2CP3:Ubiquitin conjugating enzyme E2 C pseudogene 3
UCA1:Urothelial carcinoma-associated 1.

Data Availability

The data that support the findings of this study are not publicly available due to privacy and ethical restrictions but are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

Authors’ Contributions

Jinlan Huang and Yansong Zheng contributed equally to this work.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (81802087), the youth project of Fujian Provincial Health and Family Planning Commission (2017-1-47), the Natural Science Foundation of Fujian Province (2018J05130), and the Qihang Project of Fujian Medical University (2017XQ1053). We thank the biomedical laboratory scientists for the excellent technical assistance with handling of blood samples. We also thank the patients for their participation in this study.